Communication: Direct determination of triple-point coexistence through cell model simulation.

In simulations of fluid-solid coexistence, the solid phase is modeled as a constrained system of Wigner-Seitz cells with one particle per cell. This model, commonly referred to as the constrained cell model, is a limiting case of a more general cell model, which is formed by considering a homogeneous external field that controls the number of particles per cell and, hence, the relative stability of the solid against the fluid phase. The generalized cell model provides a link that connects the disordered, fluid phase with the ordered, solid phase. In the present work, the phase diagram of this model is investigated through multicanonical simulations at constant pressure and histogram reweighting techniques for a system of 256 Lennard-Jones particles. The simulation data are used to obtain an estimate of the triple point of the Lennard-Jones system. The triple-point pressure is found to be higher compared to previous work. The likely explanation for this discrepancy is the highly compressible nature of the gas phase.

[1]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[2]  Jean-Pierre Hansen,et al.  Phase Transitions of the Lennard-Jones System , 1969 .

[3]  G. Orkoulas,et al.  Simulation of fluid–solid coexistence via thermodynamic integration using a modified cell model , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[4]  W. G. Hoover,et al.  Use of Computer Experiments to Locate the Melting Transition and Calculate the Entropy in the Solid Phase , 1967 .

[5]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[6]  J. Ilja Siepmann,et al.  Direct Gibbs Ensemble Monte Carlo Simulations for Solid−Vapor Phase Equilibria: Applications to Lennard−Jonesium and Carbon Dioxide , 2001 .

[7]  G. Orkoulas,et al.  Communication: A simple method for simulation of freezing transitions. , 2011, The Journal of chemical physics.

[8]  G. Orkoulas,et al.  Precise simulation of subcritical freezing using constrained cell models , 2012 .

[9]  Berend Smit,et al.  Chapter 3 – Monte Carlo Simulations , 2002 .

[10]  G. Orkoulas,et al.  Communication: phase transitions, criticality, and three-phase coexistence in constrained cell models. , 2012, The Journal of chemical physics.

[11]  G. Orkoulas,et al.  Precise simulation of the freezing transition of supercritical Lennard-Jones. , 2011, The Journal of chemical physics.

[12]  Borgs,et al.  Finite-size effects at asymmetric first-order phase transitions. , 1992, Physical review letters.

[13]  A. Ladd,et al.  Interfacial and co-existence properties of the Lennard-Jones system at the triple point , 1978 .

[14]  William G. Hoover,et al.  Melting Transition and Communal Entropy for Hard Spheres , 1968 .

[15]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[16]  D. Kofke,et al.  Thermodynamic and structural properties of model systems at solid-fluid coexistence: I. Fcc and bcc soft spheres , 1995 .

[17]  Michael Nayhouse,et al.  Simulation of phase boundaries using constrained cell models , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[18]  G. Orkoulas,et al.  A Monte Carlo study of the freezing transition of hard spheres , 2011, Journal of physics. Condensed matter : an Institute of Physics journal.