Hi;
I'm using amber 8 to equilibrate a 43286 atom system. After
heating the system to 300
K, 25ps (1fs timesteps, ntt=3) of NVT simulation followed (also at
300K). When I start
the pressure equillibration, I get the following error:
ASSERTion 'xtindex <= 10 .and. xtindex >= 0' failed in ew_direct.f at
line 943.
The input file I used is:
&cntrl
imin = 0,
irest=1,
ntx=5,
ntpr = 100,
ntwx = 100,
ntwe = 100,
iwrap=1,
ntwr = 1000,
nstlim =10000,
dt = 0.001,
ntb = 2,
cut = 11,
ntr = 0,
ntc = 2,
ntf = 2,
nmropt = 0,
ntt = 1,
temp0 = 300.0,
ntp=1,
pres0=1.0,
taup=10,
&end
Note I tried using taup=3 and 5 before setting this parameter to 10.
Has anyone observed
this behavior before? Some ofthe output is below:
| Flags: MPI
getting new box info from bottom of inpcrd
| peek_ewald_inpcrd: Box info found
|Largest sphere to fit in unit cell has radius = 39.009
| New format PARM file being parsed.
| Version = 1.000 Date = 06/11/09 Time = 01:52:16
NATOM = 43286 NTYPES = 17 NBONH = 41720 MBONA = 1593
NTHETH = 3642 MTHETA = 2148 NPHIH = 6709 MPHIA = 3819
NHPARM = 0 NPARM = 0 NNB = 70976 NRES = 13575
NBONA = 1593 NTHETA = 2148 NPHIA = 3819 NUMBND = 62
NUMANG = 125 NPTRA = 48 NATYP = 45 NPHB = 1
IFBOX = 1 NMXRS = 24 IFCAP = 0 NEXTRA = 0
NCOPY = 0
| Memory Use Allocated
| Real 3049515
| Hollerith 273293
| Integer 2260012
| Max Pairs 330206
| Max Rstack 1227492
| Max Istack 216430
| Total 45445 kbytes
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals
......
4. RESULTS
--------------------------------------------------------------------------------
---------------------------------------------------
APPROXIMATING switch and d/dx switch using CUBIC SPLINE INTERPOLATION
using 5000.0 points per unit in tabled values
TESTING RELATIVE ERROR over r ranging from 0.0 to cutoff
| CHECK switch(x): max rel err = 0.2738E-14 at 2.422500
| CHECK d/dx switch(x): max rel err = 0.8168E-11 at 2.732560
---------------------------------------------------
| Local SIZE OF NONBOND LIST = 192434
| TOTAL SIZE OF NONBOND LIST = 18371501
NSTEP = 100 TIME(PS) = 85.100 TEMP(K) = 300.94 PRESS
= -255.0
Etot = -106228.7798 EKtot = 26353.5248 EPtot =
-132582.3046
BOND = 551.9689 ANGLE = 1656.7456 DIHED =
1400.9006
1-4 NB = 716.6664 1-4 EEL = 9639.4064 VDWAALS =
17214.2505
EELEC = -163762.2430 EHBOND = 0.0000 RESTRAINT =
0.0000
EKCMT = 11944.6498 VIRIAL = 14825.1260 VOLUME =
523250.0270
Density =
0.8354
Ewald error estimate: 0.6977E-05
------------------------------------------------------------------------------
NSTEP = 200 TIME(PS) = 85.200 TEMP(K) = 300.69 PRESS
= -323.9
Etot = -106235.4052 EKtot = 26331.7631 EPtot =
-132567.1683
BOND = 577.8501 ANGLE = 1681.5169 DIHED =
1398.3593
1-4 NB = 698.7472 1-4 EEL = 9611.1972 VDWAALS =
17231.8852
EELEC = -163766.7241 EHBOND = 0.0000 RESTRAINT =
0.0000
EKCMT = 12004.2147 VIRIAL = 15662.5754 VOLUME =
523186.6971
Density =
0.8355
Ewald error estimate: 0.6196E-05
------------------------------------------------------------------------------
NSTEP = 300 TIME(PS) = 85.300 TEMP(K) = 300.05 PRESS
= -359.8
Etot = -106244.9948 EKtot = 26275.5856 EPtot =
-132520.5804
BOND = 552.7989 ANGLE = 1668.2645 DIHED =
1420.1295
1-4 NB = 711.0346 1-4 EEL = 9663.4828 VDWAALS =
16953.2163
EELEC = -163489.5071 EHBOND = 0.0000 RESTRAINT =
0.0000
EKCMT = 12032.4122 VIRIAL = 16096.5153 VOLUME =
523113.0068
Ewald error estimate: 0.3080E-04
------------------------------------------------------------------------------
NSTEP = 400 TIME(PS) = 85.400 TEMP(K) = 299.88 PRESS
= -347.5
Etot = -106253.9265 EKtot = 26260.5138 EPtot =
-132514.4404
BOND = 594.7058 ANGLE = 1622.2264 DIHED =
1404.7160
1-4 NB = 706.1265 1-4 EEL = 9685.9846 VDWAALS =
16848.4413
EELEC = -163376.6410 EHBOND = 0.0000 RESTRAINT =
0.0000
EKCMT = 11923.5027 VIRIAL = 15847.3449 VOLUME =
523030.9328
Density =
0.8358
Ewald error estimate: 0.5722E-04
------------------------------------------------------------------------------
(the output ends at this point).
Any ideas will be greatly appreciated.
Thanks;
-Sergio
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 06 2009 - 09:44:48 PDT
