-
PDF
- Split View
-
Views
-
Annotate
-
Cite
Cite
Fangzhi Mou, Qi Xie, Jianfeng Liu, Shengping Che, Lamya Bahmane, Ming You, Jianguo Guan, ZnO-based micromotors fueled by CO2: the first example of self-reorientation-induced biomimetic chemotaxis, National Science Review, Volume 8, Issue 11, November 2021, nwab066, https://doi.org/10.1093/nsr/nwab066
Close - Share Icon Share
Abstract
Synthetic chemotactic micro/nanomotors are envisioned to actively ‘seek out’ targets by following specific chemicals, but they are mainly powered by bioincompatible fuels and only show pseudochemotaxis (or advanced chemokinesis) due to their weak self-reorientation capabilities. Here we demonstrate that synthetic ZnO-based Janus micromotors can be powered by the alternative biocompatible fuel of CO2, and further provide the first example of self-reorientation-induced biomimetic chemotaxis using them. The ZnO-based micromotors are highly sensitive to dissolved CO2 in water, which enables the corrosion of ZnO to continuously occur by providing H+ through hydration. Thus, they can autonomously move even in water exposed to air based on self-diffusiophoresis. Furthermore, they can sense the local CO2 gradient and perform positive chemotaxis by self-reorientations under the phoretic torque. Our discovery opens a gate to developing intelligent micro/nanomotors powered by, and sensitive to, biocompatible atmospheric or endogenous gaseous chemicals for biomedical and environmental applications.
INTRODUCTION
Biological chemotaxis refers to the ability of those living organisms that are of an advanced level of evolution to move toward or away from specific chemicals [1,2]. To achieve chemotaxis, the living organisms are equipped with delicately coordinated sensing, signal processing and propelling systems, endowing them with the capability to adjust their motion direction on demand by evaluating the direction of chemical gradients at different times (e.g. the temporal comparison adopted by Escherichia coli) or perceiving signal intensity differences across their body (e.g. the spatial comparison adopted by neutrophils) [3,4]. Chemotaxis not only regulates the development of multicellular systems such as tissue development and cancer metastasis, but also enables many unicellular organisms to find nourishments, avoid toxins and coordinate collective motions [5–7].
Inspired by biological chemotaxis, researchers strive to create synthetic chemotactic micro/nanomotors (MNMs). They are envisioned to achieve self-navigation and self-targeting in complex or dynamic environments by actively following specific chemical cues, and thus may bring revolutionary changes to targeted drug delivery, microsurgeries, microfactories, etc [8]. Nonetheless, synthetic MNMs are generally too small to incorporate sophisticated sensing-processing-propelling systems and have a low sensitivity to chemical gradients. Thus, they only show a pseudochemotaxis (or macroscopic chemotaxis and advanced chemokinesis) based on a purely statistical effect rather than chemotaxis of individual MNMs because of the stochastic nature in motion for each individual [8,9]. With the pseudochemotaxis, they show negligible reorientation capabilities, and can only accumulate near (or away from) the chemoattractant (or chemorepellent) source, or shift toward (or away from) its stream in microfluidic channels due to their spatial-dependent mobility (orthokinesis) in the local chemical gradient field [8,10–13]. This statistically average movement uphill (or downhill) to a chemical gradient based on orthokinesis undermines the self-targeting performance of MNMs, and is far from the biological chemotaxis found in nature. On the other hand, chemically powered MNMs have been predicted, in theory, to achieve chemotaxis through their self-reorientation based on the spatial comparison of chemical signal intensity across their bodies, like neutrophils [14,15], but it has not been experimentally confirmed.
In this work, we demonstrate ZnO-based Janus micromotors (MMs) powered by an alternative fuel of CO2, and further provide the first example of self-reorientation-induced biomimetic chemotaxis using them. The ZnO-based MMs can autonomously move based on electrolyte self-diffusiophoresis, even in water exposed to air, as the ZnO component in the MMs is highly sensitive to CO2 dissolved in water, which continuously produces H+ through hydration to react with ZnO. In this way, the MMs can sense local CO2 signals, and exhibit a biomimetic chemotactic behavior toward the CO2 source due to their dynamic reorientations under a phoretic torque. The proposed alternative fuel of CO2, which is a constituent gas in the air and a major product of cell respiration [16,17], has remarkable advantages, including sustained remote accessibility, low cost and excellent biocompatibility. With the excellent biomedical activities of ZnO [18], the intelligent chemotactic behaviors promise that the ZnO-based MMs will ‘seek out’ specific cells or pathogen microorganisms by tracking metabolic CO2 signals emitted from them, and execute targeted biomedical and environmental operations.
RESULTS AND DISCUSSION
To experimentally achieve the biomimetic chemotaxis of artificial MNMs, we at first prepared ZnO-based MMs by sputtering a thin layer of SiO2 on the exposed surface of hexagonal wurtzite ZnO microspheres (Supplementary Fig. 1A and B), as illustrated in Supplementary Fig. 2. The scanning electron microscope (SEM) observation indicates that the prepared ZnO/SiO2 MMs have a relatively uniform spherical shape and an average diameter (d) of 2.5 μm (Fig. 1A). Energy dispersive X-ray (EDX) mapping of Zn and Si elements confirms that the ZnO/SiO2 MMs have a typical Janus structure (the inset in Fig. 1A, and Fig. 1B). The thickness of the sputtered SiO2 layer is ∼150 nm, as verified by the high-contrast SEM observation (Supplementary Fig. 3A) and line-scanning EDX analysis of Zn and Si elements (Supplementary Fig. 3B) on a typical ZnO/SiO2 MM. Zeta potential test shows that the ZnO microspheres before and after coating SiO2 have a similar negative surface charge of 24.1 and −23.7 mV, respectively.
Characterization and self-propulsions of the ZnO-based MMs. (A) SEM image of the ZnO/SiO2 MMs. Scale bar, 1 μm. The inset shows the schematic illustration of their Janus structure. (B) EDX mapping of Zn and Si elements in a ZnO/SiO2 MM. Scale bar, 500 nm. (C) Self-propulsions of ZnO-based MMs in different water media. (i and ii) Trajectories of ZnO/SiO2 MMs in 9 s in the (i) air-exposed water (Air) and (ii) CO2-eliminated water, including the freshly boiled water (Bo), Ar-gassed water (Ar), O2-gassed water (O2) and N2-gassed water (N2). Scale bars, 10 μm. (iii) Trajectories of a ZnO/Pt MM in the air-exposed water without (red curve) and with (green curve) UV irradiation. Scale bars, 5 μm. (D) MSD of the MMs in different water media. Solid color lines are fitting curves using a quadratic function. (E) Self-propulsion of a typical ZnO/SiO2 MM within a lifetime of 90 min. Red curves depict its trajectories in 9 s at different time intervals. Scale bar, 5 μm.
When the ZnO/SiO2 MMs were put into water previously exposed to air (air-exposed water, see ‘Materials’ section in the online Supplementary Data), they exhibited an autonomous motion with an average speed (v) of 7.8 μm/s (Fig. 1C(i) and Supplementary Video 1). As ZnO is poorly reactive or soluble in pure water [19], the self-propulsion of the MMs was attributed to atmospheric CO2 because it has a much higher solubility (34 mM at 25oC and 1 atm) in water compared to other major constituent gases (such as N2, O2 and Ar) [20]. To confirm the key role of CO2 in the self-propulsion of the ZnO/SiO2 MMs, we investigated their motions in CO2-eliminated water (see ‘Materials’ section in the online Supplementary Data), including freshly boiled water (pH = 7.54), Ar-gassed water (pH = 7.46), O2-gassed water (pH = 7.36) and N2-gassed water (pH = 6.67), respectively (Supplementary Video 1). As shown in Fig. 1C(ii) and D, they show an extremely slow motion in the N2-gassed water and only Brownian motions in other CO2-eliminated water, with a translational v and diffusivity (D) lower than 1.3 μm/s and 0.14 μm2/s, respectively. The translational v and D were extracted by fitting their mean square displacements (MSDs, Fig. 1D) using a quadratic function according to the previous report [21].
The SiO2 cap on the ZnO/SiO2 MMs only acts as a passive layer to break the symmetry of the surface reaction of ZnO. This can be verified by the fact that the ZnO microspheres only exhibit Brownian motions in the air-exposed water due to isotropic surface reactions (Supplementary Fig. 4), and the MM can still be self-propelled in the same condition if the SiO2 cap is replaced by a Pt cap (Fig. 1C(iii)). The ZnO/SiO2 MMs showed similar motions in low-energy red light (20.2 mW/cm2) and high-energy UV light (300 mW/cm2) (Supplementary Fig. 5), while the ZnO/Pt MMs were accelerated under the irradiation of UV light because of their high photocatalytic activity (Supplementary Video 2 and Fig. 1C(iii)). This suggests that light illumination has a negligible contribution to the self-propulsion of ZnO/SiO2 MMs. Air-exposed water has a CO2 concentration (CCO2) of 8.54 μM (pH = 5.72), which is three orders lower than the critical concentration of the H2O2 fuel (∼60 mM) required to power traditional catalytic Pt-based MNMs [22]. This suggests that the MMs are highly sensitive to dissolved CO2 in water and can be powered at an ultralow level of the chemical fuel. Moreover, the ZnO/SiO2 MMs exhibit a long lifetime of ∼90 min (Fig. 1E), which is higher than that of the active metal-based MMs (∼5 min), catalytic MnO2 MMs (≤40 min) and some enzyme-based MMs (≤40 min) [23–25].
Due to the different D of the released ionic products of HCO|$_3^{-}$| (D = 1.185 × 10–9 m2/s) and Zn2+ (D = 0.703 × 10–9 m2/s) from the ZnO/SiO2 MMs according to equation (3), HCO|$_3^{-}$| would diffuse faster than Zn2+ [32,33], leading to an uneven distribution of these ions in the area near the exposed ZnO surface (Supplementary Fig. 6A and B). To maintain electroneutrality, a local diffusioelectric field (E) arises to slow down the faster-diffusing ion, and speed up the slower ion [34], as verified by the simulated electric potential ϕ depicted in the color background of Fig. 2A and the local E shown in Supplementary Fig. 6C (black triangles). This local E then induces an electro-osmotic slip (EOS) [35] in the electric double layer of the negatively charged MM (ζ = −23.7 mV). The simulated flow field (EOF) initiated by the surface EOS around a ZnO/SiO2 MM is illustrated as black streamlines with arrows in Fig. 2A, revealing that the MMs are self-propelled based on electrolyte self-diffusiophoresis with the ZnO end forward. This propulsion direction was further experimentally verified by the close observation of a partially etched ZnO/SiO2 MM and the barely etched SiO2/ZnO MMs, both of which have a distinguishable ZnO end at a high magnification under an optical microscope (Supplementary Figs 7 and 8, Supplementary Video 3). The detailed propulsion mechanism of the ZnO/SiO2 MM fueled by the dissolved atmospheric CO2 in water is illustrated in Fig. 2B.
The self-propulsion mechanism of the ZnO/SiO2 MMs. (A) Numerical simulations of the electric potential (ϕ, the color background) and the electro-osmotic flow (EOF, black streamlines with arrows) around a ZnO/SiO2 MM. The EOF is initiated by the surface electro-osmotic slip (EOS) in its electric double layer. (B) Schematic illustration of the propulsion mechanism of the ZnO/SiO2 MMs. (C) Trajectories of the MMs in freshly boiled water at different time intervals when exposed to air. Scale bar, 5 μm. (D) Time-dependent normalized v of the MMs in a film of freshly boiled water when exposed to air (red dots) and the air-exposed water (black dots). Solid red and black curves are fitting results using exponential functions. (E) The simulated CCO2 as a function of time in a CO2-eliminated water film (160 μm in thickness) when exposed to air. (F) The v of the MMs at different CCO2 in water. The red line is the linear fitting line of the experimental data in the low CCO2 regime (0–8.54 μM), and the red curve in the inset is the fitting result according to the Michaelis-Menten kinetics in the high CCO2 regime (8.54 μM to 20.5 mM).
From the propulsion mechanism (Fig. 2B), it is expected that the ZnO-based MMs can still acquire fuel from the air even if no CO2 was pre-dissolved in water. To verify this assumption, we recorded the motion behaviors of the ZnO/SiO2 MMs in a thin film (∼160 μm in thickness) of the freshly boiled water over time when exposed to air (Fig. 2C and Supplementary Video 4). Differently to the MMs in the air-exposed water film (CCO2 = 8.54 μM), which show a decreasing v over time (the black curve with dots in Fig. 2D), those in the CO2-eliminated water display an increasing v in the first 7 min (the red curve with dots in Fig. 2D). This behavior can be verified by the different lengths of trajectories of a typical MM in 8 s (Fig. 2C) and the statistic normalized v of the MMs (Fig. 2D) at different intervals within 7 min. This increase in v reveals that atmospheric CO2 molecules can diffuse into the water as fuel to power the MMs in a sustainable and rapid manner, as verified by the numerical results of the increasing CCO2 (Fig. 2E and Supplementary Fig. 9A) and H+ concentration (|${C_{{{\rm{H}}^ + }}}$|, Supplementary Fig. 9B) over time when the water film is exposed to air. In the following 7–51 min, the v of the MMs decreases gradually due to the decreasing reaction rate of ZnO and the increasing ion concentration in the water film. With the increase of CCO2 from 0 to 8.54 μM, the v of the ZnO/SiO2 MMs increases from 2.0 to 7.8 μm/s (Fig. 2F and Supplementary Video 5). By bubbling water with pure CO2 gas, a higher CCO2 compared to that in the air-exposed water can be obtained, and the v of MMs can reach 41.5 μm/s when CCO2 increases to 20.5 mM (the inset in Fig. 2F and Supplementary Video 5). The fitting results show that the v of the MMs increases linearly in the low CCO2 regime from 0 to 8.54 μM (red line in Fig. 2F), and gradually reaches to a plateau in the high CCO2 regime (8.54 μM–20.5 mM) following Michaelis-Menten kinetics (red curve in the inset of Fig. 2F). These results reveal that the ZnO/SiO2 MMs are of a high sensitivity to CCO2, especially in the low CCO2 regime. The sensitivity can be quantified by the rate of the speed change of the ZnO/SiO2 MMs to CCO2, which has a value of 0.57 μm/(μM·s) (the slope of the red line in Fig. 2F) in the low CCO2 regime, about five orders of magnitude higher than that of common H2O2-fueled MMs (2.72 × 10–6 μm/(μM·s)) [36]. The appearance of the plateau of v at a high CCO2 of 20.5 mM indicates that the maximum speed is limited by the maximum reaction rate depending on total surface reactive sites.
Benefiting from the high sensitivity to CCO2, the ZnO/SiO2 MMs can sense a local CO2 gradient and exhibit positive chemotaxis toward a local CO2 source in water (Fig. 3A). When a micropipette filled with a CO2 solution (CCO2 = 93.7 μM) was introduced to the water medium (CCO2 = 0.2 μM), a gradient field of CCO2 (|$\nabla C$|) is formed (the color background in Fig. 3A). In this gradient field, the ZnO/SiO2 MMs swarm toward the micropipette, in analogy to the positive chemotaxis of living organisms in nature (Fig. 3B and Supplementary Video 6). The average chemotactic velocity, defined as the net displacement of the ZnO/SiO2 MMs divided by time when they are approaching the region close to (∼20 μm) the CO2 source, is 0.57 μm/s. In contrast, passive SiO2 microspheres (d = 2 μm and ζ = −18.5 mV) show no chemotactic motion in the same gradient field (Supplementary Fig. 10 and Supplementary Video 7), indicating that the MMs were actively swimming toward the CO2 source, and the possible passive diffusiophoretic attraction generated by the gradient field is negligible.
Chemotaxis of the ZnO/SiO2 MMs. (A) Schematic illustration of the positive chemotaxis of the ZnO/SiO2 MMs. The color background depicts numerical simulations of CCO2 gradient (|$\nabla C$|) near a micropipette (inner d, 20 μm) filled with a CO2 solution (CCO2 = 93.7 μM). (B) Trajectories of the MMs moving toward the CO2 source. Scale bar, 20 μm. (C) The color trajectory of a slightly etched ZnO/SiO2 MM approaching the CO2 source, in which the colors represent the normalized magnitude of v of the MM. Scale bar, 5 μm. (D) The reorientations of the ZnO/SiO2 MM shown in (C) at different times by clockwise (CW) or counterclockwise (CCW) rotations. Scale bar, 2 μm. (E) Numerical simulations of the EOSs (purple arrows) on a ZnO/SiO2 MM and the resulted flow fields (black streamlines with black arrows) when it aligns with |$\nabla C$| (θ = 0, the left panel) and deviates with it in a θ of π/4 (the right panel), respectively. The color background depicts the magnitude of the flow velocity. The green arrow represents the phoretic torque (Mp) induced by the unbalanced EOSs. (F) The calculated Mp at different θ. (G) Time-lapse images displaying the change of N near the CO2 source. Scale bar, 20 μm.
To quantify the degree of chemotaxis, we have analyzed the chemotaxis index (CI) of ZnO/SiO2 MMs, which was defined as the ratio of the displacement in the direction of the chemical gradient to the total migration distance (i.e. as the cosine of the angle θ of the MM’s motion direction with the chemical gradient) [37]. The average CI is 0.21, suggesting that they, in general, move uphill to the chemical gradient. Due to the persistent ZnO-forward configurations during the propulsion (Supplementary Figs 6 and 7), the self-reorientation dynamics of the ZnO/SiO2 MMs can be characterized by the changes of CI and the angle θ over time. When CI > 0 (i.e. |$0 \le | \theta | < {\rm{\pi }}/2$|), the ZnO/SiO2 MMs prefer to move up the gradient, while they tend to move down the gradient if CI < 0 (i.e. |${\rm{\pi }}/2 < | \theta | \le {\rm{\pi }}$|). The CI and |θ| of four typical ZnO/SiO2 MMs versus time (Supplementary Fig. 11) suggest that they move up the gradient for a longer time than moving down. The longer upgradient motion time implies that the ZnO/SiO2 MMs prefer to move with an orientation along the upgradient direction even though their motion direction was seriously randomized by the Brownian rotational diffusions.
To further demonstrate the biomimetic chemotactic behavior, the instantaneous v, orientation and rotation of a slightly etched ZnO/SiO2 MM, in which the slightly etched ZnO end was distinguishable under an optical microscope, were recorded and analyzed when it was approaching the CO2 source (Fig. 3C and D). The ZnO/SiO2 MM shows an increasing v when approaching the CO2 source because of the increasing CCO2 (i.e. orthokinesis of the motor [38]), as verified by the normalized instantaneous v in 0–15 s and 15–30 s (Fig. 3C), respectively. In addition, when the ZnO/SiO2 MM deviates with |$\nabla C$| because of the strong perturbations of Brownian randomizations, it dynamically reorientates itself by clockwise or counterclockwise rotations (Fig. 3D). When the MM moves with an orientation down |$\nabla C$|, it rotates clockwise gradually and aligns with |$\nabla C$| (0–16.5s in Fig. 3D). After a sharp turn (16.5–20 s in Fig. 3D), the MM deviates with |$\nabla C$| again, and then gradually aligns with |$\nabla C$| by a counterclockwise rotation (20–30 s in Fig. 3D). This self-reorientation behavior bears resemblance to that of chemotactic microorganisms which show a ‘memory’ of the direction of the chemical gradient and move up the gradient of chemoattractant in a ‘deviating-rectifying’ manner. This is the first experimental evidence showing that synthetic MNMs can achieve chemotaxis by dynamic self-reorientations in a gradient field [8]. After reaching the area close to the tip of the micropipette (CO2 source), the MMs were wandering in this area with no preferred directions, as depicted by the curved trajectory of the slightly etched MM in 34–163 s in Supplementary Fig. 12, so that they could accumulate near the CO2 source (Supplementary Video 6).
To decipher the chemotactic behavior, we have simulated the flow field around the ZnO-based MM when a local CO2 source is introduced, as shown in Fig. 3E. The results show that, if the axis of the MM is aligned with the CO2 gradient, EOSs (e.g. us1 and us2 in the left panel in Fig. 3E, us1 = us2) across the axis of the MM are symmetrical because of symmetric surface reactions, and the MM migrates up the gradient at this condition. However, due to the strong rotational Brownian diffusion, the symmetry axis of the MM is constantly deviating from |$\nabla C$|. When misaligned with |$\nabla C$| with an angle of θ, the ZnO surface close to the CO2 source has a higher reaction rate (due to the higher CCO2 (or |${C_{{{\rm{H}}^ + }}}$|)) than the surface further away from the source, resulting in the asymmetric distribution of chemical products and thus unbalanced EOSs across the motor axis (e.g. us1 and us2 in the right panel of Fig. 3E, us1 > us2). The unbalanced EOSs then produce a phoretic torque (Mp) to induce rotation of the MM until its axis is aligned with |$\nabla C$| (the right panel of Fig. 3E). With the increasing degree of deviation, the self-EOSs across the motor axis become more unbalanced, and thus Mp increases accordingly for reorientation, as verified by the increasing Mp with a varying θ from 0 to -π/2 and π/2 (Fig. 3F). Hence, the chemotactic behavior of the MM is attributed to dynamic self-reorientations under the phoretic torque. Due to their self-reorientations, they show a higher efficiency in chemotaxis compared to the developed (pseudo)chemotactic MNMs, which usually take hours to aggregate near the chemical source [10,39,40]. As shown in Fig. 3G, the MMs swarm to the chemical source in several minutes, as demonstrated by the increasing number (N) of MMs from 0 to 15 within 284 s in a sector region close to (80 μm) the CO2 source.
CO2 not only exists in the air and various aqueous media (such as natural waters and biological media), but also can be produced and secreted by cells or pathogenic microorganisms during their oxidative metabolism [41–44]. For instance, the CO2 tension (PCO2) of the intracellular and extracellular fluids is ∼50 and 46 mmHg, indicating that there is a high CCO2 (∼20.8 and 19.3 mM together with the chemically bound CO2) and a steep CCO2 gradient near tissue cells (e.g. ∼3.3 × 104 mol/m4 at a distance of 20 μm to the cell surface), respectively [45]. On the other hand, ZnO micro/nanomaterials are of low toxicity and good biodegradability and have shown great potential in biomedical applications due to their excellent anti-cancer, anti-diabetic, anti-bacterial and anti-inflammatory activities [18]. Thus, the chemotactic ZnO-based MMs with proper surface modifications (e.g. polyelectrolyte ion-tolerance coating [46]) are expected to act as biomimetic microrobots, ‘seeking out’ specific cells or pathogenic microorganisms by actively tracking their extracellular CO2 signals, thereby executing targeted biomedical and environmental operations.
CONCLUSION
In conclusion, we have demonstrated that ZnO-based MMs can be powered by CO2 fuels and perform intelligent positive chemotaxis by self-reorientations. We analyzed the propulsion mechanism of ZnO-based MMs and confirmed that it was the result of electrolyte self-diffusiophoresis based on the corrosion of the ZnO end by H+ continuously supplied from the dissolved CO2. Benefiting from their high sensitivity to local CO2 signals, the ZnO-based MMs can perform a biomimetic chemotactic behavior toward the CO2 source due to their dynamic self-reorientations in the chemical gradient field. The self-reorientation of the MM is attributed to the phoretic torque produced by the unbalanced EOSs across the MM axis when it is misaligned with the chemical gradient. Due to the readily accessible biocompatible CO2 fuel and intelligent chemotaxis, ZnO-based MMs have the potential to autonomously move in a wide range of aqueous media and execute targeted biomedical and environmental operations by actively tracking respiratory CO2 signals emitted from cells or pathogenic microorganisms.
METHODS
Preparation of ZnO-based MMs
The ZnO microspheres were synthesized according to the previous report [47]. To prepare ZnO/SiO2 MMs, the ZnO microspheres on the glass slide were partially coated with a SiO2 nanolayer by magnetron sputtering for 40 min (JCP500, Beijing Technol Science Co., Ltd, China). The substrate was heated to 100oC and rotated at a speed of 30 r/min during the sputtering process. In addition, ZnO/Pt MMs and SiO2/ZnO MMs were also prepared by sputtering Pt on the ZnO microspheres on the glass slide for 30 s and a ZnO layer on monodispersed SiO2 microspheres (2 μm) for 60 min, respectively.
Self-propulsion and chemotaxis
To observe the self-propulsion of ZnO/SiO2 MMs, a 20 μL aqueous suspension of the MMs (0.03 mg/mL) was dropped into a Petri dish (d = 35 mm) with 2 ml water. The self-propulsion of the MMs was observed and recorded through an inverted optical microscope (Leica DMI 3000B). All videos were analyzed using ImageJ and Video Spot Tracker V08.01 software. More than 10 MMs were analyzed to obtain the statistical result. To investigate the MM’s chemotaxis, a micropipette (inner d, 20 μm) filled with a CO2 solution (CCO2 = 93.7 μM) was fixed on the holder of a high-precision micromanipulator (Leica Microsystems, Germany) and then put in a water droplet (50 μL, pH = 6.52) with the ZnO/SiO2 MMs (0.003 mg/mL).
Numerical simulations
The simulations were performed by using the diffusions, electrostatics and creeping flow modules of COMSOL Multiphysics software [21]. Governing equations for numerical simulations are given in the online Supplementary Data. The simulation model was built up by immersing a ZnO/SiO2 MM in the middle of a cylindric cell (radius r = 50 μm, height h = 100 μm) filled with the air-exposed water (CCO2 = 8.54 μM). The fluxes (J) of Zn2+, H+ and HCO|$_3^{-}$| from the ZnO surface were set as 0.05, −0.1 and 0.1 mmol m–2 s–1, respectively. The D of Zn2+, CO2, H+ and |${{\rm HCO}_3^-}$| were set as 0.70 × 10–9, 1.77 × 10–9, 9.31 × 10–9 and 1.19 × 10–9 m2/s, respectively. The ζ of the MM was set as −23.7 mV. To simulate the reorientation of the MM during chemotaxis, a CO2 source (CCO2 = 93.7 μM) was added in the above model and placed at the left boundary with an r of 10 μm and a distance of 50 μm to the center, and J of Zn2+ and |${{\rm HCO}_3^-}$| at each point of the exposed ZnO surface was set to be proportional to the local |${C_{{{\rm{H}}^ + }}}$|. To simulate the diffusion kinetics of atmospheric CO2 (0.04 vol.% in the air) into water through the interface with the atmosphere, a simulation model of a water film (160 μm in thickness) was built and the classical two-film theory was used [48], in which the mass transfer coefficient (k) of CO2 in the water was set to be 2 μm/s [49], and the equilibrium concentration (C*) of CO2 at the interface was calculated to be 1.2 × 10–5 M according to Henry's Law equation.
Acknowledgements
The authors acknowledge Mr Kong Xiong, Dr Leilei Xu and Dr Ming Luo for their valuable discussion on the manuscript.
FUNDING
This work was supported by the National Natural Science Foundation of China (52073222, 21875175, 21975195, 21705123 and 51521001), the Natural Science Foundation of Hubei Province (2019CFA048), and the Fundamental Research Funds for the Central Universities (WUT: 2021IVA118).
AUTHOR CONTRIBUTIONS
F.M. and J.G. proposed and supervised the project. F.M. conceived and designed the experiments. F.M., Q.X., J.L. and M.Y. carried out the synthesis, structural characterization and self-propulsion tests. F.M. and S.C. carried out theoretical calculations and numerical simulations. F.M., Q.X., J.L., L.B. and J.G. co-wrote the manuscript. All authors discussed the results and participated in analyzing the experimental results.
Conflict of interest statement. None declared.
REFERENCES
Author notes
Equally contributed to this work.


