# Acoustokinetics: Crafting force landscapes from sound waves

Mohammed A. Abdelaziz    David G. Grier Department of Physics and Center for Soft Matter Research, New York University, New York, NY 10003
###### Abstract

Factoring the pressure field of a harmonic sound wave into its amplitude and phase profiles provides the foundation for an analytical framework for studying acoustic forces that not only provides novel insights into the forces exerted by specified sound waves, but also addresses the inverse problem of designing sound waves to implement desired force landscapes. We illustrate the benefits of this acoustokinetic framework through case studies of purely nonconservative force fields, standing waves, pseudo-standing waves, and tractor beams.

## I Introduction

Structured sound waves exert forces and torques that can be harnessed to transport, sort and organize insonated objects [1, 2, 3, 4]. Applications include non-contact processing of sensitive [5, 6] and hazardous [7] materials, flow focusing for materials analysis and medical diagnostics [8], and automated remote manipulation for research [9]. Rapidly growing interest in harnessing acoustic forces has inspired a fundamental reassessment of the physics of wave-mediated forces. Recent developments in the theory of acoustic forces [10, 11, 12, 13, 14] parallel the analogous theory of optical forces [15, 16]. Both offer valuable and often surprising insights into the elementary mechanisms of wave-matter interactions. The acoustokinetic framework introduced here addresses the complementary inverse problem: identifying what wave will create a desired force landscape.

The inverse problem for optical forces recently has been rendered more tractable by expressing the electromagnetic field in terms of its real-valued amplitudes and phases along each Cartesian coordinate [17, 18]. This approach is called the theory of photokinetic effects and yields useful analytic expressions for the performance of optical traps [19] including design criteria for optical tractor beams [18]. Here, we show that an analogous factorization of the pressure field in sound waves is similarly useful for understanding and implementing acoustic manipulation. We illustrate the utility of this acoustokinetic framework through case studies on nonconservative acoustic force fields, standing and pseudo-standing waves, and acoustic tractor beams.

### I.1 Light: Photokinetic analysis

We develop acoustokinetics by analogy to photokinetics and therefore briefly review the theory of optical forces. A small particle immersed in an electromagnetic wave develops an electric dipole moment proportional to the local field. This induced dipole experiences a time-averaged Lorentz force in gradients of the field that can be expressed as [16]

 $\mathbf{F}_{e}(\mathbf{r})=\frac{1}{2}\Re\left\{\alpha_{e}\sum_{j=1}^{3}E_{j}(% \mathbf{r})\nabla E_{j}^{\ast}(\mathbf{r})\right\},$ (1)

where $E_{j}(\mathbf{r})$ is the $j$-th Cartesian coordinate of the electric field and $\alpha_{e}$ is the particle’s complex dipole polarizability. Expressing the components of the electric field in terms of their real-valued amplitude and phase profiles,

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

yields the surprisingly simple expression [17],

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

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

The first term on the right-hand side of Eq. (3) is the manifestly conservative intensity-gradient force responsible for single-beam optical traps such as optical tweezers [20]. The second describes a nonconservative force [21] that is directed by phase gradients [22]. Phase-gradient forces tend to drive trapped particles out of thermodynamic equilibrium with their supporting media [22, 23], mediate the transfer of light’s orbital angular momentum [24, 25, 26, 27], and have been used to create light-driven micromachines such as pumps [28], mixers [29] and optical tractor beams [30]. Even non-absorbing dielectric particles experience nonconservative optical forces because of radiative contributions to the dipole polarizability [31, 32].

The dipole-order expression in Eq. (3) accurately describes the forces experienced by particles with radii, $a_{p}$, that are small enough to satisfy the Rayleigh criterion, $ka_{p}<1$, where $k$ is the wavenumber of light. In the Rayleigh regime, the conservative intensity-gradient force generally dominates the light-matter interaction because $\alpha_{e}^{\prime}$ scales as $(ka_{p})^{3}$, whereas $\alpha_{e}^{\prime\prime}$ scales as $(ka_{p})^{6}$.

### I.2 Sound: Acoustic radiation forces

The analogous dipole-order acoustic radiation force experienced by a small particle in a harmonic sound field may be expressed in terms of the pressure, $p(\mathbf{r},t)$, as [14]

 $\mathbf{F}(\mathbf{r})=\frac{1}{2}\Re\left\{\alpha_{a}p\nabla p^{\ast}+\beta_{% a}k^{-2}(\nabla p\cdot\nabla)\nabla p^{\ast}\right\},$ (4)

where the coefficients $\alpha_{a}$ and $\beta_{a}$ play the role of dipole and quadrupole polarizabilities, respectively. Expressing $\mathbf{F}(\mathbf{r})$ in terms of multipole polarizabilities clarifies the analogy with photokinetics. Lengths in Eq. (4) are scaled by the wavenumber, $k=\omega/c_{m}$, where $\omega$ is the sound’s frequency and $c_{m}$ is its speed in the medium. Equation (4) applies to inviscid fluids, for which the pressure satisfies the scalar wave equation

 $\nabla^{2}p=-k^{2}p.$ (5)

Our focus on traveling waves in inviscid media is inspired by our interest in developing new modalities of long-ranged non-contact manipulation. Long-range manipulation is facilitated by minimizing acoustic losses in the medium. This can be achieved in air by working at frequencies below $50\,\mathrm{kHz}$ [3, 33], for which the acoustic attenuation is less than $2\,\mathrm{dB}\,\mathrm{m}^{-1}$ under standard conditions [34] and scales as $\omega^{2}$ for lower frequencies. These conditions also minimize the influence of acoustic streaming forces, which ordinarily compete with acoustic radiation forces in viscous media and in inviscid media bounded by confining surfaces [35].

An object’s dipole and quadrupole polarizabilities generally depend on its size, shape and composition as well as the frequency of the sound and the properties of the fluid medium. For simplicity and concreteness, we will specialize to the case of a spherical scatterer of radius $a_{p}$ that is composed of a material of density $\rho_{p}$ and sound speed $c_{p}$. Such an object’s response to the sound field is characterized by the polarizabilities [14]

 $\displaystyle\alpha_{a}$ $\displaystyle=\frac{4\pi a_{p}^{3}}{3\rho_{m}c_{m}^{2}}\,f_{0}\,\left[-1+i% \frac{1}{3}(f_{0}+f_{1})(ka_{p})^{3}\right]$ (6a) $\displaystyle\beta_{a}$ $\displaystyle=\frac{2\pi a_{p}^{3}}{\rho_{m}c_{m}^{2}}\,f_{1}\,\left[1+i\frac{% 1}{6}f_{1}(ka_{p})^{3}\right],$ (6b)

where the monopole coupling coefficient,

 $f_{0}=1-\frac{\rho_{m}c_{m}^{2}}{\rho_{p}c_{p}^{2}},$ (7a) depends on the compressibility mismatch between the particle and the medium, and the dipole coupling coefficient, $f_{1}=2\frac{\rho_{p}-\rho_{m}}{2\rho_{p}+\rho_{m}},$ (7b)

gauges the density mismatch. The expressions in Eq. (6) are valid for $ka_{p}<1$ and constitute the leading-order contributions for both the real parts of the polarizabilities, $\alpha_{a}^{\prime}$ and $\beta_{a}^{\prime}$, and also the imaginary parts, $\alpha_{a}^{\prime\prime}$ and $\beta_{a}^{\prime\prime}$.

## II Acoustokinetic Framework

Drawing on the analogy with photokinetics, we express the harmonic sound wave’s pressure field in terms of its amplitude and phase profiles:

 $p(\mathbf{r},t)=u(\mathbf{r})\,e^{i\varphi(\mathbf{r})}\,e^{-i\omega t}.$ (8)

The first term from Eq. (4) then yields

 $\mathbf{F}_{\alpha}(\mathbf{r})=\frac{1}{4}\alpha_{a}^{\prime}\,\nabla u^{2}+% \frac{1}{2}\alpha_{a}^{\prime\prime}\,u^{2}\nabla\varphi,$ (9a) which is directly analogous to Eq. (3) for the dipole-order force exerted by light. These contributions to the acoustic radiation force are depicted in Fig. 1. As in the optical case, $\alpha_{a}^{\prime}$ scales as $(ka_{p})^{3}$ and $\alpha_{a}^{\prime\prime}$ scales as $(ka_{p})^{6}$, which means that the conservative force generally dominates for small particles.

The second term from Eq. (4) arises from the velocity-matching condition at the sphere’s boundary and so has no analogue in optical radiation forces. It vanishes for density-matched particles ($\beta_{\alpha}=0$), which therefore behave exactly like dielectric particles in a light field, in agreement with previous results [36]. When expressed in terms of the amplitude and phase profiles, this term separates naturally into a conservative contribution,

 $\mathbf{F}_{\beta}^{c}(\mathbf{r})=\frac{1}{4}\beta_{a}^{\prime}\nabla\left[u^% {2}+\frac{1}{2}k^{-2}\nabla^{2}u^{2}\right],$ (9b)

that augments the conservative intensity-gradient force from $\mathbf{F}_{\alpha}(\mathbf{r})$ and a nonconservative contribution,

 $\displaystyle\mathbf{F}_{\beta}^{n}(\mathbf{r})=\frac{1}{4}\beta_{a}^{\prime% \prime}k^{-2}\left[(2k^{2}u^{2}+\nabla^{2}u^{2}+2u\nabla u\cdot\nabla)\nabla% \varphi\right.\\ \displaystyle\left.-(u\nabla^{2}\varphi+2u\nabla\varphi\cdot\nabla)\nabla u% \right],$ (9c)

that is directed both by phase gradients and also by amplitude gradients. The combination,

 $\mathbf{F}_{\beta}(\mathbf{r})=\mathbf{F}_{\beta}^{c}(\mathbf{r})+\mathbf{F}_{% \beta}^{nc}(\mathbf{r}),$ (9d)

captures the sphere’s leading-order coupling to the quadrupole components of the incident field. Unlike the optical case, where quadrupolar forces generally are weaker than dipole contributions, the two terms in $\mathbf{F}_{\beta}(\mathbf{r})$ can be comparable in magnitude to their counterparts in $\mathbf{F}_{\alpha}(\mathbf{r})$ because $\beta_{a}^{\prime}$ scales as $(ka_{p})^{3}$ and $\beta_{a}^{\prime\prime}$ scales as $(ka_{p})^{6}$. These density-dependent terms therefore can be used to exert control in ways that are not possible with light.

For very small particles satisfying $ka_{p}\ll 1$, the acoustic force field is dominated by the conservative terms proportional to $\alpha_{a}^{\prime}$ and $\beta_{a}^{\prime}$. This force is described by the classic Gor’kov potential [12] which is widely used to describe acoustic trapping phenomena [3, 37]. For larger particles, and for appropriately structured sound fields, non-conservative contributions proportional to $\alpha_{a}^{\prime\prime}$ and $\beta_{a}^{\prime\prime}$ can be significant, and even dominant [38, 39]. These contributions are not accounted for by the Gor’kov potential.

The acoustokinetic framework described by Eq. (9) is the principal contribution of this work. We now demonstrate its value through case studies on realizable sound fields with exceptional properties.

## III Applications of the Acoustokinetic Framework

### III.1 Designing purely nonconservative force fields

To illustrate how the acoustokinetic framework can address the inverse problem of designing sound waves to implement desired force landscapes, we use Eq. (9) to design harmonic sound waves that exert purely nonconservative forces. This is equivalent to requiring the conservative part of the acoustic radiation force to vanish, and thus requires us to look beyond the Gor’kov potential. Equations (9a) and (9b) show that this goal can be met if the particle is not density matched, $\beta_{a}^{\prime}\neq 0$, and if the pressure intensity satisfies the inhomogeneous Helmholtz equation,

 $\nabla^{2}u^{2}+2\left(1+\frac{\alpha_{a}^{\prime}}{\beta_{a}^{\prime}}\right)% k^{2}\,u^{2}=C.$ (10)

The undetermined constant $C$ distinguishes families of non-conservative sound waves for the class of objects with compatible values of $\alpha_{a}^{\prime}/\beta_{a}^{\prime}$. Solutions to Eq. (10) must be real-valued and must be paired with real-valued phase profiles that complete the description of the pressure field and satisfy the wave equation, Eq. (5).

 One interesting set of purely nonconservative solutions has the sinusoidal amplitude profile $u(\mathbf{r})=p_{0}\cos(q(x-y)),$ (11a) with spatial frequency $q=\frac{1}{2}\sqrt{1+\frac{\alpha_{a}^{\prime}}{\beta_{a}^{\prime}}}\,k.$ (11b) The associated phase profile, $\varphi(\mathbf{r})=kz\,\cos\gamma+(x+y)\sqrt{\frac{1}{2}k^{2}\sin^{2}\gamma-q% ^{2}},$ (11c) identifies this field as the superposition of two plane waves, each oriented at angle $\gamma$ relative to $\hat{z}$ and at angle $\theta=\cos^{-1}\left(-\frac{\alpha_{a}^{\prime}}{\beta_{a}^{\prime}}\right)$ (11d) relative to one another in the $(x,y)$-plane.

Under these conditions, the in-plane component of the radiation pressure exactly cancels the conservative intensity-gradient force. The remaining scattering force is sinusoidally modulated in the transverse plane. The result is a “picket fence” of parallel force lines, which is illustrated in Fig. 2 using the rotated coordinates

 $\displaystyle x^{\prime}$ $\displaystyle=x-y\quad\text{and}$ (12a) $\displaystyle y^{\prime}$ $\displaystyle=x+y$ (12b)

for clarity. In these coordinates, the net force,

 $\mathbf{F}(\mathbf{r}^{\prime})=\frac{1}{4}kp_{0}^{2}\,f(x^{\prime})\,\hat{F},$ (13a) is directed along $\hat{F}=\sqrt{\frac{1}{2}\sin^{2}\gamma-\frac{q^{2}}{k^{2}}}\,\hat{y}^{\prime}% +\cos\gamma\,\hat{z},$ (13b) and has an amplitude that varies with the transverse coordinate as $f(x^{\prime})=C_{1}+2(C_{2}-C_{1})\cos^{2}(qx^{\prime}).$ (13c) The scale and depth of the force landscape’s modulation depend on the object’s properties through $C_{1}=\beta_{a}^{\prime\prime}(1+\alpha_{a}^{\prime}/\beta_{a}^{\prime})$ and $C_{2}=\alpha_{a}^{\prime\prime}+\beta_{a}^{\prime\prime}$.

Picket fence modes for a given type of object are distinguished by the angle of inclination, $\gamma$. The range of possible angles is limited by Eq. (11c) to

 $\sin^{2}\gamma>\frac{2q^{2}}{k^{2}}.$ (14)

With this constraint, picket fences can be created for objects that satisfy

 $-1<\frac{\alpha_{a}^{\prime}}{\beta_{a}^{\prime}}<\sin^{2}\gamma.$ (15)

Under some conditions, including those depicted in Fig. 2, the direction of the force can alternate within the fringe pattern. In terms of the standard coupling coefficients, alternating picket fences can be projected for objects satisfying

 $-\frac{3}{4}<\frac{f_{0}}{f_{1}}<-\frac{1}{2}.$ (16)

Although this shows that alternating picket fences are only possible for a small domain of material properties, picket fences in general could have practical applications for sorting objects by density or compressibility.

### III.2 Conservative forces in standing waves

The acoustokinetic framework also is useful for analyzing the force fields created by specified sound waves. Standing waves, for example, can be decomposed into superpositions of counterpropagating plane waves. The phase-dependent terms in $\mathbf{F}_{\alpha}(\mathbf{r})$ and $\mathbf{F}_{\beta}(\mathbf{r})$ vanish in such superpositions, leaving a manifestly conservative force landscape,

 $\mathbf{F}_{\text{standing}}(\mathbf{r})=\frac{1}{4}\nabla\left[\alpha_{a}^{% \prime}u^{2}+\beta_{a}^{\prime}k^{-2}(\nabla u)^{2}\right].$ (17)

Particles satisfying

 $\alpha_{a}^{\prime}-\beta_{a}^{\prime}>0$ (18)

are drawn toward antinodes of the pressure field, as illustrated in Fig. 3. This condition is satisfied by compressible, low-density particles such as bubbles in water. Particles with complementary properties are drawn toward nodes.

### III.3 Nonconservative forces in pseudo-standing waves

Not all zero-momentum waves are standing waves. Some have non-trivial phase profiles and so can exert nonconservative forces. The archetype for such pseudo-standing waves is a superposition of three plane waves with equal amplitude, $p_{0}/3$, and wave vectors

 $\mathbf{k}_{n}=-k\left[\cos\left(n\frac{2\pi}{3}\right)\,\hat{x}+\sin\left(n% \frac{2\pi}{3}\right)\,\hat{y}\right]$ (19)

that satisfy $\sum_{n=1}^{3}\mathbf{k}_{n}=0$ [40]. The pressure intensity for such a three-wave superposition is plotted in Fig. 4(a) and displays a six-fold rotational symmetry similar to that of the corresponding six-fold standing wave in Fig. 3(b). The triangular lattice of antinodes in Fig. 4(a), however, is meshed with a dual hexagonal lattice of nodes, as indicated by dotted circles in Fig. 4(a). Because nonconservative forces tend to be weaker than conservative forces by a factor of $(ka_{p})^{3}$, we focus our attention on regions where the sound field forms stable traps, and expand about these points in polar coordinates, $\mathbf{r}=(r,\theta)$ for small displacements, $kr<1$.

For the three-fold pseudo-standing wave, the conservative contributions from Eqs. (9a) and (9b) simplify to

 $\displaystyle\mathbf{F}^{c}(\mathbf{r})=\frac{1}{8}\left(2\alpha_{a}^{\prime}-% \beta_{a}^{\prime}\right)\nabla u^{2}.$ (20)

Particles satisfying

 $2\alpha_{a}^{\prime}-\beta_{a}^{\prime}>0$ (21)

therefore are drawn to pressure antinodes while complementary particles seek out nodes.

A node-seeking particle experiences amplitude and phase profiles of the form

 $\displaystyle u_{\text{node}}(\mathbf{r})$ $\displaystyle\approx\frac{1}{2}p_{0}\,kr\quad\text{and}$ (22a) $\displaystyle\varphi_{\text{node}}(\mathbf{r})$ $\displaystyle\approx\pm\theta-\frac{\pi}{2}.$ (22b) The conservative part of the associated force field exerts a Hookean restoring force, $\mathbf{F}_{\text{node}}^{c}(\mathbf{r})\approx-\frac{1}{16}kp_{0}^{2}(\beta_{% a}^{\prime}-2\alpha_{a}^{\prime})\,kr\,\hat{r},$ (22c) that keeps the particle localized near the node. The force field also includes a nonconservative component, $\displaystyle\mathbf{F}^{nc}_{\text{node}}(\mathbf{r})\approx\pm\frac{1}{8}kp_% {0}^{2}(\alpha_{a}^{\prime\prime}+\beta_{a}^{\prime\prime})\,kr\,\hat{\theta},$ (22d)

that is directed by the azimuthal phase gradient and causes the displaced particle to orbit its node. The pseudo-standing wave therefore transfers orbital angular momentum to particles moving near its nodes, with each node acting as a unit-charge acoustic vortex [1]. The sign of the orbital angular momentum alternates from site to site on the honeycomb lattice of nodes. The array of alternating acoustic vortexes in a pseudo-standing wave therefore carries no net angular momentum [41].

Nodes repel particles satisfying $2\alpha_{a}^{\prime}>\beta_{a}^{\prime}$, which instead seek out the triangular lattice of antinodes. Near an antinode, the sound field’s amplitude and phase profiles are

 $\displaystyle u_{\text{antinode}}(\mathbf{r})$ $\displaystyle\approx p_{0}\left[1-\frac{1}{4}(kr)^{2}\right]\quad\text{and}$ (23a) $\displaystyle\varphi_{\text{antinode}}(\mathbf{r})$ $\displaystyle\approx\frac{1}{24}(kr)^{3}\,\cos 3\theta.$ (23b) For small displacements, the conservative terms from $\mathbf{F}_{\alpha}(\mathbf{r})$ and $\mathbf{F}_{\beta}(\mathbf{r})$ exert a Hookean restoring force on an antinode-seeking particle: $\mathbf{F}_{\text{antinode}}^{c}(\mathbf{r})=-\frac{1}{8}kp_{0}^{2}(2\alpha_{a% }^{\prime}-\beta_{a}^{\prime})\,kr\,\hat{r}.$ (23c) The nonconservative terms create a sextupole flow that tends to drive the particle from one antinode to another: $\mathbf{F}_{\text{antinode}}^{nc}(\mathbf{r})=F_{0}\,k^{2}r^{2}(\cos 3\theta\,% \hat{r}-\sin 3\theta\,\hat{\theta}),$ (23d)

where $F_{0}=kp_{0}^{2}\,(2\alpha_{a}^{\prime\prime}-\beta_{a}^{\prime\prime})/32$.

The conservative part of the pseudo-standing wave’s force field vanishes for materials satisfying $2\alpha_{a}^{\prime}=\beta_{a}^{\prime}$, leaving a purely nonconservative force field. This contrasts with the force exerted by a standing wave, which is always conservative. Realizing this condition in practice, however, would require a precise balance of material properties.

More generally, systems satisfying $2\alpha_{a}^{\prime\prime}-\beta^{\prime\prime}_{a}>2\alpha_{a}^{\prime}-\beta% _{a}^{\prime}$ will experience nonconservative forces that rival conservative trapping forces. In the particular case of an air bubble of size $ka_{p}=0.3$ in water, for example, the nonconservative force exceeds the conservative force by three orders of magnitude. The force field for this case thus resembles the streamlines in Fig. 4.

These observations illustrate that nonconservative forces can emerge along directions where the sound field carries no net momentum. More generally, it shows that radiation pressure can be directed independently of the direction of wave propagation. This independence can be used to craft tractor beams from propagation-invariant Bessel beams [42, 30, 18].

### III.4 Bessel Beams and Tractor Beams

Both the standing wave and the pseudo-standing wave require boundary conditions that completely enclose their targets. Long-range manipulation without physical confinement is best achieved with propagation-invariant traveling waves. The natural basis for such single-beam long-range manipulators is the family of Bessel beams [43, 44, 42], which are non-diffracting solutions to Eq. (5) in cylindrical coordinates, $\mathbf{r}=(r,\theta,z)$. The amplitude and phase profiles for a Bessel beam propagating along $\hat{z}$ are

 $\displaystyle u_{\gamma,n}(\mathbf{r})$ $\displaystyle=p_{0}\,J_{n}(kr\sin\gamma)\quad\text{and}$ (24) $\displaystyle\varphi_{\gamma,n}(\mathbf{r})$ $\displaystyle=kz\cos\gamma+n\theta,$ (25)

respectively, where $J_{n}(\cdot)$ is a Bessel function of the first kind of order $n$. Bessel beams are distinguished by the convergence angle $\gamma$ that ranges from $\gamma=0$ for conventional plane waves to $\gamma=\pi/2$ for circular standing waves, and the integer $n$ that imposes a helical pitch on the beam’s wavefronts.

The conservative part of a Bessel beam’s force field is directed radially. Beams with $n=0$ have maximum intensity along the axis, $r=0$. Those with $n>0$ have zero intensity on the axis. The radial component of the acoustic force is linear in $kr$ for $\left|n\right|\leq 2$ and scales as $(kr)^{3}$, or higher, for $\left|n\right|>2$. Whether the force attracts the particle to the axis or repels it depends on the particle’s properties. For simplicity, we restrict our analysis to $\left|n\right|\leq 2$, in which case the conditions for particles to be trapped on-axis are

 $\displaystyle 4\alpha_{a}^{\prime}+\beta_{a}^{\prime}\,(1+3\cos 2\gamma)$ $\displaystyle>0,$ $\displaystyle n=0$ (26a) $\displaystyle\alpha_{a}^{\prime}+\beta_{a}^{\prime}\,\cos 2\gamma$ $\displaystyle<0,$ $\displaystyle n=\pm 1$ (26b) $\displaystyle\beta_{a}^{\prime}$ $\displaystyle<0,$ $\displaystyle n=\pm 2.$ (26c)

As noted recently [45], the condition for stable trapping by an $n=0$ Bessel beam differs qualitatively from the analogous condition in Eq. (18) for trapping at an antinode of a standing wave. For example, dense objects with large values of $\beta_{a}^{\prime}$ are repelled by the antinodes of standing waves and pseudo-standing waves, but tend to be trapped by the central antinode of a Bessel beam with $\gamma<\pi/4$. Similar reversals arise for trapping at the central node of Bessel beams with $\left|n\right|>0$.

Having established the condition for stable transverse trapping, we next analyze the axial force on a particle localized at $r=0$. Any beam satisfying $\mathbf{F}(\mathbf{r})|_{r=0}\cdot\hat{z}<0$ can be said to act as a tractor beam. It should be noted that optical Bessel beams do not act as tractor beams for small objects because the dipole-order photokinetic force is always repulsive [18].

Axial forces in propagation-invariant Bessel beams are inherently nonconservative. The relevant terms in Eq. (9) yield

 $\displaystyle\mathbf{F}(\mathbf{r})|_{\begin{subarray}{c}n=0\\ r=0\end{subarray}}$ $\displaystyle=\frac{1}{2}\left[(\alpha_{a}^{\prime\prime}+\beta_{a}^{\prime% \prime})\,u^{2}+\frac{\beta_{a}^{\prime\prime}}{2k^{2}}\nabla^{2}u^{2}\right]% \,\nabla\varphi$ (27a) $\displaystyle=\frac{1}{2}\,p_{0}^{2}k\,(\alpha_{a}^{\prime\prime}+\beta_{a}^{% \prime\prime}\,\cos^{2}\gamma)\,\cos\gamma\,\hat{z}.$ (27b)

Unlike optical Bessel beams, therefore, acoustic Bessel beams can act as tractor beams for particles satisfying both the trapping condition from Eq. (26a) and also

 $\alpha_{a}^{\prime\prime}+\beta_{a}^{\prime\prime}\,\cos^{2}\gamma<0.$ (28)

Expressed in terms of material properties, these conditions simplify to

 $\displaystyle f_{0}<\frac{3}{8}(1+3\cos 2\gamma)f_{1}\quad\text{and}$ (29a) $\displaystyle\left(\frac{f_{0}}{f_{1}}\right)^{2}+\frac{f_{0}}{f_{1}}+\frac{3}% {4}\cos^{2}\gamma<0.$ (29b)

Figure 5 shows the domain of beam shapes and particle compositions for which the $n=0$ Bessel beam acts as a tractor beam. These include the optimal condition $f_{0}/f_{1}=-1/2$ that was identified in the original discussion of acoustic tractor beams [42].

Density-matched objects ($f_{1}=0$) can only be trapped if they are compressible enough that $f_{0}<0$. This means, however, that $\mathbf{F}\cdot\hat{z}>0$, from which we conclude that Bessel beams are not tractor beams for such objects.

The analogous treatment for $n=1$ yields

 $\displaystyle\mathbf{F}(\mathbf{r})|_{\begin{subarray}{c}n=1\\ r=0\end{subarray}}$ $\displaystyle=\frac{\beta_{a}^{\prime\prime}}{4k^{2}}\,\nabla^{2}u^{2}\,\nabla\varphi$ (30a) $\displaystyle=\frac{1}{4}p_{0}^{2}k\,\beta_{a}^{\prime\prime}\,\cos\gamma\sin^% {2}\gamma\,\hat{z}.$ (30b)

Because $\beta_{a}^{\prime\prime}>0$, such beams do act as tractor beams, at least not for objects trapped on the axis. They instead drive trapped objects downstream.

The $n=2$ beam can trap small objects on the axis, but exerts no axial force at all,

 $\mathbf{F}(\mathbf{r})|_{\begin{subarray}{c}n=2\\ r=0\end{subarray}}=0.$ (31)

Such beams might serve as useful conduits for Rayleigh particles that are moved back and forth along the axis by other forces.

## IV Discussion

The theory of acoustokinetic forces presented here expresses the influence of a sound wave on a small object in terms of the amplitude and phase profiles of the pressure field. The resulting analytical framework, which is summarized in Eq. (9), offers a complementary perspective on acoustic forces to the standard development, which explicitly refers to the velocity field. Acoustokinetic expressions are particularly useful for designing acoustic force fields because they inherently account for coupling between the pressure and velocity fields and specify amplitude and phase profiles that can be controlled with transmissive, reflective or emissive elements.

The acoustokinetic approach yields expressions that closely resemble results for optical forces in the Rayleigh limit. Exploring the differences and similarities between acoustokinetic and photokinetic forces offers insights into how sound fields couple to material properties to generate useful and interesting force landscapes. The most important difference is that the particle’s polarizability includes contributions from the incident wave’s monopole component. This means that the force arising from the sound field’s dipole component includes terms at the same order of magnitude as leading quadrupole contributions. Such mixing does not arise in optical forces. It vanishes for objects that are density matched with the medium, yielding expressions for acoustic forces that are exactly analogous to their optical counterparts.

The acoustokinetic framework naturally accounts for the conservative nature of the force field exerted by standing waves and differentiates it from the influence of pseudo-standing waves that also carry no net momentum yet still exert nonconservative forces. When applied to acoustic Bessel beams, the acoustokinetic framework provides clear guidance on the limits for axial trapping and the conditions under which a zeroth-order Bessel beam can act as a tractor beam. This result contrasts with the equivalent analysis of optical Bessel beams, which do not to act as tractor beams in the Rayleigh regime.

These examples illustrate the value of the acoustokinetic framework for analyzing acoustic force fields. More generally, the expressions from Eq. (9) can be solved either numerically, or in some cases analytically, for the phase and amplitude profiles that correspond to a specified force field acting on a particle of a given size, compressibility and density. As in the optical case, the framework can be extended beyond the dipole approximation. For objects substantially smaller than the wavelength of sound, however, the dipole-order terms are simple enough to yield analytical expressions. This framework also can be extended to handle interactions among insonated objects, an effect usually referred to as secondary Bjerknes forces.

## V Acknowledgments

This work was supported by the MRSEC program of the National Science Foundation through Award Number DMR-1420073. The authors acknowledge helpful conversations with Philip Marston (Washington State University) and Konstantin Bliokh (The Australian National University and RIKEN, Japan).

## References

• [1] L. Zhang and P. L. Marston. “Angular momentum flux of nonparaxial acoustic vortex beams and torques on axisymmetric objects.” Phys. Rev. E 84, 065601 (2011).
• [2] W. Wang, L. A. Castro, M. Hoyos and T. E. Mallouk. “Autonomous motion of metallic microrods propelled by ultrasound.” ACS Nano 6, 6122–6132 (2012).
• [3] A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long and S. Subramanian. ‘‘Holographic acoustic elements for manipulation of levitated objects.” Nat. Commun. 6, 8661 (2015).
• [4] A. Marzo and B. W. Drinkwater. “Holographic acoustic tweezers.” Proc. Natl. Acad. Sci. U.S.A. 116, 84–89 (2019).
• [5] W. J. Xie, C. D. Cao, Y. J. Lü, Z. Y. Hong and B. Wei. “Acoustic method for levitation of small living animals.” Appl. Phys. Lett. 89, 214102 (2006).
• [6] H.-L. Cao, D.-C. Yin, Y.-Z. Guo, X.-L. Ma, J. He, W.-H. Guo, X.-Z. Xie and B.-R. Zhou. “Rapid crystallization from acoustically levitated droplets.” J. Acoust. Soc. Am. 131, 3164–3172 (2012).
• [7] M. A. B. Andrade, A. L. Bernassau and J. C. Adamowski. “Acoustic levitation of a large solid sphere.” Appl. Phys. Lett. 109, 044101 (2016).
• [8] A. Lenshof, C. Magnusson and T. Laurell. “Acoustofluidics 8: Applications of acoustophoresis in continuous flow microsystems.” Lab Chip 12, 1210–1223 (2012).
• [9] C. R. Courtney, B. W. Drinkwater, C. E. Demore, S. Cochran, A. Grinenko and P. D. Wilcox. “Dexterous manipulation of microparticles using Bessel-function acoustic pressure fields.” Appl. Phys. Lett. 102, 123508 (2013).
• [10] V. Bjerknes. Fields of Force (Columbia University Press, 1906).
• [11] F. Blake Jr. “Bjerknes forces in stationary sound fields.” J. Acoust. Soc. Am. 21, 551–551 (1949).
• [12] L. P. Gor’kov. “On the forces acting on a small particle in an acoustical field in an ideal fluid.” Sov. Phys. Dokl. 6, 773–775 (1962).
• [13] M. Settnes and H. Bruus. “Forces acting on a small particle in an acoustical field in a viscous fluid.” Phys. Rev. E 85, 016327 (2012).
• [14] G. T. Silva. “Acoustic radiation force and torque on an absorbing compressible particle in an inviscid fluid.” J. Acoust. Soc. Am. 136, 2405–2413 (2014).
• [15] J. P. Gordon. “Radiation forces and momenta in dielectric media.” Phys. Rev. A 8, 14–21 (1973).
• [16] 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).
• [17] 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).
• [18] A. Yevick, D. B. Ruffner and D. G. Grier. “Tractor beams in the Rayleigh limit.” Phys. Rev. A 93, 043807 (2016).
• [19] A. Yevick, D. J. Evans and D. G. Grier. “Photokinetic analysis of the forces and torques exerted by optical tweezers carrying angular momentum.” Phil. Trans. Roy. Soc. A 375, 20150432 (2017).
• [20] 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).
• [21] 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).
• [22] 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).
• [23] B. Sun, J. Lin, E. Darby, A. Y. Grosberg and D. G. Grier. “Brownian vortexes.” Phys. Rev. E 80, 010401(R) (2009).
• [24] 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).
• [25] H. He, M. E. J. Friese, N. R. Heckenberg and H. Rubinsztein-Dunlop. “Direct observation of transfer of angular momentum to absorptive particles from a laser beam with a phase singularity.” Phys. Rev. Lett. 75, 826–829 (1995).
• [26] K. T. Gahagan and G. A. Swartzlander. “Optical vortex trapping of particles.” Opt. Lett. 21, 827–829 (1996).
• [27] N. B. Simpson, L. Allen and M. J. Padgett. “Optical tweezers and optical spanners with Laguerre-Gaussian modes.” J. Mod. Opt. 43, 2485–2491 (1996).
• [28] K. Ladavac, K. Kasza and D. G. Grier. “Sorting by periodic potential energy landscapes: Optical fractionation.” Phys. Rev. E 70, 010901(R) (2004).
• [29] K. Ladavac and D. G. Grier. “Colloidal hydrodynamic coupling in concentric optical vortices.” Europhys. Lett. 70, 548–554 (2005).
• [30] S.-H. Lee, Y. Roichman and D. G. Grier. “Optical solenoid beams.” Opt. Express 18, 6988–6993 (2010).
• [31] B. T. Draine. “The discrete-dipole approximation and its application to interstellar graphite grains.” Astrophys. J. 333, 848–872 (1988).
• [32] S. Albaladejo, R. Gómez-Medina, L. S. Froufe-Pérez, H. Marinchio, R. Carminati, J. F. Torrado, G. Armelles, A. García-Martin and J. J. Sáenz. “Radiative corrections to the polarizability tensor of an electrically small anisotropic dielectric particle.” Opt. Express 18, 3556–3567 (2010).
• [33] M. X. Lim, K. A. Murphy and H. M. Jaeger. “Edges control clustering in levitated granular matter.” arXiv preprint arXiv:1901.08656 (2019).
• [34] H. E. Bass, L. C. Sutherland, A. J. Zuckerwar, D. T. Blackstock and D. M. Hester. “Atmospheric absorption of sound: Further developments.” J. Acoust. Soc. Am. 97, 680–683 (1995).
• [35] J. Lighthill. “Acoustic streaming.” J. Sound Vibration 61, 391–418 (1978).
• [36] I. D. Toftul, K. Y. Bliokh, M. I. Petrov and F. Nori. “Acoustic Radiation Force and Torque on Small Particles as Measures of the Canonical Momentum and Spin Densities.” arXiv preprint arXiv:1905.12216 (2019).
• [37] D. Bazou, A. Castro and M. Hoyos. “Controlled cell aggregation in a pulsed acoustic field.” Ultrasonics 52, 842–850 (2012).
• [38] K. Melde, A. G. Mark, T. Qiu and P. Fischer. “Holograms for acoustics.” Nature 537, 518 (2016).
• [39] C. E. M. Démoré, P. M. Dahl, Z. Yang, P. Glynne-Jones, A. Melzer, S. Cochran, M. P. MacDonald and G. C. Spalding. “Acoustic tractor beam.” Phys. Rev. Lett. 112, 174302 (2014).
• [40] J. Masajada and B. Dubik. “Optical vortex generation by three plane wave interference.” Opt. Commun. 198, 21–27 (2001).
• [41] K. Y. Bliokh and F. Nori. “Spin and orbital angular momenta of acoustic beams.” Phys. Rev. B 99, 174310 (2019).
• [42] 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).
• [43] J. Durnin. ‘‘Exact-solutions for nondiffracting beams. 1. The scalar theory.” J. Opt. Soc. Am. A 4, 651–654 (1987).
• [44] J. Durnin, J. J. Miceli, Jr and J. H. Eberly. “Diffraction-free beams.” Phys. Rev. Lett. 58, 1499–1501 (1987).
• [45] X.-D. Fan and L. Zhang. “Trapping force of acoustical Bessel beams on a sphere and stable tractor beams.” Phys. Rev. Appl. 11, 014055 (2019).