[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