AMBER Archive (2005)Subject: AMBER: Amber Poisson-Boltzmann MD :Equilibration Problem
From: praveena moldy (bc_praveena_at_yahoo.co.in)
Date: Sun Feb 27 2005 - 23:36:36 CST
Hello Amber Group,
I want to run Poisson-Botzmann MD in NVE ensemble for a small peptide Ace-Ala3-NMe with the length of all bonds fixed. I have to equilibrate it to 300K slowly within 100ps.The time step of the simulaiton which I kept is 1fs. In the first 65ps the molecule is made to heat gradually to 300K, and in the successive 25ps the temperature is restrained to 300K and no temperature scaling is adopted in the successive 10ps.
The input I have used is the one below,
&cntrl
imin=0,
ntpr=50,
ntwx=50,
ntb=0,
cut=12.0, igb=10,
nstlim=100000,dt=0.001,
ntt=0,ntc=3,ntf=3,
temp0=300.0, tempi=0.0,
/
&wt
type='TEMP0', istep1=0, istep2=65000,
value1=0.0, value2=300.0,
/
&wt
type='TEMP0', istep1=66000, istep2=90000,
value1=300.0, value2=300.0,
/
&wt
type='END',
/
hold peptide fixed
1.0
RES 1 42
END
/
While equilibrating the molecule, I got something like (after 0.15ps),
vlimit exceeded for step 157 ; vmax = 5540034186820002.
vlimit exceeded for step 158 ; vmax = 5.4060661364929577E+18
vlimit exceeded for step 159 ; vmax = 5.1700189427506561E+19
vlimit exceeded for step 160 ; vmax = 2.84653933347553812E+19
vlimit exceeded for step 161 ; vmax = 3.85341586976328188E+20
PB warning in pb_miccg(): CG maxitn exceeded!
PB warning in pb_miccg(): CG maxitn exceeded!
PB warning in pb_miccg(): CG maxitn exceeded!
PB warning in pb_miccg(): CG maxitn exceeded!..........
and this has been continued till 3.4ps and the equilibration stops with the message,
vlimit exceeded for step 3400 ; vmax = 47.43705100474962
vlimit exceeded for step 3401 ; vmax = 139.07727246991954
vlimit exceeded for step 3402 ; vmax = 33.706126147486799
vlimit exceeded for step 3403 ; vmax = 65.636635184900626
vlimit exceeded for step 3404 ; vmax = 31.068175282916044
vlimit exceeded for step 3405 ; vmax = 40.0656845829796
vlimit exceeded for step 3406 ; vmax = 92.181793757499122
vlimit exceeded for step 3407 ; vmax = 31.354599848532565
vlimit exceeded for step 3408 ; vmax = 33.413529233998702
vlimit exceeded for step 3409 ; vmax = 29.615107718328723
vlimit exceeded for step 3410 ; vmax = 25.571031269519956
vlimit exceeded for step 3411 ; vmax = 26.054697642189538
vlimit exceeded for step 3412 ; vmax = 22.508337242911296
Frac coord min, max: -9.84445917080887817E-4, 0.52496741432105076
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.
SANDER BOMB in subroutine Routine: map_coords (ew_force.f)
Atom out of bounds. If a restart has been written,
restarting should resolve the error
I do not know what went wrong. Is there any misatake in the input given.Could somebody help me regarding this.
Suggest me an idea to rectify this error.
Thank you,
Praveena.G
Yahoo! India Matrimony: Find your life partneronline.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|