The interaction between an acoustically driven microbubble and a surface is of interest for a variety of applications, such as ultrasound imaging and therapy. Prior investigations have mainly focused on acoustic effects of a rigid boundary, where it was generally observed that the wall increases inertia and reduces the microbubble resonance frequency. Here we investigate the response of a lipid-coated microbubble adherent to a rigid wall. Firm adhesion between the microbubble and a glass surface was achieved through either specific (biotin/avidin) or nonspecific (lipid/glass) interactions. Total internal reflection fluorescence microscopy was used to verify conditions leading to either adhesion or non-adhesion of the bubble to a glass or rigid polymer surface. Individual microbubbles were driven acoustically to sub-nanometer-scale radial oscillations using a photoacoustic technique. Remarkably, adherent microbubbles were shown to have a higher resonance frequency than non-adherent microbubbles resting against the wall. Analysis of the resonance curves indicates that adhesion stiffens the bubble by an apparent increase in the shell elasticity term and decrease in the shell viscosity. Based on these results, we conclude that surface adhesion is dominant over acoustic effects for low-amplitude microbubble oscillations.

Microbubbles are currently used in biomedical ultrasonics for contrast-enhanced ultrasound imaging, super-resolution imaging, molecular imaging, and image-guided therapy. Biomedical microbubbles comprise a 1–10 μm diameter gas core suspended in an aqueous medium and commonly coated with a phospholipid monolayer shell.1 The resonance frequency of a lipid-coated microbubble is typically between 1 and 10 MHz, depending on the diameter and shell stiffness.2 Microbubble resonance fortuitously is in the range of most clinical ultrasound imaging and therapy devices. Optimal sensing and actuation of the microbubble requires matching of the transducer bandwidth and the ultrasound pulse center frequency to the microbubble resonance frequency.3 

Prior research has primarily focused on characterizing microbubble resonance in an isotropic medium.4,5 However, the interaction between an acoustically driven microbubble and a surface is important for many applications. There have been some studies comparing the acoustic response and stability of individual microbubbles near a boundary.6–10 In each of these studies, microbubbles were driven with sufficiently large acoustic pressure to induce a nonlinear response and microbubble instability. The response of adherent microbubbles under stable, linear oscillations, although fundamental to understanding the physical process, has received comparatively little attention. In this letter, we investigate the dynamics of microbubbles adherent to a rigid substrate under stable, linear oscillation using a photoacoustic technique. A modulated laser source was used to generate ultrasound in the liquid adjacent to a microbubble to drive oscillations. The resulting radial displacement, on the order of tens of picometers, was detected using forward light scattering. We show that adhesion to a rigid substrate significantly increases the resonance frequency over that of non-adherent microbubbles resting against the wall. The response is consistent with an increase in stiffness, rather than an increase in mass, and can be captured quantitatively through an increase in the apparent shell elasticity. Bound microbubbles also show a decrease in both the damping ratio (friction) and the maximum radial displacement.

In this study, phospholipid-coated microbubbles were generated by shaking.11 The microbubble shells comprised 90 mol. % of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 10 mol. % of 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[methoxy(polyethylene glycol)-2000] (DSPE-PEG2000). Biotinylated microbubbles were prepared by substituting 10 mol. % of their PEGylated lipid (1 mol% overall) with 1,2-distearoyl-sn-glycerol-3-phosphoethanolamine-N-[biotinyl-(polyethylene glycol)-2000] (DSPE-PEG2000-B). These lipids were suspended in phosphate buffered saline (PBS) at a lipid concentration of 2 mg/ml and were mixed and heated to 45 °C by probe sonication. A 3-ml glass serum vial was filled with lipid suspension and sealed with a gas headspace containing perfluorobutane (PFB). Microbubbles were generated by vigorously shaking the serum vial for 40 s using a dental amalgamator, and the vial was quickly quenched to room temperature in an ice bath. A centrifugation technique12 was used to wash and separate the microbubbles by removing undesired bubble sizes and excess lipid in the solution. Microbubbles were immediately diluted with PBS (approximately 5 × 105 bubbles per ml and 4–10 μm in diameter) and pipetted between a microscope slide and a modified glass coverslip with a thin gasket separating them.

We explored adherent and non-adherent microbubbles, where adhesion resulted from either specific (biotin/avidin) or nonspecific (lipid/glass) interactions. For specific binding, biotinylated microbubbles were brought into contact with a streptavidin-coated glass coverslip, which had been incubated with 33 μg/ml streptavidin for one hour, rinsed with PBS, and then air-dried immediately prior to the experiment. This condition was compared to a non-adherent microbubble condition, where non-biotinylated microbubbles were brought into contact with the streptavidin-coated glass coverslip surface. For nonspecific binding, non-biotinylated microbubbles were brought into contact with a clean glass coverslip. This condition was compared to a non-adherent condition, where non-biotinylated microbubbles were brought into contact with a glass coverslip surface where a supported lipid bilayer was deposited to block adherence. The glass coverslip had been incubated with 2 mg/ml DPPC for one hour, rinsed with PBS and then air-dried immediately prior to the experiment. Microbubbles were allowed to equilibrate with the glass coverslips for 30 min prior to testing.

Total internal reflection fluorescence (TIRF) and epifluorescence microscopy were used to characterize the adhesion between the microbubbles and the glass coverslips. Here, microbubbles were fluorescently labeled with DiO (ex/em = 484/501 nm) during preparation. Figure 1(a) shows a schematic of the prism-TIRF microscopy setup, which incorporated a Nikon inverted microscope with a 60X objective. A 491-nm excitation laser source was directed through an index matched prism to illuminate the sample surface at an incident angle greater than the critical angle for total internal reflection at the glass-water interface. This generated an evanescent field extending from the glass-water interface (∼100 nm deep) where only the fluorophores in the microbubble shell near the surface (i.e., in the region associated with adhesion contact) were excited. The emitted light from the microbubble contact region was then directed with a dichroic mirror (505 nm) and captured by an electron-multiplying CCD camera. For epifluorescence microscopy, the microbubbles were illuminated through the microscope objective using a 488/10 nm bandpass excitation filter, and fluorescent images were captured with the focal plane at the microbubble mid-plane (i.e., the widest projection of the bubble) to quantify the resting radius.

FIG. 1.

(a) Diagram of the prism-TIRF microscope setup showing fluorescently labeled microbubbles at the glass coverslip surface. A totally internally reflected excitation laser source (491 nm) generates an evanescent field extending from the coverslip surface, and the emitted light from the fluorescent microbubbles (dashed green lines) is captured with a camera. TIRF and epifluorescence microscopy images of typical microbubbles for the (b) specific (biotin/avidin) and (c) nonspecific (lipid/glass) conditions. Scale bars are 2 μm. Also shown are contact radius vs resting radius plots derived from the fluorescence images. Solid curves show best linear fit to the data with the intercept fixed at the origin.

FIG. 1.

(a) Diagram of the prism-TIRF microscope setup showing fluorescently labeled microbubbles at the glass coverslip surface. A totally internally reflected excitation laser source (491 nm) generates an evanescent field extending from the coverslip surface, and the emitted light from the fluorescent microbubbles (dashed green lines) is captured with a camera. TIRF and epifluorescence microscopy images of typical microbubbles for the (b) specific (biotin/avidin) and (c) nonspecific (lipid/glass) conditions. Scale bars are 2 μm. Also shown are contact radius vs resting radius plots derived from the fluorescence images. Solid curves show best linear fit to the data with the intercept fixed at the origin.

Close modal

Figures 1(b) and 1(c) show representative epifluorescence and TIRF microscopy images of typical microbubbles for each of the different adhesion conditions. The bubble diameter was quantified using the averaged outer full-width-at-half-maximum taken from the vertical and horizontal pixel intensity line profiles across the center of the epifluorescence image. Similarly, the diameter of the contact region was quantified using the averaged full-width-at-80%-maximum from the line intensity profiles for TIRF images. Figures 1(b) and 1(c) show measurements of adhesion contact radius vs the bubble resting radius for the different conditions. It is interesting to note the irregular and less consistent contact regions formed for the biotin/avidin condition, possibly owing to heterogeneity in the physisorption of streptavidin to the glass, or streptavidin coating some of the microbubbles.

Kooiman et al.13 and Lajoinie et al.14 reported microbubble adhesion deformation into a spherical cap geometry and quantified this deformation by calculating a contact ratio. Here, the contact ratio is defined as Rcontact/R0, where Rcontact is the radius of the adhesion contact region and R0 is the resting radius of the microbubble. Biotinylated microbubbles resting against a streptavidin surface had a median contact ratio of 0.28, compared to non-biotinylated microbubbles that had a median contact ratio of 0.13. This result confirmed our assumptions that the biotinylated microbubbles were adhered via specific interactions and that non-biotinylated microbubbles were not adhered to the streptavidin surface. Similarly, microbubbles resting against a clean glass surface had a median contact ratio of 0.56, compared to those resting against a pre-coated lipid-on-glass surface that had a median contact ratio of 0.15. Again, this result confirmed nonspecific adhesion and non-adhesion for the two conditions, respectively. Brownian motion of microbubbles in the non-adherent conditions was also observed and supports our findings. A Shapiro–Wilk test for normality showed that the contact ratios for each case were not all normally distributed, so nonparametric statistical testing was conducted using OriginPro software and was used for all the data throughout this study. Using the Mann–Whitney U-test, the contact ratios were significantly greater for the adherent microbubbles than for the non-adherent microbubbles for both specific and nonspecific interactions (p < 0.01).

The resonant behavior of single microbubbles bound to a solid surface were quantified using a photoacoustic technique.11 Briefly, a microbubble was driven into small-amplitude radial oscillations by ultrasound waves generated by an amplitude-modulated continuous wave laser source. This 1550-nm infrared generation laser source illuminated water in the sample (3 μm 1/e diameter and 5.5 mW) approximately 100 μm away from the bubble. Absorption of the laser energy by water caused localized thermal expansion and ultrasound generation through the photoacoustic effect. We previously estimated that the laser source also produces a local temperature rise of approximately 2.5 °C at the microbubble position, and all measurements were conducted at room temperature. A second continuous wave laser source with a wavelength of 532 nm (9.5 μm 1/e diameter and 1.05 mW) was used to illuminate the microbubble and track its oscillations by means of forward scattered light. A photodetector and lock-in amplifier were used to track the light scattering from the microbubble at a single driving frequency. The frequency of the amplitude modulation for the generation laser was swept across a range of at least 2 MHz in 50 kHz steps, producing a resonance curve for a single microbubble. A microscope image of the microbubble was captured prior to measurement using a CCD camera to determine the bubble resting radius.11 

Figure 2(a) shows resonance curves for a biotinylated microbubble and a non-biotinylated microbubble resting against streptavidin-coated glass, both with a resting radius of 2.7 μm. The adherent microbubble had a resonance frequency of 3.4 MHz, which was much higher than the non-adherent microbubble resonance frequency of 1.5 MHz. Figures 2(b) and 2(c) show plots of eigenfrequency vs bubble radius for the different conditions. As with the contact ratio measurements, more variability was seen in the frequency response for the biotin/avidin system than for the lipid/glass system.

FIG. 2.

(a) Resonance curves for two microbubbles with resting radii of 2.7 μm for the biotin/avidin condition. The solid curves are Lorentzian fits to the data. Plots of eigenfrequency vs bubble resting radius for individual microbubbles measured for (b) biotin/avidin and (c) lipid/glass conditions. Solid curves are calculated from the linearized Rayleigh–Plesset model using the median elasticity values, χ, for each condition (see the text for details). Dashed curves are for an uncoated bubble (χ = 0 N/m).

FIG. 2.

(a) Resonance curves for two microbubbles with resting radii of 2.7 μm for the biotin/avidin condition. The solid curves are Lorentzian fits to the data. Plots of eigenfrequency vs bubble resting radius for individual microbubbles measured for (b) biotin/avidin and (c) lipid/glass conditions. Solid curves are calculated from the linearized Rayleigh–Plesset model using the median elasticity values, χ, for each condition (see the text for details). Dashed curves are for an uncoated bubble (χ = 0 N/m).

Close modal

The lipid shell elasticity, χ, was calculated for each measured microbubble response using the linearized form of the modified Rayleigh–Plesset equation,3,15–17 which describes the small-amplitude radial oscillations of a microbubble. The shell elasticity is given as

χ=1.5π2f2ρLR03/12ζ23κP0R0/43κ1σ/2,
(1)

where f is the bubble resonance frequency, ρL is the density of surrounding liquid, R0 is the resting radius of the microbubble, ζ is the damping ratio, κ is the polytropic exponent for the gas core, P0 is the ambient pressure, and σ is the surface tension. The surface tension was assumed to be zero to meet the criterion for bubble stability against dissolution in a saturated medium.18 The damping ratio is given by ζ=Δf/2f, where Δf is the full-width-at-half-maximum at the peak of the resonance curve. Equation (1) includes a correction for the expected (2/3)0.5 decrease in resonance frequency from the acoustic effect of a rigid boundary.19 

The biotinylated microbubbles adherent to the streptavidin surface had a greater median elasticity of χ = 2.71 N/m, compared to the non-adherent microbubbles on the streptavidin surface where χ = 0.57 N/m (p < 0.01). Similarly, microbubbles adherent to the clean glass surface had a greater median elasticity value of χ = 2.84 N/m, compared to non-adherent microbubbles on the lipid surface where χ = 0.57 N/m (p < 0.01). The apparent elasticity of adherent microbubbles agrees with our prior photoacoustic measurements,11,19 while that for non-adherent microbubbles agrees well with a prior high-speed videomicroscopy study for similar shell compositions.20 These median elasticity values were used with Eq. (1) to plot the solid curves seen in Figs. 2(b) and 2(c). The dashed curves represent the frequency response of uncoated bubbles, where χ = 0 N/m. Interestingly, the frequency response of the adherent microbubbles was well-described by the modified Rayleigh–Plesset model even though the model assumes spherical oscillations in an unbounded medium.

Figures 3(a) and 3(b) show the damping ratio vs bubble radius for the same bubbles measured in Figs. 2(b) and 2(c), respectively. The solid lines are the median values. The adherent biotinylated microbubbles had a lower median damping ratio (ζ = 0.14) than the non-adherent microbubbles (ζ = 0.23) in contact with the streptavidin surface (p < 0.01). Similarly, the nonspecific adherent microbubbles had a lower median damping ratio (ζ = 0.10) compared to their non-adherent counterparts (ζ = 0.20) (p < 0.01).

FIG. 3.

Plots of damping ratio vs bubble resting radius for the (a) biotin/avidin and (b) lipid/glass conditions. Solid lines indicate median damping ratio for each condition. Also shown are plots of the maximum radial displacement vs bubble resting radius for the (c) biotin/avidin and (d) lipid/glass conditions. Solid curves show the predicted max radial displacement based on the linear harmonic oscillator model (see the text for details).

FIG. 3.

Plots of damping ratio vs bubble resting radius for the (a) biotin/avidin and (b) lipid/glass conditions. Solid lines indicate median damping ratio for each condition. Also shown are plots of the maximum radial displacement vs bubble resting radius for the (c) biotin/avidin and (d) lipid/glass conditions. Solid curves show the predicted max radial displacement based on the linear harmonic oscillator model (see the text for details).

Close modal

The scatter plots in Figs. 3(c) and 3(d) show the maximum radial displacements vs bubble radius. The maximum radial displacements were calculated from the maximum signal amplitude of the resonance curve using the calibration procedure described previously.11 In the conversion from the measured response to the radial displacement, we assume spherically symmetric oscillations for both adherent and non-adherent bubble populations. The displacement was also found from the measured damping ratio and resonance frequency using the harmonic oscillator model. The maximum radial displacement, ΔRmax, is given by

ΔRmax=Pa12ζ2/8π2ζf2ρLR01ζ2,
(2)

where Pa is the acoustic driving pressure. The parameters ζ and f were set to the corresponding median values determined experimentally. Pa was estimated to be 2.8 Pa for all cases from a best fit to the experimental data [solid curves in Figs. 3(c) and 3(d)]. Adherent microbubbles had lower maximum radial displacements than non-adherent microbubbles for both the specific and nonspecific interactions (p < 0.01).

To summarize, we report on the behavior of individual adherent microbubbles undergoing small-amplitude oscillations. While it is difficult to compare our findings directly with past literature due to differences in experimental parameters, we can compare our results to a similar study by Rooij et al.,10 who investigated the dynamics of individual biotinylated microbubbles with lipid coatings composed primarily of DPPC lipid that were bound to a streptavidin surface. They also reported higher resonance frequencies for adherent DPPC microbubbles compared to non-adherent microbubbles driven by 50 kPa ultrasound although the increase was not statistically significant. However, Rooij et al. also reported an increase in maximum radial displacements for adherent microbubbles compared to non-adherent microbubbles, while we observed a decrease in maximum radial displacements for adherent microbubbles. This disagreement could involve how Rooij et al. monitored the microbubble oscillations optically in one plane-parallel to the contact surface using a high-speed camera while our measurements are sensitive to any dynamic changes in the optical scattering cross sections of the microbubbles. We note that our reported maximum radial displacements could also be influenced by the contact ratio and non-spherical oscillations of microbubbles in the adherent bubble populations. The disagreement with maximum radial displacements could also be due to the differences in amplitudes of oscillation, where our study involved very small amplitudes of oscillation (tens of picometers) and Rooij et al.'s microbubbles were driven into larger amplitudes of oscillation (hundreds of nanometers).

Additionally, Rooij et al. used a polymer substrate (OptiCell) instead of glass. For comparison, we repeated our nonspecific binding experiment using the OptiCell (Fig. 4). TIRF microscopy confirmed adhesion and non-adhesion for the uncoated OptiCell vs lipid-coated OptiCell, respectively. Consistent with our results for glass, microbubbles adherent to the rigid polymer substrate showed an increase in resonance frequency, decrease in damping ratio, and decrease in maximum radial displacement. Changes in shell elasticity, damping ratio, and maximum radial displacement were all statistically significant (p < 0.01). Thus, our results are consistent for both glass and rigid polymer surfaces.

FIG. 4.

Results for microbubbles in contact with an OptiCell polymer surface for the nonspecific adhesion condition. (a) Representative TIRF and epifluorescence microscopy images and plot of contact radius vs resting radius (median contact ratios are 0.58 for adherent and 0.18 for non-adherent). (b) Plots of eigenfrequency vs bubble radius (χ = 2.40 N/m for adherent and χ = 0.55 N/m for non-adherent). The dashed curve is for an uncoated bubble (χ = 0 N/m). (c) Plot of damping ratio vs bubble radius (ζ = 0.15 for adherent and ζ = 0.12 for non-adherent). A nonparametric Mann–Whitney U test confirmed that damping ratios for the adherent microbubbles on OptiCell were statistically greater than the values for nonadherent microbubbles. (d) Plot of the maximum radial displacement vs bubble radius. Scale bars are 2 μm.

FIG. 4.

Results for microbubbles in contact with an OptiCell polymer surface for the nonspecific adhesion condition. (a) Representative TIRF and epifluorescence microscopy images and plot of contact radius vs resting radius (median contact ratios are 0.58 for adherent and 0.18 for non-adherent). (b) Plots of eigenfrequency vs bubble radius (χ = 2.40 N/m for adherent and χ = 0.55 N/m for non-adherent). The dashed curve is for an uncoated bubble (χ = 0 N/m). (c) Plot of damping ratio vs bubble radius (ζ = 0.15 for adherent and ζ = 0.12 for non-adherent). A nonparametric Mann–Whitney U test confirmed that damping ratios for the adherent microbubbles on OptiCell were statistically greater than the values for nonadherent microbubbles. (d) Plot of the maximum radial displacement vs bubble radius. Scale bars are 2 μm.

Close modal

The surfaces used here (glass and OptiCell) are sufficiently rigid compared to the microbubble that their acoustic response is negligible. This was confirmed in our experiments without a microbubble, where no resonance peak was observed. Additionally, the non-adherent microbubble measurements for both substrates (coated-glass or coated-OptiCell) had very similar resonance frequencies and responses, leading us to conclude that the forward light scattering technique was independent of the substrate mechanics.

The mechanism for the observed stiffening is currently unknown. One possible explanation is that the contact line of the adhesion zone imparts an additional tension in the lipid shell that manifests as an increase in the apparent elasticity. For liposomes, adhesion is known to induce a significant force normal to the contact line of the adhesion zone.21 This adhesive tension was shown to deform a soft substrate (hence it was called a “traction force”), pulling up a rim of soft substrate material along the contact line.21 One can imagine a similar situation for microbubbles, where an additional tension in the shell induced by adhesion increases the overall system stiffness and apparent microbubble shell elasticity. Future work combining modeling and experimentation may further elucidate the mechanism.

In conclusion, we verified microbubble adhesion and non-adhesion under specific or nonspecific surface interactions by use of TIRF microscopy. We then measured stable, linear oscillations using a photoacoustic technique to show that adherent microbubbles exhibit an increase in resonance frequency, decrease in damping ratio, and decrease in maximum radial displacement compared to non-adherent microbubbles. Comparison to the linearized Rayleigh–Plesset model with shell viscoelastic terms indicates that adhesion stiffens the system, instead of increasing its inertia or friction. While the apparent shell stiffness is observed to increase, the stiffening effect occurs over the system as a whole. Our working hypothesis is that the adherent microbubble is effectively pinned at the microbubble/substrate contact perimeter, causing a restoring force to be applied across the rest of the microbubble that is in contact with the aqueous medium. In our prior photoacoustic measurements,11,19 it is likely that the microbubbles were adherent as they gave similar shell elasticities. Future studies seeking to characterize the mechanics of free lipid-coated microbubbles should avoid adhesion by pre-incubating the substrate with phospholipid. The observed microbubble stiffening may lead to deeper understanding of the dynamics of microbubbles adherent to different surfaces and shows that surface forces dominate the acoustic image-bubble effect for small-amplitude oscillations.

This work was supported by the National Institutes of Health (No. R01CA195051). J.L. acknowledges support from the CU Boulder Mechanical Engineering Research Innovation Fellowship. D.F.K. acknowledges support from the U.S. Defense Threat Reduction Agency (Award No. HDTRA1-16-1-0045). The authors are grateful to Gazendra Shakya and Hendrik Vos for insightful discussions about the experimental results.

1.
M. A.
Borden
and
K.-H.
Song
,
Adv. Colloid Interface Sci.
262
,
39
(
2018
).
2.
A. A.
Doinikov
,
J. F.
Haac
, and
P. A.
Dayton
,
Ultrasonics
49
,
263
(
2009
).
3.
N.
de Jong
,
L.
Hoff
,
T.
Skotland
, and
N.
Bom
,
Ultrasonics
30
,
95
(
1992
).
4.
B.
Helfield
,
Ultrasound Med. Biol.
45
,
282
(
2019
).
5.
B.
Dollet
,
P.
Marmottant
, and
V.
Garbin
,
Annu. Rev. Fluid Mech.
51
,
331
355
(
2019
).
6.
S.
Zhao
,
K. W.
Ferrara
, and
P. A.
Dayton
,
Appl. Phys. Lett.
87
,
134103
(
2005
).
7.
M.
Overvelde
,
V.
Garbin
,
B.
Dollet
,
N.
de Jong
,
D.
Lohse
, and
M.
Versluis
,
Ultrasound Med. Biol.
37
,
1500
(
2011
).
8.
B. L.
Helfield
,
E.
Cherin
,
F. S.
Foster
, and
D. E.
Goertz
,
Ultrasound Med. Biol.
39
,
345
(
2013
).
9.
J.
Casey
,
C.
Sennoga
,
H.
Mulvana
,
J. V.
Hajnal
,
M.-X.
Tang
, and
R. J.
Eckersley
,
Ultrasound Med. Biol.
39
,
903
(
2013
).
10.
T.
van Rooij
,
I.
Beekers
,
K. R.
Lattwein
,
A. F. W.
van der Steen
,
N.
de Jong
, and
K.
Kooiman
,
IEEE Trans. Ultrason., Ferroelectr., Freq. Control
64
,
785
(
2017
).
11.
J. S.
Lum
,
D. M.
Stobbe
,
M. A.
Borden
, and
T. W.
Murray
,
Appl. Phys. Lett.
112
,
111905
(
2018
).
12.
J. A.
Feshitan
,
C. C.
Chen
,
J. J.
Kwan
, and
M. A.
Borden
,
J. Colloid Interface Sci.
329
,
316
(
2009
).
13.
K.
Kooiman
,
T. J. A.
Kokhuis
,
T.
van Rooij
,
I.
Skachkov
,
A.
Nigg
,
J. G.
Bosch
,
A. F. W.
van der Steen
,
W. A.
van Cappellen
, and
N.
de Jong
,
Eur. J. Lipid Sci. Technol.
116
,
1217
(
2014
).
14.
G.
Lajoinie
,
Y.
Luan
,
E.
Gelderblom
,
B.
Dollet
,
F.
Mastik
,
H.
Dewitte
,
I.
Lentacker
,
N.
de Jong
, and
M.
Versluis
,
Commun. Phys.
1
,
22
(
2018
).
15.
L.
Hoff
,
P. C.
Sontum
, and
J. M.
Hovem
,
J. Acoust. Soc. Am.
107
,
2272
(
2000
).
16.
K.
Sarkar
,
W. T.
Shi
,
D.
Chatterjee
, and
F.
Forsberg
,
J. Acoust. Soc. Am.
118
,
539
(
2005
).
17.
P.
Marmottant
,
S.
van der Meer
,
M.
Emmer
,
M.
Versluis
,
N.
de Jong
,
S.
Hilgenfeldt
, and
D.
Lohse
,
J. Acoust. Soc. Am.
118
,
3499
(
2005
).
18.
P. S.
Epstein
and
M. S.
Plesset
,
J. Chem. Phys.
18
,
1505
(
1950
).
19.
J. S.
Lum
,
J. D.
Dove
,
T. W.
Murray
, and
M. A.
Borden
,
Langmuir
32
,
9410
(
2016
).
20.
S. M.
van der Meer
,
B.
Dollet
,
M. M.
Voormolen
,
C. T.
Chin
,
A.
Bouakaz
,
N.
de Jong
,
M.
Versluis
, and
D.
Lohse
,
J. Acoust. Soc. Am.
121
,
648
(
2007
).
21.
M. P.
Murrell
,
R.
Voituriez
,
J.-F.
Joanny
,
P.
Nassoy
,
C.
Sykes
, and
M. L.
Gardel
,
Nat. Phys.
10
,
163
(
2014
).