When mixtures of hard-sphere-like colloids do not behave as mixtures of hard spheres

Glass transition and reversible gelation in asymmetric binary mixtures: A study by mode coupling theory and molecular dynamics

The glass transition and the binodals of asymmetric binary mixtures are investigated from the effective fluid approach in the mode coupling theory and by molecular dynamics. Motivated by previous theoretical predictions, the hard-sphere mixture and the Asakura-Oosawa models are used to analyze experimental results from the literature, relative to polystyrene spheres mixed either with linear polymers or with dense microgel particles. In agreement with the experimental observations, the specificity of the depletant particles is shown to favor lower density gels. It further favors equilibrium gelation by reducing also the tendency of the system to phase separate. These results are confirmed by a phenomenological modification of the mode coupling theory in which the vertex functions are computed at an effective density lower than the actual one. A model effective potential in asymmetric mixtures of hard particles is used to further check this phenomenological modification against molecular dynamics simulation.

Glass transition in hard-core fluids and beyond, using an effective static structure in the mode coupling theory

The dynamical arrest in classical fluids is studied using a simple modification of the mode coupling theory (MCT) aimed at correcting its overestimation of the tendency to glass formation while preserving its overall structure. As in previous attempts, the modification is based on the idea of tempering the static pair correlations used as input. It is implemented in this work by computing the static structure at a different state point than the one used to solve the MCT equation for the intermediate scattering function, using the pure hard-sphere glass for calibration. The location of the glass transition predicted from this modification is found to agree with simulations data for a variety of systems --pure fluids and mixtures with either purely repulsive interaction potentials or ones with attractive contributions. Besides improving the predictions in the long-time limit, and so reducing the non-ergodicity domain, the same modification works as well for the time-dependent correlators.

Self-assembly and crystallisation of indented colloids at a planar wall

We report experimental and simulation studies of the structure of a monolayer of indented ("lock and key") colloids, on a planar surface. On adding a non-absorbing polymer with prescribed radius and volume fraction, depletion interactions are induced between the colloids, with controlled range and strength. For spherical particles, this leads to crystallisation, but the indented colloids crystallise less easily than spheres, in both simulation and experiment. Nevertheless, simulations show that indented colloids do form plastic (rotator) crystals. We discuss the conditions under which this occurs, and the possibilities of lower-symmetry crystal states. We also comment on the kinetic accessibility of these states.

Bond lifetime and diffusion coefficient in colloids with short-range interactions

We use molecular dynamics simulations to study the influence of short-range structures in the interaction potential between hard-sphere-like colloidal particles. Starting from model potentials and effective potentials in binary mixtures computed from the Ornstein-Zernike equations, we investigate the influence of the range and strength of a possible tail beyond the usual core repulsion or the presence of repulsive barriers. The diffusion coefficient and mean "bond" lifetimes are used as indicators of the effect of this structure on the dynamics. The existence of correlations between the variations of these quantities with the physical parameters is discussed to assess the interpretation of dynamics slowing down in terms of long-lived bonds. We also discuss the question of a universal behaviour determined by the second virial coefficient B ((2)) and the interplay of attraction and repulsion. While the diffusion coefficient follows the B ((2)) law for purely attractive tails, this is no longer true in the presence of repulsive barriers. Furthermore, the bond lifetime shows a dependence on the physical parameters that differs from that of the diffusion coefficient. This raises the question of the precise role of bonds on the dynamics slowing down in colloidal gels.

How soft repulsion enhances the depletion mechanism

We investigate binary mixtures of large colloids interacting through soft potentials with small, ideal depletants. We show that softness has a dramatic effect on the resulting colloid-colloid effective potential when the depletant-to-colloid size ratio q is small, with significant consequences on the colloidal phase behaviour. We provide an exact relationship that allows us to obtain the effective pair potential for any type of colloid-depletant interaction in the case of ideal depletants, without having to rely on complicated and expensive full-mixture simulations. We also show that soft repulsion among depletants further enhances the tendency of colloids to aggregate. Our theoretical and numerical results demonstrate that--in the limit of small q--soft mixtures cannot be mapped onto hard systems and hence soft depletion is not a mere extension of the widely used Asakura-Oosawa potential.

Free-energy landscape for cage breaking of three hard disks

We investigate cage breaking in dense hard-disk systems using a model of three Brownian disks confined within a circular corral. This system has a six-dimensional configuration space, but can be equivalently thought to explore a symmetric one-dimensional free-energy landscape containing two energy minima separated by an energy barrier. The exact free-energy landscape can be calculated as a function of system size by a direct enumeration of states. Results of simulations show the average time between cage breaking events follows an Arrhenius scaling when the energy barrier is large. We also discuss some of the consequences of using a one-dimensional representation to understand dynamics through a multidimensional space, such as diffusion acquiring spatial dependence and discontinuities in spatial derivatives of free energy.

Depletion potentials in highly size-asymmetric binary hard-sphere mixtures: comparison of simulation results with theory

We report a detailed study, using state-of-the-art simulation and theoretical methods, of the effective (depletion) potential between a pair of big hard spheres immersed in a reservoir of much smaller hard spheres, the size disparity being measured by the ratio of diameters q ≡ σ(s)/σ(b). Small particles are treated grand canonically, their influence being parameterized in terms of their packing fraction in the reservoir η(s)(r). Two Monte Carlo simulation schemes--the geometrical cluster algorithm, and staged particle insertion--are deployed to obtain accurate depletion potentials for a number of combinations of q ≤ 0.1 and η(s)(r). After applying corrections for simulation finite-size effects, the depletion potentials are compared with the prediction of new density functional theory (DFT) calculations based on the insertion trick using the Rosenfeld functional and several subsequent modifications. While agreement between the DFT and simulation is generally good, significant discrepancies are evident at the largest reservoir packing fraction accessible to our simulation methods, namely, η(s)(r) = 0.35. These discrepancies are, however, small compared to those between simulation and the much poorer predictions of the Derjaguin approximation at this η(s)(r). The recently proposed morphometric approximation performs better than Derjaguin but is somewhat poorer than DFT for the size ratios and small-sphere packing fractions that we consider. The effective potentials from simulation, DFT, and the morphometric approximation were used to compute the second virial coefficient B(2) as a function of η(s)(r). Comparison of the results enables an assessment of the extent to which DFT can be expected to correctly predict the propensity toward fluid-fluid phase separation in additive binary hard-sphere mixtures with q ≤ 0.1. In all, the new simulation results provide a fully quantitative benchmark for assessing the relative accuracy of theoretical approaches for calculating depletion potentials in highly size-asymmetric mixtures.

Gelation and phase coexistence in colloidal suspensions with short-range forces: generic behavior versus specificity

The interplay between physical gelation and equilibrium phase transitions in asymmetric binary mixtures is analyzed from the effective fluid approach, in which the big particles interact via a short-range effective attraction beyond the core due to the depletion mechanism. The question of the universality of the scenario for dynamical arrest is then addressed. The comparison of the phase diagrams of the hard-sphere mixture and the Asakura-Oosawa models at various size ratios shows that strong specificity is observed for nonideal depletants. In particular, equilibrium gelation, without the competition with fluid-fluid transition is possible in mixtures of hard-sphere colloids. This is interpreted from the specificities of the effective potential, such as its oscillatory behavior and its complex variation with the physical parameters. The consequences on the dynamical arrest and the fluid-fluid transition are then investigated by considering in particular the role of the well at contact and the first repulsive barrier. This is done for the actual effective potential in the hard-sphere mixture and for a square well and shoulder model, which allows a separate discussion of the role of the different parameters, in particular on the localization length and the escape time. This study is next extended to mixtures of "hard-sphere-like" colloids with residual interactions. It confirms the trends relative to equilibrium gelation and illustrates a diversity of the phase behavior well beyond the scenarios expected from simple models.

Equilibrium route to colloidal gelation: mixtures of hard-sphere-like colloids

The binodals and the nonergodicity lines of a binary mixture of hard-sphere-like particles with a large size ratio are computed for studying the interplay between dynamic arrest and phase separation in depletion-driven colloidal mixtures. Contrary to the case of hard core plus short-range effective attraction, physical gelation without competition with the fluid-phase separation can occur in such mixtures. This behavior due to the oscillations in the depletion potential should concern all simple mixtures with a nonideal depletant, justifying further studies of their dynamic properties.

Equilibrium and glassy states of the Asakura-Oosawa and binary hard sphere mixtures: effective fluid approach

Motivated by recent experimental results on model binary colloidal mixtures, especially for the glass transition, we investigate the phase diagram of two models of asymmetric binary mixtures: the hard sphere and the Asakura-Oosawa mixtures. This includes the binodals and the glass transition line, computed in the effective one-component representation using the corresponding potentials of mean force at infinite dilution. The reference hypernetted chain approximation is used for computing the static properties and the glass transition line is computed in the mode coupling approximation. The similarities and the differences between the two models are discussed for different size ratios. It is shown that while both models follow a universal behavior at large asymmetry, the hard sphere mixture model leads to more original results at moderate size ratio. These results show that a modeling beyond generic effective potentials might be necessary for an appropriate description of the complete phase diagram.

Phase transitions in highly asymmetric binary hard-sphere fluids: Fluid-fluid binodal from a two-component mixture theory

Fluid-fluid binodals of binary hard-sphere mixtures are computed from the recently proposed fundamental measure functional-mean spherical approximation closure of the two-component Ornstein-Zernike equation. The results, especially in the dense fluid region that was not accessible by previous theoretical methods, are compared with the corresponding ones for the one-component fluid of big spheres with effective potential obtained from the same closure. The general trends are those expected for hard-sphere potentials but small difference are detectable. The overall agreement found validates the equivalence of the two descriptions for size ratios R = 8.5 or greater.

Influence of polydispersity on the critical parameters of an effective-potential model for asymmetric hard-sphere mixtures

We report a Monte Carlo simulation study of the properties of highly asymmetric binary hard-sphere mixtures. This system is treated within an effective fluid approximation in which the large particles interact through a depletion potential [R. Roth, Phys. Rev. E 62 5360 (2000)] designed to capture the effects of a virtual sea of small particles. We generalize this depletion potential to include the effects of explicit size dispersity in the large particles and consider the case in which the particle diameters are distributed according to a Schulz form having a degree of polydispersity 14%. The resulting alteration (with respect to the monodisperse limit) of the metastable fluid-fluid critical point parameters is determined for two values of the ratio of the diameters of the small and large particles: q(triple bond)sigma(s)/(-)sigma(b)=0.1 and q=0.05. We find that the inclusion of polydispersity moves the critical point to lower reservoir volume fractions of the small particles and high volume fractions of the large ones. The estimated critical point parameters are found to be in good agreement with those predicted by a generalized corresponding states argument which provides a link to the known critical adhesion parameter of the adhesive hard-sphere model. Finite-size scaling estimates of the cluster percolation line in the one phase fluid region indicate that inclusion of polydispersity moves the critical point deeper into the percolating regime. This suggests that phase separation is more likely to be preempted by dynamical arrest in polydisperse systems.

Potential of mean force in confined colloids: Integral equations with fundamental measure bridge functions

The potential of mean force for uncharged macroparticles suspended in a fluid confined by a wall or a narrow pore is computed for solvent-wall and solvent-macroparticle interactions with attractive forces. Bridge functions taken from Rosenfeld's density-functional theory are used in the reference hypernetted chain closure of the Ornstein-Zernike integral equations. The quality of this closure is assessed by comparison with simulation. As an illustration, the role of solvation forces is investigated. When the "residual" attractive tails are given a range appropriate to "hard sphere-like" colloids, the unexpected role of solvation forces previously observed in bulk colloids is confirmed in the confinement situation.