CHARMM Development Project
Posted By: ca4930 Coor Inertia Axes - 09/19/18 09:00 PM
I want to compute the angles between the principal axis (i, j, k) and the atomic axis (x, y, z). So I'm using coor inertia to get the principal axis but I'm a bit confused as to what the secondary and tertiary axes describe. Could someone please elaborate?

Posted By: rmv Re: Coor Inertia Axes - 09/19/18 09:38 PM
Computing the inertial axes is a 3D eigenvalue problem, with the principal axis corresponding to the smallest eigenvalue, i.e. the axis with least inertia wrt. rotation around the axis.

The secondary and tertiary axes are orthogonal to the primary axis, with successively larger inertia (eigenvalues), as more mass is rotated.

One thing to note is that the direction of the principal axis is arbitrary, in terms of the sign of the vector components. For example, the principal inertial axis of a helix could point in either direction.
Posted By: ca4930 Re: Coor Inertia Axes - 09/20/18 06:41 PM
I see. So the components of the principal, secondary and tertiary axes are essentially the cosine of the angles of rotation.
Posted By: rmv Re: Coor Inertia Axes - 09/20/18 06:54 PM
The inertial axes are unit vectors, so yes the components correspond to cosines wrt. to the lab frame axes.
© CHARMM forums