Windmilling clusters of active quadrupoles

Active matter has thrived in recent years, driven both by the insight that it underlies fundamental processes in nature, and by its vast potential for applications. This allows for innovation both inspired by experimental observations, and by construction of novel systems with desired properties. In this paper, we develop a novel system in the search for a new kind of pattern formation: microstructural motifs with orthogonal alignment. Taking a simple active Brownian particle (ABP) model applied to dumbbell-shaped particles, we add a quadrupolar interaction by positioning two antiparallel magnetic dipolar moments on each particle. We find that the phase behavior is determined by the competition between active motion and the orthogonal alignment favored by quadrupolar attraction. By varying these quantities, we are able to tune both the internal structure of the aggregates, and find a surprising stability of triangular aggregates, to the point of clusters of size $N=3$ being strongly overrepresented. Although none of the component particles are chiral, the resulting structures spin in a random, fixed direction due to combination of the polarity of the active motion. This results in an ensemble of windmilling (randomly spinning in a circular motion) aggregates with windmill-like shape (due to the three- or four core component dumbbells). Ultimately, this simple model shows an interesting range of microstructural motifs, with great potential for experimental implementations.

💡 Research Summary

In this work the authors introduce a novel active matter system in which each particle is a dumbbell‑shaped active Brownian particle (ABP) bearing two antiparallel magnetic dipoles, thereby forming a magnetic quadrupole. The particles interact via a short‑range Weeks‑Chandler‑Andersen (WCA) repulsion that enforces a rigid 2:1 aspect ratio, and via a long‑range dipole‑dipole potential that, because of the antiparallel arrangement, yields a quadrupolar interaction. The strength of this interaction is quantified by the dimensionless coupling constant λ = μ₀μ²/(4πk_BTσ³), while activity is characterized by the Péclet number Pe = v₀γ_r, where v₀ is the self‑propulsion speed and γ_r the rotational friction coefficient.

Brownian‑Dynamics simulations were performed in two dimensions with N = 1000 particles at three area fractions (φ = 0.05, 0.15, 0.3). The overdamped Langevin equations were integrated with a time step δt = 0.001, and data were collected after an equilibration period of 5 × 10⁴ δt. The authors systematically varied λ and Pe to map out the collective behavior.

A key finding is that, despite the pairwise quadrupolar potential having its minimum at orthogonal alignment (U_orth ≈ ‑1.56 λ) and being strongly repulsive for parallel alignment (U_par ≈ +2.35 λ), the energetically optimal configuration for three particles is a compact equilateral triangle with each pair at 60°. This yields a per‑pair energy U_tri ≈ ‑1.34 λ and a total cluster energy ≈ ‑4.02 λ, lower than any arrangement that preserves orthogonal contacts for only two of the three pairs. For four particles a square lattice with orthogonal contacts becomes favorable (≈ ‑6.74 λ), but at low N the triangular motif dominates.

When activity is turned on, each particle experiences a constant propulsion force along its long axis. Because the propulsion direction coincides with the dipole orientation, the self‑propulsion introduces a net torque on any cluster that lacks perfect symmetry. Consequently, the over‑represented three‑particle triangles acquire a fixed sense of rotation and behave as “windmilling” clusters: they spin persistently around a common center, reminiscent of a wind‑mill blade. Larger aggregates also rotate, but the triangular clusters are statistically the most abundant.

The authors construct a phase diagram in the (λ, Pe) plane for each φ. Four macroscopic states are identified:

1. Active‑gas (low λ, low Pe): particles remain largely dispersed, with only occasional contacts.
2. Magnetically‑dominated (high λ, low Pe): clusters form with a predominance of orthogonal contacts, similar to passive dipolar systems but with the quadrupolar twist.
3. Triangular‑active (intermediate λ, intermediate Pe): a striking excess of three‑particle triangular clusters that rotate, giving rise to the windmilling phenomenon.
4. Triangular‑magnetic (high λ, intermediate Pe): triangular clusters coexist with larger magnetic aggregates, still showing a peak at N = 3.

Statistical analysis of the cluster‑size distribution p(n) confirms that n = 3 clusters are over‑represented across a wide range of parameters, far beyond what would be expected from simple energetic arguments for larger aggregates. The authors attribute this to a competition between the quadrupolar energy landscape (which favors orthogonal contacts) and the geometric constraint that prevents three particles from simultaneously satisfying orthogonal pairwise minima. The activity then selects the compact triangular arrangement because it allows a coherent torque to develop, stabilizing the rotating state.

Methodologically, the paper provides a clear description of the simulation protocol, including the choice of reduced units (σ = 1, k_BT = 1, D_rot = 0.1), the mapping of λ to physical magnetic moments, and the implementation of the dipolar interaction within the ESPResSo package. The authors also discuss the importance of a sufficiently long equilibration for high λ, low Pe runs, where magnetic relaxation is slow.

The discussion highlights the experimental relevance of the model. Fabricating micron‑scale dumbbells with two oppositely magnetized ends (e.g., by coating each lobe with a different ferromagnetic material) would realize the quadrupolar interaction. Activity could be introduced via catalytic self‑propulsion, light‑driven phoretic motion, or magnetic field‑induced rotation. The ability to tune λ (by varying magnetic moment or particle size) and Pe (by adjusting fuel concentration or field strength) offers a practical route to explore the predicted windmilling clusters in the laboratory.

In conclusion, the study demonstrates that combining a quadrupolar magnetic interaction with active propulsion yields a previously unobserved microstructural motif: rotating three‑particle “windmills”. This challenges the conventional view that non‑chiral active particles only form static aggregates or flocking states, and opens new avenues for designing self‑assembled micromachines, magnetic fluids with controllable rheology, and testbeds for non‑equilibrium statistical physics.