Simulation of phase boundaries using constrained cell models.

Despite impressive advances, precise simulation of fluid-fluid and fluid-solid phase transitions still remains a challenging task. The present work focuses on the determination of the phase diagram of a system of particles that interact through a pair potential, φ(r), which is of the form φ(r) = 4ε[(σ/r)(2n) - (σ/r)(n)] with n = 12. The vapor-liquid phase diagram of this model is established from constant-pressure simulations and flat-histogram techniques. The properties of the solid phase are obtained from constant-pressure simulations using constrained cell models. In the constrained cell model, the simulation volume is divided into Wigner-Seitz cells and each particle is confined to moving in a single cell. The constrained cell model is a limiting case of a more general cell model which is constructed by adding a homogeneous external field that controls the relative stability of the fluid and the solid phase. Fluid-solid coexistence at a reduced temperature of 2 is established from constant-pressure simulations of the generalized cell model. The previous fluid-solid coexistence point is used as a reference point in the determination of the fluid-solid phase boundary through a thermodynamic integration type of technique based on histogram reweighting. Since the attractive interaction is of short range, the vapor-liquid transition is metastable against crystallization. In the present work, the phase diagram of the corresponding constrained cell model is also determined. The latter is found to contain a stable vapor-liquid critical point and a triple point.

Simulation of fluid-solid coexistence via thermodynamic integration using a modified cell model.

Despite recent advances, precise simulation of fluid-solid transitions still remains a challenging task. Thermodynamic integration techniques are the simplest methods to study fluid-solid coexistence. These methods are based on the calculation of the free energies of the fluid and the solid phases, starting from a state of known free energy which is usually an ideal-gas state. Despite their simplicity, the main drawback of thermodynamic integration techniques is the large number of states that must be simulated. In the present work, a thermodynamic integration technique, which reduces the number of simulated states, is proposed and tested on a system of particles interacting via an inverse twelfth-power potential energy function. The simulations are implemented at constant pressure and the solid phase is modeled according to the constrained cell model of Hoover and Ree. The fluid and the solid phases are linked together by performing constant-pressure simulations of a modified cell model. The modified cell model, which was originally proposed by Hoover and Ree, facilitates transitions between the fluid and the solid phase by tuning a homogeneous external field. This model is simulated on a constant-pressure path for a series of progressively increasing values of the field, thus allowing for direct determination of the free energy difference between the fluid and the solid phase via histogram reweighting. The size-dependent results are analyzed using histogram reweighting and finite-size scaling techniques. The scaling analysis is based on studying the size-dependent behavior of the second- and higher-order derivatives of the free energy as well as the dimensionless moment ratios of the order parameter. The results clearly demonstrate the importance of accounting for size effects in simulation studies of fluid-solid transitions.

Precise simulation of the freezing transition of supercritical Lennard-Jones.

The fluid-solid transition of the Lennard-Jones model is analyzed along a supercritical isotherm. The analysis is implemented via a simulation method which is based on a modification of the constrained cell model of Hoover and Ree. In the context of hard-sphere freezing, Hoover and Ree simulated the solid phase using a constrained cell model in which each particle is confined within its own Wigner-Seitz cell. Hoover and Ree also proposed a modified cell model by considering the effect of an external field of variable strength. High-field values favor configurations with a single particle per Wigner-Seitz cell and thus stabilize the solid phase. In previous work, a simulation method for freezing transitions, based on constant-pressure simulations of the modified cell model, was developed and tested on a system of hard spheres. In the present work, this method is used to determine the freezing transition of a Lennard-Jones model system on a supercritical isotherm at a reduced temperature of 2. As in the case of hard spheres, constant-pressure simulations of the fully occupied constrained cell model of a system of Lennard-Jones particles indicate a point of mechanical instability at a density which is approximately 70% of the density at close packing. Furthermore, constant-pressure simulations of the modified cell model indicate that as the strength of the field is reduced, the transition from the solid to the fluid is continuous below the mechanical instability point and discontinuous above. The fluid-solid transition of the Lennard-Jones system is obtained by analyzing the field-induced fluid-solid transition of the modified cell model in the high-pressure, zero-field limit. The simulations are implemented under constant pressure using tempering and histogram reweighting techniques. The coexistence pressure and densities are determined through finite-size scaling techniques for first-order phase transitions which are based on analyzing the size-dependent behavior of susceptibilities and dimensionless moment ratios of the order parameter.

A Monte Carlo study of the freezing transition of hard spheres.

A simulation method for fluid-solid transitions, which is based on a modification of the constrained cell model of Hoover and Ree, is developed and tested on a system of hard spheres. In the fully occupied constrained cell model, each particle is confined in its own Wigner-Seitz cell. Constant-pressure simulations of the constrained cell model for a system of hard spheres indicate a point of mechanical instability at a density which is about 64% of the density at the close packed limit. Below that point, the solid is mechanically unstable since without the confinement imposed by the cell walls it will disintegrate to a disordered, fluid-like phase. Hoover and Ree proposed a modified cell model by introducing an external field of variable strength. High values of the external field variable favor configurations with one particle per cell and thus stabilize the solid phase. In this work, the modified cell model of a hard-sphere system is simulated under constant-pressure conditions using tempering and histogram reweighting techniques. The simulations indicate that as the strength of the field is reduced, the transition from the solid to the fluid phase is continuous below the mechanical instability point and discontinuous above. The fluid-solid transition of the hard-sphere system is determined by analyzing the field-induced fluid-solid transition of the modified cell model in the limit in which the external field vanishes. The coexistence pressure and densities are obtained through finite-size scaling techniques and are in good accord with previous estimates.