Intra-droplet acoustic particle focusing: simulations and experimental observations

The aim of this paper is to study resonance conditions for acoustic particle focusing inside droplets in two-phase microfluidic systems. A bulk acoustic wave microfluidic chip was designed and fabricated for focusing microparticles inside aqueous droplets (plugs) surrounded by a continuous oil phase in a 380-μm-wide channel. The quality of the acoustic particle focusing was investigated by considering the influence of the acoustic properties of the continuous phase in relation to the dispersed phase. To simulate the system and study the acoustic radiation force on the particles inside droplets, a simplified 3D model was used. The resonance conditions and focusing quality were studied for two different cases: (1) the dispersed and continuous phases were acoustically mismatched (water droplets in fluorinated oil) and (2) the dispersed and continuous phases were acoustically matched (water droplets in olive oil). Experimentally, we observed poor acoustic particle focusing inside droplets surrounded by fluorinated oil while good focusing was observed in droplets surrounded by olive oil. The experimental results are supported qualitatively by our simulations. These show that the acoustic properties (density and compressibility) of the dispersed and continuous phases must be matched to generate a strong and homogeneous acoustic field inside the droplet that is suitable for high-quality intra-droplet acoustic particle focusing.


Introduction
Manipulation of particles in microfluidic channels using acoustic forces, acoustophoresis, has shown to be a useful tool in a number of lab-on-a-chip applications including the separation of rare cells (Nordin and Laurell 2012;Li et al. 2015), concentration of bacteria (Carugo et al. 2014) and cell trapping (Hammarström et al. 2010;Christakou et al. 2015;Collins et al. 2015). Recently, acoustic particle manipulation has also been implemented in droplet-based microfluidic systems for manipulation of whole droplets (Schmid et al. 2014;Leibacher et al. 2015;Sesen et al. 2015) and focusing of particles and cells inside droplets (Fornell et al. 2015Park et al. 2017).
The interest in droplet microfluidics has emerged because droplets are perfectly suited as miniaturized reaction chambers for encapsulation and investigation of cells at the singlecell level (Schneider et al. 2013;Shembekar et al. 2016), and the technology allows for improved genome engineering and screening of cells (Agresti et al. 2010;Sjostrom et al. 2014). In most of these applications, having the possibility to precisely handle the particles inside the droplets expands the complexity of the assays that can be performed on-chip. Particle manipulation allows, for example, intra-droplet particle enrichment and buffer exchange which are often required in multistep assays .
To control the position of particles inside droplets both active and passive methods have been implemented including magnetophoresis (Lombardi and Dittrich 2011;Lee et al. 2014;Brouzes et al. 2015), dielectrophoresis (Han et al. 2017), acoustophoresis (Fornell et al. 2015Park et al. 2017) and hydrodynamic methods (Kurup and Basu 2012;Sun et al. 2012;Hein et al. 2015). Of these, acoustophoresis has the advantages of being label free, gentle, and operated in non-contact mode. However, in our previous work, there have been experimental indications that the acoustic particle focusing inside droplets is impaired by a difference in acoustic properties between the droplet and the continuous phase (Fornell et al. 2015). Typically, water-based solutions are used as the dispersed phase as living cells are often encapsulated inside the droplets, and fluorinated oils such as Novec HFE-7500 or FC-40 with the addition of fluorosurfactants are used as the continuous phase (also known as the carrier oil). These oils are preferred in droplet-based applications as they allow for the generation of stable water droplets and low cross-contamination between individual droplets, as well as having high gas solubility which is essential for cell viability and proliferation (Baret 2012). However, these oils have very different acoustic properties compared with water, and that is thought to have a negative impact on intra-droplet acoustic particle focusing in two-phase microfluidic systems.
The interaction of sound waves and liquid interfaces has been studied for a long time in both micro-and macrosystems. Almost 80 years ago, the deformation of the interface between two immiscible liquids with different acoustic properties caused by a directed ultrasonic beam was studied (Hertz and Mende 1939). Additionally, the acoustic forces can also act on liquid interfaces in miscible systems, and it has been shown that laminated liquids with different acoustic properties exposed to a standing wavefield can be relocated or stabilized within a microfluidic channel (Deshmukh et al. 2014).
Motivated by these observations and studies, this paper investigates acoustic particle focusing in two-phase systems and in particular how the acoustic properties of the continuous phase influence the acoustic radiation force on particles encapsulated inside droplets. The study is restricted to two cases: water droplets (plugs) generated in fluorinated oil and vegetable olive oil, respectively. Fluorinated oil is chosen since this oil type is commonly used in droplet microfluidic systems, and olive oil is chosen since it has similar acoustic properties as water. First, a simplified 3D model is set up and the acoustic pressure field and the acoustic radiation force on the encapsulated particles is calculated. Second, acoustic particle focusing is analyzed experimentally. This study provides a wider understanding of the fundamental physical principles of intra-droplet acoustic particle focusing, which is essential for the development of optimal system design for acoustic particle manipulation in two-phase microfluidic systems.

Acoustophoresis
Ultrasonic standing waves can be used to position particles in microfluidic channels. For particles larger than 2 µm suspended in water, the acoustic boundary layer and viscosity in the MHz range can be neglected. In this case, the acoustic velocity field v 1 is given by the gradient of the acoustic pressure field p 1 as v 1 = −i wa ∇p 1 , and the acoustic radiation force F rad acting on the microparticle is a gradient force given by (Bruus 2012) Here is the angular frequency, f 0 and f 1 are the monopole and dipole coefficients, a P , P , and P are the radius, density and compressibility of the microparticle, and WA and WA are the density and compressibility of water, respectively. As seen in Eq. (1), the acoustic properties of the particle relative to the surrounding water determine the direction of the radiation force on the particle: a particle with high density and low compressibility is focused to the pressure node whereas a particle with low density and high compressibility is focused to the pressure anti-node. Both plastic particles and cells in water have such acoustic properties that they are focused to the pressure node in an acoustic standing wavefield. The focusing time t foc , defined as the time it takes a particle to move from y = − 7 16 W to y = − 1 16 W or ∼ 40% of the half channel width W , is given in terms of the particle radius a p , viscosity , and the acoustic energy density E ac as (Bruus 2012)

Model system used in the numerical simulations
The system is modelled by considering a water-in-oil droplet confined in a rectangular channel of height H = 100 µm and width W = 380 µm in an acoustically hard silicon-Pyrex chip to which the ultrasound is supplied by oscillating the bottom plane with a frequency near 2 MHz, see Fig. 1.
While the width W and height H of the channel as well as the height of the silicon and Pyrex layers equal the experimental values, the width W s of the solid chip and the length L s of the system has been reduced to enable the demanding 3D simulation on the available 128-GB RAM computer. A further reduction is achieved by limiting the study to systems with rectangular channels that are symmetric around the vertical y-z-plane, such that only half of the system needs to be simulated. Moreover, we neglect the thin film of the continuous phase between the droplet and the channel walls (Baroud (1) , and f 1 = 2( P − WA ) 2 P + WA .
(2) t foc = 9 . The material and geometrical parameters used in the simulations are listed in Table 1.

Numerical simulations: method
Simulations of the systems are performed using the finiteelement solver COMSOL Multiphysics (version 5.2) using the method of Ley and Bruus (Ley and Bruus 2017) for Pyrex-water systems, extended to include the single-crystal silicon base using the method of Dual and Schwarz (Dual and Schwarz 2012). The acoustic field is calculated numerically using the linear elastic equations for the displacement field u 1 in the solid wall with density s , damping coefficient Γ and stress tensor 1 coupled to the pressure field p 1 governed by the Helmholtz wave equation in the two fluids (a) and (b), To mimic the non-symmetric actuation of the experiment on the bottom of the silicon base, a harmonically oscillating vertical displacement of the form u 1 = d 1 .4 nm and using complex-time notation. Here = 2 f , where and f are the angular frequency and frequency, respectively. On the vertical y-z-planes at x = 0 and x = L, we use symmetry boundary conditions (hard wall conditions), and all other outer boundaries have the no-stress condition. On the internal solid-liquid and liquid-liquid interfaces, the boundary conditions are continuity of the normal component of the stress and of the displacement velocities. The corresponding tangential stresses are set to zero. In summary, where n is the interface normal vector. We performed a standard numerical convergence study (Ley and Bruus 2017) and found the fields to have converged within a relative deviation of 1%.

(4a)
Solid − fluid interface ∶ 1 ⋅ n = −p 1 n (stress) and Solid − air boundary ∶ 1 ⋅ n = 0,  (beige). The device is assumed to be symmetric around the y-z-plane, such that only half of the system ( x > 0 ) needs to be simulated. b The detailed geometry of the device is shown with the silicon layer opened up to make the water droplet (cyan) visible in the long, straight, rectangular oil-filled channel (blue). c The detailed geometry of the flattened water droplet. d End view in the symmetry plane of the water droplet (cyan) and surrounding oil (blue). The bottom plane is actuated vertically as

frequency near 2 MHz
Because the size of the chip in the numerical model is smaller than the actual chip, we first characterize the system with air in the channel to locate the resonances determined solely by the geometry of the solid. In the range from 1.5 to 2.1 MHz we find only two resonances, namely at 1.60 and at 2.01 MHz, as shown in Fig. 2, where the black line shows the average acoustic energy density E solid ac versus frequency f . In the subsequent analysis, with liquid in the channel, we then make sure to avoid these resonances by working in the range from 1.64 to 1.96 MHz. Our numerical results are thus independent of the specific smaller chip geometry of the model and can be applied to the larger actual chip geometry. The optimal frequency for acoustophoresis inside the droplet is found by locating resonance peaks in plots of the acoustic radiation force F rad versus frequency f for a 10-µm-diameter polystyrene test particle averaged over its position in the water droplet, see Eq. (1). Ideally, the vertical component should vanish, F rad z = 0, while the horizontal component should be anti-symmetric and pointing towards the x-z-plane with as big as magnitude as possible, sign(y) F rad y < 0 , where the sign of the y-coordinate, sign(y) , is inserted to pick up the anti-symmetry of the radiation force component in the spatial average. vs. f in the solid walls for the empty chip (black line), and the averaged y-and z-components F rad y,z of the acoustic radiation force on a 10-µm-diameter polystyrene test particle in the water droplet in fluorinated oil ( F rad y,HFE blue and F rad z,HFE cyan), and olive oil ( F rad y,OLV red and F rad z,OLV magneta). The resonances are found to be at f =1.718 MHz for fluorinated oil, and at f = 1.865, 1.878, 1.905, and 1.950 MHz for olive oil. The latter resonance series corresponds to the number of axial nodes being 0, 1, 2, and 3

Numerical simulations: results and discussion
The resonance curves for a water droplet surrounded by fluorinated oil (marked as HFE) and olive oil (marked as OLV) are shown in Fig. 2. Strong resonances are found at f = 1.718 MHz with sign(y) F rad y = − 1.9 pN for a water droplet in fluorinated oil, and at f = 1.865 MHz with sign(y) F rad y = − 5.4 pN for a water droplet in olive oil. The vertical components are negligible in the entire frequency range.
A physical interpretation of the resonances and their ability to facilitate acoustophoresis is obtained by studying the acoustic pressure p in the fluids and the displacement u in the solids as well as the acoustic radiation force F rad in the water droplet for the main resonances at f = 1.718 MHz (fluorinated oil) and at f = 1.865 MHz (olive oil) shown in Fig. 3a, b, respectively.
We first note that for an ideal hard-walled system of width W = 380 µm, a standing one-half wavelength resonance should occur at f WA = c WA ∕(2W) = 1.970 MHz in water and at f OLV = c OLV ∕(2W) = 1.907 MHz in olive oil where c is the sound speed, while a standing two-half wavelength resonance should occur at f HFE = c HFE ∕W = 1.734 MHz in fluorinated oil. We have tested this simple prediction with simulations on the system without a water droplet but filled uniformly with either water, fluorinated oil, or olive oil and found the following resonance frequencies numerically: f num WA = 1.918 MHz, f num HFE = 1.722 MHz, and f num OLV = 1.864 MHz. For a water droplet in fluorinated oil, Fig. 3a, we first of all note that the numerically obtained resonance frequency f = 1.718 MHz is close to the pure fluorinated oil two-half wavelength resonance of f num HFE = 1.722 MHz. However, we also observe a clear mismatch between the approximate one-half and two-half standing pressure wavelengths in the water droplet and in the fluorinated oil, respectively. This mismatch leads to a zero pressure and a zero radiation force near the front end of the droplet. Not only is the average magnitude of the radiation force lowered by the acoustic mismatch between the water and the surrounding oil, but the suppression of the radiation force near the front end implies that the radiation force cannot counterbalance the drag forces from the internal circulation flow roll that occurs as the droplet moves inside that channel, forcing the particles away from the center plane ). In conclusion, the simulation predicts low-quality acoustophoresis in water droplet surrounded by fluorinated oil. In contrast, the acoustic mismatch between water and olive oil is minute, and we see in Fig. 3b, a nearperfect standing one-half wavelength resonance throughout the fluid domain at f = 1.865 MHz, close to the pure olive oil one-half wave resonance f num OLV = 1.864 MHz, when the water droplet is surrounded by olive oil. As a consequence, there is no suppression of the radiation force near the droplet front end, and the internal flow rolls cannot disperse the particles away from the center axis. In conclusion, the simulation predicts high quality of acoustophoresis in a water droplet surrounded by olive oil.

Fabrication of the device
A silicon-glass microfluidic chip for droplet generation, particle encapsulation, and intra-droplet acoustic particle focusing is fabricated, and in Fig. 4a the chip design is shown schematically.
The microchannels are fabricated using standard photolithography and deep reactive-ion etching on a silicon wafer (500 µm thick, ⟨100⟩-orientation). Deep reactive-ion etching Fig. 3 a The chip with a water droplet in fluorinated oil at resonance f = 1.718 MHz. The acoustic pressure p from − 156 kPa (blue) to 156 kPa (red) and the magnitude of the acoustic displacement u from 0 nm (black) to 0.35 nm (white) in the solids. The inset shows the acoustic radiation force F rad inside the water droplet with a magnitude from 0 pN (black) to 4.5 pN (white). b Same but for olive oil at resonance f = 1.865 MHz and with amplitudes for p from − 448 kPa (blue) to 447 kPa (red), for u from 0 nm (black) to 0.35 nm (white), and for F rad from 0 pN (black) to 10.7 pN (white) is a highly anisotropic etching method that yields channels with rectangular cross sections. The width of the channels is 380 µm and the height is 100 µm. Holes for fluid inlets and outlets are drilled through the silicon wafer. The channels are sealed by anodic bonding of a glass wafer (700 µm thick, Borofloat-33), and the silicon-glass wafer is diced into individual chips. Short pieces of silicone tubings are glued to the chip as fluid connectors. The channels are silanized by flushing the channels with Repel-Silane (Pharmacia Biotech) to make the channels hydrophobic, which is a requirement for generating water-in-oil droplets. To supply the ultrasound, a 12 mm × 12 mm piezoelectric transducer (1 mm thick yielding 2 MHz fundamental resonance frequency, Pz26, Ferroperm Piezoceramics A/S) is glued on the silicon side of the chip using cyanoacrylate glue (Loctite 420, Henkel AG & Co.), and electrical wires are soldered on the piezoelectric transducer.

Experimental method
The experimental setup is shown in Fig. 4b. The piezoelectric transducer is actuated by a sinusoidal signal from a function generator (33220A, Agilent Technologies Inc.) after amplification of the signal (75A250, Amplifier Research). The output voltage from the function generator is kept constant in the experiments and the voltage over the transducer is monitored using a digital oscilloscope (TDS 1002, Tektronix). The voltage in the experiments over the transducer is between 10 and 25 V peak−peak depending on the exact frequency of the ultrasound. The frequency is manually swept between 1.50 and 2.50 MHz in steps of 0.05 MHz to find the strongest particle focusing at the fundamental resonance for each system. The focusing strength is determined by the operator by visually inspecting the system and identifying when the particles are focused using an optical microscope (BX51W1, Olympus) equipped with a 4× objective. Images are acquired with a camera (XM10, Olympus) mounted on the microscope.
The fluid flows are controlled by three syringe pumps (NEMESYS, Cetoni GmbH) mounted with plastic syringes and connected to the inlet channels via Teflon tubing whereas the outlet channel is kept open. In the experiments, the dispersed phase is Milli-Q water. In the acoustically mismatched system, fluorinated oil (Novec HFE-7500, 3M) containing 2% surfactant (Krytox FSH-157, Dupont) is used as the continuous phase, and in the acoustically matched droplet system vegetable olive oil (Di Luca and Di Luca) is used as the continuous phase. Microparticles (10-µm polystyrene particles, Sigma-Aldrich) are suspended in the water phase, and in the droplet generation process the particles are encapsulated inside the droplets. The total flow rate is kept constant in the experiments (5 µL/min) to ensure the residence time of the droplet in the channel is the same. However, the length of the droplets is known to be influenced by the material properties of the fluid phases (Utada et al. 2007;Chen et al. 2014), and in the experiment this effect is partly adjusted for by changing the water:oil flow ratio. For generation of water droplets in olive oil, the flow ratio is set to 2:3 compared with 1:4 for generation of water droplets in fluorinated oil.

Experimental results and discussion
In the experiments, water droplets were generated in fluorinated oil and olive oil, respectively, and polystyrene particles were encapsulated inside the droplets. Without actuation of the ultrasound, the particles were positioned in the entire droplet, see Fig. 5a, c. To focus the encapsulated particles, ultrasound at the fundamental resonance frequency was applied. For each system, the frequency was tuned manually to find the strongest particle focusing frequency. To minimize the influence of variations in the acoustic field arising from external factors such as the gluing of the transducer and differences from the microfabrication process, the same microfluidic chip and the same part of the channel was studied throughout the experimental series.
For water droplets in fluorinated oil no strong focusing was observed when the ultrasound was applied; however, it was noted that at f = 1.70 − 1.75 MHz the particles were weakly affected by the ultrasound, see Fig. 5b. In contrast, when the water droplets were surrounded by olive oil, it was observed that the encapsulated particles were strongly  Fig. 5d. The strongest particle focusing was seen at f = 1.70 − 1.90 MHz for water droplets in olive oil, but some focusing could be observed at higher frequencies, see Table 2. The fact that the observed frequency range for particle focusing in water droplets surrounded by olive oil was wide indicates that this system is robust and not particularly sensitive to small fabrication and temperature variations, and thus easy and stable to operate.

Concluding discussion
The numerical simulations and experiments show the importance of matching the acoustic properties of the dispersed and the continuous phase to achieve good quality intra-droplet acoustic particle focusing. Experimentally good particle focusing was only observed for water droplets in olive oil whereas for water droplets in fluorinated oil the significant difference in the acoustic properties between the two phases resulted in weak particle focusing. This can be explained qualitatively from the simulations, which clearly show that strong and uniform focusing occurs in systems with good acoustic matching between the droplet and the continuous phase (olive oil), while focusing is only maintained near the center of droplets having a poor acoustic matching with the continuous phase (fluorinated oil). The lack of focusing forces near the droplet ends is crucial in explaining the poor overall acoustophoretic behavior: as droplets move in microfluidic channels flow rolls are induced inside them (Kinoshita et al. 2007;Ma et al. 2014), and these flow rolls counteract the acoustic focusing force in the front of the droplet . Consequently, in acoustically mismatched droplets, the acoustic focusing is too weak to balance the defocusing Stokes' drag force from the flow rolls, and such droplets are unlikely to exhibit good acoustic particle focusing. Furthermore, in actual experiments, the acoustic fields are never completely perfect and homogeneous, but instead there are areas in the channel with a higher or lower acoustic energy ("hot spots" and "weak spots") due to periodical variations along the length of the channel (Barnkob et al. 2010) as well as bad acoustic coupling, and here the internal flow rolls can have a larger impact and spread the particles more. To quantify the influence from the flow rolls, we note that the fluid flow rolls affect the particles through the Stokes drag force F drag = 6 a (v d − v p ) , where v d is the fluid velocity in the droplet and v p is the particle velocity. In our experiments, the flow rate was 5 µL/min, which corresponds to an average fluid velocity of 2 mm/s. In Fig. 5 Photographs of acoustic focusing of microparticles inside water droplets in fluorinated oil and olive oil, respectively. At actuation of the ultrasound at the fundamental resonance, a standing wavefield is created between the channel walls, and the particles are affected by the radiation force. The particles are most strongly focused in the water droplet surrounded by olive oil. The droplets are flowing towards the right in the photographs  (Baroud et al. 2010). We, therefore, estimate the drag force on the particles in our experiment to be of the order F drag ≈ 2 − 20 pN . Following Barnkob et al. (2010), the magnitude of the radiation force in our experiments is determined by measuring the typical focusing time t foc = 1 s , which by Eq.
(2) gives a typical energy density E ac = 6 Pa . In the simulations, this value of the energy density is obtained by choosing the actuation amplitude to be d 1 = 0.4 nm, a value for which the resulting radiation force becomes F rad ≈ 11 pN , or somewhat smaller, as shown in Figs. 2 and 3. Thus, it is within the experimental control of the drag force through the flow rate and of the radiation force through the actuation voltage on the piezocrystal to cross over from a radiation force-dominated to drag force-dominated behavior. However, for a given setting of the control parameters, the focusing ability of the acoustically matched olive oil system is better than for the mismatched fluorinated oil system. Because the typical timescales for the acoustic focusing of the particles, for the passage of the droplets through the channel, and for the sedimentation of the particles in our experiments are 1, 5, and 32 s, respectively, it follows that the influence of gravity is negligible.
The purpose of this study was to investigate how the acoustic properties of the continuous phase affect particle focusing in acoustofluidic two-phase systems; however, from the application viewpoint in addition to suitable acoustic properties, care must also be taken to find fluid phases that provide monodisperse droplet generation, high droplet stability, no cross-contamination and are biocompatible.

Conclusion
In this work, acoustic focusing of microparticles encapsulated in water-in-oil droplets was studied. Two different systems were considered: one where there was a significant difference in the acoustic properties of the dispersed and continuous phase (water droplets in fluorinated oil) and one where the phases had similar acoustic properties (water droplets in olive oil). The simulations and experimental results show that to obtain high-quality acoustic focusing inside the droplets the acoustic properties of the dispersed and continuous phase should be matched. We believe that these findings will provide important information in designing experiments regarding acoustofluidics in two-phase systems.