求助：拟合二面角参数，怎么用Origin程序拟合二面角扭转项参数？

灵芝5

你好！

我按照Sob老师的帖子拟合一个二面角的参数，帖子步骤如下：
1 用量子化学程序扫描二面角扭转势能曲线。参考《详谈使用Gaussian做势能面扫描》（http://sobereva.com/474）
2 对当前体系照常产生拓扑文件，但是把那个二面角项去除，然后算(1)的扫描曲线上每个点的MM能量
3 将(1)的能量减去(2)的能量，剩下的部分就是纯粹应当由二面角扭转项描述的势能曲线。将这个作为被拟合对象，用Origin程序或自己写的程序拟合扭转项参数即可。


请问第3步的“用Origin程序拟合扭转项参数”，怎么做操作？
这样得到的二面角参数适合CHARMM力场吗？

sobereva

自己在origin建立个新函数，拟合时选择相应函数。不清楚的话google搜或者看帮助
对于各种力场都是共通的。搞清楚相应力场的函数形式、单位这些细节

灵芝5

谢谢卢老师答复！