acoustic tweezers acoustic tweezers acoustic tweezers (or acoustical tweezers) are used to manipulate the position and movement of very small objects with sound waves. strictly speaking, only a single-beam based configuration can be called acoustical tweezers. however, the broad concept of acoustical tweezers involves two configurations of beams: single beam and standing waves. the technology works by controlling the position of acoustic pressure nodes that draw objects to specific locations of a standing acoustic field. the target object must be considerably smaller than the wavelength of sound used, and the technology is typically used to manipulate microscopic particles. acoustic waves have been proven safe for biological objects, making them ideal for biomedical applications. recently, applications for acoustic tweezers have been found in manipulating sub-millimetre objects, such as flow cytometry, cell separation, cell trapping, single-cell manipulation, and nanomaterial manipulation. the use of one-dimensional standing waves to manipulate small particles was first reported in the 1982 research article "ultrasonic inspection of fiber suspensions". method in a standing acoustic field, objects experience an acoustic-radiation force that moves them to specific regions of the field. depending on an object's properties, such as density and compressibility, it can be induced to move to either acoustic pressure nodes (minimum pressure regions) or pressure antinodes (maximum pressure regions). as a result, by controlling the position of these nodes, the precise movement of objects using sound waves is feasible. acoustic tweezers do not require expensive equipment or complex experimental setups. fundamental theory particles in an acoustic field can be moved by forces originating from the interaction among the acoustic waves, fluid, and particles. these forces (including acoustic radiation force, secondary field force between particles, and stokes drag force) create the phenomena of acoustophoresis, which is the foundation of the acoustic tweezers technology. acoustic radiation force when a particle is suspended in the field of a sound wave, an acoustic radiation force that has risen from the scattering of the acoustic waves is exerted on the particle. this was first modeled and analyzed for incompressible particles in an ideal fluid by louis king in 1934. yosioka and kawasima calculated the acoustic radiation force on compressible particles in a plane wave field in 1955. gorkov summarized the previous work and proposed equations to determine the average force acting on a particle in an arbitrary acoustical field when its size is much smaller than the wavelength of the sound. recently, bruus revisited the problem and gave detailed derivation for the acoustic radiation force. fig. 1: acoustic radiation force on a small particle as shown in figure 1, the acoustic radiation force on a small particle results from a non-uniform flux of momentum in the near-field region around the particle, f r a d = − ∇ u =-\nabla , which is caused by the incoming acoustic waves and the scattering on the surface of the particle when acoustic waves propagate through it. for a compressible spherical particle with a diameter much smaller than the wavelength of acoustic waves in an ideal fluid, the acoustic radiation force can be calculated by f r a d = − ∇ u =-\nabla , where u is a given quantity, also called acoustic potential energy. the acoustic potential energy is expressed as: u = v 0 ( p i n 2 ¯ 2 ρ f c f 2 f 1 − 3 ρ f v i n 2 ¯ 4 f 2 ) (^ \over c_^-^ \over 4})} where v 0 is the particle volume, p i n is the acoustic pressure, v i n is the velocity of acoustic particles, ρ f _ is the fluid mass density, c f is the speed of sound of the fluid, < ∗ > ¯ is the time-average term, 1 t ∫ 0 t ( ∗ ) d t \int \limits _^dt} the coefficients f 1 and f 2 can be calculated by f 1 = 1 − ρ f c f 2 ρ p c p 2 =1-c_^ \over c_^ and f 2 = 2 ( ρ p − ρ f ) 2 ρ p + ρ f =-)} \over + where ρ p _ is the mass density of the particle, c p _ is the speed of sound of the particle. acoustic radiation force in standing waves the standing waves can form a stable acoustic potential energy field, so they are able to create stable acoustic radiation force distribution, which is desirable for many acoustic tweezers applications. for one-dimension planar standing waves, the acoustic fields are given by: a i n ( x , t ) = − p 0 ρ f ω sin ⁡ ( k x ) sin ⁡ ( ω t ) (x,t)= \over \omega }}\sin(kx)\sin(\omega t)} , p i n ( x , t ) = p 0 cos ⁡ ( k x ) sin ⁡ ( ω t ) (x,t)=\cos(kx)\sin(\omega t)} , v in ( x , t ) = − p 0 ρ f c f sin ⁡ ( k x ) cos ⁡ ( ω t ) ) e x _(x,t)=c_\sin(kx)\cos(\omega t))_ , where a i n is the displacement of acoustic particle, p 0 is the acoustic pressure amplitude, ω is the angular velocity, k is the wave number. with these fields, the time-average terms can be obtained. these are: p i n 2 ¯ = 1 2 p 0 2 cos 2 ( k x ) ^=p_^(kx)} , v i n 2 ¯ = p 0 2 2 ρ f 2 c f 2 sin 2 ( k x ) ^=^ \over ^c_^(kx)} , thus, the acoustic potential energy is: u = v 0 p 0 2 4 ρ f c f 2 [ cos 2 ( k x ) f 1 − 3 2 sin 2 ( k x ) f 2 ] ^ \over c_^ , then, the acoustic radiation force is found by differentiation: f x r a d = − ∂ x u = v 0 k e a c sin ⁡ ( 2 k x ) φ ^=-u=k\sin(2kx)\phi } , e a c = p 0 2 4 ρ f c f 2 =^ \over c_^ , φ = f 1 + 3 2 f 2 = 5 ρ p − 2 ρ f 2 ρ p + ρ f − ρ f c f 2 ρ p c p 2 +=-2 \over +-c_^ \over c_^ , fig. 2: positions of pressure nodes (pn) and antinodes (an) along acoustic pressure waveform where e a c is the acoustic energy density, and φ is acoustophoretic contrast factor. the term sin ⁡ ( 2 k x ) )} shows that the radiation force period is one-half of the pressure period. also, the contrast factor can be positive or negative depending on the properties of particles and fluid. for positive value of φ , the radiation force points from the pressure antinodes to the pressure nodes, as shown in figure 2, and the particles will be pushed to the pressure nodes. secondary acoustic forces when multiple particles in a suspension are exposed to a standing wave field, they will not only experience acoustic radiation force, but also secondary acoustic forces caused by waves scattered by other particles. the inter-particle forces are sometimes called bjerknes forces. a simplified equation for the inter-particle forces of identical particles is: f b ( x ) = 4 π a 6 [ ( ρ p − ρ f ) 2 ( 3 cos 2 θ − 1 ) 6 ρ f d 4 v i n 2 ( x ) − ω 2 ρ f 9 d 2 ( 1 ρ p c p 2 − 1 ρ f c f 2 ) 2 p i n 2 ( x ) ] (x)=4\pi where a is the radius of the particle, d is the distance between the particles, θ is the angle between the central line of the particles and the direction of propagation of the incident acoustic wave. the sign of the force represents its direction: a negative sign for an attractive force, and a positive sign for a repulsive force. the left side of the equation depends on the acoustic particle velocity amplitude v i n ( x ) _()} and the right side depends on the acoustic pressure amplitude p i n ( x ) _()} . the velocity-dependent term is repulsive when particles are aligned with wave propagation (θ=0°), and negative when perpendicular to wave propagation (θ=90°). the pressure-dependent term is unaffected by the particle orientation and is always attractive. in the case of a positive contrast factor, the velocity-dependent term diminishes as particles are driven to the velocity node (pressure antinode), as in the case of air bubbles and lipid vesicles. in a similar way, the pressure-dependent term diminishes as particles are driven towards the pressure node (velocity antinode), as are most solid particles in aqueous solutions. in addition to the scattering-related secondary acoustic forces, the flow field resulting from the interactions of the various acoustic streaming fields, generated by the acoustic boundary layer of each particle (sometimes called microstreaming), can induce additional viscous shear forces on each of the particles' surfaces, which then results in an additional contribution to the secondary acoustic forces in its fully viscous formulations. the viscous effects on the secondary acoustic force can become significant when compared to the perfect fluid formulation exemplified above, and even dominant in certain limit cases, yielding both quantitatively and qualitatively different results than what is predicted by inviscid theory. the relevance of the viscous contributions varies greatly depending on the specific case being investigated, and thus important care needs to be taken in selecting an appropriate secondary acoustic force model for the given scenario. the influence of the secondary forces is usually very weak, and only has an effect when the distance between particles is very small. it becomes important in aggregation and sedimentation applications, where particles are initially gathered in nodes by the acoustic radiation force. as inter-particle distances become smaller, the secondary forces assist in further aggregation until the clusters become heavy enough for sedimentation to begin. acoustic streaming acoustic streaming is a steady flow generated by the nonlinear component of the oscillations in an acoustic field. depending on the mechanisms, the acoustic streaming can be categorized into two general types, eckart streaming and rayleigh streaming. eckart streaming is driven by a time-average momentum flux created when high-amplitude acoustic waves propagate and attenuate in a fluid. rayleigh streaming, also called "boundary driven streaming", is forced by reynolds stresses in the viscous boundary layer. both of the driven mechanisms come from a time-average nonlinear effect. a perturbation approach is used to analyze the phenomenon of nonlinear acoustic streaming. the governing equations for this problem are mass conservation and navier-stokes equations： ∂ t ρ = − ∇ ⋅ ( ρ v ) =-\nabla \cdot ()} , ρ [ ∂ t v + ( v ⋅ ∇ ) v ] = − ∇ p + μ ∇ 2 v + β μ ∇ ( ∇ ⋅ v ) ) =-\nabla +\nabla ^+\nabla (\nabla \cdot ))} where ρ is the density of fluid, v is the velocity of fluid particle, p is the pressure, μ is the dynamic viscosity of fluid, β is the viscosity ratio. the perturbation series can be written as p = p 0 + p 1 + p 2 ++ , v = 0 + v 1 + v 2 =+_+_ , ρ = ρ 0 + ρ 1 + ρ 2 =_+_+_ , which are diminishing series with the higher-order terms much smaller than the lower-order ones. the liquid is quiescent and homogeneous at its zero-order state. substituting the perturbation series into the mass conservation and navier-stokes equation and using the relation of p 1 = c f 2 ρ 1 =c_^ , the first-order equations can be obtained by collecting terms in first-order, ∂ t p 1 = − ρ 0 c f 2 ∇ ⋅ v 1 _=-__^\nabla \cdot _ , ρ 0 ∂ t v 1 = − ∇ p 1 + μ ∇ 2 v 1 + β μ ∇ ( ∇ ⋅ v 1 ) _\partial __=-\nabla _+\nabla ^_+\nabla (\nabla \cdot _)} . similarly, the second-order equations can be found as well,