24小时热门版块排行榜    

查看: 3407  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 13:44 by winsaint
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 5/250 2026-08-28 12:41 by Tide man
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +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 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见