Abstract
Sound can levitate objects of different sizes and materials through air, water and tissue. This allows us to manipulate cells, liquids, compounds or living things without touching or contaminating them. However, acoustic levitation has required the targets to be enclosed with acoustic elements or had limited manoeuvrability. Here we optimize the phases used to drive an ultrasonic phased array and show that acoustic levitation can be employed to translate, rotate and manipulate particles using even a singlesided emitter. Furthermore, we introduce the holographic acoustic elements framework that permits the rapid generation of traps and provides a bridge between optical and acoustical trapping. Acoustic structures shaped as tweezers, twisters or bottles emerge as the optimum mechanisms for tractor beams or containerless transportation. Singlebeam levitation could manipulate particles inside our body for applications in targeted drug delivery or acoustically controlled micromachines that do not interfere with magnetic resonance imaging.
Introduction
Acoustic waves can exert radiation forces^{1} and form acoustic traps at points where these forces converge permitting the levitation of particles of a wide range of materials and sizes^{2} through air^{3}, water^{4} or biological tissues^{5}. This is of paramount importance for crystallography^{6}, cell manipulation^{7}, labonachip scenarios^{8}, biomaterials^{9}, containerless transportation^{3,10} and even the levitation of living things^{11}.
With previous acoustic levitators, the trapped particles had to be enclosed by acoustic elements^{3,10,12,13,14}. Singlesided (or singlebeam) levitators only exerted lateral trapping forces^{15,16}, pulling forces^{17,18} or required the use of an acoustic lens^{19}. Furthermore, translation^{3,12,13,14} and rotation^{20} of the traps were limited.
Singleaxis levitators^{3,6,10,13} are a common arrangement for generating acoustic traps. They consist of an acoustic transducer and a reflector or another transducer above it. This generates a standing wave between the two elements and the nodes of the wave act as trap. By changing the phase difference between the transducers, the traps move in a single dimension without mechanical actuation. Various configurations for twodimensional manipulation have been explored, for example, a flat array of transducers and a parallel reflector provides movement within the plane of the array^{3,13}. Alternatively, an inwardfacing circular array of transducers can translate^{4,12} and rotate^{20} a particle within the circle. Threedimensional (3D) translation is possible with four arrays placed forming a square^{14} and recently with two opposed arrays^{21}.
Recent progress has seen custommade piezoelectric elements being used to create traps with a singlesided device (acoustic tweezers)^{15,16}. However, these traps only exert lateral forces and thus the particles have to rest on a surface. Pulling forces acting counter to the propagation direction (tractor beams) have been measured in water using triangularshaped particles^{18} and in air using acoustic bottle beams^{17}. Full 3D trapping with a singlesided device has been shown theoretically^{22,23} and a static underwater 3D trap has recently been reported^{19}. Nonetheless, a physical acoustic lens was required, introducing considerable energy loss^{16} and fixing the position of the trap to the focal point.
Controlled 3D trapping, translation and rotation with a singlesided array would enable acoustic tweezers to become the largerscale counterparts of optical tweezers^{24}, opening up applications in materials processing, microscale manufacturing and biomedicine.
Here we demonstrate simultaneous 3D acoustic trapping, translation and rotation of levitated particles using a singlesided array operating in air. This is achieved by optimally adjusting the phase delays used to drive an array of transducers; in this way unprecedented acoustic structures are generated without resorting to physical lenses, custom transducers or mechanical actuation. Our approach generates optimum traps at the target positions with any spatial arrangement of transducers and significantly enhances previous manipulators^{3,12,14}. We report three optimum acoustic traps: tweezerlike twin traps, a novel acoustic phenomenon with the ability to also rotate objects; twisterlike vortex traps, whose levitation capabilities were shown theoretically^{22,23} and recently observed experimentally using a fixed acoustic lens^{19}; and bottleshaped traps, never proven or suggested to levitate objects before^{17}. We also introduce the holographic acoustic element framework based on interpreting the phase delays as a holographic plate that combines the encoding of identifiable acoustic elements. The framework permits the analysis and efficient generation of acoustic traps as well as comparisons with optical traps. This work brings the advantages of optical tweezing (that is, singlebeam, rotation, holographic control and multiple particles)^{24} to the efficiency and versatility of acoustic levitation and could lead to the development of powerful tractor beams, 3D physical displays or acoustically controlled in vivo micromachines that do not interfere with magnetic resonance imaging.
Results
Universal optimizer
We characterize a 3D trap as a point towards which the forces converge from all directions. More explicitly, the Gor’kov potential^{1} defines a field, the gradient of which gives the forces exerted on small spheres; therefore, the Laplacian operator applied to the Gor’kov potential represents the trapping strength at a certain point. The Gor’kov Laplacian function at one position in space can be expressed as a nonlinear infinitely differentiable function with the phase delays (modulations) applied to the transducers as the only variables. With this function and the gradient of its variables, we employ a Broyden–Fletcher–Goldfarb–Shanno (BFGS) optimizer^{25} to obtain the phase modulations for the transducers so that when driven with a reference signal the generated acoustic field exerts maximum trapping forces on a particle situated at the target point. Our formulations of the Gor’kov Laplacian and its gradient enable realtime optimization.
Maximizing the Gor’kov Laplacian at a point sets the phase modulation of the transducers to generate a focal point at that position. In theory, a focal point can trap a particle exactly at its centre, where all the amplitude gradient forces cancel each other and the velocity gradient forces drag the particle in; amplitude gradients push the particles from highamplitude regions to lowamplitude ones, whereas velocity gradients displace particles towards regions with high complex gradients of the acoustic field (see Methods, equation (3)). However, a focal point is only a theoretical solution; experimentally, dense particles are repelled by the focal point^{22}and it is not possible to levitate particles around a focal point in a stable manner (Supplementary Fig. 1 and Supplementary Movie 1). Consequently, our optimizer uses an objective function that simultaneously maximizes the Gor’kov Laplacian and minimizes the pressure amplitude at the target point. These silent acoustic traps are the counterpart of dark optical traps^{26}. In addition, weights are applied to each component of the Gor’kov Laplacian to control the trapping strength in each dimension (see Methods, equation (9)).
This optimization method can be applied to scenarios with reflectors and any spatial arrangement of acoustic elements. Therefore, we can use it to control and enhance previously suggested manipulators. The improvements on the working volume for some arrangements from the literature^{3,12,14} are presented in Fig. 1 as a comparative qualitative representation (Supplementary Movie 2). This illustrates the benefit of using an optimization approach over the current positioning algorithms. More importantly, we show here that the optimization method can for the first time trap, translate and rotate particles using singlesided arrays (Fig. 2). Depending on the spatial arrangement of the array and the weights selected for each dimension, different acoustic traps are created. For a detailed description of the arrangements, see Supplementary Figs 2 and 3 and Supplementary Note 1.
Optimal singlebeam acoustic traps
The three optimal traps that emerge as optimum solutions for singlesided arrays are twin (Fig. 3), vortex (Fig. 4) and bottle (Fig. 5) traps. Experimental measures of these traps are presented in Supplementary Figs 4–8 and Supplementary Notes 2 and 3.
Acoustic traps can be analysed in terms of the origin of the exerted forces; namely, radiation forces are generated by amplitude gradients or velocity gradients^{1} (see Methods, equation (3)). In addition, phase singularities can be used to characterize the traps. Phase singularities are regions with zero amplitude and thus where the phase is not defined^{27}.
As a novel method to analyse traps, we introduce the concept of holographic acoustic elements. The phase modulation applied to the transducers is interpreted as a holographic plate that when driven with a reference signal renders an acoustic field. In our case, the traps are encoded as the combination of two holographic acoustic elements: a holographic acoustic lens that generates a focal point at the trap position and an extra element dependent on the type of trap (Fig. 6). The lens is obtained by making all the emitted waves coincide in phase at the focal point. By subtracting this lens from the optimized total plate, the holographic signature of the trap is obtained. The signature is an interesting feature for analysing the traps as to some extent it is invariant to the levitation position and can be compared with existing holographical optical traps^{24,26}.
Twin traps emerge when equal weights are specified in vshape arrangements or a large x axis weight is used for other arrangements. These traps have two fingerlike cylindrical regions of high amplitude, which tweeze the particle with amplitude gradients in the x direction. Velocity gradients constrain in the other two axes. A plane phase singularity (that is, twodimensional) occurs between the cylinders in the x plane. The holographic signature has a πphase difference between the two halves of the array. By rotation of the reference coordinate system or the holographic signature, the tweezer structure and the clamped particle can be rotated. Twin traps are shown in operation in Fig. 2 and have never been reported theoretically or experimentally in acoustics or optics.
Vortex traps emerge when equal weights are used in a hemispherical cap or a flat array. The xy section of the trap shows a highamplitude ring that generates lateral trapping forces with amplitude gradients. Along the z axis, the trapping force is due to velocity gradients and the phase is a 3D corkscrew spiralling around a line phase singularity (that is, onedimensional). The holographic signature is a helicoidal pattern. A particle trapped in this vortex trap spins around its own axis following the signature direction due to transfer of angular momentum^{28,29}. In our experiments, only small particles could be trapped (diameter <0.12λ=1 mm), see Supplementary Fig. 9, Supplementary Note 4 and Supplementary Movie 3 for further details. Vortex traps in acoustics have been shown theoretically^{22,23} and recently experimentally using a fixed acoustic lens^{19}. We note that the acoustic vortex trap that emerges from our optimizer is equivalent to an optical vortex^{26} (Supplementary Note 5 and Supplementary Fig. 10).
Bottle traps emerge in all the arrangements when large weights are applied to the direction of propagation (z axis). These traps create a highamplitude cage around the levitation point and all the forces result from amplitude gradients. A point phase singularity (that is, zerodimensional) is found at the trap centre. Here the holographic signature is a circular region of πphase difference. Bottle traps have been reported in acoustics^{17} but their ability to levitate particles was never suggested or proved; in optics, they have been generated by intersecting two laser beams with different modes^{30}.
Twin traps and vortex traps have a similar working volume for the same arrangements (Supplementary Fig. 11 and Supplementary Note 6) and this is comparable to the working volume of a standard singleaxis levitator. The z axis range with the tested singlesided arrays was up to 40 mm. This range is sufficient for many applications and could be increased by using more powerful transducers, a different host medium or if it not were necessary to defy gravity (for example, underwater applications). Bottle traps were limited in working volume since the lateral forces were weak, and the bottle shape was not maintained when the trap was generated offcentre.
Particles could be transported horizontally at up to 26 cm s^{−1} (Supplementary Tables 1 and 2 and Supplementary Note 7), this speed being limited by the update rate of our custom electronics (Supplementary Note 1). Vortex and twin traps achieved similar horizontal transport speeds and much higher vertical speeds than a traditional standing wave created with a twosided device. Twin traps are not symmetric around the z axis leading to faster transport speeds in the y direction than in the x direction. Bottle traps were limited in transport speed because of their relatively weak lateral forces. Other systems have reported maintained speeds of up to 4.9 (ref. 3), 3.2 (ref. 23) and 7 cm s^{−1} (ref. 20). The accuracy of particle repositioning (Supplementary Tables 3 and 4 and Supplementary Note 8) was at least 0.4 mm (λ/21) and up to 0.05 mm (λ/171) depending on the trap and axis, which is comparable to that achieved in previous devices^{12}.
Forces of the order of μN were generated on particles due to the trapping forces, this is comparable to previous levitation systems^{12,17}. The trapping forces are presented as spring constants in Supplementary Fig. 12, Supplementary Tables 5 and 6, and Supplementary Note 9. Twin and vortex traps had lateral forces comparable to a standing wave generated with a traditional singleaxis device but the zdirection forces were around 30 times weaker. Bottle traps had zdirection forces seven times stronger than twin and vortex traps, potentially leading to a greater z range, but their low lateral forces diminished their manoeuvrability. That is, it was possible to levitate particles at relatively large distances form the array, but not to reach those positions moving the particle from the central position. In general, pressure amplitude gradients were seen to generate much stronger forces than velocity gradients. The strength of the traps affected negatively the speed of transport but this was caused by the update rate of the electronics. Strong trapping forces require smaller step sizes for stable transport (Supplementary Tables 7 and 8 and Supplementary Note 7) and thus faster updates of the phases.
Holographic acoustic framework
The holographic acoustic element framework can also be used as a fast method to generate traps at different positions, rotate and spin particles as well as to create multiple levitation points (Fig. 7). Any trap can be generated at different locations by adding its signature to the phase delays that generate a focal point at the desired position, thus moving the trap is like refocusing the holographic acoustic lens. Rotating the holographic signature of a twin trap makes the trap structure and the trapped particles to rotate. Vortex traps transfer angular orbital momentum to the levitated particle with the same direction as the signature. And, when a holographic signature is added to a plate that generates multiple focal points, these points get transformed into traps that are of the same type as the added signature.
Discussion
Until now, only standing waves^{3,10,14,20,21} or Bessel beams^{4,12} were capable of translating levitated particles. On the other hand, singlesided arrays required an acoustic lens and generated static traps^{19}. Here we have presented an optimization method that creates optimal traps at the desired positions with different array geometries. It can directly control previous manipulators offering better results in terms of working volume. More importantly, the method can be applied to singlesided arrays and generates some unprecedented acoustic structures (that is, twin traps).
The introduction of three acoustic structures for the translation and rotation of levitated particles will find applications in tractor beams, containerless handling of matter and tangible displays. Our systems use inexpensive lowpower transducers but highpower versions could enable longer range 3D transportation, orientation and assembly of heavier objects. Singlesided devices potentially enable in vivo manipulation since the device could be applied directly onto the skin with the manipulation taking place inside the body; similar to an ultrasound scanner but for manipulating particles (that is, drug capsules, kidney stones or microsurgical instruments). This is a significant advantage over twosided opposed arrangements, which require the target area to be sandwiched by the arrays; also, singlebeam traps do not have repeated patterns that could accidentally trap other particles or generate undesired secondary maxima.
We also introduced the holographic acoustic framework that allows the traps to be generated without iterative methods. A direct link between optical and acoustic trapping has now been established and we expect this to yield further advances in both fields.
Methods
Characterizing a Levitation Point
The acoustic radiation force (F) exerted on a small spherical particle can be calculated from the gradient of the Gor’kov potential^{1} U:
We characterize a levitation point as a maximum of the Laplacian operator (convergence of the gradient) applied to the Gor’kov potential, that is, a point towards which all the forces converge.
where , and a=x, y, z are the Cartesian axes.
The Gor’kov potential, U, in terms of the complex acoustic pressure (p) and its spatial derivatives is given by:
where V is the volume of the spherical particle, ω is the frequency of the emitted waves, ρ is the density and c is the speed of sound (with the subscripts 0 and p referring to the host medium and the particle material, respectively). In equation (3), the first term relates to the amplitude gradient and establishes that particles are moved from regions with high amplitude towards regions with low amplitude; the second term relates to the velocity gradient and establishes that particles are dragged into regions with large modulus of the gradient of the complex field.
If an acoustic transducer emits with a constant frequency and amplitude, then the complex pressure that the jth transducer creates at a point can be expressed as:
where is the phase delay of the transducer and is a complex number that is constant for a given transducer and point in space. Owing to linearity, this also holds true for the spatial derivatives of the pressure, for instance .
To predict the constants and their spatial derivatives, several methods can be used; namely, the matrix method, finite differences or experimental measures. We employ a farfield model of a circular piston source:
where P_{0} is a constant defined by the transducers power, J_{0} is a zerothorder Bessel function of the first kind, k is the wavenumber , r is the radius of the piston, d_{j} is the distance between the transducer and the point, and θ_{j} is the angle between the transducer normal and the point. The piston model was adequate as the simulations matched the experiments for both the complex acoustic field (Supplementary Figs 4–8 and Supplementary Note 3) and the predicted levitation positions (Supplementary Tables 9 and 10 and Supplementary Note 10).
The total acoustic field (p) generated by N transducers is the addition of the individual fields, that is, . This also holds true for its spatial derivatives, for instance .
Therefore, the Laplacian of the Gor’kov potential can be expressed as a function of the phase delay of the transducers, .
Objective function
For a stable levitation trap, the Laplacian of the Gor’kov potential (which we termed the Gor’kov Laplacian) must be maximized and the modulus of the pressure (amplitude) must be minimized. Therefore, the objective function to minimize is:
which can be expanded with the addition of individual weights for the Cartesian axes:
where w_{x}, w_{y} and w_{z} are weights used to accentuate or damp the trapping forces in particular directions; large weights are proportions of 1,000 to 1. w_{p} is used to specify the balance between maximizing the Gor’kov Laplacian and minimizing the amplitude; here a value of 1 is used.
Efficient evaluation of the objective function and its gradient
Each term of the objective function can be explicitly expressed in terms of the pressure and its spatial derivatives:
where and the operator ‘·’ is defined as:
The objective function can be differentiated with respect to the phase of the jth transducer by applying the following formula:
where f, g∈{, x, y, z, xy, xz, yz, xx, yy, zz, xxy, xxz, xyy, xzz, yzz, xxx, yyy, zzz}. That is, p_{f} and p_{g} can be the complex pressure or any of its spatial derivatives.
Using equations (10)–(13), , , we can assemble the objective function (equation (9)) and its gradient for a given point in space (that is, the required levitation point). Our algorithm is highly efficient as evaluating the target function (equation (9)) or its gradient at one point has time complexity O(N) where N is the number of transducers. In addition, the employed optimizer (BFGS) presents superlinear convergence.
Once the complex constants for the desired levitation point have been calculated, the Gor’kov Laplacian and its gradient can be evaluated as follows: First, calculate the pressure and its spatial derivatives that each transducer creates at the target point given the transducer phase delay, . Second, calculate the total pressure and its spatial derivatives, . Third, evaluate the objective function (equation (9)) using equations (10)–(12), , and the previously calculated p_{f}. And finally, calculate the gradient of the objective function by taking the derivative of the objective function over each of its variables . This can be done by applying equation (13) to each term of equations (10) and (11) using the previously calculated p_{f} and per transducer .
BFGS optimization
The BFGS algorithm is an iterative method for unconstrained nonlinear optimization. Contrary to Newton’s method, it does not require the timeconsuming calculation of the inverse Hessian matrix. At every step, the optimizer needs to evaluate the function to minimize and its gradient.
For the linear search strategy, we employed the Armijo–Camino rule with α=1, β=0.5 and σ=0.0001. We also tried Basin Hop (temperature=0.1; step size=120) to ensure that the global minimum was found but it was not necessary. In our experiments with 400 transducers, ∼9,000 iterations were sufficient to converge on a solution starting from a random set of phase delays.
BFGS optimizers^{25,31} with Basin Hop^{32} have been successfully used before to solve the structure of condensed matter^{33}, proteins^{34,35} or atoms configurations^{32,36}. Now, this powerful and versatile approach also brings marked improvements to acoustic levitation.
Additional information
How to cite this article: Marzo, A. et al. Holographic acoustic elements for manipulation of levitated objects. Nat. Commun. 6:8661 doi: 10.1038/ncomms9661 (2015).
References
 1
Bruus, H. Acoustofluidics 7: the acoustic radiation force on small particles. Lab Chip 12, 1014–1021 (2012).
 2
Brandt, E. H. Acoustic physics: suspended by sound. Nature 413, 474–475 (2001).
 3
Foresti, D., Nabavi, M., Klingauf, M., Ferrari, A. & Poulikakos, D. Acoustophoretic contactless transport and handling of matter in air. Proc. Natl Acad. Sci. USA 110, 12549–12554 (2013).
 4
Courtney, C. R. et al. Dexterous manipulation of microparticles using Besselfunction acoustic pressure fields. Appl. Phys. Lett. 102, 123508 (2013).
 5
Kang, S. T. & Yeh, C. K. Potentialwell model in acoustic tweezers. IEEE Trans. Ultrason. Ferroelect. Freq. Control 57, 1451–1459 (2010).
 6
Hong, Z. Y., Xie, W. J. & Wei, B. Acoustic levitation with selfadaptive flexible reflectors. Rev. Sci. Instrum. 82, 074904 (2011).
 7
Laurell, T., Petersson, F. & Nilsson, A. Chip integrated strategies for acoustic separation and manipulation of cells and particles. Chem. Soc. Rev. 36, 492–506 (2007).
 8
Ding, X. et al. Onchip manipulation of single microparticles, cells, and organisms using surface acoustic waves. Proc. Natl Acad. Sci. USA 109, 11105–11109 (2012).
 9
Weber, R. J. et al. Acoustic levitation: recent developments and emerging opportunities in biomaterials research. Eur. Biophys. J. 41, 397–403 (2012).
 10
Whymark, R. R. Acoustic field positioning for containerless processing. Ultrasonics 13, 251–261 (1975).
 11
Xie, W. J., Cao, C. D., Lü, Y. J., Hong, Z. Y. & Wei, B. Acoustic method for levitation of small living animals. Appl. Phys. Lett. 89, 214102 (2006).
 12
Seah, S., Drinkwater, B. W., Carter, T., Malkin, R. & Subramanian, S. Dexterous ultrasonic levitation of millimetersized objects in air. IEEE Trans. Ultrason. Ferroelect. Freq. Control 61, 1233–1236 (2014).
 13
GlynneJones, P. et al. Arraycontrolled ultrasonic manipulation of particles in planar acoustic resonator. IEEE Trans. Ultrason. Ferroelect. Freq. Control 59, 1258–1266 (2012).
 14
Ochiai, Y., Hoshi, T. & Rekimoto, J. Pixie dust: graphics generated by levitated and animated objects in computational acousticpotential field. ACM Trans. Graph. 33, 85 (2014).
 15
Lee, J. et al. Single beam acoustic trapping. Appl. Phys. Lett. 95, 073701 (2009).
 16
Lam, K. H. et al. Ultrahigh frequency lensless ultrasonic transducers for acoustic tweezers application. Biotechnol. Bioeng. 110, 881–886 (2013).
 17
Zhang, P. et al. Generation of acoustic selfbending and bottle beams by phase engineering. Nat. Commun. 5, 4316 (2014).
 18
Démoré, C. E. et al. Acoustic tractor beam. Phys. Rev. Lett. 112, 174302 (2014).
 19
Baresch, D., Thomas, J. L. & Marchiano, R. Observation of a singlebeam gradient force acoustical trap for elastic particles: acoustical tweezers. Preprint at http://arxiv.org/abs/1411.1912 (2014).
 20
Foresti, D. & Poulikakos, D. Acoustophoretic contactless elevation, orbital transport and spinning of matter in air. Phys. Rev. Lett. 112, 024301 (2014).
 21
Omirou, T., Marzo, A., Seah, S. A. & Subramanian, S. in Proceedings of the 33rd Annual ACM Conference on Human Factors in Computing Systems, 309–312 (New York, NY, USA, 2015).
 22
Baresch, D., Thomas, J. L. & Marchiano, R. Spherical vortex beams of high radial degree for enhanced singlebeam tweezers. J. Appl. Phys. 113, 184901 (2013).
 23
Silva, G. T. & Baggio, A. L. Designing singlebeam multitrapping acoustical tweezers. Ultrasonics 56, 449–455 (2015).
 24
Grier, D. G. A revolution in optical manipulation. Nature 424, 810–816 (2013).
 25
Liu, D. C. & Nocedal, J. On the limited memory BFGS method for large scale optimization. Math. Program. 45, 503–528 (1989).
 26
Curtis, J. E. & Grier, D. G. Structure of optical vortices. Phys. Rev. Lett. 90, 133901 (2003).
 27
Nye, J. F. & Berry, M. V. Dislocations in wave trains. Proc. R. Soc. Lond. A 336, 165–190 (1974).
 28
Hefner, B. T. & Marston, P. L. An acoustical helicoidal wave transducer with applications for the alignment of ultrasonic and underwater systems. J. Acoust. Soc. Am. 106, 3313–3316 (1999).
 29
VolkeSepúlveda, K., Santillán, A. O. & Boullosa, R. R. Transfer of angular momentum to matter from acoustical vortices in free space. Phys. Rev. Lett. 100, 024302 (2008).
 30
Arlt, J. & Padgett, M. J. Generation of a beam with a dark focus surrounded by regions of higher intensity: the optical bottle beam. Opt. Lett. 25, 191–193 (2000).
 31
Nash, S. G. & Nocedal, J. A numerical study of the limited memory BFGS method and the truncatedNewton method for large scale optimization. SIAM J. Optim. 1, 358–372 (1991).
 32
Wales, D. J. & Doye, J. P. Global optimization by basinhopping and the lowest energy structures of LennardJones clusters containing up to 110 atoms. J. Phys. Chem. A 101, 5111–5116 (1997).
 33
Goedecker, S. Minima hopping: an efficient search method for the global minimum of the potential energy surface of complex molecular systems. J. Chem. Phys. 120, 9911–9917 (2004).
 34
Verma, A., Schug, A., Lee, K. H. & Wenzel, W. Basin hopping simulations for allatom protein folding. J. Chem. Phys. 124, 044515 (2006).
 35
Prentiss, M. C., Wales, D. J. & Wolynes, P. G. Protein structure prediction using basinhopping. J. Chem. Phys. 128, 225106 (2008).
 36
Leary, R. H. & Doye, J. P. Tetrahedral global minimum for the 98atom LennardJones cluster. Phys. Rev. E 60, 6320–6322 (1999).
Acknowledgements
This work has been supported by the EC within the seventh framework programme through both the European Research Council (Starting Grant Agreement 278576) and the FET Open Scheme (under grant agreement no. 309191). At the University of Bristol we thank Matt Sutton for his help in the figures and video; and Tom LlewellynJones for assisting with the 3D pressure measurements.
Author information
Affiliations
Contributions
A.M. and B.W.D. designed, developed and implemented the algorithms and simulations; A.M. and S.A.S. measured the acoustic slices; A.M and D.R.S. measured the spring constants; A.M. conducted the rest of the experiments and wrote the paper; all the authors contributed to the discussion and edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 112, Supplementary Tables 110 and Supplementary Notes 110. (PDF 1209 kb)
Supplementary Movie 1
2_LevitationFocal. A particle levitating unstably around the focal point until it drops. The particle levitates 19cm. above a flat array. ·3_SpinningOutOfVortex. A particle levitates in a Vortex trap, starts to orbit and then gets ejected. (MOV 3305 kb)
Supplementary Movie 2
1_MainVideo. Main video that summarizes the paper. From 0:09 to 0:36, singlesided levitation, rotation and spinning of particles are shown simultaneously in the simulation and experiments. From 0:39 to 0:47, full acoustic trapping with levitators placed at different angles. From 0:49 to 1:54, explanation of the optimization method. Amplitude isosurface, amplitude slices, phase slices, potential and holographic decomposition of Twin traps (2:01 to 2:21), Vortex traps (2:23 to 2:45) and Bottle traps (2:45 to 3:03). From 3:05 to 3:29, dynamic holographic decomposition of moving traps. From 3:30 to 3:49, holographic method for multiple traps. From 3:50 to 4:10, enhancing previous manipulators. (MOV 47867 kb)
Supplementary Movie 3
3_SpinningOutOfVortex. A particle levitates in a Vortex trap, starts to orbit and then gets ejected. (MOV 2557 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Marzo, A., Seah, S., Drinkwater, B. et al. Holographic acoustic elements for manipulation of levitated objects. Nat Commun 6, 8661 (2015). https://doi.org/10.1038/ncomms9661
Received:
Accepted:
Published:
Further reading

Three dimensional acoustic tweezers with vortex streaming
Communications Physics (2021)

Acoustic hologram optimisation using automatic differentiation
Scientific Reports (2021)

Emitting longdistance spiral airborne sound using lowprofile planar acoustic antenna
Nature Communications (2021)

Magnetic levitation using diamagnetism: Mechanism, applications and prospects
Science China Technological Sciences (2021)

Generation of Nonaliased Twodimensional Acoustic Vortex with Enclosed Metasurface
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.