[gmx-users] diffusion, g_msd and trestart
chees at nus.edu.sg
Mon Jul 27 11:05:44 CEST 2009
I am trying to calculate diffusion coefficients for a system of (30)
glycine molecules in a box with (1000) water.
g_msd -f md.xtc -s md.tpr -mol -n molindex.ndx
I get a list of diffusion coefficients - one for each glycine molecule.
I then calculate the average and the stdev:
My problem is that the stdev is very great! More than +/- 75% of the
Why is that?
Also, what is the reason for the trestart=10 default value? Shouldn't it
be better to have trestart=1, if my xtc file is saved every 1 ps?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users