Re: AMBER: Molecular Frame

From: Byungchan Kim <>
Date: Sat, 3 Jan 2004 14:00:20 -0500 (EST)

On Fri, 2 Jan 2004, Bill Ross wrote:

> > I want to calculate order parameter (S2) of backbone N-H.
> > To calculate rotation correlation function of N-H bond vector,
> > I need to calculate the bond vector respect to molecular frame.
> > Is there any pre-set way or program to do this easily?
> If you want to measure the angle of the bond vector with some
> vector derived from the overall molecule, you will need to
> define the latter vector. One possibility might be the line
> drawn between centers of mass of two groups of atoms. Another
> possibility might be to use the vector N-{molecule center of mass}
> as the measurement vector. These two could be done with carnal
> and possibly ptraj. Or you could take the angle of N-H with
> respect to the principle axis of the molecule, which might
> require some programming to accomplish. I am not familiar
> with the mathematical definition of "molecular frame", so
> am just throwing out some ideas.
> In any case, it sounds like you are interested in per-snapshot
> 'internal' measurements, so superimposing snapshots may not be
> relevant except to perform an approximate measurement whereby
> the N-H vector of the positioned molecule could be used for
> further calculations.

Actually, my favorate is to use principal axis of a molecule as the
molecular frame. If the molecule is rigid, I would not hesitate to use
the principal axis. But if the molecule is flexible during the simulation,
principal axis would also change and that is I am wondering if I can use
principal axis.

> Bill Ross
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to
> To unsubscribe, send "unsubscribe amber" to

Thank You So Much
Byungchan Kim

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Jan 14 2004 - 15:53:13 PST
Custom Search