24小时热门版块排行榜    

查看: 1890  |  回复: 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)。
请问这样理解对吗??
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by zx04 at 2016-10-31 23:36:26
你用的是那个立场文件?

使用的是AmberTools+ACPYPE+Gaussian创建的小分子GAFF力场
3楼2016-11-01 13:45:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

zx04

木虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
jianjunjiang: 金币+5, 谢谢您的答复 2016-11-01 13:46:08
你用的是那个立场文件?
2楼2016-10-31 23:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zx04

木虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianjunjiang

铁杆木虫 (正式写手)

引用回帖:
4楼: Originally posted by zx04 at 2016-11-01 22:20:51
看看这里
http://ambermd.org/antechamber/gaff.html#function

GAFF力场里面二面角的公式应该是
sum(Vn/2  )
Vn are force constants; n is multiplicity and gamma is phase angle for torsional angle para ...

谢谢!
5楼2016-11-01 22:57:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +5 @taotao 2026-03-19 5/250 2026-03-24 21:07 by greychen00
[考研] 求调剂323材料与化工 +3 1124361 2026-03-24 3/150 2026-03-24 20:24 by peike
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 0805 316求调剂 +4 大雪深藏 2026-03-18 4/200 2026-03-24 17:34 by dick_runner
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
信息提示
请填处理意见