关于GAFF力场参数优化的一些问题

KiritsuguPapa · 发表于 Post on 2017-12-16 10:03:09

smutao 发表于 2017-12-6 00:05
你们老板好像来我们这找工作了

what他才跳槽到香港刚刚快满一年……

KiritsuguPapa · 发表于 Post on 2017-12-16 10:06:15

本帖最后由 KiritsuguPapa 于 2017-12-16 10:30 编辑

beefly 发表于 2017-12-6 07:25
你的键角力常数是如何定义的？力场程序给的公式是d^2E/dA^2，但是一些很老很经典的书上，还要把上面的公 ...

是用的能量对角度(rad)的二阶导。Amber力场文件里(parm99.dat, gaff.dat)也是同样的单位。
你说的那种做法，考虑的是 A-B-C键里A/C原子沿垂直A-B/B-C键方向的位移 (dL)，此时的力常数定义是能量对该位移的二阶导(d^2E/dL^2)。它和对角度的二阶导的换算差了一个r^2（因为r*d\theta = dL）。
另外这种做法因为有A/C两边同时要考虑，要求一个平均值。我见过Seminario方法里用的是调和平均值，大概就是这个吧。

smutao · 发表于 Post on 2018-4-26 21:27:29

有兴趣可以关注一下这个工作
doi:10.1021/acs.jctc.7b01171

KiritsuguPapa · 发表于 Post on 2018-5-26 04:46:53

smutao 发表于 2018-4-26 21:27
有兴趣可以关注一下这个工作
doi:10.1021/acs.jctc.7b01171

谢谢拜读大作

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

[Amber] 关于GAFF力场参数优化的一些问题

浏览过的版块