[gmx-users] order parameters is more than one
Karim Mahnam
karimmahnam93 at gmail.com
Fri Dec 23 19:24:45 CET 2016
Hello GMX users
I want to calculate order parameters (S2) for NH bonds of backbone of one
protein. I did 20 ns MD simulation and performed these commands according
to Proof Van der Spoel comments (
https://extras.csc.fi/chem/courses/gmx2007/analysis/):
1-gmx trjconv -f md.trr -s md.tpr -fit rot+trans -o norot.trr -nice -20
Select 'Protein' group
2-gmx make_ndx -f md.tpr -o n1.ndx
> a N H & !r 1 *
##* is number of proline residues
>name 19 NH
>q
3-gmx rotacf -s md.tpr -d -n n1.ndx -P 2 -f norot.trr -noaver -o
rotacf.xvg -nice -20 -dt 400
Select 'NH' group
4- ./process_rotacf.pl rotacf.xvg
But autocorrelation numbers for each amino acid and order parameters for
each amino acid was more than one.
Would you please tell me what is my error and whow can I solve this problem.
I attached S2.xvg file to this mail.
Thank you very much
More information about the gromacs.org_gmx-users
mailing list