Search

The phase behavior and structural properties of hard anisotropic particles (prisms and dumbbells) are examined in one-dimensional channels using the Parsons–Lee (PL) theory, and the transfer-matrix and neighbor-distribution methods. The particles are allowed to move freely along the channel, while their orientations are constrained such that one particle can occupy only two or three different lengths along the channel. In this confinement setting, hard prisms behave as an additive mixture, while hard dumbbells behave as a non-additive one. We prove that all methods provide exact results for the phase properties of hard prisms, while only the neighbor-distribution and transfer-matrix methods are exact for hard dumbbells. This shows that non-additive effects are incorrectly included into the PL theory, which is a successful theory of the isotropic-nematic phase transition of rod-like particles in higher dimensions. In the one-dimensional channel, the orientational ordering develops continuously with increasing density, i.e., the system is isotropic only at zero density, while it becomes perfectly ordered at the close-packing density. We show that there is no orientational correlation in the hard prism system, while the hard dumbbells are orientationally correlated with diverging correlation length at close packing. On the other hand, positional correlations are present for all the systems, the associated correlation length diverging at close packing. [ABSTRACT FROM AUTHOR]

The structural properties of confined single-file hard-disk fluids are studied analytically by means of a mapping of the original system onto a one-dimensional mixture of non-additive hard rods, the mapping being exact in the polydisperse limit. Standard statistical-mechanical results are used as a starting point to derive thermodynamic and structural properties of the one-dimensional mixture, where the condition that all particles have the same chemical potential must be taken into account. Analytical results are then provided for the nth neighbor probability distribution function, the radial distribution function, and the structure factor, a very good agreement being observed upon comparison with simulation data from the literature. Moreover, we have analyzed the scaling form for the disappearance of defects in the zigzag configuration for high pressure and have obtained the translational correlation length and the structural crossover in the oscillation frequency for asymptotically large distances. [ABSTRACT FROM AUTHOR]

The exact transfer-matrix solution for the longitudinal equilibrium properties of the single-file hard-disk fluid is used to study the limiting low- and high-pressure behaviors analytically as functions of the pore width. In the low-pressure regime, the exact third and fourth virial coefficients are obtained, which involve single and double integrals, respectively. Moreover, we show that the standard irreducible diagrams do not provide a complete account of the virial coefficients in confined geometries. The asymptotic equation of state in the high-pressure limit is seen to present a simple pole at the close-packing linear density, as in the hard-rod fluid, but, in contrast to the latter case, the residue is 2. Since, for an arbitrary pressure, the exact transfer-matrix treatment requires the numerical solution of an eigenvalue integral equation, we propose here two simple approximations to the equation of state, with different complexity levels, and carry out an extensive assessment of their validity and practical convenience vs the exact solution and available computer simulations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

A detailed comparison between the Boublík–Mansoori–Carnahan–Starling–Leland (BMCSL) equation of state of hard-sphere mixtures is made with Molecular Dynamics (MD) simulations of the same compositions. The Labík and Smith simulation technique [S. Labík and W. R. Smith, Mol. Simul. 12, 23–31 (1994)] was used to implement the Widom particle insertion method to calculate the excess chemical potential, β μ 0 ex , of a test particle of variable diameter, σ0, immersed in a hard-sphere fluid mixture with different compositions and values of the packing fraction, η. Use is made of the fact that the only polynomial representation of β μ 0 ex which is consistent with the limits σ0 → 0 and σ0 → ∞ has to be of the cubic form, i.e., c 0 ( η ) + c ¯ 1 ( η ) σ 0 / M 1 + c ¯ 2 ( η ) ( σ 0 / M 1 ) 2 + c ¯ 3 ( η ) ( σ 0 / M 1 ) 3 , where M1 is the first moment of the distribution. The first two coefficients, c0(η) and c ¯ 1 ( η ) , are known analytically, while c ¯ 2 ( η ) and c ¯ 3 ( η ) were obtained by fitting the MD data to this expression. This in turn provides a method to determine the excess free energy per particle, βaex, in terms of c ¯ 2 , c ¯ 3 , and the compressibility factor, Z. Very good agreement between the BMCSL formulas and the MD data is found for β μ 0 ex , Z, and βaex for binary mixtures and continuous particle size distributions with the top-hat analytic form. However, the BMCSL theory typically slightly underestimates the simulation values, especially for Z, differences which the Boublík–Carnahan–Starling–Kolafa formulas and an interpolation between two Percus–Yevick routes capture well in different ranges of the system parameter space. [ABSTRACT FROM AUTHOR]

The interest in the description of the properties of fluids of restricted dimensionality is growing for theoretical and practical reasons. In this work, we have firstly developed an analytical expression for the Helmholtz free energy of the two-dimensional square-well fluid in the Barker–Henderson framework. This equation of state is based on an approximate analytical radial distribution function for d-dimensional hard-sphere fluids (1 ≤ d ≤ 3) and is validated against existing and new simulation results. The so-obtained equation of state is implemented in a discrete perturbation theory able to account for general potential shapes. The prototypical Lennard-Jones and Yukawa fluids are tested in its two-dimensional version against available and new simulation data with semiquantitative agreement. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The Labík and Smith Monte Carlo simulation technique to implement the Widom particle insertion method is applied using Molecular Dynamics (MD) instead to calculate numerically the insertion probability, P0(η, σ0), of tracer hard-sphere (HS) particles of different diameters, σ0, in a host HS fluid of diameter σ and packing fraction, η, up to 0.5. It is shown analytically that the only polynomial representation of - ln P0(η, σ0) consistent with the limits σ0 → 0 and σ0 → ∞ has necessarily a cubic form, c0(η)+c1(η)σ0=σ +c2(η)(σ0=σ)² +c3(η)(σ0=σ)³. Our MD data for - ln P0(η, σ0) are fitted to such a cubic polynomial and the functions c0(η) and c1(η) are found to be statistically indistinguishable from their exact solution forms. Similarly, c2(η) and c3(η) agree very well with the Boublík-Mansoori- Carnahan-Starling-Leland and Boublík-Carnahan-Starling-Kolafa formulas. The cubic polynomial is extrapolated (high density) or interpolated (low density) to obtain the chemical potential of the host fluid, or σ0 → σ, as βαex = c0 + c1 + c2 + c3. Excellent agreement between the Carnahan-Starling and Carnahan-Starling-Kolafa theories with our MD data is evident. [ABSTRACT FROM AUTHOR]

The relevance of neglecting three- and four-body interactions in the coarse-grained version of the Asakura-Oosawa model is examined. A mapping between the first few virial coefficients of the binary nonadditive hard-sphere mixture representative of this model and those arising from the coarsegrained (pairwise) depletion potential approximation allows for a quantitative evaluation of the effect of such interactions. This turns out to be especially important for large size ratios and large reservoir polymer packing fractions. [ABSTRACT FROM AUTHOR]

A general class of nonadditive sticky-hard-sphere binary mixtures, where small and large spheres represent the solvent and the solute, respectively, is introduced. The solute-solute and solvent-solvent interactions are of hard-sphere type, while the solute-solvent interactions are of sticky-hard-sphere type with tunable degrees of size nonadditivity and stickiness. Two particular and complementary limits are studied using analytical and semi-analytical tools. The first case is characterized by zero nonadditivity, lending itself to a Percus-Yevick approximate solution from which the impact of stickiness on the spinodal curves and on the effective solute-solute potential is analyzed. In the opposite nonadditive case, the solvent-solvent diameter is zero and the model can then be reckoned as an extension of the well-known Asakura-Oosawa model with additional sticky solute-solvent interaction. This latter model has the property that its exact effective one-component problem involves only solute-solute pair potentials for size ratios such that a solvent particle fits inside the interstitial region of three touching solutes. In particular, we explicitly identify the three competing physical mechanisms (depletion, pulling, and bridging) giving rise to the effective interaction. Some remarks on the phase diagram of these two complementary models are also addressed through the use of the Noro-Frenkel criterion and a first-order perturbation analysis. Our findings suggest reentrance of the fluid-fluid instability as solvent density (in the first model) or adhesion (in the second model) is varied. Some perspectives in terms of the interpretation of recent experimental studies of microgels adsorbed onto large polystyrene particles are discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The mutual entropic depletion force felt by two solute "big" hard spheres immersed in a binary mixture solvent of nonadditive "small" hard spheres is calculated as a function of the surface-to-surface distance by means of canonical Monte Carlo simulations and through a recently proposed rational-function approximation [R. Fantoni and A. Santos, Phys. Rev. E 84, 041201 (2011)]. Four representative scenarios are investigated: symmetric solute particles and the limit where one of the two solute spheres becomes a planar hard wall, in both cases with symmetric and asymmetric solvents. In all cases, the influence on the depletion force due to the nonadditivity in the solvent is determined in the mixed state. Comparison between results from the theoretical approximation and from the simulation shows a good agreement for surface-to-surface distances greater than the smallest solvent diameter. [ABSTRACT FROM AUTHOR]

The fourth virial coefficient of additive hard-sphere mixtures, as predicted by the Percus-Yevick (PY) and hypernetted-chain (HNC) theories, is derived via the compressibility, virial, and chemicalpotential routes, the outcomes being compared with exact results. Except in the case of the HNC compressibility route, the other five expressions exhibit a common structure involving the first three moments of the size distribution. In both theories, the chemical-potential route is slightly better than the virial one and the best behavior is generally presented by the compressibility route. Moreover, the PY results with any of the three routes are more accurate than any of the HNC results. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The question of whether the known virial coefficients are enough to determine the packing fraction η∞ at which the fluid equation of state of a hard-sphere fluid diverges is addressed. It is found that the information derived from the direct Padé approximants to the compressibility factor constructed with the virial coefficients is inconclusive. An alternative approach is proposed which makes use of the same virial coefficients and of the equation of state in a form where the packing fraction is explicitly given as a function of the pressure. The results of this approach both for hard-disk and hard-sphere fluids, which can straightforwardly accommodate higher virial coefficients when available, lends support to the conjecture that η∞ is equal to the maximum packing fraction corresponding to an ordered crystalline structure. [ABSTRACT FROM AUTHOR]

The one-dimensional penetrable-square-well fluid is studied using both analytical tools and specialized Monte Carlo simulations. The model consists of a penetrable core characterized by a finite repulsive energy combined with a short-range attractive well. This is a many-body one-dimensional problem, lacking an exact analytical solution, for which the usual van Hove theorem on the absence of phase transition does not apply. We determine a high-penetrability approximation complementing a similar low-penetrability approximation presented in previous work. This is shown to be equivalent to the usual Debye–Hückel theory for simple charged fluids for which the virial and energy routes are identical. The internal thermodynamic consistency with the compressibility route and the validity of the approximation in describing the radial distribution function is assessed by a comparison against numerical simulations. The Fisher–Widom line separating the oscillatory and monotonic large-distance behaviors of the radial distribution function is computed within the high-penetrability approximation and compared with the opposite regime, thus providing a strong indication of the location of the line in all possible regimes. The high-penetrability approximation predicts the existence of a critical point and a spinodal line, but this occurs outside the applicability domain of the theory. We investigate the possibility of a fluid-fluid transition by the Gibbs ensemble Monte Carlo techniques, not finding any evidence of such a transition. Additional analytical arguments are given to support this claim. Finally, we find a clustering transition when Ruelle’s stability criterion is not fulfilled. The consequences of these findings on the three-dimensional phase diagrams are also discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

As is well known, approximate integral equations for liquids, such as the hypernetted chain (HNC) and Percus–Yevick (PY) theories, are in general thermodynamically inconsistent in the sense that the macroscopic properties obtained from the spatial correlation functions depend on the route followed. In particular, the values of the fourth virial coefficient B4 predicted by the HNC and PY approximations via the virial route differ from those obtained via the compressibility route. Despite this, it is shown in this paper that the value of B4 obtained from the virial route in the HNC theory is exactly three halves the value obtained from the compressibility route in the PY theory, irrespective of the interaction potential (whether isotropic or not), the number of components, and the dimensionality of the system. This simple relationship is confirmed in one-component systems by analytical results for the one-dimensional penetrable-square-well model and the three-dimensional penetrable-sphere model, as well as by numerical results for the one-dimensional Lennard-Jones model, the one-dimensional Gaussian core model, and the three-dimensional square-well model. [ABSTRACT FROM AUTHOR]

The Kirkwood–Buff integrals for two-component mixtures in one-dimensional systems are calculated directly. The results are applied to square-well particles and found to agree with those obtained by the inversion of the Kirkwood–Buff theory of solutions. [ABSTRACT FROM AUTHOR]

We study structural and thermophysical properties of a one-dimensional classical fluid made of penetrable spheres interacting via an attractive square-well potential. Penetrability of the spheres is enforced by reducing from infinite to finite the repulsive energy barrier in the pair potentials As a consequence, an exact analytical solution is lacking even in one dimension. Building upon previous exact analytical work in the low-density limit [A. Santos, R. Fantoni, and A. Giacometti, Phys. Rev. E 77, 051206 (2008)], we propose an approximate theory valid at any density and in the low-penetrable regime. By comparison with specialized Monte Carlo simulations and integral equation theories, we assess the regime of validity of the theory. We investigate the degree of inconsistency among the various routes to thermodynamics and explore the possibility of a fluid-fluid transition. Finally we locate the dependence of the Fisher–Widom line on the degree of penetrability. Our results constitute the first systematic study of penetrable spheres with attractions as a prototype model for soft systems. [ABSTRACT FROM AUTHOR]

Using the first seven known virial coefficients and forcing it to possess two branch-point singularities, a new equation of state for the hard-sphere fluid is proposed. This equation of state predicts accurate values of the higher virial coefficients, a radius of convergence smaller than the close-packing value, and it is as accurate as the rescaled virial expansion and better than the Padé [3/3] equations of state. Consequences regarding the convergence properties of the virial series and the use of similar equations of state for hard-core fluids in d dimensions are also pointed out. [ABSTRACT FROM AUTHOR]

The depletion force and depletion potential between two in principle unequal “big” hard spheres embedded in a multicomponent mixture of “small” hard spheres are computed using the rational function approximation method for the structural properties of hard-sphere mixtures [S. B. Yuste, A. Santos, and M. López de Haro, J. Chem. Phys. 108, 3683 (1998)]. The cases of equal solute particles and of one big particle and a hard planar wall in a background monodisperse hard-sphere fluid are explicitly analyzed. An improvement over the performance of the Percus–Yevick theory and good agreement with available simulation results are found. [ABSTRACT FROM AUTHOR]

The aim of this work is to investigate to what extent the structural properties of a short-range square-well (SW) fluid of range λ at a given packing fraction η and reduced temperature T*=kBT/ε can be represented by those of a sticky-hard-sphere (SHS) fluid at the same packing fraction and an effective stickiness parameter τ(T*,λ). Such an equivalence cannot hold for the radial distribution function g(r) since this function has a delta singularity at contact (r=σ) in the SHS case, while it has a jump discontinuity at r=λσ in the SW case. Therefore, the equivalence is explored with the cavity function y(r), i.e., we assume that ySW(r|η,T*;λ)≈ySHS(r|η,τ(T*,λ)). Optimization of the agreement between ySW and ySHS to first order in density suggests the choice τ(T*,λ)=[12(e1/T*-1)(λ-1)]-1. We have performed Monte Carlo (MC) simulations of the SW fluid for λ=1.05, 1.02, and 1.01 at several densities and temperatures T* such that τ(T*,λ)=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 η and τ, some local discrepancies between ySW and ySHS exist (especially for λ=1.05), the SW data converge smoothly toward the SHS values as λ-1 decreases. In fact, precursors of the singularities of ySHS at certain distances due to geometrical arrangements are clearly observed in ySW. The approximate mapping ySW→ySHS is exploited to estimate the internal energy and structure factor of the SW fluid from those of the SHS fluid. Taking for ySHS 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. [ABSTRACT FROM AUTHOR]

Bounded potentials are good models to represent the effective two-body interaction in some colloidal systems, such as the dilute solutions of polymer chains in good solvents. The simplest bounded potential is that of penetrable spheres, which takes a positive finite value if the two spheres are overlapped, being 0 otherwise. Even in the one-dimensional case, the penetrable-rod model is far from trivial, since interactions are not restricted to nearest neighbors and so its exact solution is not known. In this paper the structural properties of one-dimensional penetrable rods are studied. We first derive the exact correlation functions of the penetrable-rod fluids to second order in density at any temperature, as well as in the high-temperature and zero-temperature limits at any density. It is seen that, in contrast to what is generally believed, the Percus-Yevick equation does not yield the exact cavity function in the hard-rod limit. Next, two simple analytic theories are constructed: a high-temperature approximation based on the exact asymptotic behavior in the limit T→∞ and a low-temperature approximation inspired by the exact result in the opposite limit T→0. Finally, we perform Monte Carlo simulations for a wide range of temperatures and densities to assess the validity of both theories. It is found that they complement each other quite well, exhibiting a good agreement with the simulation data within their respective domains of applicability and becoming practically equivalent on the borderline of those domains. A comparison with numerical solutions of the Percus-Yevick and the hypernetted-chain approximations is also carried out. Finally, a perspective on the extension of our two heuristic theories to the more realistic three-dimensional case is provided. [ABSTRACT FROM AUTHOR]

The contact values g(σ,σ′) of the radial distribution functions of a fluid of (additive) hard spheres with a given size distribution f(σ) are considered. A “universality” assumption is introduced, according to which, at a given packing fraction η,g(σ,σ′)=G(z(σ,σ′)), where G is a common function independent of the number of components (either finite or infinite) and z(σ,σ′)=[2σσ′/(σ+σ′)]μ2/μ3 is a dimensionless parameter, μ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 G(z) 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. [ABSTRACT FROM AUTHOR]

We report molecular dynamics results for the contact values of the radial distribution functions of binary additive mixtures of hard disks. The simulation data are compared with theoretical predictions from expressions proposed by Jenkins and Mancini [J. Appl. Mech. 54, 27 (1987)] and Santos et al. [J. Chem. Phys. 117, 5785 (2002)]. Both theories agree quantitatively within a very small margin, which renders the former still a very useful and simple tool to work with. The latter (higher-order and self-consistent) theory provides a small qualitative correction for low densities and is superior especially in the high-density domain.© 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular-dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented. [ABSTRACT FROM AUTHOR]

The linearized Debye–Hückel theory for liquid state is shown to provide thermodynamically consistent virial and energy routes for any potential and for any dimensionality. The importance of this result for bounded potentials is discussed. [ABSTRACT FROM AUTHOR]

The direct correlation function and the (static) structure factor for a seven-dimensional hard-sphere fluid are considered. Analytical results for these quantities are derived within the Percus-Yevick [Phys. Rev.110, 1 (1958)] theory. [ABSTRACT FROM AUTHOR]

The article offers information on the equation of hard-sphere (HS) model. Among the key topics that were discussed include computer simulations, presentation of equation of state (EOS) as graph in compressibility factor, and scaling parameters of HS absolute temperatures. Other topics mentioned are the Monte Carlo (MC) method, coexistence pressure, and freezing point determination.

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. [ABSTRACT FROM AUTHOR]

It is proven that, for any soft potential φ(r) characterized by a finite Fourier transform φ(k), the virial and energy thermodynamic routes are equivalent if the Fourier transform of the total correlation function divided by the density ρ, h(k)/ρ, is an arbitrary function of ρβφ(k), where β is the inverse temperature. This class of approximations includes the mean spherical approximation as a particular case. [ABSTRACT FROM AUTHOR]

The energy route to the equation of state of hard-sphere fluids is ill defined since the internal energy is just that of an ideal gas, and thus it is independent of density. It is shown that this ambiguity can be avoided by considering a square-shoulder interaction and taking the limit of vanishing shoulder width. The resulting hard-sphere equation of state coincides exactly with the one obtained through the virial route. Therefore, the energy and virial routes to the equation of state of hard-sphere fluids can be considered as equivalent. [ABSTRACT FROM AUTHOR]

A class of binary mixtures of Janus fluids formed by colloidal spheres with the hydrophobic hemispheres constrained to point either up or down are studied by means of Gibbs ensemble Monte Carlo simulations and simple analytical approximations. These fluids can be experimentally realized by the application of an external static electrical field. The gas-liquid and demixing phase transitions in five specific models with different patch-patch affinities are analyzed. It is found that a gas-liquid transition is present in all the models, even if only one of the four possible patch-patch interactions is attractive. Moreover, provided the attraction between like particles is stronger than between unlike particles, the system demixes into two subsystems with different composition at sufficiently low temperatures and high densities. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The problem of demixing in a binary fluid mixture of highly asymmetric additive hard spheres is revisited. A comparison is presented between the results derived previously using truncated virial expansions for three finite size ratios with those that one obtains with the same approach in the extreme case in which one of the components consists of point particles. Since this latter system is known not to exhibit fluid-fluid segregation, the similarity observed for the behavior of the critical constants arising in the truncated series in all instances, while not being conclusive, may cast serious doubts as to the actual existence of a demixing fluid-fluid transition in disparate-sized binary additive hard-sphere mixtures. [ABSTRACT FROM AUTHOR]

We study thermophysical properties of a Janus fluid with constrained orientations, using analytical techniques and numerical simulations. The Janus character is modeled by means of a Kern-Frenkel potential where each sphere has one hemisphere of square-well and the other of hard-sphere character. The orientational constraint is enforced by assuming that each hemisphere can only point either North or South with equal probability. The analytical approach hinges on a mapping of the above Janus fluid onto a binary mixture interacting via a 'quasi' isotropic potential. The anisotropic nature of the original Kern-Frenkel potential is reflected by the asymmetry in the interactions occurring between the unlike components of the mixture. A rational-function approximation extending the corresponding symmetric case is obtained in the sticky limit, where the square-well becomes infinitely narrow and deep, and allows a fully analytical approach. Notwithstanding the rather drastic approximations in the analytical theory, this is shown to provide a rather precise estimate of the structural and thermodynamical properties of the original Janus fluid. [ABSTRACT FROM AUTHOR]

A theoretical model for polydisperse systems of hard spheres is said to be truncatable when the excess free energy depends on the size distribution through a finite number K of moments. This Note proves an exact scaling relation for truncatable free energies, which allows to reduce the effective degrees of freedom to K - 2 independent combinations of the moments. [ABSTRACT FROM AUTHOR]

A correction to the article "Equation of state of a seven-dimensional hard-sphere fluid. Percus-Yevick theory and molecular-dynamics simulations," that was published in the 2004 issue is presented.

The equation of state for five-dimensional hard hyperspheres arising as a weighted average of the Percus-Yevick compressibility ([FRACTION SHAPE="CASE"][NUM]3[/NUM][DEN]5[/DEN][/FRACTION]) and virial ([FRACTION SHAPE="CASE"][NUM]2[/NUM][DEN]5[/DEN][/FRACTION]) equations of state is considered. This Carnahan-Starling-type equation turns out to be extremely accurate, despite the fact that both Percus-Yevick equations of state are rather poor. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

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.

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.

A recently derived method [R. D. Rohrmann and A. Santos, Phys. Rev. E 76, 051202 (2007)] to obtain the exact solution of the Percus-Yevick equation for a fluid of hard spheres in (odd) d dimensions is used to investigate the convergence properties of the resulting virial series. This is done both for the virial and compressibility routes, in which the virial coefficients B(j) are expressed in terms of the solution of a set of (d-1)/2 coupled algebraic equations which become nonlinear for d>/=5. Results have been derived up to d=13. A confirmation of the alternating character of the series for d>/=5, due to the existence of a branch point on the negative real axis, is found and the radius of convergence is explicitly determined for each dimension. The resulting scaled density per dimension 2eta(1/d), where eta is the packing fraction, is wholly consistent with the limiting value of 1 for d-->infinity. Finally, the values for B(j) predicted by the virial and compressibility routes in the Percus-Yevick approximation are compared with the known exact values [N. Clisby and B. M. McCoy, J. Stat. Phys. 122, 15 (2006)].