Hard-sphere fluid adsorbed in an annular wedge: the depletion force of hard-body colloidal physics

Many important issues of colloidal physics can be expressed in the context of inhomogeneous fluid phenomena. When two large colloids approach one another in solvent, they interact at least partly by the response of the solvent to finding itself adsorbed in the annular wedge formed between the two colloids. At shortest range, this fluid mediated interaction is known as the depletion force/interaction because solvent is squeezed out of the wedge when the colloids approach closer than the diameter of a solvent molecule. An equivalent situation arises when a single colloid approaches a substrate/wall. Accurate treatment of this interaction is essential for any theory developed to model the phase diagrams of homogeneous and inhomogeneous colloidal systems. The aim of our paper is a test of whether or not we possess sufficient knowledge of statistical mechanics that can be trusted when applied to systems of large size asymmetry and the depletion force in particular. When the colloid particles are much larger than a solvent diameter, the depletion force is dominated by the effective two-body interaction experienced by a pair of solvated colloids. This low concentration limit of the depletion force has therefore received considerable attention. One route, which can be rigorously based on statistical mechanical sum rules, leads to an analytic result for the depletion force when evaluated by a key theoretical tool of colloidal science known as the Derjaguin approximation. A rival approach has been based on the assumption that modern density functional theories (DFT) can be trusted for systems of large size asymmetry. Unfortunately, these two theoretical predictions differ qualitatively for hard sphere models, as soon as the solvent density is higher than about 23 that at freezing. Recent theoretical attempts to understand this dramatic disagreement have led to the proposal that the Derjaguin and DFT routes represent opposite limiting behavior, for very large size asymmetry and molecular sized mixtures, respectively. This proposal implies that nanocolloidal systems lie in between the two limits, so that the depletion force no longer scales linearly with the colloid radius. That is, by decreasing the size ratio from mesoscopic to molecular sized solutes, one moves smoothly between the Derjaguin and the DFT predictions for the depletion force scaled by the colloid radius. We describe the results of a simulation study designed specifically as a test of compatibility with this complex scenario. Grand canonical simulation procedures applied to hard-sphere fluid adsorbed in a series of annular wedges, representing the depletion regime of hard-body colloidal physics, confirm that neither the Derjaguin approximation, nor advanced formulations of DFT, apply at moderate to high solvent density when the geometry is appropriate to nanosized colloids. Our simulations also allow us to report structural characteristics of hard-body solvent adsorbed in hard annular wedges. Both these aspects are key ingredients in the proposal that unifies the disparate predictions, via the introduction of new physics. Our data are consistent with this proposed physics, although as yet limited to a single colloidal size asymmetry.