Mode coupling theory for nonequilibrium glassy dynamics of thermal self-propelled particles

An elastoplastic model approach for the relaxation dynamics of active glasses

How activity affects the glassy dynamics is crucial for several biological processes. Furthermore, active glasses offer fascinating phenomenologies, extend the scope of equilibrium glass-forming liquids, and can provide novel insights into the original problem. We introduce a family of novel approaches to investigating the relaxation dynamics of active glasses via an active elastoplastic model (EPM). These approaches describe the relaxation dynamics via local plastic yielding and can provide improved insights as we can study various aspects of the system separately. Activity enters the model via three crucial features: activity-mediated plastic yielding, activated barrier crossing, and persistent rotational dynamics of the yielding direction. We first consider a minimal active EPM that adds the effect of active yielding to a thermal EPM. We show that this active EPM captures the known results of active glasses within a reasonable parameter space. The results also agree well with the analytical results for active glasses when activity is small. The minimal model breaks down at very low temperatures where other effects become important. Looking at the broader model class, we demonstrate that whereas active yielding primarily dominates the relaxation dynamics, the persistence of the yielding direction governs the dynamic heterogeneity in active glasses.

Emergent mesoscale correlations in active solids with noisy chiral dynamics

We present the linear response theory for an elastic solid composed of active Brownian particles with intrinsic individual chirality, deriving both a normal mode formulation and a continuum elastic formulation. Using this theory, we compute analytically the velocity correlations and energy spectra under different conditions, showing an excellent agreement with simulations. We generate the corresponding phase diagram, identifying chiral and achiral disordered regimes (for high chirality or noise levels), as well as chiral and achiral states with mesoscopic-range order (for low chirality and noise). The chiral ordered states display mesoscopic spatial correlations and oscillating time correlations, but no wave propagation. In the high chirality regime, we find a peak in the elastic energy spectrum that leads to a non-monotonic behavior with increasing noise strength that is consistent with the emergence of the 'hammering state' recently identified in chiral glasses. Finally, we show numerically that our theory, despite its linear response nature, can be applied beyond the idealized homogeneous solid assumed in our derivations. Indeed, by increasing the level of activity, we show that it remains a good approximation of the system dynamics until just below the melting transition. In addition, we show that there is still an excellent agreement between our analytical results and simulations when we extend our results to heterogeneous solids composed of mixtures of active particles with different intrinsic chirality and noise levels. The derived linear response theory is therefore robust and applicable to a broad range of real-world active systems. Our work provides a thorough analytical and numerical description of the emergent states in a densely packed system of chiral self-propelled Brownian disks, thus allowing a detailed understanding of the phases and dynamics identified in a minimal chiral active system.

Motility driven glassy dynamics in confluent epithelial monolayers

As wounds heal, embryos develop, cancer spreads, or asthma progresses, the cellular monolayer undergoes a glass transition between solid-like jammed and fluid-like flowing states. During some of these processes, the cells undergo an epithelial-to-mesenchymal transition (EMT): they acquire in-plane polarity and become motile. Thus, how motility drives the glassy dynamics in epithelial systems is critical for the EMT process. However, no analytical framework that is indispensable for deeper insights exists. Here, we develop such a theory inspired by a well-known glass theory. One crucial result of this work is that the confluency affects the effective persistence time-scale of active force, described by its rotational diffusivity, Deffr. Deffr differs from the bare rotational diffusivity, Dr, of the motile force due to cell shape dynamics, which acts to rectify the force dynamics: Deffr is equal to Dr when Dr is small and saturates when Dr is large. We test the theoretical prediction of Deffr and how it affects the relaxation dynamics in our simulations of the active Vertex model. This novel effect of Deffr is crucial to understanding the new and previously published simulation data of active glassy dynamics in epithelial monolayers.

Extremely persistent dense active fluids

We study the dynamics of dense three-dimensional systems of active particles for large persistence times τp at constant average self-propulsion force f. These systems are fluid counterparts of previously investigated extremely persistent systems, which in the large persistence time limit relax only on the time scale of τp. We find that many dynamic properties of the systems we study, such as the mean-squared velocity, the self-intermediate scattering function, and the shear-stress correlation function, become τp-independent in the large persistence time limit. In addition, the large τp limits of many dynamic properties, such as the mean-square velocity and the relaxation times of the scattering function, and the shear-stress correlation function, depend on f as power laws with non-trivial exponents. We conjecture that these systems constitute a new class of extremely persistent active systems.

Rotational inertia-induced glassy transition in chiral particle systems

The dense active matter exhibits characteristics reminiscent of traditional glassy phenomena, yet the role of rotational inertia in glass dynamics remains elusive. In this study, we investigate the glass dynamics of chiral active particles influenced by rotational inertia. Rotational inertia endows exponential memory to particle orientation, restricting its alteration and amplifying the effective persistence time. At lower spinning frequencies, the diffusion coefficient exhibits a peak function relative to rotational inertia for shorter persistence times, while it steadily increases with rotational inertia for longer persistence times. In the realm of high-frequency spinning, the impact of rotational inertia on diffusion behavior becomes more pronounced, resulting in a nonmonotonic and intricate relationship between the diffusion coefficient and rotational inertia. Consequently, the introduction of rotational inertia significantly alters the glassy dynamics of chiral active particles, allowing for the control over transitions between fluid and glassy states by modulating rotational inertia. Moreover, our findings indicate that at a specific spinning temperature, there exists an optimal spinning frequency at which the diffusion coefficient attains its maximum value.

Glassy Dynamics in Chiral Fluids

Chiral active matter is enjoying a rapid increase of interest, spurred by the rich variety of asymmetries that can be attained in, e.g., the shape or self-propulsion mechanism of active particles. Though this has already led to the observance of so-called chiral crystals, active chiral glasses remain largely unexplored. A possible reason for this could be the naive expectation that interactions dominate the glassy dynamics and the details of the active motion become increasingly less relevant. Here, we show that quite the opposite is true by studying the glassy dynamics of interacting chiral active Brownian particles. We demonstrate that when our chiral fluid is pushed to glassy conditions, it exhibits highly nontrivial dynamics, especially compared to a standard linear active fluid such as common active Brownian particles. Despite the added complexity, we are still able to present a full rationalization for all identified dynamical regimes. Most notably, we introduce a new "hammering" mechanism, unique to rapidly spinning particles in high-density conditions, that can fluidize a chiral active solid.

Mode-coupling theory for the dynamics of dense underdamped active Brownian particle system

We present a theory to study the inertial effect on glassy dynamics of the underdamped active Brownian particle (UABP) system. Using the assumption of the nonequilibrium steady-state, we obtain an effective Fokker-Planck equation for the probability distribution function (PDF) as a function of positions and momentums. With this equation, we achieve the evolution equation of the intermediate scattering function through the Zwanzig-Mori projection operator method and the mode-coupling theory (MCT). Theoretical analysis shows that the inertia of the particle affects the memory function and corresponding glass transition by influencing the structure factor and a velocity correlation function. The theory provides theoretical support and guidance for subsequent simulation work.

Active glassy dynamics is unaffected by the microscopic details of self-propulsion

Recent years have seen a rapid increase of interest in dense active materials, which, in the disordered state, share striking similarities with the conventional passive glass-forming matter. For such passive glassy materials, it is well established (at least in three dimensions) that the details of the microscopic dynamics, e.g., Newtonian or Brownian, do not influence the long-time glassy behavior. Here, we investigate whether this still holds true in the non-equilibrium active case by considering two simple and widely used active particle models, i.e., active Ornstein-Uhlenbeck particles (AOUPs) and active Brownian particles (ABPs). In particular, we seek to gain more insight into the role of the self-propulsion mechanism on the glassy dynamics by deriving a mode-coupling theory (MCT) for thermal AOUPs, which can be directly compared to a recently developed MCT for ABPs. Both theories explicitly take into account the active degrees of freedom. We solve the AOUP- and ABP-MCT equations in two dimensions and demonstrate that both models give almost identical results for the intermediate scattering function over a large variety of control parameters (packing fractions, active speeds, and persistence times). We also confirm this theoretical equivalence between the different self-propulsion mechanisms numerically via simulations of a polydisperse mixture of active quasi-hard spheres, thereby establishing that, at least for these model systems, the microscopic details of self-propulsion do not alter the active glassy behavior.

Motile Topological Defects Hinder Dynamical Arrest in Dense Liquids of Active Ellipsoids

Recent numerical studies have identified the persistence time of active motion as a critical parameter governing glassy dynamics in dense active matter. Here we studied dynamics in liquids of granular active ellipsoids with tunable persistence and velocity. We show that increasing the persistence time at moderate supercooling is equivalent to increasing the strength of attraction in equilibrium liquids and results in reentrant dynamics not just in the translational degrees of freedom, as anticipated, but also in the orientational ones. However, at high densities, motile topological defects, unique to active liquids of elongated particles, hindered dynamical arrest. Most remarkably, for the highest activity, we observed intermittent dynamics due to the jamming-unjamming of these defects for the first time.

Cage Length Controls the Nonmonotonic Dynamics of Active Glassy Matter

Dense active matter is gaining widespread interest due to its remarkable similarity with conventional glass-forming materials. However, active matter is inherently out of equilibrium and even simple models such as active Brownian particles (ABPs) and active Ornstein-Uhlenbeck particles (AOUPs) behave markedly differently from their passive counterparts. Controversially, this difference has been shown to manifest itself via either a speedup, slowdown, or nonmonotonic change of the glassy relaxation dynamics. Here we rationalize these seemingly contrasting views on the departure from equilibrium by identifying the ratio of the short-time length scale to the cage length, i.e., the length scale of local particle caging, as a vital and unifying control parameter for active glassy matter. In particular, we explore the glassy dynamics of both thermal and athermal ABPs and AOUPs upon increasing the persistence time. We find that for all studied systems there is an optimum of the dynamics; this optimum occurs when the cage length coincides with the corresponding short-time length scale of the system, which is either the persistence length for athermal systems or a combination of the persistence length and a diffusive length scale for thermal systems. This new insight, for which we also provide a simple physical argument, allows us to reconcile and explain the manifestly disparate departures from equilibrium reported in many previous studies of dense active materials.

Generalized mode-coupling theory for mixtures of Brownian particles

Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the nontrivial similarity between Brownian and Newtonian MCT is maintained for our multicomponent GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.

Mode-coupling theory for tagged-particle motion of active Brownian particles

We derive a mode-coupling theory (MCT) to describe the dynamics of a tracer particle that is embedded in a dense system of active Brownian particles (ABPs) in two spatial dimensions. The ABP undergo translational and rotational Brownian motion and are equipped with a fixed self-propulsion speed along their orientational vector that describes their active motility. The resulting equations of motion for the tagged-particle density-correlation functions describe the various cases of tracer dynamics close to the glass transition: that of a single active particle in a glass-forming passive host suspensions, that of a passive colloidal particle in a suspension of ABP, and that of active tracers in a bath of active particles. Numerical results are presented for these cases assuming hard-sphere interactions among the particles. The qualitative and quantitative accuracy of the theory is tested against event-driven Brownian dynamics (ED-BD) simulations of active and passive hard disks. Simulation and theory are found in quantitative agreement, provided one adjusts the overall density (as known from the passive description of glassy dynamics), and allows for a rescaling of self-propulsion velocities in the active host system. These adjustments account for the fact that ABP-MCT generally overestimates the tendency for kinetic arrest. We confirm in the simulations a peculiar feature of the transient and stationary dynamical density-correlation functions regarding their lack of symmetry under time reversal, demonstrating the nonequilibrium nature of the system and how it manifests itself in the theory.

Transport coefficients in dense active Brownian particle systems: mode-coupling theory and simulation results

We discuss recent advances in developing a mode-coupling theory of the glass transition (MCT) of two-dimensional systems of active Brownian particles (ABPs). The theory describes the structural relaxation close to the active glass in terms of transient dynamical density correlation functions. We summarize the equations of motion that have been derived for the collective density-fluctuation dynamics and those for the tagged-particle motion. The latter allow to study the dynamics of both passive and active tracers in both passive and active host systems. In the limit of small wave numbers, they give rise to equations of motion describing the mean-squared displacements (MSDs) of these tracers and hence the long-time diffusion coefficients as a transport coefficient quantifying long-range tracer motion. We specifically discuss the case of a single ABP tracer in a glass-forming passive host suspension, a case that has recently been studied in experiments on colloidal Janus particles. We employ event-driven Brownian dynamics (ED-BD) computer simulations to test the ABP-MCT and find good agreement between the two for the MSD, provided that known errors in MCT already for the passive system (i.e., an overestimation of the glassiness of the system) are accounted for by an empirical mapping of packing fractions and host-system self-propulsion forces. The ED-BD simulation results also compare well to experimental data, although a peculiar non-monotonic mapping of self-propulsion velocities is required. The ABP-MCT predicts a specific self-propulsion dependence of the Stokes-Einstein relation between the long-time diffusion coefficient and the host-system viscosity that matches well the results from simulation. An application of ABP-MCT within the integration-through transients framework to calculate the density-renormalized effective swim velocity of the interacting ABP agrees qualitatively with the ED-BD simulation data at densities close to the glass transition and quantitatively for the full density range only after the mapping of packing fractions employed for the passive system.

The "isothermal" compressibility of active matter

We demonstrate that the mechanically defined "isothermal" compressibility behaves as a thermodynamic-like response function for suspensions of active Brownian particles. The compressibility computed from the active pressure-a combination of the collision and unique swim pressures-is capable of predicting the critical point for motility induced phase separation, as expected from the mechanical stability criterion. We relate this mechanical definition to the static structure factor via an active form of the thermodynamic compressibility equation and find the two to be equivalent, as would be the case for equilibrium systems. This equivalence indicates that compressibility behaves like a thermodynamic response function, even when activity is large. Finally, we discuss the importance of the phase interface when defining an active chemical potential. Previous definitions of the active chemical potential are shown to be accurate above the critical point but breakdown in the coexistence region. Inclusion of the swim pressure in the mechanical compressibility definition suggests that the interface is essential for determining phase behavior.

Escape of a passive particle from an activity-induced energy landscape: emergence of slow and fast effective diffusion

Spontaneous persistent motions driven by active processes play a central role in maintaining living cells far from equilibrium. In the majority of research studies, the steady state dynamics of an active system has been described in terms of an effective temperature. By contrast, we have examined a prototype model for diffusion in an activity-induced rugged energy landscape to describe the slow dynamics of a tagged particle in a dense active environment. The expression for the mean escape time from the activity-induced rugged energy landscape holds only in the limit of low activity and the mean escape time from the rugged energy landscape increases with activity. The precise form of the active correlation will determine whether the mean escape time will depend on the persistence time or not. The activity-induced rugged energy landscape approach also allows an estimate of the non-equilibrium effective diffusivity characterizing the slow diffusive motion of the tagged particle due to activity. On the other hand, in a dilute environment, high activity augments the diffusion of the tagged particle. The enhanced diffusion can be attributed to an effective temperature higher than the ambient temperature and this is used to calculate the Kramers' mean escape time, which decreases with activity. Our results have direct relevance to recent experiments on tagged particle diffusion in condensed phases.

Non-monotonic dependence of polymer chain dynamics on active crowder size

Configuration dynamics of flexible polymer chains is of ubiquitous importance in many biological processes. Here, we investigate a polymer chain immersed in a bath of size-changed active particles in two dimensional space using Langevin dynamics simulations. Particular attention is paid to how the radius of gyration R_g of the polymer chain depends on the size σ_c of active crowders. We find that R_g shows nontrivial non-monotonic dependence on σ_c: The chain first swells upon increasing σ_c, reaching a fully expanded state with maximum R_g, and then, R_g decreases until the chain collapses to a compact coil state if the crowder is large enough. Interestingly, the chain may oscillate between a collapse state and a stretched state at moderate crowder size. Analysis shows that it is the competition between two effects of active particles, one stretching the chain from inside due to persistence motion and the other compressing the chain from outside, that leads to the non-monotonic dependence. Besides, the diffusion of the polymer chain also shows nontrivial non-monotonic dependence on σ_c. Our results demonstrate the important interplay between particle activity and size associated with polymer configurations in active crowding environments.

Complex modulus and compliance for airway smooth muscle cells

A cell can be described as a complex viscoelastic material with structural relaxations that is modulated by thermal and chemically nonequilibrium processes. Tissue morphology and function rely upon cells' physical responses to mechanical force. We measured the frequency-dependent mechanical relaxation response of adherent human airway smooth muscle cells under adenosine triphosphate (ATP) depletion and normal ATP conditions. The frequency dependence of the complex compliance J^{*} and modulus G^{*} was measured over the frequencies 10^{-1}<f<10^{3} Hz at selected temperatures between 4<T<54^{∘}C. Our results show characteristic relaxation features which can be interpreted by the mode-coupling theory (MCT) of viscoelastic liquids. We analyze the shape of the spectra in terms of a so-called A_{4} scenario with logarithmic scaling laws. Characteristic timescales τ_{β} and τ_{α} appear with corresponding energy barriers E_{β}≈(10-20)k_{B}T and E_{α}≈(20-30)k_{B}T. We demonstrate that cells are close to a glass transition. We find that the cell becomes softer around physiological temperatures, where its surface structure is more liquid-like with a plateau modulus around 0.1-0.8 kPa compared with the more solid-like interior cytoskeletal structures with a plateau modulus 1-15 kPa. Corresponding values for the viscosity are 10^{2}-10^{3} Pa s for the surface structures closer to the membrane and 10^{4}-10^{6} Pa s for the core cytoskeletal structures.

Effective one-component model of binary mixture: molecular arrest induced by the spatially correlated stochastic dynamics

Spatially correlated noise (SCN), i.e. the thermal noise that affects neighbouring particles in a similar manner, is ubiquitous in soft matter systems. In this work, we apply the over-damped SCN-driven Langevin equations as an effective, one-component model of the dynamics in dense binary mixtures. We derive the thermodynamically consistent fluctuation-dissipation relation for SCN to show that it predicts the molecular arrest resembling the glass transition, i.e. the critical slow-down of dynamics in the disordered phases. We show that the mechanism of singular dissipation is embedded in the dissipation matrix, accompanying SCN. We are also able to identify the characteristic length of collective dissipation, which diverges at critical packing. This novel physical quantity conveniently describes the difference between the ergodic and non-ergodic dynamics. The model is fully analytically solvable, one-dimensional and admits arbitrary interactions between the particles. It qualitatively reproduces several different modes of arrested disorder encountered in binary mixtures, including e.g. the re-entrant arrest. The model can be effectively compared to the mode coupling theory.

Active glasses

Active glassy matter has recently emerged as a novel class of non-equilibrium soft matter, combining energy-driven, active particle movement with dense and disordered glass-like behavior. Here we review the state-of-the-art in this field from an experimental, numerical, and theoretical perspective. We consider both non-living and living active glassy systems, and discuss how several hallmarks of glassy dynamics (dynamical slowdown, fragility, dynamical heterogeneity, violation of the Stokes-Einstein relation, and aging) are manifested in such materials. We start by reviewing the recent experimental evidence in this area of research, followed by an overview of the main numerical simulation studies and physical theories of active glassy matter. We conclude by outlining several open questions and possible directions for future work.

Disordered hyperuniform obstacles enhance sorting of dynamically chiral microswimmers

Disordered hyperuniformity, a brand new type of arrangement with novel physical properties, provides various practical applications in extensive fields. To highlight the great potential of applying disordered hyperuniformity to active systems, a practical example is reported here by an optimal sorting of dynamically chiral microswimmers in disordered hyperuniform obstacle environments in comparison with regular or disordered ones. This optimal chirality sorting stems from a competition between advantageous microswimmer-obstacle collisions and disadvantageous trapping of microswimmers by obstacles. Based on this mechanism, optimal chirality sorting is also realized by tuning other parameters including the number density of obstacles, the strength of driven force and the noise intensity. Our findings may open a new perspective on both theoretical and experimental investigations for further applications of disordered hyperuniformity in active systems.

Mode-coupling theory for the steady-state dynamics of active Brownian particles

We present a theory for the steady-state dynamics of a two-dimensional system of spherically symmetric active Brownian particles. The derivation of the theory consists of two steps. First, we integrate out the self-propulsions and obtain a many-particle evolution equation for the probability distribution of the particles' positions. Second, we use the projection operator technique and a mode-coupling-like factorization approximation to derive an equation of motion for the density correlation function. The nonequilibrium character of the active system manifests itself through the presence of a steady-state correlation function that quantifies spatial correlations of microscopic steady-state currents of the particles. This function determines the dependence of the short-time dynamics on the activity. It also enters into the expression for the memory matrix and thus influences the long-time glassy dynamics.

Mode-coupling theory for active Brownian particles

We present a mode-coupling theory (MCT) for the slow dynamics of two-dimensional spherical active Brownian particles (ABPs). The ABPs are characterized by a self-propulsion velocity v_{0} and by their translational and rotational diffusion coefficients D_{t} and D_{r}, respectively. Based on the integration-through-transients formalism, the theory requires as input only the equilibrium static structure factors of the passive system (where v_{0}=0). It predicts a nontrivial idealized-glass-transition diagram in the three-dimensional parameter space of density, self-propulsion velocity, and rotational diffusivity that arise because at high densities, the persistence length of active swimming ℓ_{p}=v_{0}/D_{r} interferes with the interaction length ℓ_{c} set by the caging of particles. While the low-density dynamics of ABPs is characterized by a single Péclet number Pe=v_{0}^{2}/D_{r}D_{t}, close to the glass transition the dynamics is found to depend on Pe and ℓ_{p} separately. At fixed density, increasing the self-propulsion velocity causes structural relaxation to speed up, while decreasing the persistence length slows down the relaxation. The active-MCT glass is a nonergodic state that is qualitatively different from the passive glass. In it, correlations of initial density fluctuations never fully decay, but also an infinite memory of initial orientational fluctuations is retained in the positions.

Nonequilibrium mode-coupling theory for dense active systems of self-propelled particles

The physics of active systems of self-propelled particles, in the regime of a dense liquid state, is an open puzzle of great current interest, both for statistical physics and because such systems appear in many biological contexts. We develop a nonequilibrium mode-coupling theory (MCT) for such systems, where activity is included as a colored noise with the particles having a self-propulsion force f₀ and a persistence time τ_p. Using the extended MCT and a generalized fluctuation-dissipation theorem, we calculate the effective temperature T_eff of the active fluid. The nonequilibrium nature of the systems is manifested through a time-dependent T_eff that approaches a constant in the long-time limit, which depends on the activity parameters f₀ and τ_p. We find, phenomenologically, that this long-time limit is captured by the potential energy of a single, trapped active particle (STAP). Through a scaling analysis close to the MCT glass transition point, we show that τ_α, the α-relaxation time, behaves as τ_α ∼ f₀^-2γ, where γ = 1.74 is the MCT exponent for the passive system. τ_α may increase or decrease as a function of τ_p depending on the type of active force correlations, but the behavior is always governed by the same value of the exponent γ. Comparison with the numerical solution of the nonequilibrium MCT and simulation results give excellent agreement with scaling analysis.