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

HAL
open science

# Non-uniform magnetic fields for collective behavior of self-assembled magnetic pillars

Juan J Huaroto, Franco N Piñan Basualdo, Dionne Lisa Roos Ariëns, Sarthak Misra

► To cite this version:

Juan J Huaroto, Franco N Piñan Basualdo, Dionne Lisa Roos Ariëns, Sarthak Misra. Non-uniform magnetic fields for collective behavior of self-assembled magnetic pillars. Swarm Intelligence, 2024, 10.1007/s11721-024-00240-z. http://hal-04753522

HAL Id: hal-04753522

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

Submitted on 25 Oct 2024

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.

Public Domain
# Non-uniform magnetic fields for collective behavior of self-assembled magnetic pillars

Juan J. Huaroto¹ · Franco N. Piñan Basualdo¹ · Dionne Lisa Roos Ariëns¹ · Sarthak Misra¹,²

Received: 15 July 2023 / Accepted: 10 June 2024
© The Author(s) 2024

## Abstract

Programmable and self-assembled magnetic pillars are essential to expanding the application domain of magnetic microparticle collectives. Typically, the collective behavior of self-assembled magnetic pillars is carried out by generating uniform and time-varying magnetic fields. However, magnetic field-shaping capabilities employing non-uniform fields have not been explored for magnetic pillars. In this study, we generate non-uniform magnetic fields using a nine-coil electromagnetic system to achieve object manipulation, upstream/downstream locomotion, and independent actuation. We begin analyzing the static magnetic self-assembly of reduced iron microparticles and experimentally derive the average dimensions (height and diameter) of the resulting pillars. Subsequently, we delve into the collective dynamic response under non-uniform and time-varying magnetic fields, unveiling four distinct modalities. In order to demonstrate the versatility of our approach, we extend our study to the two-dimensional manipulation of a millimeter-sized glass bead using a precessing magnetic field describing a Lissajous curve. Moreover, we showcase the ability of magnetic pillars to adapt to confined and dynamic conditions within fluidic tubes. We finally present a noteworthy case where the nine-coil electromagnetic system independently actuates two clusters of magnetic pillars. Our study shows the potential of using non-uniform magnetic fields to actuate self-assembled magnetic pillars, enabling morphology reconfiguration capabilities, object manipulation, locomotion, and independent actuation.

**Keywords** Collective behavior · Magnetic pillars · Electromagnetic actuation · Magnetic gradients

## 1 Introduction

Collective behavior is a widespread phenomenon observed in nature across all scales (Whitesides & Grzybowski, 2002; Vicsek & Zafeiris, 2012), where groups of individual agents exhibit different coordinated/collective actions (Okubo, 1986). This phenomenon has been studied extensively in fields such as biology (Deneubourg et al., 1986) and physics (Mori, 1965). The existence of coordinated behavior can lead to emergent properties not present

<footnote>Extended author information available on the last page of the article</footnote>
in the individual agents themselves (Liljeström et al., 2019), which has inspired the devel-
opment of multiple robotic systems (Rubenstein et al., 2012; Werfel et al., 2014). These
collective properties can be harnessed to accomplish complex tasks and adapt to dynamic
environments robustly (Gardi et al., 2022).

In microrobotics, collective systems are advantageous over individual micro-agents since the latter are limited to performing specific tasks (Elgeti et al., 2015). Collective systems can be reorganized in situ and accomplish complex actions such as manipulation and cargo transport (Ahmed et al., 2021; Yang et al., 2022; Xu & Xu, 2022). Thus, understanding the interaction and self-organization of micro-agents is fundamental to developing novel collaborative applications. Furthermore, contactless actuation methods driven by magnetic fields (Yang et al., 2021; Xu et al., 2023), acoustic waves (Zhou et al., 2016; Liu et al., 2022), and light (Hu et al., 2022) can be potentially implemented under in vitro and ex vivo conditions, which enables their use in biomedical scenarios.

Magnetic actuation is a promising method for controlling collective systems due to its transparency and biocompatibility (Sitti & Wiersma, 2020). Magnetic collective systems can perform multimodal locomotion (Xie et al., 2019), morphology reconfiguration (Dong & Sitti, 2020), object/droplet manipulation (Kim et al., 2015; Sun et al., 2021), and cargo transportation (Akter et al., 2022). Inspired from cilia structures in nature, recent studies have proposed magnetic micropillars for generating programmable metachronal waves (Dong et al., 2020), enhanced object manipulation (Demirörs et al., 2021), and liquid/solid transport (Miao et al., 2022; Sohn et al., 2023). Magnetic micropillars are primarily fabricated using replica molding (Ren et al., 2022; Cui et al., 2023) and self-assembly of magnetic soft composites (i.e., a blend of silicone and magnetic microparticles) (Miao et al., 2022; Demirörs et al., 2021). Mobile permanent magnets or static arrangements of electromagnetic coils have been used to actuate magnetic micropillar arrays (Ni & Wang, 2023). Although the aforementioned fabrication methods permitted the engineering of the magnetization direction of micropillars, the pillars were attached to a substrate, restricting the morphology reconfiguration of pillars (Sohn et al., 2023). Finding a strategy to fabricate programmable pillars by magnetic self-assembly techniques is a current challenge in the development process (Ni & Wang, 2023).

Self-assembly is the spontaneous organization of individual agents into a stable con-
figuration. In particular, magnetic self-assembly has been applied to control the collective
behavior of magnetic droplets (Timonen et al., 2013; Wang et al., 2018), soft ferromagnetic
beads (Collard et al., 2022), and generate programmable pillars (Xu & Xu, 2022). Self-
assembled magnetic pillars use uniform magnetic fields for actuation, demonstrating effi-
cient deployment on biological tissue, navigation through uneven surfaces, and cargo trans-
portation (Xu & Xu, 2022; Xu et al., 2023). However, the homogeneous characteristics
of the magnetic field can limit the field-shaping capabilities and generation of gradients.
Prior work in electromagnetic actuation exploits non-uniform magnetic fields generated
by arrangements of multiple electromagnets for independent manipulation of microparti-
cles (Diller et al., 2013; Ongaro et al., 2018b; Piñan Basualdo & Misra, 2023). Therefore,
expanding the application domain of non-uniform magnetic fields to magnetic self-assem-
bled magnetic pillars has the potential to achieve morphology reconfiguration and inde-
pendent actuation of self-assembled magnetic pillar clusters.

To generate programmable and fabrication-free pillars based on the self-assembly of
magnetic microparticles, we exploit the field-shaping capabilities of a nine-coil electro-
magnetic system (Ongaro et al., 2018b). This electromagnetic system can generate time-
varying non-uniform magnetic fields to explore novel actuation modalities, which are used
for object manipulation and locomotion within fluidic tubes. In order to show the potential
of non-uniform magnetic fields, we exploit the spatially selective actuation of the elec-
tromagnetic system to perform independent actuation of magnetic pillar collectives. We
start analyzing the static magnetic self-assembly of the microparticles by experimentally
obtaining the average dimensions of the magnetic pillars. Additionally, we study the pil-
lars' collective dynamic response to non-uniform and time-varying magnetic fields, iden-
tifying four primary modalities. Furthermore, we demonstrate that a precessing magnetic
field describing a Lissajous curve (i.e., parametric curve composed of two perpendicular
sinusoidal signals with different phases and frequencies) can be utilized for the two-dimen-
sional manipulation of a millimeter-sized glass bead. We study the behavior of the pillars in
confined and dynamic environments by actuating them within fluidic tubes. We showcase
the capability of the nine-coil electromagnetic system by performing independent actuation
of two clusters of magnetic pillars. Our results show that non-uniform magnetic fields are a
promising approach for actuating self-assembled magnetic pillars, enabling reconfiguration
capabilities, object manipulation, locomotion, and independent actuation.

## 2 Results and discussion

The generation of magnetic self-assembly and collective behavior is carried out in a nine-coil electromagnetic system (Fig. 1A). Microparticles of reduced iron are placed into a polymethyl methacrylate (PMMA) box (base dimensions of 22 × 22 mm²) along with pure water (i.e., Milli-Q water) and surfactant (details in Sect. 4). The PMMA box containing the microparticles and suspension is located within the electromagnetic system workspace. A constant current is used to power the coil #9 as it generates the magnetic self-assembly of microparticles (i.e., magnetic pillars). Thereon, the coils #1–8 generate non-uniform time-varying magnetic fields and gradients, which trigger the collective behavior of magnetic pillars (Fig. 1B). The magnetization curve of reduced iron in Fig. 1C shows negligible remanence (magnetically soft material) and a linear behavior if the applied magnetic field is below 250 mT. The morphology of the microparticles is acquired from scanning electron microscopy, which shows particles with non-uniform diameters (1–5 µm) (Fig. 1D). In order to study the collective behavior of magnetic pillars, we divide our analysis into five sections: (1) static self-assembly of microparticles of reduced iron, (2) collective dynamic behavior of magnetic pillars, (3) object manipulation, (4) upstream and downstream locomotion, and (5) independent actuation of two clusters of magnetic pillars.

## 2.1 Static magnetic self-assembly

We begin the study by analyzing the self-assembly of microparticles of reduced iron to
create magnetic pillars. The microparticles are magnetized by an external magnetic field
generated by coil #9, as depicted in Fig. 2A. Using the current-to-field map of the sys-
tem, the magnetic field generated by powering coil #9 with 1 A is shown for the XY-/-
YZ-planes, with black arrows representing the magnetic field direction (Fig. 2B, C). In
order to assess the uniformity of the magnetic field, we compute the curvature of the
magnetic field at the center of the workspace. Unlike previous studies, which use per-
manent magnets at varying distances to modify the curvature (Xu & Xu, 2022; Ni &
Wang, 2023), we investigate the capabilities of an iron-cored electromagnet to create
magnetic pillars while maintaining a constant distance from the workspace (25 mm).
Considering the magnetic field is cylindrical symmetric along the coil axis, the
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Self-assembled magnetic pillars generated from microparticles of reduced iron. **A** Nine-coil electromagnetic actuation system. **B** Generation of self-assembled magnetic pillars and collective behavior. The states of the microparticles of reduced iron are indicated as follows: ①, base of the workspace surrounded the electromagnetic actuation system; ②, magnetic microparticles of reduced iron located within the workspace along with Milli-Q water and surfactant; ③, the coil #9 is activated with a constant current to generate the self-assembled magnetic pillars; ④, the coils #1–8 are activated to generate non-uniform and time-varying magnetic fields, which produce the collective behavior of magnetic pillars. **C** Magnetization curve of reduced iron. **D** Scanning electron microscopy (SEM) of microparticles of reduced iron</figcaption></figure>

curvature along the X and Y-axis exhibit identical values for a static magnetic field. We compute the curvature along the Y-axis ($c_Y = -\frac{d^2|\mathbf{B}(0,0,0)|}{dY^2}$) as a function of the magnetic field ($\mathbf{B}(0,0,0)$). Our results demonstrate that the curvature ($c_Y$) varies linearly with the current of coil #9 and the magnitude of $\mathbf{B}(0,0,0)$ (Fig. 2D).

By varying the current of coil #9 ($I_9 = 2–5$ A) and mass of microparticles of reduced iron ($m_{ri} = 10–40$ mg), we estimate the magnetic field at the center of the workspace ($\mathbf{B}(0,0,0)$) and compute the average height ($\bar{h}$) and base diameter ($\bar{d}$) of magnetic pillars (Fig. 2E). The cameras located at the side and above the workspace permit the acquisition of images, which are subsequently processed to compute $\bar{h}$ and $\bar{d}$, respectively (S1 Text, Supplementary information). Figure 2F, G shows that $\bar{h}$ increases according to the magnetic field while $\bar{d}$ remains approximately constant. On the other hand, $\bar{h}$ and $\bar{d}$ show an increment according to the mass of reduced iron microparticles within the workspace. It is worth noting that the increment of current through coil #9 is realized abruptly, and the coil is powered off ($I_9 = 0$ A) between measurements. This way, we reduce the influence between adjacent experiments.
<figure><img src="image_3.png" /><figcaption>Fig. 2 Characterization of static magnetic self-assembly of microparticles of reduced iron. A Workspace configuration with respect to coil #9. Magnetic field maps in the B XY and C YZ planes (powering coil #9 with 1 A). D Confining curvature of the magnetic field ($c_Y$) as a function of the magnetic field at the center of the workspace (B(0, 0, 0)). E Illustration of a magnetic pillar and average geometric parameters ($\bar{h}$ and $\bar{d}$). Experimental characterization of F the average height ($\bar{h}$) and G base diameter ($\bar{d}$) varying B(0, 0, 0) and mass of microparticles of reduced iron</figcaption></figure>

## 2.2 Dynamic collective behavior

The study of magnetic self-assembly of microparticles of reduced iron provides valuable insights into the formation of magnetic pillars under non-uniform and static magnetic fields. In this section, we focus on exploring the collective behavior of these magnetic pillars when exposed to non-uniform and time-varying magnetic fields. To this end, we represent all magnetic fields at the center of the workspace ($\mathbf{p}_0 = [0, 0, 0]^T$) (Fig. 3A). The magnetic pillars are generated from the static magnetic field created by coil #9 ($\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)$) (Fig. 3A) while coils #1–8 are utilized to create time ($t \in \mathbb{R}^+$) varying magnetic fields ($\tilde{\mathbf{B}}(\mathbf{p}_0, t)$). By considering the spatial configuration of unitary vectors ($\hat{\mu}$ and $\hat{\omega}$), which are perpendicular to the base of the workspace and the oscillation plane, respectively. Prior research has explored precessing and oscillating fields to perform multimodal locomotion of self-assembled magnetic pillars (Yigit et al., 2019; Xu & Xu, 2022; Xu et al., 2023). Drawing inspiration from the previous body of literature, we define two distinct cases for the actuation modalities (Fig. 3A). The resulting magnetic field at the center of the workspace ($\mathbf{B}(\mathbf{p}_0)$) for cases 1 ($\hat{\omega} \parallel \hat{\mu}$) and 2 ($\hat{\omega} \perp \hat{\mu}$) are mathematically defined as follows:

$$ \left\{ \begin{array}{l} \text{Case 1: } \mathbf{B}(\mathbf{p}_0) = \tilde{\mathbf{B}}(\mathbf{p}_0, I_9) + \tilde{\mathbf{B}}_X(\mathbf{p}_0, t) + \tilde{\mathbf{B}}_Y(\mathbf{p}_0, t) \\ \text{Case 2: } \mathbf{B}(\mathbf{p}_0) = \tilde{\mathbf{B}}(\mathbf{p}_0, I_9) + \tilde{\mathbf{B}}_X(\mathbf{p}_0, t) + \tilde{\mathbf{B}}_Y(\mathbf{p}_0, t) + \tilde{\mathbf{B}}_Z(\mathbf{p}_0, t), \end{array} \right. \quad (1) $$

where the magnetic fields ($\tilde{\mathbf{B}}_X(\mathbf{p}_0, t)$, $\tilde{\mathbf{B}}_Y(\mathbf{p}_0, t)$, and $\tilde{\mathbf{B}}_Z(\mathbf{p}_0, t)$) are the components of $\tilde{\mathbf{B}}(\mathbf{p}_0, t)$, and are represented by sinusoidal functions with amplitudes ($A_X$, $A_Y$, $A_Z$), frequencies ($f_X$, $f_Y$, $f_Z$), and phases ($\varphi_X$, $\varphi_Y$, $\varphi_Z$), respectively. According to the aforementioned cases and the observed collective behavior generated by the magnetic pillars, we propose the following actuation modalities: convergence, stretching, creep, and homogenization. The actuation modalities are experimentally tested using 30 mg of microparticles of
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Collective behavior of magnetic pillars. **A** Two cases of actuation modalities according to the spatial configuration of unitary vectors ($\hat{\mu}$ and $\hat{\omega}$) perpendicular to the base of the workspace and oscillation plane, respectively. The point ($\mathbf{p}_0 = [0, 0, 0]^T$) is defined to illustrate the components of the magnetic field ($\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)$ and $\tilde{\mathbf{B}}(\mathbf{p}_0, t)$). **B** Four actuation modalities are proposed to generate different collective behavior of magnetic pillars. The non-uniformity of the magnetic field is analyzed using directional curvatures ($\kappa = [c_X, c_Y]$) at $\mathbf{p}_0$. In addition, the magnetic field ($\mathbf{B}(X, Y, 0)$), gradients ($\nabla|\mathbf{B}(X, Y, 0)|$), and an image frame of the collective behavior of magnetic pillars are provided. The red dots indicate the center of the workspace ($\mathbf{p}_0$). Please refer to the accompanying video (S1 Video, Supplementary information). **C** Illustration of the external ($\mathbf{F}(\mathbf{p}_i)$, $\mathbf{F}(\mathbf{p}_{i-1})$, $\mathbf{F}(\mathbf{p}_i, I_9)$, $\mathbf{F}(\mathbf{p}_{i-1}, I_9)$), internal forces ($\mathbf{p}_{i-1} \cdot \mathbf{F}_{\mathbf{p}_i}$ and $\mathbf{p}_i \cdot \mathbf{F}_{\mathbf{p}_{i-1}}$), and torques ($\boldsymbol{\tau}(\mathbf{p}_i)$ and $\boldsymbol{\tau}(\mathbf{p}_{i-1})$) applied to adjacent magnetic pillars (Color figure online)</figcaption></figure>

reduced iron and using the parameters presented in Table 1 (Fig. 3B). For providing a metric on the non-uniformity of the magnetic field, each modality is analyzed through the directional curvatures of the magnetic field at the center of the workspace
<table><caption>**Table 1** Actuation modalities for self-assembled magnetic pillars</caption><thead><tr><th>Modality</th><th><i>l</i><sub>9</sub> [A]</th><th><i>A</i><sub>X</sub> [mT]</th><th><i>A</i><sub>Y</sub> [mT]</th><th><i>A</i><sub>Z</sub> [mT]</th><th><i>f</i><sub>X</sub> [Hz]</th><th><i>f</i><sub>Y</sub> [Hz]</th><th><i>f</i><sub>Z</sub> [Hz]</th><th><i>φ</i><sub>X</sub> [rad]</th><th><i>φ</i><sub>Y</sub> [rad]</th><th><i>φ</i><sub>Z</sub> [rad]</th></tr></thead><tbody><tr><td>Convergence</td><td>2</td><td>10</td><td>10</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Stretching</td><td>2</td><td>0</td><td>10</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Creep</td><td>2</td><td>0</td><td>8</td><td>8</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>π/2</td></tr><tr><td>Homogenization</td><td>2</td><td>0</td><td>0</td><td>8</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>
($\kappa = [c_X = -\frac{d^2|\mathbf{B}_{(0,0,0)}|}{dX^2}, c_Y = -\frac{d^2|\mathbf{B}_{(0,0,0)}|}{dY^2}]$). Moreover, the magnetic field ($\mathbf{B}(X, Y, 0)$), gradient ($\nabla|\mathbf{B}(X, Y, 0)|$), and image frame are shown at a specific percentage of the oscillation period (yellow points) (Fig. 3B).

Prior to discussing each actuation modality, we provide an illustration that depicts all the forces and torques acting on two adjacent pillars (identified by the subscripts *i* and *i* − 1) exposed to a non-uniform magnetic field. We assume that each magnetic pillar is composed of ($N \in \mathbb{N}$) microparticles of reduced iron. By powering the electromagnetic coil #9 with a constant current, we have observed the generation of magnetic pillars within the workspace. Each pillar behaves as a magnetic dipole since the microparticles of reduced iron are magnetized in response to the external magnetic field. The magnetic gradients generated by powering the coil #9 create forces ($\mathbf{F}(\mathbf{p}_i, I_0)$ and $\mathbf{F}(\mathbf{p}_{i-1}, I_0)$) that are oriented toward the base of the workspace to maintain the shape and structural integrity of the magnetic pillars. Furthermore, each pillar interacts with its adjacent pillar through a repulsive force ($\mathbf{p}_{i-1} \cdot \mathbf{p}_i = -\mathbf{p}_i \cdot \mathbf{p}_{i-1}$), which keeps the pillars separated by a distance ($b \in \mathbb{R}^+$) (Fig. 3C). Using coils #1–8, we generate non-uniform and time-varying magnetic fields and gradients. These magnetic fields and gradients exert torques ($\boldsymbol{\tau}(\mathbf{p}_i)$) and ($\boldsymbol{\tau}(\mathbf{p}_{i-1})$) as well as forces ($\mathbf{F}(\mathbf{p}_i)$) and ($\mathbf{F}(\mathbf{p}_{i-1})$) on the magnetic pillars, making them susceptible to tilting and slipping, respectively (Fig. 3C). It is important to note that, for the purpose of analysis, we are neglecting friction and viscous forces.

By understanding the influence of the magnetic field and gradients, we can gain further insights into the collective behavior of magnetic pillars under non-uniform and time-varying magnetic fields. We begin by discussing the actuation modalities corresponding to case 1. The magnetic field oscillates in a parallel plane with respect to the workspace. Thus, the forces ($\mathbf{F}(\mathbf{p}_i, I_0)$) remain constant during a complete period (Fig. 3C). Moreover, the torques ($\boldsymbol{\tau}(\mathbf{p}_i)$) and forces ($\mathbf{F}(\mathbf{p}_i)$) determine the motion of the magnetic pillar based on the time-varying magnetic fields and gradients, respectively. In particular, the convergence modality exhibits a rotating magnetic field, causing the curvatures along the X and Y-axes to oscillate reciprocally (Fig. 3B1). The magnetic pillars align themselves according to the magnetic field (black arrows), while the gradients (red arrows) generate forces, which tend to gather the magnetic pillars toward the center of the workspace, forming a circumference (S1 Video: 00:07–00:34 s, Supplementary information). The circularity of this circumference can be controlled by adjusting the values of the magnetic field amplitudes ($A_X$ and $A_Y$). On the other hand, the stretching modality can be considered a sub-case of the convergence modality, where the magnetic field along the X-axis is set to zero ($A_X = 0$ mT). As a consequence, the collective of magnetic pillars tends to stretch along the X-axis over time (S1 Video: 00:35–01:01 s, Supplementary information). Typically, the collective can be stretched along different axes by adjusting the values of $A_X$ and $A_Y$.

For actuation modalities corresponding to case 2, the effective magnetic field along the Z-axis varies according to the amplitude of the oscillating magnetic field along the Z-axis ($A_Z$). Therefore, we impose the following ratio to maintain the magnetic self-assembly of microparticles of reduced iron: $\frac{A_Z}{|\tilde{\mathbf{B}}(\mathbf{p}_0, I_0)|} < 1$. In our experiments, we utilize a maximum ratio ($\approx 0.34$) such that the magnitude of the static magnetic field predominates over the oscillating one. Since the effective magnetic field along the Z-axis oscillates over time. The curvatures along the X and Y-axes reach a minimum value when $\tilde{\mathbf{B}}(\mathbf{p}_0, I_0)$ and $\tilde{\mathbf{B}}(\mathbf{p}_0, t)$ have opposite directions. In particular, the creep modality uses a rotating magnetic field with components along the Y and Z-axes ($A_Y$ and $A_Z$) (Fig. 3B3). This way, the magnetic pillars are displaced in response to the resultant magnetic field. (S1 Video, Supplementary
information: 01:02–01:26 s). The predefined maximum ratio (≈ 0.34) prevents the magnetic pillars from rotating as their base remains in contact with the base of the workspace. Furthermore, the creep orientation can be modified by introducing a component along the X-axis ($f_X = f_Y$ and $\varphi_X = \varphi_Y$). This actuation modality is subsequently used for experiments involving the motion of magnetic pillars in fluidic tubes. The homogenization modality is a particular case of the creep modality, generating only an oscillating magnetic field along the Z-axis. The simulation and experimental results show that the magnetic pillars tend to move toward the boundaries and the center of the workspace (S1 Video, Supplementary information: 01:27–01:56 s s). This collective behavior leads to a reorganization (i.e., homogenization) of magnetic pillars within the workspace. We utilize the homogenization modality to reorganize the magnetic pillars between experiments.

## 2.3 Object manipulation

The dynamic collective behavior of magnetic pillars, achieved through the self-assembly of iron microparticles, presents a concept primarily explored for the assembly/disassembly of various spherical particles (diameter ≤ 1 mm and mass < 1 mg) (Xu & Xu, 2022). However, the challenge remains in achieving two-dimensional and dexterous manipulation of millimeter-sized objects while simultaneously preserving the positional integrity of the collective of magnetic pillars. In this study, we conduct the manipulation of a 2 mm glass bead (mass of 9.3 mg) using magnetic pillars within a workspace containing Milli-Q water and surfactant. The microparticles of reduced iron (20 mg mass) and the bead are positioned within the workspace while powering coil #9 ($\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)$) to generate the magnetic pillars. Next, we employ a time-varying magnetic field ($\tilde{\mathbf{B}}(\mathbf{p}_0, t)$) utilizing coils #1–8 to orient the magnetic pillars. The resulting precessing magnetic field is illustrated in Fig. 4A and mathematically represented as follows:

$$ \mathbf{B}(\mathbf{p}_0) = |\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)|\hat{\mathbf{k}} + 10 \sin(2\pi f_X t + \varphi_X)\hat{\mathbf{i}} + 10 \sin(2\pi f_Y t + \varphi_Y)\hat{\mathbf{j}} \text{[mT]}, \quad (2) $$

where $I_9 = 2$ A, $|\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)| = 21.1$ mT, and the unitary vectors ($\hat{\mathbf{i}}, \hat{\mathbf{j}}$ and $\hat{\mathbf{k}}$) corresponds to the X-/Y-/Z-axis, respectively. The precessing magnetic field ($\tilde{\mathbf{B}}(\mathbf{p}_0, t) + \mathbf{B}(\mathbf{p}_0, I_9)$) describes a Lissajous curve depending on the parameters of frequency ($f_X$ and $f_Y$) and phase ($\varphi_X$ and $\varphi_Y$) along the X-/Y-axis. We arrange these parameters as components of a matrix ($M = \begin{bmatrix} f_X & \varphi_X \\ f_Y & \varphi_Y \end{bmatrix}$), which is used to encode the following fundamental motion directions: forward, backward, left, and right (Fig. 4A). Figure 4B depicts the Lissajous curve used to move the bead toward the left direction ($M = \begin{bmatrix} 2f_X & 0 \\ f_X & 0 \end{bmatrix}$). For the analysis, we divide the Lissajous curve into four segments (Fig. 4B). We experimentally demonstrate that the Lissajous curve permits the relative displacement of the glass bead with respect to the magnetic pillars collective (details in S1 Text, Supplementary information). The intervals (④ → ① and ② → ③) allow for the effective stroke of the bead. Previous research using precessing magnetic fields demonstrates that the Lissajous curve can generate non-reciprocal motion of soft magnetic microparticles (Collard et al., 2020). Based on the experimental results depicted in Figure S2 (S1 Text, Supplementary information), we demonstrate that the trajectory described by the glass bead can be decomposed into translation
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Manipulation of a passive object (2 mm glass bead). **A** The components of the magnetic field: Constant magnetic field generated by the coil #9 ($\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)$), and oscillating magnetic field ($\tilde{\mathbf{B}}(\mathbf{p}_0, t)$) with amplitudes $A_X$ and $A_Y$ along the X and Y-axis, respectively. The matrix ($M \in \mathbb{R}^{2 \times 2}$), which encompasses the frequencies ($f_X$, $f_Y$) and phases ($\varphi_X$, $\varphi_Y$), is used to encode the direction of motion of the bead. **B** Illustration of the periodic motion performed by the magnetic pillars to move the 2 mm glass bead towards the left. The precessing magnetic field ($\tilde{\mathbf{B}}(\mathbf{p}_0, t) + \tilde{\mathbf{B}}(\mathbf{p}_0, I_9)$) describes a Lissajous curve according to the corresponding matrix $M$. The effective stroke is achieved within the intervals (② → ③) and (④ → ①), whereas the recovery stroke is performed within (① → ②) and (③ → ④). **C** Manipulation of a 2 mm glass bead (encircled in yellow). The images show the trajectory (in red) of the glass bead as it describes the alphabet I: “S,” II: “R,” and III: “L,” respectively. The insets show the initial and final time stamps. *Please refer to the accompanying video* (S2 Video, Supplementary information) (Color figure online)</figcaption></figure>

and a cyclic motion following the input Lissajous curve. During an actuation cycle, the bead follows the pillars' motion. However, at the end of the cycle, the bead does not return to its initial position, exhibiting a net translation. Since this displacement is always in the same direction relative to the input Lissajous curve, it suffices to change the input signal orientation to select the advancement direction. In the object manipulation experiments, the four motion directions are tested to replicate trajectories that depict the alphabet “S,” “R,” and “L” (representing the initials of Surgical Robotics Laboratory) (Fig. 4C).

## 2.4 Upstream and downstream locomotion

The capability of magnetic pillars to withstand and move within hollow cavities transporting fluid holds importance for studying potential biomedical applications under ex vivo conditions. We start testing the capability of the magnetic pillars to withstand fluid flow. To this end,
we place 10 mg of microparticles of reduced iron inside a glass tube (8 mm inner diameter), connected by one of its ends to a syringe pump. The magnetic pillars are formed by circulating a constant current through coil #9 ($I_9$ = 2 A). Milli-Q water is automatically pumped at the maximum flow rate capacity (100 ml/min) to evaluate the displacement of the magnetic pillars (S3 Video, Supplementary information: 00:07–00:40 s). Our results show that the center of mass of magnetic pillars moves approximately 1.4 mm (for a duration of 15 s). Despite the susceptibility of the magnetic pillars to be carried away by the fluid flow, we observe that the collective of pillars withstands the drag forces acting upon it. Furthermore, experiments with 80 ml/min and 60 ml/min show a displacement of less than 0.5 mm for a duration of 15 s.

Following our study, we test the locomotion of magnetic pillars upstream and downstream using the creep modality ($f_Y = f_Z = 1.5$ Hz and $A_Y = A_Z = 8$ mT) within a fluidic tube (Fig. 5A). The experiments are carried out for a duration of 50 s (Fig. 5B, C). The center of mass of the magnetic pillars is tracked while performing periodic motion according to the actuation frequency. The black lines depicted in Fig. 5D, E represent the displacement in still water. As the flow rate increases (20 ml/min and 40 ml/min), we observe that the upstream

<figure><img src="image_6.png" /><figcaption>**Fig. 5** Motion of magnetic pillars within a tube transporting Milli-Q water. A Magnetic pillars formed from 10 mg of microparticles of reduced iron exposed to a constant magnetic field ($|\vec{B}(p_0, t_0)|$ = 21.1 mT, with $t_0$ = 2 A). The creep modality, which employs a rotating magnetic field ($\vec{B}(p_0, t)$), is utilized to move the pillars upstream (A1) and downstream (A2) the fluid flow. The yellow and pink dots represent the starting and final point, respectively (for a duration of 50 s). Top views of the starting and final frames of magnetic pillars displacing B upstream and C downstream using flow rates of 20 ml/min and 40 ml/min. Displacement of the center of mass of magnetic pillars in D upstream and E downstream conditions. The frames corresponding to the times $t = 0$ s and 50 s show the magnetic pillar states when the electromagnetic system is deactivated and activated, respectively. Please refer to the accompanying video (S3 Video, Supplementary information) (Color figure online)</figcaption></figure>
displacement of the center of mass decreases by 47.3% and 64.0%, respectively (Fig. 5D). Conversely, for downstream experiments, the center of mass displacement increases with the flow rate by 13.3% and 29.0%, respectively (Fig. 5E). Our results demonstrate that the magnetic pillars generated using 10 mg of reduced iron can navigate upstream and downstream, utilizing a maximum flow rate of 40 ml/min. Furthermore, the magnetic pillars can structurally withstand a fluid flow of 100 ml/min. It is worth noting that the saturation tendency on the displacement observed in Fig. 5D, E is due to the gradient generated by coil #9, which attracts the magnetic pillars toward the center of the workspace. In order to reach a more significant displacement, the electromagnetic coil arrangement can be assembled at the end effector of robotic platforms to constantly relocate the magnetic pillars collective (Yang et al., 2019; Wang et al., 2021).

## 2.5 Independent actuation of two clusters of magnetic pillars

The previous experiments are carried out using non-uniform and time-varying magnetic fields, which are the result of imposing a magnetic field at the center of the workspace ($\mathbf{p}_0 = [0, 0, 0]^T$). Here, we leverage the capabilities of nine iron-cored coils to perform the independent motion of two clusters of magnetic pillars. To this end, we impose two rotating magnetic fields ($\tilde{\mathbf{B}}(\mathbf{p}_1, t)$ and $\tilde{\mathbf{B}}(\mathbf{p}_2, t)$) defined at two points ($\mathbf{p}_1 = [-8, 0, 0]^T$ mm and $\mathbf{p}_2 = [8, 0, 0]^T$ mm) within the workspace (Fig. 6A). The magnetic fields exhibit opposite directions, characterized by amplitudes of 8 mT and frequencies of 1 Hz. The resultant magnetic fields ($\mathbf{B}(\mathbf{p}_1)$ and $\mathbf{B}(\mathbf{p}_2)$) are defined as follows:

<figure><img src="image_7.png" /><figcaption>**Fig. 6** A Illustration of the independent actuation of two clusters of magnetic pillars. The magnetic pillars are formed from 10 mg of microparticles of reduced iron exposed to a constant magnetic field ($|\tilde{\mathbf{B}}(\mathbf{p}_0, I_9)| = 21.1$ mT, with $I_9 = 2$ A). Two points ($\mathbf{p}_1$ and $\mathbf{p}_2$) within the workspace are used to generate two different oscillating magnetic fields ($\tilde{\mathbf{B}}(\mathbf{p}_1, t)$ and $\tilde{\mathbf{B}}(\mathbf{p}_2, t)$). **B** and **C** show the magnetic field and gradient within the workspace at $t = 0.25$ s (25% of one actuation period), respectively. The pink and purple dots represent the points $\mathbf{p}_1$ and $\mathbf{p}_2$, respectively. **D** Time frames show the independent motion of two clusters of magnetic pillars. Please refer to the accompanying video (S4 Video, Supplementary information) (Color figure online)</figcaption></figure>
$$ \left\{ \begin{array}{l} \mathbf{B}(\mathbf{p}_1) = |\bar{\mathbf{B}}(\mathbf{p}_1, I_9)|\hat{\mathbf{k}} - 8 \sin(2\pi t)\hat{\mathbf{j}} + 8 \cos(2\pi t)\hat{\mathbf{k}} \, [\text{mT}] \\ \mathbf{B}(\mathbf{p}_2) = |\bar{\mathbf{B}}(\mathbf{p}_2, I_9)|\hat{\mathbf{k}} + 8 \sin(2\pi t)\hat{\mathbf{j}} + 8 \cos(2\pi t)\hat{\mathbf{k}} \, [\text{mT}], \end{array} \right. \qquad (3) $$

where $\hat{i}$, $\hat{j}$ and $\hat{k}$ are the unitary vectors along the X-/Y-/Z-axis, respectively. Moreover, $|\bar{\mathbf{B}}(\mathbf{p}_1, I_9)| \approx |\bar{\mathbf{B}}(\mathbf{p}_2, I_9)| = 19.3 \, \text{mT}$ (using $I_9 = 2 \, \text{A}$), since the magnetic field generated by the coil #9 is cylindrical symmetric along the coil axis (Fig. 2B, C). Using the current-to-field map of the system, we obtain the magnetic field and gradients, where the pink and purple points represent $\mathbf{p}_1$ and $\mathbf{p}_2$, respectively (Fig. 6B, C). For the experiments, two transparent tubes (6 mm inner diameter) are filled with Milli-Q water and 10 mg of micro-particles each. The tubes are located within the workspace using two plastic dome nuts (Fig. 6D). The microparticles of reduced iron are approximately located at the center of each tube. Thereon, the electromagnetic system is activated according to (3) as the magnetic pillars move toward opposed directions.

## 3 Conclusions

We introduce non-uniform magnetic fields to explore the collective behavior of self-assembled magnetic pillars formed from microparticles of reduced iron. Our study utilizes an electromagnetic system composed of nine iron-cored coils, which provide field-shaping and spatially selective actuation capabilities within the workspace. We experimentally determine the dimensions of self-assembled magnetic pillars and investigate four different actuation modalities. Through experimental validation, we demonstrate that a precessing magnetic field describing a Lissajous curve can be used for in-plane manipulation of a millimeter-sized glass bead. Furthermore, we examine the performance of magnetic pillars within a fluid-filled tube, providing insights into their ability to withstand and move under varying fluid flow conditions. Lastly, we leverage the potential of the nine-coil electromagnetic system to achieve the independent actuation of two clusters of magnetic pillars. Our results open up new avenues for further exploration and applications of self-assembled magnetic pillars exposed to non-uniform and time-varying magnetic fields.

## 4 Experimental section

*Instrumentation:* The mass of microparticles of reduced iron utilized in this study is measured on a high-precision scale (SECURA225D-1S, Sartorius Lab Instruments GmbH & Co. KG, Germany). The microparticles are visualized using a scanning electron microscope (JSM-7200F, JEOL, Japan) to reveal further morphological details. In addition, the magnetization curve of microparticles of reduced iron is obtained using a vibrating-sample magnetometer (GMW 3474-140, GMW Associates, USA). Two cameras (FLIR, Wilsonville, USA), positioned above and to the side of the workspace, are equipped with microscope objectives (Optem FUSION, Qioptiq Ltd., UK) to capture high-resolution images (2048 × 2048 pixels) with a pixel resolution of approximately 10 µm. For experiments involving fluid flow, a syringe pump (NE-4000, New Era Pump Systems, USA) generates constant flow rates of Milli-Q water within a transparent tube.

*Simulations and image processing:* The simulations of the magnetic field and gradients are obtained from a current-to-field map (Ongaro et al., 2018a) that is implemented
in MATLAB (version 2023a, MathWorks, USA) for visualization. The algorithms for tracking and computing geometric parameters (average height and diameter) of magnetic pillars are implemented in a custom script of MATLAB. *Please refer to the accompanying text for details* (S1 Text, Supplementary information).

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s11721-024-00240-z.

**Acknowledgements** We want to thank Mr. Nick Helthuis for his assistance with the scanning electron microscope and Mr. Yu-Hsiang Lin for acquiring the magnetization curve of microparticles of reduced iron.

**Author contributions** J.H. conceived and managed the project and took the lead in writing the manuscript with contributions from all authors. D.L.R.A. conducted experiments on the static magnetic self-assembly of microparticles of reduced iron. Additionally, J.H. performed and recorded the remaining experiments presented in the manuscript. J.H. also compiled the figures and videos. F.N.P.B. developed the graphical user interface to conduct the experiments. All authors thoroughly reviewed the manuscript and actively participated in discussions regarding the selection of experimental demonstrations. S.M. secured the funding for the project and provided supervision throughout its course.

**Funding** This work is supported by the European Research Council (ERC) under the European Union's Horizon 2020 Research and Innovation programme under Grant #866494 (project - MAESTRO).

**Data availability** The authors confirm that the data supporting the findings of this study are available within the article and its Supplementary information. Raw data that support the findings of this study are available from the corresponding author, upon reasonable request.

## Declarations

**Conflict of interest** The Authors declare that there is no conflict of interest.

**Ethics approval and consent to participate** Not applicable.

**Consent for publication** Not applicable.

**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/.

## References

Ahmed, D., Sukhov, A., Hauri, D., Rodrigue, D., Maranta, G., Harting, J., & Nelson, B. J. (2021). Bioinspired acousto-magnetic microswarm robots with upstream motility. *Nature Machine Intelligence, 3*(2), 116-124.

Akter, M., Keya, J. J., Kayano, K., Kabir, A. M., Inoue, D., Hess, H., Sada, K., Kuzuya, A., Asanuma, H., & Kakugo, A. (2022). Cooperative cargo transportation by a swarm of molecular machines. *Science Robotics, 765*, eabm0677.

Collard, Y., Grosjean, G., & Vandewalle, N. (2020). Magnetically powered metachronal waves induce locomotion in self-assemblies. *Communications Physics, 3*(1), 112.

Collard, Y., Piñan Basualdo, F. N., Bolopion, A., Gauthier, M., Lambert, P., & Vandewalle, N. (2022). Controlled transitions between metastable states of 2d magnetocapillary crystals. *Scientific Reports, 12*(1), 16027.
Cui, Z., Wang, Y., Zhang, S., Wang, T., & den Toonder, J. M. (2023). Miniaturized metachronal magnetic artificial cilia. *Proceedings of the National Academy of Sciences, 120*(35), e2304519120.

Demirörs, A. F., Aykut, S., Ganzekoorn, S., Meier, Y. A., Hardeman, R., de Graaf, J., Mathijssen, A. J., Poloni, E., Carpenter, J. A., Ünlü, C., & Zenhäusern, D. (2021). Amphibious transport of fluids and solids by soft magnetic carpets. *Advanced Science, 8*(21), 2102510.

Deneubourg, J. L., Aron, S., Goss, S. A., Pasteels, J. M., & Duerinck, G. (1986). Random behaviour, amplification processes and number of participants: How they contribute to the foraging properties of ants. *Physica D: Nonlinear Phenomena, 22*(1–3), 176–186.

Diller, E., Giltinan, J., Jena, P., & Sitti, M. (2013). Three dimensional independent control of multiple magnetic microrobots. In *Proceedings of the 2013 IEEE international conference on robotics and automation* (pp. 2576–2581). IEEE.

Dong, X., Lum, G. Z., Hu, W., Zhang, R., Ren, Z., Onck, P. R., & Sitti, M. (2020). Bioinspired cilia arrays with programmable nonreciprocal motion and metachronal coordination. *Science Advances, 6*(45), eabc9323.

Dong, X., & Sitti, M. (2020). Controlling two-dimensional collective formation and cooperative behavior of magnetic microrobot swarms. *The International Journal of Robotics Research, 39*(5), 617–638.

Elgeti, J., Winkler, R. G., & Gompper, G. (2015). Physics of microswimmers-single particle motion and collective behavior: A review. *Reports on Progress in Physics, 78*(5), 056601.

Gardi, G., Ceron, S., Wang, W., Petersen, K., & Sitti, M. (2022). Microrobot collectives with reconfigurable morphologies, behaviors, and functions. *Nature Communications, 13*(1), 2239.

Hu, D., Wang, Y., Liu, J., Mao, Y., Chang, X., & Zhu, Y. (2022). Light-driven sequential shape transformation of block copolymer particles through three-dimensional confined self-assembly. *Nanoscale, 14*(17), 6291–6298.

Kim, J. H., Kang, S. M., Lee, B. J., Ko, H., Bae, W. G., Suh, K. Y., Kwak, M. K., & Jeong, H. E. (2015). Remote manipulation of droplets on a flexible magnetically responsive film. *Scientific Reports, 5*(1), 17843.

Liljeström, V., Chen, C., Dommersnes, P., Fossum, J. O., & Gröschel, A. H. (2019). Active structuring of colloids through field-driven self-assembly. *Current Opinion in Colloid & Interface Science, 40*, 25–41.

Liu, W., Gao, H., Liu, K., Lei, D., Pei, K., & Hu, G. (2022). A review on particle assembly in standing wave acoustic field. *Journal of Nanoparticle Research, 24*(4), 81.

Miao, J., Zhang, T., Li, G., Shang, W., & Shen, Y. (2022). Magnetic artificial cilia carpets for transport, mixing, and directional diffusion. *Advanced Engineering Materials, 24*(7), 2101399.

Mori, H. (1965). Transport, collective motion, and Brownian motion. *Progress of Theoretical Physics, 33*(3), 423–455.

Ni, K., & Wang, Z. (2023). Recent progress on the development of magnetically-responsive micropillars: Actuation, fabrication, and applications. *Advanced Functional Materials, 33*(14), 2213350.

Okubo, A. (1986). Dynamical aspects of animal grouping: Swarms, schools, flocks, and herds. *Advances in Biophysics, 22*, 1–94.

Ongaro, F., Heunis, C. M., & Misra, S. (2018). Precise model-free spline-based approach for magnetic field mapping. *IEEE Magnetics Letters, 10*, 1–5.

Ongaro, F., Pane, S., Scheggi, S., et al. (2018). Design of an electromagnetic setup for independent three-dimensional control of pairs of identical and nonidentical microrobots. *IEEE Transactions on Robotics, 35*(1), 174–183.

Piñan Basualdo, F. N., & Misra, S. (2023). Collaborative magnetic agents for 3d microrobotic grasping. *Advanced Intelligent Systems, 5*, 2300365.

Ren, Z., Zhang, M., Song, S., et al. (2022). Soft-robotic ciliated epidermis for reconfigurable coordinated fluid manipulation. *Science Advances, 8*(34), eabq2345.

Rubenstein, M., Ahler, C., & Nagpal, R. (2012). Kilobot: A low cost scalable robot system for collective behaviors. In *Proceedings of IEEE international conference on robotics and automation* (pp. 3293–3298), IEEE.

Sitti, M., & Wiersma, D. S. (2020). Pros and cons: Magnetic versus optical microrobots. *Advanced Materials, 32*(20), 1906766.

Sohn, S., Lee, H., Kee, H., et al. (2023). Reprogrammable magnetically actuated self-assembled cilia array. *Advanced Intelligent Systems, 5*(1), 2200227.

Sun, M., Fan, X., Tian, C., et al. (2021). Swarming microdroplets to a dexterous micromanipulator. *Advanced Functional Materials, 31*(19), 2011193.

Timonen, J. V., Latikka, M., Leibler, L., et al. (2013). Switchable static and dynamic self-assembly of magnetic droplets on superhydrophobic surfaces. *Science, 341*(6143), 253–257.

Vicsek, T., & Zafeiris, A. (2012). Collective motion. *Physics Reports, 517*(3–4), 71–140.
Wang, L., Wang, J., Hao, J., et al. (2021). Guiding drug through interrupted bloodstream for potentiated thrombolysis by c-shaped magnetic actuation system in vivo. *Advanced Materials*, **33**(51), 2105351.

Wang, Q., Yang, L., Wang, B., et al. (2018). Collective behavior of reconfigurable magnetic droplets via dynamic self-assembly. *ACS Applied Materials & Interfaces*, **11**(1), 1630–1637.

Werfel, J., Petersen, K., & Nagpal, R. (2014). Designing collective behavior in a termite-inspired robot construction team. *Science*, **343**(6172), 754–758.

Whitesides, G. M., & Grzybowski, B. (2002). Self-assembly at all scales. *Science*, 295(5564), 2418–2421.

Xie, H., Sun, M., Fan, X., et al. (2019). Reconfigurable magnetic microrobot swarm: Multimode transformation, locomotion, and manipulation. *Science Robotics*, **4**(28), eaav8006.

Xu, Z., Wu, Z., Yuan, M., et al. (2023). Multiple cilia-like swarms enable efficient microrobot deployment and execution. *Cell Reports Physical Science*, **4**(3), 101329.

Xu, Z., & Xu, Q. (2022). Collective behaviors of magnetic microparticle swarms: From dexterous tentacles to reconfigurable carpets. *ACS Nano*, **16**(9), 13728–13739.

Yang, L., Du, X., Yu, E., et al. (2019). Deltamag: An electromagnetic manipulation system with parallel mobile coils. In *Proceedings of the 2019 International Conference on Robotics and Automation (ICRA)* (pp. 9814–9820). IEEE

Yang, L., Yu, J., Yang, S., et al. (2021). A survey on swarm microrobotics. *IEEE Transactions on Robotics*, **38**(3), 1531–1551.

Yang, S., Wang, Q., Jin, D., et al. (2022). Probing fast transformation of magnetic colloidal microswarms in complex fluids. *ACS Nano*, **16**(11), 19025–19037.

Yigit, B., Alapan, Y., & Sitti, M. (2019). Programmable collective behavior in dynamically self-assembled mobile microrobotic swarms. *Advanced Science*, **6**(6), 1801837.

Zhou, Q., Sariola, V., Latifi, K., et al. (2016). Controlling the motion of multiple objects on a chladni plate. *Nature Communications*, **7**(1), 12764.

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Authors and Affiliations

Juan J. Huaroto<fnref n="1" /> · Franco N. Piñan Basualdo<fnref n="1" /> · Dionne Lisa Roos Ariëns<fnref n="1" /> · Sarthak Misra<fnref n="1" />,<fnref n="2" />

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

Juan J. Huaroto
j.j.j.huarotosevilla@utwente.nl

Franco N. Piñan Basualdo
f.n.pinanbasualdo@utwente.nl

Dionne Lisa Roos Ariëns
d.l.r.ariens@student.utwente.nl

Sarthak Misra
s.misra@utwente.nl

<footnote n="1">¹ Surgical Robotics Laboratory, Department of Biomechanical Engineering, University of Twente, 7522 NB Enschede, The Netherlands</footnote>

<footnote n="2">² Surgical Robotics Laboratory, Department of Biomaterials and Biomedical Technology, University Medical Centre Groningen and University of Groningen, 9713 GZ Groningen, The Netherlands</footnote>