Modifications of the SP-MC method for the computer simulation of chemical potentials: ternary mixtures of fused hard sphere fluids

Analysis of probability of inserting a hard spherical particle with small diameter in hard-sphere fluid

The probability of inserting, without overlap, a hard spherical particle of diameter σ in a hard-sphere fluid of diameter σ0 and packing fraction η determines its excess chemical potential at infinite dilution, μex(σ, η). In our previous work [R. L. Davidchack and B. B. Laird, J. Chem. Phys. 157, 074701 (2022)], we used Widom's particle insertion method within molecular dynamics simulations to obtain high precision results for μex(σ, η) with σ/σ0 ≤ 4 and η ≤ 0.5. In the current work, we investigate the behavior of this quantity at small σ. In particular, using the inclusion-exclusion principle, we relate the insertion probability to the hard-sphere fluid distribution functions and thus derive the higher-order terms in the Taylor expansion of μex(σ, η) at σ = 0. We also use direct evaluation of the excluded volume for pairs and triplets of hard spheres to obtain simulation results for μex(σ, η) at σ/σ0 ≤ 0.2247 that are of much higher precision than those obtained earlier with Widom's method. These results allow us to improve the quality of the small-σ correction in the empirical expression for μex(σ, η) presented in our previous work.

Chemical potential and surface free energy of a hard spherical particle in hard-sphere fluid over the full range of particle diameters

The excess chemical potential $\mu^\mathrm{ex}(\sigma,\eta)$ of a test hard spherical particle of diameter $\sigma$ in a fluid of hard spheres of diameter $\sigma_0$ and packing fraction $\eta$ can be computed with high precision using Widom's particle insertion method [J.~Chem.~Phys.~{\bf 39}, 2808 (1963)] for $\sigma$ between 0 and just larger than 1 and/or small $\eta$. Heyes and Santos [J.~Chem.~Phys.~{\bf 145}, 214504 (2016)] showed analytically that the only polynomial representation of $\mu^\mathrm{ex}$ consistent with the limits of $\sigma$ at zero and infinity has a cubic form. On the other hand, through the solvation free energy relationship between $\mu^\mathrm{ex}$ and the surface free energy $\gamma$ of hard-sphere fluid at a hard spherical wall, we can obtain precise measurements of $\mu^\mathrm{ex}$ for large $\sigma$, extending up to infinity (flat wall) [J.~Chem. Phys.~{\bf 149}, 174706 (2018)]. Within this approach, the cubic polynomial representation is consistent with the assumptions of Morphometric Thermodynamics. In this work, we present measurements of $\mu^\mathrm{ex}$ that combine the two methods to obtain high-precision results for the full range of $\sigma$ values from zero to infinity, which show statistically significant deviations from the cubic polynomial form. We propose an empirical functional form for $\mu^\mathrm{ex}$ dependence on $\sigma$ and $\eta$ which better fits the measurement data while remaining consistent with the analytical limiting behaviour at zero and infinite $\sigma$.

Field induced gradient simulations: a high throughput method for computing chemical potentials in multicomponent systems

We present a simulation method for direct computation of chemical potentials in multicomponent systems. The method involves application of a field to generate spatial gradients in the species number densities at equilibrium, from which the chemical potential of each species is theoretically estimated. A single simulation yields results over a range of thermodynamic states, as in high throughput experiments, and the method remains computationally efficient even at high number densities since it does not involve particle insertion at high densities. We illustrate the method by Monte Carlo simulations of binary hard sphere mixtures of particles with different sizes in a gravitational field. The results of the gradient Monte Carlo method are found to be in good agreement with chemical potentials computed using the classical Widom particle insertion method for spatially uniform systems.