计算化学公社

标题: 求助：amber里的LJ参数与gromacs中的转化 [打印本页]

作者
Author: 1900 时间: 2022-3-15 18:28
标题: 求助：amber里的LJ参数与gromacs中的转化
如题，如何将amber力场中的LJ势转化为gromacs中top文件中的LJ势。我查看了这两个手册，但是对于σ是如何转换的我还是没看懂。和acpype给出的结果差一个2倍关系。

作者
Author: sobereva 时间: 2022-3-16 03:19
Sobtop（http://sobereva.com/soft/Sobtop）里的子程序，用来把AMBER力场文件里读的epsilon和r0转化为GROMACS里用AMBER/GAFF力场时用的epsilon和sigma参数，注释写得很清楚

!--------- Convert LJ parameters (epsilon in kcal/mol, r0 in Angstrom) to GROMACS convention (epsilon in kJ/mol, sigma in nm)
subroutine LJ_convert(eps,r0,gmx_eps,gmx_sig)
real*8 eps,r0,gmx_eps,gmx_sig
gmx_eps=eps*4.184D0 !kcal/mol -> kJ/mol
gmx_sig=r0/2**(1D0/6D0)/10 !r0 -> sigma (Angstrom), then to nm
end subroutine

作者
Author: 1900 时间: 2022-3-25 11:16
谢谢sob老师

作者
Author: angervlf 时间: 2023-10-11 15:42
兄弟，我用了sob老师的公示，还是存在一个两倍的关系，我还是不太懂sigma的转换关系，你懂了嘛？为什么会存在这个两倍的关系呢？

作者
Author: youknowdcf 时间: 2024-10-28 16:43

angervlf 发表于 2023-10-11 15:42
兄弟，我用了sob老师的公示，还是存在一个两倍的关系，我还是不太懂sigma的转换关系，你懂了嘛？为什么会存 ...

要乘2，https://blog.sciencenet.cn/home. ... =blog&id=718384

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