Riccardo Fantoni

Penetrable square-well fluids: Exact results in one dimension

We introduce a model of attractive penetrable spheres by adding a short-range attractive square well outside a penetrable core, and we provide a detailed analysis of structural and thermodynamical properties in one dimension using the exact impenetrable counterpart as a starting point. The model is expected to describe star polymers in regimes of good and moderate solvent under dilute conditions. We derive the exact coefficients of a low-density expansion up to second order for the radial distribution function and up to fourth order in the virial expansion. These exact results are used as a benchmark to test the reliability of approximate theories (Percus-Yevick and hypernetted chain). Notwithstanding the lack of an exact solution for arbitrary densities, our results are expected to be rather precise within a wide range of temperatures and densities. A detailed analysis of some limiting cases is carried out. In particular, we provide a complete solution of the sticky penetrable-sphere model in one dimension up to the same order in density. The issue of Ruelles thermodynamics stability is analyzed and the region of a well-defined thermodynamic limit is identified.

Multicomponent adhesive hard sphere models and short-ranged attractive interactions in colloidal or micellar solutions.

We investigate the dependence of the stickiness parameters tij=1/(12tauij)--where the tauij are the conventional Baxter parameters--on the solute diameters sigmai and sigmaj in multicomponent sticky hard sphere (SHS) models for fluid mixtures of mesoscopic neutral particles. A variety of simple but realistic interaction potentials, utilized in the literature to model short-ranged attractions present in real solutions of colloids or reverse micelles, is reviewed. We consider: (i) van der Waals attractions, (ii) hard-sphere-depletion forces, (iii) polymer-coated colloids, and (iv) solvation effects (in particular hydrophobic bonding and attractions between reverse micelles of water-in-oil microemulsions). We map each of these potentials onto an equivalent SHS model by requiring the equality of the second virial coefficients. The main finding is that, for most of the potentials considered, the size-dependence of tij(T,sigmai,sigmaj) can be approximated by essentially the same expression, i.e., a simple polynomial in the variable sigmaisigmaj/sigmaij2, with coefficients depending on the temperature T, or--for depletion interactions--on the packing fraction eta0 of the depletant particles.

Thermodynamic instabilities of a binary mixture of sticky hard spheres

The thermodynamic instabilities of a binary mixture of sticky hard spheres (SHS) in the modified mean spherical approximation (mMSA) and the Percus-Yevick (PY) approximation are investigated using an approach devised by Chen and Forstmann [corrected] [J. Chem. Phys. [corrected] 97, 3696 (1992)]. This scheme hinges on a diagonalization of the matrix of second functional derivatives of the grand canonical potential with respect to the particle density fluctuations. The zeroes of the smallest eigenvalue and the direction of the relative eigenvector characterize the instability uniquely. We explicitly compute three different classes of examples. For a symmetrical binary mixture, analytical calculations, both for mMSA and for PY, predict that when the strength of adhesiveness between like particles is smaller than the one between unlike particles, only a pure condensation spinodal exists; in the opposite regime, a pure demixing spinodal appears at high densities. We then compare the mMSA and PY results for a mixture where like particles interact as hard spheres (HS) and unlike particles as SHS, and for a mixture of HS in a SHS fluid. In these cases, even though the mMSA and PY spinodals are quantitatively and qualitatively very different from each other, we prove that they have the same kind of instabilities. Finally, we study the mMSA solution for five different mixtures obtained by setting the stickiness parameters equal to five different functions of the hard sphere diameters. We find that four of the five mixtures exhibit very different type of instabilities. Our results are expected to provide a further step toward a more thoughtful application of SHS models to colloidal fluids.

Penetrable-square-well fluids: Analytical study and Monte Carlo simulations

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.

Phase behavior of weakly polydisperse sticky hard spheres: perturbation theory for the Percus-Yevick solution.

We study the effects of size polydispersity on the gas-liquid phase behavior of mixtures of sticky hard spheres. To achieve this, the system of coupled quadratic equations for the contact values of the partial cavity functions of the Percus-Yevick solution [R. J. Baxter, J. Chem. Phys. 49, 2770 (1968)] is solved within a perturbation expansion in the polydispersity, i.e., the normalized width of the size distribution. This allows us to make predictions for various thermodynamic quantities which can be tested against numerical simulations and experiments. In particular, we determine the leading order effects of size polydispersity on the cloud curve delimiting the region of two-phase coexistence and on the associated shadow curve; we also study the extent of size fractionation between the coexisting phases. Different choices for the size dependence of the adhesion strengths are examined carefully; the Asakura-Oosawa model [J. Chem. Phys. 22, 1255 (1954)] of a mixture of polydisperse colloids and small polymers is studied as a specific example.

A cluster theory for a Janus fluid

Abstract Recent Monte Carlo simulations on the Kern and Frenkel model of a Janus fluid have revealed that in the vapour phase there is the formation of preferred clusters made up of a well-defined number of particles: the micelles and the vesicles. A cluster theory is developed to approximate the exact clustering properties stemming from the simulations. It is shown that the theory is able to reproduce the micellisation phenomenon.

Phase diagram of the penetrable-square-well model

We study a system formed by soft colloidal spheres attracting each other via a square-well potential, using extensive Monte Carlo simulations of various nature. The softness is implemented through a reduction of the infinite part of the repulsive potential to a finite one. For sufficiently low values of the penetrability parameter we find the system to be Ruelle stable with square-well–like behavior. For high values of the penetrability the system is thermodynamically unstable and collapses into an isolated blob formed by a few clusters each containing many overlapping particles. For intermediate values of the penetrability the system has a rich phase diagram with a partial lack of thermodynamic consistency.

Thermodynamic consistency of energy and virial routes: An exact proof within the linearized Debye–Hückel theory

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.

Wertheim perturbation theory: thermodynamics and structure of patchy colloids

We critically discuss the application of the Wertheim’s theory to classes of complex associating fluids that can be today engineered in the laboratory as patchy colloids and to the prediction of their peculiar gas–liquid phase diagrams. Our systematic study, stemming from perturbative version of the theory, allows us to show that, even at the simplest level of approximation for the inter-cluster correlations, the theory is still able to provide a consistent and stable picture of the behaviour of interesting models of self-assembling colloidal suspension. We extend the analysis of a few cases of patchy systems recently introduced in the literature. In particular, we discuss for the first time in detail the consistency of the structural description underlying the perturbative approach and we are able to prove a consistency relationship between the valence as obtained from thermodynamics and from the structure for the one-site case. A simple analytical expression for the structure factor is proposed.

Monte Carlo simulation of the nonadditive restricted primitive model of ionic fluids: phase diagram and clustering.

We report an accurate Monte Carlo calculation of the phase diagram and clustering properties of the restricted primitive model with nonadditive hard-sphere diameters. At high density the positively nonadditive fluid shows more clustering than in the additive model and the negatively nonadditive fluid shows less clustering than in the additive model; at low density the reverse scenario appears. A negative nonadditivity tends to favor the formation of neutrally charged clusters starting from the dipole. A positive nonadditivity favors the pairing of like ions at high density. The critical point of the gas-liquid phase transition moves at higher temperatures and higher densities for a negative nonadditivity and at lower temperatures and lower densities for a positive nonadditivity. The law of corresponding states does not seem to hold strictly. Our results can be used to interpret recent experimental works on room temperature ionic liquids.