计算化学公社

标题: OPLSAA力场，使用ligpargen生成itp时可以将itp中的电荷替换成高斯计算的吗 [打印本页]

作者
Author: wangzihan2026 时间: 2026-3-22 23:26
标题: OPLSAA力场，使用ligpargen生成itp时可以将itp中的电荷替换成高斯计算的吗
各位老师和同学好:
有个问题想请教大家我在查阅文献时发现我的领域大多数是使用oplsaa力场，一些人在补充资料里面写电荷是通过高斯拟合的，请问这个高斯拟合出来的电荷可以直接把ligpagen自动生成的电荷替换吗

作者
Author: student0618 时间: 2026-3-22 23:51
计算适用于OPLS-AA力场做模拟的1.2*CM5原子电荷的懒人脚本
http://bbs.keinsci.com/thread-21462-1-1.html
(出处: 计算化学公社)

作者
Author: wangzihan2026 时间: 2026-3-23 00:02

student0618 发表于 2026-3-22 23:51
计算适用于OPLS-AA力场做模拟的1.2*CM5原子电荷的懒人脚本
http://bbs.keinsci.com/thread-21462-1-1.html ...

谢谢回复，另外我还是想请教上面我询问的问题，这个itp中的电荷可以被替换吗，我看了sobtop软件里面默认的有gaff和amber这些力场输出的itp文件，这些文件里面的电荷数据我可以粘贴到oplsaa的itp里面吗十分感谢

作者
Author: student0618 时间: 2026-3-23 00:08

wangzihan2026 发表于 2026-3-23 00:02
谢谢回复，另外我还是想请教上面我询问的问题，这个itp中的电荷可以被替换吗，我看了sobtop软件里面默认 ...

Carefully read the post http://bbs.keinsci.com/thread-21462-1-1.html or search related disucssion in this forum

根据J. Phys. Chem. B., 121, 3864 (2017)中的测试，基于OPLS-AA力场的模拟很适合结合1.2*CM5原子电荷。虽然用1.2*CM5的时候水合自由能的计算精度稍逊于文中测试的另一种原子电荷1.14*CM1A-LBCC（因为这种原子电荷本来就是针对水合自由能训练的校正参数），但在计算其它属性上，如蒸发焓、密度，用1.2*CM5时的误差都更小，而且经验性更低、更有普适性。所以当大家用OPLS-AA力场研究新的小分子的时候，我比较推荐用1.2*CM5电荷。1.2*CM5电荷就是在Truhlar等人提出的原始的CM5电荷基础上乘上1.2，这相当于增大了原子电荷的数量级，等效地体现出了溶剂环境对溶质的极化作用。

作者
Author: sobereva 时间: 2026-3-23 06:58
自己把Multiwfn算出来的1.2*CM5电荷粘贴到itp的[atoms]原子电荷那一列就完了

作者
Author: wangzihan2026 时间: 2026-3-23 10:19

sobereva 发表于 2026-3-23 06:58
自己把Multiwfn算出来的1.2*CM5电荷粘贴到itp的[atoms]原子电荷那一列就完了

十分感谢回复

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