A method for generating stable ultrasonic levitation of physical matter in air using single beams (also known as tractor beams) is demonstrated. The method encodes the required phase modulation in passive unit cells into which the ultrasonic sources are mounted. These unit cells use waveguides such as straight and coiled tubes to act as delay-lines. It is shown that a static tractor beam can be generated using a single electrical driving signal, and a tractor beam with one-dimensional movement along the propagation direction can be created with two signals. Acoustic tractor beams capable of holding millimeter-sized polymer particles of density 1.25 g/cm 3 and fruit-flies (Drosophila) are demonstrated. Based on these design concepts, we show that portable tractor beams can be constructed with simple components that are readily available and easily assembled, enabling applications in industrial contactless manipulation and biophysics.

generate acoustic radiation forces. 1–4 147(861), 212 (1934). 1. L. V. King, Proc. R. Soc. London(861), 212 (1934). https://doi.org/10.1098/rspa.1934.0215 2. L. P. Gor'Kov, Sov. Phys. Doklady 6, 773 (1962). 100(2), 1231 (1996). 3. A. A. Doinikov, J. Acous. Soc. Am.(2), 1231 (1996). https://doi.org/10.1121/1.415969 12(6), 1014 (2012). 4. H. Bruus, Lab Chip(6), 1014 (2012). https://doi.org/10.1039/c2lc21068a 5 413(6855), 474 (2001). 5. E. H. Brandt, Nature(6855), 474 (2001). https://doi.org/10.1038/35097192 6 109(28), 11105 (2012). 6. X. Ding, S. C. S. Lin, B. Kiraly, H. Yue, S. Li, I. K. Chiang, and T. J. Huang, Proc. Natl. Acad. Sci.(28), 11105 (2012). https://doi.org/10.1073/pnas.1209288109 7 78(8), 4845 (1995). 7. H. M. Hertz, J. Appl. Phys.(8), 4845 (1995). https://doi.org/10.1063/1.359770 8 41(4), 397 (2012). 8. R. J. Weber, C. J. Benmore, S. K. Tumber, A. N. Tailor, C. A. Rey, L. S., Taylor, and S. R. Byrn, Eur. Biophys. J.(4), 397 (2012). https://doi.org/10.1007/s00249-011-0767-3 9–11 7(9), 1125 (2007). 9. L. Puskar, R. Tuckermann, T. Frosch, J. Popp, V. Ly, D. McNaughton, and B. R. Wood, Lab Chip(9), 1125 (2007). https://doi.org/10.1039/b706997a 6, 20023 (2016). 10. T. Vasileiou, D. Foresti, A. Bayram, D. Poulikakos, and A. Ferrari, Sci. Rep., 20023 (2016). https://doi.org/10.1038/srep20023 5, 13596 (2015). 11. M. Sundvik, H. J. Nieminen, A. Salmi, P. Panula, and E. Hæggström, Sci. Rep., 13596 (2015). https://doi.org/10.1038/srep13596 12 89(21), 214102 (2006). 12. W. J. Xie, C. D. Cao, Y. J. Lü, Z. Y. Hong, and B. Wei, Appl. Phys. Lett.(21), 214102 (2006). https://doi.org/10.1063/1.2396893 Sound is a mechanical wave and as such it carries momentum that canacoustic radiation forces.When these forces are strong enough and converge from all directions, particles can be levitated against gravity.Acoustic levitation is becoming a fundamental tool in lab-on-a-chip scenarios,microscopy,pharmaceuticals,the levitation of biological samples,and even small animals.

generated using phased arrays. 13–15 112(17), 174302 (2014). 13. C. E. Démoré, P. M. Dahl, Z. Yang, P. Glynne-Jones, A. Melzer, S. Cochran, and G. C. Spalding, Phys. Rev. Lett.(17), 174302 (2014). https://doi.org/10.1103/PhysRevLett.112.174302 6, 8661 (2015). 14. A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long, and S. Subramanian, Nat. Commun., 8661 (2015). https://doi.org/10.1038/ncomms9661 116(2), 024301 (2016). 15. D. Baresch, J. L. Thomas, and R. Marchiano, Phys. Rev. Lett.(2), 024301 (2016). https://doi.org/10.1103/PhysRevLett.116.024301 tractor beams. Acoustic tractor beams have the potential to revolutionize contactless manipulation due to their high exerted force to input power ratio 15 116(2), 024301 (2016). 15. D. Baresch, J. L. Thomas, and R. Marchiano, Phys. Rev. Lett.(2), 024301 (2016). https://doi.org/10.1103/PhysRevLett.116.024301 phased array systems are currently necessary in order to generate and amplify dozens of independent electric signals making tractor beams not generally accessible due to cost, space, or complexity of operation. Recently, single-beam acoustic levitators have beenusingWe adopt terminology common in optics and term these single-beams that can trap objects in three dimensions. Acoustic tractor beams have the potential to revolutionize contactless manipulation due to their high exerted force to input power ratioand the wide variety of supported particle materials as well as sizes. However,systems are currently necessary in order toand amplify dozens of independent electric signals making tractor beams not generally accessible due to cost, space, or complexity of operation.

Phased arrays are a collection of elements that emit or receive with specific phases or time delays. They are in widespread use in radar, 16 16. A. J. Fenn, D. H. Temme, W. P. Delaney, and W. E. Courtney, Lincoln Lab. J. 12(2), 321 (2000). 17 17. J. H. G. Ender and A. R. Brenner, IEEE Proc.-Radar, Sonar Navig. 150(3), 165 (2003). 18 39(7), 525 (2006). 18. B. W. Drinkwater and P. D. Wilcox, NDT&E Int.(7), 525 (2006). https://doi.org/10.1016/j.ndteint.2006.03.006 19 19. T. L. Gereld, U.S. patent 2,387,783 (30 October 1945). are a collection of elements that emit or receive with specific phases or time delays. They are in widespread use in radar,sonar,and ultrasonic imagingsince they can dynamically change the direction or shape of the beam. At present, the phase delays are programmable and achieved using bespoke electronics, but when these systems were first developed the delays were static and introduced with passive elements. For instance, electrical signals were made to travel through long cables to introduce a time delay proportional to its length.

generate an acoustic tractor beam. The sonic devices are composed of several unit cells, and the emitters are mounted into them. Each unit cell introduces a specific phase delay. Our approach parallels developments in 2D planar metamaterials 20–25 108(11), 114301 (2012). 20. Z. Liang and J. Li, Phys. Rev. Lett.(11), 114301 (2012). https://doi.org/10.1103/PhysRevLett.108.114301 3, 2546 (2013). 21. Y. Li, B. Liang, Z. M. Gu, X. Y. Zou, and J. C. Cheng, Sci. Rep., 2546 (2013). https://doi.org/10.1038/srep02546 2(6), 064002 (2014). 22. Y. Li, X. Jiang, R. Q. Li, B. Liang, X. Y. Zou, L. L. Yin, and J. C. Cheng, Phys. Rev. Appl.(6), 064002 (2014). https://doi.org/10.1103/PhysRevApplied.2.064002 4(2), 024003 (2015). 23. Y. Li, X. Jiang, B. Liang, J. C. Cheng, and L. Zhang, Phys. Rev. Appl.(2), 024003 (2015). https://doi.org/10.1103/PhysRevApplied.4.024003 117(2), 024503 (2015). 24. K. Tang, C. Qiu, J. Lu, M. Ke, and Z. Liu, J. Appl. Phys.(2), 024503 (2015). https://doi.org/10.1063/1.4905910 117(3), 034301 (2016). 25. X. Jiang, Y. Li, B. Liang, J. C. Cheng, and L. Zhang, Phys. Rev. Lett.(3), 034301 (2016). https://doi.org/10.1103/PhysRevLett.117.034301 Inspired by these early passive delay lines, we developed devices that modulate a simple wavefront in order toan acoustic tractor beam. The sonic devices are composed of several unit cells, and the emitters are mounted into them. Each unit cell introduces a specific phase delay. Our approach parallels developments in 2D planar metamaterialswhere a layer of sub-wavelength unit-cells is used to shape and steer incident fields. In our sonic devices, the phase-shifting unit cells are not limited to the size of the wavelength and the emitters are designed as part of the structure. This approach minimises operational cost and complexity, making tractor beams portable and more readily deployable, hence enabling applications in industrial contactless manipulation and biophysics.

P at point r due to a collection of N transducers emitting through a layer of passive delay lines can be written as P ( r ) = ∑ n = 1 N P n L n D f ( θ n ) d n e i ( φ n + k Δ n ) , (1) where P n is a constant that defines the transducer output power, L n is a loss factor associated with each delay line, and D f is a far-field directivity function that depends on the angle θ n between the effective source normal and r. Here, D f = 2 J 1 ( k a sin θ n ) / k a sin θ n , the directivity function of a piston source, where J 1 is a first order Bessel function of the first kind and a is the radius of the piston. The term 1 / d n accounts for divergence, where d n is the propagation distance in free space. The phase of the wave at point r is φ n + k Δ n , where φ n is the applied electronic phase, Δ n is the propagation distance between r and the emitters, k = 2 π / λ is the wavenumber, and λ is the wavelength. The acoustic pressureat pointdue to a collection ofemitting through a layer of passive delay lines can be written aswhereis a constant that defines theoutput power,is a loss factor associated with each delay line, andis a far-field directivity function that depends on the anglebetween the effective source normal and. Here,, the directivity function of a piston source, whereis a first order Bessel function of the first kind andis the radius of the piston. The termaccounts for divergence, whereis the propagation distance in free space. The phase of the wave at pointis, whereis the applied electronic phase,is the propagation distance betweenand the emitters,is the wavenumber, andis the wavelength.

r can be controlled either by the electronic phase modulation ( φ n ) or by the distance that the wave travels ( Δ n ). Fig. 1 phased array into passive sonic devices. In the electronic phased array, the circuitry controls φ n whereas in the passive devices φ n is a constant for all the elements and the phase is controlled by changing Δ n . It should be noted that the phase of a wave arriving atcan be controlled either by the electronic) or by the distance that the wave travels (). Fig.(Multimedia view) illustrates the three concepts that we used to encode the phases of an electronicinto passive sonic devices. In the electronicthe circuitry controlswhereas in the passive devicesis a constant for all the elements and the phase is controlled by changing

Electronic phased array. An electronic phased array in which the elements are all located in a horizontal plane (z = 0) is shown in Fig. 1(a) Δ n = d n = | r − s n | and θ n = A z m ∠ ( r − s n ) , where s n is the position of a source and A z m ∠ denotes the azimuthal angle. Typically, the element outputs are uniform (i.e., P n = P o ) and there is no insertion loss ( L n = 1 ). For instance, to focus the array at a point f = ( 0 , 0 , f l ) , the phase modulation is found from simple time reversal as φ n = k | f − s n | . An electronicin which the elements are all located in a horizontal plane (0) is shown in Fig.. Here,and, whereis the position of a source anddenotes the azimuthal angle. Typically, the element outputs are uniform (i.e.,) and there is no insertion loss (). For instance, to focus the array at a point, theis found from simple time reversal as

Straight tubes. In Fig. 1(b) Δ n = d n + h n , where d n is the propagation distance in free space (i.e., once the sound waves have exited the tubes) and h n is the tube length. If t n is the position of the tube exit then: d n = | r − t n | , θ n = A z m ∠ ( r − t n ) , and h n = | t n − s n | . The tube length for introducing a phase modulation of φ is h = φ / k ; this length can be seen as a static alternative to the electronic delay. The tubes have the same diameter as the transducers to facilitate their attachment; wider (or narrower) tubes may cause an increase in the insertion loss and change the linear relation between phase modulation and length. 26 39(1), 1 (1975). 26. H. Tijdeman, J. Sound Vib.(1), 1 (1975). https://doi.org/10.1016/S0022-460X(75)80206-9 h n modulo λ ) where modulo operation gives the remainder of the integer division. In Fig., the emitters are attached below vertical rigid tubes of different heights; this sonic device has a flat upper surface. In Eq. (1) , whereis the propagation distance in free space (i.e., once the sound waves have exited the tubes) andis the tube length. Ifis the position of the tube exit then:, and. The tube length for introducing aofis; this length can be seen as a static alternative to the electronic delay. The tubes have the same diameter as theto facilitate their attachment; wider (or narrower) tubes may cause an increase in the insertion loss and change the linear relation betweenand length.Since continuous wave emission was used, the lengths of the tubes were trimmed according to () where modulo operation gives the remainder of the integer division.

Sculpted surface. The tubes and emitters from the previous device have been aligned with the lines that join the target point and the top of the tubes as illustrated in Fig. 1(c) Δ n remains the same even when the tubes are removed. The emitters are held in place by slots in a shaped surface, which also acts as a supporting structure. In Eq. θ n and L n are calculated as in the electronic phased array. The emitters are pointing at the focal point so θ n = 0 and hence D f = 1 , which means that the maximum intensity of the transducers is directed towards the target making this the most efficient possible configuration. The tubes and emitters from the previous device have been aligned with the lines that join the target point and the top of the tubes as illustrated in Fig.. Therefore, the propagation distanceremains the same even when the tubes are removed. The emitters are held in place by slots in a shaped surface, which also acts as a supporting structure. In Eq. (1) andare calculated as in the electronicThe emitters are pointing at the focal point soand hence, which means that the maximum intensity of theis directed towards the target making this the most efficient possible configuration.

Coils. In Fig. 1(d) 27,28 5, 5553 (2014). 27. Y. Xie, W. Wang, H. Chen, A. Konneker, B. I. Popa, and S. A. Cummer, Nat. Commun., 5553 (2014). https://doi.org/10.1038/ncomms6553 7, 11731 (2016). 28. X. Zhu, K. Li, P. Zhang, J. Zhu, J. Zhang, C. Tian, and S. Liu, Nat. Commun., 11731 (2016). https://doi.org/10.1038/ncomms11731 Δ n = d n + c n , where c n is the length of the coil, note that the heights of the coils ( h n ) are equal and that c n depends on the coil revolutions. Since the coiling requires extra space inside the unit cell, the coil diameter is reduced in the inner region. We used horns 29 192(11), 1533 (2003). 29. E. Bängtsson, D. Noreland, and M. Berggren, Comput. Methods Appl. Mech. Eng.(11), 1533 (2003). https://doi.org/10.1016/S0045-7825(02)00656-4 coil to transition between the diameter at the entrance and exit (transducer diameter) and the inner reduced diameter. The model for the coil follows closely that of the tube for d n and θ n . In Fig., the tubes are coiled to fit into a flat rectangular slab. This sonic device was stimulated by the active research in the field of space-coiling metamaterials in which conduits are coiled in different ways.In Eq. (1) , whereis the length of thenote that the heights of the) are equal and thatdepends on therevolutions. Since the coiling requires extra space inside the unit cell, thediameter is reduced in the inner region. We used hornsat the entrance and exit of theto transition between the diameter at the entrance and exitdiameter) and the inner reduced diameter. The model for thefollows closely that of the tube forand

generated by focusing the array at a target position and then adding a static phase signature that only depends on the type of trap to be generated; three distinct phase signatures have been shown to generate functional tractor beams. 14 6, 8661 (2015). 14. A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long, and S. Subramanian, Nat. Commun., 8661 (2015). https://doi.org/10.1038/ncomms9661 generated; and if a helical phase ramp is used, then a Vortex trap results. We focused our designs on the Twin trap since the Vortex trap can only seize small particles without ejecting them and the Bottle trap has relatively weak lateral forces. 14 6, 8661 (2015). 14. A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long, and S. Subramanian, Nat. Commun., 8661 (2015). https://doi.org/10.1038/ncomms9661 generate any type of trap. An acoustic tractor beam isby focusing the array at a target position and then adding a static phase signature that only depends on the type of trap to bethree distinct phase signatures have been shown tofunctional tractor beams.Briefly, if half the array is excited out-of-phase, then a Twin trap results; if a central circular aperture is driven out-of-phase, then a Bottle trap isand if a helical phase ramp is used, then atrap results. We focused our designs on the Twin trap since thetrap can only seize small particles without ejecting them and the Bottle trap has relatively weak lateral forces.However, it should be noted that the presented devices could easily be altered toany type of trap.

To excite the out-of-phase elements needed for a Twin trap there are two options: introducing the π-phase offset in the delay lines by lengthening or shortening them by a distance of λ/2 or by reversing the polarity of the emitters. Introducing the signatures directly to the delay lines creates an integrated simpler structure but introduces small dissimilarities in the attenuation per element as different length tubes or coils are required. Switching the polarity requires a specific cable-transducer attachment and so is not such a flexible solution, but it leads to simple and more uniform surfaces. For the Straight tubes and Coils, we used the former method whereas for the Sculpted surface we used the latter leading to a simpler surface with all the emitters having the same distance to the target.

transducers operating at 40 kHz ( λ = 8.5 mm): 10 mm diameter transducers (MA40S4S) for the Straight Tubes and Coils, and 16 mm diameter transducers for the Sculpted Surface (MCUSD16P40B12RO). With this configuration, the Straight tubes had a directivity ( D f ) similar to that of the piston source (see supplementary Figs. 1 and 2), the linear relation ( h = φ / k ) was found to hold ( R 2 = 0.997 ) (see supplementary Fig. 3), and the loss factor was L n = 0.8 on average (see supplementary Fig. 4). For the coils, the relation between coil revolutions and phase modulation is presented in supplementary Fig. 5. The coils were 10 mm diameter at the entrance and exit and 2 mm in the inner region; this configuration reduced the insertion loss (see supplementary Fig. 6). The coil horns produced a directivity function similar to a piston source (see supplementary Figs. 1 and 2). The insertion loss for the coils depending on the revolutions and on average was L n = 0.7 (see supplementary Fig. 7). The devices developed to create a Twin trap using a single driving signal are illustrated in Fig. 2 For the practical realization of the devices, we usedoperating at 40 kHz (= 8.5 mm): 10 mm diameter(MA40S4S) for the Straight Tubes andand 16 mm diameterfor the Sculpted Surface (MCUSD16P40B12RO). With this configuration, the Straight tubes had a directivity () similar to that of the piston source (see supplementary Figs. 1 and 2), the linear relationwas found to hold () (see supplementary Fig. 3), and the loss factor was= 0.8 on average (see supplementary Fig. 4). For thethe relation betweenrevolutions andis presented in supplementary Fig. 5. Thewere 10 mm diameter at the entrance and exit and 2 mm in the inner region; this configuration reduced the insertion loss (see supplementary Fig. 6). Thehorns produced a directivity function similar to a piston source (see supplementary Figs. 1 and 2). The insertion loss for thedepending on the revolutions and on average was(see supplementary Fig. 7). The devices developed to create a Twin trap using a single driving signal are illustrated in Fig.(Multimedia view). These devices are portable as well as created with components that can be bought and put together by anyone (A step by step video is included in the supplementary material , and the source code and 3D files can be found online).

generate a tractor beam. In Figs. 3(a)–3(d) generated with a Coil device. Figs. 3(e)–3(h) Vortex trap generated with a Straight tubes device. The agreement between experiment and prediction is good, especially around the trapping point. Also, full three-dimensional Finite-Difference Time-Domain (FDTD) simulations were performed and found to be in good agreement with the simple model expressed by Eq. When the three passive devices are powered with a single driving signal, theya tractor beam. In Figs.(Multimedia view), we show the modelled and experimental acoustic pressure amplitude (see supplementary material ) for a Twin trapwith adevice. Figs.show atrapwith a Straight tubes device. The agreement between experiment and prediction is good, especially around the trapping point. Also, full three-dimensionalsimulations were performed and found to be in good agreement with the simple model expressed by Eq. (1) (see supplementary Fig. 8).

R p = 2 mm, which given a density ρ E P S = 29 kg/m3 represents a levitation force of 1 μN. The size limit is imposed as the devices levitate in the Rayleigh scattering regime in which the particle has to be small in comparison to the wavelength. The EPS particles were levitated vertically (i.e., above an array located in the x-y plane) and laterally (i.e., to the side of an array located in the x-z plane); in the case of the Sculpted surface devices also levitation upside down was achieved (i.e., below an array located in the x-y plane). Additionally, it was possible to levitate a particle of polymer (Polylactic Acid, ρ P L A = 1250 kg/m3, dia. 1.85 mm, thickness 0.5 mm) as well as fruit-flies (Drosophila, 0.3 mg each). A collection of these scenarios is shown in Fig. 4 The three devices were able to levitate expanded polystyrene (EPS) spheres of2 mm, which given a density29 kg/mrepresents a levitation force of 1N. The size limit is imposed as the devices levitate in the Rayleigh scattering regime in which the particle has to be small in comparison to the wavelength. The EPS particles were levitated vertically (i.e., above an array located in the x-y plane) and laterally (i.e., to the side of an array located in the x-z plane); in the case of the Sculpted surface devices also levitation upside down was achieved (i.e., below an array located in the x-y plane). Additionally, it was possible to levitate a particle of polymer (Polylactic Acid,1250 kg/m, dia. 1.85 mm, thickness 0.5 mm) as well as fruit-flies (Drosophila, 0.3 mg each). A collection of these scenarios is shown in Fig.(Multimedia view).

R p = 1 mm EPS sphere with a Twin trap 1.8 cm above the surface and finding the minimum required voltage. This ranked the Sculpted surface as the most effective device, followed by the Straight tubes and then the Coil device. In brief, a drive signal of 14Vpp (5.6 W) was required for levitation using the Sculpted surface, 34Vpp (27.2 W) for the Straight tubes, and 43Vpp (51.6 W) for the Coils. In the normal regime of operation, the pressure generated by an emitter is proportional to the voltage of the driving signal. The trapping stiffness is the Laplacian of the Gor'kov Potential, 14 6, 8661 (2015). 14. A. Marzo, S. A. Seah, B. W. Drinkwater, D. R. Sahoo, B. Long, and S. Subramanian, Nat. Commun., 8661 (2015). https://doi.org/10.1038/ncomms9661 generated by our devices can be expressed as K = ε V 2 , where ϵ is a coefficient that defines the device efficiency and V is the applied voltage. The ε was 2.35 for the sculpted surface, 0.39 for the Straight tubes, and 0.24 for the Coils. These differences can be explained by the relatively high loss factors of the Coils (on average L n = 0.7) and because the flat arrangement of transducers incurs losses due to directivity effects (i.e., the transducers are not directed to the focal point). Good agreement was found with Eq. ε = 2.25, 0.39 and 0.31, respectively. Device performance was measured experimentally by levitating a= 1 mm EPS sphere with a Twin trap 1.8 cm above the surface and finding the minimum required voltage. This ranked the Sculpted surface as the most effective device, followed by the Straight tubes and then thedevice. In brief, a drive signal of 14Vpp (5.6 W) was required for levitation using the Sculpted surface, 34Vpp (27.2 W) for the Straight tubes, and 43Vpp (51.6 W) for theIn the normal regime of operation, the pressureby an emitter is proportional to the voltage of the driving signal. The trapping stiffness is the Laplacian of the Gor'kov Potential,which is proportional to the square of the pressure. Therefore, the trapping stiffnessby our devices can be expressed as, where ϵ is a coefficient that defines the device efficiency andis the applied voltage. Thewas 2.35 for the sculpted surface, 0.39 for the Straight tubes, and 0.24 for theThese differences can be explained by the relatively high loss factors of the(on average= 0.7) and because the flat arrangement ofincurs losses due to directivity effects (i.e., theare not directed to the focal point). Good agreement was found with Eq. (1) (RMSD = 0.06), which gave predicted stiffness coefficients of2.25, 0.39 and 0.31, respectively.

transducers in a given row is the same; therefore, a single signal can be used to drive all the transducers in that row, swapping the polarity in the opposite half of that row to get a Twin trap. We then achieve dynamic trapping by following the technique used in concentric-ring phased arrays to refocus the beam at different heights. 30 34(5), 542 (1986). 30. C. A. Cain and S. I. Umemura, IEEE Trans. Microwave Theory Tech.(5), 542 (1986). https://doi.org/10.1109/TMTT.1986.1133390 2(b) transducers removed, leaving only the top and middle rows; the electronic circuit is presented in supplementary Fig. 9 and the arrangement in supplementary Fig. 10. To explore the functionality of this two channel arrangement, the particle was placed in the natural focus of the bowl (1.8 cm above the top) and the relative phase of the two 25Vpp excitation signals was changed to move the particle upwards in increments of 0.5 mm until the particle fell out of the trap under gravity. With this setup, the maximum height reached with an EPS bead Rp = 1 mm was 4.5 cm above the top of the spherical cap. Axisymmetric designs of the other sonic device concepts would enable them to use this feature. The three devices were excited using a single electrical drive signal and consequently they produced a static trap. However, it is possible to exploit geometric symmetries to achieve dynamic focusing of the trap using a reduced number of signals. We first note that in an axisymmetric arrangement the distance between any point on the revolution axis and thein a given row is the same; therefore, a single signal can be used to drive all thein that row, swapping the polarity in the opposite half of that row to get a Twin trap. We then achieve dynamic trapping by following the technique used in concentric-ringto refocus the beam at different heights.We used the Sculpted Surface device presented in Fig.with the bottom row ofremoved, leaving only the top and middle rows; the electronic circuit is presented in supplementary Fig. 9 and the arrangement in supplementary Fig. 10. To explore the functionality of this two channel arrangement, the particle was placed in the natural focus of the bowl (1.8 cm above the top) and the relative phase of the two 25Vpp excitation signals was changed to move the particle upwards in increments of 0.5 mm until the particle fell out of the trap under gravity. With this setup, the maximum height reached with an EPS bead Rp = 1 mm was 4.5 cm above the top of the spherical cap. Axisymmetric designs of the other sonic device concepts would enable them to use this feature.

Three types of passive devices capable of generating acoustic tractor beams have been demonstrated. They employ a small number of electric driving signals, making them a significant simplification over existing tractor beams which require electronic phased arrays with dozens of independent signals. The Sculpted surface sonic device is the most efficient device but Straight tubes and Coils have the advantage of presenting a flat emitting surface, making them more widely deployable. Additionally, Coil devices could use a single source since they are also flat on the lower surface. The ability of phased arrays to steer the beam has been replaced here in favour of a more portable, affordable, simple, and compact solution. This trade-off is desirable for scenarios that require levitation at fixed locations (which could potentially be coupled with mechanical motion) or limited movement with the addition of a small number of independent channels. The presented sonic devices have the potential to make acoustic tractor beams a widespread solution in contactless manipulation.

phase modulation for the tubes (3), insertion loss of tubes (4), phase modulation for different coil revolutions (5), attenuation of different horns (6), insertion loss for different coil revolutions (7), FDTD simulations (8), Electronic circuit (9), and dynamic refocusing (10). The clip “Dynamic Trap” shows a simulation of the trap being refocused dynamically using 2 channels. The clip “Step by Step assembly” shows how to assemble a portable tractor beam step by step. See supplementary material for information about the construction of the sonic devices, the driver board, as well as the simulated and experimental fields. It also contains experimental measures and model fits for the directivity of each delay line (supplementary Figures 1 and 2),for the tubes (3), insertion loss of tubes (4),for differentrevolutions (5), attenuation of different horns (6), insertion loss for differentrevolutions (7),simulations (8), Electronic circuit (9), and dynamic refocusing (10). The clip “Dynamic Trap” shows a simulation of the trap being refocused dynamically using 2 channels. The clip “Step by Step assembly” shows how to assemble a portable tractor beam step by step.