计算化学公社

标题: 通过autoff得出的原子 sigma， epsilon 参数与文献不一样 [打印本页]

作者
Author: huangzi 时间: 2025-12-30 17:47
标题: 通过autoff得出的原子 sigma， epsilon 参数与文献不一样
本帖最后由 huangzi 于 2025-12-30 18:09 编辑

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

作者
Author: student0618 时间: 2025-12-30 21:10
文献是按什么优化这参数的？

作者
Author: huangzi 时间: 2025-12-30 21:29
本帖最后由 huangzi 于 2025-12-30 21:39 编辑

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

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

作者
Author: yuechenghua 时间: 2025-12-31 09:46
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，还是一样，找跟你结构环学环境更相似的。

作者
Author: huangzi 时间: 2025-12-31 23:18
本帖最后由 huangzi 于 2025-12-31 23:25 编辑

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

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

作者
Author: yuechenghua 时间: 2026-1-1 09:48

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

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)