Dear LAMMPS users,
I am using the serial version of LAMMPS (Jan 23, 2009) on cygwin on Windows
XP with additional packages: opt, molecule, kspace, class2, manybody.
I am simulating 125 Lennard Jones particles (epsilon = sigma = 1) in LJ
units at the number density of 0.85 at reduced temperature of 1.2 in NVE
ensemble. The volume of my simulation box is 147.059.
I'd appreciate if you can help me resolve the the following unexpected
observations:
1. I used compute stress/atom command to calculate the pressure as described in
the documentation. However, the pressure calculated by LAMMPS is not same as the
pressure calculated using the 'compute' command, see column #3 and #4 in the
output. As mentioned in the documentation, these two values should match. Here
is the link to this documentation page:
http://lammps.sandia.gov/doc/compute_stress_atom.html
2. The 'sum' mode of 'compute reduce' outputs the sum divided by the number of
particles for the global stress components, see the column #4, #5 and #6 in the
output. This is not the expected behavior as explained the manual. The manual
states that LAMMPS should output the extensive value for 'sum' mode, i.e. #4 =
-(#5+#6+#7)/(3*vol), but this is not true in the output. Here is the link to
this documentation page:
http://lammps.sandia.gov/doc/compute_reduce.html
Here is the thermo output: