------------------------------------------------------- Amber 10 SANDER 2008 ------------------------------------------------------- | PMEMD implementation of SANDER, Release 10 | Run on 04/09/2010 at 16:15:58 [-O]verwriting output File Assignments: | MDIN: wrap.in | MDOUT: 151_E461K_MGATP_prod46.out | INPCRD: 151_E461K_MGATP_prod45.rst | PARM: E461K_MGATP_sol.prmtop | RESTRT: 151_E461K_MGATP_prod46.rst | REFC: refc | MDVEL: mdvel | MDEN: mden | MDCRD: 151_E461K_MGATP_prod46.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_EFS | DIRFRC_NOVEC | MPI | PUBFFT | BINTRAJ | Largest sphere to fit in unit cell has radius = 85.730 | New format PARM file being parsed. | Version = 1.000 Date = 09/27/09 Time = 02:40:15 | Duplicated 0 dihedrals | Duplicated 0 dihedrals -------------------------------------------------------------------------------- 1. RESOURCE USE: -------------------------------------------------------------------------------- getting new box info from bottom of inpcrd NATOM = 711968 NTYPES = 20 NBONH = 657760 MBONA = 54607 NTHETH = 127652 MTHETA = 73780 NPHIH = 238210 MPHIA = 173474 NHPARM = 0 NPARM = 0 NNB = 1409494 NRES = 207974 NBONA = 54607 NTHETA = 73780 NPHIA = 173474 NUMBND = 54 NUMANG = 108 NPTRA = 66 NATYP = 38 NPHB = 1 IFBOX = 2 NMXRS = 43 IFCAP = 0 NEXTRA = 0 NCOPY = 0 | Coordinate Index Table dimensions: 31 31 31 | Direct force subcell size = 6.7740 6.7740 6.7740 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 = 209.995 Box Y = 209.995 Box Z = 209.995 Alpha = 109.471 Beta = 109.471 Gamma = 109.471 NFFT1 = 216 NFFT2 = 216 NFFT3 = 216 Cutoff= 10.000 Tol =0.100E-05 Ewald Coefficient = 0.31234 Interpolation order = 4 | PMEMD ewald parallel performance parameters: | block_fft = 0 | fft_blk_y_divisor = 4 | excl_recip = 0 | excl_master = 0 | atm_redist_freq = 320 -------------------------------------------------------------------------------- 3. ATOMIC COORDINATES AND VELOCITIES -------------------------------------------------------------------------------- begin time read from input coords = 46310.000 ps Number of triangulated 3-point waters found: 200624 Sum of charges from parm topology file = ************ Assuming uniform neutralizing plasma | Dynamic Memory, Types Used: | Reals 44354169 | Integers 26441394 | Nonbonded Pairs Initial Allocation: 8694908 | Running AMBER/MPI version on 64 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.: -209.9950312 0.E+0 0.E+0 wrapping first mol.: -209.9950312 0.E+0 0.E+0 NSTEP = 1 TIME(PS) = 46310.001 TEMP(K) = 302.71 PRESS = 0.0 Etot = -1641721.5167 EKtot = 444578.8156 EPtot = -2086300.3323 BOND = 19761.1694 ANGLE = 52203.0058 DIHED = 75748.0769 1-4 NB = 25030.3891 1-4 EEL = 334702.1664 VDWAALS = 224843.4719 EELEC = -2818588.6119 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.4136E-04 ------------------------------------------------------------------------------ A V E R A G E S O V E R 1 S T E P S NSTEP = 1 TIME(PS) = 46310.001 TEMP(K) = 302.71 PRESS = 0.0 Etot = -1641721.5167 EKtot = 444578.8156 EPtot = -2086300.3323 BOND = 19761.1694 ANGLE = 52203.0058 DIHED = 75748.0769 1-4 NB = 25030.3891 1-4 EEL = 334702.1664 VDWAALS = 224843.4719 EELEC = -2818588.6119 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.4136E-04 ------------------------------------------------------------------------------ R M S F L U C T U A T I O N S NSTEP = 1 TIME(PS) = 46310.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.44 43.19 | Nonbond 0.48 46.54 | Bond 0.00 0.01 | Angle 0.00 0.10 | Dihedral 0.00 0.27 | Shake 0.00 0.46 | RunMD 0.10 9.42 | Other 0.00 0.00 | ------------------------------ | Total 1.02 | PME Nonbond Pairlist CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | Set Up Cit 0.11 11.23 | Build List 0.15 14.67 | --------------------------------- | Total 0.26 25.90 | PME Direct Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | NonBonded Calc 0.15 14.96 | Exclude Masked 0.00 0.29 | Other 0.01 0.92 | --------------------------------- | Total 0.17 16.17 | PME Reciprocal Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | 1D bspline 0.00 0.23 | Grid Charges 0.00 0.40 | Scalar Sum 0.02 1.55 | Gradient Sum 0.01 0.54 | FFT 0.02 1.74 | --------------------------------- | Total 0.05 4.47 | 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: 5.66 seconds | Master NonSetup CPU time: 4.71 seconds | Master Total CPU time: 10.37 seconds 0.00 hours | Master Setup wall time: 6 seconds | Master NonSetup wall time: 5 seconds | Master Total wall time: 11 seconds 0.00 hours