Performance of optical traps with geometric aberrations
Abstract.
We assess the influence of geometric aberrations on optical traps' inplane performance by studying the dynamics of trapped colloidal spheres in deliberately distorted holographic optical tweezers. The traps' lateral stiffness turns out to be insensitive to moderate amounts of coma, astigmatism and spherical aberration. Moreover, holographic aberration correction enables us to compensate inherent shortcomings in the optical train, thereby adaptively improving its performance. We also demonstrate the effects of geometric aberrations on the intensity profiles of optical vortices, whose readily measured deformations suggest a method for rapidly estimating and correcting geometric aberrations in holographic trapping systems.
An optical tweezer exploits forces exerted by intensity gradients in a single strongly focused beam of light to trap mesoscopic objects in three dimensions (1). The ability to noninvasively exert and monitor forces in microscopic systems has led to exciting new insights in a wide range of fields, including the functioning of biological molecular motors and the statistical mechanics of thermally driven systems. Increasingly, these applications require optical traps with highly optimized and precisely characterized properties. The threedimensional potential energy well associated with an optical tweezer is determined by the intensity distribution near the focal point, and thus on the detailed structure of the light's wavefronts. The sharpest focus, which yields the strongest intensity gradients and thus the best traps, requires a converging beam with spherical wavefronts, and is said to be diffraction limited. Deviations from this ideal structure are dubbed aberrations and degrade the trap's performance.
It has bean shown both experimentally (2); (3); (4); (5); (6) and theoretically (7); (18) that optical tweezers' axial trapping efficiency is sensitive to the presence of spherical aberration. Less has been published about spherical aberration's influence on lateral trapping efficiency (8); (25), and less still about the influence of other geometric aberrations. This is a growing concern as optical trapping systems become increasingly sophisticated and are combined with complementary optical techniques.
Here we report measurements of lateral trap stiffness under the influence of spherical aberration, coma and astigmatism, using computergenerated holograms to introduce controlled aberrations and digital particle tracking (9); (10) to characterize the traps' performance. Not only does this approach provide a means to measure aberrations' influence on optical trapping efficiency, it also provides the basis for realtime adaptive optimization. Holographic optical tweezers (HOT) (11); (12); (13); (14); (15); (16); (10) are especially suited for such studies because their use of computergenerated holograms to structure the wavefronts of light affords a quick and flexible way to introduce aberrations without otherwise modifying the optical train. The degraded beam is subsequently used to create an array of traps whose performance can be assessed in parallel.
Our HOT implementation, shown schematically in Fig. 1, uses light from a frequencydoubled Nd:YVO laser (Coherent Verdi) expanded by a Galilean telescope to fill the 25 mm diagonal face of a reflective liquid crystal spatial light modulator (SLM). (Hamamatsu X826716 PPM). The SLM can impose phase shifts, , at each position in a array. The phasemodified beam is relayed by a second telescope and a wavelengthselective polarizing beam splitter to the input pupil of a , NA 1.4, Plan Apo oil immersion objective mounted in a Nikon TE2000U inverted optical microscope, which focuses the light into optical traps. Trapped samples are imaged onto a CCD camera using the microscope's standard optical train, with the imaging illumination passing through the beam splitter.
To assess the holographic traps' performance under the influence of geometric aberrations, we studied the dynamics of silica beads 1.53 in diameter (Duke Scientific Lot 5238) dispersed in water and trapped in a array of optical tweezers. We used standard methods of digital video microscopy (9) to track the trapped particles' motions with spatial resolution at intervals. The resulting trajectories were analyzed with statistically optimal methods (10) to extract the traps' stiffness along and perpendicular to the video scan lines while simultaneously measuring the spheres' hydrodynamic radii.
The five primary geometric aberrations in a beam of light can be characterized with Zernike's orthogonal polynomials (17):
spherical aberration  
coma  
astigmatism  
curvature of field  
distortion  (1) 
where is the radial coordinate in units of the input pupil's radius, , is the polar angle around the optical axis, and where the coefficients are measured in units of the wavelength of light. The three angles, , and , establish the orientation of coma, astigmatism and distortion, respectively. All five primary geometric aberrations can be introduced independently by adding their phase representations to the hologram encoding the array of traps, as shown in Fig. 1.
Figure 2 shows the effect of coma on optical traps' structure and performance. Coma is typically introduced by failure to center lenses along the optical axis. This is a concern, therefore, for optical trapping systems that translate their traps by displacing one or more lenses. The nine axial crosssections through holographically projected traps in Fig. 2(a) demonstrate that coma ruins the beam's radial symmetry at large distances () from the focus, but does not have an obvious effect on the intensity distribution in the plane of best focus. Even so, a typical trajectory in an aberratd trap such as the inset to Fig. 2(b) shows a slight diagonal elongation that is not evident for acomatic traps.
To quantify the traps' performance, we model each as a harmonic potential well
(2) 
characterized by distinct stiffnesses and along and transverse to the camera's scan direction, respectively. A measured trajectory consisting of time intervals then yields the associated trap's spring constants and the sphere's viscous drag coefficient as (10)
(3)  
(4) 
where
(5)  
(6) 
and where
(7) 
is the autocorrelation at lag of the trajectory's th component. Error estimates for and also follow from the methods of Ref. (10). Equations (3) through (7) yield accurate results provided the sampling interval is comparable to the viscous relaxation time . With the unaberrated traps powered by 6.25 mW per trap, we obtain and and . The variances in these results include the measured (10) 10% variation in the holographically projected traps' intensities and the rated 10% polydispersity in the spheres' radii. This yields . In all cases, the thermally driven excursions at temperature exceed the variance due to measurement error, , so that the approximations (10) used in deriving Eqs. (3) and (4) are satisfied. With these conditions satisfied, video frames obtained over 60 suffice to measure the traps' stiffnesses to within a relative error of 8%. As a final consistency check, we ensure that the spheres' viscous drag coefficients measured along different direction agree even if the traps themselves are anisotropic.
To best compile results from the entire array of traps, we normalize results for each trap by their values measured without imposed aberrations and then average over the array. The trapping light is linearly polarized, with its plane of polarization aligned with the direction to within . The direction of polarization could influence the traps' performance in our highNA optical train. The potential wells' measured curvatures along the and directions, however, indistinguishable in the absence of applied aberrations. This is consistent with recent advances in the theory of optical trapping, which suggest that polarization effects in highNA optical trapping systems only are important for smaller objects (18).
The results in Fig. 2(b) demonstrate that coma measurably degrades the traps' lateral stiffness despite its comparatively modest influence on the individual traps' intensity distributions. Coma's influence can be made more apparent by employing the SLM also as a mode converter. For example, imposing the phase profile transforms a conventional pointlike optical tweezer into a ringlike optical vortex (19); (20); (21) whose radius is determined (22); (23) by the integer winding number . In the absence of aberrations, a wellfocused optical vortex appears as a uniformly bright ring of light. By distorting the helical phase profile that defines an optical vortex, geometric aberrations distort its shape and redistribute its intensity. The structure of coma, for instance, displaces an optical vortex's center away from the optical axis, while the twofold modulation of astigmatism extends it into an elliptical shape (24). These effects are accentuated by these distortion's radial dependence. The comatic optical vortices in Fig. 2(c) and 2(e) therefore are not just displaced relative to the optical axis, but also are elongated. The example with no imposed distortions in Fig. 2(d) reflects the inherent aberrations in our optical train. The distortions in such images can serve as a metric for estimating the magnitude, , and orientation, of intrinsic geometric aberrations. Once measured, these can be removed by inverting the signs of the coefficients in Eqs. (Performance of optical traps with geometric aberrations) and adding the resulting phase modulation to any trapforming holograms.
A close examination of the data in Fig. 2(b) suggests that the traps' performance might be improved by deliberately adding a small amount of negative coma. This is consistent with the observed distortion of the optical vortex in Fig. 2(d).
Astigmatism, which can result from having lenses canted in the optical train, has a more obvious effect on optical tweezers' structure and function. The crosssections through an astigmatic trap in Fig. 3(a) demonstrate that the beam focuses along the and axes in different planes. The result is that a sphere encounters a substantially different lateral confining potential as it diffuses axially. This gives rise to the crossshaped trajectory in the inset to Fig. 3(b). Although the derivation of Eq. (3) does not account for this complicated structure, the results obtained for and are consistent with astigmatism's degrading lateral trapping performance. Furthermore, because the spheres are displaced downstream by radiation pressure, the traps appear softer along the direction in which the downstream lobe is longer. This changes with the sign of the coefficient , as shown by the data in Fig. 3(b).
The inherent astigmatism in our system can be estimated by locating the symmetry axis of Fig. 3(b), and is roughly . This is consistent also with the images of astigmatic optical vortices in Figs. 3(c), (d) and (e). These clearly show biaxial extension due to astigmatism's twofold modulation.
The other primary geometric aberrations have more subtle effects on optical traps' efficacy and on optical vortices' structure. Spherical aberration, for example, is characterized by phase terms proportional to and . The quadratic term serves mostly to displace the trap's plane of best focus along the direction. To emphasize spherical aberration's effect on the optical tweezers' lateral stiffness, we imposed only the quartic term. This enables us to adjust the degree of aberration without shifting the focal plane, which frees our measurement from several possible artifacts including changes in the optical train's inherent spherical aberration (8); (25) and changes in the spheres' hydrodynamic coupling to the wall (26). The data in Fig. 4 reveal that even large amounts of this quartic distortion have very little influence on the traps' inplane performance. This contrasts with spherical aberration's demonstrated degradation of axial trapping stiffness (2); (3); (4); (5); (6). Even a wellaligned optical train can suffer from this defect if the sample's index of refraction differs from the design value. In this case, deliberately adding negative spherical aberration can restore the traps' performance (8); (5); (25), and might account for the slight improvement we observe.
Optimal statistical analysis of optically trapped particles' thermally driven trajectories provides rapid characterization of the traps' performance (10). When combined with direct visualization of the optical train's aberrations through the structure of projected optical vortices, this provides a basis for adaptively assessing and optimizing holographic optical trapping. Optimization takes the form of an estimate for the phase modulation required to correct wavefront distortions that can be added to trapforming holograms. Although aberrations are known to substantially degrade axial trapping performance, we have demonstrated that their effect on lateral trapping is less pronounced, with astigmatism and coma apparently requiring the most attention.
This work was supported by the National Science Foundation through Grant Number DBI0233971, through the Research Experiences for Undergraduates program of the NSF, and by the donors of the Petroleum Research Fund of the American Chemical Society.
References

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

(2)
S. Hell, G. Reiner, C. Cremer, and E. H. K. Stelzer, J. Microsc. 169, 391 (1993).

(3)
H. Felgner, O. Muller, and M. Schliwa, Appl. Opt. 34, 977 (1995).

(4)
M. J. Booth, M. A. A. Neil, and T. Wilson, J. Microsc. 192, 90 (1998).

(5)
T. Ota, T. Sugiura, S. Kawata, M. J. Booth, M. A. A. Neil, R. Juškaitis, and T. Wilson, Jpn. J. Appl. Phys. 42, L701 (2003).

(6)
K. C. Neuman, E. A. Abbondanzieri, and S. M. Block, Opt. Let. 30, 1318 (2005).

(7)
E. Fällman and O. Axner, Appl. Opt. 42, 3915 (2003).

(8)
P. C. Ke and M. Gu, J. Mod. Opt. 45, 2159 (1998).

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

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

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

(12)
M. Reicherter, T. Haist, E. U. Wagemann, and H. J. Tiziani, Opt. Lett. 24, 608 (1999).

(13)
J. Liesener, M. Reicherter, T. Haist, and H. J. Tiziani, Opt. Comm. 185, 77 (2000).

(14)
E. R. Dufresne, G. C. Spalding, M. T. Dearing, S. A. Sheets, and D. G. Grier, Rev. Sci. Instr. 72, 1810 (2001).

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

(16)
D. G. Grier, Nature 424, 810 (2003).

(17)
M. Born and E. Wolf, Principles of Optics (Cambridge University Press, Cambridge, 1980), chap. 9, pp. Table XXII, page 470.

(18)
D. Ganic, X. Gan, and M. Gu, Opt. Express 12, 2670 (2004).

(19)
H. He, N. R. Heckenberg, and H. RubinszteinDunlop, J. Mod. Opt. 42, 217 (1995).

(20)
N. B. Simpson, L. Allen, and M. J. Padgett, J. Mod. Opt. 43, 2485 (1996).

(21)
K. T. Gahagan and G. A. Swartzlander, Opt. Lett. 21, 827 (1996).

(22)
J. E. Curtis and D. G. Grier, Phys. Rev. Lett. 90, 133901 (2003a).

(23)
S. Sundbeck, I. Gruzberg, and D. G. Grier, Opt. Lett. 30, 477 (2005).

(24)
J. E. Curtis and D. G. Grier, Opt. Lett. 28, 872 (2003b).

(25)
E. Theofanidou, L. Wilson, W. J. Hossack, and J. Arlt, Opt. Comm. 236, 145 (2004).

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