[gmx-users] msd and diffusion coefficient and getting unfolded trajectory

Mark Abraham Mark.Abraham at anu.edu.au
Fri May 9 10:10:10 CEST 2008


Ananya Debnath wrote:
>  Hi,
> 
>   I want to calculate MSD (mean square displacement) for DPPC-WATER
> bilayer system, without using g_msd tool in gromacs, but using any
> other code, I
> need a trajcetory which is unfolded.
> 
>>From mails in mailing list I understand that people use trjconv -pbc
> nojump, to remove pbc.
> So I used pbc = xyz in .mdp file, and then
> trjconv -f traj.trr -s topol.tpr -o trajout.trr -pbc nojump
> 
> Then if I calculate MSD (by my own code) from the trajout.trr , the
> result is same if I calculate from the original traj.trr file.
> 
> Seeing the result I am confused, because -pbc nojump sholud stop the
> molecule from jumping across the box, then how can the MSD be same if
> it is with pbc?

Maybe your trajectory does not have your relevant RMSD groups crossing 
the PBC, so trjconv is not doing anything relevant?

Mark



More information about the gromacs.org_gmx-users mailing list