24小时热门版块排行榜    

查看: 3409  |  回复: 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的回帖

小清新王虫虫

新虫 (初入文坛)


小木虫: 金币+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的回帖
查看全部 3 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +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
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +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
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见