AMBER: comparing versions 6.0 and 7.0 of AMBER

From: Therese Malliavin <Therese.Malliavin.ibpc.fr>
Date: Mon, 27 Oct 2003 10:30:21 +0100

Dear AMBER-Netters,

I am now testing the version 7.0 of AMBER and comparing it to version
6.0

1) I was really happy that, on a cluster of Power4 IBM, the CPU time was
30% shorter for 7.0 than for 6.0 in a MD simulation.

2) I realized that the ndfmin and ntcm input parameters are no more
existing in AMBER 7.0, which is to my sense an improvement. Are there
other modifications of the input between the 2 versions, and where are
they listed (I looked in the Web site, but did not find anything) ?

3) I compared the results obtained with AMBER 6.0 and 7.0 starting with
the same input. The mean values of energies, temperature and pressure,
on 300 ps of MD simulation are of the same order, but slightly different
(see below). The differences are already visible after the first 250
steps of the simulation. Did you try similar comparisons ?

Thank you for your answers,

Thérèse Malliavin
Laboratoire de Biochimie Théorique
IBPC, Paris
France


AMBER 7.0
 NSTEP = 250 TIME(PS) = 35.500 TEMP(K) = 296.68 PRESS =
44.2
 Etot = -51368.9961 EKtot = 11952.8837 EPtot =
-63321.8797
 BOND = 250.1802 ANGLE = 803.0923 DIHED =
535.9939
 1-4 NB = 330.3660 1-4 EEL = 1450.1109 VDWAALS =
8012.4293
 EELEC = -74704.0523 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 5591.7403 VIRIAL = 5401.0442 VOLUME =
200033.5357
                                                Density =
1.0057
 Ewald error estimate: 0.1345E-03

------------------------------------------------------------------------------
 
      A V E R A G E S O V E R 150000 S T E P S
 
 
 NSTEP = 150000 TIME(PS) = 335.000 TEMP(K) = 297.84 PRESS =
1.1
 Etot = -51675.2585 EKtot = 11999.7269 EPtot =
-63674.9854
 BOND = 286.8138 ANGLE = 828.2833 DIHED =
492.1864
 1-4 NB = 319.9926 1-4 EEL = 1560.9264 VDWAALS =
7965.0759
 EELEC = -75128.2637 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 5487.5470 VIRIAL = 5482.8556 VOLUME =
199336.7664
                                                Density =
1.0093
 Ewald error estimate: 0.8209E-04

------------------------------------------------------------------------------
 
 
      R M S F L U C T U A T I O N S
 
 
 NSTEP = 150000 TIME(PS) = 335.000 TEMP(K) = 1.64 PRESS =
108.4
 Etot = 58.0769 EKtot = 65.8940 EPtot =
87.0540
 BOND = 14.2855 ANGLE = 20.7062 DIHED =
13.3836
 1-4 NB = 7.7735 1-4 EEL = 40.9783 VDWAALS =
87.1298
 EELEC = 135.5214 EHBOND = 0.0000 RESTRAINT =
0.0000
 EKCMT = 51.5177 VIRIAL = 472.8587 VOLUME =
504.3893
                                                Density =
0.0026
 Ewald error estimate: 0.6058E-04

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

AMBER 6.0
 
 NSTEP = 250 TIME(PS) = 35.500 TEMP(K) = 299.69 PRESS =
32.84
 Etot = -51294.4939 EKtot = 12075.1683 EPtot =
-63369.6622
 BOND = 253.2897 ANGLE = 805.5442 DIHED =
537.8691
 1-4 NB = 330.0391 1-4 EEL = 1452.1585 VDWAALS =
7948.0344
 EELEC = -74696.5973 EHBOND = 0.0000 CONSTRAINT =
0.0000
 EKCMT = 5557.3655 VIRIAL = 5415.6249 VOLUME =
199878.1326
                                                Density =
1.0065
 Ewald error estimate: 0.1301E-04

------------------------------------------------------------------------------
 
      A V E R A G E S O V E R ***** S T E P S
 
 
 NSTEP =150000 TIME(PS) = 335.000 TEMP(K) = 298.00 PRESS =
0.99
 Etot = -51900.6808 EKtot = 12006.9748 EPtot =
-63907.6556
 BOND = 286.7658 ANGLE = 827.3234 DIHED =
483.6268
 1-4 NB = 320.9856 1-4 EEL = 1548.5629 VDWAALS =
8024.9745
 EELEC = -75399.8945 EHBOND = 0.0000 CONSTRAINT =
0.0000
 EKCMT = 5438.0421 VIRIAL = 5433.7916 VOLUME =
198855.8801
                                                Density =
1.0117
 Ewald error estimate: 0.7820E-04

------------------------------------------------------------------------------
 
 
      R M S F L U C T U A T I O N S
 
 
 NSTEP =150000 TIME(PS) = 335.000 TEMP(K) = 1.61 PRESS =
118.65
 Etot = 86.3695 EKtot = 65.0466 EPtot =
102.2181
 BOND = 14.3055 ANGLE = 20.6284 DIHED =
15.7263
 1-4 NB = 7.9256 1-4 EEL = 35.5736 VDWAALS =
86.0490
 EELEC = 149.9105 EHBOND = 0.0000 CONSTRAINT =
0.0000
 EKCMT = 50.7803 VIRIAL = 515.9662 VOLUME =
548.7427
                                                Density =
0.0028
 Ewald error estimate: 0.5846E-04

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

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Mon Oct 27 2003 - 09:53:01 PST
Custom Search