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
-- gmx-users mailing list [email protected] http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected]. Can't post? Read http://www.gromacs.org/mailing_lists/users.php

