Re: [AMBER] rotational diffusion error

From: David A Case <>
Date: Wed, 2 Oct 2013 07:34:36 -0400

On Wed, Oct 02, 2013, Vijay Manickam Achari wrote:
> rotdif nvecs 1000 ref [avg] :DDM.C15,C16 ncorr 90 ti 0.0 tf 0.80 dt
> 0.002 deffout deffs.dat itmax 500 tol 0.000001 d0 0.03 order 2 outfile
> rotdif.out

Do you really only have two atoms in your mask (":DDM.C15,C16")? I think this
might be a degenerate case where the rotation matrices are not well-defined.

Also, have you visualized the trajectory? Is the system undergoing rotational
tumbling? Be sure you can do a "standard" calculation (e.g. where the mask
covers the entire backbone, or the entire backbone that is in secondary
structure), and can understand the results. As I pointed out in an earlier
email, the rotdif command is not suited to analyze rotational behavior of just
part of system.


AMBER mailing list
Received on Wed Oct 02 2013 - 05:00:04 PDT
Custom Search