Weak long-ranged Casimir attraction in colloidal crystals
We investigate the influence of geometric confinement on the free energy of an idealized model for charge-stabilized colloidal suspensions. The mean-field Poisson-Boltzmann formulation for this system predicts pure repulsion among macroionic colloidal spheres. Fluctuations in the simple ions' distribution provide a mechanism for the macroions to attract each other at large separations. Although this Casimir interaction is long-ranged, it is too weak to influence colloidal crystals' properties.
Experimental evidence collected over 20 years (1) suggests that similarly-charged colloidal spheres dispersed in water need not simply repel each other. Under some circumstances they instead experience an unexpected long-ranged attraction. For example, like-charge attractions are implicated in the cohesion of metastable superheated colloidal crystals (2); (3) even though isolated pairs of the constituent spheres are observed to repel each other (4); (5). Comparable attractions have been measured for pairs of spheres confined by two (5); (6) charged planar walls. Recent calculations (7); (8) reveal that such confinement-induced attractions cannot be accounted for by local density theory nor by electrohydrodynamic coupling (9); (10). Such anomalous effects in charge-stabilized colloid therefore challenge our general understanding of interactions and dynamics in macroionic systems.
This Letter addresses fluctuations' contribution to the free energy of highly charged colloidal spheres surrounded by a neutralizing cloud of small singly-charged counterions. Highly symmetric monopolar fluctuations in the counterion distribution increase the system's free energy. We demonstrate that their suppression by boundary conditions at the spheres' surfaces introduces a long-range attraction into the crystal's free energy analogous to the Casimir force in quantum electrodynamics, but that it is too weak to account for anomalous behavior in charge-stabilized suspensions.
Our treatment is based on the Wigner-Seitz cell model introduced by Wennerström, Jönsson and Linse (11) which has been studied extensively (12) both theoretically and through Monte Carlo simulation. It consists of a single spherical macroion of radius carrying a uniformly distributed surface charge and surrounded by a thermal cloud of point-like counterions at temperature , each carrying a single charge . The macroion and counterions are confined by a concentric conducting spherical shell of radius . This outer shell plays a role analogous to the Wigner-Seitz cell boundary in a colloidal crystal. More generally, it models the crowding or geometric confinement characteristic of colloidal crystals (13).
Previous investigations of this and related models (12); (14) have found short-ranged correlation-driven attractions between the bounding surfaces under some conditions, particularly when the counterions are polyvalent. They have not found evidence for long-ranged attractions in monovalent electrolytes (13).
Our method for evaluating the counterions' partition function allows us to investigate much higher macroion charges than have been considered before. The outer boundary's suppression of counterion fluctuations induces a long-ranged Casimir-like attraction (14) between the macroion and its neighbors across the Wigner-Seitz cell boundary. Although this cell model is far too simple to describe the behavior of real charge-stabilized suspensions, it highlights a previously unexplored mechanism for long-ranged confinement-induced like-charge colloidal attractions.
We adopt the path integral formalism reviewed in (14) and write the counterions' canonical partition function as a functional integral over all possible counterion distributions, :
where is the thermal energy scale at temperature , the prime indicates that the number of charges is conserved (), and (15) The potential energy functional describes the counterions' interaction with the local electric potential . The system's Helmholtz free energy is then .
One ionic distribution, , minimizes , and thus has the greatest statistical weight in . We factor the partition function into the saddle point contribution where , and a term accounting for fluctuations, , away from . Expressing as a series expansion in yields as the lowest-order non-vanishing term
where is the second-order change in due to .
Terminating this expansion at Gaussian order is justified if corrections at higher order in contribute negligibly to . This condition is met if is large (16) and itself changes negligibly over the mean radial counterion separation:
Selecting sets the potential's reference point conveniently without loss of generality. Eq. (4) accounts for the solvent's influence in the so-called primitive model through its dielectric constant .
Following conventional practice, we introduce the Bjerrum length and an effective screening length . Taking the system's radial symmetry into account leads to
We solve Eq. (5) subject to two boundary conditions: Gauss' theorem at the macroion's surface gives
and electroneutrality requires
Figure 1 shows as a function of , calculated numerically from the solution of the full PB equation. Since decreases monotonically with , the mean-field theory predicts pure repulsion.
An equivalent result obtained by solving the linearized PB equation for the same system does have a local minimum for slightly bigger than . In the context of colloidal interactions, the linearized mean-field description yields the electrostatic component of the pair potential due to Derjaguin, Verwey, Landau and Overbeek (DLVO) (19). The local minimum in the free energy suggests an effective electrostatic attraction between like-charged colloidal spheres within the DLVO theory. Comparison with the full calculation, however, shows this to be an artifact of linearization.
In addition to confirming the absence of attractions in the mean-field description, our numerical results also satisfy the condition in Eq. (3). Thus, we are justified in using Eq. (2) to calculate . The second-order change in due to is
So, from Eq. (2),
We evaluate the functional integral by partitioning the system into concentric shells,
where is the volume of cell , the number density of counterions therein, and the equivalent number density for the mean-field case. This highly symmetric partition of the system is appropriate for radial or monopolar density fluctuations. A more complete description including multipole fluctuations is not necessary for our purposes because recent calculations reveal that these yield only short-ranged attractions (20). Cross terms cancel at the Gaussian level of approximation, so we can consider monopole and multipole contributions separately.
We define , and divide the system so that every cell has the same number of counterions in the mean-field distribution, i.e. . Rescaled in this way,
Then, numbering the cells from the center outwards so that if , we obtain
The prime indicates that the integral is constrained to maintain electroneutrality: . Because of this condition, we can subtract a constant from each element of without changing the value of . Let us subtract from each element of to obtain matrix . All of the elements in the last row and column of vanish except which is . Separating mode from the other modes in this way allows us to remove the explicit constraint:
where is now a vector of dimension , and is a matrix of dimension obtained by adding to each element in the first rows and columns of matrix .
may be expressed as the sum of two matrices, and , whose components are and
where is the greater of and . , so that where is the identity matrix.
Evaluating is greatly facilitated if the components of are much smaller than 1. differs little from the identity matrix. The components of , on the other hand, are bounded above by . We previously assumed in deriving Eq. (2). But for the micron-sized spheres in experimental observations, so that we may reasonably assume . Even if this were not the case, we would be justified in formally taking the limit at this point because the final result cannot depend on . Consequently,
In this approximation, the fluctuation contribution to the free energy is
Rewriting the sums over shell indices as integrals over radii, we obtain
Most of the boundaries, , between cells in the mean-field distribution are clustered near . Consequently,
where is the fraction of counterions within radius , in the mean field approximation. Integrating by parts then yields
Unlike , decreases with decreasing because the outer boundary condition suppresses fluctuations as approaches . The resulting attraction therefore is reminiscent of the Casimir attractions previously identified in confined electrolytes as well as other systems (14). It is interesting to note that this mechanism does not yield an attractive contribution in all geometries; the second term in the exponent of Eq. (13) vanishes for unbounded systems such as parallel plates and concentric cylinders. This is consistent with the absence of long-ranged like-charge attractions in measurements (21), theoretical treatments (12); (14), and simulations (12); (14) of unbounded systems. Such attractions, therefore, are peculiar to closed systems, such as the Wigner-Seitz cells of colloidal crystals. If there is an at which the attractive force
has larger absolute value than the mean-field repulsive force
then would have a minimum at that . Such a minimum would correspond to a fluctuation-mediated bound state for a colloidal crystal of nearest neighbor spacing . However, numerical results shown in Fig. 2 reveal over the entire range of conditions studied experimentally. On this basis, we conclude that monopolar fluctuations are not responsible for the strong and long-ranged attractions reported in measurements on charge-stabilized colloid (2); (3); (5); (6).
To understand why the ratio is small we estimate it analytically by assuming that most of the counterions are clustered close to the macroion surface and that only a small fraction of counterions are affected by a change in the radius of the confining shell, and that this fraction is uniformly distributed in the volume. We then evaluate and to . Truncating to is justified by numerical investigation of the mean-field solution which indicates in the region of interest. In this approximation
where , is independent of and has a maximium value at . The location and magnitude of the maximum value agrees with the full numerical solution to within factors of 1.2 and 4, respectively. This value of would correspond to half the nearest-neighbor separation in a fluctuation-stabilized colloidal crystal; it also is consistent with the lattice constants observed in superheated metastable colloidal crystallites (2); (3). However, water has a dielectric constant which corresponds to a Bjerrum length nm at room temperature. Consequently, for all reasonable colloidal radii. Changing solvents would not affect this conclusion because of the limited range of accessible values of .
We also investigated the possible influence of multivalent counterions carrying charge . Ignoring the relatively weak dependence of yields . Despite the relative strength of the attraction being larger for higher valency counterions, the effect is still too weak to induce measurable attractions for physically plausible values of .
We have demonstrated that suppression of monopolar ionic fluctuations by surfaces induces a long ranged attraction remniscent of Casimir attractions (14). This interaction is distinct from and complementary to attractions driven by multipolar fluctuations which have been studied elsewhere (20). Neither mechanism, however, accounts for the strong and long ranged attractions observed experimentally between highly charged colloidal spheres.
The long-ranged like-charge attractions observed in confined colloid are not consistent with mean field theories for electrolyte structure. Possible explanations must incorporate mechanisms such as fluctuations and high-order correlations not captured by mean field theory. While multipole fluctuations in the distribution of simple ions induce strong attractions (20), they are short ranged. The present study demonstrates that radially symmetric fluctuations can induce long-ranged attractions, but that they are too weak to influence colloidal behavior. Consequently, the explanation must lie in another mechanism not yet considered and thus remains an important outstanding challenge.
Acknowledgements.We are grateful to Tom Witten, Stuart Rice, Phil Pincus, Mehran Kardar, Ramin Golestanian, An-Chang Shi and Adrian Parsegian for enlightening conversations. This work was supported by the MRSEC Program of the National Science Foundation under Award Number DMR-9808595.
Ise N., Ito K., Matsuoka H. & Yoshida H., Ordering and Phase Transitions in Charged Colloids, A. K. Arora & B. V. R. Tata eds. (VCH Publishers, New York, 1996) .
Larsen A. E. & Grier D. G., Phys. Rev. Lett. 76, 3862 (1996).
Larsen A. E. & Grier D. G., Nature 385, 230 (1997).
Crocker J. C. & Grier D. G., Phys. Rev. Lett. 73, 352 (1994); Sugimoto T. et al., Langmuir 13, 5528 (1997).
Crocker J. C. & Grier D. G., Phys. Rev. Lett. 77, 1897 (1996).
Kepler G. M. & Fraden S., Phys. Rev. Lett. 73, 356 (1994); Carbajal-Tinoco M. D., Castro-Román F. & Arauz-Lara J. L., Phys. Rev. E 53, 3745 (1996).
Neu J. C., Phys. Rev. Lett. 82, 1072 (1999); Sader J. E. & Chan D. Y. C., J. Colloid Interface Sci. 213, 268 (1999).
Trizac E., Phys. Rev. E 62, R1465 (2000).
Squires T. M. & Brenner M. P., Phys. Rev. Lett. 85, 4976 (2000).
Behrens S. H. & Grier D. G., unpublished.
Wennerström H., Jönsson B. & Linse P., J. Chem. Phys. 76, 4665 (1982).
Attard P., Adv. Chem. Phys. XCII, I. Prigogine & S. A. Rice eds. (Wiley, New York, 1996) .
Grier D. G., J. Phys.: Cond. Matt. 12, A85 (2000).
Kardar M. & Golestanian R., Rev. Mod. Phys 71, 1233 (1999).
Zhou T., Phys. Lett. A 267, 389 (2000).
An-Chang Shi, private communication.
van Roij R. & Hansen J.-P., Phys. Rev. Lett. 79, 3082 (1997).
Netz R. R. & Orland H., Europhys. Lett. 6, 726 (1999).
Derjaguin B. V. & Landau L., Acta Phys. (URSS) 14, 633 (1941); Verwey E. J. & Overbeek J. T. G., Theory of the Stability of Lyophobic Colloids, (Elsevier, Amsterdam, 1948).
Levin Y., Physica A 265, 432 (1999).
Israelachvili J., Intermolecular and Surface Forces, (Academic, London, 1992).