Source author record

Andrés Santos

Andrés Santos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

54works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2023arXiv

A heuristic approach for the densest packing fraction of hard-sphere mixtures

In a previous work, a simple approach to derive the jamming packing fraction of a hard-sphere mixture from the knowledge of the random close-packing fraction of the monocomponent system was proposed. Now, an extension of that approach is applied to provide an approximate formula for the densest packing fraction of a given hard-sphere mixture in terms of the fcc close-packing fraction of a monocomponent hard-sphere system and of a single parameter encapsulating the dependence on the size ratios and the number of spheres in the unit cell. Comparison with recent results for such densest packing fraction of binary and ternary systems is performed and reasonable agreement is obtained.

preprint2022arXiv

Granular gas of inelastic and rough Maxwell particles

The most widely used model for granular gases is perhaps the inelastic hard-sphere model (IHSM), where the grains are assumed to be perfectly smooth spheres colliding with a constant coefficient of normal restitution. A much more tractable model is the inelastic Maxwell model (IMM), in which the velocity-dependent collision rate is replaced by an effective mean-field constant. This simplification has been taken advantage of by many researchers to find a number of exact results within the IMM. On the other hand, both the IHSM and IMM neglect the impact of roughness -- generally present in real grains -- on the dynamic properties of a granular gas. This is remedied by the inelastic rough hard-sphere model (IRHSM), where, apart from the coefficient of normal restitution, a constant coefficient of tangential restitution is introduced. In parallel to the simplification carried out when going from the IHSM to the IMM, we propose in this paper an inelastic rough Maxwell model (IRMM) as a simplification of the IRHSM. The tractability of the proposed model is illustrated by the exact evaluation of the collisional moments of first and second degree, and the most relevant ones of third and fourth degree. The results are applied to the evaluation of the rotational-to-translational temperature ratio and the velocity cumulants in the homogeneous cooling state.

preprint2022arXiv

Structural and thermodynamic properties of fluids whose molecules interact via one-, two-, and three-step potentials

The structural and thermodynamic properties of fluids whose molecules interact via potentials with a hard-core plus a square well, a square shoulder, and a second square well, are considered. Those properties are derived by using a (semi-analytical) rational-function approximation method as a particular case of the more general formulation provided earlier involving potentials with a hard-core plus $n$ piecewise constant sections. Comparison of the results with recent simulation data confirms the usefulness of the approach.

preprint2022arXiv

Thermal versus entropic Mpemba effect in molecular gases with nonlinear drag

Loosely speaking, the Mpemba effect appears when hotter systems cool sooner or, in a more abstract way, when systems further from equilibrium relax faster. In this paper, we investigate the Mpemba effect in a molecular gas with nonlinear drag, both analytically (by employing the tools of kinetic theory) and numerically (direct simulation Monte Carlo of the kinetic equation and event-driven molecular dynamics). The analysis is carried out via two alternative routes, recently considered in the literature: first, the kinetic or thermal route, in which the Mpemba effect is characterized by the crossing of the evolution curves of the kinetic temperature (average kinetic energy), and, second, the stochastic thermodynamics or entropic route, in which the Mpemba effect is characterized by the crossing of the distance to equilibrium in probability space. In general, a nonmutual correspondence between the thermal and entropic Mpemba effects is found, i.e., there may appear the thermal effect without its entropic counterpart or vice versa. Furthermore, a nontrivial overshoot with respect to equilibrium of the thermal relaxation makes it necessary to revise the usual definition of the thermal Mpemba effect, which is shown to be better described in terms of the relaxation of the local equilibrium distribution. Our theoretical framework, which involves an extended Sonine approximation in which not only the excess kurtosis but also the sixth cumulant is retained, gives an excellent account of the behavior observed in simulations.

preprint2021arXiv

Mpemba effect in inertial suspensions

The Mpemba effect (a counterintuitive thermal relaxation process where an initially hotter system may cool down to the steady state sooner than an initially colder system) is studied in terms of a model of inertial suspensions under shear. The relaxation to a common steady state of a suspension initially prepared in a quasi-equilibrium state is compared with that of a suspension initially prepared in a nonequilibrium sheared state. Two classes of Mpemba effect are identified, the normal and the anomalous one. The former is generic, in the sense that the kinetic temperature starting from a cold nonequilibrium sheared state is overtaken by the one starting from a hot quasi-equilibrium state, due to the absence of initial viscous heating in the latter, resulting in a faster initial cooling. The anomalous Mpemba effect is opposite to the normal one since, despite the initial slower cooling of the nonequilibrium sheared state, it can eventually overtake an initially colder quasi-equilibrium state. The theoretical results based on kinetic theory agree with those obtained from event-driven simulations for inelastic hard spheres. It is also confirmed the existence of the inverse Mpemba effect, which is a peculiar heating process, in these suspensions. More particularly, we find the existence of a mixed process in which both heating and cooling can be observed during relaxation.

preprint2020arXiv

Enskog kinetic theory of rheology for a moderately dense inertial suspension

The Enskog kinetic theory for moderately dense inertial suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the background fluid on suspended particles is modeled via a viscous drag force plus a Langevin-like term defined in terms of the background temperature. In a previous paper [Hayakawa et al., Phys. Rev. E 96, 042903 (2017)], Grad's moment method with the aid of a linear shear-rate expansion was employed to obtain a theory which gave good agreement with the results of event-driven Langevin simulations of hard spheres for low densities and/or small shear rates. Nevertheless, the previous approach had a limitation of applicability to the high shear-rate and high density regime. Thus, in the present paper, we extend the previous work and develop Grad's theory including higher order terms in the shear rate. This improves significantly the theoretical predictions, a quantitative agreement between theory and simulation being found in the high-density region (volume fractions smaller than or equal to $0.4$).

preprint2020arXiv

Equation of State of Four- and Five-Dimensional Hard-Hypersphere Mixtures

New proposals for the equation of state of four- and five-dimensional hard-hypersphere mixtures in terms of the equation of state of the corresponding monocomponent hard-hypersphere fluid are introduced. Such proposals (which are constructed in such a way so as to yield the exact third virial coefficient) extend, on the one hand, recent similar formulations for hard-disk and (three-dimensional) hard-sphere mixtures and, on the other hand, two of our previous proposals also linking the mixture equation of state and the one of the monocomponent fluid but unable to reproduce the exact third virial coefficient. The old and new proposals are tested by comparison with published molecular dynamics and Monte Carlo simulation results and their relative merit is evaluated

preprint2020arXiv

Mpemba effect in molecular gases under nonlinear drag

We look into the Mpemba effect---the initially hotter sample cools sooner---in a molecular gas with nonlinear viscous drag. Specifically, the gas particles interact among them via elastic collisions and also with a background fluid at equilibrium. Thus, within the framework of kinetic theory, our gas is described by an Enskog--Fokker--Planck equation. The analysis is carried out in the first Sonine approximation, in which the evolution of the temperature is coupled to that of the excess kurtosis. This coupling leads to the emergence of the Mpemba effect, which is observed in an early stage of the relaxation and when the initial temperatures of the two samples are close enough. This allows for the development of a simple theory, linearizing the temperature evolution around a reference temperature---namely the initial temperature closer to the asymptotic equilibrium value. The linear theory provides a semiquantitative description of the effect, including expressions for the crossover time and the maximum temperature difference. We also discuss the limitations of our linearized theory.

preprint2019arXiv

Large Mpemba-like effect in a gas of inelastic rough hard spheres

We report the emergence of a giant Mpemba effect in the uniformly heated gas of inelastic rough hard spheres: The initially hotter sample may cool sooner than the colder one, even when the initial temperatures differ by more than one order of magnitude. In order to understand this behavior, it suffices to consider the simplest Maxwellian approximation for the velocity distribution in a kinetic approach. The largeness of the effect stems from the fact that the rotational and translational temperatures, which obey two coupled evolution equations, are comparable. Our theoretical predictions agree very well with molecular dynamics and direct simulation Monte Carlo data.

preprint2016arXiv

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 \leq d \leq 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 fractal dimension [M. Heinen et al., Phys. Rev. Lett. \textbf{115}, 097801 (2015)], a good agreement being observed.

preprint2015arXiv

Bridging and depletion mechanisms in colloid-colloid effective interactions: A reentrant phase diagram

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.

preprint2015arXiv

Equation of state for five-dimensional hyperspheres from the chemical-potential route

We use the Percus-Yevick approach in the chemical-potential route to evaluate the equation of state of hard hyperspheres in five dimensions. The evaluation requires the derivation of an analytical expression for the contact value of the pair distribution function between particles of the bulk fluid and a solute particle with arbitrary size. The equation of state is compared with those obtained from the conventional virial and compressibility thermodynamic routes and the associated virial coefficients are computed. The pressure calculated from all routes is exact up to third density order, but it deviates with respect to simulation data as density increases, the compressibility and the chemical-potential routes exhibiting smaller deviations than the virial route. Accurate linear interpolations between the compressibility route and either the chemical-potential route or the virial one are constructed.

preprint2015arXiv

Multi-particle critical correlations

We study the role of multi-particle spatial correlations in the appearance of a liquid-vapour critical point. Our analysis is based on the exact infinite hierarchy of equations relating spatial integrals of $(k+1)$-particle correlations to the $k$-particle ones and their derivatives with respect to the density. Critical exponents corresponding to generalized compressibility equations resulting from the hierarchy are shown to grow linearly with the order of correlations. We prove that the critical behaviour requires taking into account correlation functions of arbitrary order. It is however only a necessary condition. Indeed, approximate closures of the hierarchy obtained by expressing higher order correlations in terms of lower order ones (Kirkwood's superposition approximation and its generalizations) turn out to be inconsistent with the critical behaviour.

preprint2015arXiv

Steady state in a gas of inelastic rough spheres heated by a uniform stochastic force

We study here the steady state attained in a granular gas of inelastic rough spheres that is subject to a spatially uniform random volume force. The stochastic force has the form of the so-called white noise and acts by adding impulse to the particle translational velocities. We work out an analytical solution of the corresponding velocity distribution function from a Sonine polynomial expansion that displays energy non-equipartition between the translational and rotational modes, translational and rotational kurtoses, and translational-rotational velocity correlations. By comparison with a numerical solution of the Boltzmann kinetic equation (by means of the Direct Simulation Monte Carlo method) we show that our analytical solution provides a good description that is quantitatively very accurate in certain ranges of inelasticity and roughness. We also find three important features that make the forced granular gas steady state very different from the homogeneous cooling state (attained by an unforced granular gas). First, the marginal velocity distributions are always close to a Maxwellian. Second, there is a continuous transition to the purely smooth limit (where the effects of particle rotations are ignored). And third, the angular translational-rotational velocity correlations show a preference for a quasiperpendicular mutual orientation (which is called "lifted-tennis-ball" behavior).

preprint2015arXiv

The effective colloid interaction in the Asakura-Oosawa model. Assessment of non-pairwise terms from the virial expansion

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 coarse-grained (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.

preprint2015arXiv

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.

preprint2014arXiv

A 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 $ϕ_\text{J}$ of the polydisperse system as a simple function of $M_1M_3/M_2^2$, where $M_k$ is the $k$th moment of the size distribution. An analysis of experimental and simulation data of $ϕ_\text{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.

preprint2014arXiv

Depletion force in the infinite-dilution limit in a solvent of nonadditive hard spheres

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 [Phys. Rev. E \textbf{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.

preprint2014arXiv

Depletion potential in the infinite dilution limit

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. {\bf 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

preprint2014arXiv

Equation of state of sticky-hard-sphere fluids in the chemical-potential route

The coupling-parameter method, whereby an extra particle is progressively coupled to the rest of the particles, is applied to the sticky-hard-sphere fluid to obtain its equation of state in the so-called chemical-potential route ($μ$ route). As a consistency test, the results for one-dimensional sticky particles are shown to be exact. Results corresponding to the three-dimensional case (Baxter's model) are derived within the Percus-Yevick approximation by using different prescriptions for the dependence of the interaction potential of the extra particle on the coupling parameter. The critical point and the coexistence curve of the gas-liquid phase transition are obtained in the $μ$ route and compared with predictions from other thermodynamics routes and from computer simulations. The results show that the $μ$ route yields a general better description than the virial, energy, compressibility, and zero-separation routes.

preprint2014arXiv

Fourth virial coefficient of additive hard-sphere mixtures in the Percus-Yevick and hypernetted-chain approximations

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 chemical-potential 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.

preprint2014arXiv

Hydrodynamic Burnett equations for inelastic Maxwell models of granular gases

The hydrodynamic Burnett equations and the associated transport coefficients are exactly evaluated for generalized inelastic Maxwell models. In those models, the one-particle distribution function obeys the inelastic Boltzmann equation, with a velocity-independent collision rate proportional to the $γ$ power of the temperature. The pressure tensor and the heat flux are obtained to second order in the spatial gradients of the hydrodynamic fields with explicit expressions for all the Burnett transport coefficients as functions of $γ$, the coefficient of normal restitution, and the dimensionality of the system. Some transport coefficients that are related in a simple way in the elastic limit become decoupled in the inelastic case. As a byproduct, existing results in the literature for three-dimensional elastic systems are recovered, and a generalization to any dimension of the system is given. The structure of the present results is used to estimate the Burnett coefficients for inelastic hard spheres.

preprint2014arXiv

Transport coefficients of a granular gas of inelastic rough hard spheres

The Boltzmann equation for inelastic and rough hard spheres is considered as a model of a dilute granular gas. In this model, the collisions are characterized by constant coefficients of normal and tangential restitution and hence the translational and rotational degrees of freedom are coupled. A normal solution to the Boltzmann equation is obtained by means of the Chapman-Enskog method for states near the homogeneous cooling state. The analysis is carried out to first order in the spatial gradients of the number density, the flow velocity, and the granular temperature. The constitutive equations for the momentum and heat fluxes and for the cooling rate are derived and the associated transport coefficients are expressed in terms of the solutions of linear integral equations. For practical purposes, a first Sonine approximation is used to obtain explicit expressions of the transport coefficients as nonlinear functions of both coefficients of restitution and the moment of inertia. Known results for purely smooth inelastic spheres and perfectly elastic and rough spheres are recovered in the appropriate limits.

preprint2013arXiv

Chemical-potential route for multicomponent fluids

The chemical potentials of multicomponent fluids are derived in terms of the pair correlation functions for arbitrary number of components, interaction potentials, and dimensionality. The formally exact result is particularized to hard-sphere mixtures with zero or positive nonadditivity. As a simple application, the chemical potentials of three-dimensional additive hard-sphere mixtures are derived from the Percus-Yevick theory and the associated equation of state is obtained. This Percus-Yevick chemical-route equation of state is shown to be more accurate than the virial equation of state. An interpolation between the chemical-potential and compressibility routes exhibits a better performance than the well-known Boublík-Mansoori-Carnahan-Starling-Leland equation of state.

preprint2013arXiv

Communication: Virial coefficients and demixing in highly asymmetric binary additive hard-sphere mixtures

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.

preprint2013arXiv

Janus fluid with fixed patch orientations: theory and simulations

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.

preprint2013arXiv

Multicomponent fluid of nonadditive hard spheres near a wall

A recently proposed rational-function approximation [Phys. Rev. E \textbf{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.

preprint2013arXiv

Phase diagrams of Janus fluids with up-down constrained orientations

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.

preprint2013arXiv

Playing with Marbles: Structural and Thermodynamic Properties of Hard-Sphere Systems

These lecture notes present an overview of equilibrium statistical mechanics of classical fluids, with special applications to the structural and thermodynamic properties of systems made of particles interacting via the hard-sphere potential or closely related model potentials. The exact statistical-mechanical properties of one-dimensional systems, the issue of thermodynamic (in)consistency among different routes in the context of several approximate theories, and the construction of analytical or semi-analytical approximations for the structural properties are also addressed.

preprint2013arXiv

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 [Condens. Matter Phys. {\bf 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.

preprint2012arXiv

Chemical-Potential Route: A Hidden Percus-Yevick Equation of State for Hard Spheres

The chemical potential of a hard-sphere fluid can be expressed in terms of the contact value of the radial distribution function of a solute particle with a diameter varying from zero to that of the solvent particles. Exploiting the explicit knowledge of such a contact value within the Percus--Yevick (PY) theory, and using standard thermodynamic relations, a hitherto unknown PY equation of state, $p/ρk_BT=-(9/η)\ln(1-η)-(16-31η)/2(1-η)^2$, is unveiled. This equation of state turns out to be better than the one obtained from the conventional virial route. Interpolations between the chemical-potential and compressibility routes are shown to be more accurate than the widely used Carnahan--Starling equation of state. The extension to polydisperse hard-sphere systems is also presented.

preprint2012arXiv

Class of consistent fundamental-measure free energies for hard-sphere mixtures

In fundamental-measure theories the bulk excess free-energy density of a hard-sphere fluid mixture is assumed to depend on the partial number densities ${ρ_i}$ only through the four scaled-particle-theory variables ${ξ_α}$, i.e., $Φ({ρ_i})\toΦ({ξ_α})$. By imposing consistency conditions, it is proven here that such a dependence must necessarily have the form $Φ({ξ_α})=-ξ_0\ln(1-ξ_3)+Ψ(y)ξ_1ξ_2/(1-ξ_3)$, where $y\equiv {ξ_2^2}/{12πξ_1 (1-ξ_3)}$ is a scaled variable and $Ψ(y)$ is an arbitrary dimensionless scaling function which can be determined from the free-energy density of the one-component system. Extension to the inhomogeneous case is achieved by standard replacements of the variables ${ξ_α}$ by the fundamental-measure (scalar, vector, and tensor) weighted densities ${n_α(\mathbf{r})}$. Comparison with computer simulations shows the superiority of this bulk free energy over the White Bear one.

preprint2012arXiv

Collisional rates for the inelastic Maxwell model: application to the divergence of anisotropic high-order velocity moments in the homogeneous cooling state

The collisional rates associated with the isotropic velocity moments $<V^{2r}>$ and the anisotropic moments $<V^{2r}V_i>$ and $<V^{2r}(V_iV_j-d^{-1}V^2δ_{ij})>$ are exactly derived in the case of the inelastic Maxwell model as functions of the exponent $r$, the coefficient of restitution $α$, and the dimensionality $d$. The results are applied to the evolution of the moments in the homogeneous free cooling state. It is found that, at a given value of $α$, not only the isotropic moments of a degree higher than a certain value diverge but also the anisotropic moments do. This implies that, while the scaled distribution function has been proven in the literature to converge to the isotropic self-similar solution in well-defined mathematical terms, nonzero initial anisotropic moments do not decay with time. On the other hand, our results show that the ratio between an anisotropic moment and the isotropic moment of the same degree tends to zero.

preprint2012arXiv

Fourth virial coefficients of asymmetric nonadditive hard-disc 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\leq q\leq 0.95$) and nonadditivity ($-0.5\leq Δ\leq 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 $Δ=\pm 0.2$.

preprint2012arXiv

Note: An exact scaling relation for truncatable free energies of polydisperse hard-sphere mixtures

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.

preprint2012arXiv

Unsteady non-Newtonian hydrodynamics in granular gases

The temporal evolution of a dilute granular gas, both in a compressible flow (uniform longitudinal flow) and in an incompressible flow (uniform shear flow), is investigated by means of the direct simulation Monte Carlo method to solve the Boltzmann equation. Emphasis is laid on the identification of a first "kinetic" stage (where the physical properties are strongly dependent on the initial state) subsequently followed by an unsteady "hydrodynamic" stage (where the momentum fluxes are well-defined non-Newtonian functions of the rate of strain). The simulation data are seen to support this two-stage scenario. Furthermore, the rheological functions obtained from simulation are well described by an approximate analytical solution of a model kinetic equation.

preprint2011arXiv

A Bhatnagar-Gross-Krook-like Model Kinetic Equation for a Granular Gas of Inelastic Rough Hard Spheres

The Boltzmann collision operator for a dilute granular gas of inelastic rough hard spheres is much more intricate than its counterpart for inelastic smooth spheres. Now the one-body distribution function depends not only on the translational velocity of the center of mass but also on the angular velocity of the particle. Moreover, the collision rules couple both velocities, involving not only the coefficient of normal restitution but also the coefficient of tangential restitution. The aim of this paper is to propose an extension to inelastic rough particles of a Bhatnagar-Gross-Krook-like kinetic model previously proposed for inelastic smooth particles. The Boltzmann collision operator is replaced by the sum of three terms representing: (i) the relaxation to a two-temperature local equilibrium distribution, (ii) the action of a nonconservative drag force proportional to the peculiar velocity, and (iii) the action of a nonconservative torque equal to a linear combination of the angular velocity and its mean value. The three coefficients in the force and torque are fixed to reproduce the Boltzmann collisional rates of change of the mean angular velocity and of the two granular temperatures (translational and rotational). A simpler version of the model is also constructed in the form of two coupled kinetic equations for the translational and rotational velocity distributions. The kinetic model is applied to the simple shear flow steady state and the combined influence of the two coefficients of restitution on the shear and normal stresses and on the translational velocity distribution function is analyzed.

preprint2011arXiv

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.

preprint2011arXiv

Exact solution of the Percus-Yevick integral equation for fluid mixtures of hard hyperspheres

Structural and thermodynamic properties of multicomponent hard-sphere fluids at odd dimensions have recently been derived in the framework of the rational function approximation (RFA) [Rohrmann and Santos, Phys. Rev. E \textbf{83}, 011201 (2011)]. It is demonstrated here that the RFA technique yields the exact solution of the Percus-Yevick (PY) closure to the Ornstein-Zernike (OZ) equation for binary mixtures at arbitrary odd dimensions. The proof relies mainly on the Fourier transforms $\hat{c}_{ij}(k)$ of the direct correlation functions defined by the OZ relation. From the analysis of the poles of $\hat{c}_{ij}(k)$ we show that the direct correlation functions evaluated by the RFA method vanish outside the hard core, as required by the PY theory.

preprint2011arXiv

Homogeneous Free Cooling State in Binary Granular Fluids of Inelastic Rough Hard Spheres

In a recent paper [A. Santos, G. M. Kremer, and V. Garzó, \emph{Prog. Theor. Phys. Suppl.} \textbf{184}, 31-48 (2010)] the collisional energy production rates associated with the translational and rotational granular temperatures in a granular fluid mixture of inelastic rough hard spheres have been derived. In the present paper the energy production rates are explicitly decomposed into equipartition rates (tending to make all the temperatures equal) plus genuine cooling rates (reflecting the collisional dissipation of energy). Next the homogeneous free cooling state of a binary mixture is analyzed, with special emphasis on the quasi-smooth limit. A previously reported singular behavior (according to which a vanishingly small amount of roughness has a finite effect, with respect to the perfectly smooth case, on the asymptotic long-time translational/translational temperature ratio) is further elaborated. Moreover, the study of the time evolution of the temperature ratios shows that this dramatic influence of roughness already appears in the transient regime for times comparable to the relaxation time of perfectly smooth spheres.

preprint2011arXiv

Multicomponent fluids of hard hyperspheres in odd dimensions

Mixtures of hard hyperspheres in odd space dimensionalities are studied with an analytical approximation method. This technique is based on the so-called Rational Function Approximation and provides a procedure for evaluating equations of state, structure factors, radial distribution functions, and direct correlations functions of additive mixtures of hard hyperspheres with any number of components and in arbitrary odd-dimension space. The method gives the exact solution of the Ornstein--Zernike equation coupled with the Percus--Yevick closure, thus extending to arbitrary odd dimension the solution for hard-sphere mixtures [J. L. Lebowitz, Phys.\ Rev.\ \textbf{133}, 895 (1964)]. Explicit evaluations for binary mixtures in five dimensions are performed. The results are compared with computer simulations and a good agreement is found.

preprint2011arXiv

Non-Newtonian Couette-Poiseuille flow of a dilute gas

The steady state of a dilute gas enclosed between two infinite parallel plates in relative motion and under the action of a uniform body force parallel to the plates is considered. The Bhatnagar-Gross-Krook model kinetic equation is analytically solved for this Couette-Poiseuille flow to first order in the force and for arbitrary values of the Knudsen number associated with the shear rate. This allows us to investigate the influence of the external force on the non-Newtonian properties of the Couette flow. Moreover, the Couette-Poiseuille flow is analyzed when the shear-rate Knudsen number and the scaled force are of the same order and terms up to second order are retained. In this way, the transition from the bimodal temperature profile characteristic of the pure force-driven Poiseuille flow to the parabolic profile characteristic of the pure Couette flow through several intermediate stages in the Couette-Poiseuille flow are described. A critical comparison with the Navier-Stokes solution of the problem is carried out.

preprint2011arXiv

Nonadditive hard-sphere fluid mixtures: A simple analytical theory

We construct a non-perturbative fully analytical approximation for the thermodynamics and the structure of nonadditive hard-sphere fluid mixtures. The method essentially lies in a heuristic extension of the Percus-Yevick solution for additive hard spheres. Extensive comparison with Monte Carlo simulation data shows a generally good agreement, especially in the case of like-like radial distribution functions.

preprint2011arXiv

On the relation between virial coefficients and the close-packing of hard disks and hard spheres

The question of whether the known virial coefficients are enough to determine the packing fraction $η_\infty$ 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 $η_\infty$ is equal to the maximum packing fraction corresponding to an ordered crystalline structure.

preprint2011arXiv

Sonine approximation for collisional moments of granular gases of inelastic rough spheres

We consider a dilute granular gas of hard spheres colliding inelastically with coefficients of normal and tangential restitution $α$ and $β$, respectively. The basic quantities characterizing the distribution function $f(\mathbf{v},\bmω)$ of linear ($\mathbf{v}$) and angular ($\bmω$) velocities are the second-degree moments defining the translational ($T^\text{tr}$) and rotational ($T^\text{rot}$) temperatures. The deviation of $f$ from the Maxwellian distribution parameterized by $T^\text{tr}$ and $T^\text{rot}$ can be measured by the cumulants associated with the fourth-degree velocity moments. The main objective of this paper is the evaluation of the collisional rates of change of these second- and fourth-degree moments by means of a Sonine approximation. The results are subsequently applied to the computation of the temperature ratio $T^\text{rot}/T^\text{tr}$ and the cumulants of two paradigmatic states: the homogeneous cooling state and the homogeneous steady state driven by a white-noise stochastic thermostat. It is found in both cases that the Maxwellian approximation for the temperature ratio does not deviate much from the Sonine prediction. On the other hand, non-Maxwellian properties measured by the cumulants cannot be ignored, especially in the homogeneous cooling state for medium and small roughness. In that state, moreover, the cumulant directly related to the translational velocity differs in the quasi-smooth limit $β\to -1$ from that of pure smooth spheres ($β=-1$). This singular behavior is directly related to the unsteady character of the homogeneous cooling state and thus it is absent in the stochastic thermostat case.

preprint2010arXiv

Energy production rates in fluid mixtures of inelastic rough hard spheres

The aim of this work is to explore the combined effect of polydispersity and roughness on the partial energy production rates and on the total cooling rate of a granular fluid mixture. We consider a mixture of inelastic rough hard spheres of different number densities, masses, diameters, moments of inertia, and mutual coefficients of normal and tangential restitution. Starting from the first equation of the BBGKY hierarchy, the collisional energy production rates associated with the translational and rotational temperatures ($T_i^\text{tr}$ and $T_i^\text{rot}$) are expressed in terms of two-body average values. Next, those average values are estimated by assuming a velocity distribution function based on maximum-entropy arguments, allowing us to express the energy production rates and the total cooling rate in terms of the partial temperatures and the parameters of the mixture. Finally, the results are applied to the homogeneous cooling state of a binary mixture and the influence of inelasticity and roughness on the temperature ratios $T_1^\text{tr}/T_1^\text{rot}$, $T_2^\text{tr}/T_1^\text{tr}$, and $T_2^\text{rot}/T_1^\text{rot}$ is analyzed.

preprint2010arXiv

Molecular dynamics simulation study of self-diffusion for penetrable-sphere model fluids

Molecular dynamics simulations are carried out to investigate the diffusion behavior of penetrable-sphere model fluids characterized by a finite energy barrier $ε$. The self-diffusion coefficient is evaluated from the time-dependent velocity autocorrelation function and mean-square displacement. Detailed insights into the cluster formation for penetrable spheres are gained from the Enskog factor, the effective particle volume fraction, the mean free path, and the collision frequency for both the soft-type penetrable and the hard-type reflective collisions. The simulation data are compared to theoretical predictions from the Boltzmann kinetic equation and from a simple extension to finite $ε$ of the Enskog prediction for impenetrable hard spheres ($ε\to\infty$). A reasonable agreement between theoretical and simulation results is found in the cases of $ε^*\equiv ε/k_BT=0.2$, $0.5$, and $1.0$. However, for dense systems (packing fraction $ϕ>0.6$) with a highly repulsive energy barrier ($ε^*=3.0$), a poorer agreement was observed due to metastable static effects of clustering formation and dynamic effects of correlated collision processes among these cluster-forming particles.

preprint2010arXiv

Non-Newtonian granular hydrodynamics. What do the inelastic simple shear flow and the elastic Fourier flow have in common?

We describe a special class of steady Couette flows in dilute granular gases admitting a non-Newtonian hydrodynamic description for strong dissipation. The class occurs when viscous heating exactly balances inelastic cooling, resulting in a uniform heat flux. It includes the Fourier flow of ordinary gases and the simple or uniform shear flow (USF) of granular gases as special cases. The rheological functions have the same values as in the USF and generalized thermal conductivity coefficients can be identified. These points are confirmed by molecular dynamics simulations, Monte Carlo simulations of the Boltzmann equation, and analytical results from Grad's 13-moment method.

preprint2010arXiv

Simple relationship between the virial-route hypernetted-chain and the compressibility-route Percus--Yevick values of the fourth virial coefficient

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 $B_4$ 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 $B_4$ 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.

preprint2010arXiv

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.

preprint2001arXiv

Critical Behavior of a Heavy Particle in a Granular Fluid

Behavior analogous to a second order phase transition is observed for the homogeneous cooling state of a heavy impurity particle in a granular fluid. The order parameter $ϕ$ is the ratio of impurity mean square velocity to that of the fluid, with a conjugate field $h$ proportional to the mass ratio. A parameter $β$, measuring the fluid cooling rate relative to the impurity--fluid collision rate, is the analogue of the inverse temperature. For $β<1$ the fluid is ``normal'' with $ϕ=0$ at $h=0$, as in the case of a system with elastic collisions. For $β>1$ an ``ordered'' state with $ϕ\neq 0$ occurs at $h=0$, representing an extreme breakdown of equipartition. Critical slowing and qualitative changes in the velocity distribution function for the impurity particle near the transition are noted

preprint2000arXiv

Non-Newtonian Poiseuille flow of a gas in a pipe

The Bhatnagar-Gross-Krook kinetic model of the Boltzmann equation is solved for the steady cylindrical Poiseuille flow fed by a constant gravity field. The solution is obtained as a perturbation expansion in powers of the field (through fourth order) and for a general class of repulsive potentials. The results, which are hardly sensitive to the interaction potential, suggest that the expansion is only asymptotic. A critical comparison with the profiles predicted by the Navier-Stokes equations shows that the latter fail over distances comparable to the mean free path. In particular, while the Navier-Stokes description predicts a monotonically decreasing temperature as one moves apart from the cylinder axis, the kinetic theory description shows that the temperature has a local minimum at the axis and reaches a maximum value at a distance of the order of the mean free path. Within that distance, the radial heat flows from the colder to the hotter points, in contrast to what is expected from the Fourier law. Furthermore, a longitudinal component of the heat flux exists in the absence of gradients along the longitudinal direction. Non-Newtonian effects, such as a non uniform hydrostatic pressure and normal stress differences, are also present.