I am running a dynamics calculations using amber7 on a protein neutrilized
by one ion. I have minimized the system with DRMS 0.00001 with distance
dependent dielectric. I have equilibrated the system for 500 ps and
raising the temperature (from 0K to 300K) in smaller steps of 30K/50 ps.
Visual inspection of protein at the end of equilibration shows that IB has
moved a bit far away from the protein (nearest atom of prtein to IB is 64
A distance). While subjecting the system to production run, following
error occured.
-----------------
Frac coord min, max: -8.20798303E-07 0.879446721
The system has extended beyond
the extent of the virtual box.
Restarting sander will recalculate
a new virtual box with 30 Angstroms
extra on each side, if there is a
restart file for this configuration.
EWALD BOMB in subroutine Routine: map_coords (ew_force.f)
Atom out of bounds. If a restart has been written,
restarting should resolve the error
-------------
The above error also appears to be related to that only whereby the IB is
moving further away from the protein and resulting in atoms out of bounds.
The protein is intact. I am sure that restarting it is not going to solve
the problem as IB may move further away as simulation progresses.
What are the possible reasons for IB drifting away? How to overcome this
problem?
Latest input of equilibration step (It was performed in steps):
&cntrl
imin=0, irest=1, ntx=5, ntb=0, cut=50,
ntr=0, tempi=300, temp0=300, ntt=1,
tautp=1.0, dielc=4,
ntpr= 500, ntwx=1000, ntwr=1000,
nstlim=50000,
&end
&ewald
eedmeth=5,
&end
Input for production run:
&cntrl
imin=0, irest=1, ntx=5, ntb=0, cut=50,
ntr=0, tempi=300, temp0=300, ntt=1,
tautp=0.5, dielc=4
ntpr= 5000, ntwx=10000, ntwr=10000,
nstlim=9500000,
&end
&ewald
eedmeth=5,
&end
----------------------------
-Balvinder
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Oct 07 2004 - 12:53:00 PDT
