# Coated microbubbles swim via shell buckling

Georges Chabouh, Marcel Mokbel, Benjamin van Elburg, Michel Versluis,
Tim Segers, Sebastian Aland, Catherine Quilliet, Gwennou Coupier

► To cite this version:

Georges Chabouh, Marcel Mokbel, Benjamin van Elburg, Michel Versluis, Tim Segers, et al. Coated microbubbles swim via shell buckling. Communications Engineering, 2023, 2 (1), pp.63. 10.1038/s44172-023-00113-z. hal-04200310

HAL Id: hal-04200310

https://hal.science/hal-04200310v1

Submitted on 8 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Coated microbubbles swim via shell buckling

Georges Chabouh ID <fnref n="1" />✉, Marcel Mokbel², Benjamin van Elburg³, Michel Versluis ID ³, Tim Segers ID ⁴,
Sebastian Aland ID ², Catherine Quilliet<fnref n="1" /> & Gwennou Coupier ID <fnref n="1" />✉

Engineered microswimmers show great promise in various biomedical applications. However, their application is hindered by the slow mobility, limited maneuverability and poor bio-compatibility. Lipid coated microbubbles have high compressibility and are already approved for clinical use as diagnostic ultrasound contrast agents. Here we experimentally investigate the swimming motion of these microbubbles under external cyclic overpressure. A net displacement was generated via reproducible and non-destructive cycles of deflation and re-inflation of the microbubble. We also propose a numerical model which allows a maximum swimming speed on the order of meters per second, which falls in the range of blood flow velocity in large vessels. Unlike the acoustic radiation force technique, where the displacement is always directed along the axis of ultrasound propagation, here, the direction of propulsion is controlled in the shell reference frame. This provides a solution toward controlled steering for ultrasound molecular imaging and drug delivery.

<footnote n="1">¹CNRS/Université Grenoble-Alpes, LIPhy UMR 5588, Grenoble F-38401, France. ²Technische Universität Bergakademie Freiberg, Akademiestrasse, 609599 Freiberg, Germany. ³Physics of Fluids Group, Technical Medical (TechMed) Center and MESA+ Institute for Nanotechnology, University of Twente, 7500 AE Enschede, The Netherlands. ⁴BIOS/Lab-on-a-Chip Group, Max Planck Center Twente for Complex Fluid Dynamics, MESA+ Institute for Nanotechnology, University of Twente, Enschede, The Netherlands. ✉email: georges.chabouh@sorbonne-universite.fr; gwennou.coupier@univ-grenoble-alpes.fr</footnote>
The creation of a programmable machine of several microns, capable of carrying out automatically a complex series of actions in our body, was the famous scientific challenge given by Richard Feynman back in 1959¹. It was achievable only in science fiction movies (Fantastic Voyage 1966 & Inner Space 1987) until the rapid development of microrobots research in the past two decades. Soft microrobots in particular are receiving increasing attention. Their softness, elasticity and adaptivity in their configuration is opening ways for a better interaction with the complex vascularity of the human body, compared to their rigid counterparts². Softness allows for an easy shape changes, similar to what is observed in biological swimmers performing either tail or flagella beating or full deformation cycles of their body, leading to what can be called ameboid swimming³⁻⁹, by analogy with the crawling-type ameboid motion of cells adhering on a substrate. However, bio-inspired artificial microswimmers share with their natural counterparts an essential flaw: they are slow. Micron-sized swimmers performing shape deformations hardly overcome the 100 µm/s limit¹⁰⁻¹³. Though inspiration emerging from nature has led to elegant motility strategies at microscale, this approach contains by itself intrinsic limitations set by the fluid mechanics at this scale, compromising faster displacements. This hinders their use in complex, crowded and quickly flowing media such as blood. More generally, finding the compromise between biocompatibility and swimming performance is an ongoing quest¹⁴.

When a microrobot or a microorganism performs a shape deformation, it generates a fluid flow around it. The flow is characterized by the Reynolds number $Re = UL/v_\beta$ with U the typical deformation velocity, L the characteristic length of the microswimmer, and $v_\beta$ the kinematic viscosity of the surrounding fluid. In the low Reynolds number regime ("Stokes flow"), thrust (but also drag) is mediated by viscous friction with the fluid. A net displacement can only be produced if the cyclic sequence of deformations is not reversible, a principle known as the “scallop theorem”¹⁵. Deflation and re-inflation of a spherical shell usually involves buckling, providing the much desired feature¹⁶ to the deformation cycle¹⁷⁻²¹ (Fig. 1a), i.e., an hysteretic path. Furthermore, since buckling as an instability may correspond to a high Reynolds episode, it provides an additional opportunity for self propulsion, mediated this time by inertia. Swimming under such amphibolic efficiency was reported by Djellouli et al., who studied the displacement of air-filled macroscopic shells made of a rubber-like material²² on a large range of Reynolds numbers.

Soft elastic shells also exist at the microscale. Thin shells (1–10 nm) encapsulating tiny gas bubbles (1–5 µm in radius) to avoid gas dissolution in the blood plasma are used daily in clinical

practice as ultrasound contrast agents (UCAs). UCAs have shells with different visco-elastic properties depending on their material (protein, lipids or polymers)²³⁻²⁵. Coated microbubbles are highly compressible hence they present high echogenicity which is mandatory for ultrasound imaging.

The use of UCAs in the medical domain goes beyond imaging; in fact, UCAs can be used as drug carriers²⁶,²⁷ whose cargo is released on-demand at the target site, thereby limiting adverse systemic effects. The release of drugs is triggered by an ultrasonic pulse localized in space and time above a threshold amplitude high enough to break the shell. Despite the huge potential of bubble-induced drug delivery systems, the circulation of the drug carriers in the bloodstream is still not controlled, which hinders their translation to clinical practice. The absence of specificity and sensitivity with the passively circulating microbubbles jeopardizes the accurate delivery of drugs.

Here, we activate UCAs by cyclic low amplitude variations of the external pressure, which is achieved in the experiments by pressurizing a closed chamber. We recently demonstrated that the anisotropy of lipid shells induces a deflation-inflation pattern, distinct from the one observed in homogeneous shell materials²¹. We examine the resulting displacements due to sequences of these non-symmetric deformations, as depicted in Fig. 1.

## Results

**Slow activation.** We first applied sinusoidal pressure cycles which varied from atmospheric pressure $P_{atm} = 101.3$ kPa to a maximum value $P_{max}$ between 110 and 160 kPa at two different frequencies $f=1$ Hz and $f=2$ Hz (Fig. 1). We systematically observed that when the maximal pressure $P_{max}$ was set to be slightly above a given buckling pressure $P_b$, the microbubbles would loose their spherical symmetry at threshold $P_b$ and collapse into a bowl-like shape. This deformation was reversible and could be repeated several times through slow pressure variations between $P_{atm}$ and $P_{max}$ with no apparent damage as can be seen on Fig. 1b, that shows the shape evolution within one pressure cycle.

**Ballasted microswimmer.** The fortuitous self-assembly of a three-body swimmer composed of a Sonovue microbubble ballasted by two non-buckling beads (insets in Fig. 2a, b) allows us to introduce the swimming mechanism. It was aided by the action of buoyancy that allowed for a long observation of the microswimmer over 25 cycles. The ballasted microswimmer was activated at 2 Hz with a sinusoidal pressure of amplitude 20 kPa. While the dominant motion of the swimmer, which is denser than the surrounding water, is set by its interaction with the gravity field, its trajectory shows oscillations with a frequency of 2 Hz (Fig. 2a).

A zoom on a single period reported in Fig. 2b shows that four steps can be identified, where motion is related to deformation. Upon deflation, the ballasted microswimmer undergoes a buckling event and swims up against gravity for a duration $t_1$, which results in a quasi stationary height.

After buckling, the microswimmer sediments for a duration $t_2$ while keeping a buckled shape. After that, the reinflation phase is marked by a reverse buckling event — often improperly called debuckling phase—for a duration $t_3$ during which the microswimmer sinks faster. It is followed by a slow sedimentation phase in the purely spherical shape, of duration $t_4$.

The quasi-stagnation during the buckling phase and the enhanced downward motion during the debuckling phase are strong indicators of a swimming motion that points in the center-to-buckling spot direction upon deflation and in the opposite direction upon inflation. These motions are all the more

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Deflation and re-inflation cycle of a coated microbubble.</figcaption></figure>

a Schematic of the time-asymmetric deformation cycle of an hollow shell upon an increase and decrease of external pressure P. b Snapshots of a coated microbubble, where a buckling spot appears on the top, grows then relaxes. The scale bar represents 4 µm.
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Trajectory analysis of the ballasted microswimmer. **a** Part of the trajectory of the ballasted microswimmer with and without excitation (see Supplementary Video S2). The sedimentation trajectory without excitation was used to set the gravity direction **y**. **b** Zoom on the evolution of the **y** position over one cycle (period *T* = 0.5 s), together with the evolution of the aspect ratio *λ* of the shell that quantifies its deformation (see Methods section). Insets show two selected configurations, where the central shell (a UCA of radius 4.5 µm) is buckled and not buckled. *t*₁ is the duration of the buckling phase, *t*₂ of the post-buckling phase, *t*₃ of the debuckling phase and *t*₄ of the spherical phase. Vertical velocity *v*<sub>y</sub> (not corrected by gravity) of the ballasted microswimmer for one period *T* = 0.5 s. **c** stack of successive 25 cycles (different blue lines), **d** debuckling phase with time shift correction overlaying all 25 profiles (different blue lines). The scale bar represents 4 µm. The letter **g** denotes the gravity constant.</figcaption></figure>

significant as the variations in the buoyancy during these two phases should enhance downward motion in the first case and counter-act this downward motion in the second case. We estimated the contribution of gravity (weight+buoyancy) on this microswimmer of varying volume but constant mass (see Supplementary information: Note 1) and determined that the mean net displacement in the y direction *d*<sub>net</sub> is positive and bounded by: 0.16 µm/cycle *<d*<sub>net</sub>*<0.94* µm/cycle.

As far as the dynamics is concerned, both buckling and debuckling phases are marked by important oscillations of frequency slightly larger than 20 Hz (Fig. 2c). For all the 25 cycles, these oscillations start at exactly the same time for the post-buckling ones, unlike those observed upon reinflation. However, when each of the 25 profiles is individually shifted, the first peaks of these profiles overlap perfectly (Fig. 2d). The existence of a reproducible buckling threshold and of a well defined oscillation frequency illustrate the robustness of the deformation process and opens the way toward a better steering control by appropriate activation by sound waves.

**Free contrast agents.** Isolated contrast agents that float by buoyancy do not present a constant orientation $\theta_b$ (Fig. 3a) of their buckling spot, in contrast with the ballasted microswimmer; in addition to the variability between each observed microbubble, this orientation evolves in time as the microbubble rises. We showed in ref. 21 that microbubbles have a single buckling spot, which may be due to the presence of a defect or be dictated by the locus of the first buckling event. The orientation of this buckling spot may differ at each cycle because of rotational Brownian motion of the microbubbles.

Buoyancy driven rise velocity is of order 20 µm/s and varies with the volume of the microbubble: it makes it difficult to measure the component of the swimming motion in the direction

of gravity. However, the contribution $v_x$ of the swimming motion in the transverse direction is not affected by buoyancy effects. If the buckling spot has an orientation $\theta_b$ compared to x axis, $v_x/\cos\theta_b$ yields, for each deformation cycle, a measure of the swimming motion in the center-to-buckling-spot direction.

The gravity direction may be potentially biased by micro-streaming due to the presence of trapped air bubbles or material deformation. Therefore, in order to detect the gravity direction with more precision, we mixed the microbubble suspension with hollow glass beads (HGMS, Cospheric, USA) of 5 µm average diameter with a density of 0.1–0.7 g/cm³. These beads did not buckle when a pressure was applied and could be used as reference particles to detect the main sedimentation direction $\mathbf{y}$.

*Correlation between deformation and motion.* The trajectory of a rising microbubble is represented in Fig. 3a. A marked deviation from the $\mathbf{y}$ direction is visible. In this selected case where the buckling angle is constant along the cycles, $\theta_b = 150^\circ$, it is interesting to directly quantify the correlation between deformation and motion. In Fig. 3b, the aspect ratio $\lambda$ of the bubble, averaged over 25 cycles, is plotted as a function of time, together with the velocity $v_x$ in the direction orthogonal to the main rise direction $\mathbf{y}$. One can clearly see that upon buckling (drop in $\lambda$), the velocity $v_x$ along x peaks to a negative value with a postbuckling oscillation. Inversely, upon re-inflation, where $\lambda$ goes back to its initial value, $v_x$ peaks to a positive value. However, the absolute value of $v_x$ is smaller than upon deflation. Hence over a full cycle of deflation and re-inflation the net displacement is non-zero.

*Swimming direction and intensity.* In the above example, where $\cos(\theta_b)<0$, a net motion in the direction of the buckling pole is observed (mean $v_x$ is negative). More generally, the buckling angle $\theta_b$ can change from cycle to cycle and the information of interest is the velocity $v_b$ along the buckling axis (of angle $\theta_b$). We deduce this
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Swimming motion estimation. **a** Example of microbubble track obtained by minimum intensity projection (see Supplementary Video S3), red arrow presents the direction of the corrected gravity through the reference bead, the scale bar represents 4 µm, **b** averaged aspect ratio λ (left side) and swimming velocity $v_x$ (right side) on 25 cycles, for the same microbubble (dashed lines present the error ± standard deviation). Net displacements per cycle for SonoVue microbubbles in function of the radius at (c), $f=1$ Hz, (d), $f=2$ Hz, and (e) for in-house coated microbubbles at $f=1$ Hz. The data shows the average net displacement on all the cycles (at least 7 cycles per shell) with the error bar being the standard deviation. The orange points indicate the net displacement per cycle (based on $v_x$) of the reference beads. The letter **g** denotes the gravity constant.</figcaption></figure>

velocity from the measurement of its component along the x axis: $v_b = v_x / \cos(\theta_b)$, which is calculated for each cycle. This velocity is then integrated over one cycle and averaged over all cycles, to obtain the net displacement. Figure 3c, d shows the net displacement per cycle for SonoVue microbubbles for an external driving frequency of, respectively, 1 Hz and 2 Hz. Figure 3e shows the net displacement per cycle for the microbubbles made in-house, which are stiffer than Sonovue microbubbles. The displacements are compared to that of the non-buckling reference beads. A prominent feature in this plot is the non-zero net displacement of the microbubbles. If we consider a single point of the data, it has a high standard deviation that indicates differences in the behavior from cycle to cycle. This can result from an inaccuracy of the measurement of the buckling direction, and from the loss of the contribution in the out-of-plane direction.

In the following, we refer to a “forward” propulsion when the net displacement is positive (i.e., in the center-to-buckling-spot direction) and to a “backward” propulsion when the net displacement with respect to the buckling spot orientation is negative. For SonoVue microbubbles at $f=1$ Hz (Fig. 3c), we notice that every coated microbubble is executing a backward motion. While for the SonoVue at a frequency of 2 Hz and the in-house microbubbles, we actually have both scenarios, forward and backward motion. Different swimming directions may be due to heterogeneities in shell properties. In fact, coated microbubbles of the same size are known to have different acoustic response²⁸. The overlaid boxplot in Fig. 4a summarizes all the swimming experiments on SonoVue microbubbles for two different frequencies $f=1$ Hz and $f=2$ Hz and for the in-house microbubbles at one frequency $f=1$ Hz. By comparison with the reference beads, a statistical significance in the absolute net displacement of microbubbles is shown through 2 way ANOVA test. In-house microbubbles show higher net displacement compared to SonoVue, which can be understood from their higher elastic in-plane compression modulus ($\chi=2$ N/m for the in-house microbubbles versus ~0.5 N/m for SonoVue): the volume set in motion upon buckling increases with the elastic modulus¹⁸.

## Towards tunable high velocities

While a typical motion of the micron per cycle is clearly demonstrated here, the intrinsic variability in the behavior between each coated microbubble prevents a systematic prediction for the motion of a given, unknown, UCA. Nevertheless, we can address the question of the average behavior of UCAs by means of numerical simulations. This can be useful for predicting how a colony of similar microswimmers will behave, or how a given UCA will behave after ‘calibrating’ the simulations, having identified the basic characteristics of its behavior.

We consider zero-thickness elastic microbubbles of initial radius $R_0$ (stress-free configuration) characterized by an in-plane compression modulus $\chi$, a Poisson ratio $\nu$ and a bending modulus $\kappa$. The characteristic distance for elasticity $d_{eff} = \sqrt{6(1+\nu)\kappa/\chi}$, equal to the thickness $d$ for a thin shell of an isotropic material, drives the deformation of shells in the buckled configuration²¹.

In a recent study on the buckling threshold and final shape of UCAs, we estimated that $d_{eff}/R_0 = 10\%$²¹. The characteristic distance $d_{eff}$ would then be much larger than the typical thickness $d \approx 5$ nm of the shell. This result, which is in line with another study on gel-phase lipid vesicles²⁹, is also confirmed here: by considering the thick shell, simulations show a good agreement with the observed experimental displacement. The microbubbles are immersed in a Newtonian fluid for which the Navier-Stokes equation is solved numerically. They are filled with a gas at pressure $P$ which is assumed to be instantaneously set by the shell volume according to an adiabatic process, a reasonable hypothesis considering the high velocities encountered in this problem: $PV^\kappa = P_0 V_0^\kappa$, where $V$ is the microbubble volume, $V_0$ its initial value, $P_0$ the initial pressure, and $\kappa = 1.4$. In the simulations, forward motion is always observed at slow frequency, which points to the necessity to explore the small scale to understand the backward motion observed in some experiments. Figure 4b shows that, upon slow cyclic pressurization, swimming intensity predicted by simulations is in line with experimental results, providing an effective thickness over radius ratio of 10% is chosen. This agreement leads us to confidently examine what would
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Buckling induced swimming motion through experiments and modeling. **a** The absolute net displacement of SonoVue microbubbles (two different frequencies), in-house microbubbles and control (glass beads). Each symbol represents net displacement averaged on at least eight cycles. The overlaid boxplots show the median, interquartile range, mean, and the minimum to maximum values. Statistical significance of microbubbles compared to the beads. It is indicated with * (*p* < 0.01) and **** (*p* < 0.0001) using a 2 way ANOVA test on GraphPad Prism 9 software. **b** The minimum and maximum displacement per cycle from the numerical simulations at low frequency with the two assumptions *d/R* of 1% and 10%, and initial radius in the range 1 to 5 microns. **c** Displacement per cycle in the case (*d/R* = 10%, *χ* = 2 N/m, *R*₀ = 5 µm) for different activation frequencies and pressure wave of amplitude 800 mbar. **d** The estimated swimming velocity in function of the varying excitation frequency for three different compression moduli *χ* = 0.5, 2 and 4 N/m. (obtained through a linear fit of displacement (*R*² > 0.95), see also Fig. 5 in Supplementary information: Note 2).</figcaption></figure>

be the effect of a faster activation of the coated microbubble. Figure 4c shows the displacement for sub MHz frequency activation of a shell with elastic modulus $\chi = 2$ N/m (see Supplementary Video S1). For this shell, the swimming velocity increases up to 1 m/s at $f = 400$ kHz (Fig. 4d). The quasi linear increase of the velocity with the frequency indicates a regime of independent successive buckling events. At higher frequencies, an interference with the own dynamics of the shell (whose post-buckling resonance intrinsic resonance frequency is close to 400 kHz³⁰) partly inhibits the propulsion.

By variation of the shell elasticity, and for similar amplitude of the acoustic wave, the optimal swimming velocity can be reached at different activation frequencies as seen on Fig. 4d. As already explained, softer shells swim more slowly but interestingly, in the shell-wave interaction regime reached at high frequencies, they can also swim backward (Fig. 5a in Supplementary information: Note 2). More rigid shells illustrate here the subtlety of the interactions with the wave: while they swim faster at low frequency (10–50 kHz) (Fig. 5b), they are not able to buckle at higher frequencies (Fig. 5c), since they are very close to their buckling threshold (see Fig. 5d). Indeed, as pointed out by Pelekasis et al., the time needed for buckling depends both on the activating frequency and on its amplitude³¹.

Thus, when considering an assembly of coated microbubbles with different elastic moduli, that are activated by a given pressure amplitude or a fixed input voltage of a medical echograph, an adequate choice of the pressure amplitude with respect to the buckling threshold of each shell appears to tune the relative

contribution of each shell. In other words: one can simply control the steering of the cargo by varying the frequency.

The proposed swimming mechanism allows to reach velocities up to 1 m/s for a micron-size object performing shape deformations. This is above the performances of partially shelled microbubbles, rapidly propelled thanks to acoustic streaming³²,³³, but these latter significantly influence each other³⁴, which impedes the realization of a multidirectional swimmer. Meanwhile, the possibility of “activating” only a set of coated microbubbles enables fine steering as long as the orientation of the different actors can be monitored in real time. The latter point is made possible thanks to a recent breakthrough in ultrasound imaging, due to in-vivo combination of ultrasound and microbubbles³⁵,³⁶ in an ultrafast and super-resolved technique known as Ultrasound Localization Microscopy (ULM). ULM allows localization and tracking of microbubbles within 10 microns in resolution, and accurate 3D measurements of their velocities (~mm/s) deep in organs³⁷. Thus, a dual-frequency ultrasound probe can achieve simultaneously propulsion and real time imaging³⁸,³⁹. Additionally, the type of powering we studied in this article preserves the displacement's degrees of freedom: (i) the swimmer is activated by a scalar parameter, namely pressure, which precludes any reorientation effect that may affect, e.g., swimmers activated by magnetic fields (ii) the swimming direction is completely independent of the direction of the activation wave, by contrast with microbubble control mechanisms based on the primary radiation force (“Bjerkness forces”)⁴⁰. In a complex 3D environment such as the human microvasculature where important
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Displacement from numerical simulation. a $\chi$ = 0.5 N/m, b $\chi$ = 2 N/m and c $\chi$ = 4 N/m. The amplitude of pressure is set to 800 mbar and the frequency is varied. d Time-evolution of the aspect ratio $\lambda$ of the shell for the latter case, for some selected frequencies.</figcaption></figure>

exchanges with the vessel walls take place, having a swimmer whose propelling direction is independent of the actuator features is a key point. Last but not least, a coated bubble is a robust object, biocompatible and already approved as clinical agent. Thanks to its complete coating, it has proven its ability to survive to the constraints met in the blood stream, thus making it a complete and trustful candidate for biomedical application.

UCAs were already excellent candidates for cargos enabling drug release through shell destruction and blood brain barrier opening⁴¹. We proved here that the same UCAs can also act as motors. Fast propulsion, detectability, and release: these are the prerequisites for a controlled microrobot dedicated to medical applications. We anticipate that these three functions can be accomplished by a compound of multiple UCAs, one as cargo and others as motors, that would all simply be activated by a single medical echograph.

## Methods

**Numerical simulations.** In the simulations, we assume that the inner pressure is initially of the order of the atmospheric pressure $P_{atm} = 101$ kPa. We have examined coated bubbles with a radius between 1 and 5 µm and compression moduli $\chi$ in the range 0.5–2 N/m. The simulations are based on the numerical method presented in ref. ⁴², applied to the buckling of microswimmers, presented in detail and validated in ref. ³⁰, where postbuckling oscillations were studied.

**Microbubble preparation.** SonoVue® microbubbles (Bracco Spa, Milan, Italy) consist of a phospholipid shell encapsulating a sulfur hexafluoride gas core. Before the experiment, the shells are reconstituted through a mixture of the lyophilisate with 5 ml physiological saline solution, to form a suspension that contains ~2–5 × 10⁸ microbubbles per milliliter with radius ranging from 2 to 7 µm.

We also used lipid microbubbles encapsulating $C_4F_{10}$ with a more narrow size distribution that the one we produced in a flow-focusing device. The fabrication process of these DSPC/DPPE-

PEG5000 microbubbles is summarized in ref. ²¹, following protocols detailed in refs. ⁴³,⁴⁴.

**Experimental methods.** As in ref. ²¹, the microbubble suspension was gently poured into degased water and placed in a Falcon microfluidic reservoir of 15 ml (Fisher Scientific, USA) connected to an Elveflow® pressure controller (Elvesys®, France) and to a flow-through cuvette (Aireka Scientific® Co., Ltd) using PTFE tubings.

The chamber made from quartz with a square cross section (12.5 mm × 12.5 mm) was placed under an inverted microscope (Olympus®, model IX70) which was rotated 90° through three stabilizing aluminum legs. This configuration provided an observation axis $z$ perpendicular to the gravity axis $y$. After injection in the chamber, the microbubbles rose due to buoyancy. The other end of the observation chamber was connected 1) to a valve left open to allow injection of the UCAs into the chamber through a gentle increase of the pressure in the reservoir (of order 30 mPa above atmospheric pressure), then closed for pressurization of the chamber, and 2) to a pressure sensor (MPS1, Elvesys®, France).

In order to correlate the evolution of the shell shape with the driving pressure $P_{max}$, both pressure sensor and fast camera (Miro 310, Vision Research) were triggered through the Elveflow interface. Videos were taken at a rate of 100 frame/s. An automated stage (MS-2000, ASI, USA) was used to select a microshell prior to the recording.

**Image analysis.** To track the UCAs and record their shape at each time step, we developed a tracking algorithm with sub-micron resolution²¹. The intensity profile of each image was approximated by a two-dimensional elliptic Gaussian distribution characterized by a center ($x_0, y_0$), main axis $\sigma'_x$ and $\sigma'_y$ and orientation angle $\theta$. Data were processed so that $\sigma'_x$ would always denote the minor axis and $\theta$ its angle to the $x$ axis. The angle $\theta$ was also chosen so that when the aspect ratio $\lambda = \sigma'_x/\sigma'_y$ is minimal, the corresponding angle $\theta$ indicates the buckling direction, referred to as $\theta_b$ hereafter (i.e., a visual inspection can decide whether $\theta$ or $\theta + \pi$ should be used).
## Data availability

The data that support the plots within this paper and other findings of this study are available from the author upon reasonable request.

## Code availability

The code that support the plots within this paper and other findings of this study are available from the author upon reasonable request.

Received: 14 March 2023; Accepted: 24 August 2023

Published online: 07 September 2023

## References

1. Feynman, R. P. There's plenty of room at the bottom: an invitation to enter a new field of physics. *Miniaturization, Reinhold* (1961).

2. Degen, P. Self-propelling capsules as artificial microswimmers. *Curr. Opinion Colloid Interface Sci.* **19**, 611–619 (2014).

3. Barry, N. P. & Bretscher, M. S. Dictyostelium amoebae and neutrophils can swim. *Proc Natl Acad Sci* **107**, 11376–11380 (2010).

4. Farutin, A. et al. Amoeboid swimming: a generic self-propulsion of cells in fluids by means of membrane deformations. *Phys. Rev. Lett.* **111**, 228102 (2013).

5. Wu, H. et al. Amoeboid swimming in a channel. *Soft Matter* **12**, 7470–7484 (2016).

6. Heintzelman, M. B. Cellular and molecular mechanics of gliding locomotion in eukaryotes. *Int. Rev. Cytol.* **251**, 79–129 (2006).

7. Othmer, H. G. Eukaryotic cell dynamics from crawlers to swimmers. *WIREs Comput. Mol. Sci.* **9**, e1376 (2019).

8. Noselli, G., Beran, A., Arroyo, M. & DeSimone, A. Swimming euglena respond to confinement with a behavioural change enabling effective crawling. *Nat. Phys.* **15**, 496–502 (2019).

9. Aoun, L. et al. Amoeboid swimming is propelled by molecular paddling in lymphocytes. *Biophys. J.* **119**, 1157–1177 (2020).

10. Kim, S. et al. Fabrication and manipulation of ciliary microrobots with non-reciprocal magnetic actuation. *Sci. Rep.* **6**, 30713 (2016).

11. Striggow, F. et al. Sperm-driven micromotors moving in oviduct fluid and viscoelastic media. *Small* **16**, 2000213 (2020).

12. Llacer-Wintle, J. et al. Biodegradable small-scale swimmers for biomedical applications. *Adv. Mater.* **33**, 2102049 (2021).

13. Tan, L., Wang, Z., Chen, Z., Shi, X. & Cheang, U. K. Improving swimming performance of photolithography-based microswimmers using curvature structures. *Micromachines* **13** (2022).

14. Xu, Y., Bian, Q., Wang, R. & Gao, J. Micro/nanorobots for precise drug delivery via targeted transport and triggered release: a review. *Int. J. Pharm.* **616**, 121551 (2022).

15. Purcell, E. M. Life at low reynolds number. *Am. J. Phys.* **45**, 3–11 (1977).

16. Holmes, D. P. Elasticity and stability of shape-shifting structures. *Curr. Opin. Colloid Interface Sci.* **40**, 118–137 (2019).

17. Knoche, S. & Kierfeld, J. Buckling of spherical capsules. *Phys. Rev. E* **84**, 046608 (2011).

18. Quilliet, C. Numerical deflation of beach balls with various poisson's ratios: from sphere to bowl's shape. *Eur. Phys. J. E* **35**, 48 (2012).

19. Knoche, S. & Kierfeld, J. Secondary polygonal instability of buckled spherical shells. *Europhys. Lett.* **106**, 24004 (2014).

20. Coupier, G., Djellouli, A. & Quilliet, C. Let's deflate that beach ball. *Eur. Phys. J. E* **42**, 1–10 (2019).

21. Chabouh, G. et al. Buckling of lipidic ultrasound contrast agents under quasi-static load. *Philos. Trans. R. Soc. A* **381**, 20220025 (2023).

22. Djellouli, A., Marmottant, P., Djeridi, H., Quilliet, C. & Coupier, G. Buckling instability causes inertial thrust for spherical swimmers at all scales. *Phys. Rev. Lett.* **119**, 224501 (2017).

23. Helfield, B. A review of phospholipid encapsulated ultrasound contrast agent microbubble physics. *Ultrasound Med. Biol.* **45**, 282–300 (2019).

24. Versluis, M., Stride, E., Lajoinie, G., Dollet, B. & Segers, T. Ultrasound contrast agent modeling: a review. *Ultrasound Med. Biol.* **46**, 2117–2144 (2020).

25. Chabouh, G., Dollet, B., Quilliet, C. & Coupier, G. Spherical oscillations of encapsulated microbubbles: effect of shell compressibility and anisotropy. *J. Acoust. Soc. Am.* **149**, 1240–1257 (2021).

26. Couture, O., Foley, J., Kassell, N. F., Larrat, B. & Aubry, J.-F. Review of ultrasound mediated drug delivery for cancer treatment: updates from pre-clinical studies. *Transl. Cancer Res.* **3**, 494–511 (2014).

27. Kooiman, K., Vos, H. J., Versluis, M. & de Jong, N. Acoustic behavior of microbubbles and implications for drug delivery. *Adv. Drug Deliv. Rev.* **72**, 28–48 (2014).

28. Segers, T. & Versluis, M. Acoustic bubble sorting for ultrasound contrast agent enrichment. *Lab Chip* **14**, 1705–1714 (2014).

29. Quemeneur, F., Quilliet, C., Faivre, M., Viallat, A. & Pepin-Donat, B. Gel phase vesicles buckle into specific shapes. *Phys. Rev. Lett.* **108**, 108303 (2012).

30. Mokbel, M., Djellouli, A., Quilliet, C., Aland, S. & Coupier, G. Post-buckling dynamics of spherical shells. *Proc. Roy. Soc. A* **477**, 20210378 (2021).

31. Pelekasis, N., Vlachomitrou, M. & Lytra, A. Compression-only behavior: effect of prestress and shell rheology on bifurcation diagrams and parametric stability of coated microbubbles in an unbounded flow. *Phys. Rev. Fluids* **7**, 113601 (2022).

32. Bertin, N. et al. Propulsion of bubble-based acoustic microswimmers. *Phys. Rev. Appl.* **4**, 064012 (2015).

33. Ren, L. et al. 3d steerable, acoustically powered microswimmers for single-particle manipulation. *Sci. Adv.* **5**, eaax3084 (2019).

34. Ahmed, D. et al. Selectively manipulable acoustic-powered microswimmers. *Sci. Rep.* **5**, 9744 (2015).

35. Errico, C. et al. Ultrafast ultrasound localization microscopy for deep super-resolution vascular imaging. *Nature* **527**, 499–502 (2015).

36. Couture, O., Hingot, V., Heiles, B., Muleki-Seya, P. & Tanter, M. Ultrasound localization microscopy and super-resolution: a state of the art. *IEEE Trans. Ultrason. Ferroelectr. Freq. Control* **65**, 1304–1320 (2018).

37. Chavignon, A. et al. 3d transcranial ultrasound localization microscopy in the rat brain with a multiplexed matrix probe. *IEEE Trans. Biomed. Eng.* **69**, 2132–2142 (2021).

38. Martin, K. H. et al. Dual-frequency piezoelectric transducers for contrast enhanced ultrasound imaging. *Sensors* **14**, 20825–20842 (2014).

39. Myhre, O. F., Johansen, T. F. & Johan Angelsen, B. A. Analysis of acoustic impedance matching in dual-band ultrasound transducers. *J. Acoust. Soc. Am.* **141**, 1170–1179 (2017).

40. Fonseca, A. D., Kohler, T. & Ahmed, D. Navigation of ultrasound-controlled swarms under physiological flow conditions. *bioRxiv* (2022).

41. Wang, Y. & Kohane, D. S. External triggering and triggered targeting strategies for drug delivery. *Nat. Rev. Mater.* **2**, 17020 (2017).

42. Mokbel, M. & Aland, S. An ale method for simulations of axisymmetric elastic surfaces in flow. *Int. J. Numer. Methods Fluids* **92**, 1604–1625 (2020).

43. Segers, T., Gaud, E., Versluis, M. & Frinking, P. High-precision acoustic measurements of the nonlinear dilatational elasticity of phospholipid coated monodisperse microbubbles. *Soft Matter* **14**, 9550–9561 (2018).

44. Segers, T., Lassus, A., Bussat, P., Gaud, E. & Frinking, P. Improved coalescence stability of monodisperse phospholipid-coated microbubbles formed by flow-focusing at elevated temperatures. *Lab Chip* **19**, 158–167 (2019).

## Acknowledgements

G.Ch., C.Q., and G.Co. acknowledge funding from CNRS through PEPS Mécanique du futur. G.Co. acknowledges funding from Labex Tec 21 for an outgoing grant. T.S. acknowledges financial support of the Max-Planck Center Twente. Simulations were performed at the Center for Information Services and High Performance Computing (ZIH) at TU Dresden. B.V.e. acknowledges financial support from NWO BUBBLE-X & ULTRA-X-TREME.

## Author contributions

G.Ch. designed the experimental protocol and set-up for buckling and swimming experiments and carried out these experiments. B.V.e. produced and characterized the in-house microbubbles through acoustic attenuation of ultrasound. M.M. performed the buckling simulations. G.Ch., M.V., C.Q., and G.Co. performed the data analysis and interpretation. G.Ch., C.Q., and G.Co. conceived and designed the study. G.Ch., T.S., S.A., M.V., C.Q., and G.Co. drafted the manuscript. All authors read, corrected and approved the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s44172-023-00113-z.

**Correspondence** and requests for materials should be addressed to Georges Chabouh or Gwennou Coupier.

**Peer review information** *Communications Engineering* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Liangfei Tian and Rosamund Daw, Mengying Su. A peer review file is available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023