
==== Front
ACS Appl Mater Interfaces
ACS Appl Mater Interfaces
am
aamick
ACS Applied Materials & Interfaces
1944-8244
1944-8252
American Chemical Society

39143893
10.1021/acsami.4c06480
Research Article
Acoustofluidic Diversity Achieved by Multiple Modes of Acoustic Waves Generated on Piezoelectric-Film-Coated Aluminum Sheets
https://orcid.org/0000-0001-6363-3452
Wang Yong †‡§∇
Li Xianbin ∥∇
Meng Hui †
https://orcid.org/0000-0001-5461-5930
Tao Ran ⊥§
https://orcid.org/0000-0002-1292-4142
Qian Jingui ∥
Fu Chen ⊥
Luo Jingting ⊥
https://orcid.org/0000-0003-3942-3046
Xie Jin *‡
https://orcid.org/0000-0001-9797-4036
Fu Yongqing *§
† Department of Mechanical Engineering, Hangzhou City University, Hangzhou 310015, China
‡ The State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China
§ Faculty of Engineering and Environment, University of Northumbria, Newcastle upon Tyne NE1 8ST, United Kingdom
∥ Anhui Province Key Laboratory of Measuring Theory and Precision Instrument, School of Instrument Science and Optoelectronics Engineering, Hefei University of Technology, Hefei 230009, China
⊥ Shenzhen Key Laboratory of Advanced Thin Films and Applications, College of Physics and Optoelectronic Engineering, Shenzhen University Shenzhen 518060, China
* Email: xiejin@zju.edu.cn.
* Email: richard.fu@northumbria.ac.uk.
15 08 2024
28 08 2024
16 34 4511945130
20 04 2024
08 08 2024
05 08 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

Excitation of multiple acoustic wave modes on a single chip is beneficial to implement diversified acoustofluidic functions. Conventional acoustic wave devices made of bulk LiNbO3 substrates generally generate few acoustic wave modes once the crystal-cut and electrode pattern are defined, limiting the realization of acoustofluidic diversity. In this paper, we demonstrated diversity of acoustofluidic behaviors using multiple modes of acoustic waves generated on piezoelectric-thin-film-coated aluminum sheets. Multiple acoustic wave modes were excited by varying the ratios between IDT pitch/wavelength and substrate thickness. Through systematic investigation of fluidic actuation behaviors and performances using these acoustic wave modes, we demonstrated fluidic actuation diversities using various acoustic wave modes and showed that the Rayleigh mode, pseudo-Rayleigh mode, and A0 mode of Lamb wave generally have better fluidic actuation performance than those of Sezawa mode and higher-order modes of Lamb wave, providing guidance for high-performance acoustofluidic actuation platform design. Additionally, we demonstrated diversified particle patterning functions, either on two sides of acoustic wave device or on a glass sheet by coupling acoustic waves into the glass using the gel. The pattern formation mechanisms were investigated through finite element simulations of acoustic pressure fields under different experimental configurations.

Acoustofluidic diversity
Multiple modes
Thin film acoustic waves
Fluidic actuation
Particle patterning
National Natural Science Foundation of China 10.13039/501100001809 52305619 Fundamental Research Funds for the Central Universities 10.13039/501100012226 JZ2024HGTB0243 Natural Science Foundation of Anhui Province 10.13039/501100003995 2308085QF196 National Natural Science Foundation of China 10.13039/501100001809 62301209 document-id-old-9am4c06480
document-id-new-14am4c06480
ccc-price
==== Body
pmc1 Introduction

Surface acoustic wave (SAW) technology has been demonstrated as an efficient and reliable method for manipulation of fluids and particles/cells, in either a digital format (e.g., using sessile droplets) or a continuous flow format (e.g., inside a polymer- or glass-based microchannel/chamber) with its attributes of biocompatibility, contactless and precise operation, and noninvasive features.1−7 Based on acoustic streaming and acoustic radiation effects induced by acoustic waves in the liquid,8,9 various microfluidic functions have been realized, including internal flowing/streaming,10,11 pumping,12−14 splitting/separation,15,16 jetting,17−19 nebulization,20,21 and particles/cells trapping,22−24 sorting,25−27 patterning,28−30 and alignment.31−33 Piezoelectric materials, either in their bulk forms (e.g., LiNbO3 or LiTaO3)34−36 or in thin film forms (e.g., ZnO or AlN)37−39 have been used for fabricating SAW-based acoustofluidic devices. Bulk ones, especially those based on LiNbO3 substrates, have been widely used due to their large piezoelectric constant, high electromechanical coupling coefficient (5%–11%) and high energy transduction efficiency.13,37 However, they are difficult to integrate with microelectronics for embedded control and signal processing, and also brittle and prone to fracture under a high radio frequency (RF) power and large thermoelectric shock.38,39 Moreover, it is difficult to excite multiple acoustic wave modes on LiNbO3 substrates once the crystal-cut and electrode’s pattern direction are defined. For example, a Rayleigh mode is generally generated on the 128° Y-cut LiNbO3 substrate, whereas a shear-horizontal SAW (SH-SAW) mode is usually excited on the 64° Y-cut LiNbO3 substrate.7 Multiple acoustic wave mode excitation on a single substrate or chip is beneficial for realizing diversified acoustofluidic functions.40,41 For instance, the Rayleigh wave is dominant for fluidic actuation due to its effective dissipation of acoustic wave energy into the liquid and the generation of an acoustic driving force, whereas the SH-SAW is suitable for biosensing in liquid environment on account of its slight signal attenuation in liquid. Additionally, excitation of multiple acoustic wave modes expands the range of available frequencies for various acoustofluidic manipulation, promotes the efficiency for microfluidic mixing, speeds up heat and mass transfer, and creates diversified flow fields by using multiple acoustic frequencies and vibration modes.42,43

Compared to their bulk counterparts, piezoelectric thin film acoustic wave technology has distinct advantages, in terms of device design’s flexibility (if considering in-plane isotropy of piezoelectric property), high breakdown voltage, low cost and easiness of fabrication and integration with other microelectronics.44 Besides, piezoelectric thin films, such as ZnO or AlN, can be easily deposited onto various substrates, including silicon,45,46 glass,47,48 polymers,49,50 and metals,51−53 thereby achieving various acoustic speeds, acoustic modes and applications such as the integration of fluidic actuation and biosensing functions on a single chip.41,46 In the past few years, various microfluidic actuation and particle/cell manipulation functions have been achieved using thin film acoustic wave devices fabricated on silicon, glass, and aluminum substrates with the Rayleigh mode. In our previous work, we have demonstrated microfluidic actuation/pumping functions using ZnO thin film acoustic waves on 50-μm-thick Al foil substrates with the zero-order mode (i.e., A0 and S0) of Lamb wave.12,54 Apart from the Rayleigh mode and zero-order mode of Lamb wave, there are also hybrid wave modes or higher-order modes of Lamb wave generation when the ratio between interdigital transducer (IDT) pitch/wavelength and substrate thickness is within an appropriate range.42,55 In addition, the substrate thickness plays an important role in fluidic actuation performance, as a much flexible device tends to generate a macroscopical substrate deformation during the actuation, resulting in poor fluid actuation performance. Until now, acoustofluidic behaviors (including fluidic pumping and particle patterning) under multiple acoustic wave modes, IDT pitches, and substrate thicknesses have not been systemically investigated, but these are important for realizing a diversity of acoustofluidic functions. For example, the selectively excited acoustic wave modes gain a good fluid actuation performance, the balance between the device flexibility and actuation performance maintains a good fluidic actuation performance for flexible devices, and the expansion of acoustic frequencies and resonances generate various flowing behaviors and particle patterns for multiple biological applications.

In this paper, we fabricated ZnO thin film acoustic wave devices with wavelengths varied from 100 μm to 1100 μm on Al sheet substrates with two thicknesses (i.e., 200 μm and 600 μm). We have achieved various acoustic wave mode excitations on ZnO-thin-film-coated Al sheets, including the Rayleigh, Sezawa, hybrid, and Lamb waves by varying the ratios between IDT pitch and substrate thickness. We systematically investigated fluidic actuation behaviors and performance using the generated acoustic wave modes and analyzed the actuation mechanism of the microfluidic behaviors. Moreover, we demonstrated diversified particle patterning abilities either on two sides (i.e., top and backside) of the Lamb wave device or on the glass sheet (by coupling of acoustic wave into the glass using the ultrasonic gel) with different orders of Lamb waves. We analyzed the formation mechanisms of particle patterns through finite-element analysis (FEA). Our work shows great potential for establishing diversified acoustofluidic platforms using multiple modes of acoustic waves and different experimental configurations.

2 Materials and Methods

ZnO thin films with a thickness of ∼5 μm were deposited onto commercially available Al sheet substrates (Metal Sheets Limited, U.K.) with thickness of 200 and 600 μm using a direct current (DC) magnetron sputtering system (Model NS3750, Nordiko). A zinc target with a purity of 99.99% was used for film deposition. The deposition was performed with a DC power of 400 W, a chamber pressure of ∼3 mTorr, and an Ar/O2 gas flow ratio of 3/10 sccm. The distance between the zinc target and the sample holder was 20 mm, and the sample holder was rotated to obtain uniform ZnO thin films. Crystal orientation of the deposited ZnO thin film on the Al sheet substrate was analyzed by using X-ray diffraction (XRD) (Model D5000, Siemens) with Cu Kα radiation (λ = 1.5406 Å). Cross-sectional morphology of the ZnO thin film was observed using a scanning electron microscopy (SEM) system (Model S-4100, Hitachi). The characterization results of ZnO thin film on Al sheet substrate are shown in Figure S1 in the Supporting Information. IDT electrodes of acoustic wave devices were patterned by evaporating a layer of 150-nm-thick Al on a ZnO/Al sheet substrate. There were 60 pairs of finger electrodes on each IDT with an acoustic aperture of 5 mm and wavelengths that varied from 100 to 400 μm, as shown in Figures S2(a) and S2(b). For the acoustic wave device with the wavelength of 1100 μm, the number of finger pairs was 20 and the acoustic aperture was 15 mm, as shown in Figure S2(c). The reflection spectra (S11) of the acoustic wave devices were measured using an RF network analyzer (Agilent, Model E5061B).

To understand acoustic wave vibration modes and acoustic pressure fields under different experimental configurations, FEA simulations were performed using COMSOL Multiphysics (6.1) software with solid mechanics, electrostatics, and pressure acoustic modules. For simulations of acoustic wave modes, a simplified two-dimensional (2D) model with one pair of IDT electrodes and infinite boundary conditions was used. The schematic of 2D model and meshing of 2D model are shown in Figures S3(a) and S3(b), respectively. Here, the thickness of ZnO thin film was 5 μm, and the electrode thickness was 150 nm. The IDT pitch was varied from 100 to 1100 μm, and the Al sheet substrate thickness was set as 200 and 600 μm, respectively. A polarization voltage of 1 V was assigned to one of the IDT electrodes, while another was assigned to be ground. To really reflect the wave vibration mode, the bottom boundary was set to be free, whereas the boundaries of lateral walls were set as periodic, as shown in Figure S3(b). Triangular meshes was used to mesh the 2D model with the maximum element units ranging from 10 to 30 μm and the minimum element units varied from 0.1 to 0.3 μm, depending on the IDT pitch and substrate thickness.

For simulations of acoustic pressure fields, a 3D model of acoustic wave device composed of the Al electrode, Al substrate, and ZnO thin film layer was built, as shown in Figure S3(c). Here, the electrode thickness was 150 nm, the thickness of the Al substrate was 200 μm, the thickness of ZnO thin film was 5 μm, and the IDT pitch/wavelength was 1100 μm with an acoustic aperture of 8 mm. The diameter of droplet model was 3.2 mm with a height of 1.2 mm. A PDMS model with dimensions of 3.5 mm (W) × 1.5 mm (H) × 4 mm (L) and wall thickness of 0.4 mm was used. The thick of glass model was 200 μm. Free tetrahedral meshes were used for meshing the 3D model with the maximum and minimum mesh size limited to 300 and 162 μm, respectively, as illustrated in Figure S3(d). Mesh size is a compromise between computational accuracy and computational cost. Therefore, the maximum element growth rate was set as 1.6, thereby avoiding too small mesh and inverted surface mesh. For simulations of acoustic pressure fields on the device surface, the sessile droplet or PDMS chamber containing the liquid was directly placed on the device surface, then the acoustic wave was coupled into the liquid. Whereas, for the acoustic pressure field simulation on the glass sheet, the acoustic wave was first coupled into the glass, and then the propagating acoustic wave on the glass was further coupled into the liquid. For each pair of IDT electrodes, a polarization voltage of 10 V was assigned to one of the electrodes, while the other one was assigned to be ground. Finally, acoustic pressure fields under different experimental configurations were simulated in the frequency domain analysis using the corresponding frequency values obtained from the eigenfrequency analysis. The detailed material parameters were obtained from the literature.41,55

For performing acoustofluidic tests, the acoustic wave device was put on an aluminum alloy test holder to minimize potential acoustic heating effect.38 An RF signal was generated using a signal source (RIGOL, DSG 815) and amplified using a power amplifier (Aigtek, Model ATA-1222A) before being input into the IDT electrodes (Figure S4). The input RF power was measured using a feedthrough-type power meter (Model SHX-200W-3601-HUV, SHX). The used RF power range for particle patterning was varied from 0.2 to 1 W, and the RF power range for fluidic actuation was varied from 1 to 30 W. For fluidic actuation tests, a layer of 200-nm-thick CYTOP (Asahi Glass Co., Tokyo, Japan) was coated on the surface of the acoustic wave device to make the surface hydrophobic. After the hydrophobic treatment, the droplet’s contact angle on the ZnO thin film surface is ∼109°. To visualize particle patterns induced by the acoustic waves, both silica (with a diameter of 5 μm) and polystyrene (with a diameter of 10 μm) were used, and they were added into the deionized (DI) water to prepare different particle solutions, respectively. Then, the particle solution was deposited dropwise on the device’s top surface or injected into the rectangular polydimethylsiloxane (PDMS) chamber bonded onto the top or backside of the device to observe the formation of particle patterns. The PDMS chamber with dimensions of 3.5 mm × 1.5 mm × 4 mm (W × H × L) was fabricated by using the standard soft photolithography process. A piece of glass with a thickness of 50 μm was put on top of the PDMS chamber to avoid the evaporation and also induced the wave reflection. The particle solution was also dropped on a glass sheet (200 μm thick) or injected into the gap between two pieces of glass to observe the formed particle patterns under acoustic wave agitations. Here, an ultrasonic gel (AQUASONIC100) was used to effectively transfer the acoustic wave energy into the glass from the acoustic wave device. The acoustic wave mode, which was linked with a specific wavelength, could be selectively excited by applying the corresponding frequency at the resonant peak. A CCD camera (NPX-GS130UM, 100 frames/s) was used to capture the motions of the fluid and particle trajectories.

3 Results and Discussion

Previous studies have shown that the ratio r (r = λ/h) between the IDT pitch λ and substrate thickness h played an important role in exciting various acoustic wave modes.42,55 When the ratio r was ≪1, the Rayleigh mode was dominant, whereas the Lamb wave mode was dominant when the ratio r was ≫1. Both the Rayleigh mode and the Lamb wave mode would be hybridized together when the ratio r was ∼1. Figure S5 shows FEA simulation results of acoustic wave modes and experimentally measured reflection spectra for ZnO thin film acoustic wave devices fabricated on a 600-μm-thick Al sheet with wavelengths varied from 100 to 400 μm. Here, the acoustic wave mode at the resonant peak of reflection spectrum was determined by comparing the resonant peak frequency with the simulated frequencies at different acoustic wave modes. We have found that when the device wavelengths were 100 and 200 μm, in which the ratio r was much smaller than one, both the Rayleigh mode and Sezawa mode were excited. As the device wavelengths were increased to 300 and 400 μm, hybrid pseudo-Rayleigh mode and pseudo-A0 mode of Lamb wave were obtained based on the results from both FEA simulation and experiment measurement. Whereas from the experimental results, we have not found distinct resonant peak near the frequency of simulated pseudo-S0 mode. Here, the pseudo-Rayleigh mode, pseudo-A0 mode, and pseudo-S0 mode are defined as hybrid acoustic modes, which have shown the similar vibration patterns to those of the fundamental Rayleigh mode, A0 mode and S0 mode, respectively.

To further expand the generated acoustic wave modes, the acoustic wave devices were fabricated on a thinner substrate (200 μm). Figure S6 shows FEA simulation and experimental measurement of wave vibration modes for ZnO thin film acoustic wave devices fabricated on a 200-μm-thick Al sheet with wavelengths varied from 100 to 1100 μm. When the device wavelength was 200 μm and the ratio r was 1, a hybridization between the Rayleigh mode and A0 mode of the Lamb wave was observed. Additionally, both the pseudo-S0 mode and the Sezawa mode were observed from FEA simulation and experiment measurement. As the ratio was further increased to 1.5 and 2 (corresponding to device wavelengths of 300 and 400 μm), the first-order mode (i.e., pseudo-A1 and pseudo-S1) of Lamb wave was observed, apart from the generated zero-order mode (i.e., A0 and S0). Moreover, there was a hybrid S1 mode and Sezawa mode for the acoustic wave device with a wavelength of 300 μm. When the ratio r was increased to 5.5 with a wavelength of 1100 μm, the acoustic wave mode was a pure Lamb wave, i.e., zero-order and higher-harmonic modes. Therefore, as the ratio r was increased from 0.17 to 5.5, the acoustic wave mode was changed from the Rayleigh mode (Sezawa mode, also called the two-order Raleigh wave) to hybrid mode and then to Lamb wave mode.

3.1 Fluid Actuation Using Multiple Modes of Acoustic Waves

The microfluidic behaviors induced by the acoustic waves are determined by the input RF power. When the input RF power was relatively low (e.g., tens of milliwatts), acoustic streaming effects within the droplet were observed.37 As the input RF power was increased to the level of a few watts, the droplet vibration, pumping and even jetting could be achieved.56 Here, we focus on droplet pumping or transportation behaviors by using different acoustic wave modes. Figure 1a shows droplet pumping phenomena on a 600-μm-thick Al sheet substrate using various acoustic wave modes under different wavelengths or IDT pitches. For a single-pitched IDT, the specific acoustic wave mode was selectively excited by using the corresponding frequency. For the Rayleigh mode and pseudozero-order mode (i.e., pseudo-A0 and pseudo-S0) of Lamb wave, the droplet motion on the Al sheet is a combination of rolling and sliding (Movie S1). This is mainly due to a large diffraction angle (∼32°) between the acoustic wave and the liquid, which is often called Rayleigh angle (θR = sin–1(CF/CS), where CF is the sound speed in a liquid, CS is the acoustic wave propagation speed in the substrate) for the Rayleigh wave. Whereas, for the Sezawa wave mode, the droplet movement tends to be a combination of jumping and sliding (Movie S2), due to the high acoustic speed of the Sezawa wave, thus quickly dissipating acoustic energy. Also, the smaller value of diffraction angle (∼16°) could easily lift up the droplet.

Figure 1 (a) Droplet pumping images using ZnO thin film acoustic wave devices (with wavelengths varied from 100 to 400 μm) fabricated on a 600-μm-thick Al sheet substrate with different acoustic wave modes. All scale bars are 1 mm. (b) Droplet average pumping velocities for ZnO thin film acoustic wave devices fabricated on 600-μm-thick Al sheet substrates with different acoustic wave modes under different input powers. (c) Droplet average pumping velocities for ZnO thin film acoustic wave devices fabricated on 600-μm-thick Al sheet substrates using the Sezawa mode under different input powers, the inset shows threshold pumping power (the minimum power to initiate the droplet) for different wavelength acoustic wave devices using different acoustic wave modes.

The droplet’s average pumping velocities under different acoustic wave modes were measured, with results shown in Figures 1b and 1c. As the input RF power was increased, the droplet pumping velocity was increased. For acoustic wave devices with wavelengths of 100–300 μm, the droplet pumping velocity using the Rayleigh mode or pseudo-A0 mode of Lamb wave was increased with the device wavelength, because a large wavelength can result in a long acoustic wave penetration length (i.e., the interaction distance between the acoustic wave and the liquid),57 thereby improving fluidic actuation efficiency for this given droplet size. For the acoustic wave device with a wavelength of 400 μm, the droplet pumping velocity was slightly smaller than that of the 300 μm one using either pseudo-A0 mode or pseudo-S0 mode. This is mainly because when the ratio r is closer to 1, the acoustic wave mode tends to change from the Rayleigh mode to the Lamb wave mode. The Rayleigh wave are propagating more near the device’s surface, thereby having a high energy conversion efficiency. In addition, fluidic pumping performance using the Rayleigh mode or pseudo-A0 mode is distinctly better than that of the pseudo-S0 mode and Sezawa mode. For the Sezawa mode, the droplet pumping performance is improved with an increase of device wavelength, as shown in Figure 1c. This is because as the ratio r is increased, the Sezawa mode tends to convert into the A1 mode of the Lamb wave, and the pumping performance using the A1 mode is much better than that of the Sezawa mode, as demonstrated in Figure 2c.

Figure 2 (a) Droplet pumping images using ZnO thin film acoustic wave devices (with wavelengths of 100, 200, and 400 μm) fabricated on a 200-μm-thick Al sheet substrate with different acoustic wave modes. All scale bars are 1 mm. (b) Droplet average pumping velocities for ZnO thin film acoustic wave devices fabricated on a 200-μm-thick Al sheet substrate with different acoustic wave modes under different input powers. (c) Threshold pumping powers for different wavelength acoustic wave devices using different acoustic wave modes.

To further investigate effects of acoustic wave modes/orders on fluidic actuation behaviors and performance, we selected the acoustic wave devices with an Al substrate thickness of 200 μm. Figure 2a shows droplet pumping behaviors using ZnO thin film acoustic waves on a 200-μm-thick Al sheet substrate with different acoustic wave modes. For different cases of Rayleigh mode, hybrid mode (e.g., pseudo-A0 and pseudo-S0), and zero-order mode (i.e., A0 and S0) of Lamb wave, the droplet movement on the Al sheet substrate was found to be a combination of rolling and sliding (Movie S3). The droplet pumping behaviors shown in Movies S1 and S3 are similar, indicating that fluidic actuation behaviors are mainly determined by the acoustic wave modes rather than the substrate thickness. Because as long as the used acoustic wave mode is the same, the droplet motion on the Al sheet looks similar, no matter how thick the Al sheet is. Whereas for first-order mode (i.e., A1 and S1) of Lamb wave case, the droplet movement was dominated by jumping and sliding (Movie S4). The droplet pumping behaviors using the first-order mode of Lamb wave were similar to those of the Sezawa mode, which is attributed to a high acoustic speed of the first-order mode of Lamb wave and a small diffraction angle (∼16°), which tends to lift up the droplet. Another possible reason is that the surface vibration mode of S1 mode is partly similar to that of the Sezawa mode (Figure S6).

Figure 2b shows droplet average pumping velocities for different wavelength acoustic wave devices with different modes. The droplet pumping velocities using the Rayleigh mode, A0 mode, or A1 mode were much higher than those of the S0 mode and S1 mode at the same input RF power. In addition, the threshold pumping powers using the Rayleigh mode, A0 mode, or A1 mode were much smaller than those of Sezawa mode, S0 mode, and S1 mode, as shown in Figure 2c, indicating that the Rayleigh mode and antisymmetric mode (e.g., A0 and A1) of Lamb wave present better pumping performance. Moreover, the pumping performance using the zero-order mode (i.e., A0 and S0) of the Lamb wave was generally better than that of the first-order mode (i.e., A1 and S1), revealing that the lower-order mode acoustic waves generally show a better fluid actuation performance.

When the ratio r is increased to 5.5, i.e., with a wavelength of 1100 μm, a zero-order and higher-order mode (e.g., A1 and S1, A2 and S2) of Lamb waves would be excited. Here, we further investigated fluidic pumping behaviors and performance using different orders and modes of the Lamb wave. Figure 3a shows droplet pumping behaviors using a zero-order mode (i.e., A0 and S0) of the Lamb wave. For the A0 mode, there was a threshold value for the droplet volume or size (∼5 μL) to drive the droplet, above which, the droplet could be actuated forward (Movie S5). Otherwise, the droplet could only be horizontally stretched (side view; Movie S6). This is because when the device wavelength is larger than or near the droplet size, part of acoustic wave will penetrate the droplet and generate a strong wave reflection at the back of the droplet,33 thereby preventing the droplet to move forward. For the S0 mode, there was no obvious size limitation to actuate the droplet, and the droplet tended to jump or jet under the actuation of an acoustic wave. We also studied the droplet pumping using higher-order modes (e.g., A2 and S2), with the results shown in Figure 3b. The obtained threshold pumping powers using the higher-order mode were generally much higher, as shown in Figure 3c. For example, for the S2 mode of the Lamb wave, even though the input RF power was increased to tens of watts, the droplet was still difficult to be transported forward. In brief, the fluidic pumping performance using the antisymmetric modes (e.g., A1 and A2) of Lamb wave is generally better than that of symmetric modes (e.g., S1 and S2).

Figure 3 (a) Droplet pumping images with different volumes (2 to 20 μL) using ZnO thin film acoustic waves on a 200-μm-thick Al sheet substrate with the A0 mode and S0 mode. The wavelength of the acoustic wave device is 1100 μm. (b) Droplet pumping images using ZnO thin film acoustic waves on a 200-μm-thick Al sheet substrate with the A1 mode, S1 mode, A2 mode, and S2 mode. (c) Threshold pumping powers for ZnO thin film acoustic waves on a 200-μm-thick Al sheet substrate using different acoustic wave modes. All scale bars are 2 mm.

Although different acoustic wave modes may bring different fluid actuation behaviors, for practical application of acoustofluidic actuation devices, we should focus more on the fluidic actuation performance. Therefore, effects of IDT pitch/wavelength, substrate thickness, acoustic wave mode, and droplet size on fluidic actuation performance need to be carefully considered. For example, the increase of the device’s wavelength or IDT pitch can improve fluidic actuation performance to some extent, because the increased interaction distance between the acoustic wave and the liquid enhances the energy conversion efficiency. However, when the device wavelength is increased to nearly or slightly above the substrate thickness (i.e., the ratio r is near or larger than 1), there is an acoustic wave mode transition between the Rayleigh wave and Lamb wave, thereby resulting in different fluidic actuation performances. Additionally, the substrate thickness also has an important influence on fluidic actuation performance. A more flexible/bendable acoustic wave device may produce an obvious substrate deformation or significant damping effect during the actuation,54 thereby causing acoustic energy dissipating into the substrate and significant acoustic heating effect. Furthermore, for a given size droplet, when the device wavelength was increased to nearly or above the droplet size, the droplet pumping performance using the low-frequency zero-order mode of the Lamb wave may not be good due to significant wave reflection at the droplet boundary.

3.2 Particles Patterning Using Multiple Modes of Acoustic Waves

Apart from fluidic actuation functions, acoustic waves can also be used for manipulating particles within the fluid. The gravity effect on microparticles in liquid is commonly ignored due to their small sizes. Therefore, the particles suspended in the fluid mainly experience two acoustic forces: the acoustic radiation force and the drag force resulting from acoustic streaming.43,58,59 The acoustic radiation force applied to the particle can be expressed as451

2

3

where Frad corresponds to the acoustic radiation force, p0 is the acoustic pressure, V is the particle volume, βm is the compressibility of medium, βp is the compressibility of particles, λ is the wavelength, φ(β,ρ) is the acoustic contrast factor, k is the wavenumber, x is the distance of the particle from the pressure node, ρp is the particle density, ρm is density of the medium, f is the driving frequency, cL is the longitudinal wave velocity, and θ is the diffraction angle.

The acoustic streaming induced drag force acting on a particle can be calculated as follows:604

where μ is the fluid dynamic viscosity, Rp the particle radius, and v the relative velocity between the fluid and particles.

In general, the acoustic radiation force is the main driving force for particle manipulation, whereas the acoustic streaming induced drag force is mostly regarded as the disturbance. Previous studies showed that there was a threshold particle size to balance these two forces, above which, the acoustic radiation force became dominant.61 At a driving frequency of 2 MHz, the threshold particle size is ∼2.6 μm, and this size is proportional to f–1/2. For a particle patterning test, here, a ZnO-thin-film-based Lamb wave device fabricated on a 200-μm-thick Al sheet with wavelength of 1100 μm was used. The frequency range of different orders of the Lamb wave is from 1.39 to 11. Two MHz, corresponding to a threshold particle size ranging from 3.1 to 1.1 μm. Therefore, when the particle size is larger than 3.1 μm, the particle motion within the fluid will be mainly determined by the acoustic radiation force. Here, the used particle size is 5 and 10 μm, respectively. Thus, when the applied RF power is relatively low (neglecting acoustic streaming effects), under the actuation of acoustic radiation force, the particles will be accumulated at the pressure nodes corresponding to acoustic pressure fields,62 thereby generating various particle patterns.

We first investigated particle patterning within a droplet using the Lamb traveling wave. Previous study indicated that when the Lamb traveling wave was propagated into the droplet with a relatively large size (e.g., larger than 5 μL), a standing wave field was generated within the droplet due to the reflection at the droplet boundary.63 The standing wave field then produced a periodic acoustic pressure field within the droplet, thereby generating various particle patterns. Figures 4a–d show the simulated acoustic pressure field within the droplet positioned on the device surface using different orders of the Lamb traveling wave. Clearly, a checkerboard acoustic pressure field was formed under the actuation of different orders of Lamb wave. When the particle solution was placed on the device surface, under the actuation of an acoustic radiation force, the particles in the solution would accumulate at these pressure node areas, thereby forming the checkerboard patterns of particles (Movie S7), as shown in Figures 4e–h. As the acoustic wave mode was changed to a higher order or the resonant frequency of acoustic wave mode was increased, the size of the formed checkerboard patterns became much smaller, and the checkerboard patterns became irregular due to a weak acoustic field generated at higher-order modes. The results are consistent with the previous study.63Figures 4I and 4j show the simulated acoustic pressure fields inside the rectangular PDMS chamber placed on the device surface (in front of the IDT). When the Lamb traveling waves dissipated their energy into the PDMS chamber, the interference between the radiated traveling wave and the reflected acoustic wave would generate a standing wave field within the PDMS chamber,24,33 thereby forming linear pressure node patterns parallel to the PDMS channel walls as well as the IDT. When the particle solution was injected into the PDMS chamber, the particles would be accumulated on these pressure node lines and patterned into lines parallel to the IDTs (Movie S8), as shown in Figures 4k and 4l. As the frequency of the acoustic wave mode was increased, the distance between the adjacent lines was decreased. The simulation results show good agreement with the experiment ones.

Figure 4 Simulation of acoustic pressure fields within the droplet positioned on the device surface using different modes/orders of Lamb wave: (a) A0 mode (1.39 MHz), (b) S0 mode (4.81 MHz), (c) A1 mode (4.01 MHz), and (d) A2 mode (6.83 MHz). Polystyrene particle (diameter of 10 μm, red color) patterning within the droplet (50 μL) positioned on the device surface using different modes/orders of Lamb wave, (e) A0 mode (1.38 MHz, 0.2 W), (f) S0 mode (4.86 MHz, 0.3 W), (g) A1 mode (4.01 MHz, 0.3 W), and (h) A2 mode (6.81 MHz, 0.5 W). Simulation of acoustic pressure fields inside the PDMS chamber positioned on device surface (in front of IDT) using (i) the A0 mode (1.39 MHz) and (j) the S0 mode (4.81 MHz). Polystyrene particles (diameter of 10 μm, red color) patterning inside the PDMS chamber positioned on the device surface using (k) A0 mode (1.38 MHz, 0.3 W) and (l) S0 mode (4.86 MHz, 0.3 W). The scale bars in panels (e)–(h) are 500 μm, and the scale bars in panels (i)–(l) are 300 μm. The PDMS chamber is placed parallel to the IDT. The signs of minus and plus in the legend of the simulation represent the acoustic pressure direction, the termd “max” and “min” represent the maximum value of the acoustic pressure at the node and antinode.

The Lamb wave is a plate wave, propagating through the whole substrate; thus, the backside of the Lamb wave device could also generate particle patterns. Here, we further investigated particle patterning behaviors within a PDMS chamber positioned at the backside of a Lamb wave device. Figures 5a–d show the simulated acoustic pressure fields inside the PDMS chamber at the device’s backside, which are similar to those observed on the top surface. Linear patterns of pressure nodes parallel to the PDMS channel walls and the IDT were clearly formed inside the PDMS chamber placed at the device backside, because the interference between the incident Lamb traveling wave and the reflected wave generated the standing wave. When the particle solutions were injected into the PDMS chamber, linear particle patterns parallel to the IDT were formed (Movie S9) at the pressure nodes, as illustrated in Figures 5e–h. The formed particle patterns at the backside of Lamb wave device are similar to those formed on the device’s top surface, whereas the particle patterning at the device backside needs a higher RF power due to a weaker acoustic field. However, we should address that acoustofluidic manipulation at the device’s backside reduces the potential contamination to the IDT electrodes and is convenient for cleaning. The distance between the adjacent lines is decreased with the increase of frequency of acoustic wave mode, as summarized in Figure S7, which is consistent with the previous study.45

Figure 5 Simulation of acoustic pressure fields inside the PDMS chamber positioned at backside of Lamb wave device using different modes/orders of Lamb wave: (a) A0 mode (1.39 MHz), (b) S0 mode (4.81 MHz), (c) A2 mode (6.83 MHz), and (d) S1 mode (8.81 MHz). Silica particles (diameter of 5 μm, white color) patterning inside the PDMS chamber placed at the device backside using different modes/orders of Lamb traveling wave, (e) A0 mode (1.38 MHz, 0.3 W), (f) S0 mode (4.86 MHz, 0.5 W), (g) A2 mode (6.81 MHz, 0.5 W), and (h) S1 mode (8.84 MHz, 0.6 W). All scale bars are 200 μm. The PDMS chamber is positioned parallel to the IDT.

We further studied particle patterning formed by putting a particle solution on the top of the glass or injecting into the gap between two pieces of glass, as illustrated in Figures 6a–c. Here, the acoustic wave was coupled into the glass using the ultrasonic gel. Figure S8 shows that when the acoustic wave is propagated into the glass (positioned vertical to the acoustic wave device), the generated acoustic wave in the glass is also the Lamb wave. However, the acoustic wave vibration in the glass is influenced by the loaded droplet; for example, a local circular vibration region on the glass is observed at the placed position of the droplet. Figures 6d and 6e show simulated acoustic pressure fields within the droplet positioned on the glass sheet under the actuation of the A0 mode and S0 mode of the Lamb wave, respectively. Here, the glass sheet was positioned vertically to the acoustic wave device, as illustrated in Figure 6a. It has been shown that when the Lamb wave propagating on the glass was coupled into the droplet, a standing wave field was formed within the droplet, thereby generating annular pressure node patterns. Under the actuation of acoustic radiation force, the particles in solution were accumulated on these pressure node rings, thereby forming the annular patterns of particles (Movie S10), as shown in Figures 6f and 6g. Changing the acoustic wave modes or frequencies caused changes in sizes or intervals between adjacent annular patterns.

Figure 6 Schematic of experimental setup for generating particle patterns: (a) the droplet on top of the glass sheet (the glass sheet is positioned vertical to the acoustic wave device), (b) the liquid in the gap between two pieces of glass (two glass sheets are positioned vertical to the acoustic wave device), and (c) two droplets on the glass sheet (the glass sheet is placed parallel to the device surface). Simulation of acoustic pressure fields within the droplet positioned on top of the glass sheet was performed using (d) the A0 mode (1.39 MHz) and (e) the S0 mode (4.81 MHz). Polystyrene particle (diameter of 10 μm, red color) patterning within the droplet (50 μL) placed on top of the glass sheet was driven by (f) A0 mode (1.38 MHz, 1 W) and (g) S0 mode (4.86 MHz, 1 W). Simulation of acoustic pressure fields inside the liquid gap between two pieces of glass using the (h) A0 mode (1.39 MHz) and (i) S0 mode (4.81 MHz). Silica particle (diameter of 5 μm, white color) patterning inside the gap between two pieces of glass using (j) A0 mode (1.38 MHz, 1 W), (k) S0 mode (4.86 MHz, 1 W). Simulation of acoustic pressure fields within the droplet placed on the glass sheet (the glass sheet is placed parallel to the device surface) using A0 mode for (l) the sessile droplet near the acoustic wave device and (m) the sessile droplet far away from the acoustic wave device. Silica particle (diameter of 5 μm, white color) patterning within the sessile droplet placed on the glass sheet driven by the A0 mode (1.38 MHz, 1 W) for (n) the sessile droplet near the acoustic wave device and (o) the sessile droplet far away from the acoustic wave device. All scale bars are 500 μm.

When the particle solution was injected into the gap between two pieces of glass, as shown in Figure 6b, the faveolate particle patterns were formed (Movie S11), as shown in Figures 6j and 6k. Here, the white regions represent the aggregated particles, while the dark regions represent the surface of glass sheet. These have also been verified through FEA simulation of acoustic pressure fields, as shown in Figures 6h and 6i. The faveolate particle patterns could also be formed by putting a drop of particle solution on the glass sheet that was placed parallel to the device surface, as shown in Figures 6l–o. Therefore, for particle patterning functions induced by multimode acoustic waves, we can adjust the particle pattern sizes or shapes by either varying the frequency of acoustic wave mode or changing the experimental configurations. In addition, for particle patterning on the glass, this superstrate type of design avoids the contamination of the acoustic wave device by liquid samples, thereby enabling the reuse of acoustic wave devices. All of these results show that multimode acoustofluidic technologies have great potential in biomedical applications such as diversified particles/cell patterning, focusing, trapping, and solidified particle–hydrogel patterns.

4 Conclusions

In summary, we have achieved various acoustic wave mode excitations on piezoelectric ZnO-thin-film-coated aluminum sheets by varying the ratios between IDT pitch/wavelength and substrate thickness. Based on these acoustic wave modes, we systematically investigated fluidic actuation and particle patterning behaviors, demonstrating acoustofluidic diversity using multiple modes of acoustic waves. We provided a guidance for high-performance acoustofluidic actuation platform design by demonstrating that the Rayleigh mode, hybrid mode (e.g., pseudo-Rayleigh mode and A0 mode) generally have better fluidic actuation performance than those of the Sezawa mode, symmetric mode (e.g., S0 and S1) and higher-order modes of the Lamb wave. Also, we demonstrated diversified particle patterning generation either on both sides of the Lamb wave device or on the glass sheet (by coupling of acoustic wave into the glass using the gel) by varying the frequency of the acoustic wave mode and experimental configurations. The formation mechanisms of particle patterning were analyzed through the FEA simulation of acoustic pressure fields.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsami.4c06480.Droplet (1 μL) pumping using ZnO thin film acoustic wave device fabricated on the 600-μm-thick Al sheet (wavelength of 200 μm) substrate with the Rayleigh wave mode and an input RF power of 1.5 W (Movie S1) (AVI)

Droplet (1 μL) pumping using ZnO thin film acoustic wave device fabricated on a 600-μm-thick Al sheet substrate (wavelength of 200 μm) with the Sezawa mode and an input RF power of 21 W (Movie S2) (AVI)

Droplet (1 μL) pumping using ZnO thin film acoustic wave device fabricated on a 200-μm-thick Al sheet substrate (wavelength of 400 μm) with the A0 mode and an input RF power of 0.8 W (Movie S3) (AVI)

Droplet (1 μL) pumping using ZnO thin film acoustic wave device fabricated on a 200-μm-thick Al sheet substrate (wavelength of 400 μm) with the S1 mode and an input RF power of 3 W (Movie S4) (AVI)

Droplet (8 μL) pumping using ZnO thin film acoustic wave device fabricated on a 200-μm-thick Al sheet substrate (wavelength of 1100 μm) with the A0 mode and an input RF power of 10 W (Movie S5) (AVI)

Droplet (5 μL) pumping using ZnO thin film acoustic wave device fabricated on a 200-μm-thick Al sheet substrate (wavelength of 1100 μm) with the A0 mode and an input RF power of 10 W (Movie S6) (AVI)

Polystyrene particle (diameter of 10 μm) patterning within the droplet (50 μL) positioned on the device surface using the A0 mode (wavelength of 1100 μm) and an input RF power of 0.3 W (Movie S7) (AVI)

Polystyrene particle (diameter of 10 μm) patterning within the PDMS chamber positioned on the device surface using the S0 mode (wavelength of 1100 μm) and an input RF power of 0.3 W (Movie S8) (AVI)

Silica particle (diameter of 5 μm) patterning within the PDMS chamber positioned on the device backside using the A0 mode (wavelength of 1100 μm) and an input RF power of 0.5 W (Movie S9) (AVI)

Polystyrene particle (diameter of 10 μm) patterning within the droplet (50 μL) positioned on the top of glass sheet using the A0 mode (wavelength of 1100 μm) and an input RF power of 1 W (Movie S10) (AVI)

Silica particle (diameter of 5 μm) patterning inside the gap between two pieces of glass using the A0 mode (wavelength of 1100 μm) and an input RF power of 1 W (Movie S11) (AVI)

XRD pattern and cross-sectional SEM image of ZnO thin film on a 200-μm-thick Al sheet; optical images of the fabricated acoustic wave devices; schematic of the model and meshing of the model for simulation of acoustic wave mode and acoustic pressure field; schematic of experimental setup for fluid actuation; FEA simulation of acoustic wave mode and experimental measured reflection spectra for ZnO thin film acoustic wave device with wavelengths varied from 100 μm to 400 μm fabricated on 600-μm-thick Al sheets; FEA simulation of acoustic wave modes and experimental measured reflection spectra for ZnO thin film acoustic wave device with wavelengths varied from 100 μm to 1100 μm fabricated on 200-μm-thick Al sheets; distance (μm) between the adjacent lines as a function of acoustic wave mode frequency; FEA simulation of acoustic vibration on the acoustic wave device and on the glass sheet (PDF)

Supplementary Material

am4c06480_si_001.avi

am4c06480_si_002.avi

am4c06480_si_003.avi

am4c06480_si_004.avi

am4c06480_si_005.avi

am4c06480_si_006.avi

am4c06480_si_007.avi

am4c06480_si_008.avi

am4c06480_si_009.avi

am4c06480_si_010.avi

am4c06480_si_011.avi

am4c06480_si_012.pdf

Author Contributions

∇ Authors Y. Wang and X. Li contributed equally to this paper.

The authors declare no competing financial interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant Nos. 52305619 and 62301209), the Natural Science Foundation of Anhui Province (Grant No. 2308085QF196) and the Fundamental Research Funds for the Central Universities (Grant No. JZ2024HGTB0243).
==== Refs
References

Fan Y. ; Wang X. ; Ren J. ; Lin F. ; Wu J. Recent advances in acoustofluidic separation technology in biology. Microsyst. Nanoeng. 2022, 8 , 94 10.1038/s41378-022-00435-6.36060525
Yang S. ; Rufo J. ; Zhong R. ; Rich J. ; Wang Z. ; Lee L. P. ; Huang T. J. Acoustic tweezers for high-throughput single-cell analysis. Nat. Protoc. 2023, 18 , 2441–2458. 10.1038/s41596-023-00844-5.37468650
Ozcelik A. ; Rufo J. ; Guo F. ; Gu Y. ; Li P. ; Lata J. ; Huang T. J. Acoustic tweezers for the life sciences. Nat. Methods 2018, 15 , 1021–1028. 10.1038/s41592-018-0222-9.30478321
Stringer M. ; Zeng Z. ; Zhang X. ; Chai Y. ; Li W. ; Zhang J. ; Ong H. ; Liang D. ; Dong J. ; Li Y. ; Fu Y. ; Yang X. Methodologies, technologies, and strategies for acoustic streaming-based acoustofluidics. Appl. Phys. Rev. 2023, 10 , 011315 10.1063/5.0134646.
Wu M. ; Ozcelik A. ; Rufo J. ; Wang Z. ; Fang R. ; Huang T. J. Acoustofluidic separation of cells and particles. Microsyst. Nanoeng. 2019, 5 , 32 10.1038/s41378-019-0064-3.31231539
Rufo J. ; Cai F. ; Friend J. ; Wiklund M. ; Huang T. J. Acoustofluidics for biomedical applications. Nat. Rev. Methods Primers 2022, 2 , 30 10.1038/s43586-022-00109-7.
Fu Y. Q. ; Luo J. K. ; Nguyen N. T. ; Walton A. J. ; Flewitt A. J. ; Zu X. T. ; Li Y. ; McHale G. ; Matthews A. ; Iborra E. ; Du H. ; Milne W. I. Advances in Piezoelectric Thin Films for Acoustic Biosensors, Acoustofluidics and Lab-on-Chip Applications. Prog. Mater. Sci. 2017, 89 , 31–91. 10.1016/j.pmatsci.2017.04.006.
Friend J. ; Yeo L. Y. Microscale acoustofluidics: Microfluidics driven via acoustics and ultrasonics. Rev. Mod. Phys. 2011, 83 , 647 10.1103/RevModPhys.83.647.
Shi J. ; Ahmed D. ; Mao X. ; Lin S. C. S. ; Lawit A. ; Huang T. J. Acoustic tweezers: patterning cells and microparticles using standing surface acoustic waves (SSAW). Lab Chip 2009, 9 , 2890–2895. 10.1039/b910595f.19789740
Farooq U. ; Liu X. ; Zhou W. ; Hassan M. ; Niu L. ; Meng L. Cell lysis induced by nanowire collision based on acoustic streaming using surface acoustic waves. Sens. Actuators, B 2021, 345 , 130335 10.1016/j.snb.2021.130335.
Chen C. ; Zhang S. P. ; Mao Z. ; Nama N. ; Gu Y. ; Huang P. H. ; Jing Y. ; Guo X. ; Costanzo F. ; Huang T. J. Three-dimensional numerical simulation and experimental investigation of boundary-driven streaming in surface acoustic wave microfluidics. Lab Chip 2018, 18 , 3645–3654. 10.1039/C8LC00589C.30361727
Wang Y. ; Zhang Q. ; Tao R. ; Xie J. ; Canyelles-Pericas P. ; Torun H. ; Reboud J. ; McHale G. ; Dodd L. E. ; Yang X. ; Luo J. ; Wu Q. ; Fu Y. Flexible/bendable acoustofluidics based on thin-film surface acoustic waves on thin aluminum sheets. ACS Appl. Mater. Interfaces 2021, 13 , 16978–16986. 10.1021/acsami.0c22576.33813830
Wang Y. ; Tao X. ; Tao R. ; Zhou J. ; Zhang Q. ; Chen D. ; Jin H. ; Dong S. ; Xie J. ; Fu Y. Q. Acoustofluidics along inclined surfaces based on AlN/Si Rayleigh surface acoustic waves. Sens. Actuators, A 2020, 306 , 111967 10.1016/j.sna.2020.111967.
Yeo L. Y. ; Friend J. R. Ultrafast microfluidics using surface acoustic waves. Biomicrofluidics 2009, 3 , 012002 10.1063/1.3056040.19693383
Collignon S. ; Friend J. ; Yeo L. Planar microfluidic drop splitting and merging. Lab Chip 2015, 15 , 1942–1951. 10.1039/C4LC01453G.25738425
Jung J. H. ; Destgeer G. ; Ha B. ; Park J. ; Sung H. J. On-demand droplet splitting using surface acoustic waves. Lab Chip 2016, 16 , 3235–3243. 10.1039/C6LC00648E.27435869
Castro J. O. ; Ramesan S. ; Rezk A. R. ; Yeo L. Y. Continuous tuneable droplet ejection via pulsed surface acoustic wave jetting. Soft Matter 2018, 14 , 5721–5727. 10.1039/C7SM02534C.29845144
Jangi M. J. ; Luo J. T. ; Tao R. ; Reboud J. ; Wilson R. ; Cooper J. M. ; Gibson D. ; Fu Y. Q. Concentrated vertical jetting mechanism for isotropically focused ZnO/Si surface acoustic waves. Int. J. Multiphase Flow 2019, 114 , 1–8. 10.1016/j.ijmultiphaseflow.2019.02.002.
Tan M. K. ; Friend J. R. ; Yeo L. Y. Interfacial jetting phenomena induced by focused surface vibrations. Phys. Rev. Lett. 2009, 103 , 024501 10.1103/PhysRevLett.103.024501.19659210
Sun D. ; Böhringer K. F. ; Sorensen M. ; Nilsson E. ; Edgar J. S. ; Goodlett D. R. Droplet delivery and nebulization system using surface acoustic wave for mass spectrometry. Lab Chip 2020, 20 , 3269–3277. 10.1039/D0LC00495B.32760973
Cortez-Jugo C. ; Masoumi S. ; Chan P. P. ; Friend J. ; Yeo L. Nebulization of siRNA for inhalation therapy based on a microfluidic surface acoustic wave platform. Ultrason. Sonochem. 2022, 88 , 106088 10.1016/j.ultsonch.2022.106088.35797825
Guo F. ; Mao Z. ; Chen Y. ; Xie Z. ; Lata J. P. ; Li P. ; Ren L. ; Liu J. ; Yang J. ; Dao M. ; Suresh S. ; Huang T. J. Three-dimensional manipulation of single cells using surface acoustic waves. Proc. Natl. Acad. Sci. U.S.A. 2016, 113 , 1522–1527. 10.1073/pnas.1524813113.26811444
Collins D. J. ; Morahan B. ; Garcia-Bustos J. ; Doerig C. ; Plebanski M. ; Neild A. Two-dimensional single-cell patterning with one cell per well driven by surface acoustic waves. Nat. Commun. 2015, 6 , 8686 10.1038/ncomms9686.26522429
Destgeer G. ; Sung H. J. Recent advances in microfluidic actuation and micro-object manipulation via surface acoustic waves. Lab Chip 2015, 15 , 2722–2738. 10.1039/C5LC00265F.26016538
Wang K. ; Zhou W. ; Lin Z. ; Cai F. ; Li F. ; Wu J. ; Meng L. ; Niu L. ; Zheng H. Sorting of tumour cells in a microfluidic device by multi-stage surface acoustic waves. Sens. Actuators, B 2018, 258 , 1174–1183. 10.1016/j.snb.2017.12.013.
Nawaz A. A. ; Soteriou D. ; Xu C. K. ; Goswami R. ; Herbig M. ; Guck J. ; Girardo S. Image-based cell sorting using focused travelling surface acoustic waves. Lab Chip 2023, 23 , 372–387. 10.1039/D2LC00636G.36620943
Liu X. ; Chen X. ; Yang Z. ; Xia H. ; Zhang C. ; Wei X. Surface acoustic wave based microfluidic devices for biological applications. Sens. Diagn. 2023, 2 , 507–528. 10.1039/D2SD00203E.
Guo F. ; Zhou W. ; Li P. ; Mao Z. ; Yennawar N. H. ; French J. B. ; Huang T. J. Precise manipulation and patterning of protein crystals for macromolecular crystallography using surface acoustic waves. Small 2015, 11 , 2733–2737. 10.1002/smll.201403262.25641793
Wang Y. ; Xu C. ; Liu J. ; Pan H. ; Li Y. ; Mei D. Acoustic-assisted 3D printing based on acoustofluidic microparticles patterning for conductive polymer composites fabrication. Addit. Manuf. 2022, 60 , 103247 10.1016/j.addma.2022.103247.
Xu M. ; Lee P. V. ; Collins D. J. Microfluidic acoustic sawtooth metasurfaces for patterning and separation using traveling surface acoustic waves. Lab Chip 2021, 22 , 90–99. 10.1039/D1LC00711D.34860222
Li P. Q. ; Zhou W. ; Peng B. ; Zhang C. ; Zhu X. F. ; Meng L. ; Zheng H. Holographic surface-acoustic-wave tweezers for functional manipulation of solid or liquid objects. Phys. Rev. Appl. 2023, 20 , 064003 10.1103/PhysRevApplied.20.064003.
Maramizonouz S. ; Jia C. ; Rahmati M. ; Zheng T. ; Liu Q. ; Torun H. ; Wu Q. ; Fu Y. Acoustofluidic patterning inside capillary tubes using standing surface acoustic waves. Int. J. Mech. Sci. 2022, 214 , 106893 10.1016/j.ijmecsci.2021.106893.
Wang Q. ; Maramizonouz S. ; Martin M. S. ; Zhang J. ; Ong H. L. ; Liu Q. ; Yang X. ; Rahmati M. ; Torun H. ; Ng W. P. ; Wu Q. ; Binns R. ; Fu Y. Acoustofluidic patterning in glass capillaries using travelling acoustic waves based on thin film flexible platform. Ultrasonics 2024, 136 , 107149 10.1016/j.ultras.2023.107149.37703751
Fu C. ; Quan A. J. ; Luo J. T. ; Pang H. F. ; Guo Y. J. ; Wu Q. ; Ng W. P. ; Zu X. T. ; Fu Y. Q. Vertical jetting induced by shear horizontal leaky surface acoustic wave on 36° YX LiTaO3. Appl. Phys. Lett. 2017, 110 , 173501 10.1063/1.4982073.
Yeo L. Y. ; Friend J. R. Surface acoustic wave microfluidics. Annu. Rev. Fluid Mech. 2014, 46 , 379–406. 10.1146/annurev-fluid-010313-141418.
Ding X. ; Li P. ; Lin S. C. S. ; Stratton Z. S. ; Nama N. ; Guo F. ; Slotcavage D. ; Mao X. ; Shi J. ; Costanzo F. ; Huang T. J. Surface acoustic wave microfluidics. Lab Chip 2013, 13 , 3626–3649. 10.1039/c3lc50361e.23900527
Zhou J. ; Pang H. F. ; Garcia-Gancedo L. ; Iborra E. ; Clement M. ; De Miguel-Ramos M. ; Jin H. ; Luo J. K. ; Smith S. ; Dong S. ; Wang D. M. ; Fu Y. Q. Discrete microfluidics based on aluminum nitride surface acoustic wave devices. Microfluids Nanofluids 2015, 18 , 537–548. 10.1007/s10404-014-1456-1.
Wang Y. ; Zhang Q. ; Tao R. ; Chen D. ; Xie J. ; Torun H. ; Dodd L. E. ; Luo J. ; Fu C. ; Vernon J. ; Canyelles-Pericas P. ; Binns R. ; Fu Y. A rapid and controllable acoustothermal microheater using thin film surface acoustic waves. Sens. Actuators, A 2021, 318 , 112508 10.1016/j.sna.2020.112508.
Zhou J. ; Tao X. ; Luo J. ; Li Y. ; Jin H. ; Dong S. ; Luo J. ; Duan H. ; Fu Y. Nebulization using ZnO/Si surface acoustic wave devices with focused interdigitated transducers. Surf. Coat. Technol. 2019, 367 , 127–134. 10.1016/j.surfcoat.2019.03.078.
Fu Y. Q. ; Pang H. F. ; Torun H. ; Tao R. ; McHale G. ; Reboud J. ; Tao K. ; Zhou J. ; Luo J. ; Gibson D. ; Luo J. ; Hu P. Engineering inclined orientations of piezoelectric films for integrated acoustofluidics and lab-on-a-chip operated in liquid environments. Lab Chip 2021, 21 , 254–271. 10.1039/D0LC00887G.33337457
Tao R. ; Reboud J. ; Torun H. ; McHale G. ; Dodd L. E. ; Wu Q. ; Tao K. ; Yang X. ; Luo J. T. ; Todryk S. ; Fu Y. Integrating microfluidics and biosensing on a single flexible acoustic device using hybrid modes. Lab Chip 2020, 20 , 1002–1011. 10.1039/C9LC01189G.32026889
Hsu J. C. ; Chao C. L. Acoustophoretic patterning of microparticles in a microfluidic chamber driven by standing Lamb waves. Appl. Phys. Lett. 2021, 119 , 103504 10.1063/5.0056985.
Lei J. ; Zheng G. ; Cheng F. ; Li K. ; Huang Z. Acoustofluidics in an equilateral triangular channel. Exp. Therm. Fluid Sci. 2024, 151 , 111088 10.1016/j.expthermflusci.2023.111088.
Zhou J. ; Guo Y. ; Wang Y. ; Ji Z. ; Zhang Q. ; Zhuo F. ; Luo J. ; Tao R. ; Xie J. ; Reboud J. ; McHale G. ; Dong S. ; Luo J. ; Duan H. ; Fu Y. Flexible and wearable acoustic wave technologies. Appl. Phys. Rev. 2023, 10 , 021311 10.1063/5.0142470.
Tao X. ; Nguyen T. D. ; Jin H. ; Tao R. ; Luo J. ; Yang X. ; Torun H. ; Zhou J. ; Huang S. ; Shi L. ; Gibson D. ; Cooke M. ; Du H. ; Dong S. ; Luo J. ; Fu Y. Q. 3D patterning/manipulating microparticles and yeast cells using ZnO/Si thin film surface acoustic waves. Sens. Actuators, B 2019, 299 , 126991 10.1016/j.snb.2019.126991.
Pang H. F. ; Tao R. ; Luo J. ; Zhou X. ; Zhou J. ; McHale G. ; Reboud J. ; Torun H. ; Gibson D. ; Tao K. ; Chang H. ; Fu Y. Q. Rayleigh and shear-horizontal surface acoustic waves simultaneously generated in inclined ZnO films for acoustofluidic lab-on-a-chip. Surf. Coat. Technol. 2022, 442 , 128336 10.1016/j.surfcoat.2022.128336.
Huang J. ; Wang Z. ; Liu R. ; Zhu Z. ; Zhang D. ; Guo X. ZnO/glass-based SAW tweezer for dexterous particle patterning and patterned cell culturing. Microfluids Nanofluids 2023, 27 , 34 10.1007/s10404-023-02646-3.
Ji Z. ; Zhou J. ; Lin H. ; Wu J. ; Zhang D. ; Garner S. ; Gu A. ; Dong S. ; Fu Y. ; Duan H. Flexible thin-film acoustic wave devices with off-axis bending characteristics for multisensing applications. Microsyst. Nanoeng. 2021, 7 , 97 10.1038/s41378-021-00325-3.34900331
Lamanna L. Recent Progress in Polymeric Flexible Surface Acoustic Wave Devices: Materials, Processing, and Applications. Adv. Mater. Technol. 2023, 8 , 2300362 10.1002/admt.202300362.
Jin H. ; Zhou J. ; He X. ; Wang W. ; Guo H. ; Dong S. ; Wang D. ; Xu Y. ; Geng J. ; Luo J. K. ; Milne W. I. Flexible surface acoustic wave resonators built on disposable plastic film for electronics and lab-on-a-chip applications. Sci. Rep. 2013, 3 , 1–8. 10.1038/srep02140.
Tao R. ; McHale G. ; Reboud J. ; Cooper J. M. ; Torun H. ; Luo J. ; Luo J. ; Yang X. ; Zhou J. ; Canyelles-Pericas P. ; Wu Q. ; Fu Y. Hierarchical nanotexturing enables acoustofluidics on slippery yet sticky, flexible surfaces. Nano Lett. 2020, 20 , 3263–3270. 10.1021/acs.nanolett.0c00005.32233442
Zeng X. ; Yan Z. ; Lu Y. ; Fu Y. ; Lv X. ; Yuan W. ; He Y. Reduction of ice adhesion using surface acoustic waves: nanoscale vibration and interface heating effects. Langmuir 2021, 37 , 11851–11858. 10.1021/acs.langmuir.1c01852.34585928
Zeng X. ; Ong H. ; Haworth L. ; Lu Y. ; Yang D. ; Rahmati M. ; Wu Q. ; Torun H. ; Martin J. ; Hou X. ; Lv X. ; Yuan W. ; He Y. ; Fu Y. Fundamentals of Monitoring Condensation and Frost/Ice Formation in Cold Environments Using Thin-Film Surface-Acoustic-Wave Technology. ACS Appl. Mater. Interfaces 2023, 15 , 35648–35663. 10.1021/acsami.3c04854.37432769
Liu Y. ; Li Y. ; el-Hady A. M. ; Zhao C. ; Du J. F. ; Liu Y. ; Fu Y. Q. Flexible and bendable acoustofluidics based on ZnO film coated aluminium foil. Sens. Actuators, B 2015, 221 , 230–235. 10.1016/j.snb.2015.06.083.
Tao R. ; Wang W. B. ; Luo J. T. ; Hasan S. A. ; Torun H. ; Canyelles-Pericas P. ; Zhou J. ; Xuan W. P. ; Cooke M. D. ; Gibson D. ; Wu Q. ; Ng W. P. ; Luo J. K. ; Fu Y. Q. Thin film flexible/bendable acoustic wave devices: evolution, hybridization and decoupling of multiple acoustic wave modes. Surf. Coat. Technol. 2019, 357 , 587–594. 10.1016/j.surfcoat.2018.10.042.
Wang Y. ; Qian J. Femtosecond Laser Micromachining of the Mask for Acoustofluidic Device Preparation. ACS Omega 2023, 8 , 7838–7844. 10.1021/acsomega.2c07589.36873004
Guo Y. J. ; Lv H. B. ; Li Y. F. ; He X. L. ; Zhou J. ; Luo J. K. ; Zu X. T. ; Walton A. J. ; Fu Y. Q. High frequency microfluidic performance of LiNbO3 and ZnO surface acoustic wave devices. J. Appl. Phys. 2014, 116 , 024501 10.1063/1.4885038.
Nama N. ; Barnkob R. ; Mao Z. ; Kähler C. J. ; Costanzo F. ; Huang T. J. Numerical study of acoustophoretic motion of particles in a PDMS microchannel driven by surface acoustic waves. Lab Chip 2015, 15 , 2700–2709. 10.1039/C5LC00231A.26001199
Hsu J.-C. ; Chao C.-L. Full-wave modeling of micro-acoustofluidic devices driven by standing surface acoustic waves for microparticle acoustophoresis. J. Appl. Phys. 2020, 128 , 124502 10.1063/5.0017933.
Qi M. ; Dang D. ; Yang X. ; Wang J. ; Zhang H. ; Liang W. Surface acoustic wave manipulation of bioparticles. Soft Matter 2023, 19 , 4166 10.1039/D3SM00457K.37212436
Wiklund M. ; Green R. ; Ohlin M. Acoustofluidics 14: Applications of acoustic streaming in microfluidic devices. Lab Chip 2012, 12 , 2438–2451. 10.1039/c2lc40203c.22688253
Meng L. ; Cai F. ; Li F. ; Zhou W. ; Niu L. ; Zheng H. Acoustic tweezers. J. Phys. D: Appl. Phys. 2019, 52 , 273001 10.1088/1361-6463/ab16b5.
Yeo L. ; Rezk A. Lithography-free, crystal-based multiresonant lamb waves for reconfigurable microparticle manipulation. Inf. MIDEM 2016, 46 , 176–182.
