Self-reverting vortices in chiral active matter

There is currently a strong interest in the collective behavior of chiral active particles that can propel and rotate themselves. In the presence of alignment interactions for many chiral particles, chiral self-propulsion can induce vortex patterns in the velocity fields. However, these emerging patterns are non-permanent, and do not induce global vorticity. Here we combine theoretical arguments and computer simulations to predict a so-far unknown class of collective behavior. We show that, for chiral active particles, vortices with significant dynamical coherence emerge spontaneously. They originate from the interplay between attraction interactions and chirality in the absence of alignment interactions. Depending on parameters, the vortices can either feature a constant vorticity or a vorticity that oscillates periodically in time, resulting in self-reverting vortices. Our results may guide future experiments to realize customized collective phenomena such as spontaneously rotating gears and patterns with a self-reverting order.


INTRODUCTION
Chirality refers to the property of objects to be nonsuperimposable on their mirror images.The concept originated in the mid-19th century and it is attributed to the chemist Louis Pasteur, who observed that crystals of tartaric acid exist in two distinct, non-superimposable forms, which he referred to as "right-handed" and "lefthanded."Also more than a century ago, Bronn, Jennings, and others realized that shape-asymmetric motile microorganisms generically follow chiral trajectories [1,2], i.e. they do not only self-propel but they also self-rotate, showing circular trajectories.Recently, the discovery of synthetic colloidal microswimmers in the 21st century [3][4][5] has stimulated a significant interest in chiral self-propelled particles [6,7]: Like their biological counterparts, these particles also generically follow circular trajectories if they feature a shape-anisotropy [8][9][10][11][12][13][14] or are torqued by an external field [15].In addition, it is now known that chirality can emerge due to hydrodynamic interactions with walls or interfaces, as in bacteria [16][17][18], or due to memory effects in viscoelastic environments and droplet swimmers [19][20][21].
In the present work, we combine theoretical arguments and particle-based simulations to predict the existence of a so-far unknown class of structures in chiral active matter.First, and perhaps least surprisingly, for low chirality (low self-rotation frequency) we find that attractive chiral active particles (Fig. 1 a) without alignment interactions form moving rigid clusters that feature full velocity-alignment of the contained particles [52] and spatial velocity correlations [53,54] but vanishing vorticity.However, for high chirality, we observe a transition to a rotation pattern that is characterized by a persistent and time-independent vorticity (Fig. 1 b) and is termed permanent vortex state.This state can be viewed as the superposition of the translational motion characterizing the previous state and an additional collective rotation due to chirality, which transfers from the single particle to the collective level.For even higher chirality the rotation pattern again changes and the vortex starts to dynamically revert itself, exhibiting periodic transitions between vortex and antivortex configurations (Fig. 1 c).We refer to these structures as self-reverting vortices.The occurrence of this state is a consequence of the competition between chirality and isotropic interactions, which suppresses the tendency of a cluster to collectively rotate.

Model for chiral particles
To concretely investigate these states, we consider a system of N interacting active chiral particles with mass m, where each particle is governed by underdamped equations of motion for their positions, x i , and velocities, v i = ẋi .Every particle is in contact with a thermal bath at temperature T and experiences a frictional force, γv i , with friction coefficient γ.Activity is incorporated in the dynamics as a stochastic force, which imparts to each particle a constant swim velocity, v 0 , along with an orientation vector, n i = (cos θ i , sin θ i ).Here, θ i are the orientational angles and, in accordance with the active Brownian particle (ABP) model [55] describing circular swimmers [56][57][58][59][60][61][62], evolves as Brownian noise with a constant drift angular velocity, ω.The latter is also known as particle chirality and is responsible for circular trajectories [56].Thus, the system's dynamics can be expressed as (1)b Here, D r represents the rotational diffusion coefficient, and ξ i and η i denote white noises with zero average and unit variance.The particle chirality ω determines the characteristic radius of the circular trajectory displayed by a single active chiral particle, specifically v 0 /ω.In this system, the absence of torques between particles results in their sole interaction through the force The shape of the interaction potential U (r) is obtained by truncating and shifting an attractive Lennard-Jones potential, U LJ (r) = 4ϵ σ r 12 − σ r 6 .The potential U (r) is therefore defined as U (r) = U LJ (r)−U LJ (3σ) for r ≤ 3σ, and zero otherwise.Here, σ signifies the nominal particle diameter, while ϵ stands for the energy scale of the interactions.The interparticle attraction is sufficiently large to guarantee that the cluster structure remains stable.The system is characterized by three primary time scales: the inertial time, τ I = m/γ, determining the velocity relaxation; the persistence time, τ = 1/D r , which dictates the duration required for active particles to randomize their orientations; and the time 1/ω necessary for the orientation to complete a full rotation due to chirality.We remark that our model considers self-propulsion and chirality, i.e. self-rotations, as two independent mechanisms.Indeed, even if these two propulsions are often related in chiral active colloids, this is not the case in other physical systems, for instance active granular particles with an intrinsic chirality and spinners where the self-propulsion is even absent.

Theoretical prediction for vortex states
To understand the collective behavior of attractive chiral active particles, we first develop a mapping showing that the overall dynamics are governed by the competi- (d)-(e): Illustrations of a chiral active dumbbell, anchored to one of the particles.The permanent vorticity state occurs when the chiral radius v0/ω (grey line centered on a diamond) is larger than the dumbbell size (black line centered on a circle), allowing the mobile particle to complete a full rotation around the other (orange star).The self-reverting vorticity state occurs in the opposite regime, such that the mobile particle cannot complete a full rotation around the other.Indeed, its self-propulsion is reversed before completion resulting in partial clockwise and counterclockwise rotations (purple star).
tion of velocity alignment and an effective Lorentz force.Specifically, Eqs.(1) can be mapped to alternative dynamics using an exact change of variables and a lattice approximation, applicable to strongly attractive active particles in the large persistence time regime (see Methods).The evolution of the particle velocity v i is effectively described by: vi Here, j encompasses neighboring particles to the i-th one, z represents the normal vector to the plane of motion, and the elements of the matrix J j are reported in the Methods.Equation (2) holds in the large persistence regime and reveals that the system's behavior is primarily governed by two distinct forces.The first one, independent of chirality, manifests as an effective alignment force emerging from the interplay between interactions and activity.It accounts for the observed velocity alignment and, indeed, is minimized in three particle configurations: i) full alignment; ii) vortex; iii) antivortex (see Methods).In the absence of chirality, there is no preference among these configurations.The second term in Eq. ( 2) is solely induced by chirality, being ∝ ω, and op-erates as an effective Lorentz force.This indicates that chirality influences the dynamics of an active particle, akin to an effective magnetic field, responsible for particle rotations.When ω is comparable to J j , the effective magnetic field selectively promotes vortex or antivortex states for negative and positive ω, respectively.Consequently, this analytical argument predicts the spontaneous emergence of substantial vorticity in the system.

Simulations unveil self-reverting vortices
To observe the predicted vortex-states, we perform simulations within a box of size L under periodic boundary conditions, ensuring that the particle packing fraction ϕ = (N/L 2 )σ 2 π/4 = 0.3 remains constant.It is crucial to emphasize that our findings pertain to the large persistence regime (τ /τ I ≫ 1), resulting in effectively overdamped dynamics (see Methods).Given this choice, the same results can be obtained by considering an overdamped dynamics if the thermal noise is sufficiently small.The condition τ /τ I ≫ 1 signifies that the persistence length v 0 τ is the dominant length scale in the system, notably larger than the cluster size , where N c is the number of particles in the cluster.We remark that in the opposite small persistence time regime τ /τ I ≪ 1, the system is close to equilibrium and the active force behaves as thermal noise (see Methods).Thus, no collective motion can be observed in this regime.Therefore, we conduct a numerical study by keeping fixed τ /τ I ≫ 1.In addition, the dynamical states shown here are obtained only in the regime of large attractions compared to thermal noise strength and activity, i.e. when the typical potential energy due to the interparticle interactions is large compared to the thermal energy and the kinetic energy associated with selfpropulsion ≈ mv 2 0 /2.Indeed, without this condition, the cluster is not stable because particles are able to leave it and therefore it is not possible to observe collective motion.Here, to investigate the influence of chirality, we vary the associated dimensionless parameter, the reduced chirality ωτ , and examine different cluster sizes L c .
We discover phenomena in active systems uniquely induced by circular motion and attractive forces.Reduced chirality ωτ fosters collective rotational motion, with the entire cluster tracing persistent circular trajectories (see Supplementary Movie 1).For further increasing values of ωτ , the cluster displays spinning dynamics rather than a circular trajectory.Indeed, the center of mass of the cluster undergoes rotations with a characteristic radius smaller than the cluster size L c (see Supplementary Movie 2).To characterize rotational motion, we monitor the evolution of the spatial average of the vorticity field ⟨Ω⟩, defined as which reads zero for particle velocities aligned in the same direction but assumes positive and negative values for antivortex and vortex configurations, respectively.In the limit of vanishing chirality ωτ (black curve in Fig. 2 b), the time-trajectory of ⟨Ω(t)⟩ fluctuates around zero, signifying the absence of a preferred vorticity.By contrast, as ωτ increases, ⟨Ω(t)⟩ exhibits minor time-fluctuations around a value greater than zero, indicating a positive spatial average vorticity aligned with the single-particle chirality (Fig. 2 b).The breaking of rotational symmetry of a single chiral particle propagates to the collective level, resulting in a non-zero global vorticity: These configurations are identified as permanent vorticity states.Conversely, when the cluster exhibits spinning dynamics, ⟨Ω(t)⟩ displays periodic-time oscillations (Fig. 2 c).This implies that particle velocities periodically switch between vortex and antivortex configurations, i.e. the cluster exhibits a self-reverting vorticity.This phenomenon is a consequence of the additional time scale introduced by chirality, as confirmed by the oscillation period which scales as ∼ 1/ω.

Mechanism behind self-reverting vorticity
Permanent vorticity and self-reverting vorticity states can be intuitively explained by considering a chiral particle anchored to a fixed point with size determined by the strong attraction, σ.Chirality enables the selfpropulsion force to persistently rotate at a frequency of 1/ω and, consequently, induces circular motion in the particle around the fixed point.If the radius of the circular trajectory is larger than the distance with the fixed point, v 0 /ω > σ (low chirality), the particle performs complete, persistent rotations, in the direction promoted by chirality (Fig. 2 d).This simple mechanism generates the permanent vorticity state at the collective level.By contrast, in the opposite regime of large chirality (v 0 /ω < σ), chirality completely reverses the direction of the active force before the particle completes a rotation of π radians around the immobile particle.Consequently, the particle moves backward and forward, effectively alternating between clockwise and counterclockwise rotations (Fig. 2 e).This explains the observed self-reverting vorticity state at the collective level.This idea can be further supported by calculating the total torque M acting on the cluster, which is dominated by the outer particle layer at distance R from the middle of the cluster (see Methods) Here r i is a vector pointing from the center of the cluster to the i-particle position.After time π/ω, each n i rotates by π and thus also its spatial average.However, if n i rotates with a period π/ω smaller than the period of r, M is continuously subject to sign changes before a full cluster rotation: the cluster displays self-reverting vorticity.By contrast, in the opposite regime, M never changes sign and the system displays a permanent vorticity state.

State diagram
Our findings are systematically explored by varying the cluster size L c and the reduced chirality ωτ on a state diagram (Fig. 2 a).We identify different states with different colors: non-permanent vorticity states (black dots) when the average vorticity is smaller than its time fluctuations; permanent vorticity states (orange dots), when the previous condition is fulfilled; self-reverting vorticity state (purple dots) when vorticity displays periodic oscillations.Consistent with our intuitive explanation, the transition line between permanent and self-reverting vorticity state occurs when the cluster size L c = σ √ N c approaches the typical radius of the chiral trajectory ∼ v 0 /ω.This argument suggests the following scaling law which fairly reproduces our numerical results (Fig. 2 a).This scaling law suggests that a larger cluster size favors the self-reverting vorticity state over the permanent vorticity state.Additionally, it is worth noting that an increase in cluster size promotes permanent vorticity states over states without vorticity.This is because the time fluctuations of the vorticity field decrease with increasing N c .We remark that the crossover between different states is not a sharp transition but occurs smoothly.Indeed, the regions in Fig. 2 a are obtained by following the threshold criterion defined in the methods section.
To quantitatively characterize the different states, we consider the time average of the spatial average vorticity as a order parameter This observable shows a non-monotonic behavior with the reduced chirality ωτ for different values of the cluster size L c (Fig. 3 a).For vanishing ωτ , the absence of permanent vortices (black state in Fig. 2 a) induces rather small values of ⟨Ω⟩.The increase of ωτ enhances the value of ⟨Ω⟩ until it becomes larger than its typical time fluctuations and the system approaches the permanent vorticity state.In this regime, ⟨Ω⟩ monotonically increases until a maximum is achieved.This maximum occurs before the system approaches the self-reverting vorticity state, for which the periodic oscillations sharply lead to vanishing values of ⟨Ω⟩.The amplitude of these oscillations is investigated by evaluating the time average of the modulus of the spatial average vorticity, ⟨|Ω|⟩ (Fig. 3 b).This observable monotonically increases with ωτ until the self-reverting vorticity state is approached when ⟨|Ω|⟩ saturates to a constant value.This implies that the amplitude of the vorticity oscillations remains constant with ωτ and does not significantly change with the cluster size.Finally, the oscillation period (Fig. 3 c) decreases with the reduced chirality as 1/(ωτ ).This scaling confirms our intuitive explanation of this phenomenon: Before completing a full rotation, the orientation of chiral active particles is reversed after a time period ∼ 2π/ω.This implies that these periodic oscillations are uniquely induced by chirality.

CONCLUSIONS
The central insight of this work is that the presence of attractions in chiral active matter without alignment interactions induces self-organized vortices involving coherent dynamics of adjacent particles.These vortices can either be persistent or show a periodically oscillating vorticity, leading to patterns that self-revert their order.
The theoretical arguments developed here (for instance Eq. ( 4)) could shed light on the link between chiral active systems and materials with odd properties [63][64][65][66], such as crystals characterized by odd elasticity [67] and liquid governed by odd viscosity [68][69][70].Indeed, living chiral crystals exhibit self-sustained chiral oscillations as well as various unconventional deformation response behaviors recently predicted for odd elastic materials [63].Our argument rationalizes these findings, suggesting that self-propulsion plays the role of the transverse neighbor forces typical of odd materials.
Even if here collective phenomena spontaneously emerge without alignment interactions, it could be interesting to evaluate the effects of explicit alignment mechanisms on chiral active particles at high density, in cluster configurations.This is a rather common scenario in selfpropelled colloids that can behave as chiral microswimmers by simply introducing a rotational asymmetry in their body [8].
This finding opens the door to the observation of customizable collective phenomena.They have the potential to inform the design and optimization of particlebased micromotors.Instead of creating asymmetric gears powered by active particles [71][72][73], spontaneous gear rotation can be achieved by harnessing chirality in active matter [74].Our study could inspire experiments across a wide range of chiral active matter experiments, such as high-density chiral active colloids [8] attracting by means of Van-der-Waals interactions, or chiral active granular particles [75][76][77] which can be connected by springs to create crystal-like configurations [78].

Derivation of the theoretical prediction, Eq. (2)
To derive Eq. ( 2), in the following we employ a similar idea as has been used in Ref. [79] for straight active particles.As we will see, accounting for chirality, leads to an additional term in the resulting equation that competes with the effective alignment that has been found in Ref. [52].This competition is at the heart of the phenomenology which we predict and observe in the present article, as discussed in the main text.

Mapping on the dynamics to an effective description
Before proceeding to the exact mapping, it is convenient to express the dynamics of the activity in Cartesian coordinates.By applying Ito calculus rules, Eq. (1)b can be expressed in Ito's convention as Here, the vector ξ i = (0, 0, ξ i ) consists only of the third component orthogonal to the plane where the particle motion takes place, namely the xy plane.In this way, the noise vector can be expressed in a compact form as ξ i = zξ i , where z is the unit vector normal to the xy plane.
Even if the dynamics (1) is underdamped, the extremely small value of the reduced inertia (i.e. of the inertial time compared to the persistence time) allows us to take the overdamped regime, vi ≈ 0, so that the equation of motion for chiral active particles is effectively given by In addition, the small value of the reduced temperature (i.e. the small value of T compared to the self-propulsion velocity square), allows us to drop the passive Brownian motion term.By applying the time-derivative to Eq. ( 8)a with T = 0, and by defining the velocity variable, v i = ẋi , we obtain where we have assumed Einstein's convention on repeated indices.By replacing ṅi by Eq. ( 8)b immediately we have (10) Now, we proceed by replacing n by the Eq. ( 8)a (again with T = 0), obtaining Dynamics ( 11) is mathematically equivalent to Eq. ( 1).
In order to proceed analytically, we consider further approximations described in the next subsections.

Lattice approximation for solid-like configurations
The strong interparticle attractive interactions induce almost-perfect solid-like configurations with an almostperfect hexagonal order.This allows us to consider the lattice approximation by fixing the particle positions on the vertices of a triangular lattice.In this way, every particle is characterized by six neighbors.This implies that we consider systems sufficiently large to neglect the contribution of the outer layer of particles, whose number scales as ≈ √ N c .In this approximation, interparticle forces are perfectly balanced because of the lattice translational invariance.As a consequence, we need only to evaluate the second derivative of the total potential in Eq. ( 11), Here, r ij is the distance between particle i-th and particle j-th, and the sum, * is restricted over the six neighbors of the target particle i.The truncation at first neighbors works if the potential is short-range, as in the Lennard-Jones potential numerically considered in the numerical simulations.To proceed further, we can calculate the spatial components of the Hessian matrix, which is a 2×2 matrix, in two dimensions.In particular, we have where we have denoted the spatial components by Greek upper indices, and r α ij = r α i − r α j , with α = x, y.Here, each prime on the potential U means a spatial derivative.
We remark that the potential depends only on the interparticle distance and, thus the following property holds: To switch to a more suitable description accounting for the lattice symmetry, it is convenient to express the Cartesian components in polar coordinates, such that r x ij /|r ij | = cos (δ j ) and α = x, y of r y ij /|r ij | = sin (δ j ).Here, δ j the angle between the r ij vector and the x-axis.
The triangular lattice structure implies the target particle i has 6 first neighbors, uniquely identified by δ j = δ 0 +jπ/3 with j = 0, 1, ..., 5.The phase δ 0 represents the orientation of the hexagon with respect to the reference frame that can be set to zero without loss of generality.In this way, by denoting |r ij | = r (the lattice constant), we can rewrite the components of the Hessian matrix as By summarizing, the left-hand-side of Eq. ( 12) can be expressed as where the matrix J j has elements Because of the following properties: we can conclude that the force (18) has the shape of an effective alignment interaction between the particle i and its 6 first neighbors.We remark that to apply our theory the potential has to be differentiable: in particular first and second derivatives of the potential should be defined.Our choice of Lennard Jones potential, truncated at 3σ as usual in numerical studies, does not represent a problem for the applicability of the theory.Indeed, we resorted to the first-neighbors approximation, which allows us simply to select the interactions between the six-neighboring particles, which are at distance ≈ σ < 3σ where the first two derivatives of the potential are well-defined.

Effect of chirality
By summarizing the results of previous sections, the behavior of chiral active particles is well-described by the following dynamics obtained after performing the lattice approximation Equation ( 23) corresponds to the dynamics (2).The first term is an effective friction force, whose friction coefficient is determined by the inverse of the persistence time τ , this term here dissipates the energy injected by the thermal bath, whose amplitude is determined by v 0 / √ τ .The unconventional shape of this noise term is due to the choice of active Brownian particle dynamics which conserves the modulus of the active force and here involves the cross product with n i .It is worth noting that by excluding chirality and interactions, the velocity scale is purely determined by v 0 while τ plays a negligible role, as expected.Both terms are vanishing in the large persistence limit τ γ → ∞, becoming subleading in the dynamics.The third term in the dynamics − 1 γ * j J j (v i − v j ) accounts for particle interactions and has the shape of an effective alignment interaction term spontaneously emerge from this analytical calculation.Indeed, the particle i feels a force proportional to the difference between the velocities of neighboring particles, −(v i − v j ), i.e. particle i tends to align its velocity to those of neighboring particles.As stated in the results, this effective alignment force is minimized in three different configurations: i) aligned velocities; ii) vortexdistributed velocities; and iii) antivortex-distributed velocities.The three configurations are illustrated in Fig. 4. In i), the tagged particle velocity v i = v is equal to any neighboring particle velocity (Fig. 4 a).Consequently, each term of the alignment force * j J j • (v j − v i ) independently vanishes.In ii) and iii), the tagged particle velocity is zero while the six neighboring particle velocities are distributed on a vortex (Fig. 4 b) and an antivortex configuration (Fig. 4 c), respectively.Thus, particles on opposite vertices of the hexagon have equal velocities with opposite directions which perfectly balance.This implies that * j or, in other words, the effective alignment interaction is not only minimized by aligned velocities but also by vortex and antivortex configurations.Finally, the last term ωv i ×z accounts for the role of chirality.Such a force term has the shape of an effective magnetic field with amplitude ω and it is responsible at the single particle level for particle rotations.Intuitively, this term selects vortex or antivortex configurations depending on the sign of the chirality ω.
We also remark that in our theory we resort to a linearization of the force between different particles.This is possible because of the solid structure.In principle, perturbation theory can be applied to mass defects [80] or non-linear potentials with a weak non-linearity, such that the force F ≈ −k 0 x − k 1 |x| 2 x, with k 1 /k 2 ≪ 1.In this case, we expect that the theory could quantitatively provide a correction to our results without changing the observation of the three dynamical states.

Small persistence time regime
In the small persistence time regime τ ≪ τ I , the cluster does not show any coherent motion and simply diffuses.Indeed, in this case, the active force γv 0 n changes fast in its direction, and it can be approximated by an effective Brownian motion.Therefore, in this regime, chirality plays a negligible role and consequently, the three states we have identified, i.e. non-permanent vorticity state, permanent vorticity state, and self-reverting vorticity state, cannot be observed.
This conclusion can be analytically derived by considering the dynamics (1)a with active force evolving in Cartesian coordinates (7).In particular, it is convenient to express the activity dynamics by resorting to a matrix formalism where B is a matrix with components In the small persistence time regime, τ ≪ τ I , τ is the faster time scale and we can take the overdamped limit in the equation for n, by setting ṅ = 0: where B −1 is the inverse of B with components By substituting Eq.( 27) in the dynamics (1)a, we obtain As a consequence, the active force simply behaves as a white noise which cannot induce the non-equilibrium collective motion observed in the regime of large persistence time.

Dimensionless dynamics and dimensional parameters
Simulations are performed by considering Eqs.(1) with rescaled variables.Particle positions are rescaled with the particle diameter σ, so that x ′ = x i /σ, while time is rescaled with the persistence time τ = 1/D r , such that t ′ = t/τ .With this choice, Eqs. ( 1) can be integrated using the Euler method with time-step dt ′ = dt/τ = 10 −6 /τ and reduces to  where i (t) are the increment of particle position and velocity after a time-step dt ′ , while δθ ′ i (t) = θ ′ i (t+dt)−θ ′ i (t) represents the time integral of the orientational angle of the particle i.In addition, dη ′ i (t ′ ) and dξ i (t ′ ) are two dimensionless Wiener processes with zero average that can be numerically generated by Gaussian numbers with unit variance and we have used the definition of the inertial time τ I = m/γ.The dynamics ( 30) is governed by five dimensionless parameters that are listed and commented on below: (i) Reduced inertial time τ I /τ = m/(γτ ) = 10 −6 which determines the velocity relaxation in units of persistence time.
(v) Reduced chirality ωτ , which governs the time-scale associated with chirality and it is varied in the simulations to address its effect.
With this choice of parameters in particular τ I /τ = 10 −6 , the dynamics (1) (or the dimensionless Eqs. ( 30)) are effectively in the overdamped regime.However, since Eqs.( 30) is an underdamped equation of motion, velocities v i are well-defined.The underdamped choice is particularly convenient to calculate velocity and vorticity fields because they remain well-defined even in the presence of thermal noise.However, the numerical results reported in this paper can be also observed with an overdamped active model if the thermal noise is sufficiently small.In addition, in numerical simulations, we explore different cluster size square σ 2 N = 113, 226, 452, 904, 1809 and packing fraction ϕ = N/L 2 πσ 2 = 0.35.The size of the box L is chosen accordingly.The system spontaneously evolves to a state characterized by a unique cluster because of attractive interactions.However, depending on the total number of particles in simulations, the system could take a long transient time to reach the steady state.Thus, when needed, simulations were directly initialized in the cluster configuration.

Details on the distinction beween the different states
In Fig. 2 a, we have distinguished between three states: i) Non-permanent states (black dots in Fig. 2 a).This state is characterized by fluctuating vorticity, from negative to positive values, which is compatible with configurations with negligible chirality.
ii) Permanent vorticity states (orange dots in Fig. 2 a).In this state, the cluster is characterized by a permanent vorticity and displays a permanent rotating trajectory aligned to the particle chirality.
iii) Vortex-antivortex state (violet dots in Fig. 2 a).This state shows the self-reverting vorticity observed and the cluster spinning dynamics.
States i), ii), and iii) are characterized by a continuous crossover rather than a sharp phase transition.This feature is already evident by evaluating the time-averaged value of the spatial average vorticity ⟨Ω⟩τ (Fig. 3 a) and the time-average of the modulus of the rescaled spatial average vorticity ⟨|Ω|⟩τ (Fig. 3 for t → ∞.Finally, configurations belonging to state iii) again satisfy the condition Eq. (31).However, at variance with state 1, ⟨Ω(t)⟩ switches from negative to positive values periodically in time.

Derivation of the theoretical argument (4)
To calculate the total torque on the cluster, let us consider the total force exerted by each microscopic active particle, F i , given by the sum of attractive interactions and active forces where the sum * j runs over the six neighbors of the ith particle.To calculate the torque due to the particle i-th, we have to apply the vector product of the relative particle position calculated from the center of the cluster r i By summing over i, the contribution of the internal force vanishes by symmetry and the total torque reads By assuming that clusters have spherical shapes, the torque can be decomposed as where r ′ is the radial coordinate with respect to the center of the cluster and R is the cluster radius.As a consequence, m(r ′ ) is the torque due to the particles at distance r ′ from the cluster center and can be expressed as (37) This expression corresponds to Eq. ( 4) after recognizing that, as a first approximation, M ≈ m(R), since the particles in the outer layer provide the larger contribution to the torque being those at the larger distance from the center.

Description of the Supplementary Movies
Supplementary Movie 1 and Supplementary Movie 2 report the time evolution of the system from two simulations in the steady-state, for reduced chirality ωτ = 10 and 2 × 10 2 , respectively.The first shows a typical configuration in the permanent vorticity state, while the second displays a typical one in the self-reverting vorticity state.In both cases, the cluster size square is given by L 2 c = 904, corresponding to a number of particles in the cluster, N c = 904.Each movie consists of two adjacent videos lasting for a total time interval 10 2 τ .In the left video, particles are plotted as circular points and colored according to their orientational angle θ i .In the right video, we plot the coarse-grained fields: black arrows are used to denote the velocity field v(r), while the color gradient shows the vorticity field Ω(r).In both cases, orientations θ i are random since there are no torques between different particles, and in general, the velocity field v(r) shows global alignment characterized by vortex-like configurations.The main differences between the two movies appear in the global cluster motion and in the time evolution of the vorticity field Ω(r).In Supplementary Movie 1, the cluster shows a global counterclockwise rotation and Ω(r) is locally and globally always larger than zero.In Supplementary Movie 2, the cluster is characterized by spinning dynamics while Ω(r) displays local oscillations between negative and positive values.

DATA AVAILABILITY
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon request, while Supplementary Movie 1 and Supplementary Movie 2 are uploaded as Supplemental Material.

Figure 1 .
Figure 1.Chirality-induced collective motion.(a): Illustrations of chiral active particles.The particle orientation is indicated by a dark blue cap with L-shape responsible for the particle chirality.The black and blue arrows serve as a schematic representation of particle velocity and chirality.We sketch the typical rotating trajectory of a chiral active particle and illustrate a chiral active cluster maintained by attractive interactions.(b): Time series of snapshot configurations showing the permanent vorticity state, as revealed by the vorticity field (colors).The cluster exhibits collective rotations, following circular trajectories and displaying counterclockwise vorticity.(c): Time series of snapshots for the cluster exhibiting spinning dynamics and revealing the self-reverting vorticity state, i.e. a periodic alternation between vortex and antivortex configurations.Arrows indicate the velocity field while colors denote the vorticity value.(b)-(c) are obtained with cluster size square L 2 c = 904 and reduced chirality ωτ = 10, 10 2 , respectively.The remaining parameters of the simulations are: τI /τ = 10 −6 , Pe = τ v0/σ = 50, τ 2 ϵ/(σ 2 m) = 5 × 10 3 , τ 2 T /(mσ 2 ) = 10 −5 .

Figure 2 .
Figure 2. Vorticity states.(a): State diagram for chiral active particles in the plane of reduced chirality ωτ and cluster size square L 2 c .Black, orange, and purple colors denote negligible vorticity, permanent vorticity, and self-reverting vorticity states, respectively.(b)-(c): Spatial average chirality, ⟨Ω⟩τ , as a function of the rescaled time t/τ .The three time-trajectories in (b) and (c) correspond to the stars in (a): Specifically, the black, orange, and purple curves are obtained for ωτ = 5 × 10 −1 , 10, 10 2 with L 2 c = 452.The remaining parameters of the simulations are: τI /τ = 10 −6 , Pe = τ v0/σ = 50, τ 2 ϵ/(σ 2 m) = 5 × 10 3 , and τ 2 T /(mσ 2 ) = 10 −5 .(d)-(e):Illustrations of a chiral active dumbbell, anchored to one of the particles.The permanent vorticity state occurs when the chiral radius v0/ω (grey line centered on a diamond) is larger than the dumbbell size (black line centered on a circle), allowing the mobile particle to complete a full rotation around the other (orange star).The self-reverting vorticity state occurs in the opposite regime, such that the mobile particle cannot complete a full rotation around the other.Indeed, its self-propulsion is reversed before completion resulting in partial clockwise and counterclockwise rotations (purple star).

Figure 4 .
Figure 4. Illustrations of configurations that minimize the effective alignment interactions.The tagged particle is in the middle while neighboring particles are placed on the vertices of a hexagon.Velocities are represented by black arrows.(a): Aligned velocities; (b) Vortex-distributed velocities; (c) Antivortex-distributed velocities.
b.In particular, the first two states are distinguished by comparing the time fluctuations and time average of the total vorticity field ⟨Ω(t)⟩.In particular, configurations that belong to state i) are characterized by a time-standard deviation of ⟨Ω(t)⟩ larger than its average, i.e. by the following