Your system is probably unstable. With shake off have you looked carefully
at the structure over time (mdcrd)? What about the energies and the
temperature?
Try running 1,000 frames or so with ntwx set to 1. Then visualise it in
something like VMD, the problem should be obvious, you either have two atoms
close together or your system is unstable. I trust you have carried out
extensive minimisation on your system before starting MD.
All the best
Ross
_____
From: owner-amber.scripps.edu [mailto:owner-amber.scripps.edu] On Behalf Of
Huang, Hai
Sent: 19 November 2004 08:25
To: amber.scripps.edu
Subject: AMBER: SHAKE probelm
Hi, all,
I am a new user of AMBER. I generated a crosslink of oligonucleotide and
peptide with amber7. I try to simulate its structrure with molecular
dynamics, it works well with SHAKE off (ntc=1, ntf=1), however, when SHAKE
is turned on (ntc=2, ntf=2), the running crashs rapidly and error message
appears:
Coordinate resetting (SHAKE) cannot be accomplished,
deviation is too large
NITER, NIT, LL, I and J are : 0 11 7 12 13
Note: This is usually a symptom of some deeper
problem with the energetics of the system.
I have carefully checked the structure. LEAP doesn't report any error or
warning. Does anyone have any idea about this problem? Thank you.
Hai Huang
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Nov 19 2004 - 17:53:00 PST
