------------------------------------------------------- Amber 10 SANDER 2008 ------------------------------------------------------- | PMEMD implementation of SANDER, Release 10 | Run on 03/16/2009 at 18:08:31 [-O]verwriting output File Assignments: | MDIN: wrap.in | MDOUT: 121_1XCK_chainA_apo_prod144.out | INPCRD: 121_1XCK_chainA_apo_prod143.rst | PARM: 1XCK_chainA_sol.prmtop | RESTRT: 121_1XCK_chainA_apo_prod144.rst | REFC: refc | MDVEL: mdvel | MDEN: mden | MDCRD: 121_1XCK_chainA_apo_prod144.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 = 41.506 | New format PARM file being parsed. | Version = 1.000 Date = 04/07/08 Time = 13:13:41 | Duplicated 0 dihedrals | Duplicated 0 dihedrals -------------------------------------------------------------------------------- 1. RESOURCE USE: -------------------------------------------------------------------------------- getting new box info from bottom of inpcrd NATOM = 81053 NTYPES = 16 NBONH = 77197 MBONA = 3884 NTHETH = 9089 MTHETA = 5245 NPHIH = 16965 MPHIA = 12349 NHPARM = 0 NPARM = 0 NNB = 140830 NRES = 24930 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: 15 15 15 | Direct force subcell size = 6.7779 6.7779 6.7779 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 = 101.668 Box Y = 101.668 Box Z = 101.668 Alpha = 109.471 Beta = 109.471 Gamma = 109.471 NFFT1 = 108 NFFT2 = 108 NFFT3 = 108 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 =143050.000 ps Number of triangulated 3-point waters found: 24406 Sum of charges from parm topology file = -18.00000061 Assuming uniform neutralizing plasma | Dynamic Memory, Types Used: | Reals 2668127 | Integers 2810439 | Nonbonded Pairs Initial Allocation: 496196 | 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.: -33.889457232219968 95.85386877949297 0.E+0 wrapping first mol.: -33.889457232219968 95.85386877949297 0.E+0 NSTEP = 1 TIME(PS) = 143050.001 TEMP(K) = 301.09 PRESS = 0.0 Etot = -195649.5861 EKtot = 49648.7468 EPtot = -245298.3329 BOND = 1415.2799 ANGLE = 3717.3327 DIHED = 5385.9396 1-4 NB = 1805.3237 1-4 EEL = 24446.9481 VDWAALS = 29738.0322 EELEC = -311807.1891 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.1485E-04 ------------------------------------------------------------------------------ A V E R A G E S O V E R 1 S T E P S NSTEP = 1 TIME(PS) = 143050.001 TEMP(K) = 301.09 PRESS = 0.0 Etot = -195649.5861 EKtot = 49648.7468 EPtot = -245298.3329 BOND = 1415.2799 ANGLE = 3717.3327 DIHED = 5385.9396 1-4 NB = 1805.3237 1-4 EEL = 24446.9481 VDWAALS = 29738.0322 EELEC = -311807.1891 EHBOND = 0.0000 RESTRAINT = 0.0000 Ewald error estimate: 0.1485E-04 ------------------------------------------------------------------------------ R M S F L U C T U A T I O N S NSTEP = 1 TIME(PS) = 143050.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.01 26.25 | Nonbond 0.03 61.06 | Bond 0.00 0.00 | Angle 0.00 0.13 | Dihedral 0.00 0.31 | Shake 0.00 0.56 | RunMD 0.01 11.69 | Other 0.00 0.00 | ------------------------------ | Total 0.05 | PME Nonbond Pairlist CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | Set Up Cit 0.01 16.00 | Build List 0.01 18.63 | --------------------------------- | Total 0.02 34.62 | PME Direct Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | NonBonded Calc 0.01 19.00 | Exclude Masked 0.00 0.25 | Other 0.00 1.25 | --------------------------------- | Total 0.01 20.50 | PME Reciprocal Force CPU Time, Average for All Tasks: | | Routine Sec % | --------------------------------- | 1D bspline 0.00 0.38 | Grid Charges 0.00 0.62 | Scalar Sum 0.00 1.63 | Gradient Sum 0.00 0.88 | FFT 0.00 2.44 | --------------------------------- | Total 0.00 5.94 | 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: 0.68 seconds | Master NonSetup CPU time: 0.50 seconds | Master Total CPU time: 1.18 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