| 查看: 1852 | 回复: 4 | ||
jianjunjiang铁杆木虫 (正式写手)
|
[求助]
一个关于gromacs二面角参数的问题 已有1人参与
|
|
各位大侠： 本人现在利用gromacs做分子动力学，在立场文件中，本人有一些关于二面角的参数不太清楚，具体如下： 比如说： [ dihedrals ] ; propers ; for gromacs 4.5 or higher, using funct 9 ; i j k l func phase kd pn 1 2 3 4 9 180.00 15.16700 2 ; 请问后面的phase，kd和pn这三个参数各代表什么势函数当中的什么参数？翻了一下手册，发现手册上有一个二面角扭转势能的公式 V(φ)=k(1+cos(nφ-φs)) 那么上面的参数文件中 phase 是不是就是φs（180），pn就是n（2），而kd就是k（15.16700）。 请问这样理解对吗？？ |
回复此楼» 猜你喜欢
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有5人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
面上项目申报
已经有3人回复
-
酰胺脱乙酰基
已经有9人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
zx04
木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 2129.8
- 散金: 200
- 帖子: 10
- 在线: 577.6小时
- 虫号: 1398620
- 注册: 2011-09-11
- 专业: 理论和计算化学
2楼2016-10-31 23:36:26
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8778.2
- 散金: 60
- 红花: 1
- 帖子: 925
- 在线: 250.3小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2016-11-01 13:45:48
zx04
木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 2129.8
- 散金: 200
- 帖子: 10
- 在线: 577.6小时
- 虫号: 1398620
- 注册: 2011-09-11
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jianjunjiang: 金币+15, ★★★★★最佳答案 2016-11-01 22:56:53
jianjunjiang: 金币+15, ★★★★★最佳答案 2016-11-01 22:56:53
|
看看这里 http://ambermd.org/antechamber/gaff.html#function GAFF力场里面二面角的公式应该是 sum(Vn/2 [1 + cos(n*theta -gamma )] ) Vn are force constants; n is multiplicity and gamma is phase angle for torsional angle parameters. 所以，你的参数差不多，但公式应该检查一下。 |
4楼2016-11-01 22:20:51
jianjunjiang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8778.2
- 散金: 60
- 红花: 1
- 帖子: 925
- 在线: 250.3小时
- 虫号: 509147
- 注册: 2008-02-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2016-11-01 22:57:17