------------------------------------------------------- Amber 10 SANDER 2008 ------------------------------------------------------- | PMEMD implementation of SANDER, Release 10 | Run on 02/12/2009 at 19:41:58 [-O]verwriting output File Assignments: | MDIN: wrap.in | MDOUT: 123_1SVT_chainA_apo_prod358.out | INPCRD: 123_1SVT_chainA_apo_prod357.rst | PARM: 1SVT_chainA_sol.prmtop | RESTRT: 123_1SVT_chainA_apo_prod358.rst | REFC: refc | MDVEL: mdvel | MDEN: mden | MDCRD: 123_1SVT_chainA_apo_prod358.mdcrd | MDINFO: mdinfo |LOGFILE: logfile Here is the input file: MD simulation NVE &cntrl imin= 0, irest= 1, ntx = 5, ntb = 1, cut = 10, ntc = 2, ntf = 2, tol = 0.000001, ntt = 0, nstlim = 1, dt = 0.001, ntpr = 1, ntwx = 1, ntwr = 5, ioutfm = 1, iwrap = 1, &ewald dsum_tol = 0.000001, / | Conditional Compilation Defines Used: | DIRFRC_COMTRANS | DIRFRC_EFS | MPI | PUBFFT | BINTRAJ | Largest sphere to fit in unit cell has radius = 48.491 | New format PARM file being parsed. | Version = 1.000 Date = 04/14/08 Time = 16:53:02 | Duplicated 0 dihedrals | Duplicated 0 dihedrals -------------------------------------------------------------------------------- 1. RESOURCE USE: -------------------------------------------------------------------------------- getting new box info from bottom of inpcrd NATOM = 128405 NTYPES = 16 NBONH = 124549 MBONA = 3884 NTHETH = 9089 MTHETA = 5245 NPHIH = 16965 MPHIA = 12349 NHPARM = 0 NPARM = 0 NNB = 203966 NRES = 40714 NBONA = 3884 NTHETA = 5245 NPHIA = 12349 NUMBND = 37 NUMANG = 73 NPTRA = 43 NATYP = 28 NPHB = 1 IFBOX = 2 NMXRS = 24 IFCAP = 0 NEXTRA = 0 NCOPY = 0 | Coordinate Index Table dimensions: 17 17 17 | Direct force subcell size = 6.9870 6.9870 6.9870 BOX TYPE: TRUNCATED OCTAHEDRON -------------------------------------------------------------------------------- 2. CONTROL DATA FOR THE RUN -------------------------------------------------------------------------------- General flags: imin = 0, nmropt = 0 Nature and format of input: ntx = 5, irest = 1, ntrx = 1 Nature and format of output: ntxo = 1, ntpr = 1, ntrx = 1, ntwr = 5 iwrap = 1, ntwx = 1, ntwv = 0, ntwe = 0 ioutfm = 1, ntwprt = 0, idecomp = 0, rbornstat= 0 Potential function: ntf = 2, ntb = 1, igb = 0, nsnb = 25 ipol = 0, gbsa = 0, iesp = 0 dielc = 1.00000, cut = 10.00000, intdiel = 1.00000 scnb = 2.00000, scee = 1.20000 Frozen or restrained atoms: ibelly = 0, ntr = 0 Molecular dynamics: nstlim = 1, nscm = 1000, nrespa = 1 t = 0.00000, dt = 0.00100, vlimit = 20.00000 SHAKE: ntc = 2, jfastw = 0 tol = 0.00000 | Intermolecular bonds treatment: | no_intermolecular_bonds = 1 | Energy averages sample interval: | ene_avg_sampling = 1 Ewald parameters: verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1 vdwmeth = 1, eedmeth = 1, netfrc = 1 Box X = 118.778 Box Y = 118.778 Box Z = 118.778 Alpha = 109.471 Beta = 109.471 Gamma = 109.471 NFFT1 = 120 NFFT2 = 120 NFFT3 = 120 Cutoff= 10.000 Tol =0.100E-05 Ewald Coefficient = 0.31234 Interpolation order = 4 | PMEMD ewald parallel performance parameters: | block_fft = 1 | fft_blk_y_divisor = 4 | excl_recip = 1 | excl_master = 1 | atm_redist_freq = 32 -------------------------------------------------------------------------------- 3. ATOMIC COORDINATES AND VELOCITIES -------------------------------------------------------------------------------- begin time read from input coords =351800.000 ps Number of triangulated 3-point waters found: 40190 Sum of charges from parm topology file = -18.00000061 Assuming uniform neutralizing plasma | Dynamic Memory, Types Used: | Reals 4136087 | Integers 4132349 | Nonbonded Pairs Initial Allocation: 786078 | Running AMBER/MPI version on 128 nodes -------------------------------------------------------------------------------- 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.8987E-11 at 2.875760 --------------------------------------------------- |--------------------------------------------------- | APPROXIMATING direct energy using CUBIC SPLINE INTERPOLATION | with 50.0 points per unit in tabled values | Relative Error Limit not exceeded for r .gt. 2.33 | APPROXIMATING direct force using CUBIC SPLINE INTERPOLATION | with 50.0 points per unit in tabled values | Relative Error Limit not exceeded for r .gt. 2.77 |--------------------------------------------------- wrapping first mol.: -79.185529277405024 55.992636583360628 96.982083586904935 wrapping first mol.: -79.185529277405024 55.992636583360628 96.982083586904935 NSTEP = 1 TIME(PS) = 351800.001 TEMP(K) = 300.35 PRESS = 0.0 Etot = -318023.3232 EKtot = 77787.9915 EPtot = -395811.3147 BOND = 1416.0408 ANGLE = 3793.4961 DIHED = 5429.2274 1-4 NB = 1798.9917 1-4 EEL = 24374.6886 VDWAALS = 52627.6899 EELEC = -485251.4493 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.1780E-04 ------------------------------------------------------------------------------ A V E R A G E S O V E R 1 S T E P S NSTEP = 1 TIME(PS) = 351800.001 TEMP(K) = 300.35 PRESS = 0.0 Etot = -318023.3232 EKtot = 77787.9915 EPtot = -395811.3147 BOND = 1416.0408 ANGLE = 3793.4961 DIHED = 5429.2274 1-4 NB = 1798.9917 1-4 EEL = 24374.6886 VDWAALS = 52627.6899 EELEC = -485251.4493 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.1780E-04 ------------------------------------------------------------------------------ R M S F L U C T U A T I O N S NSTEP = 1 TIME(PS) = 351800.001 TEMP(K) = 0.00 PRESS = 0.0 Etot = 0.0000 EKtot = 0.0000 EPtot = 0.0000 BOND = 0.0000 ANGLE = 0.0000 DIHED = 0.0000 1-4 NB = 0.0000 1-4 EEL = 0.0000 VDWAALS = 0.0000 EELEC = 0.0000 EHBOND = 0.0000 RESTRAINT = 0.0000 |E(PBS) = 0.0000 Ewald error estimate: 0.0000E+00 ------------------------------------------------------------------------------ -------------------------------------------------------------------------------- 5. TIMINGS -------------------------------------------------------------------------------- | NonSetup CPU Time in Major Routines, Average for All Tasks: | | Routine Sec % | ------------------------------ | DataDistrib 0.02 24.34 | Nonbond 0.05 62.25 | Bond 0.00 0.00 | Angle 0.00 0.04 | Dihedral 0.00 0.24 | Shake 0.00 0.08 | RunMD 0.01 13.06 | Other 0.00 0.00 | ------------------------------ | Total 0.08 | PME Nonbond Pairlist CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | Set Up Cit 0.01 18.58 | Build List 0.01 18.34 | --------------------------------- | Total 0.03 36.92 | PME Direct Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | NonBonded Calc 0.01 18.26 | Exclude Masked 0.00 0.36 | Other 0.00 1.85 | --------------------------------- | Total 0.02 20.47 | PME Reciprocal Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | 1D bspline 0.00 0.24 | Grid Charges 0.00 0.47 | Scalar Sum 0.00 1.70 | Gradient Sum 0.00 0.59 | FFT 0.00 1.85 | --------------------------------- | Total 0.00 4.85 | PME Load Balancing CPU Time, Average for All Tasks: | | Routine Sec % | ------------------------------------ | Atom Reassign 0.00 0.00 | Image Reassign 0.00 0.00 | FFT Reassign 0.00 0.00 | ------------------------------------ | Total 0.00 0.00 | Master Setup CPU time: 1.01 seconds | Master NonSetup CPU time: 0.78 seconds | Master Total CPU time: 1.78 seconds 0.00 hours | Master Setup wall time: 1 seconds | Master NonSetup wall time: 1 seconds | Master Total wall time: 2 seconds 0.00 hours