[gmx-users] Problem in energy minimization for a empty membrane system
Bing Bing
jarbing09 at gmail.com
Thu Jun 25 11:09:36 CEST 2009
Hi! Thanks a alot for all the suggestions i've got from Justin and Mark.
I found out that i actually deleted the box information from the starting
structure PDB. Once i put it back, It worked. However,
it stopped with this message.
Stepsize too small, or no change in energy.
Converged to machine precision,
but not to the requested precision Fmax < 1000
You might need to increase your constraint accuracy, or turn
off constraints alltogether (set constraints = none in mdp file)
writing lowest energy coordinates.
Back Off! I just backed up em.trr to ./#em.trr.19#
Back Off! I just backed up em.gro to ./#em.gro.11#
Steepest Descents converged to machine precision in 76 steps,
but did not reach the requested Fmax < 1000.
Potential Energy = -2.24796600790506e+05
Maximum force = 1.14539522791599e+03 on atom 6619
Norm of force = 1.62353061702859e+04
Is it mean that the minimization is done and is already converged?
Here the result for the final Step 75.
Energies (kJ/mol)
Bond Angle Proper Dih. Ryckaert-Bell. Improper Dih.
1.23491e+03 9.85552e+03 4.48086e+03 4.81663e+03 1.33894e+03
LJ-14 Coulomb-14 LJ (SR) Coulomb (SR) Coul. recip.
4.09985e+03 8.28955e+03 -1.68827e+04 -1.04035e+05 -1.38691e+05
Position Rest. Potential Kinetic En. Total Energy Temperature
6.95754e+02 -2.24797e+05 0.00000e+00 -2.24797e+05 0.00000e+00
Pressure (bar)
0.00000e+00
The position restraint is working fine as well....
Thank you so much.
Regards,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090625/b0a363d3/attachment.html>
More information about the gromacs.org_gmx-users
mailing list