[gmx-users] Fwd: regarding dihedral restraints
Sadaf Rani
sadafrani6 at gmail.com
Mon Mar 2 18:44:03 CET 2020
Dear Gromacs users
I am restraining three dihedral angles between atoms of two
different residues. when I measure it using the following command:-
gmx angle -type dihedral -f md_prod.xtc -n dihedral.ndx -ov group.xvg -of
trans.xvg
I am getting one of my dihedrals graphs as mentioned in the link.
https://www.dropbox.com/sh/1ff45is9412dld4/AACDXjV0T47P01RgCJ1Baccca?dl=0
My dihedral angles have a multiplicity of 1 as mentioned below:-
[ dihedral_restraints ]
; ai aj ak al type phiA dphiA fcA phiB dphiB
fcB
72 74 31 33 1 146.72 0.0 0 146.72 0.0
41.840
75 74 31 27 1 173.10 0.0 0 173.10 0.0 41.840
75 74 31 25 1 158.21 0.0 0 158.21 0.0 41.840
However, in gromacs manual, I read this:-
*Counting transitions only works for dihedrals with multiplicity 3*
Am I doing something wrong? How should I measure dihedral restraint over
time?
Thanks.
Sadaf
More information about the gromacs.org_gmx-users
mailing list