[gmx-users] Problem with NPT equilibration

[Bing Bing]

Dear all,
I've tried to simulate a membrane system with POPC and protein following
Justin's tutorial. The system was minimized with 2 stages. First using
restraint on everything accept water for 10000 SD steps:-
Steepest Descents converged to Fmax < 1000 in 1785 steps
Potential Energy  = -1.69457097061162e+06
Maximum force     =  9.15408957613777e+02 on atom 24360
Norm of force     =  7.53074616222725e+03
I've proceed with 2nd stage minimization without restraint,
Steepest Descents converged to Fmax < 1000 in 844 steps
Potential Energy  = -1.87072963657420e+06
Maximum force     =  9.14155621660386e+02 on atom 24360
Norm of force     =  7.84646890755063e+03
The energy curve seems ok. With this, i went on with NVT equilibration (
restraint on protein) for 100ps. Both energy and temperature plot look fine
also.
But it failed while running NPT 1ns(restraint on protein also), it stop with
Range checking error, whereby the ci value exceeded the ci cutoff. it is
suggested that the minimzation might not done properly. But the minimization
were converged in both stages. Is there anything i missed out here?

Please advice.Thanks in advance

regards
bing bing
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090904/5b7d9c8e/attachment.html>