weird energy values of LJ fluids

From: Lindong Weng (
Date: Tue Dec 11 2012 - 16:15:55 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. Is it related to the cutoff I used? I
used cutoff=3.5 in my simulation but the literature used 2.5 or even
less. But I found that the difference induced by cutoff value was minor.

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