[gmx-users] Energy drift
mark.j.abraham at gmail.com
Wed Aug 20 16:34:11 CEST 2014
On Wed, Aug 20, 2014 at 3:55 PM, Xiaobin Cao <xcao at lsu.edu> wrote:
> Thanks for your response, Mark. Can you help me to check the parameters
> for my mdp file to see if it is correct? Can you send me a mdp file for
> NVE? I will test it today.
Just turn off T and P coupling!
> By the way, my GROMACS version is single precision? Does it matter?
Generally not, if everything else is correct. But you can try double
precision to see for yourself.
> >Given that their bond constant for N-N is ~4.5 times higher than their
> >constant for C-H, and both have comparable equilibrium lengths, and
> >unconstrained C-H bonds are normally treated with either 1fs or 0.5fs, and
> >that time step should be rather shorter than the period of bond
> >oscillation, I'd have serious reservations about the quality of that
> >and your results suggest reason for that. If you've got the parameters
> >right, then you should be able to suppress (post-equilibration) NVE drift
> >from integration inaccuracy by taking a time step that is small enough.
> >small is small enough?
> >Thank you
> Gromacs Users mailing list
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users