Comment: The effect of Mie resonances on trapping in optical tweezers
Abstract.
Recently, Stilgoe, et al., (Optics Express 16, 15039 (2008)) reported calculations of the force on an optically trapped sphere performed using the “optical tweezers toolbox”. This software suffers from numerical inaccuracies that lead to qualitative errors in the state diagram for stable trapping, particularly for spheres smaller than the wavelength of light.
An optical tweezer is a singlebeam optical trap that can localize objects in three dimensions using forces exerted by a strongly focused beam of light (1). Recently, Stilgoe et al. reported calculations of the axial trapping efficiency of optical tweezers for dielectric spheres, which emphasized the important role of Mie resonances on the stability and strength of optical traps, particularly for highindex particles (2). These calculations are based on the LorenzMie theory of light scattering by spheres (3), and so should provide exact results for all sphere sizes and compositions. This approach also is numerically efficient, and can be generalized to compute the forces and torques due to arbitrary light fields (4) on objects of arbitrary shape (5); (6). The calculations reported in Ref. (2) were performed with the “optical tweezers toolbox” software suite (7) described in Ref. (5), which is formally correct, but relies on numerically inaccurate Hankel functions provided by Matlab® (versions R2007a through R2008b, The MathWorks, Natick, MA). Consequently, some of the results reported in Ref. (2) are incorrect, particularly for very small spheres.
The incident light field, , the scattered wave, , and the field within the sphere, , are all described in LorenzMie theory as expansions in vector spherical harmonics, (3)
(1)  
(2)  
(3) 
where is the wavenumber of light in a medium of refractive index .
The incident field's coefficients, and , were computed in Ref. (2) through the overdetermined pointmatching method (8) but can be derived equally well with the DebyeWolf path integral formalism (4). The scattering coefficients, and , depend on the on the size and composition of the sphere, and on the structure of the light field illuminating it. The interior field's coefficients, and , similarly are computed to satisfy continuity conditions at the sphere's surface (6):
(4)  
(5) 
where , is the particle's relative refractive index, and where and are RicattiBessel functions.
Both and contribute to the Maxwell stress tensor on the sphere's surface, whose integral yields the force. References (2) and (5) follow Crichton and Marston (9) by performing the integral analytically. The result for the axial trapping efficiency in an optical trap of total power is then (9)
(6) 
Equation (10) of Ref. (5) misstates this result, inconsistently taking the real part of the second term. The distributed software, however, is formally correct (7).
The central results of Ref. (2) rely on numerical convergence of the sums in Eq. (6). Unfortunately, the scattering coefficients, and , are notoriously difficult to compute accurately (3); (6) because the RicattiBessel functions in which they are expressed have numerically unstable recurrence relations. The expansions' accuracy can be assessed by checking the fields' continuity on the sphere's surface using Eqs. (3) through (5). In our implementation, the relative error in the computed fields is smaller than so that the relative error in the force is less than .
Figure 1 shows the maximum axial trapping efficiency, as a function of the sphere's radius, and relative refractive index, , in an optical tweezer of numerical aperture 1.2. These conditions correspond to those of Fig. 2(b) of Ref. (2). The principal qualitative difference is that Fig. 1 predicts stable trapping for spheres smaller than , while Ref. (2) does not. The present results are consistent with an analytic computation of the force on a Rayleigh particle in a Gaussian trap (10), and thus would appear to be correct. They also agree quantitatively with results obtained by numerically integrating the Maxwell stress tensor over the surface of the sphere (4), which confirms the correct implementation of Eq. (6). Although these results suggest that even vanishing small spheres can be trapped, the maximum trapping force is limited in practice by saturation of the sphere's polarization, a nonlinear effect not considered in the linear lightscattering formalism discussed here and in Ref. (2).
Explicitly computing the surface fields, such as the examples presented in Fig. 1(a) through (e), reveals smooth symmetric surface field distributions for stably trapped spheres and smooth but asymmetric distributions for unstable conditions. The results in Fig. 1 therefore appear to be free from obvious numerical artifacts.
Major qualitative and quantitative discrepancies we observe with the results of Ref. (2) can be attributed to the computational toolbox's use of the standard Matlab routine besselh(), which implements the Hankel function that is used to compute the spherical Hankel function . Taking as a test case, accurate calculation of requires the sum in Eq. (6) to be computed to at least (3); (6), with . The corresponding result obtained with Matlab is , whose real component is in error by a factor of ten. Such numerical errors diminish with increasing and are negligible for . Other quantitative differences between Fig. 1 and Fig. 2(b) of Ref. (2) are likely due to slight differences in the description of the incoming wave.
This work was supported in part by the NSF through Grant Number DMR0606415 and in part by a grant from the Keck Foundation. B.S. acknowledges support from a Kessler Family Foundation Fellowship.
References

(1)
A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, ``Observation of a singlebeam gradient force optical trap for dielectric particles,” Opt. Lett. 11(5), 288–290 (1986).

(2)
A. B. Stilgoe, T. A. Nieminen, G. Knöener, N. R. Heckenberg, and H. RubinszteinDunlop, “The effect of Mie resonances on trapping in optical tweezers,” Opt. Express 16, 15,039–15,051 (2008).

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

(4)
B. Sun, Y. Roichman, and D. G. Grier, “Theory of holographic optical trapping,” Opt. Express 16, 15,765–15,776 (2008).

(5)
T. A. Nieminen, V. L. Y. Loke, A. B. Stilgoe, G. Knoner, A. M. Branczyk, N. R. Heckenberg, and H. RubinszteinDunlop, “Optical tweezers computational toolbox,” J. Opt. A 9, S196–S203 (2007).

(6)
M. I. Mishchenko, L. D. Travis, and A. A. Lacis, Scattering, Absorption and Emission of Light by Small Particles (Cambridge University Press, Cambridge, 2001).

(7)
T. A. Nieminen, V. L. Y. Loke, A. B. Stilgoe, Y. Hu, G. Knoener, and A. M. Branczyk, “Optical tweezers toolbox 1.1,” Tech. rep., University of Queensland (2007). URL http://www.physics.uq.edu.au/people/nieminen/software.html.

(8)
T. A. Nieminen, H. RubinszteinDunlop, and N. R. Heckenberg, “Multipole expansion of strongly focussed laser beams,” J. Quant. Spec. Rad. Transfer 7980, 1005–1017 (2003).

(9)
J. H. Crichton and P. L. Marston, “The measurable distinction between the spin and orbital angular momenta of electromagnetic radiation,” Elec. J. Diff. Eq. Conf. 04, 37–50 (2000).

(10)
Y. Harada and T. Asakura, “Radiation forces on a dielectric sphere in the Rayleigh scattering regime,” Opt. Commun. 124, 529–541 (1996).