Grand canonical simulation of phase behaviour in highly size-asymmetrical binary fluids

Multilevel simulation of hard-sphere mixtures

We present a multilevel Monte Carlo simulation method for analysing multi-scale physical systems via a hierarchy of coarse-grained representations, to obtain numerically-exact results, at the most detailed level. We apply the method to a mixture of size-asymmetric hard spheres, in the grand canonical ensemble. A three-level version of the method is compared with a previously-studied two-level version. The extra level interpolates between the full mixture and a coarse-grained description where only the large particles are present -- this is achieved by restricting the small particles to regions close to the large ones. The three-level method improves the performance of the estimator, at fixed computational cost. We analyse the asymptotic variance of the estimator, and discuss the mechanisms for the improved performance.

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.

Critical casimir forces and colloidal phase transitions in a near-critical solvent: a simple model reveals a rich phase diagram

From experimental studies, it is well known that colloidal particles suspended in a near-critical binary solvent exhibit interesting aggregation phenomena, often associated with colloidal phase transitions and assumed to be driven by long-ranged solvent-mediated (SM) interactions (critical Casimir forces), set by the (diverging) correlation length of the solvent. We present the first simulation and theoretical study of an explicit model of a ternary mixture that mimics this situation. Both the effective SM pair interactions and the full ternary phase diagram are determined for Brownian disks suspended in an explicit two-dimensional supercritical binary liquid mixture. Gas-liquid and fluid-solid transitions are observed in a region that extends well away from criticality of the solvent reservoir. We discuss to what extent an effective pair-potential description can account for the phase behavior we observe. Our study provides a fresh perspective on how proximity to the critical point of the solvent reservoir might influence colloidal self-assembly.

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.