# Tractor Beams in the Rayleigh Limit

Aaron Yevick    David B. Ruffner    David G. Grier Department of Physics and Center for Soft Matter Research, New York University, New York, NY 10003
###### Abstract

A tractor beam is a travelling wave that transports illuminated objects back to its source, opposite to the wave’s direction of propagation, along its entire length. The requisite retrograde force arises when an object scatters the wave’s momentum density downstream into the direction of propagation, and then recoils upstream by conservation of momentum. Achieving this condition imposes constraints on the structure of the wave, which we elucidate in the Rayleigh limit, when the wavelength exceeds the size of the object. Continuously propagation-invariant modes such Bessel beams do not satisfy these conditions at dipole order in the multipole expansion, and so cannot serve as general-purpose long-ranged tractor beams. Modes with discrete propagation invariance, however, can act as first-order tractor beams. We demonstrate this by introducing a class of minimal solenoidal waves together with a set of design criteria that distinguish tractor beams that pull objects from repulsor beams that push them.

## I Introduction

Radiation pressure commonly is understood to push illuminated objects downstream along the direction of light’s propagation. The recent realization [1, 2, 3, 4] that appropriately structured traveling waves can pull objects upstream consequently has excited considerable interest. Retrograde forces arise when an illuminated object scatters more of the wave’s momentum into the downstream direction than otherwise would have been present [5, 6, 7]. The object then recoils upstream to conserve momentum. This happens, for example, when a small particle enters the diverging wave downstream of the focus in an optical trap, redirects light into the forward direction, and then recoils toward the focus [5, 6, 7]. The trapping beam’s divergence, however, inherently limits how far it can pull an object upstream.

Achieving long-range retrograde transport requires a non-diffracting or weakly-diffracting mode. A plane wave cannot serve in this capacity because all of its momentum is directed along the axis of propagation. Remarkably enough, some other non-diffracting traveling waves [1, 2, 3, 4, 8] have the requisite characteristics. Such modes have come to be known as tractor beams because of their resemblance to the popular science fiction trope [9].

Previously identified tractor-beam modes were discovered empirically [1, 2, 8]. Few have been demonstrated experimentally [2, 8, 10], and most of these demonstrations have relied on fine tuning of material properties [8, 10]. Here, we present an analytical theory for optical tractor beams that applies in the simplifying, yet practically important limit that the transported object is smaller than the wavelength of light, the so-called Rayleigh regime. This theoretical framework expresses the forces experienced by an illuminated object in terms of the experimentally controllable amplitude and phase profiles of a vector mode of light. It reveals that the most commonly discussed tractor-beam mode [1, 3, 4, 11] is inherently undermined by its continuous propagation invariance; such beams can pull only under exceptional circumstances, and then only weakly. These drawbacks can be avoided by creating tractor beams from modes that display discrete propagation invariance. We illustrate this by introducing a family of minimal solenoidal tractor-beam modes that exert retrograde forces as first-order photokinetic effects.

## II Photokinetic Forces

The electric and magnetic fields in a monochromatic beam of light at frequency $\omega$ may be projected into Cartesian components,

 $\displaystyle\boldsymbol{E}(\boldsymbol{r},t)$ $\displaystyle=\sum_{j=1}^{3}E_{j}(\boldsymbol{r})\,e^{-i\omega t}\,\hat{% \epsilon}_{j}\quad\text{and}$ (1a) $\displaystyle\boldsymbol{B}(\boldsymbol{r},t)$ $\displaystyle=-\frac{i}{\omega}\,\nabla\times\boldsymbol{E}(\boldsymbol{r},t)$ (1b) $\displaystyle=\sum_{j=1}^{3}B_{j}(\boldsymbol{r})\,e^{-i\omega t}\,\hat{% \epsilon}_{j},$ (1c)

where $E_{j}(\boldsymbol{r})$ and $B_{j}(\boldsymbol{r})$ are the complex-valued electric and magnetic field amplitudes along the Cartesian coordinate $\hat{\epsilon}_{j}$. For comparison with experimental implementations, it is convenient to express these field components in terms of real-valued amplitudes and phases,

 $E_{j}(\boldsymbol{r})=u_{j}(\boldsymbol{r})\,e^{i\phi_{j}(\boldsymbol{r})},$ (2)

and similarly for $B_{j}(\boldsymbol{r})$. Expressing optical forces in terms of $u_{j}(\boldsymbol{r})$ and $\phi_{j}(\boldsymbol{r})$ yields insights that are useful for designing and optimizing beams of light for optical micromanipulation, including tractor beams.

### II.1 First-order photokinetic forces

Photokinetic forces arise from light’s ability to polarize small objects. For simplicity and clarity, we will confine our discussion to optically isotropic objects whose dimensions are small enough relative to the wavelength of light that they may be treated as point-like scatterers. In this so-called Rayleigh regime, the light’s leading-order effect is to induce electric and magnetic dipole moments,

 $\displaystyle\boldsymbol{p}(\boldsymbol{r},t)$ $\displaystyle=\alpha_{e}\boldsymbol{E}(\boldsymbol{r},t)\quad\text{and}$ (3a) $\displaystyle\boldsymbol{m}(\boldsymbol{r},t)$ $\displaystyle=\alpha_{m}\boldsymbol{B}(\boldsymbol{r},t),$ (3b)

respectively, where $\alpha_{e}$ and $\alpha_{m}$ are the object’s electric and magnetic dipole polarizabilities.

The induced dipole moments described by Eq. (3) experience time-averaged forces in gradients of the electric [12] and magnetic [13] fields,

 $\displaystyle\boldsymbol{F}_{e}(\boldsymbol{r})$ $\displaystyle=\frac{1}{2}\,\Re\left\{\sum_{j=1}^{3}p_{j}(\boldsymbol{r},t)% \nabla E_{j}^{\ast}(\boldsymbol{r},t)\right\}\quad\text{and}$ (4a) $\displaystyle\boldsymbol{F}_{m}(\boldsymbol{r})$ $\displaystyle=\frac{1}{2}\,\Re\left\{\sum_{j=1}^{3}m_{j}(\boldsymbol{r},t)% \nabla B_{j}^{\ast}(\boldsymbol{r},t)\right\}.$ (4b)

Although Eq. (4) completely describes optical dipole forces, how to control these forces is more clearly revealed when they are presented in terms of the amplitudes and the phase of the fields. Expressed in these terms, the electric dipole force has the form [14, 15]

 $\boldsymbol{F}_{e}(\boldsymbol{r})=\frac{1}{4}\,\alpha_{e}^{\prime}\nabla\sum_% {j=1}^{3}u_{j}^{2}(\boldsymbol{r})+\frac{1}{2}\,\alpha_{e}^{\prime\prime}\sum_% {j=1}^{3}u_{j}^{2}(\boldsymbol{r})\nabla\phi_{j}(\boldsymbol{r}),$ (5)

where $\alpha_{e}^{\prime}$ and $\alpha_{e}^{\prime\prime}$ are the real and imaginary parts of $\alpha_{e}$, respectively.

Equation (5) and the analogous expression for the magnetic dipole force constitute the first-order photokinetic forces. The first term on the right-hand side of Eq. (5) describes the manifestly conservative intensity-gradient force responsible for trapping by single-beam optical traps such as optical tweezers [5]. The second describes a non-conservative contribution to radiation pressure that is directed by gradients of the phase [14, 16, 17]. Recognizing that these phase gradients define the local wave vector [18] and that $\alpha_{e}^{\prime\prime}$ and $\alpha_{m}^{\prime\prime}$ are positive for conventional materials, the phase-gradient term tends to drive objects downstream along the beam of light.

Light’s polarization plays remarkably little role in first-order photokinetic effects. This can be appreciated because phase differences among the Cartesian components that determine the degree of circular polarization are absent from Eq. (5). Even so, polarization-dependent photokinetic effects are well documented in the experimental literature [19, 20, 21] and emerge naturally in numerical studies [19, 22, 23, 24]. Such effects can be attributed to higher-order photokinetic forces.

### II.2 Higher-order photokinetic forces

Higher-order photokinetic effects result both from forces exerted on higher-order induced multipoles and also from interference between different types and orders of multipole scattering [3]. The multipole polarizabilities that parametrize these higher-order terms generally become rapidly smaller with increasing multipole order, particularly for objects that are smaller than the wavelength of light. For Rayleigh particles, the leading-order corrections to the first-order dipole forces arise from the interference between the scattered electric and magnetic dipole fields, and from the induced electric quadrupole moment. The dipole-interference term in particular provides a mechanism for creating propagation-invariant tractor beams.

Dipole scattering by a small isotropic particle gives rise to a second-order force of the form [13, 3]

 $\boldsymbol{F}_{em}(\boldsymbol{r})=-\frac{k^{4}}{12\pi\epsilon_{0}c}\,\Re% \left\{\alpha_{e}\alpha_{m}^{\ast}\boldsymbol{E}(\boldsymbol{r},t)\times% \boldsymbol{B}^{\ast}(\boldsymbol{r},t)\right\}.$ (6)

Expressing this in terms of the fields’ amplitudes and phases yields four distinctive contributions,

 $\displaystyle\boldsymbol{F}_{em}(\boldsymbol{r})$ $\displaystyle=\frac{\epsilon k^{3}}{24\pi\epsilon_{0}}\,\Im\left\{\alpha_{e}% \alpha_{m}^{\ast}\right\}\,\nabla\sum_{j=1}^{3}u_{j}^{2}(\boldsymbol{r})$ (7a) $\displaystyle\quad-\frac{\epsilon k^{3}}{12\pi\epsilon_{0}}\,\Re\left\{\alpha_% {e}\alpha_{m}^{\ast}\right\}\,\sum_{j=1}^{3}u_{j}^{2}(\boldsymbol{r})\nabla% \phi_{j}(\boldsymbol{r})$ (7b) $\displaystyle\quad-\frac{\epsilon k^{3}}{12\pi\epsilon_{0}}\,\Re\left\{\alpha_% {e}\alpha_{m}^{\ast}\right\}\,\nabla\cdot\mathsf{T}^{(a)}(\boldsymbol{r})$ (7c) $\displaystyle\quad-\frac{\epsilon k^{3}}{12\pi\epsilon_{0}}\,\Im\left\{\alpha_% {e}\alpha_{m}^{\ast}\right\}\nabla\cdot\mathsf{T}^{(s)}(\boldsymbol{r}),$ (7d)

where $\epsilon$ is the permeability of the medium and where we have introduced the antisymmetric $(a)$ and symmetric $(s)$ components of the Maxwell stress tensor,

 $\displaystyle T^{(a)}_{ij}(\boldsymbol{r})$ $\displaystyle=u_{i}(\boldsymbol{r})u_{j}(\boldsymbol{r})\,\sin\!\big{(}\phi_{j% }(\boldsymbol{r})-\phi_{i}(\boldsymbol{r})\big{)}\quad\text{and}$ (8a) $\displaystyle T^{(s)}_{ij}(\boldsymbol{r})$ $\displaystyle=u_{i}(\boldsymbol{r})u_{j}(\boldsymbol{r})\,\cos\!\big{(}\phi_{j% }(\boldsymbol{r})-\phi_{i}(\boldsymbol{r})\big{)}.$ (8b)

Equation (7) is the complete result for the dipole-interference contribution, including both electric and magnetic contributions. The force it describes tends to be weaker than the first-order dipole force described by Eq. (5) because $k^{3}\left|\alpha_{e}\right|\ll 1$ in the Rayleigh regime. Even so, each term in Eq. (7) has interesting implications for optical micromanipulation.

The intensity-gradient term in Eq. (7a) either strengthens or weakens the first-order trapping force depending on material parameters that affect the sign of $\Im\left\{\alpha_{e}\alpha_{m}^{\ast}\right\}$. Because $\Re\left\{\alpha_{e}\alpha_{m}^{\ast}\right\}\gtrsim 0$ for conventional materials, the phase-gradient term in Eq. (7b) describes a retrograde force that could provide the basis for a tractor beam were it not generally weaker than the first-order phase-gradient term from Eq. (5).

Equations (7c) and (7d) describe optical forces with properties distinct from those of dipole forces, and thus create new avenues for optical micromanipulation. We can interpret Eq. (7c) by noting that

 $\boldsymbol{s}(\boldsymbol{r})=\frac{\epsilon}{2\omega}\Im\left\{\boldsymbol{E% }^{\ast}(\boldsymbol{r},t)\times\boldsymbol{E}(\boldsymbol{r},t)\right\}$ (9a) is the time-averaged spin angular momentum density carried by the beam of light [25, 18, 26, 27]. From this, we obtain the general result, $\nabla\cdot\mathsf{T}^{(a)}(\boldsymbol{r})=\frac{\omega}{\epsilon}\,\nabla% \times\boldsymbol{s}(\boldsymbol{r}),$ (9b)

which identifies Eq. (7c) as a force arising from the curl of the spin angular momentum density in a nonuniform beam of light [21]. Spin-curl forces have been a subject of active research [28, 29, 30, 21, 15, 31, 27] and here may be identified as a second-order photokinetic effect. They may be directed along the axis of propagation, and thus can contribute to the retrograde forces in a tractor beam.

The remaining term in Eq. (7d) describes a complementary curl-free force field directed along the local axis of linear polarization. This term also can give rise to a retrograde force in longitudinally polarized traveling waves of the kind observed in experiments on crossed plane waves [8].

Induced-quadrupole contributions to the optical force also arise at second order in the multipole expansion and so may be comparable in strength to the dipole-interference force. This contribution has the form [3]

 $\boldsymbol{F}_{q}(\boldsymbol{r})=\frac{1}{4}\Re\left\{\sum_{j,k}Q_{jk}(% \boldsymbol{r},t)\nabla\partial_{j}E_{k}^{\ast}(\boldsymbol{r},t)\right\},$ (10)

where the induced quadrupole moment is

 $Q_{jk}(\boldsymbol{r},t)=\frac{1}{2}\gamma_{e}\left[\partial_{j}E_{k}(% \boldsymbol{r},t)+\partial_{k}E_{j}(\boldsymbol{r},t)\right]$ (11)

for a neutral particle with scalar quadrupole susceptibility, $\gamma_{e}$. The induced-quadrupole force is symmetric under exchange of the Cartesian components of $\boldsymbol{E}(\boldsymbol{r},t)$. Consequently, it does not depend on the spin angular momentum of the light, and does not contribute to spin-dependent photokinetic effects.

Second-order forces are weaker than first-order forces for Rayleigh particles. Higher multipole contributions will be weaker still. We will focus, therefore, on identifying characteristics of propagation-invariant modes of light that are useful for creating first-order tractor beams.

## III Propagation-invariant tractor beams

Achieving long-ranged retrograde transport without active intervention [32] requires a tractor beam to be based on a weakly-diffracting or propagation-invariant mode. Ideally, such a mode would carry retrograde momentum [33] along its entire length, and thus exert retrograde radiation pressure. It is sufficient, however, for the axial component of the beam’s momentum density to be smaller than the plane-wave limit. In that case, an illuminated particle can scatter additional momentum downstream, and thus recoil upstream. The archetype for such a tractor beam is the Bessel beam [34, 35, 36, 37], which originally was proposed for this application in the context of acoustic waves [1].

Bessel beams can be generated from the longitudinal vector potential

 $\displaystyle\boldsymbol{A}_{m,\alpha}(\boldsymbol{r})$ $\displaystyle=u_{m,\alpha}\,A_{m,\alpha}(\boldsymbol{r})\hat{z},\quad\text{where}$ (12a) $\displaystyle A_{m,\alpha}(\boldsymbol{r})$ $\displaystyle=J_{m}(kr\sin\alpha)\,e^{im\theta}\,e^{ikz\cos\alpha}$ (12b)

is the scalar cylindrical harmonic that is characterized by convergence angle $\alpha$ and an integer winding number $m$. The overall amplitude, $u_{m,\alpha}$, sets the intensity of the beam. Equation (12) is expressed in cylindrical coordinates, $\boldsymbol{r}=(\rho,\theta,z)$, aligned with the propagation direction, $\hat{z}$. The fields in a Bessel beam may be obtained from $\boldsymbol{A}_{m,\alpha}(\boldsymbol{r})$ as [34, 35]

 $\boldsymbol{E}_{m,\alpha}^{\text{TE}}(\boldsymbol{r},t)=\nabla\times% \boldsymbol{A}_{m,\alpha}(\boldsymbol{r})\,e^{-i\omega t},$ (13)

with the magnetic field following from Eq. (1b). Equation (13) describes a transverse electric (TE) Bessel beam that is equivalent to the form presented in [3]. The corresponding transverse magnetic (TM) mode can be obtained by exchanging the electric and magnetic fields in Eq. (13). Bessel beams with more general polarizations may be obtained as linear superpositions of these TE and TM modes. The range of non-diffracting propagation is proportional to $\csc\alpha$, and so is longer for beams with smaller convergence angles, $\alpha$ [37, 38].

Because the Bessel beams form a complete orthogonal basis, any propagation-invariant mode can be expressed as a linear superposition of Bessel beams. Those superpositions comprising Bessel beams with the same convergence angle, $\alpha$, display continuous propagation invariance. Any such non-diffracting beam lacks axial intensity gradients as a matter of definition, so that $\partial_{z}u^{2}(\boldsymbol{r})=0$. The remaining phase-gradient contribution to the axial dipole force is purely repulsive for any material with $\alpha_{e}^{\prime\prime}>0$ and $\alpha_{m}^{\prime\prime}>0$. From this, we conclude that no continuously propagation-invariant beam of light can serve as a general-purpose first-order tractor beam.

Non-diffracting modes nonetheless may act as tractor beams for particles whose scattering properties favor retrograde forces from higher-order contributions [3]. To illustrate this, we consider the photokinetic force exerted by TE Bessel beams and contrast this with the performance of TM Bessel beams. The axial component of the force then has the form

 $\boldsymbol{F}(\boldsymbol{r})\cdot\hat{z}=\frac{1}{2}k\,u_{m,\alpha}^{2}\frac% {J_{m}^{2}(kr\sin\alpha)}{r^{2}}\,\cos\alpha\,f(\alpha,x),$ (14)

where $x=kr\sin\alpha$ and where $f(\alpha,x)$ is a polarization-dependent factor describing the particle’s coupling to the light. The axial force vanishes in the strongly converging limit because the Bessel beam becomes a standing wave at $\alpha=\pi/2$.

The coupling factor for a TE-polarized Bessel beam due to dipole and dipole-interference contributions is

 $\displaystyle f^{\text{TE}}(\alpha,x)=g_{m}^{2}(x)\left(\alpha_{e}^{\prime% \prime}+\alpha_{m}^{\prime\prime}-\frac{\epsilon k^{3}}{6\pi\epsilon_{0}}\Re% \left\{\alpha_{e}\alpha_{m}^{\ast}\right\}\right)\\ \displaystyle+\sin^{2}\alpha\,[x^{2}-g_{m}^{2}(x)]\,\alpha_{m}^{\prime\prime},$ (15)

with the mode-dependent geometric factor

 $g_{m}^{2}(x)=m^{2}+x^{2}\frac{{J_{m}^{\prime}}^{2}(x)}{{J_{m}}^{2}(x)}.$ (16)

The equivalent expression for TM modes is obtained by transposing $\alpha_{e}$ and $\alpha_{m}$.

The electric and magnetic quadrupole contributions to the axial force are strictly repulsive, with the electric quadrupole repulsion being somewhat weaker in TM-polarized modes. Quadrupole contributions therefore do not help to create a second-order tractor beam.

Considering, for the moment, only those contributions to the axial force arising from dipole terms, the direction of motion depends on the sign of $f^{\text{TE}}(\alpha,x)$, and therefore on the particle’s position within the beam. Retrograde forces appear at positions $x$ that satisfy

 $\frac{g_{m}^{2}(x)}{x^{2}}<\frac{\sin^{2}\alpha\,\alpha_{m}^{\prime\prime}}{% \frac{\epsilon k^{3}}{6\pi\epsilon_{0}}\,\Re\left\{\alpha_{e}\alpha_{m}^{\ast}% \right\}-\alpha_{e}^{\prime\prime}-\cos^{2}\alpha\,\alpha_{m}^{\prime\prime}}.$ (17)

The absence of solutions in the limit of small $\alpha$ demonstrates that weakly converging Bessel beams generally act as repulsors. More strongly converging modes can act as tractor beams if the radial component of the force vanishes at positions that satisfy Eq. (17). This occurs, for example, in TE modes with $m=\pm 1$ along the principal intensity maximum at $r=0$. Such beams can stably pull an object upstream if the convergence angle satisfies

 $\sin^{2}\alpha>1+\frac{\alpha_{e}^{\prime\prime}}{\alpha_{m}^{\prime\prime}}-% \frac{\epsilon k^{3}}{6\pi\epsilon_{0}}\left(\alpha_{e}^{\prime\prime}+\alpha_% {e}^{\prime}\,\frac{\alpha_{m}^{\prime}}{\alpha_{m}^{\prime\prime}}\right).$ (18)

Because Eq. (17) neglects the quadrupole contribution, it constitutes a necessary condition for a Bessel beam to act as a second-order tractor beam, but is not a sufficient condition.

Numerical studies involving all orders of multipole scattering confirm that Bessel beams can act as tractor beams for some objects that are larger than the wavelength of light [3, 4, 11]. The present considerations demonstrate that the necessary pulling force does not arise at first order, but rather emerges from a competition between first-order repulsion and higher-order retrograde forces. Because such tractor beams require comparatively steep convergence angles, moreover, they do not lend themselves to long-range transport.

## IV First-order solenoidal tractor beams

Relaxing the condition of continuous propagation invariance in favor of discrete propagation invariance,

 $u^{2}(r,\theta,z+\Delta z)=u^{2}(r,\theta,z),$ (19)

creates opportunities for strong, long-ranged, general-purpose tractor beams that operate at first order. Axial intensity gradients in a discretely propagation invariant beam can establish retrograde forces that counteract repulsive phase-gradient forces. This is the principle by which an optical conveyor establishes an axial array of traps [32]. Achieving continuous retrograde transport in such a beam requires a continuous path through the beam along which the net optical force acts in the retrograde direction. Optical conveyors do not have this property, but waves displaying screw invariance can. These include spiral waves [39] and solenoidal beams [2]. Although solenoidal waves have been described and their efficacy as tractor beams has been demonstrated experimentally [2], the nature of their force fields has not been explored previously. Applying the theory of photokinetic effects reveals that solenoidal modes can act as first-order tractor beams.

Discretely propagation-invariant beams can be created by superposing two mutually coherent Bessel beams with different convergence angles. Because we are interested in long-range transport, we choose component Bessel beams with small convergence angles, which may be described in the paraxial approximation as TEM modes with uniform polarization, $\hat{\epsilon}$. The superposition then has the form

 $\boldsymbol{E}(\boldsymbol{r},t)=\left[u_{m,\alpha}A_{m,\alpha}(\boldsymbol{r}% ,t)+u_{n,\beta}A_{n,\beta}(\boldsymbol{r},t)\right]\hat{\epsilon}.$ (20)

The resulting axial interference pattern is spatially periodic, as required, and stationary. Superposing modes with the same winding number creates optical conveyors [40, 32, 41]. Superposing modes with different winding numbers creates a solenoidal wave. We will order the contributions with $\left|m\right|>\left|n\right|$. The original description of solenoidal waves [2] involved superpositions of large numbers of Bessel modes. Being composed of just two modes, the solenoidal waves presented here may be considered minimal realizations of the genus.

We ensure that a solenoidal mode traps particles at a radius $r=R$ by choosing the convergence angles $\alpha$ and $\beta$ so that extremum $p$ of mode $m$ overlaps extremum $q$ of mode $n$. This condition ensures $\partial_{r}\left|\boldsymbol{E}(\boldsymbol{r},t)\right|^{2}=0$ along $r=R$. Equation (5) then ensures that $\left.\hat{r}\cdot\boldsymbol{F}_{e}(\boldsymbol{r})\right|_{r=R}=0$. The intensity maximum at $r=R$ acts as a potential energy well along the radial direction for light-seeking particles with $\alpha_{e}^{\prime}>0$. It is achieved by setting

 $\frac{\sin\alpha}{j^{\prime}_{m,p}}=\frac{\sin\beta}{j^{\prime}_{n,q}}=\frac{1% }{kR},$ (21)

where $x=j_{m,p}^{\prime}$ is the $p$-th root of $J_{m}^{\prime}(x)=0$.

Interference is optimized if the two components contribute equally to the amplitude, $u_{0}$, at $r=R$. This establishes a relationship between the amplitudes

 $u_{\alpha,m}J_{m}(j_{m,p}^{\prime})=u_{\beta,n}J_{n}(j_{n,q}^{\prime})=\frac{1% }{2}u_{0}.$ (22)

The resulting superposition has the form

 $\boldsymbol{E}(\boldsymbol{r},t)=u(\boldsymbol{r})\,e^{i\phi(\boldsymbol{r})}e% ^{-i\omega t}\hat{x},$ (23a) with amplitude and phase profiles along $r=R$ of $\displaystyle\left.u(\boldsymbol{r})\right|_{r=R}$ $\displaystyle=u_{0}\,\cos\!\left(\frac{m-n}{2}\,\left[\theta-2\pi\frac{z}{% \Delta z}\right]\right)\quad\text{and}$ (23b) $\displaystyle\left.\phi(\boldsymbol{r})\right|_{r=R}$ $\displaystyle=\ell\theta+\frac{\cos\alpha+\cos\beta}{2}\,kz.$ (23c) This is a helical mode that carries orbital angular momentum [42] proportional to its winding number, $\ell=\frac{m+n}{2}.$ (23d) The axial pitch of the solenoid’s spiraling intensity distribution is $\Delta z=\frac{2\pi}{k}\frac{m-n}{\cos\alpha-\cos\beta}.$ (23e)

Interestingly, solenoidal waves with half-integer winding numbers do not differ qualitatively from those with integer winding numbers. The $2\ell$-fold sign changes in $u(\boldsymbol{r})$, as defined in Eq. (23b), eliminate the transverse vortices that otherwise would form in a helical beam with non-integer winding number [43, 44].

The transverse intensity distribution is proportional to $u^{2}(\boldsymbol{r})$ and has a number of maxima equal to $\left|m-n\right|$. These maxima spiral once around $\hat{z}$ each time the axial position $z$ advances by $\Delta z$. The spiral is right-handed for positive values of $\Delta z$.

For simplicity, we will consider this wave’s influence on dielectric Rayleigh particles, for which the electric dipole force described by Eq. (5) dominates the magnetic dipole force and all other multipole contributions. We then define the locus, $\boldsymbol{r}_{s}(z)=(R,\theta_{s},z)$, of the solenoidal trap to be the set of points along which intensity-gradient forces compensate the axial component of the radiation pressure, $F_{z}(\boldsymbol{r}_{s})=0$. This condition is satisfied along the spiral curve

 $\theta_{s}(z)=2\pi\frac{z}{\Delta z}+\theta_{0},$ (24a) with angular offset, $\theta_{0}=\frac{2}{m-n}\,\tan^{-1}\!\left(\frac{\alpha_{e}^{\prime\prime}}{% \alpha_{e}^{\prime}}\frac{\cos\alpha+\cos\beta}{\cos\alpha-\cos\beta}\right).$ (24b)

The path that a particular object takes through the beam runs parallel to this. Such solutions exist for small objects regardless of the light-scattering properties encoded in $\alpha_{e}$. Whether or not these solutions act as traps depends on whether or not the particle also is stably trapped in the radial direction. Small dielectric particles are stably trapped along the bright locus described by Eq. (23); reflecting and absorbing particles might not be.

A solenoidal trap acts as a tractor beam if the azimuthal component of the radiation pressure, $F_{\theta}(\boldsymbol{r}_{s})$, drives the trapped particle into a region where the axial component of the force, $F_{z}(\boldsymbol{r}_{s})$, points in the retrograde direction, $F_{\theta}(\boldsymbol{r}_{s})\,\left.\partial_{\theta}F_{z}(\boldsymbol{r})% \right|_{\boldsymbol{r}=\boldsymbol{r}_{s}}<0$. This condition is satisfied if

 $n\cos\alpha>m\cos\beta.$ (25)

Flipping the sense of this inequality transforms a tractor beam into a repulsor that pushes objects downstream rather than pulling them upstream.

Equation (25) is the principal design requirement for creating a general-purpose first-order tractor beam. With this constraint, Eq. (21) then requires $q>p$. Minimizing $\alpha$ and $\beta$ to maximize pulling range then requires $q=p+1$. Maximizing the pulling force by maximizing the intensity along the solenoid’s principal spiral then yields $p=1$ and $q=2$ as optimal choices.

Equations (21), (22) and (25) constitute previously lacking design criteria for general-purpose first-order tractor beams. As described, they establish conditions on the individual Bessel components’ convergence angles, $\alpha$ and $\beta$, and winding numbers, $m$ and $n$. Equivalently, these conditions may be expressed in terms of the radius, $R$, pitch, $\Delta z$, winding number, $\ell$, and non-diffracting range of the desired tractor beam. Previously, functional solutions had to be sought numerically [2].

The example of a solenoidal tractor beam presented in Fig. 1(a) was designed according to these principles. Figure 1(b) shows the computed trajectory of a $50\,\mathrm{nm}$-diameter silica sphere moving through water in its force field. Because solenoidal tractor beams operate at dipole order, they should be inherently stronger than higher-order tractor beams based on continuously invariant modes.

Accounting for the magnetic field’s contribution to the optical force shifts $\theta_{0}$ by an amount that depends on the relative magnitudes of $\alpha_{e}$ and $\alpha_{m}$. The net effect is to strengthen the pulling force, at least for conventional materials.

Polarization also can enhance a solenoidal tractor beam’s performance. For example, Eq. (7c) shows that a solenoid beam created with right-circularly polarized light has a transverse spin-curl force of the form

 $\boldsymbol{F}_{\text{sc}}(\boldsymbol{r})=\frac{\epsilon k^{3}}{24\pi\epsilon% _{0}}(\alpha_{e}^{\prime}\alpha_{m}^{\prime}+\alpha_{e}^{\prime\prime}\alpha_{% m}^{\prime\prime})\left(\frac{\hat{r}}{r}\partial_{\theta}-\hat{\theta}% \partial_{r}\right)u^{2}(\boldsymbol{r}).$ (26)

When applied to a right-handed solenoidal tractor beam, this increases $F_{\theta}(r)$ for $r and decreases it for $r>R$. By moderating the trapped particle’s speed in this way, $\boldsymbol{F}_{\text{sc}}(\boldsymbol{r})$ contributes to a kinematic trap that tends to localize an object near $r=R$ as it is transported by the beam.

## V Discussion

The two-component solenoidal tractor beams described here illustrate the process by which first-order tractor beams may be designed and optimized using the theory of photokinetic effects. More sophisticated superpositions [2] may have superior properties, but at the cost of additional complexity. For some applications, solenoidal waves based on other helical modes, such as Laguerre-Gaussian beams, may be preferable.

Beyond its immediate application to designing tractor beams, the theory of photokinetic effects offers additional benefits for optical micromanipulation. Expressing optical forces in terms of experimentally accessible quantities helps to elucidate the origin of experimentally observed phenomena. The polarization-dependent contribution from Eq. (7d), for example, explains why linearly polarized optical tweezers exert anisotropic in-plane forces [23, 16]. The circulation of optically isotropic spheres in circularly polarized optical tweezers similarly may be identified with the spin-curl force from Eq. (7c). Polarization-dependent size selectivity in the force created by crossed plane waves [8] can be explained as arising from competition between the phase-gradient terms in Eq. (5) and Eq. (7b). Such observations provide guidance for designing and optimizing new modes of optical micromanipulation, including tractor beams.

The development presented here directly addresses the forces experienced by objects that are substantially smaller than the wavelength of light. For larger particles, the forces described by Eqs. (5) and (7) are augmented and possibly dominated by higher-order multipole contributions. For Rayleigh particles, we expect the first-order forces described by Eq. (5) to be substantially stronger than any higher-order contributions, including those described by Eq. (7). Applying these considerations to the particular case of tractor beams, we conclude that Bessel beams generally push small objects whereas appropriately designed solenoidal waves can pull them.

## Acknowledgments

This work was supported by the National Science Foundation under Grant No. DMR-1305875. AY acknowledges support from a NASA Space Technology Research Fellowship under Award No. NNX15AQ40H.

## References

• [1] P. L. Marston. “Axial radiation force of a Bessel beam on a sphere and direction reversal of the force.” J. Acoust. Soc. Am. 120, 3518–3524 (2006).
• [2] S.-H. Lee, Y. Roichman and D. G. Grier. “Optical solenoid beams.” Opt. Express 18, 6988–6993 (2010).
• [3] J. Chen, J. Ng, Z. Lin and C. T. Chan. “Optical pulling force.” Nature Photonics 5, 531–534 (2011).
• [4] A. Novitsky, C.-W. Qiu and H. Wang. ‘‘Single gradientless light beam drags particles as tractor beams.” Phys. Rev. Lett. 107, 203601 (2011).
• [5] 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, 288–290 (1986).
• [6] A. Ashkin. “Forces of a single-beam gradient laser trap on a dielectric sphere in the ray optics regime.” Biophys. J. 61, 569–582 (1992).
• [7] A. Ashkin. “Optical trapping and manipulation of neutral particles using lasers.” Proc. Natl. Acad. Sci. U.S.A. 94, 4853–4860 (1997).
• [8] 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, 123–127 (2013).
• [9] E. E. Smith. “Spacehounds of IPC.” Amazing Stories 6, 1–100 (1931).
• [10] V. Shvedov, A. R. Davoyan, C. Hnatovsky, N. Engheta and W. Krolikowski. “A long-range polarization-controlled optical tractor beam.” Nature Phot. 8, 846–850 (2014).
• [11] A. Novitsky, C.-W. Qiu and A. Lavrinenko. “Material-independent and size-independent tractor beams for dipole objects.” Phys. Rev. Lett. 109, 023902 (2012).
• [12] P. C. Chaumet and M. Nieto-Vesperinas. “Time-averaged total force on a dipolar sphere in an electromagnetic field.” Opt. Lett. 25, 1065–1067 (2000).
• [13] P. C. Chaumet and A. Rahmani. “Electromagnetic force and torque on magnetic and negative-index scatterers.” Opt. Express 17, 2224–2234 (2009).
• [14] Y. Roichman, B. Sun, Y. Roichman, J. Amato-Grill and D. G. Grier. “Optical forces arising from phase gradients.” Phys. Rev. Lett. 100, 013602 (2008).
• [15] D. B. Ruffner and D. G. Grier. “Comment on ‘Scattering Forces from the Curl of the Spin Angular Momentum of a Light Field’.” Phys. Rev. Lett. 111, 059301 (2013).
• [16] Y. Roichman, B. Sun, A. Stolarski and D. G. Grier. “Influence of non-conservative optical forces on the dynamics of optically trapped colloidal spheres: The fountain of probability.” Phys. Rev. Lett. 101, 128301 (2008).
• [17] P. Y. Wu, R. X. Huang, C. Tischer, A. Jonas and E. L. Florin. “Direct measurement of the nonconservative force field generated by optical tweezers.” Phys. Rev. Lett. 103, 108101 (2009).
• [18] M. V. Berry. “Optical currents.” J. Opt. A 11, 094001 (2009).
• [19] A. Rohrbach. “Stiffness of optical traps: Quantitative agreement between experiments and electromagnetic theory.” Phys. Rev. Lett. 95, 168102 (2005).
• [20] E. Madadi, A. Samadi, M. Cheraghian and N. S. Reihani. “Polarization-induced stiffness asymmetry of optical tweezers.” Opt. Lett. 37, 3519–3521 (2012).
• [21] D. B. Ruffner and D. G. Grier. “Optical forces and torques in non-uniform beams of light.” Phys. Rev. Lett. 108, 173602 (2012).
• [22] A. R. Zahkarian, P. Polynkin, M. Mansuripur and J. V. Moloney. “Single-beam trapping of micro-beads in polarized light: Numerical simulations.” Opt. Express 14, 3660–3676 (2006).
• [23] R. S. Dutra, N. B. Viana, P. A. Maia Neto and H. M. Nussenzveig. “Polarization effects in optical tweezers.” J. Opt. A 9, S221–S227 (2007).
• [24] B. Sun, Y. Roichman and D. G. Grier. “Theory of holographic optical trapping.” Opt. Express 16, 15765–15776 (2008).
• [25] M. V. Berry and M. R. Dennis. “Polarization singularities in isotropic random vector waves.” Proc. Royal Soc. A 457, 141–155 (2001).
• [26] A. Bekshaev, K. Y. Bliokh and M. Soskin. “Internal flows and energy circulation in light beams.” J. Opt. 13, 053001 (2011).
• [27] K. Bliokh, J. Dressel and F. Nori. “Conservation of the spin and orbital angular momentum in electromagnetism.” New J. Phys. 16, 093037 (2014).
• [28] N. B. Simpson, K. Dholakia, L. Allen and M. J. Padgett. “Mechanical equivalence of spin and orbital angular momentum of light: An optical spanner.” Opt. Lett. 22, 52–54 (1997).
• [29] S. Albaladejo, M. I. Marqués, M. Laroche and J. J. Sáenz. “Scattering forces from the curl of the spin angular momentum.” Phys. Rev. Lett. 102, 113602 (2009).
• [30] X. L. Wang, J. Chen, Y. Li, J. Ding, C. S. Guo and H. T. Wang. “Optical orbital angular momentum from the curl of polarization.” Phys. Rev. Lett. 105, 253602 (2010).
• [31] A. Y. Bekshaev. “Subwavelength particles in an inhomogeneous light field: optical forces associated with the spin and orbital energy flows.” J. Opt. 15, 044004 (2013).
• [32] D. B. Ruffner and D. G. Grier. “Optical conveyors: A class of active tractor beams.” Phys. Rev. Lett. 109, 163903 (2012).
• [33] M. V. Berry. “Quantum backflow, negative kinetic energy, and optical retro-propagation.” J. Phys. A 43, 415302 (2010).
• [34] X. Zhou. “Normal and abnormal vector wave functions and their conversion relations.” Scientia Sinica 27, 1226–1232 (1984).
• [35] Z. Bouchal and M. Olvik. “Non-diffractive vector Bessel beams.” J. Mod. Opt. 42, 1555–1566 (1995).
• [36] J. Durnin. “Exact-solutions for nondiffracting beams. 1. The scalar theory.” J. Opt. Soc. Am. A 4, 651–654 (1987).
• [37] J. Durnin, J. J. Miceli, Jr and J. H. Eberly. “Diffraction-free beams.” Phys. Rev. Lett. 58, 1499–1501 (1987).
• [38] J. Arlt, V. Garcés-Chávez, W. Sibbett and K. Dholakia. “Optical micromanipulation using a Bessel light beam.” Opt. Commun. 197, 239–245 (2001).
• [39] S. Chávez-Cerda, G. S. McDonald and G. H. C. New. “Nondiffracting beams: travelling, standing, rotating and spiral waves.” Opt. Commun. 123, 225–233 (1996).
• [40] 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, 174101 (2005).
• [41] D. B. Ruffner and D. G. Grier. “Universal, strong and long-ranged trapping by optical conveyors.” Opt. Express 22, 26834–26853 (2014).
• [42] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw and J. P. Woerdman. “Orbital angular-momentum of light and the transformation of Laguerre-Gaussian laser modes.” Phys. Rev. A 45, 8185–8189 (1992).
• [43] M. V. Berry. “Optical vortices evolving from helicoidal integer and fractional phase steps.” J. Opt. A 6, 259–269 (2004).
• [44] J. Leach, E. Yao and M. J. Padgett. “Observation of the vortex structure of a non-integer vortex beam.” New J. Phys. 6, 71 (2004).