Influence of nonconservative optical forces on the dynamics of optically trapped colloidal spheres: The fountain of probability
Abstract.
We demonstrate both experimentally and theoretically that a colloidal sphere trapped in a static optical tweezer does not come to equilibrium, but rather reaches a steady state in which its probability flux traces out a toroidal vortex. This nonequilibrium behavior can be ascribed to a subtle bias of thermal fluctuations by nonconservative optical forces. The circulating sphere therefore acts as a Brownian motor. We briefly discuss ramifications of this effect for studies in which optical tweezers have been treated as potential energy wells.
Most discussions of the dynamics of optically trapped particles assume at least implicitly that the forces exerted by an optical tweezer (1) are pathindependent and therefore conserve mechanical energy. Optical forces due to gradients in the intensity are manifestly conservative in this sense (2). Radiation pressure, by contrast, is not (2); (3). The experimental studies described in this Letter demonstrate that the nonconservative component of the optical force has measurable consequences for the dynamics of optically trapped colloidal spheres. In particular, the probability density for a sphere trapped in a static optical tweezer exhibits steadystate toroidal currents, a phenomenon we call the fountain of probability. We use a FokkerPlanck formalism to explain how nonconservative forces bias random thermal fluctuations to induce circulating probability currents.
Figure 1 schematically represents the deceptively simple system. A single colloidal sphere is drawn to the focus of a converging laser beam by forces arising from gradients in the beam's intensity (1); (2). These intensitygradient forces establish a threedimensional potential energy well, , determined by the local intensity, . A particle trapped in this well also experiences radiation pressure that drives it downstream with a force proportional to . In the absence of thermal fluctuations, a trapped particle would come to rest at a stable mechanical equilibrium downstream of the focus.
Treating the displaced equilibrium point as the origin of an effective potential energy well is tempting but misleading. To appreciate the problem, consider a thermally driven trajectory such as the example shown schematically in Fig. 1(a). Were the system in thermodynamic equilibrium, forward and reverse trajectories around this loop would have equal probability. Because the light is more intense near the optical axis, however, radiation pressure biases the random walk in the forward direction. This departure from detailed balance biases the particle's fluctuations to create irreversible circulation (4).
We demonstrate this effect by observing the motions of colloidal silica spheres 2.2 in diameter (Bangs Laboratories, SS04N/7651) dispersed in water and individually trapped in discrete optical tweezers. Each trap is formed from a few milliwatts of laser light at a vacuum wavelength of 532 nm (Coherent Verdi) projected with the holographic optical trapping technique (5). The beams are brought to a diffractionlimited focus by an objective lens (Nikon Plan Apo , oil immersion) mounted in an inverted optical microscope (Nikon TE2000U). We track the particle in three dimensions with nanometer precision using video holographic microscopy (6); (7). Inline holograms, such as the example in Fig. 1(b), are created by illuminating the sphere with a collimated laser beam, in this case the 3 mm diameter beam provided by a 10 mW HeNe laser (Uniphase) operating at a vacuum wavelength of 633 nm. Light scattered by the particle interferes with the unscattered portion of the incident beam to produce an interference pattern in the microscope's focal plane. This is magnified and projected onto a CCD camera (NEC TI324AII), which records holograms at video rates. Each frame in the video stream then is fit to the predictions of LorenzMie scattering theory (8) to obtain each particle's position in three dimensions, , its radius , and its complex index of refraction (7). Figure 1(b) also shows the computed image resulting from fitting to the experimental image. A sequence of such fits yields the particle's threedimensional trajectory, as plotted in Fig. 1(c).
Were the particle to come to equilibrium within the trap, the probability to find it within of would be given by the Boltzmann distribution, , where is the thermal energy scale at absolute temperature and is a normalization constant. This is the basis for a popular method to calibrate optical traps (9), with measured particle positions being compiled into , which then is inverted to yield . In the further approximation that the trap can be modeled as a threedimensional harmonic well centered at the origin, the threedimensional probability density can be factored into onedimensional contributions, , where . The trap is then characterized by the three effective spring constants, .
Figure 2(a) shows the result of applying the harmonic approximation to 10,000 measurements of a single particle's position in an optical tweezer powered by of light. The difference between the trap's apparent stiffness, , along the axis of the light's polarization and in the perpendicular direction is consistent with previous reports of polarization effects in optical trapping (10). The axial stiffness, is a factor of 5 smaller because axial intensity gradients are correspondingly weaker (11). Results in Fig. 2(b) for the same particle at different laser powers confirms that .
To quantify nonequilibrium effects in the particle's dynamics, we consider its trajectory in cylindrical coordinates, , centered on the trap's equilibrium position. The diffusing particle explores a roughly elliptical region in the plane. A statistical bias due to radiation pressure along the axis should appear as a tendency of the particle's trajectory to wind clockwise in the plane. To quantify this, we define a measure of the instantaneous circulation rate over the interval between video frames,
(1) 
which is positive for clockwise circulation and negative for retrograde motions. The angle brackets in Eq. (1) denote an average over the trajectory. To identify trends in despite large thermal fluctuations, we define
(2) 
which measures the accumulation of clockwise circulation in the particle's trajectory. The increasing trend in , plotted in Fig. 3(a), confirms the existence of a circulating steady state in which the particle completes one net cycle in roughly 9 s with a mean drift speed of 150 nm/s.
The measured circulation rate increases with laser power, as shown in Fig. 3(b). The particle also becomes increasingly well localized as the trap stiffens, however, so that the range of experimentally accessible laser powers is limited by our instrumental resolution to .
Observing a convective flux in the trapped particles' trajectories confirms the system's departure from equilibrium. The nature of the nonequilibrium state is clarified by considering an idealized model of the system. The probability flux induced by a force acting on a Brownian particle is , where is the particle's mobility and is its diffusion coefficient. For a sphere in water, , where is the water's viscosity. We assume that the system reaches steady state so that does not depend on time. Continuity of the probability density () then yields the FokkerPlanck equation
(3) 
The divergencefree flux still can support circulation,
(4) 
provided that .
We model the trap as a radially symmetric harmonic well with radiation pressure directed along the optical axis with a strength proportional to the local intensity:
(5) 
This model is parametrized the trap's stiffness, , and the scale of radiation pressure, , both of which are proportional to the laser's power, . and also by the effective width of the trap, , which depends on the quality of the focus and the size of the trapped particle (12). For a particle whose radius is substantially larger than the wavelength of light, (12), and we set in our model for the trap. Stable trapping requires the restoring force to exceed radiation pressure, so that may be taken to be small parameter independent of laser power. In that case, has a point of stable mechanical equilibrium at .
The lighter (brown) trace in Fig. 3(a) shows the circulation observed in a fourthorder RungeKutta Brownian dynamics simulation of a particle diffusing in the force field described by Eq. (5), using parameters obtained from Fig. 2, and setting . Both the slope of and the magnitude of its fluctuations are consistent with experimental results. This agreement both supports our interpretation of our experimental observations and also indicates that Eq. (5) is sufficiently detailed for quantitative comparisons with our measurements.
Expanding to first order in , , yields an approximate solution (13) to Eq. (3),
(6)  
(7) 
Equation (6) is the result that would be obtained for a particle at equilibrium in a harmonic potential. Equation (7) describes the lowestorder correction due to spatially nonuniform radiation pressure. Because this expression for is obtained in the stabletrapping limit, , it underestimates the probability for the particle to make largescale excursions from the trap. The results that follow therefore are conservative underestimates for the influence of nonconservative optical forces on a trapped particle's dynamics.
Even if the particle were to reach equilibrium in the force field , the distortion of by would affect measurements of colloidal forces calibrated by thermal fluctuation analysis. The firstorder correction not only displaces the center of the probability distribution downstream, but also broadens it. Consequently, the effective trap stiffness obtained from a probe particle's unloaded fluctuations systematically underestimate the forces exerted by the trap under load.
Our experimental results demonstrate, furthermore, that the distorted probability density does not come to equilibrium, but rather undergoes steadystate circulation. Expanding to lowest nonvanishing order in , we obtain the mean circulation rate (13)
(8)  
(9) 
This result can be generalized for an anisotropic trap by appropriately scaling the Cartesian coordinates in Eq. (5). In that case, Eq. (9) still holds with the effective spring constant set by the harmonic mean, . Given the results from Fig. 2, the measured circulation rates in Fig. 3(b) agree well with the predictions of Eq. (9) (dashed green curve) for , including the predicted dependence on laser power.
Equation (9) also predicts that the circulation rate scales with temperature as for . This confirms that the particle would not circulate at all were it not for thermal fluctuations. The optically trapped particle therefore is an exceptionally simple example of a Brownian motor (14) whose ability to perform work relies on rectification of thermal noise. Unlike previous optical implementations of thermal ratchets (15), the fountain of probability involves no timedependent driving, but rather is biased into motion by the nonconservative component of the optical force (3).
The fountain's efficiency as a Brownian motor can be estimated by computing the power dissipated into the water by viscous drag given the particle's mean drift velocity, (17). The result (13),
(10) 
depends strongly on laser power at low powers but becomes independent of in the strongtrapping limit. The asymptotic dissipation rate, , is a fixed proportion of the particle's thermallydriven velocity fluctuations. The system transduces fluctuations into heat through the bias imposed by radiation pressure. Its efficiency is limited by the material properties that determine , and thus by the requirement that the particle remain trapped.
Although this model is consistent with our measurements, other mechanisms also might cause a particle to circulate in an optical trap. For example, localized heating due to absorption of light at the focal point could create a toroidal convection roll that would advect the particle. This has been demonstrated with infrared optical tweezers operating at 1480 nm (18). The optical absorption coefficient of water is five orders of magnitude smaller at 532 nm, however (19). Scaling the previously reported (18) convection rate accordingly suggests that the maximum drift speed due to thermal convection in our system should be no greater than at the highest laser powers in out study. Thermal convection therefore appears unlikely to account for our observations.
If, indeed, steadystate circulation is an intrinsic feature of optically trapped particles' dynamics, the nonequilibrium effects we have identified should influence any measurement based on analysis of their thermal fluctuations. How such effects might have affected previously reported measurements remains to be determined. Establishing that a colloidal particle trapped in a static optical tweezer acts as a Brownian motor also creates new opportunities for research in nonequilibrium statistical physics. For example, the nonuniform intensity in an optical tweezer also exert positiondependent torques on the trapped particle that we have not considered here. The resulting nonequilibrium behavior in the particle's rotational degrees of freedom may share features in common with the fountainlike translational bias we have identified, and the two may be coupled in interesting ways. These considerations also may be extended to account for nonequilibrium effects in more general light fields.
We are grateful to Vincent Pereira and Alexander Grosberg for enlightening conversations. This work was supported by the National Science Foundation through Grant Number DMR0451589 and by a grant from Consolidated Edison.
References

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

(2)
A. Ashkin, Biophys. J. 61, 569 (1992).

(3)
Y. Roichman, B. Sun, Y. Roichman, J. AmatoGrill, and D. G. Grier, Phys. Rev. Lett. 100, 013602 (2008).

(4)
K. Tomita and H. Tomita, Prog. Theor. Phys. 51, 1731 (1974).

(5)
E. R. Dufresne and D. G. Grier, Rev. Sci. Instrum. 69, 1974 (1998); D. G. Grier, Nature 424, 810 (2003); M. Polin, K. Ladavac, S.H. Lee, Y. Roichman, and D. G. Grier, Opt. Express 13, 5831 (2005).

(6)
J. Sheng, E. Malkiel, and J. Katz, Appl. Opt. 45, 3893 (2006); S.H. Lee and D. G. Grier, Opt. Express 15, 1505 (2007).

(7)
S.H. Lee, Y. Roichman, G.R. Yi, S.H. Kim, S.M. Yang, A. van Blaaderen, P. van Oostrum, and D. G. Grier, Opt. Express 15, 18275 (2007).

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

(9)
E. L. Florin, A. Pralle, E. H. K. Stelzer, and J. K. H. Horber, Appl. Phys. A 66, S75 (1998).

(10)
A. Rohrbach, Phys. Rev. Lett. 95, 168102 (2005); A. R. Zahkarian, P. Polynkin, M. Mansuripur, and J. V. Moloney, Opt. Express 14, 3660 (2006); R. S. Dutra, N. B. Viana, P. A. Maia Neto, and H. M. Nussenzveig, J. Opt. A 9 (2007).

(11)
K. C. Vermeulen, G. J. L. Wuite, G. J. M. Stienen, and C. F. Schmidt, Appl. Opt. 45, 1812 (2006).

(12)
K. Ladavac, K. Kasza, and D. G. Grier, Phys. Rev. E 70, 010901(R) (2004); M. Pelton, K. Ladavac, and D. G. Grier, Phys. Rev. E 70, 031108 (2004); Y. Roichman, V. Wong, and D. G. Grier, Phys. Rev. E 75, 011407 (2007).

(13)
See EPAPS Document No. [number will be inserted by publisher] (2008), for details of the derivation of an optically trapped particle's circulation. For more information on EPAPS, see http://www.aip.org/pubservs/epaps.html.

(14)
P. Reimann, R. Bartussek, R. Haussler, and P. Hänggi, Phys. Lett. A 215, 26 (1996); P. Reimann, Phys. Rep. 361, 57 (2002); P. Reimann and P. Hänggi, Appl. Phys. A 75, 169 (2002).

(15)
L. P. Faucheux, L. S. Bourdieu, P. D. Kaplan, and A. J. Libchaber, Phys. Rev. Lett. 74, 1504 (1995a); L. P. Faucheux, G. Stolovitzky, and A. Libchaber, Phys. Rev. E 51, 5239 (1995b); S.H. Lee, K. Ladavac, M. Polin, and D. G. Grier, Phys. Rev. Lett. 94, 110601 (2005); S.H. Lee and D. G. Grier, Phys. Rev. E 71, 060102(R) (2005).

(16)
S.H. Lee and D. G. Grier, J. Phys.: Condens. Matter 17, S3685 (2006).

(17)
U. Seifert, Phys. Rev. Lett. 95, 040602 (2005).

(18)
D. Braun and A. Libchaber, Phys. Rev. Lett. 89, 188103 (2002).

(19)
G. M. Hale and M. R. Querry, Appl. Opt. 12, 555 (1973).