Benchmark results for the LennardJones fluidThe purpose of these pages is to provide some explicit results from molecular dynamics and Monte Carlo simulations for the LennardJones fluid. It is intended to provide guides for testing codes. Reproducing these results is a test of the correctness of codes, either written by the user or obtained elsewhere. The explicit conditions for each of the sets of results are supplied so that meaningful comparisons of your results with the ones listed here are possible. The information presented here has been organized into eight different pages. It is available as web pages or as a single excel file.
As is usually the case, temperature, density (number density), pressure, etc., are given in reduced units (denoted by *). That is, these properties are expressed in terms of the characteristic energy, ε, and length scale, σ, defined by the LennardJones potential: Therefore, the reduced temperature T*, density ρ*, and pressure p* are k_{B}T/ε, ρσ^{3}, pσ^{3}/ε, respectively. The critical parameters for the pure LJ fluid have been determined to be T_{c}* = 1.3120(7), ρ_{c}* = 0.316(1), and p_{c}* = 0.1279(6) [3]. For computational expediency, the potential and force are usually truncated at some cutoff distance r_{c}. That is, the effective potential V(r) is
In this work, the cutoff distance is taken to be 3σ, unless noted otherwise. An approximate correction for this truncation to the energy and pressure of the system can be obtained by assuming the spatial correlations beyond the cutoff distance are unity. The reader is referred to Refs. [4] and [5] for these socalled “standard long range corrections” (sLRC). It should be noted that this not a good assumption in inhomogeneous fluids. For completeness, the working expressions for obtaining the long range corrections to the potential energy per particle and virial pressure are provided below: where ρ is the bulk number density, V_{LJ}(r) is the LennardJones potential energy, and r_{c} is the truncation (cutoff) distance. A subtle point to note is that the above pressure correction does not account for impulsive effects at r_{c} where the potential energy (and therefore the force) changes discontinuously. This is particularly important when trying to simulate directly interfacial phenomena [6]. Finally, it should be mentioned that other truncation and long range correction schemes for dealing with nonbonded interactions exist [48]. In some cases, the truncation schemes can yield critical properties that differ from those quoted above [8]. In future versions, coordinate sets will be available for which the potential energy per particle and the virial are specified. This information can be used to test energy/force routines. References1. J. R. Errington, J. Chem. Phys. 118, 9915 (2003).2. J. K. Johnson, J. A. Zollweg, and K. E. Gubbins, Mol. Phys. 78, 591 (1993). 3. J. J. Potoff and A. Z. Panagiotopoulos, J. Chem. Phys. 109, 10914 (1998). 4. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, New York, 1989). 5. D. Frenkel and B. Smit, Understanding Molecular Simulation, 2nd ed. (Academic, San Diego, 2002)., pp.3738. 6. A. Trkohymchuk and J. Alejandre, J. Chem. Phys. 111, 8510 (1999). 7. D. N. Theodorou and U. W. Suter, J. Chem. Phys. 82, 955 (1985). 8. B. Smit, J. Chem. Phys. 96, 8639 (1992).
