------------------------------------------------------- Amber 10 SANDER 2008 ------------------------------------------------------- | PMEMD implementation of SANDER, Release 10 | Run on 07/08/2008 at 23:41:17 [-O]verwriting output File Assignments: | MDIN: wrap.in | MDOUT: 128_1XCK_hepta_MGATP_prod52.out | INPCRD: 128_1XCK_hepta_MGATP_prod51.rst | PARM: 1XCK_hept_MGATP_sol.prmtop | RESTRT: 128_1XCK_hepta_MGATP_prod52.rst | REFC: refc | MDVEL: mdvel | MDEN: mden | MDCRD: 128_1XCK_hepta_MGATP_prod52.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 = 64.710 | New format PARM file being parsed. | Version = 1.000 Date = 05/19/08 Time = 13:14:30 | Duplicated 0 dihedrals | Duplicated 0 dihedrals -------------------------------------------------------------------------------- 1. RESOURCE USE: -------------------------------------------------------------------------------- getting new box info from bottom of inpcrd NATOM = 307984 NTYPES = 20 NBONH = 280768 MBONA = 27419 NTHETH = 63805 MTHETA = 37079 NPHIH = 119161 MPHIA = 86989 NHPARM = 0 NPARM = 0 NNB = 641062 NRES = 87959 NBONA = 27419 NTHETA = 37079 NPHIA = 86989 NUMBND = 54 NUMANG = 108 NPTRA = 66 NATYP = 38 NPHB = 1 IFBOX = 2 NMXRS = 43 IFCAP = 0 NEXTRA = 0 NCOPY = 0 | Coordinate Index Table dimensions: 23 23 23 | Direct force subcell size = 6.8916 6.8916 6.8916 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 = 158.507 Box Y = 158.507 Box Z = 158.507 Alpha = 109.471 Beta = 109.471 Gamma = 109.471 NFFT1 = 160 NFFT2 = 160 NFFT3 = 160 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 = 0 | excl_master = 0 | atm_redist_freq = 320 -------------------------------------------------------------------------------- 3. ATOMIC COORDINATES AND VELOCITIES -------------------------------------------------------------------------------- begin time read from input coords = 53300.000 ps Number of triangulated 3-point waters found: 84277 Sum of charges from parm topology file = ************ Assuming uniform neutralizing plasma | Dynamic Memory, Types Used: | Reals 9894098 | Integers 11900982 | Nonbonded Pairs Initial Allocation: 2873050 | Running AMBER/MPI version on 84 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.: -158.5073898 0.E+0 0.E+0 wrapping first mol.: -158.5073898 0.E+0 0.E+0 NSTEP = 1 TIME(PS) = 53300.001 TEMP(K) = 303.77 PRESS = 0.0 Etot = -702957.3505 EKtot = 194126.4554 EPtot = -897083.8060 BOND = 9903.0486 ANGLE = 26096.0764 DIHED = 38057.7134 1-4 NB = 12603.5356 1-4 EEL = 163151.0338 VDWAALS = 87773.4055 EELEC = -1234668.6192 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.8204E-04 ------------------------------------------------------------------------------ A V E R A G E S O V E R 1 S T E P S NSTEP = 1 TIME(PS) = 53300.001 TEMP(K) = 303.77 PRESS = 0.0 Etot = -702957.3505 EKtot = 194126.4554 EPtot = -897083.8060 BOND = 9903.0486 ANGLE = 26096.0764 DIHED = 38057.7134 1-4 NB = 12603.5356 1-4 EEL = 163151.0338 VDWAALS = 87773.4055 EELEC = -1234668.6192 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.8204E-04 ------------------------------------------------------------------------------ R M S F L U C T U A T I O N S NSTEP = 1 TIME(PS) = 53300.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.04 13.63 | Nonbond 0.19 72.22 | Bond 0.00 0.00 | Angle 0.00 0.18 | Dihedral 0.00 0.44 | Shake 0.00 0.51 | RunMD 0.03 12.99 | Other 0.00 0.02 | ------------------------------ | Total 0.26 | PME Nonbond Pairlist CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | Set Up Cit 0.04 16.87 | Build List 0.06 23.77 | --------------------------------- | Total 0.11 40.64 | PME Direct Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | NonBonded Calc 0.06 23.79 | Exclude Masked 0.00 0.38 | Other 0.00 1.10 | --------------------------------- | Total 0.07 25.27 | PME Reciprocal Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | 1D bspline 0.00 0.35 | Grid Charges 0.00 0.53 | Scalar Sum 0.00 1.87 | Gradient Sum 0.00 0.70 | FFT 0.01 2.87 | --------------------------------- | Total 0.02 6.31 | 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: 2.73 seconds | Master NonSetup CPU time: 1.90 seconds | Master Total CPU time: 4.64 seconds 0.00 hours | Master Setup wall time: 3 seconds | Master NonSetup wall time: 2 seconds | Master Total wall time: 5 seconds 0.00 hours