Monte Carlo simulations within the grand canonical ensemble are used to explore the liquid-vapour coexistence curve and critical point properties of the Lennard-Jones fluid. Attention is focused on the joint distribution of density and energy fluctuations at coexistence. In the vicinity of the critical point, this distribution is analysed using mixed-field finite-size scaling techniques aided by histogram reweighting methods. The analysis yields highly accurate estimates of the critical point parameters, as well as exposing the size and character of corrections to scaling. In the sub-critical coexistence region the density distribution is obtained by combining multicanonical simulations with histogram reweighting techniques. It is demonstrated that this procedure permits an efficient and accurate mapping of the coexistence curve, even deep within the two phase region.

Fluids in which the interparticle potential has a hard core, is attractive at moderate separations, and repulsive at greater separations are known to exhibit novel phase behavior, including stable inhomogeneous phases. Here we report a joint simulation and theoretical study of such a fluid, focusing on the relationship between the liquid-vapor transition line and any new phases. The phase diagram is studied as a function of the amplitude of the attraction for a certain fixed amplitude of the long ranged repulsion. We find that the effect of the repulsion is to substitute the liquidvapor critical point and a portion of the associated liquid-vapor transition line, by two first order transitions. One of these transitions separates the vapor from a fluid of spherical liquidlike clusters; the other separates the liquid from a fluid of spherical voids. At low temperature, the two transition lines intersect one another and a vapor-liquid transition line at a triple point. While most integral equation theories are unable to describe the new phase transitions, the Percus Yevick approximation does succeed in capturing the vapor-cluster transition, as well as aspects of the structure of the cluster fluid, in reasonable agreement with the simulation results.

No abstract

AbsiraeiWe develop a finite-suemaling theory describing the joint density and energy fluctuations in a near-mitical fluid. As a result of the mixing of the lemperalure and chemical potential in the WO relevant scaling fields, the energy operator features in the crilical density distribution as an antisymmetric mrrestion lo the Limiting ale-invariant form Both the limiling form and the correction am predicted to be fuaetions that are characteristic of lhe king universality dass, and am independently knom. The theory is tested with extensive Monte Carlo studies of lhe lvmdimensional LennardJones fluid, within the grand canonical ensemble. The simulations and sgling framework together are shown to prwide a powerful way of identifying lhe loCalion of lhe liquid-gas aitical point, while confirming and darimng its esenlially king charader. 'Ihe simulations also show a clearly identifiable signature of the rield-mixing responsible for the hilure of ihe law of rectilinear dismeler.

We present a method for the direct evaluation of the difference between the free energies of two crystalline structures, of different symmetry. The method rests on a Monte Carlo procedure which allows one to sample along a path, through atomic-displacement-space, leading from one structure to the other by way of an intervening transformation that switches one set of lattice vectors for another. The configurations of both structures can thus be sampled within a single Monte Carlo process, and the difference between their free energies evaluated directly from the ratio of the measured probabilities of each. The method is used to determine the difference between the free energies of the fcc and hcp crystalline phases of a system of hard spheres.