Acoustic Macroscopic Rigid Body Levitation by Responsive Boundary Hologram

Propagated acoustic waves, which generate radiation pressure, exert a non-contact force on a remote object. By suitably designing the wave field, remote tweezers are produced that stably levitate particles in the air without any mechanical contact forces[23, 3, 5]. Recent works have revealed that holographic traps can levitate particles even with a single-sided wave source[8, 7, 24, 11]. However, the levitatable objects in the previous studies were limited to particles smaller than the wavelength, or flat parts placed near a rigid wall[1]. Here, we achieve a stable levitation of a macroscopic rigid body by a holographic design of acoustic field without any dynamic control. The levitator models the acoustic radiation force and torque applied to a rigid body by discretising the body’s surface, as well as the acoustic wave sources, and optimizes the acoustic field on the body surface to achieve the Lyapunov stability so that the field can properly respond to the fluctuation of the body position and rotation. In an experiment, a 40 kHz (8.5 mm wavelength) ultrasonic phased array levitated a polystyrene sphere and a regular octahedron with a size of ∼50 mm located 200 mm away from acoustic elements in the air. This method not only expands the variety of levitatable objects but also contributes to microscopic contexts, such as in-vivo micromachines, since shorter-wavelength ultrasound than the size of target objects can be used to achieve higher controllability and stability. ∗inoue@hapis.k.u-tokyo.ac.jp 1 ar X iv :1 70 8. 05 98 8v 1 [ ph ys ic s. ap pph ] 2 0 A ug 2 01 7 Acoustic radiation pressure is a static pressure applied to an object in a propagating or standing sound wave, and can apply a force to a remote object. Acoustic tweezers trap particles at points in space by utilizing acoustic radiation pressure. Compared with optical radiation force tweezers for nanoparticles or magnetic levitation targeting ferromagnetic materials, acoustic levitation by ultrasonic waves can levitate various materials, opening up the possibility of a wide range of applications, such as crystallography[15], chemistry[17], microbiology[4] and human–computer interfaces[14, 16]. However, previous midair acoustic levitation was limited to particles smaller than the wavelength of the sonic wave, where the wave field is designed under an assumption that sound field disturbance by the target object is negligible. An exceptional example of macroscopic object levitation is squeeze-film levitation. An ultrasonic vibrator can reduce the friction and lift an object larger than the wavelength using squeeze-film effect. However, the levitation distance is much smaller than the wavelength in typical designs[21]. Another approach is to create a standing wave between wave sources and the target object [25, 1]. Although it has been theoretically predicted that this approach can levitate an object at a distance of half the wavelength or integer multiples thereof, no experiments have been reported that demonstrate stable levitation in a three-dimensional space beyond distances equal to the wavelength. Basically, this approach cannot stabilize the transverse orientation of the object to the wave direction. A typical levitator for particles uses a standing waves generated by a pair of transducers or a set of a transducers and a reflector[23, 3, 5]. A recent research has expanded its degrees of freedom of manipulation into threedimensions[14]. In addition, single-sided traps formed by travelling waves have been realized recently, by extending the technique used in optical tweezers, first for lateral trapping[8] and also for full dimensional trapping[11]. The main approach used in acoustic tweezers is to apply an approximation model, such as the Gor’kov potential, that holds on a sphere significantly smaller than the wavelength, assuming the field is not affected by the scatterer[3, 7, 5, 11, 12, 10]. This approach is often highly effective for sparse particles because it becomes a free space problem, which is analytically solvable. Another approach is to employ a cloaking technique, where the wave avoids an object and is scattered less[24]. However, it yields less force, and it has not been reported that this technique gives a force sufficient to lift a solid in the air because the acoustic radiation pressure results from the scattered energy density[22]. Levitation by shorter wavelengths than the object size will not only expand the variety of levitatable targets but will also enable precise control of small particles. Levitation of a rigid body with a large clearance while stabilizing its orientation will be of benefit in a wide range of applications, including industrial assembly such as electronics, drug delivery, micro-machine control and multiple bodies control, and will make it possible to manipulate

which generate radiation pressure, exert a non-contact force on a remote object.By suitably designing the wave field, remote tweezers are produced that stably levitate particles in the air without any mechanical contact forces[23,3,5].Recent works have revealed that holographic traps can levitate particles even with a single-sided wave source[8,7,24,11].However, the levitatable objects in the previous studies were limited to particles smaller than the wavelength, or flat parts placed near a rigid wall[1].Here, we achieve a stable levitation of a macroscopic rigid body by a holographic design of acoustic field without any dynamic control.The levitator models the acoustic radiation force and torque applied to a rigid body by discretising the body's surface, as well as the acoustic wave sources, and optimizes the acoustic field on the body surface to achieve the Lyapunov stability so that the field can properly respond to the fluctuation of the body position and rotation.In an experiment, a 40 kHz (8.5 mm wavelength) ultrasonic phased array levitated a polystyrene sphere and a regular octahedron with a size of ∼50 mm located 200 mm away from acoustic elements in the air.This method not only expands the variety of levitatable objects but also contributes to microscopic contexts, such as in-vivo micromachines, since shorter-wavelength ultrasound than the size of target objects can be used to achieve higher controllability and stability.

Acoustic radiation pressure is a st tic pressure applied to an object in a propagating or standing sound wave, and can apply a force to a remote object.Acoustic tweezers trap particles at points in space by utilizing acoustic radiation pressure.Compared with optical radiation force tweezers for nanoparticles or magnetic levitation targeting ferromagnetic materials, acoustic levitation by ultrasonic waves can levitate various materials, opening up the possibility of a wide range of applications, such as crystallography [15], chemistry [17], microbiology [4] and human-computer interfaces [14,16].

However, previous midair acoustic lev tation was limited to particles smaller than the wavelength of the sonic wave, where the wave field is designed under an assumption that sound field disturbance by the target object is negligible.An exceptional example of macroscopic object levitation is squeeze-film levitation.An ultrasonic vibrator can reduce the friction and lift an object larger than the wavelength using squeeze-film effect.However, the levitation distance is much smaller than the wavelength in typical designs [21].Another approach is to create a standing wave between wave sources and the target object [25,1].Although it has been theoretically predicted that this approach can levitate an object at a distance of half the wavelength or integer multiples thereof, no experiments have been reported that demonstrate stable levitation in a three-dimensional space beyond distances equal to the wavelength.Basically, this approach cannot stabilize the transverse orientation of the object to the wave direction.

A typical levitator for particles uses a st nding waves generated by a pair of transducers or a set of a transducers and a reflector [23,3,5].A recent research has expanded its degrees of freedom of manipulation into threedimensions [14].In addition, single-sided traps formed by travelling waves have been realized recently, by extending the technique used in optical tweezers, first for lateral trapping [8] and also for full dimensional trapping [11].

The main approach used in acoustic tweezers i to apply an approximation model, such as the Gor'kov potential, that holds on a sphere significantly smaller than the wavelength, assuming the field is not affected by the scatterer [3,7,5,11,12,10].This approach is often highly effective for sparse particles because it becomes a free space problem, which is analytically solvable.Another approach is to employ a cloaking technique, where the wave avoids an object and is scattered less [24].However, it yields less force, and it has not been reported that this technique gives a force sufficient to lift a solid in the air because the acoustic radiation pressure results from the scattered energy density [22].

Levitation by shorter wavelengths than the objec size will not only expand the variety of levitatable targets but will also enable precise control of small particles.Levitation of a rigid body with a large clearance while stabilizing its orientation will be of benefit in a wide range of applications, including industrial assembly such as electronics, drug delivery, micro-machine control and multiple bodies control, and will make it possible to manipulate the object along any path, even paths that would be impossible to follow by mechanical arms due to interference.

In this paper we present a novel approach: respon ive boundary holography that designs wave field on impedance boundaries that can provide desirable restoring force and torque.The formulation accurately models the radiation pressure for a macro-scale rigid body by discretising its surface and optimizes the wavefield on the surface to lift and stabilize the rigid body.The condition for achieving stable rigid body levitation is that both the force and torque are in mechanical equilibrium and, in addition, adequate restoring force and torque are generated against the positional and rotational fluctuation of the rigid body.Our levitator seeks an optimum acoustic field that satisfies this condition for an arbitrarily given phased array and rigid body.As examples, we demonstrate that a sphere with a diameter of 30 mm and a regular octahedron with an edge length of 35.4 mm (a diagonal length of 50 mm) can be stably levitated by single-and double-sided phased arrays, respectively without any dynamic feedback controls.

In our approach, we derive the radiation pressure on scatterer by discretising not only the wave sources but also the surface of the rigid body.According to the Helmholtz-Kirchhoff equation, surface elements on the scatterer are considered as passive acoustic elements.This principle is historically employed in the Boundary Element Method (BEM) [18].As shown in Figure 1, we describe the system with active acoustic elements (transducers) and passive acoustic elements (surface of the rigid body).Since the transducers and the rigid body are far enough apart, re-reflection on the transducers is neglected.In this formulation, the radiation force, torque and their Jacobians at the levitation point are calculated as surface integrals of the Lagrange density and are represented by a quadratic matrix form of phased array gains, which can be efficiently handled on modern parallel processors.

The conditions for rigid body levitation are: 1) the exter al forces, which are typically gravity and acoustic radiation force, are balanced; 2) a restoring force is generated against positional fluctuation; 3) the net moment is zero; and 4) a restoring torque is generated against rotational fluctuation.In previous studies on particles, condition 1) is automatically satisfied by a position shift of the particle, and 3) and 4) are out of consideration.Marzo el al. solved condition 2) by maximizing the Gor'kov Laplacian to make it positive for particles [11].However, the Gor'kov potential becomes imprecise as the target becomes larger and aspherical, and the positive Laplacian is not a sufficient condition for stabilizing a rigid body even if it is small.We formulate an optimization problem to synthesize a phased array output as eigenvalue minimization with constraints.It is solvable by gradient methods such as the BFGS method.In other words, we optimize holographic surface along the boundary of the rigid body.When ignoring drag force, the rigid body levitates if and only if all eigenvalues of its Jacobian of generalized force in Method are non-positive real numbers in the neighbourhood where the linear approximation holds, which realizes the Lyapunov stability condition known in control theory (see Method, equation ( 15)).We obtain the optimum solution by minimizing the eigenvalues keeping the force-torque equilibrium condition.It is obvious the existence of the solution is not always guaranteed and it fails to find a solution satisfying the constraints when the target is too heavy or has a too complex shape for the given phased array.In the formulation, we assume that the body boundary is sound-hard since it is valid for most rigid materials in gaseous media.The viscosity of the air is neglected in this report, though it contributes to the further stability.

Figure 2 (a) and (d) show levitation of a sphere with a diameter of 30 mm (mass 0.6 g) and a regular octahedron with a diagonal length of 50 mm (mass 0.5 g) by 40 kHz (wavelength 8.5 mm) ultrasound in the air at a position 200 mm above the bottom phased array.In the figure, three green laser lines show an orthogonal basis crossing at the levitation point.Figure 2 (b) and (e) show FEM simulation of the absolute acoustic pressure on a plane perpendicular to the phased array.The left side of each shows the acoustic field in the absence of a target rigid body, and the right side shows the one in the presence of the target rigid body.The acoustic fields are significantly changed by the rig d bodies so as to support the mass.It will be difficult for readers to recognize the acoustic fields of the left-side (without body) and right-side (with body) are originating form an identical sound source, especially in case of the octahedron.Figure 2 (c) and (f) show absolute acoustic pressure on the rigid body surface calculated by our model.We can see that the pressure spreads over the surface.

The sphere was levitated by both a travelling wave type (only by the bottom phased array, which consists of 996 transducers) and a stationary wave type (by the bottom and top phased arrays, which consist of 1992 transducers), and the latter was more stable.In the case of the octahedron, the optimizer did not reach a stable solution for the single-sided phased array while it did reach a stable solution for the double-sided phased array.In actuality, the double-sided array levitated the octahedron but the singlesided array did not in our experiments.

Figure 3 shows the estimated and measured restoring force against the displacement around the equilibrium point, the origin, for the acoustic fields and rigid objects shown in Figure 2. Figure 3 (a) and (b) show the horizontal forces for horizontal displacement applied to the spher and the octahedron.The restoring force properly acted in a direction opposite to the displacement and converged at the origin.Figure 3 (c) shows the net resultant force including the restoring force and the gravity applied to the sphere in a vertical plane including the origin.The force converged at the origin in this case too, and the gradient of the force was observed along the vertical line.

The stable region of the rigid body is the neighbourhood of the levita-tion point where restoring force acts like linear spring since our optimizer constrains the gradient of the restoring force at the equilibrium point (see Figure 3 (a)).To see this, we tracked the position of a sphere with a diameter of 30 mm from various initial points.Figure 4 shows the results of ten successive trials of the levitation.The different initial positions of the sphere were set as shown in the figure at t = 0.The five spheres that dropped within four seconds were initially at positions more distant than 5 mm from the origin.The stable region of radius 5mm corresponded to the linear region in Figure 3 (a).The other five spheres were levitated for 50 seconds as the minimum and 14 seconds at the maximum.Unfortunately, perfect asymptotic stability has not been achieved yet.The reason of this instability is that a lossy factor is not included in the system.In the previous microscopic particle levitation, the lossy factor is provided by the fluid viscosity.But in a macroscopic scenario with a large Reynolds number, the lack of the loss makes the system unstable.

In summary, we achieved acoustic levitation for macroscopic rigid bodies.We employed a BEM formulation for acoustic radiation pressure applied to an arbitrary smooth surface and formulated an optimization problem that maximizes the stability while balancing all forces and torques, at a point distant enough from the acoustic elements.This model can be applied to arbitrary shaped rigid bodies and arbitrarily arranged phased arrays and even metamaterials [12,13,10].This levitation technique is applicable to three-dimensional manipulation of multiple objects to follow complex 3D paths in the air or liquid, which would be impossible to

chieve b
mechanical arms.Another aspect of the presented me hod is it enables use of higher frequency waves for micro particle manipulation, which improves the temporal response and strengthens the holding force by steepe lographic-optimization approach is also applicable to electromagnetic fields, which improves optical tweezers [6,19].


Methods

Acoustic radiation pressure on sound-hard boundary

The force F that acts on a rigid body surface RB due to acoustic radiation pressure is expressed by a Lagrangian density + ρ 0 RB dS • uu ,(1)
where ρ 0 is the density of the medium, u is the particle velocity, and • denotes the time domain average [22].By using the relative surface admittance β and assuming rmal n:
u = − 1 jc 0 kρ 0 ∇p, n • u = −ρ 0 β ∂p ∂t = −jkβp,(2)
where d of sound in the medium, and k is the wavenumber of the ultrasonic wave.We can rewr = RB 1 2ρ 0 c 2 0 p 2 − 1 2 ρ 0 u 2 dS + ρ 0 RB u(dS • u) (3) = 1 2ρ 0 c 2 0 RB p 2 n − 1 k

∇p 2 n + 2ρ 0 c 0 β p∇p dS.(4)
Especially on a sound-hard boundar
such that β = 0, this is simply:
F = 1 2ρ 0 c 2 0 RB p 2 − 1 k 2 ∇ p 2 dS,(5)
where ∇ p means the tangential derivative at is:
T = 1 2ρ 0 c 2 0 RB p 2 − 1 k 2 ∇ p 2 (r − c) × dS,(6)
where r is the position of surface elements and c is the centre of mass of the rigid body.


Scattered acoustic wave field exerted by phased array transducers

The pressure of an incident wave at a point r is described as an integral of wavefields generated where g is Green's Function for the Helmholtz Equation, q(s) is the normal particle velocity on the phased array surface, and PAT denotes the surface of the phased array transducers.The Kirchhoff-Helmholtz integral theory states that the scattered wave from a rigid body surface RB is written in the follo

) ∂n dS.(8)
Let the total pressure p
r) be the sum of the incident and scattered waves at point r on a sufficiently smooth surface of the rigid body.It is expressed by taking the one-sided limit of the sum onto the surface [2].On a sound-hard boundary, where the normal velocity is zero, they are combined to give:
1 2 p(r) = PAT q(s)g(r, s)dS − RB p(s) ∂g(r, s) ∂n dS.(9)

Discretization of acoustic elements

The equation ( 9) can be discretized by a known protocol used in the boundary element method (BEM) [18,2].We employ a tri d body and apply a continuous piecewise linear function to it, meaning that the dimension of the discretized pressure vector p becomes the number of points on the surface.Each of the phased array ele esented as a plane so that vector q gives the phased array gain.This discretization gives the following matrix equation in strong form:
1 2 p = Gq − G n p,(10)
where matrices G and G n represent acoustic c e torque T i along an axis i ∈ {x, y, z} can also be disc

tized from equations ( 5) and (6) using a n
rmal vector n and a discretized tangential differential operator D on the surface:
N i = diag [n 1,i , n 2,i , . . . ] M i = diag [((r 1 − c) × n 1 ) i , ((r 2 − c) × n 2 ) i , . . . ] F i = p * N i − 1 k 2 D * i N i D i p(12)T i = p * M i − 1 k 2 D * i M i D i p,(13)
where n l denotes nor al vector of l-th surface element.


Stable levitation condition for rigid body

Let a small perturbat et levitation point be ∆x = [∆x, ∆y, ∆z] T and that of rotation be ∆Ω = [∆Ω x , ∆Ω y , ∆Ω z ] T .We combine those two as a generalized position X = [∆x, ∆y, ∆z, ∆Ω x , ∆Ω y ] T .

For sufficiently small X , F is expressed in a linear approximation by using a Jacobian matrix ∇F as:
F (X ) F (0) + ∇F (0)X .(14)
When all the forces and torques acting upon the object balance each other at the levitation point and resistance from the medium is ignored, the equation of motion becomes

M d 2 X dt 2 = ∇F (0)X X (t) = X (0)
√ M −1 ∇F (0)t ,(15)
where M = diag(mI; I), using m and I as the mass and inertia matrix of the rigid body, respectively.Here, X (t) is bounded for a bounded initial position X (0) if and only if all eigenvalues of M −1 ∇F (0) are non-positive real numbers.This is called the Lyapunov stability condition in control theory [20].


Objective function and optimiza q that gives the exact force and torque to the body and satisfies the stability condition stated above.This problem is a so-called quadratically constrained quadratic program (QCQP), and many algorithms have been proposed to solve this.However, we empirically employ the following naïve unconstrained objective function:
min q 6 i=1 w i F i − Fi 2 + v i [λ i ],(16)
where w, v is a po itive weight parameter, F is the external force, such as gravity, and torque applied to the rigid body at the levitation point, and λ = eig(M −1 ∇F (0)) represents the eigenvalues.Although there is no explicit term to make the imaginary part of the eigenvalues zero in this formul

ion, minimization of the real part will spontan
ously place poles on the real axis.

We employ the L the transducer are assumed to be constant.The BFGS optimizer has been succes G * B −1 * N y − 1 k 2 D * y N y D y B −1 ∂G ∂x q(19)
∂T y ∂x = 2 q * G * B −1 * M y − 1 k 2 D * y M y D y B −1 ∂G ∂x q (20)
∂F y ∂Ω x = 2 q * G * B −1 * N y − 1 k 2 D * y N y D y B −1 K x q

q * G * B −1 * N z − 1 k 2 D * y N z D
y B −1 Gq (21)∂T y ∂Ω x = 2 q * G * B −1 * M y − 1 k 2 D * y M y D y B −1 K x q − q * G * B −1 * M z − 1 k 2 D * y M z D y B −1 q, (22)
where ∂G ∂x and K x are discretized matrices of kernels ∂g(r,s) ∂x and (r×∇g(r, s)) x , respectively.These terms are analytically differentiable with respect to the phase of each transducer.


Implementation and experimental setup

We implemented the phased array with commercially available transducers (NIPPON CERAMIC T4010).Each transducer w s driven by a 40 kHz rectangular wave of V pk−pk = 24 V.The rigid bodies were made of extruded polystyrene foam.Each rigid body had a through-hole with a radius of 1 mm, and a pin inserted in the hole was driven by a mechanical arm to move the body to the initial points.The pin was silently removed after the acoustic radiation pressure was applied to the body.

The optimization parameters were w = (1, 1, 1, 0, 0, 0), v = (1, 1, 1, 0, 0, 0) for the sphere with a diameter of 30 mm, and w = (1, 1, 10, 1, 1, 1), v = (1, 1, 1, 1, 1, 1) for the regular octahedron with a diagonal length of 50 mm.

Triangular meshes with a mean edge size of 1 mm were employed to ensure sufficient precision in the experim nt.A Gauss quadrature rule of order 4 was used for generating coefficients of connectivity ma