Sorting Mesoscopic Objects with Periodic Potential Landscapes: Optical Fractionation
Abstract.
Viscously damped objects driven through a periodically modulated potential energy landscape can become kinetically locked in to commensurate directions through the landscape, and thus can be deflected away from the driving direction. We demonstrate that the threshold for an object to become kinetically locked in to an array can depend exceptionally strongly on its size. When implemented with an array of holographic optical tweezers, this provides the basis for a continuous and continuously optimized sorting technique for mesoscopic objects called “optical fractionation”.
Many natural and technologically important processes involve classical transport of small objects through modulated potential energy landscapes. While the generic behavior of modulated transport is well understood in one dimension (1), fundamental questions remain for higher dimensions. Colloidal particles flowing through arrays of optical tweezers (2); (3) provide a uniquely accessible experimental archetype for this class of problems. Experiments on transport through square arrays have revealed a Devil's staircase hierarchy of kinetically lockedin states as a function of orientation (4). Within each lockedin state, particles select commensurate paths through the array independent of the driving direction. The ability to selectively deflect one fraction out of a flowing mixture was predicted (4) to be useful for sorting and purifying mesoscopic samples. Here, we describe a practical implementation of this process, which we term optical fractionation. Examining the kinematics of optical fractionation further reveals that the threshold for kinetic lockin can depend exponentially on size.
Optical fractionation exploits a competition between optical gradient forces exerted by optical traps (2) and an externally applied force, as shown in Fig. 1. A driven particle's trajectory can be deflected enough by an encounter with one trap to pass into the domain of the next, and so on down the line. Such a trajectory is said to be kinetically lockedin to the array. Optical fractionation's selectivity emerges because objects with different sizes, shapes or compositions can experience substantially different potential energy landscapes in the same light field; periodicity emphasizes these differences. Objects that escape from the array flow away in the driving direction, while lockedin objects can be deflected. The two resulting fractions can be collected in separate microfluidic channels downstream.
To demonstrate optical fractionation in practice, we studied the transport of waterborne colloidal particles flowing past a linear array of holographic optical tweezers (3). The colloidal suspension was confined to a glass channel formed by bonding the edges of two glass cover slips. Capillary forces at the channel's inlet were used to create a flow of about along the midplane. This flow carried a mixture of monodisperse silica spheres of radius (Duke Scientific Lot No. 24169) and (Duke Scientific Lot No. 19057), which can be distinguished visually and tracked to within 30 nm in the plane at 1/60 s intervals using digital video microscopy (5).
Colloidal silica spheres are roughly twice as dense as water and settle into a monolayer just above the channel's lower wall, with the smaller spheres floating about 1 higher. Given the Poisseuille flow profile in the channel, the larger spheres travel slower, with a mean speed of , compared with the smaller spheres' . The associated driving force, , is characterized by a sizedependent drag coefficient, , modified by proximity to surfaces (6).
Twelve discrete optical tweezers were arranged in a line with centertocenter spacing oriented at with respect to the flow. Each trap was powered by of laser light at 532 nm, which slightly exceeded the empirically determined lockin threshold for the larger spheres, given and . Each trap can capture either type of sphere in the absence of flow. The trapping plane was adjusted to minimize outofplane motions, so that the system is effectively twodimensional for both populations. Although this is useful for illustrative purposes, optical fractionation also works in thick samples with threedimensional trap arrays.
The trajectories in Fig. 2(a) and (b) demonstrate that the larger spheres are indeed systematically deflected by the array of traps, while the smaller spheres are not. Consequently, the array creates a shadow in the distribution of large spheres into which the small spheres can flow. Because the purification of small spheres and concentration of large results from lateral deflection, this optical fractionation process can proceed continuously, in contrast to most competing techniques (7).
Figures 2(c) and (d) show statistics compiled from tens of thousands of trajectories. Here, we plot the two populations' timeaveraged areal densities normalized by their means. The separation's quality is assessed in Fig. 3 with , which reaches unity in regions containing only large spheres, and minus one in regions with only small spheres. A transverse section taken along line in Fig. 3(a) and plotted as squares in Fig. 3(b) reveals a thoroughly mixed sample with approaching the traps. A similar section along line downstream of the array demonstrates roughly 40 percent purification of both large and small spheres. Much of the background can be attributed to large particles escaping from the weakest traps in our array (4). In denser suspensions, this escape rate is increased by collisions. Both processes can be mitigated by projecting multiple lines of traps, with perfect deflection of the lockedin fraction having been demonstrated (4) at throughputs exceeding 10,000 particles/sec. Larger arrays, thicker samples, and faster flows would facilitate much higher throughputs.
Optical fractionation's ability to distinguish objects arises as a general and previously unappreciated feature of transport through periodically structured environments. Analyzing such transport not only provides insights into optimizing practical sorting, but also sheds new light on a range of analogous processes.
The potential energy landscape presented by an optical trap array is a convolution of the traps' intensity profile with an object's optical form factor : . The total applied force then is , where is the driving force. An overdamped particles' trajectory, , has one component, , along the row of traps and another, , perpendicular. Although this problem is well understood in one dimension (1), few analytic results are available for the inclined line, and fewer still incorporate thermal or quenched randomness. Consequently, we focus on the kinematic limit in which the driving and trapping forces dominate and trajectories may be treated deterministically. We then estimate the threshold for an object to escape from an array of optical traps, and thereby establish the selectivity of optical fractionation.
Any trajectory entrained by the traps, such as the example in Fig. 1(b), is characterized by turning points where . Conversely, any trajectory without such turning points must be unbounded. This establishes as the maximum possible lockedin deflection angle
(1) 
At a given , objects are either deflected or not with a selectivity set by the dependence of on material properties. To estimate , we model the array as a periodically modulated line of light with intensity :
(2) 
where , the dimensionless transverse distribution ) is symmetrically peaked around , and the coefficients account for the tweezers' detailed structure with .
Convolving first along the direction by applying the Fourier convolution theorem to each term in the sum, and then noting that yields
(3) 
where is the potential wells' depth and is the form factor's Fourier transform along the direction. The array's periodicity thus selects a discrete set of wavenumbers from the continuous whose dependence on endows optical fractionation with exceptional size selectivity. This is most clearly demonstrated if can be factored into inline and transverse components, . In this case,
(4) 
with and . Equivalent results can be obtained when is not separable.
We turn our attention first to the transverse contribution. If a particle is comparable in size to the optical tweezers' width, , then depends no more strongly on size than . For example, if and are Gaussians of widths and , respectively, then . Similarly, the potential depth and driving force generally scale as simple powers of . Comparable algebraic sensitivity to size and material properties is offered by other techniques such as gel electrophoresis and field flow fractionation (7), and would be obtained with an unmodulated line of light ().
The inline contribution is more interesting. Because a particle's form factor vanishes outside the interval , its Fourier transform depends very strongly on wavenumber. For example, a uniform dielectric cube aligned with the array has a separable form factor, for whose Fourier transform, , is bounded by the leadingorder cumulant expansion, . The equivalent result for a sphere (8) with for is , and . All of these bounding approximations surpass exponential selectivity for , with the actual form factors depending even more strongly on . Any smooth, bounded, positivedefinite on would behave similarly. Applying this insight to Eq. (4) establishes the lockin transition's exponential size sensitivity for
Comparably strong dependence on control parameters is observed in analogous transitions between subharmonic steps in driven charge density waves (9) and between kinetically lockedin states in twodimensional optical trap arrays (4). Similar results also can be obtained for arrays of potential barriers, suggesting that arrays of optical tweezers also should be effective for sorting absorbing, reflecting and lowdielectric particles that are repelled by laser light. This analysis also carries over to filtration by arrays of micromachined posts (10), which therefore should be able to resolve objects substantially smaller than the interpost separation.
Both and the coefficients fall off rapidly with index . Consequently, the sum in Eq. (4) is dominated by the first term, . This improves the approximations used in deriving Eq. (4) and suggests that the result may be considered an estimate for rather than simply a bound.
To demonstrate this, we apply this analysis to our present experimental data, modeling the individual optical traps as Gaussian potential wells
(5) 
with (11). In this model,
(6) 
where . The weakest trap's occupancy, , is inversely proportional to the particles' minimum speed, , as they pass through. Consequently, we can estimate the relative trap strength from the data in Fig. 2 as . Similarly, the separation between the depleted region ahead of the traps and the position of maximum occupancy is . From the histograms in Fig. 2(c) and (d), we obtain and , and and for the large and small spheres, respectively (11). These results suggest for the large spheres and for the small, which is consistent with the observation that only the large spheres are systematically deflected at .
The threshold, , depends only linearly on and . Thus, imperfections in practical trap arrays and fluctuations in the driving force can be largely compensated for by the substantially stronger dependence on particle size. Indeed, Figs. 2 and 3 demonstrate robust size separation despite more than 20 percent variation in flow velocity over the course of the experiment.
Equation (6) also offers insights into applying optical fractionation to nanometerscale objects. Stokes drag scales linearly with , and the optical trapping potential for Rayleigh particles scales as , so that . Sorting proteins or nanoclusters, therefore, will require enhancing by four orders of magnitude. This might be accomplished by increasing the light's intensity and reducing its wavelength (12). Even then, only algebraic size sensitivity should be expected for objects with because in this limit. Exploiting resonances might overcome this limitation.
We have focused on effects due to deflection transverse to the optical axis. Multidimensional separations could take advantage of Bessel beams' ability to exert controlled axial forces (13) to distribute objects both transverse to and along the optical axis.
In summary, we have demonstrated optical fractionation for a model system of bidisperse colloidal spheres. This approach lends itself to continuous, rather than batchmode fractionation, with continuous tuning and dynamic optimization over the entire accessible size range, i.e. nanometers to micrometers. The abrupt transition from free flow to kinetically lockedin transport should offer exponential size selectivity for objects larger than roughly 100 nm. Separation on the basis of other characteristics also can be optimized, although exponential sensitivity should not be expected in general. Our analysis focuses on the kinematic limit, , which is both tractable and appropriate for weaklytrapped micrometerscale colloid. Stronger trapping would require a more detailed treatment of thermally assisted hopping (1). A substantially more sophisticated analysis also would be required for higherdimensional arrays. Similarly abrupt transitions should occur with variation of driving or trapping strength for vortex creep through patterned typeII superconductors (14), electron transport through twodimensional electron gases (15), and electromigration on crystal surfaces, with potentially useful applications resulting in each case.
Since submitting this article, we have become aware of a subsequently submitted independent experimental study of colloidal sorting in an optical lattice (16). We are grateful to Paul Chaikin and Matthew Pelton for enlightening discussions. This research was supported primarily by the National Science Foundation through Grants Number DMR0304906 and DBI0233971 and in part by the MRSEC and REU programs of the NSF through Grant Number DMR0213745.
References

(1)
H. Risken, The FokkerPlanck Equation (SpringerVerlag, Berlin, 1989), 2nd ed.

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

(3)
E. R. Dufresne and D. G. Grier, Rev. Sci. Instr. 69, 1974 (1998); E. R. Dufresne, G. C. Spalding, M. T. Dearing, S. A. Sheets, and D. G. Grier, Rev. Sci. Instr. 72, 1810 (2001a); J. E. Curtis, B. A. Koss, and D. G. Grier, Opt. Comm. 207, 169 (2002).

(4)
P. T. Korda, M. B. Taylor, and D. G. Grier, Phys. Rev. Lett. 89, 128301 (2002).

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

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

(7)
I. D. Wilson, ed., Encyclopedia of Separation Science (Academic Press, San Diego, 2000).

(8)
C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley Interscience, New York, 1983).

(9)
S. E. Brown, G. Mozurkewich, and G. Grüner, Phys. Rev. Lett. 52, 2277 (1984).

(10)
T. A. J. Duke and R. H. Austin, Phys. Rev. Lett. 80, 1552 (1997).

(11)
This form applies in both Rayleigh () and rayoptics () limits (12), and is consistent with our measurements, assuming .

(12)
P. A. M. Neto and H. M. Nussenzveig, Europhys. Lett. 50, 702 (2000).

(13)
J. Arlt, V. GarcesChavez, W. Sibbett, and K. Dholakia, Opt. Comm. 197, 239 (2001).

(14)
C. Reichhardt and F. Nori, Phys. Rev. Lett. 82, 414 (1999).

(15)
J. Wiersig and K.H. Ahn, Phys. Rev. Lett. 87, 026803 (2001).

(16)
M. P. MacDonald, G. C. Spalding and K. Dholakia, Nature 426, 421 (2003). See K. Ladavac, M. Pelton and D. G. Grier, condmat/0404242 for a critical discussion of the sorting mechanism proposed in this paper.