[gmx-users] christoph.freudenberger help me please!
    PeiQuan Chen 
    gromacs at 163.com
       
    Wed Apr 16 12:22:01 CEST 2003
    
    
  
Dear christoph.freudenberger:
  I have downloaded the Acetonitrile solvent box that you make in the www.gromacs.org,
and run a 3000ps MD on it.At last I have following your instructions to use the g_energy to 
calculate the potential enenergy! The result is as follow:
Statistics over 3000001 steps [ 0.0000 thru 3000.0002 ps ], 1 data sets
Energy                      Average       RMSD     Fluct.      Drift  Tot-Drift
-------------------------------------------------------------------------------
Potential                  -9147.37    61.3934    61.3661  0.0021147    6.34409
gcq#94: "Once Again Let Me Do This" (Urban Dance Squad)
so the DHvap can be caculated as the following:
DHvap=RT-Epot=[8.314*300-(-9147.37)] (J/mol)= 11.641 KJ/mol.
And that you result is DHvap=32.82 KJ/mol.
Can you tell me where is my wrong?
Sincerely yours
---- 
PeiQuan Chen
________________________________________________________________________
PeiQuan Chen 	Department of Chemistry, Nankai University
WeiJin Road 94, TianJin, 300071, P.R.China
phone:	(86)2223506728		Mobile Phone: (+86)13920675030
gromacs at 163.com	
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    
    
More information about the gromacs.org_gmx-users
mailing list