AMBER: using sander without SHAKE

From: Wayne Dawson <dawson.nih.go.jp>
Date: Mon, 22 Dec 2003 15:38:14 +0900

Dear list,

I am trying to study solvent effects on an electrostatic protein.
I have a structure I obtained, and I have done simulations on
this thing out to ns using SHAKE. But I want to look at some
simulations without SHAKE (i.e., setting ntc=1).

I simply cannot seem to run sander successfully with ntc=1. I've
already equilibrated this protein under NPT conditions out to 0.250 ns
using SHAKE, yet I still keep failing because of problems with
exceeding the vlimit. I've tried time-steps as short as 0.1
fs. Nothing seems to stabilize the calculation.

I can definitely calculate this same kind of simulation on Discover
(Amber 6) with comparable conditions (below). I want to check these
Discover results using TIP3P water and amber7. What am I doing wrong?

Here is the complete AMBER7 output and the stderr output is also
listed at the end. (A minor note: I don't know why this was run in
1952, maybe I've suffered from an AIX/ENIAC twin-paradox.)
------------------------------------------

          -------------------------------------------------------
          Amber 7 SANDER Scripps/UCSF 2002
          -------------------------------------------------------

| Run on Tue Apr 15 00:27:43 1952

  [-O]verwriting output

File Assignments:
| MDIN: md2-ntc1.in
| MDOUT: test_Boxmdx1.out
|INPCRD: test_Boxmdx0.rst
| PARM: test_Box.top
|RESTRT: test_Boxmdx1.rst
| REFC: refc
| MDVEL: mdvel
| MDEN: test_Boxmdx1.mden
| MDCRD: test_Boxmdx1.mdcrd
|MDINFO: mdinfo
|INPDIP: inpdip
|RSTDIP: rstdip


 Here is the input file:

MD simulation, 10 ps
&cntrl
   imin=0, irest=1, ntx=5,
   nstlim=1000,
   nmropt=1,
   ntwr=200,
   ntwe=100,
   ntwx=100,
   ntpr=10,
   dt=0.00025,
   ntt=1,
   ntp=1,
   pres0=1.015,
   cut=15.0,
   ntb=2,
   ntc=1,
   ntf=2,
   vlimit=20,
   tempi=310.0,
   temp0=310.0,
   tautp=1.0,
 /
&ewald
   eedmeth=1,
 /
&wt type='END'
 /
DISANG=RST

#

--------------------------------------------------------------------------------
   1. RESOURCE USE:
--------------------------------------------------------------------------------

| Flags: MEM_ALLOC MPI RLE ROWAT
 getting new box info from bottom of inpcrd
| peek_ewald_inpcrd: Box info found
| New format PARM file being parsed.
| Version = 1.000 Date = 11/25/03 Time = 15:27:54
 NATOM = 17461 NTYPES = 18 NBONH = 17032 MBONA = 436
 NTHETH = 913 MTHETA = 590 NPHIH = 1745 MPHIA = 1092
 NHPARM = 0 NPARM = 0 NNB = 26713 NRES = 5602
 NBONA = 436 NTHETA = 590 NPHIA = 1092 NUMBND = 42
 NUMANG = 86 NPTRA = 39 NATYP = 31 NPHB = 1
 IFBOX = 1 NMXRS = 24 IFCAP = 0 NEXTRA = 0


| Memory Use Allocated
| Real 1258602
| Hollerith 110370
| Integer 862094
| Max Pairs 6127563
| Max Rstack 2000000
| Max Istack 2000000
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals

     BOX TYPE: RECTILINEAR

--------------------------------------------------------------------------------
   2. CONTROL DATA FOR THE RUN
--------------------------------------------------------------------------------



General flags:
     imin = 0, nmropt = 1

Nature and format of input:
     ntx = 5, irest = 1, ntrx = 1

Nature and format of output:
     ntxo = 1, ntpr = 10, ntrx = 1, ntwr = 200
     iwrap = 0, ntwx = 100, ntwv = 0, ntwe = 100
     ioutfm = 0, ntwprt = 0, idecomp = 0, rbornstat= 0

Potential function:
     ntf = 2, ntb = 2, igb = 0, nsnb = 25
     ipol = 0, gbsa = 0
     dielc = 1.00000, cut = 15.00000, intdiel = 1.00000
     scnb = 2.00000, scee = 1.20000

Frozen or restrained atoms:
     ibelly = 0, ntr = 0

Molecular dynamics:
     nstlim = 1000, nscm = 1000, nrespa = 1
     t = 0.00000, dt = 0.00025, vlimit = 20.00000

Temperature regulation:
     ig = 71277, ntt = 1, vrand = 0
     temp0 = 310.00000, tempi = 310.00000, heat = 0.00000
     dtemp = 5.00000, tautp = 1.00000

Pressure regulation:
     ntp = 1, plevel = 1
     pres0 = 1.01500, comp = 44.60000, taup = 0.20000

NMR refinement options:
     iscale = 0, noeskp = 1, ipnlty = 1, mxcub = 1
     scalm = 100.00000, pencut = 0.10000, tausw = 0.10000

Ewald parameters:
     verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1
     vdwmeth = 1, eedmeth = 1, frc_int = 0, netfrc = 1
     Box X = 55.892 Box Y = 56.422 Box Z = 56.350
     Alpha = 90.000 Beta = 90.000 Gamma = 90.000
     NFFT1 = 60 NFFT2 = 60 NFFT3 = 60
     Cutoff= 15.000 Tol =0.100E-04
     Ewald Coefficient = 0.17871
     Interpolation order = 4
| PLEVEL = 1: runmd parallelization, no EKCMR

--------------------------------------------------------------------------------
   3. ATOMIC COORDINATES AND VELOCITIES
--------------------------------------------------------------------------------


 begin time read from input coords = 232.600 ps



           Begin reading energy term weight changes/NMR restraints
 WEIGHT CHANGES:
                         ** No weight changes given **

 RESTRAINTS:
 Requested file redirections:
  DISANG = RST
 Restraints will be read from file: RST
 Here are comments from the DISANG input file:
 #
 # 25 CYS SG 30 CYS SG 2.8 3.2

                       Number of restraints read = 1

                  Done reading weight changes/NMR restraints


| Reorder water Setup:
| ----- Reordering water will NOT be done
| NMROPT not equal 0
|
 Number of triangulated 3-point waters found: 5542
| Atom division among processors:
| 0 4366 8731 13096 17461
| Running AMBER/MPI version on 4 nodes


     Sum of charges from parm topology file = -0.00000005
     Forcing neutrality...

--------------------------------------------------------------------------------
   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.3388E-14 at 2.505100
| CHECK d/dx switch(x): max rel err = 0.7986E-11 at 2.679240
 ---------------------------------------------------
| Local SIZE OF NONBOND LIST = 4496835
| TOTAL SIZE OF NONBOND LIST = 17628514

 NSTEP = 10 TIME(PS) = 232.602 TEMP(K) = 298.10 PRESS = -2185.8
 Etot = -42809.5148 EKtot = 15514.4310 EPtot = -58323.9458
 BOND = 204.2093 ANGLE = 520.6311 DIHED = 530.3931
 1-4 NB = 181.2874 1-4 EEL = 2788.4967 VDWAALS = 7116.6741
 EELEC = -69666.2954 EHBOND = 0.0000 RESTRAINT = 0.6579
 EAMBER (non-restraint) = -58324.6037
 EKCMT = 5095.8433 VIRIAL = 13479.5162 VOLUME = 177644.1534
                                                Density = 0.9913
 Ewald error estimate: 0.6407E-05
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 10
 Energy (this step): Bond = 0.658 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.662 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.779 0.779 0.781 0.781 0.579 0.579 0.581 0.581
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================

 NSTEP = 20 TIME(PS) = 232.605 TEMP(K) = 627.70 PRESS =-10173.1
 Etot = -42878.7961 EKtot = 32668.2956 EPtot = -75547.0917
 BOND = 191.8645 ANGLE = 501.2709 DIHED = 529.0825
 1-4 NB = 181.3232 1-4 EEL = 2782.9497 VDWAALS = 7302.6889
 EELEC = -87036.9179 EHBOND = 0.0000 RESTRAINT = 0.6465
 EAMBER (non-restraint) = -75547.7382
 EKCMT = 5075.0875 VIRIAL = 43981.8937 VOLUME = 177130.4398
                                                Density = 0.9941
 Ewald error estimate: 0.8882E-05
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 20
 Energy (this step): Bond = 0.646 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.657 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.773 0.773 0.778 0.778 0.573 0.573 0.578 0.578
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================

 NSTEP = 30 TIME(PS) = 232.607 TEMP(K) = 1369.17 PRESS =-28768.3
 Etot = -43463.4825 EKtot = 71258.0206 EPtot = -114721.5031
 BOND = 166.7530 ANGLE = 443.3562 DIHED = 522.7411
 1-4 NB = 179.4392 1-4 EEL = 2778.0989 VDWAALS = 7972.7864
 EELEC = -126785.3098 EHBOND = 0.0000 RESTRAINT = 0.6319
 EAMBER (non-restraint) = -114722.1351
 EKCMT = 5078.1413 VIRIAL = 114047.4140 VOLUME = 175433.4875
                                                Density = 1.0038
 Ewald error estimate: 0.5301E-04
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 30
 Energy (this step): Bond = 0.632 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.651 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.766 0.766 0.775 0.775 0.566 0.566 0.575 0.575
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================

 NSTEP = 40 TIME(PS) = 232.610 TEMP(K) = 3142.52 PRESS =-81131.9
 Etot = -47182.2867 EKtot = 163551.1041 EPtot = -210733.3908
 BOND = 167.3398 ANGLE = 423.6836 DIHED = 514.9278
 1-4 NB = 178.1821 1-4 EEL = 2777.0986 VDWAALS = 10130.2256
 EELEC = -224925.4634 EHBOND = 0.0000 RESTRAINT = 0.6152
 EAMBER (non-restraint) = -210734.0060
 EKCMT = 5474.8584 VIRIAL = 304847.4455 VOLUME = 170900.1706
                                                Density = 1.0304
 Ewald error estimate: 0.2535E-04
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 40
 Energy (this step): Bond = 0.615 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.644 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.758 0.758 0.772 0.772 0.558 0.558 0.572 0.572
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================
 vlimit exceeded for step 42 ; vmax = 27.5954012646802589
 vlimit exceeded for step 43 ; vmax = 138.959560144215487
 vlimit exceeded for step 44 ; vmax = 119.776795686669814
 vlimit exceeded for step 45 ; vmax = 165.223684329653508
 vlimit exceeded for step 46 ; vmax = 180.570936430453798
 vlimit exceeded for step 47 ; vmax = 229.954197612266995
 vlimit exceeded for step 48 ; vmax = 749.379462585204124
 vlimit exceeded for step 49 ; vmax = 809.343960320155134

 NSTEP = 50 TIME(PS) = 232.612 TEMP(K) =100064.41 PRESS =********
 Etot = 4290807.7799 EKtot = 5207807.0642 EPtot = -916999.2843
 BOND = 179.1112 ANGLE = 452.0682 DIHED = 509.5124
 1-4 NB = 178.7702 1-4 EEL = 2779.2682 VDWAALS = 23717.5414
 EELEC = -944816.1529 EHBOND = 0.0000 RESTRAINT = 0.5969
 EAMBER (non-restraint) = -916999.8813
 EKCMT = 225135.2696 VIRIAL = 1273894.2971 VOLUME = 156000.2398
                                                Density = 1.1288
 Ewald error estimate: 0.1103E-04
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 50
 Energy (this step): Bond = 0.597 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.636 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.748 0.748 0.768 0.768 0.548 0.548 0.568 0.568
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================
 vlimit exceeded for step 50 ; vmax = 698.993868856018707
 vlimit exceeded for step 51 ; vmax = 414.910311906575089
 vlimit exceeded for step 52 ; vmax = 1628.37707913765939
 vlimit exceeded for step 53 ; vmax = 1547.51242049349730
 vlimit exceeded for step 54 ; vmax = 2698.52571209905454
 vlimit exceeded for step 55 ; vmax = 2696.23180467852671
 vlimit exceeded for step 56 ; vmax = 2295.21036036854639
 vlimit exceeded for step 57 ; vmax = 3185.33984517002500
 vlimit exceeded for step 58 ; vmax = 4707.94695676927222
 vlimit exceeded for step 59 ; vmax = 1693.37577817678016

 NSTEP = 60 TIME(PS) = 232.615 TEMP(K) =227552.14 PRESS =********
 Etot = 9367982.8928 EKtot = ************ EPtot = ************
 BOND = 620.6349 ANGLE = 677.5709 DIHED = 552.3036
 1-4 NB = 185.8595 1-4 EEL = 2781.7005 VDWAALS = 64956.2605
 EELEC = ************ EHBOND = 0.0000 RESTRAINT = 0.5781
 EAMBER (non-restraint) = ************
 EKCMT = 3048945.8319 VIRIAL = 5276292.3048 VOLUME = 165545.6053
                                                Density = 1.0637
 Ewald error estimate: 0.2075E-05
 ------------------------------------------------------------------------------

===============================================================================
                      NMR restraints for step 60
 Energy (this step): Bond = 0.578 Angle = 0.000 Torsion = 0.000
 Energy (tot. run) : Bond = 0.628 Angle = 0.000 Torsion = 0.000

 DEVIATIONS: Target=(r2+r3)/2 Target = closer of r2/r3
            This step Entire run This step Entire run
           ave. rms ave. rms ave. rms ave. rms
 Bond 0.739 0.739 0.764 0.764 0.539 0.539 0.564 0.564
 Angle 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 Torsion 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
===============================================================================
 vlimit exceeded for step 60 ; vmax = 1823.84265923123826
 vlimit exceeded for step 61 ; vmax = 525.756187843766384
 vlimit exceeded for step 62 ; vmax = 2541.60983804967918
 vlimit exceeded for step 63 ; vmax = 925.816197967856510
 EWALD BOMB in subroutine ewald_list
  volume of ucell too big, too many subcells
  list grid memory needs to be reallocated, restart sander

********************************************
********************************************
Here is the stderr output:
------------------------------
AIX 224% EWALD BOMB in subroutine ewald_list
  volume of ucell too big, too many subcells
  list grid memory needs to be reallocated, restart sander
 EWALD BOMB in subroutine ewald_list
 EWALD BOMB in subroutine ewald_list
  volume of ucell too big, too many subcells
  list grid memory needs to be reallocated, restart sander
  volume of ucell too big, too many subcells
  list grid memory needs to be reallocated, restart sander
ERROR: 0031-250 task 0: Terminated
ERROR: 0031-250 task 1: Terminated
ERROR: 0031-250 task 2: Terminated
ERROR: 0031-250 task 3: Terminated


********************************************
********************************************

Thank you,
Sincerely
Wayne

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 14 2004 - 15:53:11 PST
Custom Search