[AMBER] Etot positive

From: Jio M <jiomm.yahoo.com>
Date: Sun, 15 Nov 2009 21:06:22 -0800 (PST)

Dear Amber users;



1) Is it necessary that Etot should be negative .... I have seen in
most of the querries and tutorial output, that Etot is negative.

I used the following input for equilibrium phase after heating from 0 - 300K

here is equilibrium input:

 &cntrl

  imin = 0, irest = 1, ntx = 7,

  ntb = 2, pres0 = 1.0, ntp = 1,

  taup = 2.0,

  cut = 10.0, ntr = 0,

  ntc = 2, ntf = 2,

  tempi = 300.0, temp0 = 300.0,

  ntt = 3, gamma_ln = 1.0,

  nstlim = 100000, dt = 0.002,

  ntpr = 100, ntwx = 100, ntwr = 1000



these last steps energy output values



NSTEP =    99900   TIME(PS)
=     249.800  TEMP(K) =  
295.30  PRESS =   -56.5

 Etot   =      6133.5908 
EKtot   =      3811.4443 
EPtot      =     
2322.1465

 BOND   =      
635.3000  ANGLE   =     
1956.2672  DIHED     
=       405.6718

 1-4 NB =       195.6267  1-4
EEL =      4000.9903 
VDWAALS    =     -4121.3236

 EELEC  =      -750.3858 
EHBOND  =        
0.0000  RESTRAINT 
=         0.0000

 EKCMT  =       673.1555 
VIRIAL  =       788.8184 
VOLUME     =     94811.4804

                                                   
Density   
=         1.2381

 Ewald error estimate:   0.5658E-03

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





 NSTEP =   100000   TIME(PS)
=     250.000  TEMP(K) =  
299.31  PRESS =    59.6

 Etot   =      6178.0066 
EKtot   =      3863.0849 
EPtot      =     
2314.9217

 BOND   =      
674.1848  ANGLE   =     
1878.4011  DIHED     
=       425.5195

 1-4 NB =       189.5570  1-4
EEL =      4010.5437 
VDWAALS    =     -4097.2245

 EELEC  =      -766.0599 
EHBOND  =        
0.0000  RESTRAINT 
=         0.0000

 EKCMT  =       685.0976 
VIRIAL  =       563.0999 
VOLUME     =     94823.4042

                                                   
Density   
=         1.2380

 Ewald error estimate:   0.1530E-03





Is there anything wrong in starting structures or input files??



please suggest





thanks and regards;

Jiomm




_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Nov 15 2009 - 21:30:02 PST
Custom Search