Remarks:
Uncertainties were obtained
from 5 independent simulations and represent 67% confidence limits (one
standard deviation).
Liquid-vapor coexistence was determined by adjusting the activity such
that the pressures of the liquid and vapor phases were
equal. Following, identification of state points below the coexistence pressure was accomplished by adjusting the activity until the pressure of the vapor phase was equal to the desired fraction of psat[6]. 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 = Λ-3 exp(βμ)
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.
Phase-coexistence energies were obtained by determining the
mean potential energy at a given value of N for an additional 8 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, linear force shifted at 2.5σ, the critical properties are Tc*=0.937, ρc*=0.320, and pc*=0.0820 [7].
[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).
[4] J. R. Errington and A. Z. Panagiotopoulos,
J. Chem. Phys. 109, 1093
(1998).
[5] A. Z. Panagiotopoulos, J. Phys.: Condens.
Matter 12, 25 (2000).
[6] J. R. Errington and V. K. Shen, J. Chem. Phys. 123, 164103 (2005).
[7] J. R. Errington P. G. Debenedetti, and S. Torquato, J. Chem. Phys. 118, 2256 (2003).
[8] V. K. Shen and D. W. Siderius, J. Chem. Phys., 140, 244106, (2014).
[9] V. K. Shen and J. R. Errington, J. Phys. Chem. B 108, 19595, (2004).
[10] V. K. Shen and J. R. Errington, J. Chem. Phys. 122, 064508, (2005).
[11] V. K. Shen, R. D. Mountain, and J. R. Errington, J. Phys. Chem. B 111, 6198, (2007).
[12] D. W. Siderius and V. K. Shen, J. Phys. Chem. 117, 5681, (2013).
|