Anomalous vibrational dispersion in holographically trapped colloidal arrays
Abstract.
Colloidal spheres localized in an array of harmonic wells form a thermally excited, viscously damped dynamical system capable of supporting propagating elastic waves. Experimentally realized with micrometerscale polystyrene spheres localized in a line of holographic optical traps, the hydrodynamically coupled array's behavior is quantitatively explained by a model based on the Oseen superposition approximation. The spheres' purely dissipative coupling is predicted to mediate a crossover to a regime of underdamped propagating elastic waves with uniformly negative group velocities that has yet to be verified experimentally.
The manybody dynamics of driven dissipative systems often are governed by hydrodynamic coupling. Even in the simplest systems, this can lead to interesting and unexpected behavior. For example, Meiners and Quake (1) demonstrated that the thermally driven fluctuations of two noninteracting colloidal spheres individually localized in the harmonic potential wells of two optical tweezers are strongly anticorrelated at short times and small separations. Cooperative motion in this case is a consequence of the hydrodynamic forces the spheres exert on each other as they relax into their respective potential wells. Its appearance suggests that an extended system of independently trapped hydrodynamically coupled spheres might behave as an elastic medium, and that the spheres' predilection for anticorrelated motion might give rise to unusual timedependent properties.
This Letter describes an experimental and theoretical study of thermally driven vibrations in a onedimensional hydrodynamically coupled colloidal crystal created with the holographic optical trapping technique (2); (3); (4). Each sphere is localized in an individual optical tweezer (5), and the spheres interact only through their influence on the surrounding viscous medium. We show experimentally that the normal modes of such a dissipatively coupled solid are consistent with a LangevinOseen model. This theory further predicts that elastic waves can propagate in this system despite the absence of conservative interactions. These waves exhibit anomalous dispersion and negative group velocities as a generic feature.
Related phenomena have been studied in dusty plasmas whose highly charged constituent particles can be harmonically constrained in one or two dimensions by external electric fields (6); (7); (8); (9), and in chargestabilized colloidal crystals (10); (11). The present system differs from these and most other solids in that cooperativity arises from the spatial dependence of a purely dissipative interaction. Negative group velocities also distinguish the optical properties of cold atomic gases and certain semiconducting systems, where they result from resonances.
Our experimental setup is sketched in Fig. 1(a). We create arrays of uniformly strong optical tweezers by imprinting computergenerated holograms (4) onto the wavefronts of a laser beam (Coherent Verdi) with a liquid crystal spatial light modulator (SLM) (Hamamatsu X7550 PALSLM). The modified beam is focused into an array of traps with a 100, NA 1.4, SPlan Apo oil immersion objective lens mounted in a Zeiss S100TV Axiovert inverted optical microscope. Placing a frontsurface mirror on the microscope's stage allows us to photograph the focused trap array, as in Fig. 1(b). Photometry on such images confirms that the traps' intensities vary by less than 5 percent from the mean, and thus that their trapping characteristics should be comparably uniform (4).
The traps are projected into a sample of colloidal polystyrene spheres in diameter (Bangs Labs Lot number 6064) dispersed in a 75 percent (w/w) mixture of glycerol (Fisher Scientific Lot number 04377) in deionized water. Allowing the dispersion to equilibrate in air reduces the screening length to less than 300 (12); (13) and effectively eliminates electrostatic interactions among the spheres. Adding 1 mM NaCl reduces the screening length to 12 without affecting the results that follow. The sample is sealed within a slit pore thick formed by bonding the edges of a #1 cover glass to a microscope slide, and the traps are focused to the midplane, , to minimize hydrodynamic coupling of the spheres to the walls.
We regulate the sample's temperature, and thus its viscosity by heating from above with a transparent indiumtinoxide heater whose temperature is fixed to by a Lake Shore 330 temperature controller. Thermal contact with the microscope objective at establishes a vertical temperature gradient, with the temperature at the trapping plane estimated to be . Under these conditions, the solution should have a viscosity of , which corresponds to a freesphere viscous drag coefficient of . Evaporation of water during sample preparation, however, can increase this substantially, and we instead measure in situ.
Figure 1(c) shows a typical video frame of ten spheres trapped in the array. Experiments were performed with centertocenter separations ranging from to , which are large enough to avoid artifacts due to overlapping images (14). Standard methods of digital video analysis (15) allow us to track the spheres' centroids in the plane with resolution at 1/30 intervals. With an estimated 3 mW of light per trap, the particles' thermally driven root meansquared (RMS) displacements are around . Tracking individual spheres over 5 minutes allows us to characterize both the traps' potential wells and the particles' viscous drag coefficients (4). The traps are accurately described as symmetric harmonic wells in the plane with a stiffness of . The uncertainty here is comparable to the singletrap measurement error (4) of , with additional contributions coming from a combination of the measured differences in the traps' intensities and the 2.2 percent variation in these spheres' diameters inferred from measurements of . The extracted value, , is consistent with an 81.6 percent (w/w) glycerol concentration. The associated viscous relaxation time, , exceeds the video sampling interval of 33 sufficiently that the array's dynamics can be captured by conventional digital video analysis.
The particles' measured trajectories were decomposed into the normal modes anticipated for hydrodynamically coupled spheres. We compute these using Faxén's law (16) in the approximation that a particle at location is advected by the flow field generated by its neighbors. In the absence of other interactions, a system of identical spheres of radius at positions responds to a set of forces by moving with velocities dictated by the body Oseen tensor
(1) 
Assuming that the spheres are well localized in traps at , with , we may approximate by the pairwise superposition of stokeslets (17),
(2) 
where , , and is the identity tensor. This approximation has proved effective in previous experimental studies of colloidal hydrodynamics (1); (18); (19); (20). If we further assume that the spheres' fluctuations about their equilibrium positions are much smaller than their mean separation, , we may replace with in Eq. (2) to obtain a timeindependent approximation to . The Langevin equation for the th sphere is then
(3) 
where is the sphere's mass and describes the thermal forces acting on the sphere. For a system in equilibrium at temperature , and
(4) 
where is Boltzmann's constant. Retaining the inertial term in Eq. (3) enables us investigate conditions under which the optically trapped array admits propagating modes.
The Oseen tensor can be decomposed into longitudinal () and transverse () components relative to the array's axis, . The associated normal modes, and , are eigenvectors of the matrices
(5)  
(6) 
Measured trajectories are combined accordingly.
The eigenvalue associated with the normal mode in a particular polarization describes that mode's time evolution through the body Langevin equation
(7) 
where the thermal forcing terms are Gaussian random variables with zero mean and covariance
(8) 
Accordingly, solutions to Eq. (7) in the limit of strong damping satisfy
(9) 
where the variance, , and the relaxation time, , completely characterize the motion.
The success of this description can be gauged from comparisons with data in Fig. 2. Normal mode fluctuations were recorded for 20 minutes at each of seven different lattice constants ranging from to . Both the normal modes and their associated eigenvalues are determined solely by the geometric ratio , and the data were analyzed accordingly. As expected for a thermally equilibrated system obeying equipartition, the normal modes' normalized covariances, , plotted in Fig. 2(a) collapse to unity to within 10 percent, independent of lattice constant. The common mode (mode 0) departs from unity due to its extreme sensitivity to slow instrumental drifts. Its offset from the other modes in all seven runs is consistent with a random drift of over 10 minutes. The other normal modes are more robust against such artifacts, and their variances are consistent with the measured differences in the traps' intensities. The inverse autocorrelation times plotted in Fig. 2(b) also agree with calculated eigenvalues, thereby confirming the correct decomposition of the particles' motions into normal modes. Successful collapse of the data over a range of interparticle separations also helps to rule out any effect of electrostatic interactions.
Having verified the model's accuracy in this range of parameters, we can now comment more generally on the nature of vibrational waves' propagation in this class of systems. From Eq. (7), the frequency associated with the th normal mode is
(10) 
where and . For the present system, . We recast Eq. (10) as a dispersion relation for elastic waves by associating a wavenumber with each mode. For an infinite array, translation invariance implies eigenmodes of the form
(11) 
with wavenumbers , amplitudes and polarization vectors or . The associated eigenvalues are
(12)  
(13) 
Substituting these into Eq. (10) yields the dispersion relations, and , for longitudinal and transverse elastic waves, respectively. Wavenumbers may be associated with the eigenmodes of the finite system by fitting sinusoids to the eigenvectors of .
Figure 3 shows representative data from Fig. 2, obtained at , replotted as the imaginary branch of a dispersion relation. The results agree quantitatively with Eqs. (10), (12) and (13), which are shown as dashed curves. The absence of propagating modes in such a highly damped system is not surprising. Diverging eigenvalues at longer wavelengths, however, guarantee a crossover to realvalued eigenfrequencies and thus to propagating modes for any values of , and , even when the singlesphere Reynolds number is small. The crossover to real dispersion can be moved into the experimentally accessible realm by trapping larger spheres more rigidly in a lowerviscosity medium, as shown by the solid curves in Fig. 3. Under all conditions, decreases monotonically with so that the propagating excitations all have negative group velocities. Both the crossover to underdamped wave propagation and also these modes' anomalous dispersion arise not from direct interparticle interactions as in other systems, but rather from the separation dependence of the viscous damping.
In practice, hydrodynamic coupling to bounding surfaces tends to suppress longwavelength propagating modes. Accounting for these in the stokeslet approximation (17); (21) (dotted curves in Fig. 3) yields no measurable change in the predicted dispersion for our experimental conditions, and no qualitative change for the predicted crossover to propagating longitudinal modes.
A crossover from overdamped to underdamped dynamics at large wavelengths has been predicted (10) and recently observed (11) for the transverse lattice vibrations of colloidal crystals. This crossover, however, is governed by direct interparticle interactions rather than by hydrodynamic coupling. Colloidal crystals consequently exhibit normal dispersion for transverse waves in the propagating regime, and their longitudinal modes are always overdamped. The hydrodynamically coupled lattice, by contrast, supports propagating modes with anomalous dispersion in both polarizations.
Underdamped elastic waves propagate in dusty plasmas, with no crossover to overdamped dynamics (8); (9). The plasma medium has low enough viscosity that hydrodynamic coupling between neighboring particles is negligible. Here again, confining the spheres to one or two dimensions with externally imposed potentials induces anomalous dispersion for waves polarized along the confined directions (8); (9). In this case, however, negative group velocities result from the interplay of confining potentials and the particles' electrostatic interactions. Spheres in the present system, by contrast, interact only through their influence on the viscous medium, with propagating modes vanishing at low viscosity. Optically organized arrays of noninteracting spheres in viscous fluids therefore constitute a distinct class of elastic media whose ability to transmit elastic waves is a consequence of dissipation instead of conservative interactions.
We thank Paul Chaikin and John Goree for helpful discussions. This work was supported by the National Science Foundation under Grant Number DMR0451589.
References

(1)
J.C. Meiners and S. R. Quake, Phys. Rev. Lett. 82, 2211 (1999).

(2)
E. R. Dufresne and D. G. Grier, Rev. Sci. Instr. 69, 1974 (1998).

(3)
J. E. Curtis, B. A. Koss, and D. G. Grier, Opt. Comm. 207, 169 (2002).

(4)
M. Polin, K. Ladavac, S.H. Lee, Y. Roichman, and D. G. Grier, Opt. Express 13, 5831 (2005).

(5)
A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, Opt. Lett. 11, 288 (1986).

(6)
S. V. Vladimirov, P. V. Shevchenko, and N. F. Cramer, Phys. Rev. E 56, R74 (1997).

(7)
G. Piacente, F. M. Peeters, and J. J. Betouras, Phys. Rev. E 70, 036406 (2004).

(8)
T. Misawa, N. Ohno, K. Asano, M. Sawai, S. Takamura, and P. K. Kaw, Phys. Rev. Lett. 86, 1219 (2001).

(9)
B. Liu, K. Avinash, and J. Goree, Phys. Rev. Lett. 91, 255003 (2003).

(10)
A. J. Hurd, N. A. Clark, R. C. Mockler, and W. J. O'Sullivan, Phys. Rev. A 26, 2869 (1982).

(11)
B. V. R. Tata, P. S. Mohanty, M. C. Valsakumar, and J. Yamanaka, Phys. Rev. Lett. 93, 268303 (2004).

(12)
C. S. Miner and N. N. Dalton, Glycerol (Reinhold, New York, 1953).

(13)
J. N. Butler Carbon Dioxide Equilibria and Their Application (AddisonWesley, London, 1982).

(14)
J. Baumgartl and C. Bechinger, Europhys. Lett. 71, 487 (2005).

(15)
J. C. Crocker and D. G. Grier, J. Colloid Interface Sci. 179, 298 (1996).

(16)
J. Happel and H. Brenner, Low Reynolds Number Hydrodynamics (Kluwer, Dordrecht, 1991).

(17)
C. Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow (Cambridge University Press, New York, 1992).

(18)
J. Crocker, J. Chem. Phys. 106, 2837 (1997).

(19)
E. R. Dufresne, T. M. Squires, M. P. Brenner, and D. G. Grier, Phys. Rev. Lett. 85, 3317 (2000).

(20)
E. R. Dufresne, D. Altman, and D. G. Grier, Europhys. Lett. 53, 264 (2001).

(21)
A. J. Hurd, N. A. Clark, R. C. Mockler, and W. J. O'Sullivan, J. Fluid Mech. 153, 401 (1985).