Structure of hard-sphere metastable fluids

Structural and thermodynamic properties of hard-sphere fluids

This Perspective article provides an overview of some of our analytical approaches to the computation of the structural and thermodynamic properties of single-component and multicomponent hard-sphere fluids. For the structural properties, they yield a thermodynamically consistent formulation, thus improving and extending the known analytical results of the Percus-Yevick theory. Approximate expressions linking the equation of state of the single-component fluid to the one of the multicomponent mixtures are also discussed.

Structural properties of additive binary hard-sphere mixtures

An approach to obtain the structural properties of additive binary hard-sphere mixtures is presented. Such an approach, which is a nontrivial generalization of the one recently used for monocomponent hard-sphere fluids [S. Pieprzyk, A. C. Brańka, and D. M. Heyes, Phys. Rev. E 95, 062104 (2017)2470-004510.1103/PhysRevE.95.062104], combines accurate molecular-dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation. A comparison of the direct correlation functions obtained with the present scheme with those derived from theoretical results stemming from the Percus-Yevick (PY) closure and the so-called rational-function approximation (RFA) is performed. The density dependence of the leading poles of the Fourier transforms of the total correlation functions and the decay of the pair correlation functions of the mixtures are also addressed and compared to the predictions of the two theoretical approximations. A very good overall agreement between the results of the present scheme and those of the RFA is found, thus suggesting that the latter (which is an improvement over the PY approximation) can safely be used to predict reasonably well the long-range behavior, including the structural crossover, of the correlation functions of additive binary hard-sphere mixtures.

Radial distribution function for hard spheres in fractal dimensions: A heuristic approximation

Analytic approximations for the radial distribution function, the structure factor, and the equation of state of hard-core fluids in fractal dimension d (1≤d≤3) are developed as heuristic interpolations from the knowledge of the exact and Percus-Yevick results for the hard-rod and hard-sphere fluids, respectively. In order to assess their value, such approximate results are compared with those of recent Monte Carlo simulations and numerical solutions of the Percus-Yevick equation for a fractal dimension [M. Heinen et al., Phys. Rev. Lett. 115, 097801 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.097801], a good agreement being observed.

Equation of state and critical point behavior of hard-core double-Yukawa fluids

A theoretical study on the equation of state and the critical point behavior of hard-core double-Yukawa fluids is presented. Thermodynamic perturbation theory, restricted to first order in the inverse temperature and having the hard-sphere fluid as the reference system, is used to derive a relatively simple analytical equation of state of hard-core multi-Yukawa fluids. Using such an equation of state, the compressibility factor and phase behavior of six representative hard-core double-Yukawa fluids are examined and compared with available simulation results. The effect of varying the parameters of the hard-core double-Yukawa intermolecular potential on the location of the critical point is also analyzed using different perspectives. The relevance of this analysis for fluids whose molecules interact with realistic potentials is also pointed out.

Multicomponent fluid of nonadditive hard spheres near a wall

A recently proposed rational-function approximation [Phys. Rev. E 84, 041201 (2011)] for the structural properties of nonadditive hard spheres is applied to evaluate analytically (in Laplace space) the local density profiles of multicomponent nonadditive hard-sphere mixtures near a planar nonadditive hard wall. The theory is assessed by comparison with NVT Monte Carlo simulations of binary mixtures with a size ratio 1:3 in three possible scenarios: a mixture with either positive or negative nonadditivity near an additive wall, an additive mixture with a nonadditive wall, and a nonadditive mixture with a nonadditive wall. It is observed that, while the theory tends to underestimate the local densities at contact (especially in the case of the big spheres) it captures very well the initial decay of the densities with increasing separation from the wall and the subsequent oscillations.

Use of the Ornstein-Zernike Percus-Yevick equation to extract interaction potentials from pair correlation functions

In this work, we test the ability of the Ornstein-Zernike equation in the Percus-Yevick approximation (OZPY) to generate interaction potentials from pair correlation functions (PCFs) of monatomic and diatomic Lennard-Jones fluids. The PCFs are generated by solving OZPY equation (monatomic fluid) and molecular-dynamics (MD) simulations (diatomic fluid). Since the interaction potentials are inputs in the OZPY method and the MD simulation, the extraction of the potential from the PCFs using OZPY is a test of self-consistency. This test is necessary if the procedure is to be used to generate coarse-grained (CG) potentials from PCFs. We find that the procedure is completely self-consistent for the monatomic fluid in the whole range of densities studied (reduced density up to 0.55, under reduced temperature of 2.0). In the diatomic case, we find that the procedure is generally self-consistent under both low and high densities, although there is a systematic deviation at high densities. The method is able to reproduce the two parameters (ε and σ) of the input Lennard-Jones potential model to within about 1%. This CG potential generating procedure can be straightforwardly extended to more complicated molecules.

Estimation of the solute diffusion coefficient of a dilute liquid alloy: static structure factor and isothermal compressibility estimates obtained using the rational function approximation of the radial distribution

A simple method for estimating the mass diffusion coefficient of a dilute binary liquid alloy that sequentially uses experimental data for the static structure factor and isothermal susceptibility of the solvent is presented, as well as another using the static structure factor alone and a method using the isothermal susceptibility alone. A fourth method that simultaneously uses the static structure factor and isothermal susceptibility is also noted. Of significance is the fact that these methods do not require information about the interatomic potential. Stability with respect to weights in the optimization process employed has been established and is reported, as well as some indication of the upper limits on the applicable solute concentration. Comparisons are made with results from a high quality capillary experiment for Pb 1 wt% Au liquid alloy performed in microgravity, and with velocity autocorrelation estimates derived from molecular dynamics simulation. The results suggest that the capillary experiments are influenced by reverse diffusion of the solvent, and actually measure an average of the mass diffusion coefficients, D(ij), weighted by the equilibrium concentrations of the solvent, x(1), and solute, x(2), defined by [Formula: see text] The three methods are required to provide upper and lower estimates for the mixed solvent-solute diffusion coefficient, which is not directly accessible from the experimental data, and demonstrate agreement with the experiment via D(tot).

Structure of hard-hypersphere fluids in odd dimensions

The structural properties of single component fluids of hard hyperspheres in odd space dimensionalities d are studied with an analytical approximation method that generalizes the rational function approximation earlier introduced in the study of hard-sphere fluids [S. B. Yuste and A. Santos, Phys. Rev. A 43, 5418 (1991)]. The theory makes use of the exact form of the radial distribution function to first order in density and extends it to finite density by assuming a rational form for a function defined in Laplace space, the coefficients being determined by simple physical requirements. Fourier transform in terms of reverse Bessel polynomials constitute the mathematical framework of this approximation, from which an analytical expression for the static structure factor is obtained. In its most elementary form, the method recovers the solution of the Percus-Yevick closure to the Ornstein-Zernike equation for hyperspheres at odd dimensions. The present formalism allows one to go beyond by yielding solutions with thermodynamic consistency between the virial and compressibility routes to any desired equation of state. Excellent agreement with available computer simulation data at d=5 and d=7 is obtained.

Thermodynamic properties of van der Waals fluids from Monte Carlo simulations and perturbative Monte Carlo theory

Computer simulations have been performed for fluids with van der Waals potential, that is, hard spheres with attractive inverse power tails, to determine the equation of state and the excess energy. On the other hand, the first- and second-order perturbative contributions to the energy and the zero- and first-order perturbative contributions to the compressibility factor have been determined too from Monte Carlo simulations performed on the reference hard-sphere system. The aim was to test the reliability of this "exact" perturbation theory. It has been found that the results obtained from the Monte Carlo perturbation theory for these two thermodynamic properties agree well with the direct Monte Carlo simulations. Moreover, it has been found that results from the Barker-Henderson [J. Chem. Phys. 47, 2856 (1967)] perturbation theory are in good agreement with those from the exact perturbation theory.

Volume terms for charged colloids: a grand-canonical treatment

We present a study of thermodynamic properties of suspensions of charged colloids on the basis of linear Poisson-Boltzmann theory. We calculate the effective Hamiltonian of the colloids by integrating out the ionic degrees of freedom grand canonically. This procedure not only yields the well-known pairwise screened-Coulomb interaction between the colloids, but also additional volume terms that affect the phase behavior and the thermodynamic properties, such as the osmotic pressure. These calculations are greatly facilitated by the grand-canonical character of our treatment of the ions and allow for relatively fast computations compared to earlier studies in the canonical ensemble. Moreover, the present derivation of the volume terms are relatively simple, make a direct connection with Donnan equilibrium, yield an explicit expression for the effective screening constant, and allow for extensions to include, for instance, nonlinear effects.

Simple analytic equations of state for hard-core single and double Yukawa fluids and mixtures based on second-order Barker-Henderson perturbation theory

A simple analytic expression with high precision for the radial distribution function of hard spheres is proposed. The form of the expression has been carefully selected to combine the well-known Camahan-Starling equation of state in it and satisfy the limit condition at low density, its simplicity and precision is superior to the well-known Percus-Yevick expression. The coefficients contained in the expression have been determined by fitting the Monte Carlo data for the first coordination shell, and by fitting both the Monte Carlo data and the numerical results of the Percus-Yevick expression for the second coordination shell. The expression has been applied to develop simple analytic equations of state for the hard-core single, double Yukawa fluids, and the hard-core Yukawa mixtures. The comparisons show that the agreement of our model with the computer simulation data is slightly better than the mean spherical approximation and other analytic models.

Generalized van der Waals theory for the thermodynamic properties of square-well fluids

A theory previously developed for the coordination number of square-well fluids is used within the context of a generalized van der Waals theory to obtain the compressibility factor and the internal energy of these fluids. Results are compared with computer simulations for several densities, temperatures, and potential widths, which are also reported.

Heterogeneities in supercooled liquids: a density-functional study

A metastable state, characterized by a low degree of mass localization, is identified using density-functional theory (DFT). This free energy minimum, located through the proper evaluation of competing terms in the free energy functional, is independent of the specific form of the DFT used. Computer simulation results on particle motion indicate that this heterogeneous state corresponds to the deeply supercooled state.