[gmx-users] a question in energy minimization step

[Banafsheh Mehrazma]

Hi everyone;
I am about to minimize a system in 3 steps;
first, I want my molecule to be minimized, so I used ''define =
-DPOSRES_WATER " in my *.mdp file., Then after using g_energy command to
make potential.xvg, it gave me a nice diagram. Next I wanted to perform
minimization with my protein to be constained and then a
whole_system_minimization without any constraints. However, the next
diagrams are not appropriate; one is linear, and another is increasing,
instead of decreasing!
I should mention that they all converged before 1340, 201, 47 steps,
respectively, and I set emtol option to be 100 and nsteps = 500000.
Do you know what is the problem?

yours sincerely
B.Mehrazma
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120125/47909595/attachment.html>