Re: SHAKE troubles

From: Lepsa <>
Date: Wed 29 May 2002 16:02:15 +0200

Dear Arvid,
this happened to me as well as my system was equilibrated badly. You might
consider slower warming using the &wt option (which helped me) as is
discussed in the DNA tutorial or run short (several ps) MD of the watbox
alone prior to the whole-system MD.

Hope this helps.


Martin Lepsik, PhD student            Phone:  +420/2/20183-540, Fax:
Dept. of Theor Chem & Center for Complex molecular Systems and Biomolecules
Institute of Organic Chemistry and Biochemistry (IOCB)
Flemingovo nam 2,
Czech Academy of Sciences,
166 10, Prague 6,
Czech Rep.
----- Original Message -----
From: "Arvid Soederhaell" <>
To: "amber" <>
Sent: Wednesday, May 29, 2002 12:07 PM
Subject: SHAKE troubles
> Dear all,
> I have some problems with SHAKE that I do not understand. I get the error
> message
>      NITER, NIT, LL, I AND J ARE :  0 0 1963 3195 3196
> no matter what i do???!!! The error message comes from the shake.f
> subprogram. The message implyes that i take a too long timestep, but even
> a timestep of 0.5 fs is too long! (Normally one should be able to use a 2
> fs timestep.) The simulation is stable if i turn off SHAKE and use a
> 0.2 fs timestep. The temperature in the simulation never exceeds 350 K.
> Any suggestions to what could be the problem?
> (The system is a DPPC lipid membrane with TIP3 water, and I use amber6)
> Arvid
Received on Wed May 29 2002 - 07:02:15 PDT
Custom Search