Preparation and characterization of the sperm-like nanorobot. The fabrication process of sperm-like nanorobots was briefly shown in Fig. 1a. Firstly, we prepared flexible Au/PPy composite nanowires with the Au ends modified using streptavidin, which could act as artificial flagella of the nanorobots to generate undulatory wave propagation. Here, multi-step electrochemical deposition and monolayer modification were combined for fabrication. Au segments of 3 µm and PPy tails of given lengths were deposited in sequence. After mechanical polishing of the sacrificial gold layer, the Au tips of the nanowires were exposed and coated with monolayer of DTSP molecules, which could increase the contact area and improve the assembly efficiency for the following process. Then wet-etching of the Al2O3 templates was conducted and the nanowires were released to be modified with streptavidin. After that, the nanowires were mixed with biotin-coated Fe3O4 nanoparticles for bonding and formed the final sperm-like nanorobots. In this case, the Fe3O4 head endows the nanorobot with magnetic response and the flexible flagellum ensures undulatory propagation for effective propulsion.
Morphology and elemental analysis of the artificial flagella were conducted, and the SEM and EDS results were shown in Fig. 1b and 1c. Distinct two-stage structure can be observed on the nanowire, which corresponds to the Au tip and the PPy body, respectively. Flexibility of the artificial flagellum is presented with slightly bending in this view. Schematics of the self-assembled nanorobot’s head-to-tail structure and the corresponding SEM image were shown in Fig. 1d and 1e, respectively. It can be observed that an intact self-assembled microstructure has been successfully achieved between the magnetic head and the flexible flagellum, to form the sperm-like nanorobot.
Theoretical analysis of the nanorobot’s propulsion. The nanorobot in this study is fabricated via self-assembly using a superparamagnetic nanobead and an artificial flexible ultra-thin nanorod. Assuming the head has a magnetic moment M in the direction under the magnetic field, the flexible flagellum has a body length L with a given elasticity modulus. For an arbitrary point P along the tail, we mark it with a generalized coordinate q respective to the head and describe its position at time t with vector p(q, t), which can also be transferred to the head’s frame as
\(p\left(q,t\right)=(q\hspace{1em}{\varphi }_{y}(q,t)\text{}\hspace{1em}{\varphi }_{z}(q,t\left)\right)\) | (1) |
where ϕy(q,t) and ϕz(q,t) represent the corresponding deformation along Y and Z axis, respectively. |
Then the moving Frenet-Serret frame can be established on point P along local tangent, normal, and bi-normal directions:
\(t=\frac{dp(q,t)/dq}{\left|\right|dp(q,t)/dq\left|\right|}, n=\frac{dt/dq}{\left|\right|dt/dq\left|\right|}, b=t\times n\) | (2) |
The head can orient itself to the time-varying magnetic field due to the magnetic torque exerted on the dipole. Due to the head’s rigid body rotation, bending waves are generated and can be propagating along the flexible flagellum. According to the typical resistive force theory, the force and torque balance based on magnetic and fluidic fields can be expressed as:
\(\left[\begin{array}{c}{F}_{m}+{F}_{d}\\ {T}_{m}+{T}_{d}\end{array}\right]=0\) | (3) |
Here, gravity is neglected, Fm and Tm are magnetic force and torque that can be expressed as: |
\(\left[\begin{array}{c}{F}_{m}\\ {T}_{m}\end{array}\right]=\left[\begin{array}{c}V(M\cdot \nabla )B\\ VM\times B\end{array}\right]\) | (4) |
F d and Td are fluidic drag force and torque, which can be decomposed into contributions of the head and the tail:
\(\left[\begin{array}{c}{F}_{d, head}\\ {T}_{d, head}\end{array}\right]=\left[\begin{array}{cc}{D}_{h}& -{D}_{h}{S}_{h}\\ {S}_{h}{D}_{h}& {N}_{h}\end{array}\right]\left[\begin{array}{c}V\\ \varOmega \end{array}\right]\) | (5) |
\(\left[\begin{array}{c}{F}_{d, tail}\\ {T}_{d, tail}\end{array}\right]={\int }_{0}^{L}\left[\begin{array}{cc}RC{R}^{T}& -RC{R}^{T}{S}_{t}\\ {S}_{t}RC{R}^{T}& -{S}_{t}RC{R}^{T}{S}_{t}\end{array}\right]\hspace{0.33em}dl\left[\begin{array}{c}V\\ \varOmega \end{array}\right]\) | (6) |
where V and Ω are propulsion and rotation velocity matrices, Dh and Nh are resistance matrices for the head, and C represents the resistance coefficient matrix for the tail. Besides, R is the rotation matrix between the local Frenet-Serret coordinates and the nanorobot’s frame, Sh and St are the corresponding position transformation matrices for the head and the tail.
The above equations describe the elastohydrodynamics of the sperm-shaped nanorobot swimming in viscous fluids 27. In the experiments, the nanorobot was actuated to swim in the horizontal plane (XY) and was investigated from the top view. Apart form the elastohydrodynamics, a simplified undulatory propagation modelling could also be used 33,34. Swimming of the nanorobot in the X direction can be simplified into two orthogonal components in the XY and XZ planes. Motion of the nanorobot in each plane can be treated as a bending wave propagation across the two distal ends, and relative phase lag between the head and tail in two planes are noted as ϕxy and ϕxz. Here, flagellar propulsion of the sperm-like nanorobot along X direction can be ascribed to the overall oscillation along the body induced by two separate oscillations in both XY and XZ planes as:
\(y(x,t)={a}_{1}{sin}(2\pi ft)\left(\frac{x}{L}\right)+{a}_{2}{sin}(2\pi ft+{\phi }_{xy}){\left(\frac{x}{L}\right)}^{m}\) | (7) |
\(z(x,t)={a}_{3}{sin}(2\pi ft)\left(\frac{x}{L}\right)+{a}_{4}{sin}(2\pi ft+{\phi }_{xz}){\left(\frac{x}{L}\right)}^{m}\) | (8) |
where a1 and a3 are amplitudes of magnetic oscillation that correspond to the head, a2 and a4 are amplitudes of fluidic oscillation that correspond to the tail. Besides, f is the actuation frequency, and m is the curvature of bending deformation.
The propulsive forces in two planes can be expressed as
\({F}_{x}\left(xy\right)={\int }_{0}^{L}{C}_{\parallel }\left[\left(\frac{{C}_{\perp }}{{C}_{\parallel }}-1\right)\frac{dy}{dt}\frac{dy}{dx}-{V}_{x}\left(xy\right)\right]dl\) | (9) |
\({F}_{x}\left(xz\right)={\int }_{0}^{L}{C}_{\parallel }\left[\left(\frac{{C}_{\perp }}{{C}_{\parallel }}-1\right)\frac{dz}{dt}\frac{dz}{dx}-{V}_{x}\left(xz\right)\right]dl\) | (10) |
where dl is an infinitesimal section along the body, \({C}_{\parallel }\) and \({C}_{\perp }\) are the drag coefficients in the tangential and normal directions. Besides, \({V}_{x}\left(xy\right)\) and \({V}_{x}\left(xz\right)\) are two contributions of the X-direction velocity resulted from oscillations in XY and XZ planes, respectively. When the nanorobot reaches a steady swimming state, the total force equals zero, and the velocity components can be calculated to be
\({V}_{x}\left(xy\right)=\pi \left(\frac{{C}_{\perp }}{{C}_{\parallel }}-1\right)\left(\frac{m-1}{m+1}\right)\frac{{a}_{1}{a}_{2}}{Lf}{sin}\left({\phi }_{xy}\right)\) | (11) |
\({V}_{x}\left(xz\right)=\pi \left(\frac{{C}_{\perp }}{{C}_{\parallel }}-1\right)\left(\frac{m-1}{m+1}\right)\frac{{a}_{3}{a}_{4}}{Lf}{sin}\left({\phi }_{xz}\right)\) | (12) |
And the resultant velocity along X-direction can be expressed as |
\({V}_{x}={V}_{x}\left(xy\right)+{V}_{x}\left(xz\right)\) | (13) |
Hence, we can deduce that a phase difference in oscillation is requisite to generate effective locomotion with a nonzero velocity. The oscillations in two perpendicular planes undergo different fluidic drag conditions due to the substrate-induced asymmetry, and diverse amplitudes as well as phase lags are stimulated. The resultant motion of the nanorobot is determined by combination of the decomposed oscillations, which can result in a balance between propulsion and retarding in the form of forward or backward locomotion. In this case, dynamic magnetic field with a precessing angle is essential for swimming direction reversal since spatial oscillations become accessible to lead to bidirectional resultant.
Propulsion performance of the sperm-like nanorobot. To actuate the sperm-like nanorobot and test its propulsion performance, a precessing magnetic field or a so-called conical rotating magnetic field was applied using a custom-made triaxial Helmholtz coil system. The externally actuated magnetic field can be expressed as
\(\overrightarrow{B}={B}_{yz}\left[{cos}(2\pi ft\right)\overrightarrow{{e}_{y}}+{sin}(2\pi ft\left)\overrightarrow{{e}_{z}}\right]+{B}_{x}\overrightarrow{{e}_{x}}\) | (14) |
where Byz and f are the amplitude and rotating frequency of the circular polarized component of the magnetic field. And Bx is the static component of the magnetic field. Besides, \(\overrightarrow{{e}_{x}}\), \(\overrightarrow{{e}_{y}}\), \(\overrightarrow{{e}_{z}}\) are the unit vectors along the corresponding coordinate axis.
The resultant field vector at the head of the nanorobot was rotating in a cone-like path, and the angle between the field vector and the cone axis was noted as the precessing angle (θ) of the actuation field. For flagellar propulsion of sperm-like nanorobots, undulatory motion theory has been developed as described above, which is focused on the bending wave propagating process along the body. Actuated by the magnetic torque exerted on the head, continuous rotation is triggered to drive the body to fluctuate. The stimulated wave propagates along the flagellum yet different oscillating amplitudes are generated from the head to the tail. Flexibility of the artificial flagellum contributes to effective undulatory propagation and a constant phase lag exists compared to the magnetic head. Here, time-dependent deformation of the flexible flagellum was measured over one whole period and morphing changing process of the body was also illustrated. The nanorobot was actuated to oscillate by a 2 Hz precessing field (B = 100 Gs, θ = 30°) (Movie 1). The results indicate the whole propagating process of the generated wave during helical flagellar propulsion, and a distinct oscillating angle can be observed from the superimposed image as shown in Fig. 2a. Such nanorobot could be actuated to precess synchronously with the low-frequency alternating magnetic field.
A planar coordinate system was defined to facilitate position recording and calculation. In this view, the upward tail tip of the nanorobot corresponds to a positive value of oscillating angle yet the downward tail tip corresponds to a negative one. As shown in Fig. 2b, we measured the oscillating angle change of the nanorobot in a given time period, which exhibited a sinusoidal pattern with time and kept consistent with the precessing field’s oscillation. To investigate the relative position relations between the nanorobot’s head and flagellum, we further recorded the displacements of the head and the tail tip along the Y axis, with respect to the original position (Y = 0) at t = 0 s. As shown in Fig. 2c, the displacement waveforms near the distal end as well as the head were recorded. Both the head and the flagellum undergo a sinusoidal waving locomotion in the direction perpendicular to the cone axis, and a fixed phase diversity exists between them. The head and the flagellum exhibit diverse oscillation amplitudes during propulsion, in which the tail tip oscillates with a larger amplitude compared with the head. In this case, the nanosphere head dominates the magnetic response and acts as an oscillation source. It is actuated to precess around and drives the flagellum to rotate in the same pattern with a phase lag under the sinusoidal alternating magnetic field. It also proves that asymmetrical shape deformation occurs to cater the scallop theorem and achieve effective locomotion as a result. At the same time, a negative positional offset of the swimmer’s body can be observed, which is related to the minor lateral drift during swimming.
Bidirectional locomotion property of the nanorobot. For magnetic actuation of the sperm-like nanorobot, four types of spatially oriented precessing magnetic fields were defined according to the corresponding directions (Fig. 3a). When the field vector exhibited an acute angle with the flagellum-head direction, the field was defined as a “Head-pointing Field” (HF), and the rotating direction was further classified into counter-clockwise (CCW) or clockwise (CW) types on the basis of the right-hand rule. Schematics of the sperm-like nanorobot actuated by a HF type field with a precessing angle θ was shown in Fig. 3b. Similarly, when the field vector was pointing in the head-flagellum direction, it could be defined as a “Flagellum-pointing Field” (FF) with a CCW or CW rotating directions. In our experiments, the nanorobotic sperms would not easily turn around for coincident magnetic alignment when the precessing axis of the high-frequency field (f > 5 Hz) was abruptly turned in an opposite direction. This could be ascribed to the relatively low magnetic torque as well as the encountered fluidic drag of the slender nanorobot near the substrate surface, which induced the insensitive magnetic response under high-frequency dynamic fields. This was also confirmed using magnetic nanowires actuated under 3 Hz and 5 Hz precessing fields for comparison (Movie 5 and 6). Both HF and FF type fields could be used to actuate the nanorobot and the flexible flagellum could be driven to oscillate around with the precessing magnetic head. Effective locomotion can be observed near the substrate due to periodic nonreciprocation that breaks the time symmetry, which is depended on integrated effects of elastic flagellum, hydrodynamic resistance as well as fluidic flows generated by the moving head. Due to continuous interactions with the viscous fluids, helical flagellar propulsion was generated for the nanorobot.
To study the motion directionality, diverse precessing angles and directions were adopted for experimental comparison (Movie 2). As shown in Fig. 3c, the field strength and frequency were set to be constant at 100 Gs, 40 Hz, yet the spatial directions were different. In Fig. 3c-1, when a field of HF-CCW at θ = 10° was applied, backward locomotion in the flagellum-pointing direction could be observed, which corresponded to a negative velocity (V < 0). However, when the field was set to be precessing in the CW direction, the nanorobot changed to swim forward on the contrary. In this case, the nanorobot was capable to move in an opposite direction without a U-turn trajectory, which was different from other one-tail magnetic microswimmers. This property could be attributed to the chirality of the sperm-like nanorobot, which was a combined effect of the nonideal self-assembled structural configuration as well as the time-dependent undulatory dynamics. Thus magnetic fields of diverse precessing directions were capable to drive the nanorobot to rotate and achieve bidirectional locomotion. In Fig. 3c-2, the magnetic field (θ = 10°) was tuned to be a FF type, the nanorobot kept the original pointing direction without turning around. As for magnetic field transformation, a HF-CW type field could be tuned to FF-CW type as just turning the precessing axis 180°, yet the nanorobot could be driven to oscillate and move in an opposite direction as a result. Therefore, a FF-CW and FF-CCW type of precessing magnetic field could actuate the nanorobot to move in the same direction with the results of HF-CCW and HF-CW cases, respectively.
In addition, the influence of the actuation field’s precession angle was also investigated as shown in Fig. <link rid="fig3">3</link>c-3. Under a HF type field at θ = 30°, the nanorobot could be effectively actuated forward in both CCW and CW directions, which was quite different with the results demonstrated in Fig. 3c-1. It could be explained that actuation field of a large precessing angle induced relatively violent precession with higher oscillation frequency, and the chirality was no longer valid in this case. The results confirmed that flexible switching between swimming forward and backward could be achieved via tuning the actuation field parameters, including precessing direction and angle as well as field frequency. In the future, such bidirectional propulsion property can be applied to control multiple sperm-like nanorobots with distinguished strategies and complete cooperative tasks.
Based on the preliminary actuation study above, actuation experiments under fields of various parameters were conducted, and locomotion directions as well as resultant velocities were systematically measured. As shown in Fig. 4a-1 to 3, the field strength was set to be constant at 50 Gs, yet the precessing angles were 10°, 30°, 60°, respectively. Four types of fields (HF-CCW, HF-CW, FF-CCW, FF-CW) at a given frequency ranging from 5 Hz to 40 Hz were applied to actuate the nanorobot. The positive velocity corresponded to swimming forward yet the negative velocity represented backward propulsion. At a lower value of precession angle (θ = 10°), the nanorobot was actuated to swim backward over the frequency range regardless of spatial directions of the fields. However, the results at higher precession angles (θ = 30° and 60°) were completely different since both forward and backward locomotion could be observed. Specifically, the nanorobot was actuated forward under HF-CW or FF-CCW fields, yet backward under HF-CCW or FF-CW fields. A highly effective forward velocity of 2.11 µm/s and backward velocity of 2.77 µm/s could be easily obtained at θ = 30° via tuning the precession directions of the field, indicating distinct bidirectional propulsion property which was different from previously reported sperm-shaped microswimmers.
As shown in Fig. 4b-1 to 3, under precessing magnetic fields of 70 Gs, the propulsion results were similar with the cases in 50 Gs when other parameters kept the same. The nanorobots still kept swimming backward at θ = 10°, yet locomotion differentiation occurred when θ = 30° or 60°. However, when the field strength was large enough (100 Gs), it turned to be much easier for the nanorobots to move forward. Specifically, both HF-CW and FF-CCW fields could successfully drive the nanorobots forward at θ = 10°, despite backward locomotion was still destined under HF-CCW or FF-CW fields. A large forward velocity of 2.96 µm/s and backward velocity of 4.26 µm/s could both be achieved at 40 Hz using the opposite precessing directions (HF-CW and FF-CW fields, respectively). As increasing precessing angle of the 100 Gs fields (θ = 30°, 60°), robust forward propulsion turned to be dominated regardless of the spatial directions of the actuation fields. In this case, a significantly high forward velocity could be obtained over a wide frequency range.
Thereon, we measured and summarized locomotion directions of the nanorobot actuated by precession magnetic fields of diverse parameters, and the results were listed in Table 1. From the observations in our experiments, it can be deduced that locomotion direction of the sperm-like nanorobot is determined by multiple actuation parameters including field strength, frequency, precessing angle as well as direction. In the experiments, the sperm-like nanorobots generally tended to move forward under precessing magnetic fields of HF-CW or FF-CCW types. To propel the nanorobots forward, precessing magnetic fields of an intensified strength as well as a large precessing angle were desired for actuation, which was also consistent with the velocity measurement results in Fig. 4.
Table 1
Locomotion directions of the sperm-like nanorobot actuated by a precessing magnetic field of given parameters over a frequency range (5–40 Hz). Here, “+” represents forward locomotion yet “-” represents backward locomotion.
θ | B (Gs) | HF-CW/HF-CCW | FF-CW/FF-CCW |
10° | 50 | -/- | -/- |
70 | -/- | -/- |
100 | +/- | -/+ |
30° | 50 | +/- | -/+ |
70 | +/- | -/+ |
100 | +/+ | +/+ |
60° | 50 | +/- | -/+ |
70 | +/- | -/+ |
100 | +/+ | +/+ |
Moreover, we also fabricated nanorobots with shorter flagella, which were synthesized via controlling the electrochemical deposition time. In Figure S2, SEM image of the short-tailed nanorobot (tail length about 8 µm) was demonstrated. Actuation experiments of such nanorobot under precessing magnetic fields of 100 Gs were carried out (Movie 3 and 4) and locomotion velocities were calculated as shown in Fig. 5. Similarly, four directional types of fields and precessing angles of 10°, 30°, 60° were applied successively for flagellar propulsion. Bidirectional locomotion could be distinctly observed under fields of any precession angles, which was quite different compared with the long-tailed nanorobots that tended to move forward under the same actuation condition. Under the 100 Gs magnetic field of a moderate precessing angle (θ = 30°), typical bidirectional property could be observed and relatively high locomotion velocity exceeded 2 µm/s could be achieved both in forward and backward directions. The nanorobot could be actuated to move forward under precessing magnetic fields of HF-CW or FF-CCW types, and move backward under HF-CCW or FF-CW fields, which were still similar with the experimental results of long-tailed ones actuated under 70 Gs fields. With diversity in body length, the sperm-like nanorobots exhibited different sensitivity to the dynamic fields yet the locomotion directions kept basicly consistent. It can also be predicted that short-tailed nanorobots tend to move forward under precessing fields of an enhanced intensity.