LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Angular velocities in principal axes frame
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Angular velocities in principal axes frame

From: Maurice de Koning <dekoning@...527...>
Date: Sun, 15 Apr 2018 17:04:41 -0300

Just to clarify, what I’m trying to do is to express a given “global” vector (i.e., the COM angular velocity in the inertial frame of the simulation box described by its fixed x, y and z directions)
in terms of the system of axes formed by the principal axes of the molecule. In other words, I’m just looking for a change of basis: going from the fixed simulation-box basis (centered at the COM) 
to that formed by the principal axes of the molecule.  

On 15 Apr 2018, at 14:39, Giacomo Fiorin <giacomo.fiorin@...24...> wrote:

Can you clarify what you mean by "principal axes frame"?

If you define a rotating frame of reference that is always aligned with the principal axes of a rigid molecule, the angular velocity of the molecule should be zero in that frame by definition.


On Sat, Apr 14, 2018 at 12:43 PM, Maurice de Koning <dekoning@...527...> wrote:
I’m interested in computing so-called rotational diffusion constants for rigid water molecules (say TIP4P, using shake to maintain the constraints).
One of the definitions for the rotational mean-squared displacement requires access to
molecular angular velocities expressed in the instantaneous principal axes frame. Searching
in the manual it seems that such angular velocities are available only with respect to the simulation
box frame directions. On the other hand, I’ve seen that for certain calculations (body particles), lammps
does compute the principal axes internally.

Would you have any suggestions as to how to compute angular velocities in the principal axes frame?



Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list

Giacomo Fiorin
Associate Professor of Research, Temple University, Philadelphia, PA
Contractor, National Institutes of Health, Bethesda, MD