24小时热门版块排行榜    

查看: 3410  |  回复: 2

zhuimeng08

银虫 (正式写手)

[交流] gromacs 计算二面角 已有2人参与

我用的命令;g_angle -f md.xtc  -n angle.ndx -oc dihan.xvg
报的错误:Fatal error:
number of index elements not multiple of 3, these can not be angle triplets
求帮助!!!!!

但是说明:g_angle computes the angle distribution for a number of angles or dihedrals.

With option -ov, you can plot the average angle of a group of angles as a
function of time. With the -all option, the first graph is the average and
the rest are the individual angles.

With the -of option, g_angle also calculates the fraction of trans dihedrals
(only for dihedrals) as function of time, but this is probably only fun for a
select few.

With option -oc, a dihedral correlation function is calculated.

It should be noted that the index file must contain atom triplets for angles
or atom quadruplets for dihedrals. If this is not the case, the program will
crash.

With option -or, a trajectory file is dumped containing cos and sin of
selected dihedral angles, which subsequently can be used as input for a
principal components analysis using g_covar.

Option -ot plots when transitions occur between dihedral rotamers of
multiplicity 3 and -oh records a histogram of the times between such
transitions, assuming the input trajectory frames are equally spaced in time.

应该没有问题呀,奇怪了就是算不了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

Cong2014

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
估计是你的index文件不对,计算angle的话,index格式如下:
#####
[ angle ]
a1 b1 c1
a2 b2 c2
........
........
an-1 bn-1 cn-2
an bn cn
#####
g_angle只能计算上述N个angle的平均值随时间的变化,或者其distribution.
2楼2014-11-07 10:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小清新王虫虫

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
gmx angle -f md.xtc -n index.ndx -type dihedral -ov dihedral.xvg ,之前和你情况一样,用该命令已解决问题
3楼2017-11-09 11:02:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuimeng08 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见