[gmx-users] The RMSD problem
qing.zhu at ttu.edu
Thu Sep 1 16:33:58 CEST 2005
For my understanding, what you said happens all the time. Because our
simulation systems are very small and simulation time-scale is kind of
short, system fluctuations and statistical errors may be big; even you
run the simulations using the same initial conditions but only different
random seeds, the simulation results could be significantly different.
So my suggestion is run replicate simulation (5 or 6 might be enough),
average the results, do statistical analysis, and find the real
statistical average of your RMSD.
Tanping Li wrote:
>Another problem bothering me for the
>protein(myoglobin) is the RMSD. Because it is a very
>compact protein made of 7 helices, so I expect the
>RMSD should be stable and low. I follow the steps EM,
>PR MD and full MD. But I still find out the RMSD is
>very sensitive. It converges to 0.15nm in one of the
>gound state I made, but if the initial condition is a
>little bit different, it converges to 0.25 sometimes.
>Could I believe that 0.15nm is the correct one and
>throw away the 0.25nm one? Where does this come from,
>force field or something else? Thank you very much to
>share your experience and view.
>gmx-users mailing list
>gmx-users at gromacs.org
>Please don't post (un)subscribe requests to the list. Use the
>www interface or send it to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users