Search

Flocking systems are known to be strongly out of equilibrium. Energy input occurs at the individual level to ensure self-propulsion, and the individual motility in turn contributes to ordering, enhancing information propagation and strengthening collective motion. However, even beyond ordering, a crucial feature of natural aggregations is response. How, then, do off-equilibrium features affect the response of the system? In this work, we consider a minimal model of flocking and investigate response behavior under directional perturbations. We show that equilibrium dynamical fluctuation-dissipation relations between response and correlations are violated, both at the local and at the global level. The amount of violation peaks at the ordering transition, exactly as for the entropy production rate. Entropy is always produced locally and connected to the local fluctuation-dissipation violation via Harada-Sasa relationships. However, cooperative mechanisms close to the transition spread off-equilibrium effects to the whole system, producing an out of equilibrium response on the global scale. Our findings elucidate the role of activity and interactions in the cost repartition of collective behavior and explain what observed in experiments on natural living groups., Comment: 15 pages, 6 figures

Experiments on bird flocks and midge swarms reveal that these natural systems are well described by an active theory in which conservation laws play a crucial role. By building a symplectic structure that couples the particles' velocities to the generator of their internal rotations (spin), the Inertial Spin Model (ISM) reinstates a second-order temporal dynamics that captures many phenomenological traits of flocks and swarms. The reversible structure of the ISM predicts that the total spin is a constant of motion, the central conservation law responsible for all the novel dynamical features of the model. However, fluctuations and dissipation introduced in the original model to make it relax, violate the spin conservation law, so that the ISM aligns with the biophysical phenomenology only within finite-size regimes, beyond which the overdamped dynamics characteristic of the Vicsek model takes over. Here, we introduce a novel version of the ISM, in which the irreversible terms needed to relax the dynamics strictly respect the conservation of the spin. We perform a numerical investigation of the fully conservative model, exploring both the fixed-network case, which belongs to the equilibrium class of Model G, and the active case, characterized by self-propulsion of the agents and an out-of-equilibrium reshuffling of the underlying interaction network. Our simulations not only capture the correct spin wave phenomenology of the ordered phase, but they also yield dynamical critical exponents in the near-ordering phase that agree very well with the theoretical predictions., Comment: 14 pages, 6 figures

Numerical analysis of conserved field dynamics has been generally performed with pseudo spectral methods. Finite differences integration, the common procedure for non-conserved field dynamics, indeed struggles to implement a conservative noise in the discrete spatial domain. In this work, we present a novel method to generate a conservative noise in the finite differences framework, which works for any discrete topology and boundary conditions. We apply it to numerically solve the conserved Kardar-Parisi-Zhang (cKPZ) equation, widely used to describe surface roughening when the number of particles is conserved. Our numerical simulations recover the correct scaling exponents $\alpha$, $\beta$, and $z$ in $d=1$ and in $d=2$. To illustrate the potentiality of the method, we further consider the cKPZ equation on different kinds of non-standard lattices and on the random Euclidean graph. This is the first numerical study of conserved field dynamics on an irregular topology, paving the way to a broad spectrum of possible applications., Comment: Multiplicative noise case added, 6 figures

Living groups move in complex environments and are constantly subject to external stimuli, predatory attacks and disturbances. An efficient response to such perturbations is vital to maintain the group's coherence and cohesion. Perturbations are often local, i.e. they are initially perceived only by few individuals in the group, but can elicit a global response. This is the case of starling flocks, that can turn very quickly to evade predators. In this paper, we investigate the conditions under which a global change of direction can occur upon local perturbations. Using minimal models of self-propelled particles, we show that a collective directional response occurs on timescales that grow with the system size and it is, therefore, a finite-size effect. The larger the group is, the longer it will take to turn. We also show that global coherent turns can only take place if i) the mechanism for information propagation is efficient enough to transmit the local reaction undamped through the whole group; and if ii) motility is not too strong, to avoid that the perturbed individual leaves the group before the turn is complete. No compliance with such conditions results in the group's fragmentation or in a non-efficient response., Comment: 23 pages, 7 figures

A well-established numerical technique to study the dynamics of spin systems in which symmetries and conservation laws play an important role is to microcanonically integrate their reversible equations of motion, obtaining thermalization through initial conditions drawn with the canonical distribution. In order to achieve a more realistic relaxation of the magnetic energy, numerically expensive methods that explicitly couple the spins to the underlying lattice are normally employed. Here we introduce a stochastic conservative thermostat that relaxes the magnetic energy while preserving the constant of motions, thus turning microcanonical spin dynamics into a conservative canonical dynamics, without actually simulating the lattice. We test the thermostat on the Heisenberg antiferromagnet in d=3 and show that the method reproduces the exact values of the static and dynamic critical exponents, while in the low-temperature phase it yields the correct spin wave phenomenology. Finally, we demonstrate that the relaxation coefficient of the new thermostat is quantitatively connected to the microscopic parameters of the spin-lattice coupling., Comment: Revised version. 30 pages, 7 figures

Flocking in $d=2$ is a genuine non-equilibrium phenomenon for which irreversibility is an essential ingredient. We study a class of minimal flocking models whose only source of irreversibility is self-propulsion and use the entropy production rate (EPR) to quantify the departure from equilibrium across their phase diagrams. The EPR is maximal in the vicinity of the order-disorder transition, where reshuffling of the interaction network is fast. We show that signatures of irreversibility come in the form of asymmetries in the steady state distribution of the flock's microstates. They occur as consequences of the time reversal symmetry breaking in the considered self-propelled systems, independently of the interaction details. In the case of metric pairwise forces, they reduce to local asymmetries in the distribution of pairs of particles. This study suggests a possible use of pair asymmetries both to quantify the departure from equilibrium and to learn relevant information about aligning interaction potentials from data., Comment: 8 pages + Appendix; 6 figures