On measuring colloidal volume fractions Wilson C. K. Poon Scottish Universities Physics Alliance (SUPA) and The School of Physics, University of Edinburgh, Kings Buildings, Mayfield Road, Edinburgh EH9 3JZ, U. K. Eric R. Weeks Department of Physics, Emory University, Atlanta, GA 30322 U.S.A. C. Patrick Royall School of Chemistry, University of Bristol, Bristol, BS8 1TS, U.K Hard-sphere colloids are popular as models for testing fundamental theories in condensed matter and statistical physics, from crystal nucleation to the glass transition. A single parameter, the volume fraction (φ), characterizes an ideal, monodisperse hard-sphere suspension. In comparing experiments with theories and simulation, researchers to date have paid little attention to likely uncertainties in experimentally-quoted φ values. We critically review the experimental measurement of φ in hard-sphere colloids, and show that while statistical uncertainties in comparing relative values − of φ can be as low as 10 4, systematic errors of 3-6% are probably unavoidable. The consequences of this are illustrated by way of a case study comparing literature data sets on hard-sphere viscosity and diffusion. I. INTRODUCTION spheres is controlled by one parameter, the volume frac- tion φ, i.e. the fraction of the total volume V that is filled by N spheres, each of radius a, Excluded volume effects dominate the behaviour of liquids around the triple point [1] and play a key role 4 3 N φ = πa . (1) in structuring crystalline [2] and amorphous solids [3]. 3 V Thus, hard spheres have long functioned as a reference Since φ is precisely known in theory or simulations, a system for theoretical and simulational studies of con- comparison with experiments is straightforward provided densed matter. In 1986, Pusey and van Megen demon- that this quantity is also accurately measurable for real strated [4] that suspensions of sterically-stabilised poly- suspensions. Much of the literature has indeed pro- methylmethacrylate (PMMA) colloids showed nearly- ceeded on this basis, assuming that φ is unproblemati- perfect hard-sphere equilibrium phase behaviour, under- cally known from experiments. going a first-order phase transition from a fluid to a crys- However, as Pusey and van Megen pointed out in a talline state at concentrations around those predicted symposium article [17] following their Nature paper,[4] some time ago by computer simulations for hard spher- the experimental determination of φ is emphatically not ical particles [5, 6]. Soon afterward, the same authors unproblematic, because: (1) no real colloid is truly showed [7] that PMMA colloids underwent a glass tran- ‘hard’, since there is always some softness in the inter- sition at even higher concentrations. Subsequently, in particle potential; and (2) real colloids always have a a 1991 review, Pusey enunciated the ‘colloids as atoms’ finite size distribution, i.e. they are polydisperse. Thus, paradigm [8] — Brownian suspensions can be used as Pusey and van Megen calculated an experimental ‘ef- ‘test tube simulations’ of many generic condensed mat- fective’ hard-sphere volume fraction φ , and found that ter phenomena such as crystallization and vitrification. E the freezing and melting volume fractions of their system Since then, the use of colloids as models has become very were 0.494 and 0.535 respectively, compared to 0.494 and popular, especially since the addition of non-adsorbing arXiv:1106.2566v1 [cond-mat.soft] 13 Jun 2011 0.545 in simulations.[5, 6] Their careful conclusion reads: polymers can be used to induce an inter-particle attrac- ‘Despite ambiguities . in the experimental determina- tion ‘tuneable’ separately in its range and depth [9, 10]. tion of the coexistence region this difference is probably In such colloid-polymer mixtures, gas-liquid coexistence significant.’ can be studied [11, 12], including interfaces and critical- The same degree of caution has not characterized the ity [13, 14], as well as novel modes of arrest [15]. Due literature since. Experimental reports typically do not to their slow intrinsic time scales, colloids are also ideal discuss in any detail the method used for arriving at φ. models for studying phase transition kinetics [16]. On the other hand, theory or simulations almost always While a range of inter-particle interactions are now take experimental reports of φ at face value and proceed available in model colloids, hard spheres remain an im- to use the data on this basis. This is unsatisfactory, par- portant reference system for which very direct compari- ticularly in situations where theory testing demands a son between experiments and theoretical calculations or degree of accuracy and certainty in the experimental φ computer simulations is in principle possible. The be- that is probably unattainable. In this paper, we critically haviour of a single-sized, or monodisperse, system of hard review a plethora of methods for the experimental deter- 2 mination of φ in hard-sphere suspensions, evaluate the (a) (b) degree of accuracy attainable in each case, comment on the potential discrepancies between methods, and give a case study showing how different experiments should be compared taking into account possible difference in φ determination. With the increasing popularity of confocal microscopy, direct counting of particles is becoming a standard method for determining φ (see Sec. V D). This method acs=ac acs ac δ depends on knowing the particle size. Thus, after in- troducing model colloids (Section II), we review particle sizing (Section III). The ‘classic’ method for determin- FIG. 1: Definitions of particle radii. (a) Non-sterically- ing φ is via the crystallization phase behaviour, which stabilized particles, e.g. silica, with a core radius ac. (b) changes with polydispersity[18, 19]. So we review poly- Sterically-stabilized particle, with surface ‘hairs’ (not drawn ¯ dispersity measurements (Section IV) before turning to to scale), where additionally the average hair thickness δ and ¯ consider the determination of φ in detail (Section V). We the core-shall radius aCS = aC + δ are needed for a full char- acterisation. finish with a case study (Section VI) and a Conclusion. of what we say will also apply to silica and other model II. THE PARTICLES systems. We focus on suspensions of nearly-perfect hard spheres. A system that potentially behaves most like perfect III. MEASURING SIZE hard spheres is charge-stabilised colloidal silica in wa- ter. When the charges are sufficiently screened out by The basic parameter characterizing spherical colloidal the addition of salt [20–22], the resulting suspension is particles is their radius, a. Some methods of determining very close to being hard-sphere-like. However, a signifi- φ, most obviously by ‘counting’ from confocal microscopy cant drawback of silica as a model system is that these images, depend directly on measuring a. In this section, particles have a density ρ 2.2 g/cm3, and it has proved we critically review the measurement of particle size. impossible to find solvents≈ that match this density. The sedimentation problem can be alleviated by using smaller particles, but at the expense of increasing polydisper- A. Definition of particle radius sity. Charged-stabilised polystyrene spheres are also po- tentially very hard-sphere-like, and are easier to density At first sight, defining ‘the particle radius’ should be match, though (unlike silica) almost impossible to refrac- simple in a nearly-monodisperse colloid. But this is de- tive index match. ceptive, Fig. 1. For sterically-stabilized particles like A more popular model hard sphere system is PMMA, PMMA, we can usefully define as least four different sterically stabilized by a δ . 10 nm layer of PHSA (poly- radii. First, the hydrodynamic radius, aH , occurs in the 12-hydroxystearic acid) [23]. This layer confers a degree relationship between the drag, f, on a particle moving of softness to the interparticle potential on the scale of at velocity v in a fluid of viscosity η at low Reynolds δ. PMMA particles can be dispersed in solvent mixtures numbers, f = ξv, where the friction coefficient that match both the particles’ density and index of re- fraction [24, 25]. Unfortunately, particle swelling by sol- ξ =6πηaH . (2) vents is endemic [26, 27]. Furthermore, the swelling pro- cess can take several weeks, so the particle size changes Relating aH to φ is a non-trivial problem in hydrody- over time, though heat shock may speed up the process namics. Somewhat more directly related to φ is the core to taking only a few hours [28]. As swelling is poorly radius, aC, which is the radius of the sterically-stabilized characterized, in situ measurement is the only reliable particles minus the stabilizing hairs. Thirdly, if we can means of ensuring that it is complete before the parti- determine the average hair thickness δ¯, then the core- cles are used in experiments. This is particularly impor- shell radius, aCS = aC + δ¯. For ‘hairless’ particles, such tant at high φ, where many properties are steep func- as charge-stabilized silica, aCS = aC. Finally, we may tions of the concentration: an x% increase in the particle assign an effective hard sphere radius, aeff , to the par- radius translates into & 3x% in φ. Thus, e.g., an index- ticles to obtain the best fit to theory or simulations of matching mixture of cis-decalin and tetrachloroethylene hard-sphere behaviour in a certain range of φ, so that causes 20% swelling, which has a drastic effect upon the aeff is inevitably dependent on the chosen property and colloid volume fraction [27]. Finally, batch-to-batch vari- φ range. ations generate further uncertainties. We will not review established sizing methods in any Below, we focus on PMMA particles, although much detail, but will reference existing literature and note cau- 3 tionary points. Then we will introduce a number of newer polydispersity is a significant complication[36].