[AMBER] Performance drop for large systems

From: Andreas Tosstorff <andreas.tosstorff.cup.uni-muenchen.de>
Date: Tue, 8 Nov 2016 11:11:46 +0100

Hi all,

 

I am running some simulations for a 240.000 atoms system on a Titan X
Pascal.

 

My 240.000 atoms: 16.01 ns/day.

Ambermd 408609 atoms Cellulose: 18 ns/day.

 

When running the smaller DHFR simulation, I achieve however comparable
results. With my system even outperforming the ambermd benchmark:

 

My DHFR: 552.81 ns/day

Ambermd DHFR: 542.58 ns/day

 

So I was wondering if there are any parameters I'm missing that have an
impact on the performance for my larger system.

 

 

Below are my output files for both runs.

 

-----------------------DHFR----------------------------------

 

 

          -------------------------------------------------------

          Amber 16 PMEMD 2016

          -------------------------------------------------------

 

| PMEMD implementation of SANDER, Release 16

 

| Run on 11/07/2016 at 18:20:59

 

| Executable path: pmemd.cuda

| Working directory: /storage/home/WIN/andtos

| Hostname: Unknown

  [-O]verwriting output

 

File Assignments:

| MDIN: ./Dropbox/PhD/Benchmark/2c2t_md.in


| MDOUT: ./Dropbox/PhD/Benchmark/2c2t_md.out


| INPCRD: ./Dropbox/PhD/Benchmark/2c2t_min.rst


| PARM: ./Dropbox/PhD/Benchmark/2c2t.parm7


| RESTRT: ./Dropbox/PhD/Benchmark/2c2t_md.rst


| REFC: refc


| MDVEL: mdvel


| MDEN: mden


| MDCRD: mdcrd


| MDINFO: mdinfo


| MDFRC: mdfrc


 

 

Here is the input file:

 

Typical Production MD NVE with


 GOOD energy conservation, HMR, 4fs.


 &cntrl


   ntx=1, irest=0,


   ntc=2, ntf=2, tol=0.000001,


   nstlim=75000,


   ntpr=1000, ntwx=1000,


   ntwr=10000,


   dt=0.004, cut=8.,


   ntt=0, ntb=1, ntp=0,


   ioutfm=1,


 /


 &ewald


   dsum_tol=0.000001,


 /


 

 

Note: ig = -1. Setting random seed to 754551 based on wallclock time in

      microseconds.

|--------------------- INFORMATION ----------------------

| GPU (CUDA) Version of PMEMD in use: NVIDIA GPU IN USE.

| Version 16.0.0

|

| 02/25/2016

|

| Implementation by:

| Ross C. Walker (SDSC)

| Scott Le Grand (nVIDIA)

|

| Precision model in use:

| [SPFP] - Single Precision Forces, 64-bit Fixed Point

| Accumulation. (Default)

|

|--------------------------------------------------------

|----------------- CITATION INFORMATION -----------------

|

| When publishing work that utilized the CUDA version

| of AMBER, please cite the following in addition to

| the regular AMBER citations:

|

| - Romelia Salomon-Ferrer; Andreas W. Goetz; Duncan

| Poole; Scott Le Grand; Ross C. Walker "Routine

| microsecond molecular dynamics simulations with

| AMBER - Part II: Particle Mesh Ewald", J. Chem.

| Theory Comput., 2013, 9 (9), pp3878-3888,

| DOI: 10.1021/ct400314y.

|

| - Andreas W. Goetz; Mark J. Williamson; Dong Xu;

| Duncan Poole; Scott Le Grand; Ross C. Walker

| "Routine microsecond molecular dynamics simulations

| with AMBER - Part I: Generalized Born", J. Chem.

| Theory Comput., 2012, 8 (5), pp1542-1555.

|

| - Scott Le Grand; Andreas W. Goetz; Ross C. Walker

| "SPFP: Speed without compromise - a mixed precision

| model for GPU accelerated molecular dynamics

| simulations.", Comp. Phys. Comm., 2013, 184

| pp374-380, DOI: 10.1016/j.cpc.2012.09.022

|

|--------------------------------------------------------

|------------------- GPU DEVICE INFO --------------------

|

| CUDA_VISIBLE_DEVICES: not set

| CUDA Capable Devices Detected: 2

| CUDA Device ID in use: 1

| CUDA Device Name: TITAN X (Pascal)

| CUDA Device Global Mem Size: 12189 MB

| CUDA Device Num Multiprocessors: 28

| CUDA Device Core Freq: 1.53 GHz

|

|--------------------------------------------------------

 

| Conditional Compilation Defines Used:

| PUBFFT

| BINTRAJ

| CUDA

| EMIL

 

| Largest sphere to fit in unit cell has radius = 29.217

 

| New format PARM file being parsed.

| Version = 1.000 Date = 11/07/16 Time = 18:20:23

 

| Note: 1-4 EEL scale factors are being read from the topology file.

 

| Note: 1-4 VDW scale factors are being read from the topology file.

| Duplicated 0 dihedrals

 

| Duplicated 0 dihedrals

 

----------------------------------------------------------------------------
----
   1.  RESOURCE   USE: 
----------------------------------------------------------------------------
----
 
getting box info from netcdf restart file
NATOM  =   24304 NTYPES =      16 NBONH =   22802 MBONA  =    1537
NTHETH =    3465 MTHETA =    2075 NPHIH =    6892 MPHIA  =    6536
NHPARM =       0 NPARM  =       0 NNB   =   45123 NRES   =    7283
NBONA  =    1537 NTHETA =    2075 NPHIA =    6536 NUMBND =      71
NUMANG =     161 NPTRA  =     201 NATYP =      36 NPHB   =       1
IFBOX  =       2 NMXRS  =      24 IFCAP =       0 NEXTRA =       0
NCOPY  =       0
 
| Coordinate Index Table dimensions:    14   14   14
| Direct force subcell size =     5.1119    5.1119    5.1119
 
     BOX TYPE: TRUNCATED OCTAHEDRON
 
----------------------------------------------------------------------------
----
   2.  CONTROL  DATA  FOR  THE  RUN
----------------------------------------------------------------------------
----
 
default_name
 
General flags:
     imin    =       0, nmropt  =       0
 
Nature and format of input:
     ntx     =       1, irest   =       0, ntrx    =       1
 
Nature and format of output:
     ntxo    =       2, ntpr    =    1000, ntrx    =       1, ntwr    =
10000
     iwrap   =       0, ntwx    =    1000, 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     =   8.00000, intdiel =   1.00000
 
Frozen or restrained atoms:
     ibelly  =       0, ntr     =       0
 
Molecular dynamics:
     nstlim  =     75000, nscm    =      1000, nrespa  =         1
     t       =   0.00000, dt      =   0.00400, vlimit  =  -1.00000
 
SHAKE:
     ntc     =       2, jfastw  =       0
     tol     =   0.00000
 
| Intermolecular bonds treatment:
|     no_intermolecular_bonds =       1
 
| Energy averages sample interval:
|     ene_avg_sampling =    1000
 
Ewald parameters:
     verbose =       0, ew_type =       0, nbflag  =       1, use_pme =
1
     vdwmeth =       1, eedmeth =       1, netfrc  =       1
     Box X =   71.567   Box Y =   71.567   Box Z =   71.567
     Alpha =  109.471   Beta  =  109.471   Gamma =  109.471
     NFFT1 =   72       NFFT2 =   72       NFFT3 =   72
     Cutoff=    8.000   Tol   =0.100E-05
     Ewald Coefficient =  0.39467
     Interpolation order =    4
 
----------------------------------------------------------------------------
----
   3.  ATOMIC COORDINATES AND VELOCITIES
----------------------------------------------------------------------------
----
 
default_name
 begin time read from input coords =     0.000 ps
 
 Number of triangulated 3-point waters found:     7097
 
     Sum of charges from parm topology file =  -0.00000007
     Forcing neutrality...
 
| Dynamic Memory, Types Used:
| Reals              745601
| Integers           812199
 
| Nonbonded Pairs Initial Allocation:     4059375
 
| GPU memory information (estimate):
| KB of GPU memory in use:     66637
| KB of CPU memory in use:     34224
 
----------------------------------------------------------------------------
----
   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.52
| 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.92
|---------------------------------------------------
 
NSTEP =        0   TIME(PS) =       0.000  TEMP(K) =     0.00  PRESS =
0.0
Etot   =    -95978.3214  EKtot   =         0.0000  EPtot      =
-95978.3214
BOND   =       111.6690  ANGLE   =       374.5261  DIHED      =
2117.9406
1-4 NB =       578.1100  1-4 EEL =      8180.9609  VDWAALS    =
13099.2755
EELEC  =   -120440.8034  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000487     0.00(Removed)
 
NSTEP =     1000   TIME(PS) =       4.000  TEMP(K) =    72.69  PRESS =
0.0
Etot   =    -83628.9618  EKtot   =      3618.7891  EPtot      =
-87247.7508
BOND   =       125.7101  ANGLE   =       461.9804  DIHED      =
2177.3242
1-4 NB =       603.1699  1-4 EEL =      8174.6980  VDWAALS    =
12930.3395
EELEC  =   -111720.9728  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000103     0.00(Removed)
 
NSTEP =     2000   TIME(PS) =       8.000  TEMP(K) =    77.59  PRESS =
0.0
Etot   =    -83642.1628  EKtot   =      3863.0061  EPtot      =
-87505.1689
BOND   =       133.4382  ANGLE   =       495.6591  DIHED      =
2171.3048
1-4 NB =       606.1061  1-4 EEL =      8165.8796  VDWAALS    =
13249.4005
EELEC  =   -112326.9571  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000136     0.00(Removed)
 
NSTEP =     3000   TIME(PS) =      12.000  TEMP(K) =    81.02  PRESS =
0.0
Etot   =    -83651.0804  EKtot   =      4033.5830  EPtot      =
-87684.6634
BOND   =       150.2517  ANGLE   =       514.5180  DIHED      =
2189.5612
1-4 NB =       602.7806  1-4 EEL =      8171.7382  VDWAALS    =
13168.0116
EELEC  =   -112481.5248  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000345     0.00(Removed)
 
NSTEP =     4000   TIME(PS) =      16.000  TEMP(K) =    82.62  PRESS =
0.0
Etot   =    -83653.3115  EKtot   =      4113.2148  EPtot      =
-87766.5263
BOND   =       158.1691  ANGLE   =       529.9045  DIHED      =
2192.2796
1-4 NB =       607.0153  1-4 EEL =      8183.1536  VDWAALS    =
13385.0126
EELEC  =   -112822.0610  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000514     0.00(Removed)
 
NSTEP =     5000   TIME(PS) =      20.000  TEMP(K) =    82.90  PRESS =
0.0
Etot   =    -83657.1927  EKtot   =      4127.2090  EPtot      =
-87784.4017
BOND   =       168.8873  ANGLE   =       560.4574  DIHED      =
2201.6659
1-4 NB =       614.7597  1-4 EEL =      8162.8854  VDWAALS    =
13379.4005
EELEC  =   -112872.4580  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000333     0.00(Removed)
 
NSTEP =     6000   TIME(PS) =      24.000  TEMP(K) =    83.94  PRESS =
0.0
Etot   =    -83659.4541  EKtot   =      4179.2251  EPtot      =
-87838.6792
BOND   =       177.1570  ANGLE   =       587.6517  DIHED      =
2182.1998
1-4 NB =       614.2575  1-4 EEL =      8184.5689  VDWAALS    =
13341.1927
EELEC  =   -112925.7068  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000259     0.00(Removed)
 
 
NSTEP =    49000   TIME(PS) =     196.000  TEMP(K) =    91.15  PRESS =
0.0
Etot   =    -83676.4700  EKtot   =      4538.2090  EPtot      =
-88214.6790
BOND   =       260.0450  ANGLE   =       756.7518  DIHED      =
2207.6491
1-4 NB =       622.8510  1-4 EEL =      8174.0375  VDWAALS    =
13893.8624
EELEC  =   -114129.8757  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000257     0.00(Removed)
 
NSTEP =    50000   TIME(PS) =     200.000  TEMP(K) =    91.67  PRESS =
0.0
Etot   =    -83678.6813  EKtot   =      4563.7769  EPtot      =
-88242.4581
BOND   =       257.0876  ANGLE   =       781.7082  DIHED      =
2200.4814
1-4 NB =       609.4546  1-4 EEL =      8172.0605  VDWAALS    =
13958.7950
EELEC  =   -114222.0454  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000291     0.00(Removed)
 
NSTEP =    51000   TIME(PS) =     204.000  TEMP(K) =    91.43  PRESS =
0.0
Etot   =    -83682.7109  EKtot   =      4551.9438  EPtot      =
-88234.6548
BOND   =       271.5436  ANGLE   =       752.8684  DIHED      =
2197.5209
1-4 NB =       608.3235  1-4 EEL =      8181.0471  VDWAALS    =
14027.6596
EELEC  =   -114273.6179  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000220     0.00(Removed)
 
NSTEP =    52000   TIME(PS) =     208.000  TEMP(K) =    92.20  PRESS =
0.0
Etot   =    -83683.0532  EKtot   =      4590.3027  EPtot      =
-88273.3559
BOND   =       275.4873  ANGLE   =       763.1347  DIHED      =
2197.0612
1-4 NB =       615.9724  1-4 EEL =      8162.1206  VDWAALS    =
14101.1722
EELEC  =   -114388.3043  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000292     0.00(Removed)
 
NSTEP =    53000   TIME(PS) =     212.000  TEMP(K) =    91.95  PRESS =
0.0
Etot   =    -83678.1396  EKtot   =      4578.0439  EPtot      =
-88256.1836
BOND   =       270.3665  ANGLE   =       771.0835  DIHED      =
2205.8472
1-4 NB =       618.9741  1-4 EEL =      8181.4440  VDWAALS    =
13962.6542
EELEC  =   -114266.5530  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000308     0.00(Removed)
 
NSTEP =    54000   TIME(PS) =     216.000  TEMP(K) =    93.04  PRESS =
0.0
Etot   =    -83682.5249  EKtot   =      4632.1533  EPtot      =
-88314.6783
BOND   =       253.2044  ANGLE   =       749.8910  DIHED      =
2211.2431
1-4 NB =       620.9964  1-4 EEL =      8170.8711  VDWAALS    =
13991.6644
EELEC  =   -114312.5488  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000228     0.00(Removed)
 
NSTEP =    55000   TIME(PS) =     220.000  TEMP(K) =    92.51  PRESS =
0.0
Etot   =    -83684.1123  EKtot   =      4605.5918  EPtot      =
-88289.7041
BOND   =       271.0046  ANGLE   =       740.1745  DIHED      =
2197.9832
1-4 NB =       611.5845  1-4 EEL =      8155.6705  VDWAALS    =
13958.9234
EELEC  =   -114225.0449  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000350     0.00(Removed)
 
NSTEP =    56000   TIME(PS) =     224.000  TEMP(K) =    92.42  PRESS =
0.0
Etot   =    -83682.5195  EKtot   =      4601.4458  EPtot      =
-88283.9653
BOND   =       266.7120  ANGLE   =       747.8794  DIHED      =
2224.5047
1-4 NB =       607.2874  1-4 EEL =      8184.1817  VDWAALS    =
13901.9492
EELEC  =   -114216.4797  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000338     0.00(Removed)
 
NSTEP =    57000   TIME(PS) =     228.000  TEMP(K) =    92.14  PRESS =
0.0
Etot   =    -83680.5615  EKtot   =      4587.1587  EPtot      =
-88267.7202
BOND   =       260.9190  ANGLE   =       749.1269  DIHED      =
2190.8955
1-4 NB =       618.5479  1-4 EEL =      8169.6077  VDWAALS    =
13973.1872
EELEC  =   -114230.0043  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000236     0.00(Removed)
 
NSTEP =    58000   TIME(PS) =     232.000  TEMP(K) =    92.26  PRESS =
0.0
Etot   =    -83681.0792  EKtot   =      4593.1919  EPtot      =
-88274.2710
BOND   =       263.3553  ANGLE   =       758.1997  DIHED      =
2215.2101
1-4 NB =       608.4316  1-4 EEL =      8168.6263  VDWAALS    =
14031.7889
EELEC  =   -114319.8829  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000170     0.00(Removed)
 
NSTEP =    59000   TIME(PS) =     236.000  TEMP(K) =    92.44  PRESS =
0.0
Etot   =    -83681.9664  EKtot   =      4602.2891  EPtot      =
-88284.2554
BOND   =       259.6166  ANGLE   =       767.2233  DIHED      =
2199.2495
1-4 NB =       615.4153  1-4 EEL =      8165.5614  VDWAALS    =
13984.7174
EELEC  =   -114276.0389  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000164     0.00(Removed)
 
NSTEP =    60000   TIME(PS) =     240.000  TEMP(K) =    91.98  PRESS =
0.0
Etot   =    -83682.7769  EKtot   =      4579.3384  EPtot      =
-88262.1153
BOND   =       249.8541  ANGLE   =       768.3124  DIHED      =
2206.7481
1-4 NB =       615.0716  1-4 EEL =      8182.0497  VDWAALS    =
13988.6436
EELEC  =   -114272.7946  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000189     0.00(Removed)
 
NSTEP =    61000   TIME(PS) =     244.000  TEMP(K) =    93.15  PRESS =
0.0
Etot   =    -83682.8932  EKtot   =      4637.6934  EPtot      =
-88320.5866
BOND   =       254.8020  ANGLE   =       772.0134  DIHED      =
2210.5500
1-4 NB =       620.6896  1-4 EEL =      8160.4999  VDWAALS    =
14049.1907
EELEC  =   -114388.3322  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000267     0.00(Removed)
 
NSTEP =    62000   TIME(PS) =     248.000  TEMP(K) =    92.73  PRESS =
0.0
Etot   =    -83679.8375  EKtot   =      4616.7720  EPtot      =
-88296.6095
BOND   =       261.2831  ANGLE   =       758.9755  DIHED      =
2211.3710
1-4 NB =       614.7107  1-4 EEL =      8175.4674  VDWAALS    =
13967.2000
EELEC  =   -114285.6172  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000235     0.00(Removed)
 
NSTEP =    63000   TIME(PS) =     252.000  TEMP(K) =    92.93  PRESS =
0.0
Etot   =    -83680.9919  EKtot   =      4626.6831  EPtot      =
-88307.6751
BOND   =       250.3440  ANGLE   =       757.7667  DIHED      =
2205.4394
1-4 NB =       609.2522  1-4 EEL =      8149.6342  VDWAALS    =
14000.9999
EELEC  =   -114281.1114  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000184     0.00(Removed)
 
NSTEP =    64000   TIME(PS) =     256.000  TEMP(K) =    92.83  PRESS =
0.0
Etot   =    -83683.0711  EKtot   =      4621.5508  EPtot      =
-88304.6219
BOND   =       239.5525  ANGLE   =       745.6270  DIHED      =
2220.8738
1-4 NB =       617.6087  1-4 EEL =      8172.3025  VDWAALS    =
13844.9444
EELEC  =   -114145.5307  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000173     0.00(Removed)
 
NSTEP =    65000   TIME(PS) =     260.000  TEMP(K) =    93.01  PRESS =
0.0
Etot   =    -83682.5638  EKtot   =      4630.7598  EPtot      =
-88313.3236
BOND   =       260.0190  ANGLE   =       746.4905  DIHED      =
2208.5661
1-4 NB =       624.4461  1-4 EEL =      8176.1744  VDWAALS    =
13898.9726
EELEC  =   -114227.9924  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000038     0.00(Removed)
 
NSTEP =    66000   TIME(PS) =     264.000  TEMP(K) =    92.64  PRESS =
0.0
Etot   =    -83681.2938  EKtot   =      4612.2617  EPtot      =
-88293.5555
BOND   =       257.1811  ANGLE   =       769.9406  DIHED      =
2196.7013
1-4 NB =       626.1907  1-4 EEL =      8157.3730  VDWAALS    =
14069.7699
EELEC  =   -114370.7122  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000097     0.00(Removed)
 
NSTEP =    67000   TIME(PS) =     268.000  TEMP(K) =    92.61  PRESS =
0.0
Etot   =    -83684.7753  EKtot   =      4610.5615  EPtot      =
-88295.3369
BOND   =       253.8543  ANGLE   =       753.2212  DIHED      =
2207.1275
1-4 NB =       621.2360  1-4 EEL =      8179.1783  VDWAALS    =
13991.5329
EELEC  =   -114301.4870  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000053     0.00(Removed)
 
NSTEP =    68000   TIME(PS) =     272.000  TEMP(K) =    92.25  PRESS =
0.0
Etot   =    -83679.0085  EKtot   =      4592.6260  EPtot      =
-88271.6344
BOND   =       247.6602  ANGLE   =       771.3300  DIHED      =
2216.5672
1-4 NB =       616.8963  1-4 EEL =      8194.8695  VDWAALS    =
14052.3802
EELEC  =   -114371.3380  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000173     0.00(Removed)
 
NSTEP =    69000   TIME(PS) =     276.000  TEMP(K) =    92.74  PRESS =
0.0
Etot   =    -83685.3891  EKtot   =      4617.1226  EPtot      =
-88302.5116
BOND   =       252.6687  ANGLE   =       765.1628  DIHED      =
2211.9533
1-4 NB =       608.5757  1-4 EEL =      8163.2522  VDWAALS    =
13906.7188
EELEC  =   -114210.8431  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000274     0.00(Removed)
 
NSTEP =    70000   TIME(PS) =     280.000  TEMP(K) =    93.10  PRESS =
0.0
Etot   =    -83678.7184  EKtot   =      4635.1396  EPtot      =
-88313.8581
BOND   =       257.8200  ANGLE   =       766.9005  DIHED      =
2206.5389
1-4 NB =       621.9698  1-4 EEL =      8183.3024  VDWAALS    =
13974.4145
EELEC  =   -114324.8042  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000164     0.00(Removed)
 
NSTEP =    71000   TIME(PS) =     284.000  TEMP(K) =    93.49  PRESS =
0.0
Etot   =    -83684.8987  EKtot   =      4654.3032  EPtot      =
-88339.2019
BOND   =       249.9093  ANGLE   =       725.3916  DIHED      =
2210.9032
1-4 NB =       620.0760  1-4 EEL =      8210.5792  VDWAALS    =
13841.0249
EELEC  =   -114197.0862  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000265     0.00(Removed)
 
NSTEP =    72000   TIME(PS) =     288.000  TEMP(K) =    93.76  PRESS =
0.0
Etot   =    -83686.1725  EKtot   =      4667.8823  EPtot      =
-88354.0548
BOND   =       257.4239  ANGLE   =       737.9259  DIHED      =
2207.0769
1-4 NB =       610.9468  1-4 EEL =      8190.1325  VDWAALS    =
13924.2696
EELEC  =   -114281.8304  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000179     0.00(Removed)
 
NSTEP =    73000   TIME(PS) =     292.000  TEMP(K) =    93.11  PRESS =
0.0
Etot   =    -83683.4873  EKtot   =      4635.7412  EPtot      =
-88319.2285
BOND   =       259.7581  ANGLE   =       759.5119  DIHED      =
2202.1605
1-4 NB =       603.1771  1-4 EEL =      8186.0657  VDWAALS    =
14086.0266
EELEC  =   -114415.9283  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000222     0.00(Removed)
 
NSTEP =    74000   TIME(PS) =     296.000  TEMP(K) =    93.53  PRESS =
0.0
Etot   =    -83681.9285  EKtot   =      4656.5522  EPtot      =
-88338.4808
BOND   =       262.9754  ANGLE   =       769.8954  DIHED      =
2196.6080
1-4 NB =       611.0025  1-4 EEL =      8164.7565  VDWAALS    =
14009.2751
EELEC  =   -114352.9936  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
check COM velocity, temp:        0.000137     0.00(Removed)
 
NSTEP =    75000   TIME(PS) =     300.000  TEMP(K) =    93.57  PRESS =
0.0
Etot   =    -83684.5976  EKtot   =      4658.5552  EPtot      =
-88343.1528
BOND   =       256.6035  ANGLE   =       777.0333  DIHED      =
2207.7610
1-4 NB =       623.0649  1-4 EEL =      8176.3008  VDWAALS    =
13948.7081
EELEC  =   -114332.6242  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
 
      A V E R A G E S   O V E R      75 S T E P S
 
 
NSTEP =    75000   TIME(PS) =     300.000  TEMP(K) =    89.78  PRESS =
0.0
Etot   =    -83676.6057  EKtot   =      4469.9454  EPtot      =
-88146.5511
BOND   =       239.2836  ANGLE   =       718.2113  DIHED      =
2204.9564
1-4 NB =       614.5756  1-4 EEL =      8173.7806  VDWAALS    =
13799.7408
EELEC  =   -113897.0993  EHBOND  =         0.0000  RESTRAINT  =
0.0000
----------------------------------------------------------------------------
--
 
 
      R M S  F L U C T U A T I O N S
 
 
NSTEP =    75000   TIME(PS) =     300.000  TEMP(K) =     3.75  PRESS =
0.0
Etot   =        10.2132  EKtot   =       186.6720  EPtot      =
196.5398
BOND   =        31.3975  ANGLE   =        68.9747  DIHED      =
9.9859
1-4 NB =         5.9357  1-4 EEL =        13.2575  VDWAALS    =
238.0087
EELEC  =       530.3662  EHBOND  =         0.0000  RESTRAINT  =
0.0000
|E(PBS) =        44.6263
----------------------------------------------------------------------------
--
 
----------------------------------------------------------------------------
----
   5.  TIMINGS
----------------------------------------------------------------------------
----
 
|  NonSetup CPU Time in Major Routines:
|
|     Routine           Sec        %
|     ------------------------------
|     Nonbond          10.34   22.72
|     Bond              0.00    0.00
|     Angle             0.00    0.00
|     Dihedral          0.00    0.00
|     Shake             0.86    1.90
|     RunMD            34.16   75.06
|     Other             0.15    0.32
|     ------------------------------
|     Total            45.51
 
|  PME Nonbond Pairlist CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     Set Up Cit           0.00    0.00
|     Build List           0.00    0.00
|     ---------------------------------
|     Total                0.00    0.00
 
|  PME Direct Force CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     NonBonded Calc       0.00    0.00
|     Exclude Masked       0.00    0.00
|     Other                0.08    0.17
|     ---------------------------------
|     Total                0.08    0.17
 
|  PME Reciprocal Force CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     1D bspline           0.00    0.00
|     Grid Charges         0.00    0.00
|     Scalar Sum           0.00    0.00
|     Gradient Sum         0.00    0.00
|     FFT                  0.00    0.00
|     ---------------------------------
|     Total                0.00    0.00
 
|  Final Performance Info:
|     -----------------------------------------------------
|     Average timings for last   75000 steps:
|     Elapsed(s) =      46.89 Per Step(ms) =       0.63
|         ns/day =     552.81   seconds/ns =     156.29
|
|     Average timings for all steps:
|     Elapsed(s) =      46.89 Per Step(ms) =       0.63
|         ns/day =     552.81   seconds/ns =     156.29
|     -----------------------------------------------------
 
|  Setup CPU time:            1.00 seconds
|  NonSetup CPU time:        45.51 seconds
|  Total CPU time:           46.51 seconds     0.01 hours
 
|  Setup wall time:           1    seconds
|  NonSetup wall time:       47    seconds
|  Total wall time:          48    seconds     0.01 hours
 
 
--------------------------------------240000
System--------------------------
 
 
 
          -------------------------------------------------------
          Amber 16 PMEMD                              2016
          -------------------------------------------------------
 
| PMEMD implementation of SANDER, Release 16
 
| Run on 11/07/2016 at 21:33:39
 
|   Executable path: pmemd.cuda
| Working directory: /storage/home/WIN/andtos
|          Hostname: Unknown
  [-O]verwriting output
 
File Assignments:
|   MDIN:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/mdin_prod.
|  MDOUT:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/TRZ_2016_9
| INPCRD:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/TRZ_2016_9
|   PARM:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/TRZ_2016_9
| RESTRT:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/TRZ_2016_9
|   REFC: refc
|  MDVEL: mdvel
|   MDEN: mden
|  MDCRD:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/TRZ_2016_9
| MDINFO: mdinfo
|  MDFRC: mdfrc
 
 
Here is the input file:
 
200 ps NPT production for 63 A
Constant pressure constant temperature production run
 &cntrl
  nstlim=500000, dt=0.002, ntx=5, irest=1, ntpr=5000, ntwr=10000, ntwx=5000,
  temp0=300.0, ntt=1, tautp=2.0, cut=8.,
 
  ntb=2, ntp=1,
 
  ntc=2, ntf=2, tol=0.000001,
 
  nmropt = 1, ioutfm = 1,
 
&end
 
&ewald
  dsum_tol=0.000001,
 /
 
&wt
  type='DUMPFREQ', istep1=50,
 &end
 &wt
  type='END',
 &end
DISANG=./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/disang.63
DUMPAVE=./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/distral_
63.
 
 
Note: ig = -1. Setting random seed to   402342 based on wallclock time in 
      microseconds.
|--------------------- INFORMATION ----------------------
| GPU (CUDA) Version of PMEMD in use: NVIDIA GPU IN USE.
|                    Version 16.0.0
| 
|                      02/25/2016
| 
| Implementation by:
|                    Ross C. Walker     (SDSC)
|                    Scott Le Grand     (nVIDIA)
| 
| Precision model in use:
|      [SPFP] - Single Precision Forces, 64-bit Fixed Point
|               Accumulation. (Default)
| 
|--------------------------------------------------------
|----------------- CITATION INFORMATION -----------------
|
|    When publishing work that utilized the CUDA version
|    of AMBER, please cite the following in addition to
|    the regular AMBER citations:
|
|  - Romelia Salomon-Ferrer; Andreas W. Goetz; Duncan
|    Poole; Scott Le Grand; Ross C. Walker "Routine
|    microsecond molecular dynamics simulations with
|    AMBER - Part II: Particle Mesh Ewald", J. Chem.
|    Theory Comput., 2013, 9 (9), pp3878-3888,
|    DOI: 10.1021/ct400314y.
|
|  - Andreas W. Goetz; Mark J. Williamson; Dong Xu;
|    Duncan Poole; Scott Le Grand; Ross C. Walker
|    "Routine microsecond molecular dynamics simulations
|    with AMBER - Part I: Generalized Born", J. Chem.
|    Theory Comput., 2012, 8 (5), pp1542-1555.
|
|  - Scott Le Grand; Andreas W. Goetz; Ross C. Walker
|    "SPFP: Speed without compromise - a mixed precision
|    model for GPU accelerated molecular dynamics
|    simulations.", Comp. Phys. Comm., 2013, 184
|    pp374-380, DOI: 10.1016/j.cpc.2012.09.022
|
|--------------------------------------------------------
|------------------- GPU DEVICE INFO --------------------
|
|            CUDA_VISIBLE_DEVICES: not set
|   CUDA Capable Devices Detected:      2
|           CUDA Device ID in use:      1
|                CUDA Device Name: TITAN X (Pascal)
|     CUDA Device Global Mem Size:  12189 MB
| CUDA Device Num Multiprocessors:     28
|           CUDA Device Core Freq:   1.53 GHz
|
|--------------------------------------------------------
 
| Conditional Compilation Defines Used:
| PUBFFT
| BINTRAJ
| CUDA
| EMIL
 
| Largest sphere to fit in unit cell has radius =    42.715
 
| New format PARM file being parsed.
| Version =    1.000 Date = 11/01/16 Time = 14:20:36
 
| Note: 1-4 EEL scale factors are being read from the topology file.
 
| Note: 1-4 VDW scale factors are being read from the topology file.
| Duplicated    0 dihedrals
 
| Duplicated    0 dihedrals
 
----------------------------------------------------------------------------
----
   1.  RESOURCE   USE: 
----------------------------------------------------------------------------
----
 
getting box info from netcdf restart file
NATOM  =  237560 NTYPES =      17 NBONH =  230892 MBONA  =    6826
NTHETH =   14656 MTHETA =    9302 NPHIH =   30538 MPHIA  =   29664
NHPARM =       0 NPARM  =       0 NNB   =  371832 NRES   =   75696
NBONA  =    6826 NTHETA =    9302 NPHIA =   29664 NUMBND =      70
NUMANG =     160 NPTRA  =     203 NATYP =      35 NPHB   =       1
IFBOX  =       1 NMXRS  =      24 IFCAP =       0 NEXTRA =       0
NCOPY  =       0
 
| Coordinate Index Table dimensions:    18   36   38
| Direct force subcell size =     4.7462    4.6011    4.5090
 
     BOX TYPE: RECTILINEAR
 
----------------------------------------------------------------------------
----
   2.  CONTROL  DATA  FOR  THE  RUN
----------------------------------------------------------------------------
----
 
default_name
 
General flags:
     imin    =       0, nmropt  =       1
 
Nature and format of input:
     ntx     =       5, irest   =       1, ntrx    =       1
 
Nature and format of output:
     ntxo    =       2, ntpr    =    5000, ntrx    =       1, ntwr    =
10000
     iwrap   =       0, ntwx    =    5000, ntwv    =       0, ntwe    =
0
     ioutfm  =       1, ntwprt  =       0, idecomp =       0, rbornstat=
0
 
Potential function:
     ntf     =       2, ntb     =       2, igb     =       0, nsnb    =
25
     ipol    =       0, gbsa    =       0, iesp    =       0
     dielc   =   1.00000, cut     =   8.00000, intdiel =   1.00000
 
Frozen or restrained atoms:
     ibelly  =       0, ntr     =       0
 
Molecular dynamics:
     nstlim  =    500000, nscm    =      1000, nrespa  =         1
     t       =   0.00000, dt      =   0.00200, vlimit  =  -1.00000
 
Berendsen (weak-coupling) temperature regulation:
     temp0   = 300.00000, tempi   =   0.00000, tautp   =   2.00000
 
Pressure regulation:
     ntp     =       1
     pres0   =   1.00000, comp    =  44.60000, taup    =   1.00000
 
SHAKE:
     ntc     =       2, jfastw  =       0
     tol     =   0.00000
 
NMR refinement options:
     iscale  =       0, noeskp  =       1, ipnlty  =       1, mxsub   =
1
     scalm   = 100.00000, pencut  =   0.10000, tausw   =   0.10000
 
| Intermolecular bonds treatment:
|     no_intermolecular_bonds =       1
 
| Energy averages sample interval:
|     ene_avg_sampling =    5000
 
Ewald parameters:
     verbose =       0, ew_type =       0, nbflag  =       1, use_pme =
1
     vdwmeth =       1, eedmeth =       1, netfrc  =       1
     Box X =   85.431   Box Y =  165.641   Box Z =  171.343
     Alpha =   90.000   Beta  =   90.000   Gamma =   90.000
     NFFT1 =   96       NFFT2 =  168       NFFT3 =  180
     Cutoff=    8.000   Tol   =0.100E-05
     Ewald Coefficient =  0.39467
     Interpolation order =    4
 
----------------------------------------------------------------------------
----
   3.  ATOMIC COORDINATES AND VELOCITIES
----------------------------------------------------------------------------
----
 
default_name
 begin time read from input coords =   500.000 ps
 
 
 
           Begin reading energy term weight changes/NMR restraints
WEIGHT CHANGES:
DUMPFREQ     50      0    0.000000    0.000000      0      0
                         ** No weight changes given **
 
RESTRAINTS:
Requested file redirections:
  DISANG    =
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/disang.63
  DUMPAVE   =
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/distral_63.d
Restraints will be read from file:
./Dropbox/PhD/Trastuzumab/Her_Modeller_2_FabHer/TRZ_2016_93/disang.63
Here are comments from the DISANG input file:
 
                       Number of restraints read =     1
 
                  Done reading weight changes/NMR restraints
 
 
 Number of triangulated 3-point waters found:    74808
 
     Sum of charges from parm topology file =  -0.00000063
     Forcing neutrality...
 
| Dynamic Memory, Types Used:
| Reals             9342683
| Integers          6974732
 
| Nonbonded Pairs Initial Allocation:    39678459
 
| GPU memory information (estimate):
| KB of GPU memory in use:    545960
| KB of CPU memory in use:    255089
 
----------------------------------------------------------------------------
----
   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.52
| 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.92
|---------------------------------------------------
check COM velocity, temp:        0.000709     0.00(Removed)
check COM velocity, temp:        0.000007     0.00(Removed)
check COM velocity, temp:        0.000023     0.00(Removed)
check COM velocity, temp:        0.000012     0.00(Removed)
check COM velocity, temp:        0.000014     0.00(Removed)
 
NSTEP =     5000   TIME(PS) =     510.000  TEMP(K) =   301.00  PRESS =
-8.9
Etot   =   -583426.5852  EKtot   =    144088.3125  EPtot      =
-727514.8977
BOND   =      2411.3904  ANGLE   =      6575.2942  DIHED      =
10924.3112
1-4 NB =      2986.6906  1-4 EEL =     37777.8614  VDWAALS    =
101938.5885
EELEC  =   -890129.0621  EHBOND  =         0.0000  RESTRAINT  =
0.0282
EAMBER (non-restraint)  =   -727514.9258
EKCMT  =     67518.3748  VIRIAL  =     67984.6864  VOLUME     =
2425373.0425
                                                    Density    =
0.9878
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.028   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000003     0.00(Removed)
check COM velocity, temp:        0.000009     0.00(Removed)
check COM velocity, temp:        0.000019     0.00(Removed)
check COM velocity, temp:        0.000007     0.00(Removed)
check COM velocity, temp:        0.000011     0.00(Removed)
 
NSTEP =    10000   TIME(PS) =     520.000  TEMP(K) =   299.79  PRESS =
56.0
Etot   =   -583429.1779  EKtot   =    143508.9062  EPtot      =
-726938.0841
BOND   =      2543.0911  ANGLE   =      6528.6244  DIHED      =
10924.5879
1-4 NB =      2924.9081  1-4 EEL =     37597.3950  VDWAALS    =
101766.8935
EELEC  =   -889224.2185  EHBOND  =         0.0000  RESTRAINT  =
0.6345
EAMBER (non-restraint)  =   -726938.7187
EKCMT  =     67254.8340  VIRIAL  =     64320.4999  VOLUME     =
2424852.8794
                                                    Density    =
0.9880
----------------------------------------------------------------------------
--
 
 
NMR restraints: Bond =    0.790   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000024     0.00(Removed)
check COM velocity, temp:        0.000008     0.00(Removed)
check COM velocity, temp:        0.000022     0.00(Removed)
check COM velocity, temp:        0.000043     0.00(Removed)
check COM velocity, temp:        0.000017     0.00(Removed)
 
NSTEP =    40000   TIME(PS) =     580.000  TEMP(K) =   300.87  PRESS =
-28.7
Etot   =   -583156.3506  EKtot   =    144027.8594  EPtot      =
-727184.2100
BOND   =      2428.8020  ANGLE   =      6373.5315  DIHED      =
10866.5979
1-4 NB =      2961.7757  1-4 EEL =     37626.7203  VDWAALS    =
101893.9364
EELEC  =   -889335.5749  EHBOND  =         0.0000  RESTRAINT  =
0.0012
EAMBER (non-restraint)  =   -727184.2112
EKCMT  =     67879.2284  VIRIAL  =     69383.9613  VOLUME     =
2425363.1708
                                                    Density    =
0.9878
----------------------------------------------------------------------------
--
 
 
NMR restraints: Bond =    0.000   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000030     0.00(Removed)
check COM velocity, temp:        0.000027     0.00(Removed)
check COM velocity, temp:        0.000016     0.00(Removed)
check COM velocity, temp:        0.000031     0.00(Removed)
check COM velocity, temp:        0.000027     0.00(Removed)
 
NSTEP =   445000   TIME(PS) =    1390.000  TEMP(K) =   299.28  PRESS =
-20.7
Etot   =   -583367.8274  EKtot   =    143267.2344  EPtot      =
-726635.0618
BOND   =      2460.4063  ANGLE   =      6568.2509  DIHED      =
10909.4730
1-4 NB =      2971.8806  1-4 EEL =     37555.3302  VDWAALS    =
101599.8766
EELEC  =   -888701.0002  EHBOND  =         0.0000  RESTRAINT  =
0.7209
EAMBER (non-restraint)  =   -726635.7827
EKCMT  =     67239.0759  VIRIAL  =     68324.6566  VOLUME     =
2425529.1213
                                                    Density    =
0.9877
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.721   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000029     0.00(Removed)
check COM velocity, temp:        0.000026     0.00(Removed)
check COM velocity, temp:        0.000028     0.00(Removed)
check COM velocity, temp:        0.000016     0.00(Removed)
check COM velocity, temp:        0.000024     0.00(Removed)
 
NSTEP =   450000   TIME(PS) =    1400.000  TEMP(K) =   299.75  PRESS =
-32.3
Etot   =   -583485.5216  EKtot   =    143488.4219  EPtot      =
-726973.9434
BOND   =      2425.2339  ANGLE   =      6439.1236  DIHED      =
10883.5475
1-4 NB =      2929.6868  1-4 EEL =     37556.9613  VDWAALS    =
101519.2403
EELEC  =   -888727.7467  EHBOND  =         0.0000  RESTRAINT  =
0.0099
EAMBER (non-restraint)  =   -726973.9533
EKCMT  =     67535.1148  VIRIAL  =     69225.2684  VOLUME     =
2424452.5393
                                                    Density    =
0.9881
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.010   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000006     0.00(Removed)
check COM velocity, temp:        0.000028     0.00(Removed)
check COM velocity, temp:        0.000005     0.00(Removed)
check COM velocity, temp:        0.000024     0.00(Removed)
check COM velocity, temp:        0.000012     0.00(Removed)
 
NSTEP =   455000   TIME(PS) =    1410.000  TEMP(K) =   300.09  PRESS =
2.9
Etot   =   -583560.6779  EKtot   =    143652.5000  EPtot      =
-727213.1779
BOND   =      2446.8756  ANGLE   =      6449.9512  DIHED      =
10857.9126
1-4 NB =      2945.6408  1-4 EEL =     37626.2206  VDWAALS    =
101714.3762
EELEC  =   -889254.2061  EHBOND  =         0.0000  RESTRAINT  =
0.0512
EAMBER (non-restraint)  =   -727213.2291
EKCMT  =     67352.9578  VIRIAL  =     67202.6006  VOLUME     =
2422655.2974
                                                    Density    =
0.9889
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.051   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000037     0.00(Removed)
check COM velocity, temp:        0.000016     0.00(Removed)
check COM velocity, temp:        0.000018     0.00(Removed)
check COM velocity, temp:        0.000013     0.00(Removed)
check COM velocity, temp:        0.000014     0.00(Removed)
 
NSTEP =   460000   TIME(PS) =    1420.000  TEMP(K) =   299.57  PRESS =
-9.4
Etot   =   -583579.4939  EKtot   =    143402.8438  EPtot      =
-726982.3376
BOND   =      2560.1740  ANGLE   =      6402.3817  DIHED      =
10851.2438
1-4 NB =      2953.5230  1-4 EEL =     37665.3074  VDWAALS    =
101068.2431
EELEC  =   -888484.0765  EHBOND  =         0.0000  RESTRAINT  =
0.8660
EAMBER (non-restraint)  =   -726983.2036
EKCMT  =     67038.9020  VIRIAL  =     67530.9966  VOLUME     =
2421554.1060
                                                    Density    =
0.9893
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.866   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000005     0.00(Removed)
check COM velocity, temp:        0.000030     0.00(Removed)
check COM velocity, temp:        0.000014     0.00(Removed)
check COM velocity, temp:        0.000016     0.00(Removed)
check COM velocity, temp:        0.000033     0.00(Removed)
 
NSTEP =   465000   TIME(PS) =    1430.000  TEMP(K) =   300.08  PRESS =
-15.7
Etot   =   -583582.6762  EKtot   =    143647.4219  EPtot      =
-727230.0981
BOND   =      2549.6391  ANGLE   =      6499.1456  DIHED      =
10852.2084
1-4 NB =      2938.9084  1-4 EEL =     37535.9829  VDWAALS    =
101133.8505
EELEC  =   -888740.3070  EHBOND  =         0.0000  RESTRAINT  =
0.4741
EAMBER (non-restraint)  =   -727230.5723
EKCMT  =     67698.5799  VIRIAL  =     68518.1682  VOLUME     =
2425240.1684
                                                    Density    =
0.9878
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.474   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000030     0.00(Removed)
check COM velocity, temp:        0.000035     0.00(Removed)
check COM velocity, temp:        0.000016     0.00(Removed)
check COM velocity, temp:        0.000026     0.00(Removed)
check COM velocity, temp:        0.000015     0.00(Removed)
 
NSTEP =   470000   TIME(PS) =    1440.000  TEMP(K) =   300.02  PRESS =
8.4
Etot   =   -583726.1146  EKtot   =    143617.6562  EPtot      =
-727343.7709
BOND   =      2435.9157  ANGLE   =      6584.2141  DIHED      =
10891.6130
1-4 NB =      2975.6154  1-4 EEL =     37579.9380  VDWAALS    =
101922.0598
EELEC  =   -889734.8197  EHBOND  =         0.0000  RESTRAINT  =
1.6929
EAMBER (non-restraint)  =   -727345.4637
EKCMT  =     67276.4848  VIRIAL  =     66838.2977  VOLUME     =
2422698.7705
                                                    Density    =
0.9888
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    1.693   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000018     0.00(Removed)
check COM velocity, temp:        0.000032     0.00(Removed)
check COM velocity, temp:        0.000022     0.00(Removed)
check COM velocity, temp:        0.000018     0.00(Removed)
check COM velocity, temp:        0.000021     0.00(Removed)
 
NSTEP =   475000   TIME(PS) =    1450.000  TEMP(K) =   300.61  PRESS =
10.3
Etot   =   -583577.7926  EKtot   =    143904.2031  EPtot      =
-727481.9958
BOND   =      2487.1492  ANGLE   =      6432.8062  DIHED      =
10856.2220
1-4 NB =      2938.9757  1-4 EEL =     37650.2697  VDWAALS    =
102196.3844
EELEC  =   -890043.8582  EHBOND  =         0.0000  RESTRAINT  =
0.0552
EAMBER (non-restraint)  =   -727482.0510
EKCMT  =     67378.7643  VIRIAL  =     66839.1539  VOLUME     =
2425083.5180
                                                    Density    =
0.9879
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.055   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000017     0.00(Removed)
check COM velocity, temp:        0.000021     0.00(Removed)
check COM velocity, temp:        0.000019     0.00(Removed)
check COM velocity, temp:        0.000023     0.00(Removed)
check COM velocity, temp:        0.000028     0.00(Removed)
 
NSTEP =   480000   TIME(PS) =    1460.000  TEMP(K) =   299.19  PRESS =
121.2
Etot   =   -583521.5779  EKtot   =    143224.1719  EPtot      =
-726745.7498
BOND   =      2474.7239  ANGLE   =      6427.4124  DIHED      =
10873.5902
1-4 NB =      2931.8030  1-4 EEL =     37629.7202  VDWAALS    =
102321.8948
EELEC  =   -889404.9305  EHBOND  =         0.0000  RESTRAINT  =
0.0363
EAMBER (non-restraint)  =   -726745.7861
EKCMT  =     67133.5011  VIRIAL  =     60792.5776  VOLUME     =
2423113.2470
                                                    Density    =
0.9887
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.036   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000029     0.00(Removed)
check COM velocity, temp:        0.000030     0.00(Removed)
check COM velocity, temp:        0.000028     0.00(Removed)
check COM velocity, temp:        0.000007     0.00(Removed)
check COM velocity, temp:        0.000013     0.00(Removed)
 
NSTEP =   485000   TIME(PS) =    1470.000  TEMP(K) =   299.74  PRESS =
-17.9
Etot   =   -583382.3932  EKtot   =    143487.7656  EPtot      =
-726870.1589
BOND   =      2455.8525  ANGLE   =      6417.7771  DIHED      =
10867.2605
1-4 NB =      2904.0877  1-4 EEL =     37722.4489  VDWAALS    =
101360.4174
EELEC  =   -888598.2112  EHBOND  =         0.0000  RESTRAINT  =
0.2083
EAMBER (non-restraint)  =   -726870.3671
EKCMT  =     67301.7343  VIRIAL  =     68239.2145  VOLUME     =
2426758.1399
                                                    Density    =
0.9872
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.208   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000017     0.00(Removed)
check COM velocity, temp:        0.000018     0.00(Removed)
check COM velocity, temp:        0.000021     0.00(Removed)
check COM velocity, temp:        0.000021     0.00(Removed)
check COM velocity, temp:        0.000014     0.00(Removed)
 
NSTEP =   490000   TIME(PS) =    1480.000  TEMP(K) =   300.29  PRESS =
13.2
Etot   =   -583514.8671  EKtot   =    143750.0312  EPtot      =
-727264.8984
BOND   =      2460.2323  ANGLE   =      6430.9124  DIHED      =
10858.7648
1-4 NB =      2922.3655  1-4 EEL =     37578.8979  VDWAALS    =
101687.4587
EELEC  =   -889203.5948  EHBOND  =         0.0000  RESTRAINT  =
0.0648
EAMBER (non-restraint)  =   -727264.9632
EKCMT  =     67555.6103  VIRIAL  =     66866.9070  VOLUME     =
2424430.6879
                                                    Density    =
0.9881
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.065   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000027     0.00(Removed)
check COM velocity, temp:        0.000037     0.00(Removed)
check COM velocity, temp:        0.000019     0.00(Removed)
check COM velocity, temp:        0.000027     0.00(Removed)
check COM velocity, temp:        0.000027     0.00(Removed)
 
NSTEP =   495000   TIME(PS) =    1490.000  TEMP(K) =   301.57  PRESS =
-99.9
Etot   =   -583482.3502  EKtot   =    144363.6406  EPtot      =
-727845.9908
BOND   =      2469.9403  ANGLE   =      6483.8895  DIHED      =
10871.7995
1-4 NB =      2960.6924  1-4 EEL =     37576.6710  VDWAALS    =
101119.1546
EELEC  =   -889328.2254  EHBOND  =         0.0000  RESTRAINT  =
0.0874
EAMBER (non-restraint)  =   -727846.0783
EKCMT  =     67754.6141  VIRIAL  =     72981.9763  VOLUME     =
2424439.7535
                                                    Density    =
0.9881
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.087   Angle =     0.000   Torsion =     0.000
============================================================================
===
check COM velocity, temp:        0.000025     0.00(Removed)
check COM velocity, temp:        0.000010     0.00(Removed)
check COM velocity, temp:        0.000023     0.00(Removed)
check COM velocity, temp:        0.000014     0.00(Removed)
check COM velocity, temp:        0.000048     0.00(Removed)
 
NSTEP =   500000   TIME(PS) =    1500.000  TEMP(K) =   299.39  PRESS =
31.7
Etot   =   -583440.6474  EKtot   =    143320.0938  EPtot      =
-726760.7412
BOND   =      2496.2454  ANGLE   =      6394.2353  DIHED      =
10850.7954
1-4 NB =      2966.9909  1-4 EEL =     37625.2657  VDWAALS    =
101684.4320
EELEC  =   -888778.8052  EHBOND  =         0.0000  RESTRAINT  =
0.0993
EAMBER (non-restraint)  =   -726760.8405
EKCMT  =     67272.7750  VIRIAL  =     65612.3020  VOLUME     =
2423361.3665
                                                    Density    =
0.9886
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.099   Angle =     0.000   Torsion =     0.000
============================================================================
===
 
      A V E R A G E S   O V E R     100 S T E P S
 
 
NSTEP =   500000   TIME(PS) =    1500.000  TEMP(K) =   300.01  PRESS =
2.2
Etot   =   -583424.8011  EKtot   =    143612.6097  EPtot      =
-727037.4108
BOND   =      2470.7500  ANGLE   =      6444.4494  DIHED      =
10888.6137
1-4 NB =      2954.5731  1-4 EEL =     37607.9120  VDWAALS    =
101754.1528
EELEC  =   -889158.1626  EHBOND  =         0.0000  RESTRAINT  =
0.3007
EAMBER (non-restraint)  =   -727037.7115
EKCMT  =     67332.8792  VIRIAL  =     67215.3502  VOLUME     =
2424878.1713
                                                    Density    =
0.9880
----------------------------------------------------------------------------
--
 
NMR restraints: Bond =    0.099   Angle =     0.000   Torsion =     0.000
============================================================================
===
 
      R M S  F L U C T U A T I O N S
 
 
NSTEP =   500000   TIME(PS) =    1500.000  TEMP(K) =     0.52  PRESS =
49.7
Etot   =        99.0740  EKtot   =       250.9695  EPtot      =
270.3225
BOND   =        41.4968  ANGLE   =        60.1593  DIHED      =
36.6356
1-4 NB =        23.3731  1-4 EEL =        69.0685  VDWAALS    =
383.9619
EELEC  =       493.1101  EHBOND  =         0.0000  RESTRAINT  =
0.3833
EAMBER (non-restraint)  =       269.9393
EKCMT  =       188.6404  VIRIAL  =      2606.6294  VOLUME     =
1245.9789
                                                    Density    =
0.0005
----------------------------------------------------------------------------
--
 
 
NMR restraints on final step:
 
----------------------------------------------------------------------------
----
   5.  TIMINGS
----------------------------------------------------------------------------
----
 
|  NonSetup CPU Time in Major Routines:
|
|     Routine           Sec        %
|     ------------------------------
|     Nonbond        4686.16   87.23
|     Bond              0.00    0.00
|     Angle             0.00    0.00
|     Dihedral          0.00    0.00
|     Shake             8.36    0.16
|     RunMD           675.67   12.58
|     Other             1.98    0.04
|     ------------------------------
|     Total          5372.16
 
|  PME Nonbond Pairlist CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     Set Up Cit           0.00    0.00
|     Build List           0.00    0.00
|     ---------------------------------
|     Total                0.00    0.00
 
|  PME Direct Force CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     NonBonded Calc       0.00    0.00
|     Exclude Masked       0.00    0.00
|     Other                1.20    0.02
|     ---------------------------------
|     Total                1.20    0.02
 
|  PME Reciprocal Force CPU Time:
|
|     Routine              Sec        %
|     ---------------------------------
|     1D bspline           0.00    0.00
|     Grid Charges         0.00    0.00
|     Scalar Sum           0.00    0.00
|     Gradient Sum         0.00    0.00
|     FFT                  0.00    0.00
|     ---------------------------------
|     Total                0.00    0.00
 
|  Final Performance Info:
|     -----------------------------------------------------
|     Average timings for last    5000 steps:
|     Elapsed(s) =      55.33 Per Step(ms) =      11.07
|         ns/day =      15.61   seconds/ns =    5533.40
|
|     Average timings for all steps:
|     Elapsed(s) =    5396.43 Per Step(ms) =      10.79
|         ns/day =      16.01   seconds/ns =    5396.43
|     -----------------------------------------------------
 
|  Setup CPU time:            2.31 seconds
|  NonSetup CPU time:      5372.16 seconds
|  Total CPU time:         5374.47 seconds     1.49 hours
 
|  Setup wall time:           3    seconds
|  NonSetup wall time:     5396    seconds
|  Total wall time:        5399    seconds     1.50 hours
 
 
 
------------------------------------------------------------------
 
 
 
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Nov 08 2016 - 02:30:02 PST
Custom Search