Hard sphere compressibility factors for equation of state development

Ex Machina Determination of Structural Correlation Functions

Determining the structural properties of condensed-phase systems is a fundamental problem in theoretical statistical mechanics. Here we present a machine learning method that is able to predict structural correlation functions with significantly improved accuracy in comparison with traditional approaches. The usefulness of this ex machina (from the machine) approach is illustrated by predicting the radial distribution functions of two paradigmatic condensed-phase systems, a Lennard-Jones fluid and a hard-sphere fluid, and then comparing those results to the results obtained using both integral equation methods and empirically motivated analytical functions. We find that application of the developed ex machina method typically decreases the predictive error by more than an order of magnitude in comparison with traditional theoretical methods.

Equations of the state of hard sphere fluids based on recent accurate virial coefficients B5–B12

A review on the numerical virial coefficients, compressibility factor, fluid–solid phase transition point and equations of the state of hard sphere fluids.

Extension of the BMCSL equation of state for hard spheres to the metastable disordered region: Application to the SAFT approach

A simple modification of the Boublík-Mansoori-Carnahan-Starling-Leland equation of state is proposed for an application to the metastable disordered region. The new model has a positive pole at the jamming limit and can accurately describe the molecular simulation data of pure hard in the stable fluid region and along the metastable branch. The new model has also been applied to binary mixtures hard spheres, and an excellent description of the fluid and metastable branches can be obtained by adjusting the jamming packing fraction. The new model for hard sphere mixtures can be used as the repulsive term of equations of state for real fluids. In this case, the modified equations of state give very similar predictions of thermodynamic properties as the original models, and one can remove the multiple liquid density roots observed for some versions of the Statistical Associating Fluid Theory (SAFT) at low temperature without any modification of the dispersion term.

Pseudo hard-sphere potential for use in continuous molecular-dynamics simulation of spherical and chain molecules

We present a continuous pseudo-hard-sphere potential based on a cut-and-shifted Mie (generalized Lennard-Jones) potential with exponents (50, 49). Using this potential one can mimic the volumetric, structural, and dynamic properties of the discontinuous hard-sphere potential over the whole fluid range. The continuous pseudo potential has the advantage that it may be incorporated directly into off-the-shelf molecular-dynamics code, allowing the user to capitalise on existing hardware and software advances. Simulation results for the compressibility factor of the fluid and solid phases of our pseudo hard spheres are presented and compared both to the Carnahan-Starling equation of state of the fluid and published data, the differences being indistinguishable within simulation uncertainty. The specific form of the potential is employed to simulate flexible chains formed from these pseudo hard spheres at contact (pearl-necklace model) for m(c) = 4, 5, 7, 8, 16, 20, 100, 201, and 500 monomer segments. The compressibility factor of the chains per unit of monomer, m(c), approaches a limiting value at reasonably small values, m(c) < 50, as predicted by Wertheim's first order thermodynamic perturbation theory. Simulation results are also presented for highly asymmetric mixtures of pseudo hard spheres, with diameter ratios of 3:1, 5:1, 20:1 over the whole composition range.

High-order virial coefficients and equation of state for hard sphere and hard disk systems

A very simple and accurate approach is proposed to predict the high-order virial coefficients of hard spheres and hard disks. In the approach, the nth virial coefficient B(n) is expressed as the sum of n(D-1) and a remainder, where D is the spatial dimension of the system. When n > or = 3, the remainders of the virials can be accurately expressed with Padé-type functions of n. The maximum deviations of predicted B(5)-B(10) for the two systems are only 0.0209%-0.0044% and 0.0390%-0.0525%, respectively, which are much better than the numerous existing approaches. The virial equation based on the predicted virials diverges when packing fraction eta = 1. With the predicted virials, the compressibility factors of hard sphere system can be predicted very accurately in the whole stable fluid region, and those in the metastable fluid region can also be well predicted up to eta = 0.545. The compressibility factors of hard disk fluid can be predicted very accurately up to eta = 0.63. The simulated B(7) and B(10) for hard spheres are found to be inconsistent with the other known virials and therefore they are modified as 53.2467 and 105.042, respectively.

Alternative fundamental measure theory for additive hard sphere mixtures

The purpose of this short paper is to present an alternative fundamental measure theory (FMT) for hard sphere mixtures. Keeping the main features of the original Rosenfeld's FMT [Phys. Rev. Lett. 63, 980 (1989)] and using the dimensional and the low-density limit conditions a new functional is derived incorporating Boublik's multicomponent extension [Mol. Phys. 59, 371 (1986)] of highly accurate Kolafa's equation of state for pure hard spheres. We test the theory for pure hard spheres and hard sphere mixtures near a planar hard wall and compare the results with the original Rosenfeld's FMT and one of its modifications and with new very accurate simulation data. The test reveals an excellent agreement between the results based on the alternative FMT and simulation data for density profile near a contact and some improvement over the original Rosenfeld's FMT and its modification at the contact region.

Influence of two-body and three-body interatomic forces on gas, liquid, and solid phases

Accurate molecular dynamics simulations are reported which quantify the contributions of two- and three-body interactions in the gas, liquid, and solid phases of argon at both subcritical and supercritical conditions. The calculations use an accurate two-body potential in addition to contributions from three-body dispersion interactions from third-order triple-dipole interactions. The number dependence of three-body interactions is quantified, indicating that a system size of at least five hundred atoms is required for reliable calculations. The results indicate that, although the contribution of three-body interaction to the overall energy is small, three-body interactions significantly affect the pressure at which vapor-liquid and solid-liquid transitions are observed. In particular, three-body interactions substantially increase the pressure of the freezing point. Unlike two-body interactions, which vary with both density and temperature, for a given density, three-body interactions have a near-constant 'background' value irrespective of the temperature. Both two-body interactions and kinetic energy have an important role in vapor-liquid equilibria whereas solid-liquid equilibria are dominated by two-body interactions.

Nonanalytical equation of state of the hard sphere fluid

A model based on classical nucleation theory is proposed to describe phase behavior in stable and metastable regions near a first-order phase transition. The resulting equation of state is not an analytical function at the phase transition point. The model is tested on the hard sphere fluid where it is combined with the virial expansion at low densities. This nonanalytical equation of state is able to capture the observed "anomalous increase" of pressure at high densities. Consequences to description of phase equilibria by classical equations of state are discussed.

Shape factors in equations of state : Part II. Repulsion phenomena in multicomponent chain fluids

An equation of state for the multicomponent fluid phase of nonattracting rigid particles of arbitrary shape is presented. The equation is a generalization of a previously presented equation of state for pure fluids of rigid particles; the approach describes the volumetric properties of a pure fluid in terms of a shape factor, zeta, which can be back calculated by scaling the volumetric properties of pure fluids to that of a hard sphere. The performance of the proposed equation is tested against mixtures of chain fluids immersed in a "monomeric" solvent of hard spheres of equal and different sizes. Extensive new Monte Carlo simulation data are presented for 19 binary mixtures of hard homonuclear tangent freely-jointed hard sphere chains (pearl-necklace) of various lengths (three to five segments), with spheres of several size ratios and at various compositions. The performance of the proposed equation is compared to the hard-sphere SAFT approach and found to be of comparable accuracy. The equation proposed is further tested for mixtures of spheres with spherocylinders. In all cases, the equation proved to be accurate and simple to use.

Thermodynamic properties of lattice hard-sphere models

Thermodynamic properties of several lattice hard-sphere models were obtained from grand canonical histogram- reweighting Monte Carlo simulations. Sphere centers occupy positions on a simple cubic lattice of unit spacing and exclude neighboring sites up to a distance sigma. The nearestneighbor exclusion model, sigma = radical2, was previously found to have a second-order transition. Models with integer values of sigma = 1 or 2 do not have any transitions. Models with sigma = radical3 and sigma = 3 have weak first-order fluid-solid transitions while those with sigma = 2 radical2, 2 radical3, and 3 radical2 have strong fluid-solid transitions. Pressure, chemical potential, and density are reported for all models and compared to the results for the continuum, theoretical predictions, and prior simulations when available.