Disorder-Induced Long-Ranged Correlations in Scalar Active Matter

Motility-Induced Pinning in Flocking System with Discrete Symmetry

We report a motility-induced pinning transition in the active Ising model for a self-propelled particle system with discrete symmetry. This model was known to exhibit a liquid-gas type flocking phase transition, but a recent study reveals that the polar order is metastable due to droplet excitation. Using extensive Monte Carlo simulations, we demonstrate that, for an intermediate alignment interaction strength, the steady state is characterized by traveling local domains, which renders the polar order short-ranged in both space and time. We further demonstrate that interfaces between colliding domains become pinned as the alignment interaction strength increases. A resonating back-and-forth motion of individual self-propelled particles across interfaces is identified as a mechanism for the pinning. We present a numerical phase diagram for the motility-induced pinning transition, and an approximate analytic theory for the growth and shrink dynamics of pinned interfaces. Our results show that pinned interfaces grow to a macroscopic size preventing the polar order in the regime where the particle diffusion rate is sufficiently smaller than the self-propulsion rate. The growth behavior in the opposite regime and its implications on the polar order remain unresolved and require further investigation.

Continuous symmetry breaking of low-dimensional systems driven by inhomogeneous oscillatory driving forces

The driving forces of chiral active particles and deformations of cells are often modeled by spatially inhomogeneous but temporally periodic driving forces. Such inhomogeneous oscillatory driving forces have only recently been proposed in the context of active matter, and their effects on the systems are not yet fully understood. In this work, we theoretically study the impact of spatially inhomogeneous oscillatory driving forces on continuous symmetry breaking. We first analyze the linear model for the soft modes in the ordered phase to derive the lower critical dimension of the model, and then analyze the spherical model to investigate more detailed phase behaviors. Interestingly, our analysis reveals that symmetry breaking occurs even in one and two dimensions, where the Hohenberg-Mermin-Wagner theorem prohibits continuous symmetry breaking in equilibrium. Furthermore, fluctuations of conserved quantities, such as density, are anomalously suppressed in the long-wavelength, i.e., show hyperuniformity.

Synchronized rotations of active particles on chemical substrates

Many microorganisms use chemical 'signaling' - a quintessential self-organizing strategy in non-equilibrium - that can induce spontaneous aggregation and coordinated motion. Using synthetic signaling as a design principle, we construct a minimal model of active Brownian particles (ABPs) having soft repulsive interactions on a chemically quenched patterned substrate. The interplay between chemo-phoretic interactions and activity is numerically investigated for a proposed variant of the Keller-Segel model for chemotaxis. Such competition not only results in a chemo-motility-induced phase-separated state, but also results in a new cohesive clustering phase with synchronized rotations. Our results suggest that rotational order can emerge in systems by virtue of activity and repulsive interactions alone without an explicit alignment interaction. These rotations can also be exploited by designing mechanical devices that can generate reorienting torques using active particles.

Stationary particle currents in sedimenting active matter wetting a wall

Recently it was predicted, on the basis of a lattice gas model, that scalar active matter in a gravitational field would rise against gravity up a confining wall or inside a thin capillary-in spite of repulsive particle-wall interactions [Phys. Rev. Lett. 124, 048001 (2020)0031-900710.1103/PhysRevLett.124.048001]. In this paper we confirm this prediction with sedimenting active Brownian particles (ABPs) in a box numerically and elucidate the mechanism leading to the formation of a meniscus rising above the bulk of the sedimentation region. The height of the meniscus increases with the activity of the system, algebraically with the Péclet number. The formation of the meniscus is determined by a stationary circular particle current, a vortex, centered at the base of the meniscus, whose size and strength increase with the ABP activity. The origin of these vortices can be traced back to the confinement of the ABPs in a box: already the stationary state of ideal (noninteracting) ABPs without gravitation displays circular currents that arrange in a highly symmetric way in the eight octants of the box. Gravitation distorts this vortex configuration downward, leaving two major vortices at the two side walls, with a strong downward flow along the walls. Repulsive interactions between the ABPs change this situation only as soon as motility induced phase separation (MIPS) sets in and forms a dense, sedimented liquid region at the bottom, which pushes the center of the vortex upwards towards the liquid-gas interface. Self-propelled particles therefore represent an impressive realization of scalar active matter that forms stationary particle currents being able to perform visible work against gravity or any other external field, which we predict to be observable experimentally in active colloids under gravitation.

Detecting and characterizing phase transitions in active matter using entropy

A major challenge in the study of active matter lies in quantitative characterization of phases and transitions between them. We show how the entropy of a collection of active objects can be used to classify regimes and spatial patterns in their collective behavior. Specifically, we estimate the contributions to the total entropy from correlations between the degrees of freedom of position and orientation. This analysis pin-points the flocking transition in the Vicsek model while clarifying the physical mechanism behind the transition. When applied to experiments on swarming Bacillus subtilis with different cell aspect ratios and overall bacterial area fractions, the entropy analysis reveals a rich phase diagram with transitions between qualitatively different swarm statistics. We discuss physical and biological implications of these findings.

Resolving entropy contributions in nonequilibrium transitions

We derive a functional for the entropy contributed by any microscopic degrees of freedom as arising from their measurable pair correlations. Applicable both in and out of equilibrium, this functional yields the maximum entropy which a system can have given a certain correlation function. When applied to different correlations, the method allows us to identify the degrees of freedom governing a certain physical regime, thus capturing and characterizing dynamic transitions. The formalism applies also to systems whose translational invariance is broken by external forces and whose number of particles may vary. We apply it to experimental results for jammed bidisperse emulsions, capturing the crossover of this nonequilibrium system from crystalline to disordered hyperuniform structures as a function of mixture composition. We discover that the cross-correlations between the positions and sizes of droplets in the emulsion play the central role in the formation of the disordered hyperuniform states. We discuss implications of the approach for entropy estimation out of equilibrium and for characterizing transitions in disordered systems.

Morphologies and dynamics of free surfaces of crystals composed of active particles

Active matter exhibits various collective motions and nonequilibrium phases, such as crystals; however, their surface properties have been poorly explored. Here, we use Brownian dynamics simulations to investigate the surface morphology and dynamics of two-dimensional active crystals during and after growth. For crystal growth on a substrate, the position and roughness of the crystal surface reach steady states at different times. In the steady state, the surface exhibits superdiffusive behaviour at the short time, and the roughness is insensitive to the roughening process and particle activity. We observe two-stage and three-stage surface roughening at different Péclet numbers. The result of dynamic scaling analysis shows that the surface is similar to anomalous roughening, which is distinct from the normal roughening typically found in conventional passive systems. Capillary wave theory for a thermal equilibrium system can describe the active surface fluctuations only in the long-wavelength regime, indicating that active particles mainly drive the surface out of equilibrium locally. These similarities and differences between the active and passive crystal surfaces are essential for understanding active crystals and interfaces.

Incompressible Polar Active Fluids with Quenched Random Field Disorder in Dimensions d>2

We present a hydrodynamic theory of incompressible polar active fluids with quenched random field disorder. This theory shows that such fluids can overcome the disruption caused by the quenched disorder and move coherently, in the sense of having a nonzero mean velocity in the hydrodynamic limit. However, the scaling behavior of this class of active systems cannot be described by linearized hydrodynamics in spatial dimensions between 2 and 5. Nonetheless, we obtain the exact dimension-dependent scaling exponents in these dimensions.

Packed Swarms on Dirt: Two-Dimensional Incompressible Flocks with Quenched and Annealed Disorder

We show that incompressible polar active fluids can exhibit an ordered, coherently moving phase even in the presence of quenched disorder in two dimensions. Unlike such active fluids with annealed disorder (i.e., time-dependent random white noise) only, which behave like equilibrium ferromagnets with long-range interactions, this robustness against quenched disorder is a fundamentally nonequilibrium phenomenon. The ordered state belongs to a new universality class, whose scaling laws we calculate using three different renormalization group schemes, which all give scaling exponents within 0.02 of each other, indicating that our results are quite accurate. Our predictions can be quantitatively tested in readily available artificial active systems and imply that biological systems such as cell layers can move coherently in vivo, where disorder is inevitable.

Passive objects in confined active fluids: A localization transition

We study how walls confining active fluids interact with asymmetric passive objects placed in their bulk. We show that the objects experience nonconservative long-ranged forces mediated by the active bath. To leading order, these forces can be computed using a generalized image theorem. The walls repel asymmetric objects, irrespective of their microscopic properties or their orientations. For circular cavities, we demonstrate how this may lead to the localization of asymmetric objects in the center of the cavity, something impossible for symmetric ones.

Hydrodynamic theory of two-dimensional incompressible polar active fluids with quenched and annealed disorder

We study the moving phase of two-dimensional (2D) incompressible polar active fluids in the presence of both quenched and annealed disorder. We show that long-range polar order persists even in this defect-ridden two-dimensional system. We obtain the large-distance, long-time scaling laws of the velocity fluctuations using three distinct dynamic renormalization group schemes. These are an uncontrolled one-loop calculation in exactly two dimensions, and two d=(d_{c}-ε) expansions to O(ε), obtained by two different analytic continuations of our 2D model to higher spatial dimensions: a "hard" continuation which has d_{c}=7/3, and a "soft" continuation with d_{c}=5/2. Surprisingly, the quenched and annealed parts of the velocity correlation function have the same anisotropy exponent and the relaxational and propagating parts of the dispersion relation have the same dynamic exponent in the nonlinear theory even though they are distinct in the linearized theory. This is due to anomalous hydrodynamics. Furthermore, all three renormalization schemes yield very similar values for the universal exponents, and therefore we expect the numerical values that we predict for them to be highly accurate.

Dense polar active fluids in a disordered environment

We examine the influence of quenched disorder on the flocking transition of dense polar active matter. We consider incompressible systems of active particles with aligning interactions under the effect of either quenched random forces or random dilution. The system displays a continuous disorder-order (flocking) transition, and the associated scaling behavior is described by a new universality class which is controlled by a quenched Navier-Stokes fixed point. We determine the critical exponents through a perturbative renormalization group analysis. We show that the two forms of quenched disorder, random force and random mass (dilution), belong to the same universality class, in contrast with the situation at equilibrium.

Susceptibility of Polar Flocks to Spatial Anisotropy

We study the effect of spatial anisotropy on polar flocks by investigating active q-state clock models in two dimensions. In contrast to the equilibrium case, we find that any amount of anisotropy is asymptotically relevant, drastically altering the phenomenology from that of the rotationally invariant case. All of the well-known physics of the Vicsek model, from giant density fluctuations to microphase separation, is replaced by that of the active Ising model, with short-range correlations and complete phase separation. These changes appear beyond a length scale that diverges in the q→∞ limit, so that the Vicsek-model phenomenology is observed in finite systems for weak enough anisotropy, i.e., sufficiently high q. We provide a scaling argument which explains why anisotropy has such different effects in the passive and active cases.

Scaling of entropy production under coarse graining in active disordered media

Entropy production plays a fundamental role in the study of nonequilibrium systems by offering a quantitative handle on the degree of time-reversal symmetry breaking. It depends crucially on the degree of freedom considered as well as on the scale of description. How the entropy production at one resolution of the degrees of freedom is related to the entropy production at another resolution is a fundamental question which has recently attracted interest. This relationship is of particular relevance to coarse-grained and continuum descriptions of a given phenomenon. In this work, we derive the scaling of the entropy production under iterative coarse graining on the basis of the correlations of the underlying microscopic transition rates for noninteracting particles in active disordered media. Our approach unveils a natural criterion to distinguish equilibrium-like and genuinely nonequilibrium macroscopic phenomena based on the sign of the scaling exponent of the entropy production per mesostate.

Confinement and Collective Escape of Active Particles

Active matter broadly covers the dynamics of self-propelled particles. While the onset of collective behavior in homogenous active systems is relatively well understood, the effect of inhomogeneities such as obstacles and traps lacks overall clarity. Here, we study how interacting, self-propelled particles become trapped and released from a trap. We have found that captured particles aggregate into an orbiting condensate with a crystalline structure. As more particles are added, the trapped condensates escape as a whole. Our results shed light on the effects of confinement and quenched disorder in active matter.

Transport and phase separation of active Brownian particles in fluctuating environments

In this work, we study the dynamics of a single active Brownian particle, as well as the collective behavior of interacting active Brownian particles, in a fluctuating heterogeneous environment. We employ a variant of the diffusing diffusivity model where the equation of motion of the active particle involves a time-dependent motility and diffusivities. Within our model, those fluctuations are coupled to each other. Using analytical methods, we obtain the probability distribution function of particle displacement and its moments for a single particle. We then investigate the impact of the environmental fluctuations on the collective behavior of the active Brownian particles by means of extensive numerical simulations. Our results show that the fluctuations hinder the motility-induced phase separation, accompanied by a significant change of the density dependence of particle velocities. These effects are interpreted using our analytical results for the dynamics of a single particle.

Susceptibility of Orientationally Ordered Active Matter to Chirality Disorder

We investigate the susceptibility of long-range ordered phases of two-dimensional dry aligning active matter to population disorder, taken in the form of a distribution of intrinsic individual chiralities. Using a combination of particle-level models and hydrodynamic theories derived from them, we show that while in finite systems all ordered phases resist a finite amount of such chirality disorder, the homogeneous ones (polar flocks and active nematics) are unstable to any amount of disorder in the infinite-size limit. On the other hand, we find that the inhomogeneous solutions of the coexistence phase (bands) may resist a finite amount of chirality disorder even asymptotically.

Kinetic Monte Carlo Algorithms for Active Matter Systems

We study kinetic Monte Carlo (KMC) descriptions of active particles. We show that, when they rely on purely persistent, active steps, their continuous-time limit is ill-defined, leading to the vanishing of trademark behaviors of active matter such as the motility-induced phase separation, ratchet effects, as well as to a diverging mechanical pressure. We then show how, under an appropriate scaling, mixing passive steps with active ones leads to a well-defined continuous-time limit that however differs from standard active dynamics. Finally, we propose new KMC algorithms whose continuous-time limits lead to the dynamics of active Ornstein-Uhlenbeck, active Brownian, and run-and-tumble particles.

Breakdown of Ergodicity and Self-Averaging in Polar Flocks with Quenched Disorder

We show that spatial quenched disorder affects polar active matter in ways more complex and far reaching than heretofore believed. Using simulations of the 2D Vicsek model subjected to random couplings or a disordered scattering field, we find in particular that ergodicity is lost in the ordered phase, the nature of which we show to depend qualitatively on the type of quenched disorder: for random couplings, it remains long-range ordered, but qualitatively different from the pure (disorderless) case. For random scatterers, polar order varies with system size but we find strong non-self-averaging, with sample-to-sample fluctuations dominating asymptotically, which prevents us from elucidating the asymptotic status of order.