24小时热门版块排行榜    

查看: 3277  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见