Weird energy values of LJ fluids simulation

From: Lindong Weng (
Date: Tue Dec 11 2012 - 13:43:16 CST

Hi, NAMDers,

I have figured out how to get the simulation of pure Ar or Kr fluids
started. Thanks for your help. But I got another problem afterwards, thats
about the results. The VdW energy term given by my simulation (~ -952
Kcal/mol or -3978 KJ/mol) is quite different from the literature (~ -5.4
KJ/mol). The total energy did not match neither. And the self-diffusion
coefficient of Ar or Kr atoms is also much smaller than the literature
(0.738 10e-9 m2/s in my study and 2.560 10e-9 m2/s in the literature). But
the density of the simulation system and the RDFs of various atom pairs are
nearly the same as the literature.

Do you have any ideas about the reasons? Thanks in advance.

Below is the force field I used, the parameters for L-J potential were
derived from =0.3405 nm and /B=119.8 K for Ar-Ar and =0.3633 nm and
/B=167.0 K for Kr-Kr. I think CHARMM has its own method to calculate

NONBONDED nbxmod 5 atom cdiel shift vatom vdistance vswitch -

cutnb 14.0 ctofnb 12.0 ctonnb 10.0 eps 1.0 e14fac 1.0 wmin 1.5

                !adm jr., 5/08/91, suggested cutoff scheme


!V(Lennard-Jones) = Eps,i,j[(Rmin,i,j/ri,j)**12 - 2(Rmin,i,j/ri,j)**6]


!epsilon: kcal/mole, Eps,i,j = sqrt(eps,i * eps,j)

!Rmin/2: A, Rmin,i,j = Rmin/2,i + Rmin/2,j


!atom ignored epsilon Rmin/2 ignored eps,1-4 Rmin/2,1-4


! ions, note lack of NBFIXes


AR 0.0 -0.2381 1.91099 ! Ar

KR 0.0 -0.3319 2.03895 ! Kr

HBOND CUTHB 0.5 ! If you want to do hbond analysis (only), then use

                 ! READ PARAM APPEND CARD

                 ! to append hbond parameters from the file: par_hbond.inp



UNC Charlotte

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:22:50 CST