[gmx-users] mdrun

leila karami karami.leila1 at gmail.com
Sat Nov 14 16:13:08 CET 2009


dear chris

I did EM :

Steepest Descents did not converge to Fmax < 2000 in 101 steps.
Potential Energy  = -3.2421097e+05
Maximum force     =  1.4766100e+04 on atom 1441
Norm of force     =  2.5146213e+02

I am sure that  constraints=none during EM.

how do I increase  the LINCS/SHAKE iterations to the appropriate number?
 this is my em.mdp file

cpp                 =  /usr/bin/cpp
define              =  -DFLEX_SPC
constraints         =  none
integrator          =  steep
nsteps              =  100
;
;       Energy minimizing stuff
;
emtol               =  2000
emstep              =  0.01
nstcomm             =  1
ns_type             =  grid
rlist               =  1
rcoulomb            =  1.0
rvdw                =  1.0
Tcoupl              =  no
Pcoupl              =  no
gen_vel             =  no
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20091114/75af3eee/attachment.html>


More information about the gromacs.org_gmx-users mailing list