Structural properties of additive binary hard-sphere mixtures. III. Direct correlation functions.

An analysis of the direct correlation functions c_{ij}(r) of binary additive hard-sphere mixtures of diameters σ_{s} and σ_{b} (where the subscripts s and b refer to the "small" and "big" spheres, respectively), as obtained with the rational-function approximation method and the WM scheme introduced in previous work [S. Pieprzyk et al., Phys. Rev. E 101, 012117 (2020)2470-004510.1103/PhysRevE.101.012117], is performed. The results indicate that the functions c_{ss}(r<σ_{s}) and c_{bb}(r<σ_{b}) in both approaches are monotonic and can be well represented by a low-order polynomial, while the function c_{sb}(r<1/2(σ_{b}+σ_{s})) is not monotonic and exhibits a well-defined minimum near r=1/2(σ_{b}-σ_{s}), whose properties are studied in detail. Additionally, we show that the second derivative c_{sb}^{''}(r) presents a jump discontinuity at r=1/2(σ_{b}-σ_{s}) whose magnitude satisfies the same relationship with the contact values of the radial distribution function as in the Percus-Yevick theory.

Structural properties of additive binary hard-sphere mixtures. II. Asymptotic behavior and structural crossovers.

The structural properties of additive binary hard-sphere mixtures are addressed as a follow-up of a previous paper [S. Pieprzyk et al., Phys. Rev. E 101, 012117 (2020)]2470-004510.1103/PhysRevE.101.012117. The so-called rational-function approximation method and an approach combining accurate molecular dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation are considered. The density, composition, and size-ratio dependencies of the leading poles of the Fourier transforms of the total correlation functions h_{ij}(r) of such mixtures are presented, those poles accounting for the asymptotic decay of h_{ij}(r) for large r. Structural crossovers, in which the asymptotic wavelength of the oscillations of the total correlation functions changes discontinuously, are investigated. The behavior of the structural crossover lines as the size ratio and densities of the two species are changed is also discussed.

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 the Jagla fluid.

The structural properties of the Jagla fluid are studied by Monte Carlo (MC) simulations, numerical solutions of integral equation theories, and the (semi-analytical) rational-function approximation (RFA) method. In the latter case, the results are obtained from the assumption (supported by our MC simulations) that the Jagla potential and a potential with a hard core plus an appropriate piecewise constant function lead to practically the same cavity function. The predictions obtained for the radial distribution function, g(r), from this approach are compared against MC simulations and integral equations for the Jagla model, and also for the limiting cases of the triangle-well potential and the ramp potential, with a general good agreement. The analytical form of the RFA in Laplace space allows us to describe the asymptotic behavior of g(r) in a clean way and compare it with MC simulations for representative states with oscillatory or monotonic decay. The RFA predictions for the Fisher-Widom and Widom lines of the Jagla fluid are obtained.

Virial coefficients, equation of state, and demixing of binary asymmetric nonadditive hard-disk mixtures.

Values of the fifth virial coefficient, compressibility factors, and fluid-fluid coexistence curves of binary asymmetric nonadditive mixtures of hard disks are reported. The former correspond to a wide range of size ratios and positive nonadditivities and have been obtained through a standard Monte Carlo method for the computation of the corresponding cluster integrals. The compressibility factors as functions of density, derived from canonical Monte Carlo simulations, have been obtained for two values of the size ratio (q = 0.4 and q = 0.5), a value of the nonadditivity parameter (Δ = 0.3), and five values of the mole fraction of the species with the biggest diameter (x1 = 0.1, 0.3, 0.5, 0.7, and 0.9). Some points of the coexistence line relative to the fluid-fluid phase transition for the same values of the size ratios and nonadditivity parameter have been obtained from Gibbs ensemble Monte Carlo simulations. A comparison is made between the numerical results and those that follow from some theoretical equations of state.

Equation of state of polydisperse hard-disk mixtures in the high-density regime.

A proposal to link the equation of state of a monocomponent hard-disk fluid to the equation of state of a polydisperse hard-disk mixture is presented. Event-driven molecular dynamics simulations are performed to obtain data for the compressibility factor of the monocomponent fluid and of 26 polydisperse mixtures with different size distributions. Those data are used to assess the proposal and to infer the values of the compressibility factor of the monocomponent hard-disk fluid in the metastable region from those of mixtures in the high-density region. The collapse of the curves for the different mixtures is excellent in the stable region. In the metastable regime, except for two mixtures in which crystallization is present, the outcome of the approach exhibits a rather good performance. The simulation results indicate that a (reduced) variance of the size distribution larger than about 0.01 is sufficient to avoid crystallization and explore the metastable fluid branch.

Virial coefficients and demixing in the Asakura-Oosawa model.

The problem of demixing in the Asakura-Oosawa colloid-polymer model is considered. The critical constants are computed using truncated virial expansions up to fifth order. While the exact analytical results for the second and third virial coefficients are known for any size ratio, analytical results for the fourth virial coefficient are provided here, and fifth virial coefficients are obtained numerically for particular size ratios using standard Monte Carlo techniques. We have computed the critical constants by successively considering the truncated virial series up to the second, third, fourth, and fifth virial coefficients. The results for the critical colloid and (reservoir) polymer packing fractions are compared with those that follow from available Monte Carlo simulations in the grand canonical ensemble. Limitations and perspectives of this approach are pointed out.

Simple effective rule to estimate the jamming packing fraction of polydisperse hard spheres.

A recent proposal in which the equation of state of a polydisperse hard-sphere mixture is mapped onto that of the one-component fluid is extrapolated beyond the freezing point to estimate the jamming packing fraction ϕJ of the polydisperse system as a simple function of M1M3/M22, where Mk is the kth moment of the size distribution. An analysis of experimental and simulation data of ϕJ for a large number of different mixtures shows a remarkable general agreement with the theoretical estimate. To give extra support to the procedure, simulation data for seventeen mixtures in the high-density region are used to infer the equation of state of the pure hard-sphere system in the metastable region. An excellent collapse of the inferred curves up to the glass transition and a significant narrowing of the different out-of-equilibrium glass branches all the way to jamming are observed. Thus, the present approach provides an extremely simple criterion to unify in a common framework and to give coherence to data coming from very different polydisperse hard-sphere mixtures.

Structural properties of fluids interacting via piece-wise constant potentials with a hard core.

The structural properties of fluids whose molecules interact via potentials with a hard core plus two piece-wise constant sections of different widths and heights are presented. These follow from the more general development previously introduced for potentials with a hard core plus n piece-wise constant sections [A. Santos, S. B. Yuste, and M. Lopez de Haro, Condens. Matter Phys. 15, 23602 (2012)] in which use was made of a semi-analytic rational-function approximation method. The results of illustrative cases comprising eight different combinations of wells and shoulders are compared both with simulation data and with those that follow from the numerical solution of the Percus-Yevick and hypernetted-chain integral equations. It is found that the rational-function approximation generally predicts a more accurate radial distribution function than the Percus-Yevick theory and is comparable or even superior to the hypernetted-chain theory. This superiority over both integral equation theories is lost, however, at high densities, especially as the widths of the wells and/or the barriers increase.

Fourth virial coefficients of asymmetric nonadditive hard-disk mixtures.

The fourth virial coefficient of asymmetric nonadditive binary mixtures of hard disks is computed with a standard Monte Carlo method. Wide ranges of size ratio (0.05 ≤ q ≤ 0.95) and nonadditivity (-0.5 ≤ Δ ≤ 0.5) are covered. A comparison is made between the numerical results and those that follow from some theoretical developments. The possible use of these data in the derivation of new equations of state for these mixtures is illustrated by considering a rescaled virial expansion truncated to fourth order. The numerical results obtained using this equation of state are compared with Monte Carlo simulation data in the case of a size ratio q = 0.7 and two nonadditivities Δ = ±0.2.

Communication: inferring the equation of state of a metastable hard-sphere fluid from the equation of state of a hard-sphere mixture at high densities.

A possible approximate route to obtain the equation of state of the monodisperse hard-sphere system in the metastable fluid region from the knowledge of the equation of state of a hard-sphere mixture at high densities is discussed. The proposal is illustrated by using recent Monte Carlo simulation data for the pressure of a binary mixture. It is further shown to exhibit high internal consistency.

Virial coefficients, thermodynamic properties, and fluid-fluid transition of nonadditive hard-sphere mixtures.

Different theoretical approaches for the thermodynamic properties and the equation of state for multicomponent mixtures of nonadditive hard spheres in d dimensions are presented in a unified way. These include the theory by Hamad, our previous formulation, the original MIX1 theory, a recently proposed modified MIX1 theory, as well as a nonlinear extension of the MIX1 theory proposed in this paper. Explicit expressions for the compressibility factor, Helmholtz free energy, and second, third, and fourth virial coefficients are provided. A comparison is carried out with recent Monte Carlo data for the virial coefficients of asymmetric mixtures and with available simulation data for the compressibility factor, the critical consolute point, and the liquid-liquid coexistence curves. The merits and limitations of each theory are pointed out.

Survival probability of a subdiffusive particle in a d-dimensional sea of mobile traps.

We investigate the long-time behavior of the survival probability P(t) of a mobile particle in d-dimensional continuous Euclidean media doped with noninteracting mobile traps. The particle is strictly subdiffusive, implying that its mean-square displacement grows as tgamma' with 0or=2 we find that at late times the survival probability is that of the pure target problem (the problem where the particle remains immobile) in agreement with previous studies for the d=1 case. These decay laws remain invariant over the whole gamma range as opposed to the dynamical crossover observed for the case of a purely diffusive particle (gamma'=1) where, for gamma<2/(2+d) , the survival probability becomes that of the so-called trapping problem (the problem where the particle moves in a sea of static traps). This behavior implies that for sufficiently low values of gamma(gamma<2/(2+d)) the survival probability becomes singular in the limit gamma'-->1: trappinglike for gamma'=1 and targetlike for any gamma'<1.

Simple equation of state for hard disks on the hyperbolic plane.

A simple equation of state for hard disks on the hyperbolic plane is proposed. It yields the exact second virial coefficient and contains a pole at the highest possible packing. A comparison with another very recent theoretical proposal and simulation data is presented.

Low-temperature and high-temperature approximations for penetrable-sphere fluids: comparison with Monte Carlo simulations and integral equation theories.

The two-body interaction in dilute solutions of polymer chains in good solvents can be modeled by means of effective bounded potentials, the simplest of which being that of penetrable spheres (PSs). In this paper we construct two simple analytical theories for the structural properties of PS fluids: a low-temperature (LT) approximation, that can be seen as an extension to PSs of the well-known solution of the Percus-Yevick (PY) equation for hard spheres, and a high-temperature (HT) approximation based on the exact asymptotic behavior in the limit of infinite temperature. Monte Carlo simulations for a wide range of temperatures and densities are performed to assess the validity of both theories. It is found that, despite their simplicity, the HT and LT approximations exhibit a fair agreement with the simulation data within their respective domains of applicability, so that they complement each other. A comparison with numerical solutions of the PY and the hypernetted-chain approximations is also carried out, the latter showing a very good performance, except inside the core at low temperatures.

Multicomponent fluid of hard spheres near a wall.

The rational function approximation method, density functional theory, and NVT Monte Carlo simulation are used to obtain the density profiles of multicomponent hard-sphere mixtures near a planar hard wall. Binary mixtures with a size ratio 1:3 in which both components occupy a similar volume are specifically examined. The results indicate that the present version of density functional theory yields an excellent overall performance. A reasonably accurate behavior of the rational function approximation method is also observed, except in the vicinity of the first minimum, where it may even predict unphysical negative values.

How "sticky" are short-range square-well fluids?

The aim of this work is to investigate to what extent the structural properties of a short-range square-well (SW) fluid of range lambda at a given packing fraction eta and reduced temperature T* = kBT/epsilon can be represented by those of a sticky-hard-sphere (SHS) fluid at the same packing fraction and an effective stickiness parameter tau(T*,lambda). Such an equivalence cannot hold for the radial distribution function g(r) since this function has a delta singularity at contact (r = sigma) in the SHS case, while it has a jump discontinuity at r = lambda sigma in the SW case. Therefore, the equivalence is explored with the cavity function y(r), i.e., we assume that [formula: see text]. Optimization of the agreement between y(SW) and y(SHS) to first order in density suggests the choice tau(T*,lambda) = [12(e(1/T* - 1)(lambda - 1)](-1). We have performed Monte Carlo (MC) simulations of the SW fluid for lambda = 1.05, 1.02, and 1.01 at several densities and temperatures T* such that tau(T*,lambda) = 0.13, 0.2, and 0.5. The resulting cavity functions have been compared with MC data of SHS fluids obtained by Miller and Frenkel[J. Phys.: Condens. Matter 16, S4901 (2004)]. Although, at given values of eta and tau, some local discrepancies between y(SW) and y(SHS) exist (especially for lambda = 1.05), the SW data converge smoothly toward the SHS values as lambda-1 decreases. In fact, precursors of the singularities of y(SHS) at certain distances due to geometrical arrangements are clearly observed in y(SW). The approximate mapping y(SW)-->y(SHS) is exploited to estimate the internal energy and structure factor of the SW fluid from those of the SHS fluid. Taking for y(SHS) the solution of the Percus-Yevick equation as well as the rational-function approximation, the radial distribution function g(r) of the SW fluid is theoretically estimated and a good agreement with our MC simulations is found. Finally, a similar study is carried out for short-range SW fluid mixtures.

Contact values of the particle-particle and wall-particle correlation functions in a hard-sphere polydisperse fluid.

The contact values g(sigma,sigma') of the radial distribution functions of a fluid of (additive) hard spheres with a given size distribution f(sigma) are considered. A "universality" assumption is introduced, according to which, at a given packing fraction eta,g(sigma,sigma')=G(z(sigma,sigma')), where G is a common function independent of the number of components (either finite or infinite) and z(sigma,sigma')=[2sigmasigma'/(sigma+sigma')]mu2/mu3 is a dimensionless parameter, mu n being the nth moment of the diameter distribution. A cubic form proposal for the z dependence of G is made and known exact consistency conditions for the point particle and equal size limits, as well as between two different routes to compute the pressure of the system in the presence of a hard wall, are used to express Gz in terms of the radial distribution at contact of the one-component system. For polydisperse systems we compare the contact values of the wall-particle correlation function and the compressibility factor with those obtained from recent Monte Carlo simulations.

Structure of ternary additive hard-sphere fluid mixtures.

Monte Carlo simulations on the structural properties of ternary fluid mixtures of additive hard spheres are reported. The results are compared with those obtained from a recent analytical approximation [S. B. Yuste, A. Santos, and M. López de Haro, J. Chem. Phys. 108, 3683 (1998)] to the radial distribution functions of hard-sphere mixtures and with the results derived from the solution of the Ornstein-Zernike integral equation with both the Martynov-Sarkisov and the Percus-Yevick closures. Very good agreement between the results of the first two approaches and simulation is observed, with a noticeable improvement over the Percus-Yevick predictions especially near contact.