Hello, Amber users
I'm studying a protein using MD and would like to calculate correlation functions with the iRED method in order to compare order parameters from NMR. My protein is 138 residues long and contains 10 prolines, and therefore I have 128 N-H vectors. I'm using Amber10. My input file looks like this:
trajin ../mdcrd5.gz
trajin ../mdcrd6.gz
vector v2 :2.N ired :2.H
vector v3 :3.N ired :3.H
..
vector v138 :138.N ired :138.H
matrix ired name matired order 2
analyze matrix matired vecs 128 out ired.vec
vector v2 :2.N corrired :2.H order 2 modes ired.vec beg 1 end 128 npair 1
vector v3 :3.N corrired :3.H order 2 modes ired.vec beg 1 end 128 npair 2
..
vector v138 :138.N corrired :138.H order 2 modes ired.vec beg 1 end 128 npair 128
analyze timecorr vec1 v2 tstep 10.0 tcorr 10000 out Ired/v2.out
analyze timecorr vec1 v3 tstep 10.0 tcorr 10000 out Ired/v3.out
..
analyze timecorr vec1 v138 tstep 10.0 tcorr 10000 out Ired/v138.out
(I've have also tried to break it up in two ptraj scripts, since the manual is a little bit vague if this is neccessary.) The problem is that all the correlation functions calculated, v2.out, v3.out, .. v138.out is the same, i.e. all the output files contains the same numbers. What am I doing wrong? I can hardly believe that all the correlation functions should be identical.
And when I'm at writing - what is the best way to obtain the order parameters from the correlation functions?
Best regards, Samuel
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Jun 22 2008 - 06:07:30 PDT
