Universal, strong and long-ranged trapping by optical conveyors

Transcription

1 Universal, strong and long-ranged trapping by optical conveyors David B. Ruffner and David G. Grier Department of Physics and Center for Soft Matter Research, New York University, New York, NY 10003, USA Optical conveyors are active tractor beams that selectively transport illuminated objects either upstream or downstream along their axes. Formed by the coherent superposition of coaxial Bessel beams, an optical conveyor features an axial array of equally spaced intensity maxima that act as optical traps for small objects. We demonstrate through measurements on colloidal spheres and numerical calculations based on the generalized Lorenz-Mie theory that optical conveyors interferometric structure endows them with trapping characteristics far superior to those of conventional optical tweezers. Optical conveyors form substantially stiffer traps and can transport a wider variety of materials over a much longer axial range. INTRODUCTION: OPTICAL CONVEYORS Long a staple of science fiction [1], tractor beams are traveling waves that transport illuminated objects back to their sources, opposite to the direction of energy flow. Such retrograde transport has been demonstrated in optical trapping experiments using collimated solenoid beams [], interfering plane waves [3], and optical conveyors created with coaxial Bessel beams [4 6]. Among the more surprising proposals arising from these initial studies is that optical conveyors might act as universal tractor beams, transporting small objects at uniform speed regardless of their composition [4]. This contrasts with the performance of conventional optical tweezers, which can only trap bright-seeking objects over a limited domain of size, shape and composition [7 9]. Another surprising suggestion is that optical conveyors might exceed the trapping stiffness even of diffraction-limited optical tweezers because their intensity actually vanishes at regular intervals along the optical axis [4]. Here, we report the results of a head-to-head comparison between the trapping characteristics of optical tweezers and optical conveyors involving both experimental measurements on micrometer-scale spheres and also numerical evaluations of optical forces using the generalized Lorenz-Mie theory. In confirming the superior performance of optical conveyors, these studies also enable us to establish the axial range over which optical conveyors can usefully transport material, and provide guidance for developing long-ranged tractor beams. The electric field of a monochromatic optical conveyor of angular frequency ω, linearly polarized along ˆx and propagating along ẑ through a medium of refractive index n m, is the superposition of two Bessel beams [4 6], E(r, t) = 1 [ ] E 0e iωt b 1 (kr) + e iϕ(t) b (kr), (1) each of which may be described as a conical superposition of plane waves [6], b j (kr) = π 0 ˆɛ(θ j, φ) e ik(θj,φ) r dφ, (a) where ˆɛ(θ, φ) = cos φ ˆθ + sin φ ˆφ (b) is the polarization of a plane wave incident on the optical axis at polar angle θ and azimuthal angle φ, and where k(θ, φ) = k(sin θ cos φ ˆx + sin θ sin φ ŷ + cos θ ẑ) (c) is the corresponding wave vector. In the paraxial approximation, which is appropriate for long-range tractor beams, Eq. () reduces to ( ) b j (kr) J 0 1 ηj kr e iηjkz ˆx, (3) where η j is related to the Bessel beam s cone angle by η j = cos θ j. The approximate expression in Eq. (3) differs from the exact expression in Eq. () by terms involving higher-order Bessel functions[0, 3], which vanish on the optical axis. The numerical results developed in Sec. are based on Eq. (). Analytical results developed in Sec. are more readily obtained from the approximation in Eq. (3). The two Bessel beams comprising an optical conveyor share the same amplitude E 0, frequency ω and polarization along ˆx, but differ in their relative phase, ϕ(t), and also in their axial wave numbers, η j k, which are reduced from the plane wave value, k = n m ω/c, by the dimensionless factor η j (0, 1]. Here, c is the speed of light in vacuum. The superposition is usefully characterized by the mean convergence factor η = (η 1 + η )/ and the difference η = η 1 η. The upper limit of an optical conveyor s range, R, is set by the non-diffracting range of the most strongly converging Bessel beam [10], R A cot θ 1 = A η 1 1 η 1 (4) where A is the radius of the beam s aperture. Longer ranges can be achieved with values of η approaching 1. Figures 1(a) and 1(b) show the measured [11] threedimensional intensity distribution of optical conveyors with equal values of η and differing values of η that were projected with the holographic optical trapping

2 FIG. 1. (a) Experimental reconstruction of an optical conveyor η = 0.8, η = 0.04, and (b) η = 0.8, η = (c) Measured trajectory of a 1.5 µm-diameter silica sphere trapped in one of the intensity maxima in (b). (d) Trajectory of the same particle trapped in a conventional optical tweezer projected by the same instrument with the same peak intensity. (e) Measured transverse and axial stiffness as a function η. (f) Stiffness ratio, kz /kr, as a function of η. Solid curves in (e) and (f) represent predictions of the Lorenz-Mie theory. Shaded bands represent uncertainty in the measured size and refractive index of the trapped silica sphere. Highlighted plot symbols represent results from (c) and (d). Dotted horizontal lines represent the measured performance of the optical tweezer. The dashed horizontal line in (f) shows the theoretical limit for optical tweezer performance. The shaded region above this line represents the optical conveyor s superior performance for optical micromanipulation. technique [4, 1]. These beams were powered by a linearly polarized diode-pumped solid-state laser (Coherent Verdi) operating at a vacuum wavelength of λ = 53 nm that was shaped by a liquid crystal spatial light modulator (Holoeye Pluto) before being projected with a microscope objective lens (Nikon Plan-Apo, 100, numerical aperture 1.4, oil immersion). The mode projected by this method is not the simple superposition of Bessel beams, but rather incorporates contributions from small range of axial wave numbers around ηj k [4]. Equation (1) therefore should be considered an idealized model for the actual beam. The intensity distribution, I(r, t) = 1 nm 0 c E(r, t), has maxima at axial positions zn (t) = πn + ϕ(t), η k (5) each of which can act as an optical trap for a small object. Here, 0 is the permittivity of space. Varying ϕ(t) as a function of time moves these extrema, and thus conveys trapped objects along the beam. Increasing η reduces the spacing between maxima, as shown in Figs. 1(a) and 1(b), and so provides control over the optical force profile. TRAP STIFFNESS We previously have proposed [4] that optical conveyors should make better traps than conventional optical tweezers because their intensity vanishes altogether between maxima. The data in Figs. 1(c) and 1(d) demonstrate this to be true. Figure 1(c) shows the measured trajectory of a colloidal silica sphere diffusing through water in one of the potential energy wells of the static optical conveyor from Fig. 1(b). The optical conveyor has a peak intensity of 79 mw/µm, as measured with imaging photometry. Holographic characterization [13 15] reveals the sphere s radius to be ap = ± µm, and its refractive index to be np = 1.44 ± Holographic tracking [13, 16] yields the sphere s position with 1 nm precision in-plane and 3 nm resolution axially [15, 17] at 16.7 ms intervals. Figure 1(d) shows the same sphere diffusing in a conventional optical tweezer with the same peak intensity projected by the same instrument. In both cases, the trapped particle explores the optical force landscape under the influence of random thermal forces. The optical conveyor restricts the particle s axial excursions to less than half the range of the optical

3 3 tweezer, resulting in a nearly isotropic trajectory. This suggests that an optical conveyor makes a substantially stiffer trap, even though the.8:1 aspect ratio of the optical tweezer s trajectory approaches the theoretical limit for a diffraction-limited Gaussian trap [18]. Modeling the traps as cylindrically symmetric harmonic potential energy wells, U(r) = 1 k rr + 1 k zz, (6) we may estimate the transverse and axial trap stiffness, k r and k z from the particle s trajectory using thermal fluctuation analysis [19, 1]. These measurements are performed for fixed values of the relative phase, ϕ(t), so that the optical traps do not move during the measurement. For the trajectory in Fig. 1(c), we obtain k r =. ± 0.1 pn/µm and k z = 0.89 ± 0.03 pn/µm, and an anisotropy of k r /k z = 1.57 ± The equivalent results for the optical tweezer are k r =.8 ± 0.4 pn/µm, k z = 0.41±0.08 pn/µm and k r /k z =.6±0.4. The optical conveyor performs as well as the optical tweezer in the transverse direction and is nearly twice as stiff along the axis. The performance of single beam optical traps typically is limited by their axial trapping ability. Figure 1(e) shows how the optical conveyor s stiffness varies with η for fixed η, and compares this with the performance of a diffraction-limited optical tweezer. Highlighted points correspond to the data from Figs. 1(c) and 1(d). The optical conveyor s axial stiffness exceeds that of a diffraction-limited optical tweezer by as much as a factor of two. This advantage is emphasized by the ratio between axial and transverse stiffness plotted in Fig. 1(f). Dashed horizontal dashed lines represent the measured and theoretical maximum performance of a Gaussian optical tweezer. Optical conveyors exceed this performance for η > 0.06, which corresponds to an axial period less than 6.7 µm. Generalized Lorenz-Mie theory Optical conveyors superior trapping performance is consistent with predictions of generalized Lorenz-Mie theory [6, 5], which are plotted as continuous curves in Figs. 1(e) and 1(f). For these calculations, each Bessel beam is expanded as a series, b j (r) = n n=1 m= n [ a mn (θ j ) M (1) nm(kr) ] +b mn (θ j ) N (1) nm(kr), (7) in the vector spherical harmonics, M (1) nm(kr) and N (1) nm(kr), that constitute the natural basis for transverse electric (TE) and transverse magnetic (TM) waves, respectively [6]. The expansion coefficients, a mn (θ j ) and b mn (θ j ), have been reported previously for individual Bessel beams []; they can be superimposed according to Eq. (1) to create an optical conveyor with specified values of η, η and relative phase ϕ(t). The associated expansion of the scattered field, E s (r, t) = E 0 e iωt n {[ ] r mn (θ 1 ) + e iϕ(t) r mn (θ ) n=1 m= n + [ s mn (θ 1 ) + e iϕ(t) s mn (θ ) ] M (3) nm(kr) } N (3) nm(kr), (8a) is obtained by matching boundary conditions for the electric and magnetic fields at the sphere s surface [3 5, 7 9]. The expansion coefficients, r mn (θ j ) = a n a mn (θ j ) and (8b) s mn (θ j ) = b n b mn (θ j ), (8c) are related to the expansion coefficients of the incident field by the standard Mie coefficients, a n and b n, [7], which depend on the radius and refractive index of the sphere. The superposed incident and scattered fields, E(r, t) + E s (r, t), contribute to the Maxwell stress tensor, T(r, t), whose integral over a closed surface provides an estimate for the optically-induced force on a sphere centered at r: F(r, t) = ˆn T(r, t) dr, (9) S where ˆn is the unit normal to the surface S enclosing the sphere. In practice, the integral is computed directly from the combined expansion coefficients using established techniques [7 9, 30]. This calculation must also be repeated for each value of the relative phase ϕ(t). A trap s position, r 0 (t), for a given value of ϕ(t) is calculated numerically as a solution of F(r 0 (t), t) = 0. The traps s effective stiffness along ˆr ν then is calculated as k ν = ν F ν (r, t) r=r0(t). (10) Results for the predicted transverse and axial stiffness are plotted in Fig. 1(e) for conveyor beams with η = 0.8 and η ranging up to 0.13, and for values of the sphere s radius and refractive index obtained with holographic microscopy. The curves are scaled to a maximum intensity of 41 mw/µm, which is less than the measured value presumably because holographically projected Bessel beams are not uniformly bright. This discrepancy does not affect the computed ratio of stiffnesses in Fig. 1(f), which also agree well with experimental results, with no adjustable parameters.

4 4 Dipole approximation Generalized Lorenz-Mie theory is useful for computing the forces on specific particles in particular traps. To assess trends in optical traps capabilities, we invoke the dipole approximation in which the time-averaged force [31, 3], { } F(r, t) = 1 3 R E ν (r, t) Eν(r, t), (11) α e ν=1 is proportional to the object s polarizability α e = α e + iα e. The polarizability of a dielectric sphere is related to its size and refractive index by the Clausius-Mossotti- Draine relation [33] α e = 4πɛ 0n mk a 3 p 1 i, (1) 3 Kk3 a 3 p where K = (n p n m)/(n p + n m) is the Lorentz-Lorenz factor. Absorptivity increases the imaginary part of α e. Conductivity contributes an imaginary part to K. The dipole approximation typically applies in the Rayleigh limit, for particles much smaller than the wavelength of light. An object near the axis of an optical conveyor experiences an axial force F z (r, t) E α e η sin ( Φ(z, t) ) + and a transverse force F r (r, t) E 0 α e kr α e η cos ( 1 Φ(z, t) ) (13) (1 η 14 η ) cos ( 1 Φ(z, t) ) α e η η kr sin ( Φ(z, t) ), (14) where Φ(z, t) = η kz ϕ(t), and where we have assumed kr 1. These results are obtained by averaging over times long compared with the optical cycle but short compared with variations in ϕ(t). The particle is trapped where the force vanishes, which occurs at axial positions Z n (t) that are displaced from the intensity maxima by an amount that depends on the particle s light-scattering properties, Z n (t) z n (t) = ηk tan 1 ( α e α e ) η. (15) η Equations (10) and (16) then yield the traps axial stiffness, k z = 1 4 α e ke 0 η, (16) FIG.. (a) Trap stiffness as a function of particle size for silica spheres in the optical conveyor from Fig. 1(b). Predictions from Lorenz-Mie theory are plotted as solid curves, and the corresponding results in the dipole approximation are plotted as dashed curves. Discrete points show experimental results obtained from the data in Fig. 1(c). which is strictly positive. An optical conveyor thus can trap and transport any dipolar particle, regardless of its light-scattering characteristics. This universal, materialindependent trapping capability does not require feedback [34] or fine tuning of the beam s properties [35]. Decreasing the inter-trap separation by increasing η enhances intensity gradients and thus increases trap stiffness. Picking the largest possible value of η to optimize axial trapping is not necessarily the best strategy. An optical conveyor s transverse stiffness, η α ( e 1 η 1 4 k r = k η) α e z 1 α e η + α e, (17) η can vanish or even change sign as η increases. The dependence of trap stiffness on particle size predicted by Eqs. (1), (16) and (17) is plotted as dashed curves in Fig. for silica spheres in an optical conveyor with η = 0.8 and η = Results from the dipole approximation agree well with the Lorenz-Mie results that are plotted as solid curves in Fig., at least for ka p 1. The Lorenz-Mie predictions, in turn, agree quantitatively with results from Fig. 1, which are plotted as discrete points. The dipole approximation severely underestimates the optical conveyor s stiffness for ka p > 1, and suggests that

5 5 the particle studied in Fig. 1 would not have been stably trapped in the transverse direction. Rather than displaying a single crossover to instability, the generalized Lorenz-Mie result displays limited domains of instability for particles in particular size ranges. Results from the dipole approximation therefore are useful for establishing strict lower bounds on the performance of optical conveyors. For example, requiring stable transverse trapping (kr > 0) according to Eq. (17) establishes an upper bound on η and, through Eq. (4), a lower bound on the range over which an optical conveyor can transport small objects: 00 A 0 0 < R. (18) R 00 The practically accessible range of transport may be substantially greater than R for particles larger than the wavelength of light. Comparison with optical tweezers Conventional optical tweezers do not share optical conveyors universal trapping ability. To show this, we model an optical tweezer as a focused Gaussian beam whose axial electric field profile is [36, 37] EG (z, t) = E0 p zr z + zr eikz eiζ(z) e iωt ˆ, (19) where zr R /(ka ) is the Rayleigh range of a Gaussian beam converging at distance R from an aperture of radius A, and where ζ(z) = tan 1 (z/zr ) is the Gouy phase. Because we are interested in long-ranged axial transport, we assume R > A. The axial component of the associated force, 00 FG (z) 1 z0 k z zr, (0) = zr ) E0 (z + zr forms a trap only for particles satisfying 0 > 4kzR (kzr 1). 00 (1) This condition is most easily satisfied in strongly converging beams for which zr is small. The condition on is qualitatively consistent with numerical studies [7 9] which describe the difficulty of trapping high index or absorbing particles with laser tweezers. TRANSPORT RANGE The ultimate rate of an optical conveyor is limited by the non-diffracting range of the constituent Bessel beams. FIG. 3. (a) Experimental reconstruction of an optical conveyor with η = 0.96 and η = (b) Trajectory of a 1.5 µmsilica sphere transported by the same optical conveyor over a range of 66 µm. (c) Axial stiffness as a function of transport range. Solid curves show Lorenz-Mie predictions for optical conveyors transporting large (1.5 µm-diameter) and small (0.036 µm-diameter, kap = 0.5) silica spheres through water. Dashed curves show corresponding results for optical tweezers. Discrete symbols show results for the optical tweezer in Fig. 1(c) and the optical conveyors in Figs. 1(d) and 3(a). Equation (1), by contrast, establishes an inherent upper limit on an optical tweezer s range v s u 0 Au t () RG = 00 At ranges R beyond RG, radiation pressure overwhelms trapping forces due to axial intensity gradients and ejects the particle. The upper bound of an optical tweezer s range therefore is smaller than the lower bound of an optical conveyor s. This means that optical conveyors can transport objects over substantially longer axial ranges, and at much lower numerical apertures. The price for long-ranged transport is paid in the strength of an optical conveyor s traps. Figure 3(a) is a volumetric reconstruction of an optical conveyor with η = 0.96 and η = 0.04 that was projected with a 60 objective lens (Nikon Plan-Apo, numerical aperture 1.4 oil immersion). Holographic projection limits the conveyor s range to R = 70 µm, and thus suggests an effective numerical aperture of 0.43 given the convergence angle of θ = 19. The holographically measured trajectory plotted in Fig. 3(b) shows this optical conveyor trapping and transporting a 1.5 µm-diameter silica sphere over 66 µm. The hologram at the bottom of this figure was recorded when the sphere was located at the position indicated by the sphere in Fig. 3(b). This measurement demonstrates that an optical conveyor can transport objects along its entire length, even at low nu-

6 6 merical aperture. Figure 3(c) shows how the traps stiffness falls off with range. The solid and dashed curves are Lorenz-Mie calculations of the axial trap stiffness of optical conveyors and optical tweezers, respectively, and are scaled for a peak intensity of 41 mw/µm. One set of curves is calculated for a 1.5 µmsilica sphere, and agrees reasonably well with experimental results for the optical tweezer in Fig. 1(c), the optical conveyor in Fig. 1(d) and the optical conveyor in Fig. 3(a). The other set is calculated in the Rayleigh regime at ka p = 0.5, and agrees quantitatively with the dipole result from Eq. (0). Despite the dipole prediction from Eq. (13), the force experienced by a large particle in an optical conveyor can be purely repulsive. Lorenz-Mie calculations reveal that additional radiation pressure due to off-axis scattering can overwhelm trapping forces due to axial intensity gradients for particles larger than the wavelength of light. Although a conveyor can be projected with a range exceeding this limit, it will not be able to transport large particles in the retrograde direction, and so will not act as a tractor beam. The dependence of the trapping force on range for small particles shows no such crossover from stable trapping to repulsion. While the optical tweezer is inherently limited to R/A < 9 the optical conveyor extends indefinitely, albeit with a stiffness that falls off as R 4. The maximum trapping force in this regime falls off as R. Both the force and the stiffness scale with the intensity of the beam, and therefore with the laser power. These considerations demonstrate that optical conveyors are viable candidates for long-ranged tractor beams, particularly for objects that are smaller than the wavelength of light. CONCLUSIONS Interference endows optical conveyors with trapping characteristics surpassing those of conventional singlebeam optical traps. Most notably, optical conveyors are universal traps that can hold and transport small objects regardless of their light scattering properties. Under conditions where optical tweezers also are effective, moreover, optical conveyors make stiffer traps, particularly in the axial direction. Optical conveyors act as universal traps because their intensity vanishes at points along the optical axis [4 6]. This ensures that retrograde forces arising from intensity gradients can counteract radiation pressure to hold illuminated objects in place. These traps can be positioned any placed along the propagation-invariant range of the conveyor by varying the relative phase of the constituent Bessel beams. As a consequence, optical conveyors can transport objects over substantially larger axial distances than conventional optical traps [4 6], which are inherently limited by radiation pressure. Further improvements in performance almost certainly can be realized by appropriately structuring the intensity and phase gradients [38 40] along the projected beam. The range and trapping strength also should improve in optical conveyors created with radially polarized light [41 43]. These insights apply also to other interferometrically structured beams of light, such as solenoidal waves []. Both as tractor beams, and also as ordinary optical traps, interferometrically structured beams of light offer clear benefits for optical micromanipulation. ACKNOWLEDGMENTS This work was supported principally by the National Science Foundation through grant number DMR , and in part by a grant from NASA. [1] E. E. Smith, Spacehounds of IPC, Amazing Stories, July (1931). [] S.-H. Lee, Y. Roichman, and D. G. Grier, Optical solenoid beams, Opt. Express 18(7), (010). [3] O. Brzobohatý, V. Karásek, M. Šiler, L. Chvátal, T. Čižmár, and P. Zemánek, Experimental demonstration of optical transport, sorting and self-arrangement using a tractor beam, Nature Photonics 7(), (013). [4] D. B. Ruffner and D. G. Grier, Optical conveyors: a class of active tractor beams, Phys. Rev. Lett. 109(16), (01). [5] T. Čižmar, V. Garcés-Chávez, K. Dhokalia, and P. Zemánek, Optical conveyor belt for delivery of submicron objects, Appl. Phys. Lett. 86(17), (005). [6] T. Čižmar, V. Kollárová, Z. Bouchal, and P. Zemánek, Sub-micron particle organization by self-imaging of nondiffracting beams, New J. Phys. 8, 43 (006). [7] T. A. Nieminen, V. L. Y. Loke, A. B. Stilgoe, G. Knoner, A. M. Branczyk, N. R. Heckenberg, and H. Rubinsztein- Dunlop, Optical tweezers computational toolbox, J. Opt. A 9(8), S196 S03 (007). [8] A. B. Stilgoe, T. A. Nieminen, G. Knöener, N. R. Heckenberg, and H. Rubinsztein-Dunlop, The effect of Mie resonances on trapping in optical tweezers, Opt. Express 16(19), 15,039 15,051 (008). [9] B. Sun, Y. Roichman, and D. G. Grier, Theory of holographic optical trapping, Opt. Express 16(0), 15,765 15,776 (008). [10] J. Durnin, J. J. Miceli, Jr., and J. H. Eberly, Diffractionfree beams, Phys. Rev. Lett. 58(15), (1987). [11] Y. Roichman, I. Cholis, and D. G. Grier, Volumetric imaging of holographic optical traps, Opt. Express 14(), 10,907 10,91 (006). [1] D. G. Grier, A revolution in optical manipulation, Nature 44, (003). [13] S.-H. Lee, Y. Roichman, G.-R. Yi, S.-H. Kim, S.-M. Yang, A. van Blaaderen, P. van Oostrum, and D. G. Grier, Characterizing and tracking single colloidal par-

7 7 ticles with video holographic microscopy, Opt. Express 15(6), 18,75 18,8 (007). [14] F. C. Cheong, K. Xiao, D. J. Pine, and D. G. Grier, Holographic characterization of individual colloidal spheres porosities, Soft Matter 7(15), (011). [15] B. J. Krishnatreya, A. Colen-Landy, P. Hasebe, B. A. Bell, J. R. Jones, A. Sunda-Meya, and D. G. Grier, Measuring Boltzmann s constant through holographic video microscopy of a single sphere, Am. J. Phys. 8(1), 3 31 (014). [16] F. C. Cheong, B. J. Krishnatreya, and D. G. Grier, Strategies for three-dimensional particle tracking with holographic video microscopy, Opt. Express 18(13), 13,563 13,573 (010). [17] H. Moyses, B. J. Krishnatreya, and D. G. Grier, Robustness of holographic video microscopy against defects in illumination, Opt. Express 1(5), (013). [18] A. Rohrbach, Stiffness of optical traps: quantitative agreement between experiments and electromagnetic theory, Phys. Rev. Lett. 95(16), (005). [19] E. L. Florin, A. Pralle, E. H. K. Stelzer, and J. K. H. Horber, Photonic force microscope calibration by thermal noise analysis, Appl. Phys. A 66, S75 S78 (1998). [0] Z. Bouchal and M. Olivk, Non-diffractive vector Bessel beams, J. Mod. Opt. 4(8), (1995). [1] M. Polin, K. Ladavac, S.-H. Lee, Y. Roichman, and D. G. Grier, Optimized holographic optical traps, Opt. Express 13(15), (005). [] J. M. Taylor and G. D. Love, Multipole expansion of Bessel and Gaussian beams for Mie scattering calculations, J. Opt. Soc. Am. A 6(), 78 8 (009). [3] J. Chen, J. Ng, Z. Lin, and C. T. Chan, Optical pulling force, Nature Photonics 5(9), (011). [4] M. Šiler and P. Zemánek, Optical forces in a nondiffracting vortex beam, J. Quant. Spectr. Rad. Trans. 16, (013). [5] N. Wang, J. Chen, S. Liu, and Z. Lin, Dynamical and phase-diagram study on stable optical pulling force in Bessel beams, Phys. Rev. A 87(6), (013). [6] G. Gouesbet, T-matrix formulation and generalized Lorenz-Mie theories in spherical coordinates, Opt. Comm. 83(4), (010). [7] C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley Interscience, 1987). [8] G. Gouesbet and G. Gréhan, Generalized Lorenz-Mie Theories (Springer-Verlag, 011). [9] J. P. Barton, Internal and near-surface electromagnetic fields for a spheroidal particle with arbitrary illumination, Appl. Opt. 34(4), (1995). [30] J. P. Barton, D. R. Alexander, and S. A. Schaub, Theoretical determination of net radiation force and torque for a spherical particle illuminated by a focused laser beam, J. Appl. Phys. 66(10), (1989). [31] J. P. Gordon, Radiation forces and momenta in dielectric media, Phys. Rev. A 8(1), 14 1 (1973). [3] P. C. Chaumet and M. Nieto-Vesperinas, Time-averaged total force on a dipolar sphere in an electromagnetic field, Opt. Lett. 5(15), (000). [33] B. T. Draine, The discrete-dipole approximation and its application to interstellar graphite grains, Astrophys. J. 333(), (1988). [34] S. Sukhov and A. Dogariu, Negative nonconservative forces: optical tractor beams for arbitrary objects, Phys. Rev. Lett. 107(0), 0360 (011). [35] A. Novitsky, C.-W. Qiu, and A. Lavrinenko, Materialindependent and size-independent tractor beams for dipole objects, Phys. Rev. Lett. 109(), 0390 (01). [36] A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, Observation of a single-beam gradient force optical trap for dielectric particles, Opt. Lett. 11(5), (1986). [37] K. Harada, O. Kamimura, H. Kasai, T. Matsuda, A. Tonomura, and V. V. Moshchalkov, Direct observation of vortex dynamics in superconducting films with regular arrays of defects, Science 74(590), (1996). [38] Y. Roichman, B. Sun, Y. Roichman, J. Amato-Grill, and D. G. Grier, Optical forces arising from phase gradients, Phys. Rev. Lett. 100(1), (008). [39] T. Čižmár and K. Dholakia, Tunable Bessel light modes: engineering the axial propagation, Opt. Express 17(18), 15,558 15,570 (009). [40] M. Lei, Z. Li, S. Yan, B. Yao, D. Dan, Y. Qi, J. Qian, Y. Yang, P. Gao, and T. Ye, Long-distance axial trapping with focused annular laser beams, PLoS ONE 8(3), e57984 (013). [41] H. Kawauchi, K. Yonezawa, Y. Kozawa, and S. Sato, Calculation of optical trapping forces on a dielectric sphere in the ray optics regime produced by a radially polarized beam, Opt. Lett. 3(13), (007). [4] M. Michihata, T. Hayashi, and Y. Takaya, Measurement of axial and transverse trapping stiffness of optical tweezers in air using a radially polarized beam, Appl. Opt. 48(3), (009). [43] L. Carretero, P. Acebal, and S. Blaya, Threedimensional analysis of optical forces generated by an active tractor beam using radial polarization, Opt. Express (3), (014).