Manipulation of emergent vortices in swarms of magnetic rollers


Magnetic colloidal rollers

In our experiments we sediment ferromagnetic nickel microspheres (RNi = 69 μm) dispersed in a liquid at the bottom of a glass container. A uniform vertical alternating magnetic field (Bz = B0sin(2πfBt), where t is time, B0 is the field amplitude, and fB is the frequency) is used to energize the particles (see the Methods section). In a certain range of excitation field parameters29 particles spontaneously break the symmetry of the clockwise/counterclockwise rotations experienced by a magnetic particle in a uniaxial field and start to spin20,29 creating a directionally uncorrelated rolling motions of many particles. The particles steadily spin when the following condition is satisfied:


29. Here the Mathieu characteristic exponent ν is the function of the two parameters p = αr/(ωI) and q = μB0/(ω2I) where ω = 2πfB, η is a fluid kinematic viscosity and μ, m,




are correspondingly the magnetic moment, mass, moment of inertia and the rotational drag coefficient of a roller. In the state of a steady rotation the magnetic torque on the particle, particle,


, is balanced by its viscous torque,


, and characteristic viscous time τv is comparable with a time scale of the applied field frequency τf. The frequency of the alternating magnetic field controls the speed of the rolling motion and it is employed to tune the activity in the system. Magnetic colloidal rollers demonstrate strong propensity toward the onset of a large scale collective motion. A set of dynamic phases ranging from gas to intermittent flocks and emergence of a global vortex has been observed in this system in response to changes in the rollers’ activity29. The formation of vortices, see Fig. 1a, b and Supplementary Movie 1 and 2, is observed in our system in a relatively narrow region of the field excitation frequencies (about 10 band); however, a magnetic roller vortex is a robust reproducible entity. The chirality of the vortex rotation (clock- or counterclockwise) is random from experiment to experiment. The core of the vortex is characterized almost linear tangential velocity profile (see Fig. 1d). Remarkably, the emergence of the magnetic roller vortices does not rely on the presence of geometrical boundaries or confinement. A vortex was successfully observed in a flat bottom glass container with the surface area significantly larger than the area of the self-organized roller vortex so that the rollers forming the vortex do not interact with the container’s walls (see Supplementary Movie 3).

Fig. 1
Fig. 1

A magnetic roller vortex. a An overlay of five consecutive images showing a single roller vortex in a spherical potential well (radius of curvature 52 mm, rotating clockwise (CW). A particle grayscale intensity indicates time of the frame, light-colored particles being the initial frame and dark particles the final (see Supplementary Movie 1). Scale bar length is 2 mm. b Calculated velocity amplitude field. The white star symbol is the position of the vortex eye as determined by velocity and vorticity maps. Dashed circle around a vortex eye illustrates the core of the vortex as determined from the azimuthally averaged velocity profile of the vortex shown in d. Scale bar is 2 mm. c Potential energy U of the vortex eye in the spherical well has a harmonic dependence on distance from center r. U was calculated using Eq. 1 from scatter plots of the vortex eye positions (insets). Scale bars are 1. Error bars are the standard deviation of the measurements. As indicated by the order parameter ϕR in e the vortex spontaneously switched rotation direction during the course of the experiment. For scatter plots only frames with |ϕR| > 0.3 were chosen. d Azimuthally averaged vortex tangential velocity profile versus distance from the center for the vortex shown in b. Dashed line marks the boundary of the defined vortex core where the rollers speed reaches maximum. e Time evolution of the polar order parameter ϕR in a roller vortex demonstrating a spontaneous switch of the vortex chirality. f The pair correlation function for the rollers forming the vortex. g(r) has a single peak, implying the presence of a characteristic spacing between rollers. ae The magnetic field direction is aligned with the gravity and oscillates in- and out-of-plane (amplitude B0 = 5.78 mT, frequency fB = 40 Hz)

To quantify the vortex we calculate the polar order parameter,


. Here,


are the in-plane unit vectors in angular direction,


are the in-plane velocities over the grid points and N is the number of grid points. In the ideal global polar state the order parameter will reach unity in magnitude, but in the case of the roller vortex state in our system ϕR(t) is smaller since the vortex is finite (not global) and some of the particles are in the gas phase beyond the boundaries of the vortex. As the system is dynamic by nature the polar order parameter fluctuates around a mean value with a normal distribution of fluctuations (see Supplementary Note 1 and Supplementary Figure 1).

Chirality state of a roller-vortex

Rollers from the gas perpetually join and leave the vortex and the whole system is very dynamic. While the direction of the vortex rotation is randomly selected by the system from experiment to experiment, occasionally the vortex can spontaneously change its chirality state as demonstrated in Fig. 1e. In Fig. 2a we demonstrate chirality switching events as they manifest in the polar order parameter of a roller vortex formed in a soft harmonic gravitational confinement. As one can see from the graph the vortex chirality switching is stochastic and relatively rare event on the scale of a typical experiment. The spontaneous chirality switching usually proceeds through the intermittent formation of flocks and takes of the order of 10 s (also see Supplementary Movie 4). During this intermittent process the rollers again randomly selects new chirality state that also can result in instances when vortex falls back to its previous chiral state as shown in Fig. 2b. To get an additional insight into the statistics of the vortex chirality switching in our system we extracted probability distribution functions (Fig. 2c) for the vortex to have no switching events (P0) and have one successful chirality switching (P1). Both curves suggest that on average in our system a vortex chirality switching event happens on the scale of about 170 min. As we showed previously in simulations29 phase synchronization between rollers plays an essential role in the formation of the roller vortex phase and apparent intermittencies in the vortex behavior resulting in vortex chirality switching most probably stem from the temporary loss and subsequent recovery of the global phase synchrony of the rollers comprising the vortex. Temporary local loss of synchronization can be triggered in the experiments by the collisions with neighboring rollers that have pronounced shape imperfections resulting in abnormal rotational diffusion.

Fig. 2
Fig. 2

Chirality switching statistics of a roller vortex. a Chirality switching events illustrated by the polar order parameter of a roller vortex formed in a soft harmonic gravitational confinement. B0 = 5.78 mT, fB = 40 Hz. b Events of the flock intermittency not resulted in a roller vortex chirality switch. c Probability distribution functions for the roller vortex to have no chirality switching events (P0) or have one successful chirality change (P1). Error bars are the standard deviation of the measurements

Unconfined roller vortices

In contrast to its Quincke roller counterpart where nearly all particles are concentrated along the container boundary in the vortex phase13, the magnetic roller vortex core has almost a ‘solid’ structure. Its density is approximately uniform in the core (see Supplementary Movie 2). To quantify the structure of the vortex core the pair correlation function g(r) has been calculated for the rollers forming the vortex (see Fig. 1d). It exhibits a pronounced peak at about twice the particle diameter indicative of a short-range spacing between rollers forming the vortex.

Spontaneous formation of the magnetic roller vortices at a flat surface is a nontrivial collective phenomenon. To form and maintain a roller vortex a certain local number density of rollers has to be met. In a harmonic gravitational confinement29 this is automatically reached due to the herding of rollers by the confining potential. At a flat surface the situation is different, however the system is able to spontaneously form and maintain a local vortex. The number density of particles inside a roller vortex usually exceeds the average surface number density of the system (about 11 mm−2 inside the vortex versus 6 mm−2 overall in the system), see Supplementary Movie 3. Dynamic local densifications of rollers forming the vortices at a flat surface are driven by the rollers themselves and are intrinsic property of the magnetic roller system. The size of the vortex is selected by a dynamic self-induced densification and a range of different sizes can be realized in the system at the same experimental conditions, see Supplementary Note 2 and Supplementary Figure 2. It is possible to have vortices as small as 1.2 mm and as large as 3.1 mm in diameter. However, the large vortices are less stable and may on average fall apart faster to form smaller entities; on the other hand small vortices may evaporate to a gas. The most probable size in the studied system was about 2 mm. This size could be altered by the number density of rollers, however it is insensitive to the driving field frequency or amplitude manipulations. In a concave surface case the needed density for the vortex formation is maintained and stabilized by a soft confinement and much wide range of vortex sizes can be observed. Once formed, the roller vortex is a dynamic entity and can move around the surface for minutes before it disintegrates due to interactions with obstacles or other flocks.

The main mechanism behind dynamic densifications forming the vortex in our system is analogous in nature to Vicsek flocks37 but driven in our system by a fine interplay between flows (advection forces) generated by individual magnetic rollers and magnetic interactions between rollers. In particular, each magnetic roller has a complex time-averaged interaction profile: rotation of the sphere in the fluid (Re > 1 for rollers, inertia is important) creates attractive hydrodynamic interactions in the lateral direction (along the axis of rotation, perpendicular to the direction of the roller motion) and repulsive in the direction perpendicular to the axis of rotation (along the direction of the roller motion)38,39. As a result, rollers hydrodynamicaly attract neighbors laterally and repel them if they are along the rollers direction of motion. Both forces decay as 1/r3 with a distance39,40. Time averaged magnetic interactions, on the other hand, are attractive along the rollers direction of motion and repulsive in lateral direction. Corresponding magnetic time averaged forces decay as 1/r415. Thus, hydrodynamic and magnetic interactions between rollers keep the roller vortex from falling apart (long range attractions prevent rollers from departing the vortex), or collapsing to clusters (short range repulsions keep rollers from getting to close to form chains or clusters). However, collisions with obstacles (scatterers), or other flocks of rollers may disturb the balance and phase synchrony29 of the rollers forming the vortex and it may get transformed to flocks or disintegrate.

To demonstrate the pivotal role of the induced hydrodynamic interactions in the emergence of vortices in a magnetic roller ensemble we completely eliminated the hydrodynamic effects by performing similar experiments in air. This is possible since spontaneous rotation of magnetic rollers do not rely on the presence of a liquid (in contrast to Quincke rollers). No vortices or flocks have been observed in a full range of the field parameters and number densities, where steady rolling of particles was possible in air, also see Supplementary Movie 5.

Roller vortex and inert scatterers

Once formed, the roller vortex is not a static entity. The vortex core explores the bottom surface of the experimental container (see Supplementary Movie 3), and it can be considered itself as a quasi-particle. To investigate in detail the effects of inert scatterers on dynamics of a single magnetic roller vortex we isolated the vortex in a soft harmonic gravitational well realized by a container with a concave bottom (we use a glass lens with a radius of curvature Rcurv = 52 mm). This way the vortex gets gravitationally localized around the center of the lens that prevents it from escaping the field of view or moving too close to the walls of the container that may disrupt or modify the vortex under investigation.

The roller vortex explores the vicinity of the concave surface center. We determined the position of the vortex eye (where velocity is close to zero, see Fig. 1b) for each time step and treat it as a quasi-particle exploring the potential landscape. The confining potential imposed by the concave surface on the vortex could be estimated using inverted Boltzmann equation41,



Here n(r) is a radial displacement histogram, and A is a constant dependent on total number of measurements. Our observations (see Fig. 1c) reveal a harmonic displacement dependance of the potential


, where kp is the effective potential stiffness felt by the active vortex on the concave surface. The extracted potential stiffness shows only very weak dependence on the size of the roller vortex, as demonstrated in Fig. 3a, and an active vortex almost twice the size (in number of rollers involved NNi) remains approximately as mobile as the smallest one. Indeed, the vortex speed on a flat surface does not show dependence on the size, see Supplementary Figure 2.

Fig. 3
Fig. 3

Potential stiffness manipulation with passive scatterers. a Trapping potential stiffness as a function of number of magnetic rollers comprising the vortex NNi. The number of rollers has a weak influence on the trapping potential. Black solid line is a linear fit. B0 = 5.78 mT, fB = 40 Hz. b Effect of small scatterers (dbead = 150 ± 9 μm) on the vortex trapping potential kp (blue triangles, dashed line marks the value for a pristine vortex). Vortex angular speed ω (inset) decreases with the scatterers number. Black solid lines are linear fits. B0 = 5.78 mT, fB = 40 Hz, NNi = 214. c Pinning of the roller vortex by intermediate bead sizes. Three-fold increase in stiffness of the trapping potential kp for intermediate sized scatterers (red squares, dashed line marks the value for a pristine vortex). The blue region indicates the system response as in b and the green area shows the bead sizes, where even a single scatterer destroys the vortex state. B0 = 7.54 mT, fB = 43 Hz, NNi = 271. Error bars in all panels are the standard deviation of the measurements

To explore the behavior of the magnetic roller vortex in the presence of passive scatterers we added a small amount of nonmagnetic glass spheres (see Methods). Scatterers with diameter less than about 300 μm get dispersed and incorporated inside of the rotating body of the roller vortex. At first glance the roller vortex does not ‘see’ them. Nevertheless, since rollers collide with scatterers they consume part of the vortex energy resulting in overall decrease of the effective vortex angular speed ω with the number of scatterers Nbead, see inset of Fig. 3b. Surprisingly, the stiffness of the confining potential felt by the vortex decreases (Fig. 3b), as if the core of the vortex becomes suddenly more mobile with the number of scatterers. The contradiction is resolved by the fact that the increase in the number of scatterers leads to the instability of the vortex structure so that it adjusts its position by spontaneous re-assembly of itself in a new location rather than by a continuous motion. This leads to apparent sudden vortex core shifts and results in an effective softening of the observed potential stiffness. Continuous increase of the scatterers’ number eventually jeopardizes the existence of the vortex state. Alternatively the vortex state can be destroyed by a single large particle (dbead > 710 μm, green region in Fig. 3c).

Intermediate size beads


interact with the roller vortex in a significantly different manner. These beads get pushed inside of the eye of the vortex core. Once the bead is in the eye, the roller vortex mobility is suppressed and it gets pinned by the inert bead. This is revealed by almost three-fold increase in the effective stiffness of the confining potential kp shown in Fig. 3c for a range of inert particles sizes. The stronger confining potential in the presence of those inert beads in the vortex eye is due to the scattering of the rollers close to the vortex eye on the bead. As all rollers are hydrodynamically coupled and synchronized29 within the vortex, collisions will keep a vortex eye on the bead to minimize obstructed motion of the rollers and sustain the vortex. This observed behavior is similar to pinning of Abrikosov vortices in superconductors by defects and inclusions used to improve transport properties of superconductors36. The property of the magnetic roller vortices to be pinned by certain inert defects provides us with a new tool to manipulate vortex dynamics in active roller suspensions. Conversely, if the inert particle is mobile the vortex can capture it inside of the core and transport with the roller vortex motion. This scenario is demonstrated in Fig. 4b where the inert 500 μm particle has been captured and transported for about 14 by the magnetic roller vortex generated at a flat surface, see Fig. 4 and Supplementary Movie 6. Interestingly, the capture of the passive particle proceeded through intermittent flocking of the vortex’ rollers around it. In a related approach42 the trapping and manipulation of microscopic objects has been realized with the help of an induced hydrodynamic vortex generated by a rotating magnetic micro-wire in a rotational magnetic field. There the trapping force is very local (a few microns) and directional (the trap should approach at a specific angle to the object). In contrast, while each individual roller creates a similar hydrodynamic microvortex and is capable of hydrodynamic trapping of microparticles, the self-assembled roller vortex (collection of cooperating rollers) has a significantly extended trapping range (up to a few millimeters) and allows caging and manipulation of much larger particles compared to a single roller.

Fig. 4
Fig. 4

Inert bead capture and transport by a roller vortex. ac The snapshots demonstrate the capture and transport of a passive glass bead (dbead = 500 ± 25 μm). See also Supplementary Movie 2. The red curve marks the path traveled by the passive bead and the blue circle encloses the particles that are a part of the vortex. The arrow denotes the direction of the vortex rotation. b the scatterer is captured by the vortex through intermittent flocking. The passive particle is transported by the vortex for about 14 s. The vortex spontaneously disintegrates (frame c) to release the bead. The energizing field: B0 = 6.07 mT, fB = 43 Hz. The scale bars are 2 mm

Source link


Please enter your comment!
Please enter your name here