通过autoff得出的原子 sigma， epsilon 参数与文献不一样

huangzi

老师们，我用autoff得出了tfsi-离子的力场，这是获得的原子的 sigma， epsilon 参数 ，然后我下午偶然看到一篇文献其中原子的 sigma， epsilon 参数为 我仔细对比了键角二面角发现两者一致，只有f和o的不符合，想问一下老师我应该将autoff导出来f、o的sigma， epsilon 参数修改成文献里的吗

student0618

文献是按什么优化这参数的？

huangzi

student0618 发表于 2025-12-30 21:10
文献是按什么优化这参数的？

https://pubs.acs.org/doi/10.1021/jp0476545     老师这是参考的这篇文献 ，，其中有一段这个

yuechenghua

oplsaa里面F的种类有好几种，lj参数之间有差异。
opls_786 是PF6-里的F，sigama和epsilon分别是3.11810e-01 、2.55224e-01
opls_965是氟代烷烃里的F，sigama和epsilon分别是2.95000e-01 、2.21752e-01，这个就是你给的文献里的参数。
我不知道你的tfsi-离子是什么结构，你可以自己根据化学环境判断下。
然后O的参数看起来判断成了opls_441，磷酸酯里的双键O，你给出的文献更像opls_278里的C=O双键O，还是一样，找跟你结构环学环境更相似的。

huangzi

yuechenghua 发表于 2025-12-31 09:46
oplsaa里面F的种类有好几种，lj参数之间有差异。
opls_786 是PF6-里的F，sigama和epsilon分别是3.11810e-0 ...

老师，这是我tfsi-离子的结构，f的参数我倾向于使用文献中的，o的目前我不太确定，老师可以给一下建议吗？

yuechenghua

huangzi 发表于 2025-12-31 23:18
老师，这是我tfsi-离子的结构，f的参数我倾向于使用文献中的，o的目前我不太确定，老师可以给一下建议吗 ...

这个结构N带负电的话多多少少会对S有极化影响，使得LJ势能曲线有一定变形，导致参数有变化。而在oplsaa原力场里没有一模一样的化学环境，最接近这个结构的就是砜类化合物里的S=O（opls_494），但参数和你给的文献不同。所以，如果文献中给出的化合物和你的结构相似你就用文献里的。
这可能的考量是用其他原子种类会缺少对应的二面角参数（比如F-C-S-O），如果你O的原子种类选的很偏，oplsaa中可能没有这个二面角，所以为了避免额外拟合就选了一个有参数的类型去做模拟了。

slxc920113

CL&P或者OPLS2009-IL都是为离子液体专门优化的力场参数，本身就和你文献里面的OPLS-AA不是同一种力场，你不能要求两种不同的力场给出的LJ是一样的吧。