[gmx-users] Water molecule can not be settled - mdrun error
jalemkul at vt.edu
Mon Sep 24 23:05:44 CEST 2012
On 9/24/12 4:57 PM, Lara Bunte wrote:
> How to make the energy minimization in double precision? Do this mean to use grompp twice?
No, you would have to recompile Gromacs in a different manner:
I don't know if we have to go that far yet. Your forces look very bad, such
that I don't know that a change in precision would make a difference. Before
taking the time to recompile Gromacs, there are other things to look at with the
> I changed in my energy minimization I changed emtol from 1000 to 100 and got:
> Steepest Descents:
> Tolerance (Fmax) = 1.00000e+02
> Number of steps = 5000
> Stepsize too small, or no change in energy.
> Converged to machine precision,
> but not to the requested precision Fmax < 100
> Steepest Descents converged to machine precision in 355 steps,
> but did not reach the requested Fmax < 100.
> Potential Energy = -1.3305919e+05
> Maximum force = 3.6814004e+03 on atom 19
> Norm of force = 6.2041531e+01
> Is this good or bad?
Bad. You're asking mdrun to converge to an Fmax < 100 and it reaches Fmax =
3681.4 instead. Most normal simulations need an Fmax < 1000 or so to be
feasible, better if you can converge to a lower force.
> I also changed my .mdp file for the equilibration of the water around the flavin with the right cut-offs. Now it looks like this:
> define = -DPOSRES
> integrator = md
> dt = 0.002
> emtol = 1000.0
> emstep = 0.01
Note that with the md integrator, you are not doing energy minimization, so
emtol and emstep are ignored. Not a problem per se, just pointing out
potentially confusing settings.
> nsteps = 5000
> nstlist = 1
You need not set nstlist = 1 for running MD; this is only necessary during
minimization. A value of 5 or 10 would be more appropriate.
> rlist = 1.2
> rlistlong = 1.4
> rcoulomb = 1.2
> coulombtype = pme
> vdw-type = switch
> rvdw_switch = 0.8
> rvdw = 1.2
> nstxout = 100
> nstvout = 100
> nstenergy = 100
> nstlog = 100
> tcoupl = V-rescale
> tc-grps = ISO SOL
> tau_t = 0.1 0.1
> ref_t = 300 300
> pcoupl = no
> Is this okay now? I have a further question about that. As far as I understand I want to equilibrate the water around the lumiflavin which is again a kind of minimization. Why do I have to use the integrator md and not steep like in the minimization before? I use integrator md and also DPOSRES because my tutorial says I have to use this by I asked myself why not steep?
Energy minimization and molecular dynamics are different processes. The "steep"
(steepest descents) is a minimization (non-dynamical) algorithm. It is not MD.
> Sadly with this settings I am again at the problem from the beginning of our discussion.
> step 266: Water molecule starting at atom 4416 can not be settled.
> Check for bad contacts and/or reduce the timestep if appropriate.
> Wrote pdb files with previous and current coordinates
This is a result of inadequate energy minimization. Your starting forces are
too high and lead to instability. There is no point in running any MD if your
EM does not converge adequately.
> I hope that understanding more about this minimization stuff will help me to fix this. I am interested what could be now a problem in the .mdp file.
The .mdp file looks reasonable (the nstlist setting should not be the cause of
the problem), indicating to me that your topology is unsound. I've asked about
the source of the parameters several times, but maybe I should ask it again -
from where did you get the parameters? Did you derive them yourself? Does the
simulation run in vacuo (i.e. no water, infinite cutoffs, and no PBC)?
Justin A. Lemkul, Ph.D.
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users