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.