计算化学公社

标题: 求助关于Coulomb (SR)的手动计算和gromacs计算结果不一致 [打印本页]

作者
Author: Ymw2021 时间: 2021-7-29 14:14
标题: 求助关于Coulomb (SR)的手动计算和gromacs计算结果不一致

为了弄清楚gromacs中Coulomb的计算，我手动计算了一个乙醇分子的Coulomb-14和Coulomb (SR)值，并将其和gromacs的计算结果对比。计算结果如下图所示。
gromacs计算中，盒子大小为4nm左右。coulombtype = cut-off，rcoulomb =2。
Coulomb-14手动计算结果和gromacs一致,但是Coulomb (SR)的计算值不一致。请问这是什么原因造成的？
Coulomb-14手动计算方法为（f=0.5, if i, j is 1,4）,Coulomb (SR)的计算方法为（f=0, if i, j is 1,2 or 1,3; f=0.5, if i, j is 1,4；otherwise f=1）

作者
Author: Ymw2021 时间: 2021-7-29 16:15
解决了，原来是Coulomb（SR）不包含Coulomb-14，手算时要把i,j=1,4的能量去掉。

作者
Author: sobereva 时间: 2021-7-29 16:15
看不懂你的Coulomb (SR)的计算描述，也不知道算的是什么体系，没法简单判断

注意手算的时候应只考虑2nm以内的，并且注意charge group的影响

作者
Author: Ymw2021 时间: 2021-7-29 23:56

sobereva 发表于 2021-7-29 16:15
看不懂你的Coulomb (SR)的计算描述，也不知道算的是什么体系，没法简单判断

注意手算的时候应只考虑2nm ...

谢谢sob的回复，问题已经得到了解决。

作者
Author: wenshang956 时间: 2021-8-2 16:59
请问你这个Coulomb (SR)值是怎么计算的呀？没太看懂，是计算的单个乙醇分子内不同原子之间的Coulomb值吗？能不能简单说一下你在gmx中的选择方式

作者
Author: Ymw2021 时间: 2021-8-11 15:16

wenshang956 发表于 2021-8-2 16:59
请问你这个Coulomb (SR)值是怎么计算的呀？没太看懂，是计算的单个乙醇分子内不同原子之间的Coulomb值吗？ ...

不好意思sob老师，这两天没有登录论坛，没有及时回复您的信息。我将gmx计算的mdp文件附上，系统内只有一个乙醇分子，是计算的单个乙醇分子内不同原子之间的Coulomb值。

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