Holographic acoustic elements for manipulation of levitated objects

1Acoustic waves can exert radiation forces and form acoustic traps at points where these forces converge permitting the levitation of particles of a wide range of materials and sizes through air, water or biological tissues. This is of paramount importance for crystallography, cell manipulation, lab-on-a-chip scenarios, biomaterials, containerless transportation and even the levitation of living things.

2With previous acoustic levitators, the trapped particles had to be enclosed by acoustic elements. Single-sided (or single-beam) levitators only exerted lateral trapping forces, pulling forces or required the use of an acoustic lens. Furthermore, translation and rotation of the traps were limited.

3Single-axis levitators 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 two-dimensional manipulation have been explored, for example, a flat array of transducers and a parallel reflector provides movement within the plane of the array. Alternatively, an inward-facing circular array of transducers can translate and rotate a particle within the circle. Three-dimensional (3D) translation is possible with four arrays placed forming a square and recently with two opposed arrays.

4Recent progress has seen custom-made piezoelectric elements being used to create traps with a single-sided device (acoustic tweezers). 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 triangular-shaped particles and in air using acoustic bottle beams. Full 3D trapping with a single-sided device has been shown theoretically and a static underwater 3D trap has recently been reported. Nonetheless, a physical acoustic lens was required, introducing considerable energy loss and fixing the position of the trap to the focal point.

5Controlled 3D trapping, translation and rotation with a single-sided array would enable acoustic tweezers to become the larger-scale counterparts of optical tweezers, opening up applications in materials processing, micro-scale manufacturing and biomedicine.

6Here we demonstrate simultaneous 3D acoustic trapping, translation and rotation of levitated particles using a single-sided 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. We report three optimum acoustic traps: tweezer-like twin traps, a novel acoustic phenomenon with the ability to also rotate objects; twister-like vortex traps, whose levitation capabilities were shown theoretically and recently observed experimentally using a fixed acoustic lens; and bottle-shaped traps, never proven or suggested to levitate objects before. 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, single-beam, rotation, holographic control and multiple particles) 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 micro-machines that do not interfere with magnetic resonance imaging.

7We characterize a 3D trap as a point towards which the forces converge from all directions. More explicitly, the Gor’kov potential 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 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 real-time optimization.

8Maximizing 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 high-amplitude regions to low-amplitude 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 pointand 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. 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)).

9This 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 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 single-sided 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.

10The three optimal traps that emerge as optimum solutions for single-sided 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.

11Acoustic traps can be analysed in terms of the origin of the exerted forces; namely, radiation forces are generated by amplitude gradients or velocity gradients (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.

12As 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.

13Twin traps emerge when equal weights are specified in v-shape arrangements or a large x axis weight is used for other arrangements. These traps have two finger-like 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, two-dimensional) 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 co-ordinate 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.

14Vortex traps emerge when equal weights are used in a hemispherical cap or a flat array. The xy section of the trap shows a high-amplitude 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, one-dimensional). 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. 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 and recently experimentally using a fixed acoustic lens. We note that the acoustic vortex trap that emerges from our optimizer is equivalent to an optical vortex (Supplementary Note 5 and Supplementary Fig. 10).

15Bottle traps emerge in all the arrangements when large weights are applied to the direction of propagation (z axis). These traps create a high-amplitude cage around the levitation point and all the forces result from amplitude gradients. A point phase singularity (that is, zero-dimensional) is found at the trap centre. Here the holographic signature is a circular region of π-phase difference. Bottle traps have been reported in acoustics 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.

16Twin 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 single-axis levitator. The z axis range with the tested single-sided 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 off-centre.

17Particles 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 two-sided 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.2 (ref. ) and 7 cm s−1 (ref. ). 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.

18Forces of the order of μN were generated on particles due to the trapping forces, this is comparable to previous levitation systems. 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 single-axis device but the z-direction forces were around 30 times weaker. Bottle traps had z-direction 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.

19The 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.

20Until now, only standing waves or Bessel beams were capable of translating levitated particles. On the other hand, single-sided arrays required an acoustic lens and generated static traps. 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 single-sided arrays and generates some unprecedented acoustic structures (that is, twin traps).

21The 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 low-power transducers but high-power versions could enable longer range 3D transportation, orientation and assembly of heavier objects. Single-sided 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 micro-surgical instruments). This is a significant advantage over two-sided opposed arrangements, which require the target area to be sandwiched by the arrays; also, single-beam traps do not have repeated patterns that could accidentally trap other particles or generate undesired secondary maxima.

22We 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.

23The acoustic radiation force (F) exerted on a small spherical particle can be calculated from the gradient of the Gor’kov potential U:

24We 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.

25where, and a=x, y, z are the Cartesian axes.

26The Gor’kov potential, U, in terms of the complex acoustic pressure (p) and its spatial derivatives is given by:

27where 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.

28If 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:

29where 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.

30To predict the constants and their spatial derivatives, several methods can be used; namely, the matrix method, finite differences or experimental measures. We employ a far-field model of a circular piston source:

31where P0 is a constant defined by the transducers power, J0 is a zeroth-order Bessel function of the first kind, k is the wavenumber, r is the radius of the piston, dj 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).

32The 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.

33Therefore, the Laplacian of the Gor’kov potential can be expressed as a function of the phase delay of the transducers,.

34For 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:

35which can be expanded with the addition of individual weights for the Cartesian axes:

36where wx, wy and wz are weights used to accentuate or damp the trapping forces in particular directions; large weights are proportions of 1,000 to 1. wp is used to specify the balance between maximizing the Gor’kov Laplacian and minimizing the amplitude; here a value of 1 is used.

37Each term of the objective function can be explicitly expressed in terms of the pressure and its spatial derivatives:

38where and the operator ‘·’ is defined as:

39The objective function can be differentiated with respect to the phase of the jth transducer by applying the following formula:

40where f, g∈{, x, y, z, xy, xz, yz, xx, yy, zz, xxy, xxz, xyy, xzz, yzz, xxx, yyy, zzz}. That is, pf and pg can be the complex pressure or any of its spatial derivatives.

41Using 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.

42Once 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 pf. 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 pf and per transducer.

43The BFGS algorithm is an iterative method for unconstrained nonlinear optimization. Contrary to Newton’s method, it does not require the time-consuming calculation of the inverse Hessian matrix. At every step, the optimizer needs to evaluate the function to minimize and its gradient.

44For 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.

45BFGS optimizers with Basin Hop have been successfully used before to solve the structure of condensed matter, proteins or atoms configurations. Now, this powerful and versatile approach also brings marked improvements to acoustic levitation.

46How 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).