计算化学公社

标题: 求助：关于在Gromacs中使用scaled charge [打印本页]

作者
Author: 北京的大猫 时间: 2020-7-16 09:44
标题: 求助：关于在Gromacs中使用scaled charge
大家好！

想问一下如果在gromacs中使用OPLSAA力场时，想对电荷进行0.8倍的缩放，是不是将每个原子的电荷*0.8后放入相关文件（比如rtp或者itp文件）就可以了？

谢谢各位老师！

作者
Author: sobereva 时间: 2020-7-16 09:45
是

作者
Author: 北京的大猫 时间: 2020-7-16 10:47

sobereva 发表于 2020-7-16 09:45
是

谢谢sob老师！

作者
Author: bobosiji 时间: 2020-7-17 20:37
为啥要缩放呢

作者
Author: 北京的大猫 时间: 2020-7-18 11:59
本帖最后由北京的大猫于 2020-7-18 12:02 编辑

bobosiji 发表于 2020-7-17 20:37
为啥要缩放呢

应该是能更好地复现实验性质吧,用缩放原子电荷的方法来解释极化效应，以模拟极化和电荷转移引起的平均电荷屏蔽.

作者
Author: bobosiji 时间: 2020-7-18 22:25

北京的大猫发表于 2020-7-18 11:59
应该是能更好地复现实验性质吧,用缩放原子电荷的方法来解释极化效应，以模拟极化和电荷转移引起的平均电 ...

如果遇到较真的审稿人问：这么做（缩放电荷）的理论依据？缩放电荷模拟某些物理量的结果能改善，其他的物理量呢？最关键的是：分子力场参数，比如键角、二面角参数和电荷都是耦合的，如果缩放电荷，其他力场参数也应该要相应调整吧？
如何回答呢？多谢提点~

作者
Author: 北京的大猫 时间: 2020-7-19 16:19

bobosiji 发表于 2020-7-18 22:25
如果遇到较真的审稿人问：这么做（缩放电荷）的理论依据？缩放电荷模拟某些物理量的结果能改善，其他的物 ...

我觉得您说得有道理,使用原始的OPLSAA力场中电荷应该保持原值吧,但是我的用的文献中开发力场时电荷是缩放的,它的开发只是opls力场形式的.不知我说的对不对.另外键角的参数和电荷不是耦合的吧?它是通过拟合振动频率得到的.不知我说得对不?新手一枚

作者
Author: bobosiji 时间: 2020-7-19 17:04

北京的大猫发表于 2020-7-19 16:19
我觉得您说得有道理,使用原始的OPLSAA力场中电荷应该保持原值吧,但是我的用的文献中开发力场时电荷是缩放 ...

二面角参数肯定和电荷耦合，除非不考虑1-4作用

作者
Author: sobereva 时间: 2020-7-20 02:20

bobosiji 发表于 2020-7-18 22:25
如果遇到较真的审稿人问：这么做（缩放电荷）的理论依据？缩放电荷模拟某些物理量的结果能改善，其他的物 ...

这种做法都是之前有人测试证明合理的
可以看看比如0.8 scaled OPLS-2009IL离子液体力场这篇（J. Chem. Theory Comput. 2017, 13, 6131）

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