[gmx-users] g_msd and MSD analysis
dommert at icp.uni-stuttgart.de
Sun Jun 10 10:55:53 CEST 2012
On Sat, 2012-06-09 at 15:33 -0700, Mr Bernard Ramos wrote:
> Hi everyone!
> I have a 20 ns simulation (2fs timestep but coordinates saved every
> 0.2 ps) and I was able to calculate the mean-square-displacement of
> the oxygen atoms of my water solvent. The entire MSD plot looks very
> linear to me. I need to fit this set of MSD data to obtain the
> diffusion coefficient. However, I am not sure if I will have to use
> the entire MSD points or may be just to consider a few picoseconds at
> the latter part of the simulation time. In papers I have read, they
> use up to a 100 ps, although no mention were made as to where those
> set of points came from. Please help.
have you looked at a linear graph or on a log-log representation. With
the latter representation, it is much easier to identify the different
regimes of diffusion (ballistic, subdiffusive, and linear). This should
allow you to chose an appropriate fit range. Moreover if you fit a
straight line to the log-log representation you will immediately see,
where the MSD curve deviates from the linear behaviour.
> This post seems to suggest a nanosecond window for the analysis of the
> MSD curve
> Although, this other one suggests using the latter 40 ps of his MSD
> Your help is greatly appreciated.
> gmx-users mailing list gmx-users at gromacs.org
> Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Dipl. - Phys.
Institute for Computational Physics
EMail: dommert at icp.uni-stuttgart.de
Tel.: +49 - (0)711 - 68563613
Fax.: +49 - (0)711 - 68563658
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 198 bytes
Desc: This is a digitally signed message part
More information about the gromacs.org_gmx-users