计算化学公社

标题: 求助：阴离子聚合物溶液模拟的力场、电荷及水模型的选择 [打印本页]

作者
Author: 陈欣羽 时间: 2026-3-19 11:36
标题: 求助：阴离子聚合物溶液模拟的力场、电荷及水模型的选择
本来想使用论坛sob老师推荐对聚合物的OPLS-AA＋1.2*CM5的组合，但尝试后发现因为我的聚合物长链带3个负电荷，电荷乘上1.2后电荷变成了-3.6；后继续在论坛上找到sob老师说这种需要用RESP电荷，所以我想对我模拟的阴离子聚合物溶液：阴离子聚合物和金属离子使用OPLS-AA力场，阴离子聚合物电荷使用RESP电荷，水模型采用TIP4P-ICE（需负温环境模拟）。因为没有找到使用完全相同文献，所以想来求助一下大家这种OPLS-AA + RESP + TIP4P-ICE的模拟组合是否可行，若发表论文是否会被审稿老师认可。感谢。

作者
Author: Loading0760 时间: 2026-3-19 13:34
这里引用sob老师的原话。
我的理解是1.2是一种经验的修正，CM5电荷本身与OPLS的相性会好很多，我处理类似情况的时候是直接采用CM5电荷，也就是不进行缩放。

前述J. Phys. Chem. B.文章中只是考虑了中性体系，没有考虑带电体系。带电体系用什么原子电荷结合OPLS-AA没有统一说法，一般情况下不能用1.2*CM5电荷，因为此时总电荷都不是整数了。笔者建议此时用Multiwfn算RESP2电荷，虽然与OPLS-AA兼容性没有充分测试，但原理上问题不大。
1.2*CM5电荷就是在Truhlar等人提出的原始的CM5电荷基础上乘上1.2，这相当于增大了原子电荷的数量级，等效地体现出了溶剂环境对溶质的极化作用。

作者
Author: 陈欣羽 时间: 2026-3-19 14:57

Loading0760 发表于 2026-3-19 13:34
这里引用sob老师的原话。
我的理解是1.2是一种经验的修正，CM5电荷本身与OPLS的相性会好很多，我处理类似 ...

感谢您的解答，我曾经有尝试过查找带电聚合物是否有用CM5电荷的文章，但是并没有找到较有说服力期刊的文献对这个观点进行依托，如果已发SCI文章作为模拟目标，CM5电荷是否能被审稿人承认呢？

作者
Author: sobereva 时间: 2026-3-19 16:33

sob老师推荐对聚合物的OPLS-AA＋1.2*CM5的组合

这个说法不完整。我绝对没有将这作为唯一推荐，本身我也挺不喜欢OPLS-AA力场。我很推荐用GAFF或GAFF2结合RESP/RESP2电荷算这种体系，无论是中性还是离子情况都一样处理

作者
Author: 陈欣羽 时间: 2026-3-19 17:31

sobereva 发表于 2026-3-19 16:33
这个说法不完整。我绝对没有将这作为唯一推荐，本身我也挺不喜欢OPLS-AA力场。我很推荐用GAFF或GAFF2结合 ...

非常感谢sob老师的指导，根据您的推荐并查阅相关资料，决定使用GAFF2+RESP+TIP4P-ICE的结合方式。
原因：1、GAFF2 相比 GAFF对于长链聚合物的模拟更真实平滑，会在零下环境中更自然的扭转。
2、RESP2 电荷强制要求在量子化学计算时加入隐式溶剂模型，会和我之后在体系加入TIP4P-ICE产生“双重极化”问题，聚合物不真实的亲水，导致模拟产生致命错误。而RESP电荷在真空中计算，会和TIP4P-ICE更好的结合。
再次感谢sob老师。

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