[MMTK] calculate RMSD of a collection of atoms from trajectory

Shirley Siu siu at mpi-sb.mpg.de
Tue Mar 15 18:49:18 CET 2005

Dear list,

I would like to calculate the RMSD of some atoms from the trajectory,
using the first frame as the reference. However, the result showed only
very little different and it is suspicious. I wonder if my way of
doing it is correct or not.

... load the trajectory...

collect=proteins[1]  # only calculate RMSD for chain 1

for step in self.trajectory[1:]:
    print rms

The RMSD ranges from 0 to less than 0.3 over all the frames.

Thanks in advance.


"Art is the lie that helps tell the truth..."     ^o^ Pablo Picasso

More information about the mmtk mailing list