Skip to main content
U.S. flag

An official website of the United States government

Official websites use .gov
A .gov website belongs to an official government organization in the United States.

Secure .gov websites use HTTPS
A lock ( ) or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.

SAT-TMMC: Liquid-vapor coexistence properties - Cut Potential at 5.0σ

Liquid-vapor coexistence properties obtained by grand-canonical transition-matrix Monte Carlo and histogram re-weighting over the reduced temperature range 0.60 to 1.25 at increments of 0.05. Mean values of the saturation pressure, density, potential energy per molecule, and activity (chemical potential- see below) for each phase are reported.

METHOD Grand-canonical transition-matrix Monte Carlo and histogram re-weighting [1, 8-12]
V/σ3 1000
TRUNCATION Cut potential (no tail correction) at 5.0σ
Prob. of Disp. Move 0.4
Prob. of Ins/Del Move 0.6
Biasing Function Update Frequency 1.0E6 trial moves
Simulation Length 2.0E10 trial moves















0.60 4.765E-04 6.447E-07 8.812E-01 1.015E-04 2.841E-04 3.998E-07 -6.732E-03 2.377E-05 -6.395E+00 1.417E-02 -7.661E+00 1.126E-03
0.65 1.084E-03 8.909E-07 8.609E-01 2.245E-04 6.959E-04 4.147E-07 -1.140E-02 1.374E-05 -6.218E+00 7.299E-03 -6.852E+00 6.276E-04
0.70 2.164E-03 6.955E-07 8.398E-01 3.081E-04 1.483E-03 5.282E-07 -2.626E-02 1.010E-05 -6.026E+00 3.823E-03 -6.178E+00 3.326E-04
0.73 3.120E-03 9.424E-07 8.268E-01 7.671E-05 2.213E-03 5.363E-07 -3.652E-02 3.492E-05 -5.907E+00 2.267E-03 -5.827E+00 2.725E-04
0.75 3.913E-03 8.244E-07 8.183E-01 1.304E-04 2.835E-03 7.779E-07 -4.480E-02 2.422E-05 -5.831E+00 1.846E-03 -5.611E+00 2.480E-04
0.80 6.548E-03 2.512E-06 7.958E-01 9.656E-05 4.972E-03 1.647E-06 -7.130E-02 5.459E-05 -5.637E+00 2.798E-03 -5.130E+00 2.658E-04
0.85 1.031E-02 4.335E-06 7.727E-01 8.343E-05 8.113E-03 3.008E-06 -1.075E-01 6.212E-05 -5.442E+00 2.668E-03 -4.718E+00 2.648E-04
0.90 1.548E-02 6.535E-06 7.485E-01 1.332E-04 1.257E-02 3.913E-06 -1.553E-01 7.623E-05 -5.239E+00 1.931E-03 -4.364E+00 1.557E-04
0.95 2.239E-02 1.442E-06 7.231E-01 1.541E-04 1.854E-02 7.579E-07 -2.170E-01 3.475E-05 -5.033E+00 1.610E-03 -4.056E+00 5.940E-05
1.00 3.145E-02 7.132E-06 6.958E-01 7.023E-05 2.632E-02 2.361E-06 -2.953E-01 4.548E-05 -4.815E+00 1.395E-03 -3.787E+00 1.120E-04
1.05 4.328E-02 6.936E-06 6.663E-01 1.115E-04 3.616E-02 4.671E-06 -3.946E-01 9.746E-05 -4.588E+00 6.876E-04 -3.551E+00 5.547E-05
1.10 5.873E-02 1.656E-05 6.339E-01 1.576E-04 4.834E-02 8.428E-06 -5.208E-01 1.598E-04 -4.344E+00 7.039E-04 -3.343E+00 1.107E-04
1.15 7.928E-02 3.461E-05 5.970E-01 9.314E-05 6.314E-02 1.070E-05 -6.842E-01 3.649E-04 -4.075E+00 6.823E-04 -3.158E+00 5.119E-05
1.20 1.080E-01 4.722E-05 5.527E-01 1.207E-04 8.091E-02 1.860E-05 -9.068E-01 4.098E-04 -3.768E+00 9.739E-04 -2.994E+00 1.181E-04
1.25 1.553E-01 1.144E-04 4.906E-01 1.891E-04 1.021E-01 1.420E-05 -1.276E+00 1.204E-03 -3.368E+00 1.146E-03 -2.848E+00 4.959E-05






Uncertainties were obtained from five independent simulations and represent 95% confidence limits based on a standard t statistic. Liquid-vapor coexistence was determined by adjusting the activity such that the pressures of the liquid and vapor phases were equal. Here, the pressure is not the conventional virial pressure [2,3] but is the actual thermodynamic pressure, based on the fact that the absolute free energies can be obtained from the distributions determined from simulation [4]. Alternative methods, for example Gibbs-ensemble Monte Carlo and combination grand-canonical Monte Carlo and histogram re-weighting, can be used to determine liquid-vapor coexistence. A review of standard methods of phase equilibria simulations can be found in Ref. 5.

As introduced in Refs. 2 and 3, the activity, z, is defined as

$$ z = \dfrac{ \exp\left( \beta \mu \right)}{\lambda^3}$$

where Λ is the de Broglie wavelength, β = 1/(kBT) (where kB is Boltzmann's constant), and μ is the chemical potential. It is sometimes more convenient to work with ln z in the simulations and in post-processing. (NOTE: The reported activity is dimensionless, having been scaled by the LJ length cubed.)

Phase-coexistence energies were obtained by determining the mean potential energy at a given value of N for an additional 40 billion MC trials. Combining this information with the particle number probability distribution, the mean potential energy of the coexisting phases can be calculated [6].

For the Lennard-Jones fluid, cut potential at 5.σ, the critical properties were estimated to be Tc*=1.284, ρc*=0.318, and pc*=0.118. Estimates were found via rectilinear diameter analysis of TMMC data computed with V*=1000 close to the critical point [7]. (Finite-size scaling analysis has not been completed, so these critical properties should be taken simply as estimates.)


  1. J. R. Errington, J. Chem. Phys. 118, 9915 (2003).
  2. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, New York, 1989).
  3. D. Frenkel and B. Smit, Understanding Molecular Simulation, 2nd ed. (Academic, San Diego, 2002)., pp.37-38.
  4. J. R. Errington and A. Z. Panagiotopoulos, J. Chem. Phys., 109, 1093 (1998).
  5. A. Z. Panagiotopoulos, J. Phys.: Condens. Matter, 12, R25-R52, (2000).
  6. J. R. Errington and V. K. Shen, J. Chem. Phys., 123, 164103 (2005).
  7. B. Smit and C. P. Williams, J. Phys.: Condens. Matter, 2, 4281-4288 (1990).
  8. V. K. Shen and D. W. Siderius, J. Chem. Phys., 140, 244106, 2014.
  9. V. K. Shen and J. R. Errington, JPC B 108, 19595, 2004.
  10. V. K. Shen and J. R. Errington, JCP 122, 064508, 2005.
  11. V. K. Shen, R. D. Mountain, and J. R. Errington, JPC B 111, 6198, 2007.
  12. D. W. Siderius and V. K. Shen, JPC C 117, 5681, 2013.
Created September 11, 2015, Updated March 4, 2022