[AMBER] "vlimit exceeded for step" for constant surface tension simulation

From: Wei Lin <w5lin.ucsd.edu>
Date: Thu, 20 Jun 2013 18:09:32 -0700

Dear All,

I am running a constant surface tension simulation for a lipid bilayer with
TIP3P water, after treating the system with minimization and heating.

After running 1735 steps, the simulation was terminated with error:

*vlimit exceeded for step 1735; vmax = 36.7881*

I went through amber mailing list to short out this but haven't get over
this error. Following are the output and input files for constant surface
tension simulation:

1). The last five steps of mdout:

 ------------------------------------------------------------------------------


 NSTEP = 1731 TIME(PS) = 122.462 TEMP(K) = 298.80 PRESS =
-91.2
 Etot = -97016.5662 EKtot = 27167.5502 EPtot =
-124184.1164
 BOND = 2038.0397 ANGLE = 8024.9141 DIHED =
 4573.5582
 1-4 NB = 2245.9637 1-4 EEL = -8137.7064 VDWAALS =
 8542.2872
 EELEC = -141471.1728 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 8707.8767 VIRIAL = 9568.4910 VOLUME =
 436886.7010
                                                    SURFTEN =
-88.8329
                                                    Density =
0.9749
 Ewald error estimate: 0.1524E-05
 ------------------------------------------------------------------------------


 NSTEP = 1732 TIME(PS) = 122.464 TEMP(K) = 299.44 PRESS =
-68.8
 Etot = -97011.4312 EKtot = 27225.3013 EPtot =
-124236.7326
 BOND = 2018.0758 ANGLE = 7988.6552 DIHED =
 4580.1512
 1-4 NB = 2246.0192 1-4 EEL = -8151.3022 VDWAALS =
 8552.7792
 EELEC = -141471.1110 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 8719.6867 VIRIAL = 9368.4374 VOLUME =
 436883.7151
                                                    SURFTEN =
-86.8008
                                                    Density =
0.9749
 Ewald error estimate: 0.4167E-05
 ------------------------------------------------------------------------------


 NSTEP = 1733 TIME(PS) = 122.466 TEMP(K) = 299.83 PRESS =
-47.6
 Etot = -97011.4352 EKtot = 27260.8468 EPtot =
-124272.2819
 BOND = 2020.0481 ANGLE = 7942.0973 DIHED =
 4588.1579
 1-4 NB = 2244.8143 1-4 EEL = -8170.2810 VDWAALS =
 8568.4146
 EELEC = -141465.5332 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 8718.3490 VIRIAL = 9167.1057 VOLUME =
 436881.6044
                                                    SURFTEN =
-74.1213
                                                    Density =
0.9749
 Ewald error estimate: 0.2476E-04
 ------------------------------------------------------------------------------


 NSTEP = 1734 TIME(PS) = 122.468 TEMP(K) = 299.92 PRESS =
-29.7
 Etot = -97015.7172 EKtot = 27269.0714 EPtot =
-124284.7886
 BOND = 2034.8483 ANGLE = 7918.7966 DIHED =
 4596.3835
 1-4 NB = 2241.9542 1-4 EEL = -8201.9951 VDWAALS =
 8582.1423
 EELEC = -141456.9185 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 8708.0900 VIRIAL = 8988.1628 VOLUME =
 436880.3199
                                                    SURFTEN =
-52.8906
                                                    Density =
0.9749
 Ewald error estimate: 0.1619E-04
 ------------------------------------------------------------------------------


 NSTEP = 1735 TIME(PS) = 122.470 TEMP(K) = 299.91 PRESS =
 -9.4
 Etot = -97014.3348 EKtot = 27267.9876 EPtot =
-124282.3225
 BOND = 2045.6125 ANGLE = 7931.1125 DIHED =
 4602.9317
 1-4 NB = 2239.3833 1-4 EEL = -8249.8078 VDWAALS =
 8594.3270
 EELEC = -141445.8816 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 8700.2075 VIRIAL = 8788.9502 VOLUME =
 436879.7324
                                                    SURFTEN =
-26.9472
                                                    Density =
0.9749
 Ewald error estimate: 0.1246E-04
 ------------------------------------------------------------------------------


2) The input file,

&cntrl
   imin = 0,ntx = 5,irest = 1,
   ntc = 2,ntf = 2,tol=0.000001,
   ntb = 2,ntp = 3,taup = 1.0,
   csurften = 3,gamma_ten = 10.0,ninterface = 2,
   tempi = 300.,temp0 = 300.,ntt = 3,gamma_ln = 1.,
   nstlim = 2500,dt = 0.002,ig=-1,cut = 9.,
   iwrap = 1,
   ntpr = 1,ntwr = 1,ntwx = 1,ntwe = 1,
/


Thanks,
Wei
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jun 20 2013 - 18:30:03 PDT
Custom Search