为什么二面角拟合时MM能量大于QM能量

lts！ · 发表于 Post on 2025-9-9 15:59:39

本帖最后由 lts！于 2025-9-13 13:48 编辑

我首先对初始结构进行DFT优化和频率计算，并用multiwfn得到了RESP电荷，通过sobtop得到itp和top文件（GAFF2力场），然后柔性扫描如图所示的二面角，将得到的结构保存为pdb文件后再转成gro文件，最后用cat命令生成多帧gro文件，然后rerun计算各个结构的potential能量作为MM能量，然后都减去初始结构的能量，最后作图。
首先有一个问题，就是在10-19-21-26二面角在旋转时，51-19-21-26也在同步的旋转，这在计算MM能量时如何处理呢。

图1 扫描的二面角10-19-21-26

图2 同步旋转的二面角51-19-21-26

之后我在处理itp文件时，将两个二面角其中一个的力常数设为0，和同时将两个二面角力常数都设为0，分别计算了MM能量，图中所示为MM(1)和MM(2)。发现其能量均大于QM。这是由于什么原因呢，是由于柔性扫描时还有其他二面角都发生了变化，然后力常数也不精确，导致MM能量过大吗？

图3 二面角扫描曲线

combined.gro (133.08 KB, 下载次数 Times of downloads: 0)

L8BOgujiah.itp (69.16 KB, 下载次数 Times of downloads: 1)

L8BOgujiah.top (352 Bytes, 下载次数 Times of downloads: 0)

mm.mdp (595 Bytes, 下载次数 Times of downloads: 1)

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GROMACS] 为什么二面角拟合时MM能量大于QM能量

浏览过的版块