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.

1 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 [2], interfering plane waves [3], and optical conveyors created with coaxial Bessel beams [4, 5, 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, 8, 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 ω\omega, linearly polarized along x^\hat{x} and propagating along z^\hat{z} through a medium of refractive index nmn_{m}, is the superposition of two Bessel beams [5, 6, 4],

𝐄(𝐫,t)=12E0e-iωt[𝐛1(k𝐫)+eiφ(t)𝐛2(k𝐫)],\mathbf{E}(\mathbf{r},t)=\frac{1}{2}\,E_{0}e^{-i\omega t}\,\left[\mathbf{b}_{1% }(k\mathbf{r})+e^{i\varphi(t)}\mathbf{b}_{2}(k\mathbf{r})\right], (1)

each of which may be described as a conical superposition of plane waves [6],

𝐛j(k𝐫)=02πϵ^(θj,ϕ)ei𝐤(θj,ϕ)𝐫dϕ,\mathbf{b}_{j}(k\mathbf{r})=\int_{0}^{2\pi}\hat{\epsilon}(\theta_{j},\phi)\,e^% {i\mathbf{k}(\theta_{j},\phi)\cdot\mathbf{r}}\,d\phi, (2a)
ϵ^(θ,ϕ)=cosϕθ^+sinϕϕ^\hat{\epsilon}(\theta,\phi)=\cos\phi\,\hat{\theta}+\sin\phi\,\hat{\phi} (2b)
is the polarization of a plane wave incident on the optical axis at polar angle θ\theta and azimuthal angle ϕ\phi, and where
𝐤(θ,ϕ)=k(sinθcosϕx^+sinθsinϕy^+cosθz^)\mathbf{k}(\theta,\phi)=k(\sin\theta\cos\phi\,\hat{x}+\sin\theta\sin\phi\,\hat% {y}+\cos\theta\,\hat{z}) (2c)

is the corresponding wave vector. In the paraxial approximation, which is appropriate for long-range tractor beams, Eq. (2) reduces to

𝐛j(k𝐫)J0(1-ηj2kr)eiηjkzx^,\mathbf{b}_{j}(k\mathbf{r})\approx J_{0}\!\left(\sqrt{1-\eta_{j}^{2}}kr\right)% e^{i\eta_{j}kz}\,\hat{x}, (3)

where ηj\eta_{j} is related to the Bessel beam’s cone angle by ηj=cosθj\eta_{j}=\cos\theta_{j}. The approximate expression in Eq. (3) differs from the exact expression in Eq. (2) by terms involving higher-order Bessel functions[20, 23], which vanish on the optical axis. The numerical results developed in Sec. 2.1 are based on Eq. (2). Analytical results developed in Sec. 2.2 are more readily obtained from the approximation in Eq. (3).

The two Bessel beams comprising an optical conveyor share the same amplitude E0E_{0}, frequency ω\omega and polarization along x^\hat{x}, but differ in their relative phase, φ(t)\varphi(t), and also in their axial wave numbers, ηjk\eta_{j}k, which are reduced from the plane wave value, k=nmω/ck=n_{m}\omega/c, by the dimensionless factor ηj(0,1]\eta_{j}\in(0,1]. Here, cc is the speed of light in vacuum. The superposition is usefully characterized by the mean convergence factor η=(η1+η2)/2\eta=(\eta_{1}+\eta_{2})/2 and the difference Δη=|η1-η2|\Delta\eta=\left|\eta_{1}-\eta_{2}\right|. The upper limit of an optical conveyor’s range, RR, is set by the non-diffracting range of the most strongly converging Bessel beam [10],

RAcotθ1=Aη11-η12R\leq A\cot\theta_{1}=A\frac{\eta_{1}}{\sqrt{1-\eta_{1}^{2}}} (4)

where AA is the radius of the beam’s aperture. Longer ranges can be achieved with values of η\eta approaching 1.

Figure 1: (a) Experimental reconstruction of an optical conveyor η=0.8\eta=0.8, Δη=0.04\Delta\eta=0.04, and (b) η=0.8\eta=0.8, Δη=0.086\Delta\eta=0.086. (c) Measured trajectory of a 1.5 μm\mu\text{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 Δη\Delta\eta. (f) Stiffness ratio, kz/krk_{z}/k_{r}, as a function of Δη\Delta\eta. 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.

Figures 1(a) and 1(b) show the measured [11] three-dimensional intensity distribution of optical conveyors with equal values of η\eta and differing values of Δη\Delta\eta that were projected with the holographic optical trapping technique [12, 4]. These beams were powered by a linearly polarized diode-pumped solid-state laser (Coherent Verdi) operating at a vacuum wavelength of λ=532nm\lambda=532~{}\text{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×100\times, 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 ηjk\eta_{j}k [4]. Equation (1) therefore should be considered an idealized model for the actual beam.

The intensity distribution, I(𝐫,t)=12nmϵ0c|𝐄(𝐫,t)|2I(\mathbf{r},t)=\frac{1}{2}n_{m}\epsilon_{0}c\left|\mathbf{E}(\mathbf{r},t)% \right|^{2}, has maxima at axial positions

zn(t)=2πn+φ(t)Δηk,z_{n}(t)=\frac{2\pi n+\varphi(t)}{\Delta\eta\,k}, (5)

each of which can act as an optical trap for a small object. Here, ϵ0\epsilon_{0} is the permittivity of space. Varying φ(t)\varphi(t) as a function of time moves these extrema, and thus conveys trapped objects along the beam. Increasing Δη\Delta\eta reduces the spacing between maxima, as shown in Figs. 1(a) and 1(b), and so provides control over the optical force profile.

2 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 79mW/μm279~{}\text{mW}/\mu\text{m}^{2}, as measured with imaging photometry. Holographic characterization [13, 14, 15] reveals the sphere’s radius to be ap=0.730±0.005μma_{p}=0.730\pm 0.005~{}\mu\text{m}, and its refractive index to be np=1.424±0.005n_{p}=1.424\pm 0.005. Holographic tracking [13, 16] yields the sphere’s position with 1 nm precision in-plane and 3 nm resolution axially [17, 15] 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 tweezer, resulting in a nearly isotropic trajectory. This suggests that an optical conveyor makes a substantially stiffer trap, even though the 2.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(𝐫)=12krr2+12kzz2,U(\mathbf{r})=\frac{1}{2}k_{r}r^{2}+\frac{1}{2}k_{z}z^{2}, (6)

we may estimate the transverse and axial trap stiffness, krk_{r} and kzk_{z} from the particle’s trajectory using thermal fluctuation analysis [19, 21]. These measurements are performed for fixed values of the relative phase, φ(t)\varphi(t), so that the optical traps do not move during the measurement. For the trajectory in Fig. 1(c), we obtain kr=2.2±0.1pN/μmk_{r}=2.2\pm 0.1~{}\text{pN}/\mu\text{m} and kz=0.89±0.03pN/μmk_{z}=0.89\pm 0.03~{}\text{pN}/\mu\text{m}, and an anisotropy of kr/kz=1.57±0.04\sqrt{k_{r}/k_{z}}=1.57\pm 0.04. The equivalent results for the optical tweezer are kr=2.8±0.4pN/μmk_{r}=2.8\pm 0.4~{}\text{pN}/\mu\text{m}, kz=0.41±0.08pN/μmk_{z}=0.41\pm 0.08~{}\text{pN}/\mu\text{m} and kr/kz=2.6±0.4\sqrt{k_{r}/k_{z}}=2.6\pm 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 Δη\Delta\eta for fixed η\eta, 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\Delta\eta>0.06, which corresponds to an axial period less than 6.7 μm\mu\text{m}.

2.1 Generalized Lorenz-Mie theory

Optical conveyors’ superior trapping performance is consistent with predictions of generalized Lorenz-Mie theory [22, 23, 6, 24, 25], which are plotted as continuous curves in Figs. 1(e) and 1(f). For these calculations, each Bessel beam is expanded as a series,

𝐛j(𝐫)=n=1m=-nn[amn(θj)𝐌nm(1)(k𝐫)+bmn(θj)𝐍nm(1)(k𝐫)],\mathbf{b}_{j}(\mathbf{r})=\sum_{n=1}^{\infty}\sum_{m=-n}^{n}\left[a_{mn}(% \theta_{j})\,\mathbf{M}_{nm}^{(1)}(k\mathbf{r})+b_{mn}(\theta_{j})\,\mathbf{N}% _{nm}^{(1)}(k\mathbf{r})\right], (7)

in the vector spherical harmonics, 𝐌nm(1)(k𝐫)\mathbf{M}_{nm}^{(1)}(k\mathbf{r}) and 𝐍nm(1)(k𝐫)\mathbf{N}_{nm}^{(1)}(k\mathbf{r}), that constitute the natural basis for transverse electric (TE) and transverse magnetic (TM) waves, respectively [26]. The expansion coefficients, amn(θj)a_{mn}(\theta_{j}) and bmn(θj)b_{mn}(\theta_{j}), have been reported previously for individual Bessel beams [22]; they can be superimposed according to Eq. (1) to create an optical conveyor with specified values of η\eta, Δη\Delta\eta and relative phase φ(t)\varphi(t).

The associated expansion of the scattered field,

𝐄s(𝐫,t)=E0e-iωtn=1m=-nn{[rmn(θ1)+eiφ(t)rmn(θ2)]𝐌nm(3)(k𝐫)+[smn(θ1)+eiφ(t)smn(θ2)]𝐍nm(3)(k𝐫)},\mathbf{E}_{s}(\mathbf{r},t)=E_{0}\,e^{-i\omega t}\,\sum_{n=1}^{\infty}\sum_{m% =-n}^{n}\left\{\left[r_{mn}(\theta_{1})+e^{i\varphi(t)}r_{mn}(\theta_{2})% \right]\,\mathbf{M}_{nm}^{(3)}(k\mathbf{r})+\left[s_{mn}(\theta_{1})+e^{i% \varphi(t)}s_{mn}(\theta_{2})\right]\,\mathbf{N}_{nm}^{(3)}(k\mathbf{r})\right\}, (8a)
is obtained by matching boundary conditions for the electric and magnetic fields at the sphere’s surface [23, 24, 25, 27, 28, 29]. The expansion coefficients,
rmn(θj)\displaystyle r_{mn}(\theta_{j}) =-anamn(θj)and\displaystyle=-a_{n}\,a_{mn}(\theta_{j})\quad\text{and} (8b)
smn(θj)\displaystyle s_{mn}(\theta_{j}) =-bnbmn(θj),\displaystyle=-b_{n}\,b_{mn}(\theta_{j}), (8c)

are related to the expansion coefficients of the incident field by the standard Mie coefficients, ana_{n} and bnb_{n}, [27], which depend on the radius and refractive index of the sphere.

The superposed incident and scattered fields, 𝐄(𝐫,t)+𝐄s(𝐫,t)\mathbf{E}(\mathbf{r},t)+\mathbf{E}_{s}(\mathbf{r},t), contribute to the Maxwell stress tensor, 𝐓(𝐫,t)\mathbf{T}(\mathbf{r},t), whose integral over a closed surface provides an estimate for the optically-induced force on a sphere centered at 𝐫\mathbf{r}:

𝐅(𝐫,t)=Sn^𝐓(𝐫,t)d𝐫,\mathbf{F}(\mathbf{r},t)=\oint_{S}\hat{n}\cdot\mathbf{T}(\mathbf{r}^{\prime},t% )\,d\mathbf{r}^{\prime}, (9)

where n^\hat{n} is the unit normal to the surface SS enclosing the sphere. In practice, the integral is computed directly from the combined expansion coefficients using established techniques [7, 8, 9, 30]. This calculation must also be repeated for each value of the relative phase φ(t)\varphi(t).

A trap’s position, 𝐫0(t)\mathbf{r}_{0}(t), for a given value of φ(t)\varphi(t) is calculated numerically as a solution of 𝐅(𝐫0(t),t)=0\mathbf{F}(\mathbf{r}_{0}(t),t)=0. The traps’s effective stiffness along r^ν\hat{r}_{\nu} then is calculated as

kν=-νFν(𝐫,t)|𝐫=𝐫0(t).k_{\nu}=-\left.\partial_{\nu}F_{\nu}(\mathbf{r},t)\right|_{\mathbf{r}=\mathbf{% r}_{0}(t)}. (10)

Results for the predicted transverse and axial stiffness are plotted in Fig. 1(e) for conveyor beams with η=0.8\eta=0.8 and Δη\Delta\eta 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 41mW/μm241~{}\text{mW}/\mu\text{m}^{2}, 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.

2.2 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, 32],

𝐅(𝐫,t)=12{αeν=13Eν(𝐫,t)Eν(𝐫,t)},\mathbf{F}(\mathbf{r},t)=\frac{1}{2}\,\Re\left\{\alpha_{e}\,\sum_{\nu=1}^{3}E_% {\nu}(\mathbf{r},t)\nabla E_{\nu}^{\ast}(\mathbf{r},t)\right\}, (11)

is proportional to the object’s polarizability αe=αe+iαe′′\alpha_{e}=\alpha_{e}^{\prime}+i\alpha_{e}^{\prime\prime}. The polarizability of a dielectric sphere is related to its size and refractive index by the Clausius-Mossotti-Draine relation [33]

αe=4πϵ0nm2Kap31-i23Kk3ap3,\alpha_{e}=\frac{4\pi\epsilon_{0}n_{m}^{2}K\,a_{p}^{3}}{1-i\frac{2}{3}Kk^{3}a_% {p}^{3}}, (12)

where K=(np2-nm2)/(np2+2nm2)K=(n_{p}^{2}-n_{m}^{2})/(n_{p}^{2}+2n_{m}^{2}) is the Lorentz-Lorenz factor. Absorptivity increases the imaginary part of αe\alpha_{e}. Conductivity contributes an imaginary part to KK. 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

Fz(𝐫,t)E02-14αeΔηsin(Φ(z,t))+αe′′ηcos2(12Φ(z,t))\frac{F_{z}(\mathbf{r},t)}{E_{0}^{2}}\approx-\frac{1}{4}\alpha_{e}^{\prime}% \Delta\eta\,\sin\bigl(\Phi(z,t)\bigr)+\alpha_{e}^{\prime\prime}\,\eta\,\cos^{2% }\!\left(\frac{1}{2}\Phi(z,t)\right) (13)

and a transverse force

Fr(𝐫,t)E02αekr2(1-η2-14Δη2)cos2(12Φ(z,t))+14αe′′ηΔηkrsin(Φ(z,t)),\frac{F_{r}(\mathbf{r},t)}{E_{0}^{2}}\approx\alpha_{e}^{\prime}\,\frac{kr}{2}% \left(1-\eta^{2}-\frac{1}{4}\Delta\eta^{2}\right)\cos^{2}\!\left(\frac{1}{2}% \Phi(z,t)\right)+\frac{1}{4}\alpha_{e}^{\prime\prime}\,\eta\,\Delta\eta\,kr% \sin\bigl(\Phi(z,t)\bigr), (14)

where Φ(z,t)=Δηkz-φ(t)\Phi(z,t)=\Delta\eta\,kz-\varphi(t), and where we have assumed kr1kr\ll 1. These results are obtained by averaging over times long compared with the optical cycle but short compared with variations in φ(t)\varphi(t). The particle is trapped where the force vanishes, which occurs at axial positions Zn(t)Z_{n}(t) that are displaced from the intensity maxima by an amount that depends on the particle’s light-scattering properties,

Zn(t)-zn(t)=2Δηktan-1(αe′′αe2ηΔη).Z_{n}(t)-z_{n}(t)=\frac{2}{\Delta\eta k}\,\tan^{-1}\!\left(\frac{\alpha_{e}^{% \prime\prime}}{\alpha_{e}^{\prime}}\frac{2\eta}{\Delta\eta}\right). (15)

Equations (10) and (16) then yield the traps’ axial stiffness,

kz=14|αe|kE02Δη2,k_{z}=\frac{1}{4}\left|\alpha_{e}^{\prime}\right|kE_{0}^{2}\,\Delta\eta^{2}, (16)

which is strictly positive. An optical conveyor thus can trap and transport any dipolar particle, regardless of its light-scattering characteristics. This universal, material-independent trapping capability does not require feedback [34] or fine tuning of the beam’s properties [35]. Decreasing the inter-trap separation by increasing Δη\Delta\eta enhances intensity gradients and thus increases trap stiffness.

Picking the largest possible value of Δη\Delta\eta to optimize axial trapping is not necessarily the best strategy. An optical conveyor’s transverse stiffness,

kr=kzαe2(1-η2-14Δη2)-2αe′′2η212αe2Δη2+2αe′′2η2,k_{r}=k_{z}\frac{{\alpha_{e}^{\prime}}^{2}\left(1-\eta^{2}-\frac{1}{4}\Delta% \eta^{2}\right)-2{\alpha_{e}^{\prime\prime}}^{2}\eta^{2}}{\frac{1}{2}{\alpha_{% e}^{\prime}}^{2}\Delta\eta^{2}+2{\alpha_{e}^{\prime\prime}}^{2}\eta^{2}}, (17)

can vanish or even change sign as Δη\Delta\eta increases.

Figure 2: (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).

The dependence of trap stiffness on particle size predicted by Eqs. (12), (16) and (17) is plotted as dashed curves in Fig. 2 for silica spheres in an optical conveyor with η=0.8\eta=0.8 and Δη=0.086\Delta\eta=0.086. Results from the dipole approximation agree well with the Lorenz-Mie results that are plotted as solid curves in Fig. 2, at least for kap1ka_{p}\leq 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 kap>1ka_{p}>1, and suggests that 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>0k_{r}>0) according to Eq. (17) establishes an upper bound on Δη\Delta\eta and, through Eq. (4), a lower bound on the range over which an optical conveyor can transport small objects:

RA2(αeαe′′-αe′′αe)<R.R_{\perp}\leq\frac{A}{2}\left(\frac{\alpha_{e}^{\prime}}{\alpha_{e}^{\prime% \prime}}-\frac{\alpha_{e}^{\prime\prime}}{\alpha_{e}^{\prime}}\right)<R. (18)

The practically accessible range of transport may be substantially greater than RR_{\perp} for particles larger than the wavelength of light.

2.3 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]

𝐄G(z,t)=E0zRz2+zR2eikzeiζ(z)e-iωtϵ^,\mathbf{E}_{G}(z,t)=E_{0}\frac{z_{R}}{\sqrt{z^{2}+z_{R}^{2}}}e^{ikz}e^{i\zeta(% z)}e^{-i\omega t}\,\hat{\epsilon}, (19)

where zR2R2/(kA2)z_{R}\approx 2R^{2}/(kA^{2}) is the Rayleigh range of a Gaussian beam converging at distance RR from an aperture of radius AA, and where ζ(z)=tan-1(z/zR)\zeta(z)=\tan^{-1}(z/z_{R}) is the Gouy phase. Because we are interested in long-ranged axial transport, we assume R>AR>A. The axial component of the associated force,

FG(z)E02=-12zR2zαe-k(z2-zR2)αe′′(z2+zR2)2,\frac{F_{G}(z)}{E_{0}^{2}}=-\frac{1}{2}z_{R}^{2}\,\frac{z\alpha_{e}^{\prime}-k% \left(z^{2}-z_{R}^{2}\right)\alpha_{e}^{\prime\prime}}{\left(z^{2}+z_{R}^{2}% \right)^{2}}, (20)

forms a trap only for particles satisfying

(αeαe′′)2>4kzR(kzR-1).\left(\frac{\alpha_{e}^{\prime}}{\alpha_{e}^{\prime\prime}}\right)^{2}>4kz_{R}% (kz_{R}-1). (21)

This condition is most easily satisfied in strongly converging beams for which zRz_{R} is small. The condition on αe\alpha_{e} is qualitatively consistent with numerical studies [7, 8, 9] which describe the difficulty of trapping high index or absorbing particles with laser tweezers.

3 Transport range

Figure 3: (a) Experimental reconstruction of an optical conveyor with η=0.96\eta=0.96 and Δη=0.04\Delta\eta=0.04. (b) Trajectory of a 1.5 μm\mu\text{m}silica sphere transported by the same optical conveyor over a range of 66 μm\mu\text{m}. (c) Axial stiffness as a function of transport range. Solid curves show Lorenz-Mie predictions for optical conveyors transporting large (1.5 μm\mu\text{m}-diameter) and small (0.036 μm\mu\text{m}-diameter, kap=0.5ka_{p}=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).

The ultimate rate of an optical conveyor is limited by the non-diffracting range of the constituent Bessel beams. Equation (21), by contrast, establishes an inherent upper limit on an optical tweezer’s range

RG=A21+1+(αeαe′′)2.R_{G}=\frac{A}{2}\,\sqrt{1+\sqrt{1+\left(\frac{\alpha_{e}^{\prime}}{\alpha_{e}% ^{\prime\prime}}\right)^{2}}}. (22)

At ranges RR beyond RGR_{G}, 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\eta=0.96 and Δη=0.04\Delta\eta=0.04 that was projected with a 60×60\times objective lens (Nikon Plan-Apo, numerical aperture 1.4 oil immersion). Holographic projection limits the conveyor’s range to R=70μmR=70~{}\mu\text{m}, and thus suggests an effective numerical aperture of 0.43 given the convergence angle of θ=19\theta=19^{\circ}. The holographically measured trajectory plotted in Fig. 3(b) shows this optical conveyor trapping and transporting a 1.5 μm\mu\text{m}-diameter silica sphere over 66 μm\mu\text{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 numerical 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 41mW/μm241~{}\text{mW}/\mu\text{m}^{2}. One set of curves is calculated for a 1.5 μm\mu\text{m}silica 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 kap=0.5ka_{p}=0.5, and agrees quantitatively with the dipole result from Eq. (20).

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<9R/A<9 the optical conveyor extends indefinitely, albeit with a stiffness that falls off as R-4R^{-4}. The maximum trapping force in this regime falls off as R-2R^{-2}. 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.

4 Conclusions

Interference endows optical conveyors with trapping characteristics surpassing those of conventional single-beam 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 [5, 6, 4]. 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 [5, 6, 4], 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, 39, 40] along the projected beam. The range and trapping strength also should improve in optical conveyors created with radially polarized light [41, 42, 43]. These insights apply also to other interferometrically structured beams of light, such as solenoidal waves [2]. Both as tractor beams, and also as ordinary optical traps, interferometrically structured beams of light offer clear benefits for optical micromanipulation.


This work was supported principally by the National Science Foundation through grant number DMR-1305875, and in part by a grant from NASA.


  • 1 E. E. Smith, “Spacehounds of IPC,” Amazing Stories, July (1931).
  • 2 S.-H. Lee, Y. Roichman, and D. G. Grier, “Optical solenoid beams,” Opt. Express 18(7), 6988–6993 (2010).
  • 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(2), 123–127 (2013).
  • 4 D. B. Ruffner and D. G. Grier, “Optical conveyors: a class of active tractor beams,” Phys. Rev. Lett. 109(16), 163903 (2012).
  • 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), 174101 (2005).
  • 6 T. Čižmar, V. Kollárová, Z. Bouchal, and P. Zemánek, ‘‘Sub-micron particle organization by self-imaging of non-diffracting beams,” New J. Phys. 8, 43 (2006).
  • 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–S203 (2007).
  • 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 (2008).
  • 9 B. Sun, Y. Roichman, and D. G. Grier, “Theory of holographic optical trapping,” Opt. Express 16(20), 15,765–15,776 (2008).
  • 10 J. Durnin, J. J. Miceli, Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987).
  • 11 Y. Roichman, I. Cholis, and D. G. Grier, “Volumetric imaging of holographic optical traps,” Opt. Express 14(22), 10,907–10,912 (2006).
  • 12 D. G. Grier, “A revolution in optical manipulation,” Nature 424, 810–816 (2003).
  • 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 particles with video holographic microscopy,” Opt. Express 15(26), 18,275–18,282 (2007).
  • 14 F. C. Cheong, K. Xiao, D. J. Pine, and D. G. Grier, “Holographic characterization of individual colloidal spheres’ porosities,” Soft Matter 7(15), 6816–6819 (2011).
  • 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. 82(1), 23–31 (2014).
  • 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 (2010).
  • 17 H. Moyses, B. J. Krishnatreya, and D. G. Grier, “Robustness of holographic video microscopy against defects in illumination,” Opt. Express 21(5), 5968–5973 (2013).
  • 18 A. Rohrbach, “Stiffness of optical traps: quantitative agreement between experiments and electromagnetic theory,” Phys. Rev. Lett. 95(16), 168102 (2005).
  • 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).
  • 20 Z. Bouchal and M. Olivík, “Non-diffractive vector Bessel beams,” J. Mod. Opt. 42(8), 1555-1566 (1995).
  • 21 M. Polin, K. Ladavac, S.-H. Lee, Y. Roichman, and D. G. Grier, “Optimized holographic optical traps,” Opt. Express 13(15), 5831–5845 (2005).
  • 22 J. M. Taylor and G. D. Love, “Multipole expansion of Bessel and Gaussian beams for Mie scattering calculations,” J. Opt. Soc. Am. A 26(2), 278–282 (2009).
  • 23 J. Chen, J. Ng, Z. Lin, and C. T. Chan, “Optical pulling force,” Nature Photonics 5(9), 531–534 (2011).
  • 24 M. Šiler and P. Zemánek, “Optical forces in a non-diffracting vortex beam,” J. Quant. Spectr. Rad. Trans. 126, 78–83 (2013).
  • 25 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), 063812 (2013).
  • 26 G. Gouesbet, “T-matrix formulation and generalized Lorenz-Mie theories in spherical coordinates,” Opt. Comm. 283(4), 517–521 (2010).
  • 27 C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley Interscience, 1987).
  • 28 G. Gouesbet and G. Gréhan, Generalized Lorenz-Mie Theories (Springer-Verlag, 2011).
  • 29 J. P. Barton, “Internal and near-surface electromagnetic fields for a spheroidal particle with arbitrary illumination,” Appl. Opt. 34(24), 5542–5551 (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), 4594–4602 (1989).
  • 31 J. P. Gordon, “Radiation forces and momenta in dielectric media,” Phys. Rev. A 8(1), 14–21 (1973).
  • 32 P. C. Chaumet and M. Nieto-Vesperinas, “Time-averaged total force on a dipolar sphere in an electromagnetic field,” Opt. Lett. 25(15), 1065–1067 (2000).
  • 33 B. T. Draine, “The discrete-dipole approximation and its application to interstellar graphite grains,” Astrophys. J. 333(2), 848–872 (1988).
  • 34 S. Sukhov and A. Dogariu, “Negative nonconservative forces: optical “tractor beams” for arbitrary objects,” Phys. Rev. Lett. 107(20), 203602 (2011).
  • 35 A. Novitsky, C.-W. Qiu, and A. Lavrinenko, “Material-independent and size-independent tractor beams for dipole objects,” Phys. Rev. Lett. 109(2), 023902 (2012).
  • 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), 288–290 (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 274(5290), 1167–1170 (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), 013602 (2008).
  • 39 T. Čižmár and K. Dholakia, “Tunable Bessel light modes: engineering the axial propagation,” Opt. Express 17(18), 15,558–15,570 (2009).
  • 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 (2013).
  • 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. 32(13), 1839–1841 (2007).
  • 42 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(32), 6143–6151 (2009).
  • 43 L. Carretero, P. Acebal, and S. Blaya, “Three-dimensional analysis of optical forces generated by an active tractor beam using radial polarization,” Opt. Express 22(3), 3284–3295 (2014).