[gmx-users] Normal Mode Analysis on Calphas only ??

Bert de Groot bgroot at gwdg.de
Tue May 17 18:44:54 CEST 2005


Berk is right, you need to diagonalise the full Hessian.
Not only because of the minimisation but also because
of the coupled interactions via other atoms in the Hessian,
which would be neglected by diagonalising a CA-only

Deriving CA-eigenvectors from the Hessian eigenvectors
can be done approximately by extracting CA components from
the full vectors, but if you require a fully orthogonal set
then the easiest procedure for now (awating the formula Berk
mentioned) is to generate an ensemble based on the normal
modes with g_nmens and do a CA PCA on that.


Bert de Groot, PhD

Max Planck Institute for Biophysical Chemistry
Computational biomolecular dynamics group
Am Fassberg 11
37077 Goettingen, Germany

tel: +49-551-2012308, fax: +49-551-2012302

email: bgroot at gwdg.de

More information about the gromacs.org_gmx-users mailing list