An extension of the NpT plus test particle method for the determination of the vapour-liquid equilibria of pure fluids

Predicting structural properties of fluids by thermodynamic extrapolation

We describe a methodology for extrapolating the structural properties of multicomponent fluids from one thermodynamic state to another. These properties generally include features of a system that may be computed from an individual configuration such as radial distribution functions, cluster size distributions, or a polymer's radius of gyration. This approach is based on the principle of using fluctuations in a system's extensive thermodynamic variables, such as energy, to construct an appropriate Taylor series expansion for these structural properties in terms of intensive conjugate variables, such as temperature. Thus, one may extrapolate these properties from one state to another when the series is truncated to some finite order. We demonstrate this extrapolation for simple and coarse-grained fluids in both the canonical and grand canonical ensembles, in terms of both temperatures and the chemical potentials of different components. The results show that this method is able to reasonably approximate structural properties of such fluids over a broad range of conditions. Consequently, this methodology may be employed to increase the computational efficiency of molecular simulations used to measure the structural properties of certain fluid systems, especially those used in high-throughput or data-driven investigations.

Communication: Predicting virial coefficients and alchemical transformations by extrapolating Mayer-sampling Monte Carlo simulations

Virial coefficients are predicted over a large range of both temperatures and model parameter values (i.e., alchemical transformation) from an individual Mayer-sampling Monte Carlo simulation by statistical mechanical extrapolation with minimal increase in computational cost. With this extrapolation method, a Mayer-sampling Monte Carlo simulation of the SPC/E (extended simple point charge) water model quantitatively predicted the second virial coefficient as a continuous function spanning over four orders of magnitude in value and over three orders of magnitude in temperature with less than a 2% deviation. In addition, the same simulation predicted the second virial coefficient if the site charges were scaled by a constant factor, from an increase of 40% down to zero charge. This method is also shown to perform well for the third virial coefficient and the exponential parameter for a Lennard-Jones fluid.

Magnetization and susceptibility of ferrofluids

A second-order Taylor series expansion of the free energy functional provides analytical expressions for the magnetic field dependence of the free energy and of the magnetization of ferrofluids, here modeled by dipolar Yukawa interaction potentials. The corresponding hard core dipolar Yukawa reference fluid is studied within the framework of the mean spherical approximation. Our findings for the magnetic and phase equilibrium properties are in quantitative agreement with previously published and new Monte Carlo simulation data.

Computing the starting state for Gibbs-Duhem integration

Gibbs-Duhem integration implies the numerical integration of a Clapeyron equation. To start the numerical integration, an initial coexistence point and a corresponding initial slope of the Clapeyron equation are needed. In order to apply Gibbs-Duhem integration to all kinds of systems at diverse physical conditions, one has to investigate and assess the available methods that can be used to compute these initial values. This publication focuses on vapor-liquid equilibria in binary mixtures comprising chain molecules. The initial coexistence point is either computed with the NVbeta Gibbs ensemble or with the Npbeta+test molecule method with overlapping distributions, which is introduced in this publication. Although computationally demanding, the Npbeta+test molecule method with overlapping distributions is applicable at conditions where the NVbeta Gibbs ensemble fails. We investigated three methods that can be employed to compute the initial slope of the Clapeyron equation. The Widom method and the overlapping-distributions difference method provide correct values for the initial slope. The difference method does only provide the correct answer in special cases. The possibility to judge the reliability of the results makes the overlapping-distributions difference method the safest route to the initial slope. Gibbs-Duhem integration requires the frequent computation of the slope of the Clapeyron equation. This slope depends on ensemble averages of the composition. A new bias method for efficient sampling of the composition in a semigrand-canonical simulation of chain molecules is presented. This bias method considerably enhances the composition sampling in systems comprising chain molecules of different sizes.