Scale-free chaos in the confined Vicsek flocking model

Inertial spin model of flocking with position-dependent forces

We propose an extension to the inertial spin model (ISM) of flocking and swarming. The model has been introduced to explain certain dynamic features of swarming (second sound, a lower than expected dynamic critical exponent) while preserving the mechanism for onset of order provided by the Vicsek model. The inertial spin model (ISM) has only been formulated with an imitation ("ferromagnetic") interaction between velocities. Here we show how to add position-dependent forces in the model, which allows to consider effects such as cohesion, excluded volume, confinement, and perturbation with external position-dependent field, and thus study this model without periodic boundary conditions. We study numerically a single particle with an harmonic confining field and compare it to a Brownian harmonic oscillator and to a harmonically confined active Browinian particle, finding qualitatively different behavior in the three cases.

Collective motion of chiral particles in complex noise environments

Collective motion of chiral particles in complex noise environments is investigated based on the Vicsek model. In the model, we added chirality, along with complex noise, affecting particles clustering motion. Particles can only avoid noise interference in a specific channel, and this consideration is more realistic due to the complexity of the environment. Via simulations, we find that the channel proportion, p, critically influences chiral particle synchronization. Specifically, we observe a disorder-order transition at critical [Formula: see text], only when [Formula: see text], the system can achieve global synchronization. Combined with our definition of spatial distribution parameter and observation of the model, the reason is that particles begin to escape from the noise region under the influence of complex noise. In addition, the value of [Formula: see text] increases linearly with velocity, while it decreases monotonically with the increase in chirality and interaction radius. Interestingly, an appropriate noise amplitude minimizes [Formula: see text]. Our findings may inspire novel strategies to manipulate self-propelled particles of distinct chirality to achieve desired spatial migration and global synchronization.

Power laws of natural swarms as fingerprints of an extended critical region

Collective biological systems display power laws for macroscopic quantities and are fertile probing grounds for statistical physics. Besides power laws, natural insect swarms present strong scale-free correlations, suggesting closeness to phase transitions. Swarms exhibit imperfect dynamic scaling: their dynamical correlation functions collapse into single curves when written as functions of the scaled time tξ^{-z} (ξ: correlation length, z: dynamic exponent), but only for short times. Triggered by markers, natural swarms are not invariant under space translations. Measured static and dynamic critical exponents differ from those of equilibrium and many nonequilibrium phase transitions. Here we show the following: (i) The recently discovered scale-free-chaos phase transition of the harmonically confined Vicsek model has a novel extended critical region for N (finite) insects that contains several critical lines. (ii) As alignment noise vanishes, there are power laws connecting critical confinement and noise that allow calculating static critical exponents for fixed N. These power laws imply that the unmeasurable confinement strength is proportional to the perception range measured in natural swarms. (iii) Observations of natural swarms occur at different times and under different atmospheric conditions, which we mimic by considering mixtures of data on different critical lines and N. Unlike results of other theoretical approaches, our numerical simulations reproduce the previously described features of natural swarms and yield static and dynamic critical exponents that agree with observations.

Scale-Free Chaos in the 2D Harmonically Confined Vicsek Model

Animal motion and flocking are ubiquitous nonequilibrium phenomena that are often studied within active matter. In examples such as insect swarms, macroscopic quantities exhibit power laws with measurable critical exponents and ideas from phase transitions and statistical mechanics have been explored to explain them. The widely used Vicsek model with periodic boundary conditions has an ordering phase transition but the corresponding homogeneous ordered or disordered phases are different from observations of natural swarms. If a harmonic potential (instead of a periodic box) is used to confine particles, then the numerical simulations of the Vicsek model display periodic, quasiperiodic, and chaotic attractors. The latter are scale-free on critical curves that produce power laws and critical exponents. Here, we investigate the scale-free chaos phase transition in two space dimensions. We show that the shape of the chaotic swarm on the critical curve reflects the split between the core and the vapor of insects observed in midge swarms and that the dynamic correlation function collapses only for a finite interval of small scaled times. We explain the algorithms used to calculate the largest Lyapunov exponents, the static and dynamic critical exponents, and compare them to those of the three-dimensional model.

Soliton approximation in continuum models of leader-follower behavior

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.

Mean-field theory of chaotic insect swarms

The harmonically confined Vicsek model displays qualitative and quantitative features observed in natural insect swarms. It exhibits a scale-free transition between single and multicluster chaotic phases. Finite-size scaling indicates that this unusual phase transition occurs at zero confinement [Phys. Rev. E 107, 014209 (2023)2470-004510.1103/PhysRevE.107.014209]. While the evidence of the scale-free-chaos phase transition comes from numerical simulations, here we present its mean-field theory. Analytically determined critical exponents are those of the Landau theory of equilibrium phase transitions plus dynamical critical exponent z=1 and a new critical exponent φ=0.5 for the largest Lyapunov exponent. The phase transition occurs at zero confinement and noise in the mean-field theory. The noise line of zero largest Lyapunov exponents informs observed behavior: (i) the qualitative shape of the swarm (on average, the center of mass rotates slowly at the rate marked by the winding number and its trajectory fills compactly the space, similarly to the observed condensed nucleus surrounded by vapor) and (ii) the critical exponents resemble those observed in natural swarms. Our predictions include power laws for the frequency of the maximal spectral amplitude and the winding number.