Traveling and resting crystals in active systems

Activity waves in condensed excitable phases of Quincke rollers

Traveling waves are universal in excitable systems; yet, the microscopic dynamics of wave propagation is inaccessible in conventional excitable systems. Here, we show that active colloids of Quincke rollers driven by a periodic electric field can form condensed excitable phases. Distinct from existing excitable media, condensed excitable colloids can be tuned reversibly between active liquids and active crystals in which two distinct waves can be excited, respectively. In active liquids, waves propagate by splitting and cross over each other, like sound waves, in collision. In active crystals, waves annihilate or converge, like shock waves, in collision. We show that the microscopic dynamics of sound waves is dominated by electrostatic repulsions while the dynamics of shock waves is encoded with a local density-dependent memory of propulsion. The condensed excitable colloids with tunable and controllable dynamics offer unexplored opportunities for the study of nonlinear phenomena.

Active smectics on a sphere

The dynamics of active smectic liquid crystals confined on a spherical surface is explored through an active phase field crystal model. Starting from an initially randomly perturbed isotropic phase, several types of topological defects are spontaneously formed, and then annihilate during a coarsening process until a steady state is achieved. The coarsening process is highly complex involving several scaling laws of defect densities as a function of time where different dynamical exponents can be identified. In general the exponent for the final stage towards the steady state is significantly larger than that in the passive and in the planar case, i.e. the coarsening is getting accelerated both by activity and by the topological and geometrical properties of the sphere. A defect type characteristic for this active system is a rotating spiral of evolving smectic layering lines. On a sphere this defect type also determines the steady state. Our results can in principle be confirmed by dense systems of synthetic or biological active particles.

Motility and swimming: universal description and generic trajectories

Autonomous locomotion is a ubiquitous phenomenon in biology and in physics of active systems at microscopic scale. This includes prokaryotic, eukaryotic cells (crawling and swimming) and artificial swimmers. An outstanding feature is the ability of these entities to follow complex trajectories, ranging from straight, curved (circular, helical...), to random-like ones. The non-straight nature of these trajectories is often explained as a consequence of the asymmetry of the particle or the medium in which it moves, or due to the presence of bounding walls, etc... Here, we show that for a particle driven by a concentration field of an active species, straight, circular and helical trajectories emerge naturally in the absence of asymmetry of the particle or that of suspending medium. Our proof is based on general considerations, without referring to an explicit form of a model. We show that these three trajectories correspond to self-congruent solutions. Self-congruency means that the states of the system at different moments of time can be made identical by an appropriate combination of rotation and translation of the coordinate space. We show that these solutions are exhibited by spherically symmetric particles as a result of a series of pitchfork bifurcations, leading to spontaneous symmetry breaking in the concentration field driving the particle motility. Self-congruent dynamics in one and two dimensions are analyzed as well. Finally, we present a simple explicit nonlinear exactly solvable model of fully isotropic phoretic particle that shows the transitions from a non-motile state to straight motion to circular motion to helical motion as a series of spontaneous symmetry-breaking bifurcations. Whether a system exhibits or not a given trajectory only depends on the numerical values of parameters entering the model, while asymmetry of swimmer shape, or anisotropy of the suspending medium, or influence of bounding walls are not necessary.

Time-reversal and parity-time symmetry breaking in non-Hermitian field theories

We study time-reversal symmetry breaking in non-Hermitian fluctuating field theories with conserved dynamics, comprising the mesoscopic descriptions of a wide range of nonequilibrium phenomena. They exhibit continuous parity-time (PT) symmetry-breaking phase transitions to dynamical phases. For two concrete transition scenarios, exclusive to non-Hermitian dynamics, namely, oscillatory instabilities and critical exceptional points, a low-noise expansion exposes a pretransitional surge of the mesoscale (informatic) entropy production rate, inside the static phases. Its scaling in the susceptibility contrasts conventional critical points (such as second-order phase transitions), where the susceptibility also diverges, but the entropy production generally remains finite. The difference can be attributed to active fluctuations in the wavelengths that become unstable. For critical exceptional points, we identify the coupling of eigenmodes as the entropy-generating mechanism, causing a drastic noise amplification in the Goldstone mode.

Entropy production in the nonreciprocal Cahn-Hilliard model

We study the nonreciprocal Cahn-Hilliard model with thermal noise as a prototypical example of a generic class of non-Hermitian stochastic field theories, analyzed in two companion papers [Suchanek, Kroy, and Loos, Phys. Rev. Lett. 131, 258302 (2023)10.1103/PhysRevLett.131.258302; Phys. Rev. E 108, 064123 (2023)10.1103/PhysRevE.108.064123]. Due to the nonreciprocal coupling between two field components, the model is inherently out of equilibrium and can be regarded as an active field theory. Beyond the conventional homogeneous and static-demixed phases, it exhibits a traveling-wave phase, which can be entered via either an oscillatory instability or a critical exceptional point. By means of a Fourier decomposition of the entropy production rate, we quantify the associated scale-resolved time-reversal symmetry breaking, in all phases and across the transitions, in the low-noise regime. Our perturbative calculation reveals its dependence on the strength of the nonreciprocal coupling. Surging entropy production near the static-dynamic transitions can be attributed to entropy-generating fluctuations in the longest wavelength Fourier mode and heralds the emerging traveling wave. Its translational dynamics can be mapped on the dissipative ballistic motion of an active (quasi)particle.

Entropy production and collective excitations of crystals out of equilibrium: The concept of entropons

We study the collective vibrational excitations of crystals under out-of-equilibrium steady conditions that give rise to entropy production. Their excitation spectrum comprises equilibriumlike phonons of thermal origin and additional collective excitations called entropons because each of them represents a mode of spectral entropy production. Entropons coexist with phonons and dominate them when the system is far from equilibrium while they are negligible in near-equilibrium regimes. The concept of entropons has been recently introduced and verified in a special case of crystals formed by self-propelled particles. Here we show that entropons exist in a broader class of active crystals that are intrinsically out of equilibrium and characterized by the lack of detailed balance. After a general derivation, several explicit examples are discussed, including crystals consisting of particles with alignment interactions and frictional contact forces.

Chemotactic Motility-Induced Phase Separation

Collectives of actively moving particles can spontaneously separate into dilute and dense phases-a fascinating phenomenon known as motility-induced phase separation (MIPS). MIPS is well-studied for randomly moving particles with no directional bias. However, many forms of active matter exhibit collective chemotaxis, directed motion along a chemical gradient that the constituent particles can generate themselves. Here, using theory and simulations, we demonstrate that collective chemotaxis strongly competes with MIPS-in some cases, arresting or completely suppressing phase separation, or in other cases, generating fundamentally new dynamic instabilities. We establish principles describing this competition, thereby helping to reveal and clarify the rich physics underlying active matter systems that perform chemotaxis, ranging from cells to robots.

Extreme Spontaneous Deformations of Active Crystals

We demonstrate that two-dimensional crystals made of active particles can experience extremely large spontaneous deformations without melting. Using particles mostly interacting via pairwise repulsive forces, we show that such active crystals maintain long-range bond order and algebraically decaying positional order, but with an exponent η not limited by the 1/3 bound given by the (equilibrium) KTHNY theory. We rationalize our findings using linear elastic theory and show the existence of two well-defined effective temperatures quantifying respectively large-scale deformations and bond-order fluctuations. The root of these phenomena lies in the sole time-persistence of the intrinsic axes of particles, and they should thus be observed in many different situations.

Entropons as collective excitations in active solids

The vibrational dynamics of solids is described by phonons constituting basic collective excitations in equilibrium crystals. Here, we consider a non-equilibrium active solid, formed by self-propelled particles, which bring the system into a non-equilibrium steady-state. We identify novel vibrational collective excitations of non-equilibrium (active) origin, which coexist with phonons and dominate over them when the system is far from equilibrium. These vibrational excitations are interpreted in the framework of non-equilibrium physics, in particular, stochastic thermodynamics. We call them "entropons" because they are the modes of spectral entropy production (at a given frequency and wave vector). The existence of entropons could be verified in future experiments on dense self-propelled colloidal Janus particles and granular active matter, as well as in living systems, such as dense cell monolayers.

Stationary broken parity states in active matter models

We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.

Non-reciprocity induces resonances in a two-field Cahn-Hilliard model

We consider a non-reciprocally coupled two-field Cahn-Hilliard system that has been shown to allow for oscillatory behaviour and suppression of coarsening. After introducing the model, we first review the linear stability of steady uniform states and show that all instability thresholds are identical to the ones for a corresponding two-species reaction-diffusion system. Next, we consider a specific interaction of linear modes-a 'Hopf-Turing' resonance-and derive the corresponding amplitude equations using a weakly nonlinear approach. We discuss the weakly nonlinear results and finally compare them with fully nonlinear simulations for a specific conserved amended FitzHugh-Nagumo system. We conclude with a discussion of the limitations of the employed weakly nonlinear approach. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

From a microscopic inertial active matter model to the Schrödinger equation

Active field theories, such as the paradigmatic model known as 'active model B+', are simple yet very powerful tools for describing phenomena such as motility-induced phase separation. No comparable theory has been derived yet for the underdamped case. In this work, we introduce active model I+, an extension of active model B+ to particles with inertia. The governing equations of active model I+ are systematically derived from the microscopic Langevin equations. We show that, for underdamped active particles, thermodynamic and mechanical definitions of the velocity field no longer coincide and that the density-dependent swimming speed plays the role of an effective viscosity. Moreover, active model I+ contains an analog of the Schrödinger equation in Madelung form as a limiting case, allowing one to find analoga of the quantum-mechanical tunnel effect and of fuzzy dark matter in active fluids. We investigate the active tunnel effect analytically and via numerical continuation.

Discontinuous Tension-Controlled Transition between Collective Actuations in Active Solids

The recent finding of collective actuation in active solids-solids embedded with active units-is a new promise for the design of multifunctional materials with genuine autonomy, and a better understanding of dense biological systems. Here, we combine the experimental study of centimetric model active solids, the numerical study of an agent-based model, and theoretical arguments to reveal a new form of collective actuation and how mechanical tension can serve as a general mechanism for transitioning between different collective actuation regimes. The presence of hysteresis when varying tension back and forth highlights the nontrivial selectivity of collective actuations.

Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

Spontaneous organization and phase separation of skyrmions in chiral active matter

Skyrmions are topologically protected vortex-like excitations that hold promise for applications such as information processing and electron manipulation. Here we combine theoretical analysis and numerical simulations to show that skyrmions can spontaneously emerge in chiral active matter without external confinements or regulation. Strikingly, these activity-driven skyrmions can either self-organize into a periodic, stable square lattice consisting of half Néel skyrmions and antiskyrmions, where the in-plane flows display an antiferromagnetic vortex array, or undergo phase separation between skyrmions with different topological numbers. We identify that the emerging skyrmion dynamics stems from the competition between the chiral and polar coherence length scales dictated by the interplay of intrinsic chirality, polarity, and elasticity in the system. Our results reveal unanticipated topological excitations, self-organization, and phase separation in non-equilibrium systems and also suggest a potential way towards engineering complicated bespoke skyrmionic structures through manipulating active matter.

Breaking action-reaction with active apolar colloids: emergent transport and velocity inversion

Artificial active particles are autonomous agents able to convert energy from the environment into net propulsion, breaking detailed balance and the action-reaction law, clear signatures of their out-of-equilibrium nature. Here we investigate the emergence of directed motion in clusters composed of passive and catalytically active apolar colloids. We use a light-induced chemophoretic flow to rapidly assemble hybrid self-propelling clusters composed of hematite particles and passive silica spheres. By increasing the size of the passive cargo, we observe a reversal in the transport direction of the pair. We explain this complex yet rich phenomenon using a theoretical model which accounts for the generated chemical field and its coupling with the surrounding medium. We exploit further our technique to build up more complex, chemically driven, architectures capable of carrying several passive or active species, that quickly assemble and disassemble under light control.

Active nonreciprocal attraction between motile particles in an elastic medium

We show from experiments and simulations on vibration-activated granular matter that self-propelled polar rods in an elastic medium on a substrate turn and move towards each other. We account for this effective attraction through a coarse-grained theory of a motile particle as a moving point-force density that creates elastic strains in the medium that reorient other particles. Our measurements confirm qualitatively the predicted features of the distortions created by the rods, including the |x|^{-1/2} tail of the trailing displacement field and nonreciprocal sensing and pursuit. A discrepancy between the magnitudes of displacements along and transverse to the direction of motion remains. Our theory should be of relevance to the interaction of motile cells in the extracellular matrix or in a supported layer of gel or tissue.

Repulsive torques alone trigger crystallization of constant speed active particles

We investigate the possibility for self-propelled particles to crystallize without reducing their intrinsic speed. We illuminate how, in the absence of any force, the competition between self-propulsion and repulsive torques determines the macroscopic phases of constant-speed active particles. This minimal model expands upon existing approaches for an improved understanding of crystallization of active matter.

Deformable active nematic particles and emerging edge currents in circular confinements

We consider a microscopic field theoretical approach for interacting active nematic particles. With only steric interactions the self-propulsion strength in such systems can lead to different collective behaviour, e.g. synchronized self-spinning and collective translation. The different behaviour results from the delicate interplay between internal nematic structure, particle shape deformation and particle-particle interaction. For intermediate active strength an asymmetric particle shape emerges and leads to chirality and self-spinning crystals. For larger active strength the shape is symmetric and translational collective motion emerges. Within circular confinements, depending on the packing fraction, the self-spinning regime either stabilizes positional and orientational order or can lead to edge currents and global rotation which destroys the synchronized self-spinning crystalline structure.

Morphologies and dynamics of the interfaces between active and passive phases

Active matters exhibit interesting collective behaviors and novel phases, which provide an important platform for the study of nonequilibrium physics. Mixtures of active and passive particles have been intensively studied in motility-induced phase separation, but the morphology of the active-passive interface has been poorly explored. In this work, we investigate the interface morphology in two-dimensional mixtures of active and passive particles using Brownian dynamics simulations. By systematically changing the Péclet number (Pe) and area fraction (ρ), we obtain the phase diagram of the active-passive interface, including rough sharp, rough invasive and flat interdiffusive interfaces. For a sharp interface, dynamic scaling analysis in the propagation stage shows that the roughness exponent α, the growth exponent β, the time exponent κ, and the dynamic exponent z satisfy z = α/(β - κ). Such anomalous scaling indicates that the roughening behavior does not belong to the conventional universality classes with Family-Vicsek scaling for the growth of passive interfaces. On the other hand, the interface in the middle-wavelength regime during the morphology relaxation stage can be described by capillary wave theory. The mean interface position propagates with time as t^1/2, which is robust at different ρ and Pe values in the propagation stage and exhibits superdiffusion in the morphology relaxation stage. These similarities and differences between the active-inactive interfaces and passive interfaces cast light on the interfacial growth of active matter.

Emergence and melting of active vortex crystals

Melting of two-dimensional (2D) equilibrium crystals is a complex phenomenon characterized by the sequential loss of positional and orientational order. In contrast to passive systems, active crystals can self-assemble and melt into an active fluid by virtue of their intrinsic motility and inherent non-equilibrium stresses. Currently, the non-equilibrium physics of active crystallization and melting processes is not well understood. Here, we establish the emergence and investigate the melting of self-organized vortex crystals in 2D active fluids using a generalized Toner-Tu theory. Performing extensive hydrodynamic simulations, we find rich transition scenarios. On small domains, we identify a hysteretic transition as well as a transition featuring temporal coexistence of active vortex lattices and active turbulence, both of which can be controlled by self-propulsion and active stresses. On large domains, an active vortex crystal with solid order forms within the parameter range corresponding to active vortex lattices. The melting of this crystal proceeds through an intermediate hexatic phase. Generally, these results highlight the differences and similarities between crystalline phases in active fluids and their equilibrium counterparts.

Spatial velocity correlations in inertial systems of active Brownian particles

Recently, it has been discovered that systems of active Brownian particles (APB) at high density organise their velocities into coherent domains showing large spatial structures in the velocity field. This collective behavior occurs spontaneously, i.e. is not caused by any specific interparticle force favoring the alignment of the velocities. This phenomenon was investigated in the absence of thermal noise and in the overdamped regime where inertial forces could be neglected. In this work, we demonstrate through numerical simulations and theoretical analysis that velocity alignment is a robust property of ABP and persists even in the presence of inertial forces and thermal fluctuations. We also show that a single dimensionless parameter, such as the Péclet number customarily employed in the description of self-propelled particles, is not sufficient to fully characterize this phenomenon either in the regimes of large viscosity or small mass. Indeed, the size of the velocity domains, measured through the correlation length of the spatial velocity correlation, remains constant when the swim velocity increases and decreases as the rotational diffusion becomes larger. We find that, contrary to the common belief, the spatial velocity correlation not only depends on inertia but is also non-symmetrically affected by mass and inverse viscosity variations. We conclude that in self-propelled systems, at variance with passive systems, variations in the inertial time (mass over solvent viscosity) and mass act as independent control parameters. Finally, we highlight the non-thermal nature of the spatial velocity correlations that are fairly insensitive both to solvent and active temperatures.

Two-dimensional localized states in an active phase-field-crystal model

The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations, and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.

Phase-field-crystal description of active crystallites: Elastic and inelastic collisions

The active Phase-Field-Crystal (aPFC) model combines elements of the Toner-Tu theory for self-propelled particles and the classical Phase-Field-Crystal (PFC) model that describes the transition between liquid and crystalline phases. In the liquid-crystal coexistence region of the PFC model, crystalline clusters exist in the form of localized states that coexist with a homogeneous background. At sufficiently strong activity (related to self-propulsion strength), they start to travel. We employ numerical path continuation and direct time simulations to first investigate the existence regions of different types of localized states in one spatial dimension. The results are summarized in morphological phase diagrams in the parameter plane spanned by activity and mean density. Then we focus on the interaction of traveling localized states, studying their collision behavior. As a result, we distinguish "elastic" and "inelastic" collisions. In the former, localized states recover their properties after a collision, while in the latter, they may completely or partially annihilate, forming resting bound states or various traveling states.

Dynamical Crystallites of Active Chiral Particles

One of the intrinsic characteristics of far-from-equilibrium systems is the nonrelaxational nature of the system dynamics, which leads to novel properties that cannot be understood and described by conventional pathways based on thermodynamic potentials. Of particular interest are the formation and evolution of ordered patterns composed of active particles that exhibit collective behavior. Here we examine such a type of nonpotential active system, focusing on effects of coupling and competition between chiral particle self-propulsion and self-spinning. It leads to the transition between three bulk dynamical regimes dominated by collective translative motion, spinning-induced structural arrest, and dynamical frustration. In addition, a persistently dynamical state of self-rotating crystallites is identified as a result of a localized-delocalized transition induced by the crystal-melt interface. The mechanism for the breaking of localized bulk states can also be utilized to achieve self-shearing or self-flow of active crystalline layers.

Active phase field crystal systems with inertial delay and underdamped dynamics

Active matter systems often are well approximated as overdamped, meaning that any inertial momentum is immediately dissipated by the environment. On the other hand, especially for macroscopic systems but also for many mesoscopic ones particle mass can become relevant for the dynamics. For such systems we recently proposed an underdamped continuum model which captures translationally inertial dynamics via two contributions. First, convection and second a damping time scale of inertial motion. In this paper, we ask how both of these features influence the collective behavior compared to overdamped dynamics by studying the example of the active phase field crystal model. We first focus on the case of suppressed convection to study the role of the damping time. We quantify that the relaxation process to the steady collective motion state is considerably prolonged with damping time due to the increasing occurrence of transient groups of circularly moving density peaks. Finally, we illustrate the fully underdamped case with convection. Instead of collective motion of density peaks we then find a coexistence of constant high and low density phases reminiscent of motility-induced phase separation.

Predictive local field theory for interacting active Brownian spheres in two spatial dimensions

We present a predictive local field theory for the nonequilibrium dynamics of interacting active Brownian particles with a spherical shape in two spatial dimensions. The theory is derived by a rigorous coarse-graining starting from the Langevin equations that describe the trajectories of the individual particles. For high accuracy and generality of the theory, it includes configurational order parameters and derivatives up to infinite order. In addition, we discuss possible approximations of the theory and present reduced models that are easier to apply. We show that our theory contains popular models such as Active Model B+  as special cases and that it provides explicit expressions for the coefficients occurring in these and other, often phenomenological, models. As a further outcome, the theory yields an analytical expression for the density-dependent mean swimming speed of the particles. To demonstrate an application of the new theory, we analyze a simple reduced model of the lowest nontrivial order in derivatives, which is able to predict the onset of motility-induced phase separation of the particles. By a linear stability analysis, an analytical expression for the spinodal corresponding to motility-induced phase separation is obtained. This expression is evaluated for the case of particles interacting repulsively by a Weeks-Chandler-Andersen potential. The analytical predictions for the spinodal associated with these particles are found to be in very good agreement with the results of Brownian dynamics simulations that are based on the same Langevin equations as our theory. Furthermore, the critical point predicted by our analytical results agrees excellently with recent computational results from the literature.

Bifurcations of front motion in passive and active Allen-Cahn-type equations

The well-known cubic Allen-Cahn (AC) equation is a simple gradient dynamics (or variational) model for a nonconserved order parameter field. After revising main literature results for the occurrence of different types of moving fronts, we employ path continuation to determine their bifurcation diagram in dependence of the external field strength or chemical potential. We then employ the same methodology to systematically analyze fronts for more involved AC-type models. In particular, we consider a cubic-quintic variational AC model and two different nonvariational generalizations. We determine and compare the bifurcation diagrams of front solutions in the four considered models.

Oriented Active Solids

We present a complete analysis of the linearized dynamics of active solids with uniaxial orientational order, taking into account a hitherto overlooked consequence of rotation invariance. Our predictions include a purely active response of two-dimensional orientationally ordered solids to shear, the possibility of stable active solids with quasi-long-range order in two dimensions and long-range order in three dimensions, generic instability of the solid for one sign of active forcing, and the instability of the uniaxially ordered phase in momentum-conserved systems for large active forcing irrespective of its sign.

A comparative study between two models of active cluster crystals

We study a system of active particles with soft repulsive interactions that lead to an active cluster-crystal phase in two dimensions. We use two different modelizations of the active force - Active Brownian particles (ABP) and Ornstein-Uhlenbeck particles (AOUP) - and focus on analogies and differences between them. We study the different phases appearing in the system, in particular, the formation of ordered patterns drifting in space without being altered. We develop an effective description which captures some properties of the stable clusters for both ABP and AOUP. As an additional point, we confine such a system in a large channel, in order to study the interplay between the cluster crystal phase and the well-known accumulation near the walls, a phenomenology typical of active particles. For small activities, we find clusters attached to the walls and deformed, while for large values of the active force they collapse in stripes parallel to the walls.

Traveling band formation in feedback-driven colloids

Using simulation and theory we study the dynamics of a colloidal suspension in two dimensions subject to a time-delayed repulsive feedback that depends on the positions of the colloidal particles. The colloidal particles experience an additional potential that is a superposition of repulsive potential energies centered around the positions of all the particles a delay time ago. Here we show that such a feedback leads to self-organization of the particles into traveling bands. The width of the bands and their propagation speed can be tuned by the delay time and the range of the imposed repulsive potential. The emerging traveling band behavior is observed in Brownian dynamics computer simulations as well as microscopic dynamic density functional theory. Traveling band formation also persists in systems of finite size leading to rotating traveling waves in the case of circularly confined systems.

Two types of exclusion interactions for self-propelled objects and collective motion induced by their combination

Exclusive interactions between self-driven objects may play crucial roles in their collective behavior, e.g., in collective migration of living cells. Here, such collective behavior is studied based on a simple but sufficient model taking account the exclusion effects, which incorporate the following two distinct kinds of exclusion interactions in two dimensions: The first is the mechanical exclusion wherein two objects mechanically repel each other when they overlap. The second is the scattering exclusion, wherein the directions along which each object tries to move are modulated to avoid overlapping. We propose a theoretical model based on two principles: (1) Each object maintains its own polarity with a fixed strength and attempts to move into the polarity direction and (2) objects interact with each other through the abovementioned exclusions. Based on this model, we look at the difference of consequences and combinatory effects of these two kinds of exclusions. Furthermore, we calculate the polar order of polarity directions without an external directional bias. Our results suggest that the combination of these two kinds of exclusions leads to effectively inelastic scattering of two objects, which eventually gives rise to global polar ordering. We also find that the traveling band can arise by this mechanism of alignment at the intermediate density, as generally seen in collective motion with polar alignment and investigated in various earlier works. Characteristics of transitions among disordered, traveling band, and homogeneously ordered states of the presented model are investigated, and their similarities and differences with those given by the explicit alignment interaction are discussed.

Thermodynamic phases in two-dimensional active matter

Active matter has been much studied for its intriguing properties such as collective motion, motility-induced phase separation and giant fluctuations. However, it has remained unclear how the states of active materials connect with the equilibrium phases. For two-dimensional systems, this is also because the understanding of the liquid, hexatic, and solid equilibrium phases and their phase transitions is recent. Here we show that two-dimensional self-propelled point particles with inverse-power-law repulsions moving with a kinetic Monte Carlo algorithm without alignment interactions preserve all equilibrium phases up to very large activities. Furthermore, at high activity within the liquid phase, a critical point opens up a gas–liquid motility-induced phase separation region. In our model, two-step melting and motility-induced phase separation are thus independent phenomena. We discuss the reasons for these findings to be common to a wide class of two-dimensional active systems.

Adapting statistical physics tools to study active systems is challenging due to their non-equilibrium nature. Here the authors use simulations to present a phase diagram of a 2D active system, showing a two-step melting scenario far from equilibrium along with gas-liquid motility-induced phase separation.

Resting and traveling localized states in an active phase-field-crystal model

The conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) provides a simple microscopic description of the thermodynamic transition between fluid and crystalline states. Combining it with elements of the Toner-Tu theory for self-propelled particles, Menzel and Löwen [Phys. Rev. Lett. 110, 055702 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.055702] obtained a model for crystallization (swarm formation) in active systems. Here, we study the occurrence of resting and traveling localized states, i.e., crystalline clusters, within the resulting active PFC model. Based on linear stability analyses and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of periodic and localized, resting and traveling states in a one-dimensional active PFC model. This allows us, for instance, to explore how the slanted homoclinic snaking of steady localized states found for the passive PFC model is amended by activity. A particular focus lies on the onset of motion, where we show that it occurs either through a drift-pitchfork or a drift-transcritical bifurcation. A corresponding general analytical criterion is derived.

Spontaneously Flowing Crystal of Self-Propelled Particles

We experimentally and numerically study the structure and dynamics of a monodisperse packing of spontaneously aligning self-propelled hard disks. The packings are such that their equilibrium counterparts form perfectly ordered hexagonal structures. Experimentally, we first form a perfect crystal in a hexagonal arena which respects the same crystalline symmetry. Frustration of the hexagonal order, obtained by removing a few particles, leads to the formation of a rapidly diffusing "droplet." Removing more particles, the whole system spontaneously forms a macroscopic sheared flow, while conserving an overall crystalline structure. This flowing crystalline structure, which we call a "rheocrystal," is made possible by the condensation of shear along localized stacking faults. Numerical simulations very well reproduce the experimental observations and allow us to explore the parameter space. They demonstrate that the rheocrystal is induced neither by frustration nor by noise. They further show that larger systems flow faster while still remaining ordered.

Flocking transitions in confluent tissues

Collective cell migration in dense tissues underlies important biological processes, such as embryonic development, wound healing and cancer invasion. While many aspects of single cell movements are now well established, the mechanisms leading to displacements of cohesive cell groups are still poorly understood. To elucidate the emergence of collective migration in mechanosensitive cells, we examine a self-propelled Voronoi (SPV) model of confluent tissues with an orientational feedback that aligns a cell's polarization with its local migration velocity. While shape and motility are known to regulate a density-independent liquid-solid transition in tissues, we find that aligning interactions facilitate collective motion and promote solidification, with transitions that can be predicted by extending statistical physics tools such as effective temperature to this far-from-equilibrium system. In addition to accounting for recent experimental observations obtained with epithelial monolayers, our model predicts structural and dynamical signatures of flocking, which may serve as gateway to a more quantitative characterization of collective motility.

Active crystals on a sphere

Two-dimensional crystals on curved manifolds exhibit nontrivial defect structures. Here we consider "active crystals" on a sphere, which are composed of self-propelled colloidal particles. Our work is based on a phase-field-crystal-type model that involves a density and a polarization field on the sphere. Depending on the strength of the self-propulsion, three different types of crystals are found: a static crystal, a self-spinning "vortex-vortex" crystal containing two vortical poles of the local velocity, and a self-translating "source-sink" crystal with a source pole where crystallization occurs and a sink pole where the active crystal melts. These different crystalline states as well as their defects are studied theoretically here and can in principle be confirmed in experiments.

Aggregation-fragmentation and individual dynamics of active clusters

A remarkable feature of active matter is the propensity to self-organize. One striking instance of this ability to generate spatial structures is the cluster phase, where clusters broadly distributed in size constantly move and evolve through particle exchange, breaking or merging. Here we propose an exhaustive description of the cluster dynamics in apolar active matter. Exploiting large statistics gathered on thousands of Janus colloids, we measure the aggregation and fragmentation rates and rationalize the resulting cluster size distribution and fluctuations. We also show that the motion of individual clusters is entirely consistent with a model positing random orientation of colloids. Our findings establish a simple, generic model of cluster phase, and pave the way for a thorough understanding of clustering in active matter.

Flow of colloidal solids and fluids through constrictions: dynamical density functional theory versus simulation

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

Stabilization of active matter by flow-vortex lattices and defect ordering

Active systems, from bacterial suspensions to cellular monolayers, are continuously driven out of equilibrium by local injection of energy from their constituent elements and exhibit turbulent-like and chaotic patterns. Here we demonstrate both theoretically and through numerical simulations, that the crossover between wet active systems, whose behaviour is dominated by hydrodynamics, and dry active matter where any flow is screened, can be achieved by using friction as a control parameter. Moreover, we discover unexpected vortex ordering at this wet-dry crossover. We show that the self organization of vortices into lattices is accompanied by the spatial ordering of topological defects leading to active crystal-like structures. The emergence of vortex lattices, which leads to the positional ordering of topological defects, suggests potential applications in the design and control of active materials.

Getting drowned in a swirl: Deformable bead-spring model microswimmers in external flow fields

Deformability is a central feature of many types of microswimmers, e.g., for artificially generated self-propelled droplets. Here, we analyze deformable bead-spring microswimmers in an externally imposed solvent flow field as simple theoretical model systems. We focus on their behavior in a circular swirl flow in two spatial dimensions. Linear (straight) two-bead swimmers are found to circle around the swirl with a slight drift to the outside with increasing activity. In contrast to that, we observe for triangular three-bead or squarelike four-bead swimmers a tendency of being drawn into the swirl and finally getting drowned, although a radial inward component is absent in the flow field. During one cycle around the swirl, the self-propulsion direction of an active triangular or squarelike swimmer remains almost constant, while their orbits become deformed exhibiting an "egglike" shape. Over time, the swirl flow induces slight net rotations of these swimmer types, which leads to net rotations of the egg-shaped orbits. Interestingly, in certain cases, the orbital rotation changes sense when the swimmer approaches the flow singularity. Our predictions can be verified in real-space experiments on artificial microswimmers.

Effect of self-propulsion on equilibrium clustering

In equilibrium, colloidal suspensions governed by short-range attractive and long-range repulsive interactions form thermodynamically stable clusters. Using Brownian dynamics computer simulations, we investigate how this equilibrium clustering is affected when such particles are self-propelled. We find that the clustering process is stable under self-propulsion. For the range of interaction parameters studied and at low particle density, the cluster size increases with the speed of self-propulsion (activity) and for higher activity the cluster size decreases, showing a nonmonotonic variation of cluster size with activity. This clustering behavior is distinct from the pure kinetic (or motility-induced) clustering of self-propelling particles which is observed at significantly higher activities and densities. We present an equilibrium model incorporating the effect of activity as activity-induced attraction and repulsion by imposing that the strength of these interactions depend on activity superlinearly. The model explains the cluster size dependence of activity obtained from simulations semiquantitatively. Our predictions are verifiable in experiments on interacting synthetic colloidal microswimmers.