Active glasses

Length-scale dependence of Stokes-Einstein breakdown in active glass-forming liquids

The Stokes-Einstein (SE) relation, which relates diffusion constants with the viscosity of a liquid at high temperatures in equilibrium, is violated in the supercooled temperature regime. Whether this relation is obeyed in nonequilibrium active liquids is a question of significant current interest to the statistical physics community trying to develop the theoretical framework of nonequilibrium statistical mechanics. Via extensive computer simulations of model active glass-forming liquids in three dimensions, we show that SE is obeyed at a high temperature similar to the equilibrium behavior, and it gets violated in the supercooled temperature regimes. The degree of violation increases systematically with the increasing activity which quantifies the amount the system is driven out of equilibrium. First passage-time (FPT) distributions helped us to gain insights into this enhanced breakdown from the increased short-time peak, depicting hoppers. Subsequently, we study the wave vector dependence of SE relation and show that it gets restored at a wave vector that decreases with increasing activity, and the crossover wave vector is found to be proportional to the inverse of the dynamical heterogeneity length scale in the system. Our work showed how SE violation in active supercooled liquids could be rationalized using the growth of dynamic length scale, which is found to grow enormously with increasing activity in these systems.

Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Emergence of bacterial glass

Densely packed, motile bacteria can adopt collective states not seen in conventional, passive materials. These states remain in many ways mysterious, and their physical characterization can aid our understanding of natural bacterial colonies and biofilms as well as materials in general. Here, we overcome challenges associated with generating uniformly growing, large, quasi-two-dimensional bacterial assemblies by a membrane-based microfluidic device and report the emergence of glassy states in two-dimensional suspension of Escherichia coli. As the number density increases by cell growth, populations of motile bacteria transition to a glassy state, where cells are packed and unable to move. This takes place in two steps, the first one suppressing only the orientational modes and the second one vitrifying the motion completely. Characterizing each phase through statistical analyses and investigations of individual motion of bacteria, we find not only characteristic features of glass such as rapid slowdown, dynamic heterogeneity, and cage effects, but also a few properties distinguished from those of thermal glass. These distinctive properties include the spontaneous formation of micro-domains of aligned cells with collective motion, the appearance of an unusual signal in the dynamic susceptibility, and the dynamic slowdown with a density dependence generally forbidden for thermal systems. Our results are expected to capture general characteristics of such active rod glass, which may serve as a physical mechanism underlying dense bacterial aggregates.

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.

Microscopic theory for nonequilibrium correlation functions in dense active fluids

One of the key hallmarks of dense active matter in the liquid, supercooled, and solid phases is the so-called equal-time velocity correlations. Crucially, these correlations can emerge spontaneously, i.e., they require no explicit alignment interactions, and therefore represent a generic feature of dense active matter. This indicates that for a meaningful comparison or possible mapping between active and passive liquids one not only needs to understand their structural properties, but also the impact of these velocity correlations. This has already prompted several simulation and theoretical studies, though they are mostly focused on athermal systems and thus overlook the effect of translational diffusion. Here, we present a fully microscopic method to calculate nonequilibrium correlations in two-dimensional systems of thermal active Brownian particles (ABPs). We use the integration through transients formalism together with (active) mode-coupling theory and analytically calculate qualitatively consistent static structure factors and active velocity correlations. We complement our theoretical results with simulations of both thermal and athermal ABPs which exemplify the disruptive role that thermal noise has on velocity correlations.

Activity-dependent glassy cell mechanics II: Nonthermal fluctuations under metabolic activity

The glassy cytoplasm, crowded with bio-macromolecules, is fluidized in living cells by mechanical energy derived from metabolism. Characterizing the living cytoplasm as a nonequilibrium system is crucial in elucidating the intricate mechanism that relates cell mechanics to metabolic activities. In this study, we conducted active and passive microrheology in eukaryotic cells, and quantified nonthermal fluctuations by examining the violation of the fluctuation-dissipation theorem. The power spectral density of active force generation was estimated following the Langevin theory extended to nonequilibrium systems. However, experiments performed while regulating cellular metabolic activity showed that the nonthermal displacement fluctuation, rather than the active nonthermal force, is linked to metabolism. We discuss that mechano-enzymes in living cells do not act as microscopic objects. Instead, they generate meso-scale collective fluctuations with displacements controlled by enzymatic activity. The activity induces structural relaxations in glassy cytoplasm. Even though the autocorrelation of nonthermal fluctuations is lost at long timescales due to the structural relaxations, the nonthermal displacement fluctuation remains regulated by metabolic reactions. Our results therefore demonstrate that nonthermal fluctuations serve as a valuable indicator of a cell's metabolic activities, regardless of the presence or absence of structural relaxations.

Fluid-Glass-Jamming Rheology of Soft Active Brownian Particles

We numerically study the shear rheology of a binary mixture of soft active Brownian particles, from the fluid to the disordered solid regime. At low shear rates, we find a Newtonian regime, where a Green-Kubo relation with an effective temperature provides the linear viscosity. It is followed by a shear-thinning regime at high shear rates. At high densities, solidification is signaled by the emergence of a finite yield stress. We construct a "fluid-glass-jamming" phase diagram with activity replacing temperature. While both parameters gauge fluctuations, activity also changes the exponent characterizing the decay of the diffusivity close to the glass transition and the shape of the yield stress surface. The dense disordered active solid appears to be mostly dominated by athermal jamming rather than glass rheology.

Method to probe the pronounced growth of correlation lengths in active glass-forming liquids using an elongated probe

The growth of correlation lengths in equilibrium glass-forming liquids near the glass transition is considered a critical finding in the quest to understand the physics of glass formation. These understandings helped us understand various dynamical phenomena observed in supercooled liquids. It is known that at least two different length scales exist; one is of thermodynamic origin, while the other is dynamical in nature. Recent observations of glassy dynamics in biological and synthetic systems where the external or internal driving source controls the dynamics, apart from the usual thermal noise, lead to the emergence of the field of active glassy matter. A question of whether the physics of glass formation in these active systems is also accompanied by growing dynamic and static lengths is indeed timely. In this article, we probe the growth of dynamic and static lengths in a model active glass system using rod-like elongated probe particles, an experimentally viable method. We show that the dynamic and static lengths in these nonequilibrium systems grow much more rapidly than their passive counterparts. We then offer an understanding of the violation of the Stokes-Einstein relation and Stokes-Einstein-Debye relation using these lengths via a scaling theory.

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.

Reconfiguration, Interrupted Aging, and Enhanced Dynamics of a Colloidal Gel Using Photoswitchable Active Doping

We study quasi-2D gels made of a colloidal network doped with Janus particles activated by light. Following the gel formation, we monitor both the structure and dynamics before, during, and after the activation period. Before activity is switched on, the gel is slowly aging. During the activation, the mobility of the passive particles exhibits a characteristic scale-dependent response, while the colloidal network remains connected, and the gel maintains its structural integrity. Once activity is switched off, the gel stops aging and keeps the memory of the structure inherited from the active phase. Remarkably, the motility remains larger than that of the gel, before the active period. The system has turned into a genuinely softer gel, with frozen dynamics, but with more space for thermal fluctuations. The above conclusions remain valid long after the activity period.

Activity-dependent glassy cell mechanics Ⅰ: Mechanical properties measured with active microrheology

Active microrheology was conducted in living cells by applying an optical-trapping force to vigorously fluctuating tracer beads with feedback-tracking technology. The complex shear modulus G(ω)=G'(ω)-iG″(ω) was measured in HeLa cells in an epithelial-like confluent monolayer. We found that G(ω)∝(-iω)1/2 over a wide range of frequencies (1 Hz < ω/2π < 10 kHz). Actin disruption and cell-cycle progression from G1 to S and G2 phases only had a limited effect on G(ω) in living cells. On the other hand, G(ω) was found to be dependent on cell metabolism; ATP-depleted cells showed an increased elastic modulus G'(ω) at low frequencies, giving rise to a constant plateau such that G(ω)=G0+A(-iω)1/2. Both the plateau and the additional frequency dependency ∝(-iω)1/2 of ATP-depleted cells are consistent with a rheological response typical of colloidal jamming. On the other hand, the plateau G0 disappeared in ordinary metabolically active cells, implying that living cells fluidize their internal states such that they approach the critical jamming point.

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.

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

Disordered Collective Motion in Dense Assemblies of Persistent Particles

We explore the emergence of nonequilibrium collective motion in disordered nonthermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field.

Chemotactic self-caging in active emulsions

The out-of-equilibrium dynamics of chemotactic active matter—be it animate or inanimate—is closely coupled to the environment, a chemical landscape shaped by secretions from the motile agents, fuel uptake, or autochemotactic signaling. This gives rise to complex collective effects, which can be exploited by the agents for colony migration strategies or pattern formation. We study such effects using an idealized experimental system: self-propelled microdroplets that communicate via chemorepulsive trails. We present a comprehensive experimental analysis that involves direct probing of the diffusing chemical trails and the trail–droplet interactions and use it to construct a generic theoretical model. We connect these repulsive autochemotactic interactions to the collective dynamics in emulsions, demonstrating a state of dynamical arrest: chemotactic self-caging.

A common feature of biological self-organization is how active agents communicate with each other or their environment via chemical signaling. Such communications, mediated by self-generated chemical gradients, have consequences for both individual motility strategies and collective migration patterns. Here, in a purely physicochemical system, we use self-propelling droplets as a model for chemically active particles that modify their environment by leaving chemical footprints, which act as chemorepulsive signals to other droplets. We analyze this communication mechanism quantitatively both on the scale of individual agent–trail collisions as well as on the collective scale where droplets actively remodel their environment while adapting their dynamics to that evolving chemical landscape. We show in experiment and simulation how these interactions cause a transient dynamical arrest in active emulsions where swimmers are caged between each other’s trails of secreted chemicals. Our findings provide insight into the collective dynamics of chemically active particles and yield principles for predicting how negative autochemotaxis shapes their navigation strategy.

Aging or DEAD: Origin of the non-monotonic response to weak self-propulsion in active glasses

Among amorphous states, glass is defined by relaxation times longer than the observation time. This nonergodic nature makes the understanding of glassy systems an involved topic, with complex aging effects or responses to further out-of-equilibrium external drivings. In this respect, active glasses made of self-propelled particles have recently emerged as a stimulating systems, which broadens and challenges our current understanding of glasses by considering novel internal out-of-equilibrium degrees of freedom. In previous experimental studies we have shown that in the ergodicity broken phase, the dynamics of dense passive particles first slows down as particles are made slightly active, before speeding up at larger activity. Here, we show that this nonmonotonic behavior also emerges in simulations of soft active Brownian particles and explore its cause. We refute that the deadlock by emergence of active directionality model we proposed earlier describes our data. However, we demonstrate that the nonmonotonic response is due to activity enhanced aging and thus confirm the link with ergodicity breaking. Beyond self-propelled systems, our results suggest that aging in active glasses is not fully understood.

Multiscale heterogeneous dynamics in two-dimensional glassy colloids

On approaching the glass transition, a dense colloid exhibits a dramatic slowdown with minute structural changes. Most microscopy experiments directly follow the motion of individual particles in real space, whereas scattering experiments typically probe the collective dynamics in reciprocal space, at variable wavevector q. Multiscale studies of glassy dynamics are experimentally demanding and thus seldom performed. By using two-dimensional hard-sphere colloids at various area fractions φ, we show here that Differential Dynamic Microscopy (DDM) can be effectively used to measure the collective dynamics of a glassy colloid in a range of q within a single experiment. As φ is increased, the single decay of the intermediate scattering functions is progressively replaced by a more complex relaxation that we fit to a sum of two stretched-exponential decays. The slowest process, corresponding to the long-time particle escapes from caging, has a characteristic time τs = 1/(DLq2 ) with diffusion coefficient DL ∼ (φc −φ)2.8 , and φc ≈ 0.81. The fast process exhibits, instead, a non-Brownian scaling of the characteristic time τf(q) and a relative amplitude a(q) that monotonically increases with q. Despite the non-Brownian nature of τf(q), we succeed in estimating the short-time diffusion coefficient Dcage, whose φ-dependence is practically negligible compared to the one of DL. Finally, we extend DDM to measure the q-dependent dynamical susceptibility χ4(q,t), a powerful yet hard-to-access multiscale indicator of dynamical heterogeneities. Our results show that DDM is a convenient tool to study the dynamics of colloidal glasses over a broad range of time and length-scales.

Controlling the volume fraction of glass-forming colloidal suspensions using thermosensitive host "mesogels"

The key parameter controlling the glass transition of colloidal suspensions is φ, the fraction of the sample volume occupied by the particles. Unfortunately, changing φ by varying an external parameter, e.g., temperature T as in molecular glass formers, is not possible, unless one uses thermosensitive colloidal particles, such as the popular poly(N-isopropylacrylamide) (PNiPAM) microgels. These, however, have several drawbacks, including high deformability, osmotic deswelling, and interpenetration, which complicate their use as a model system to study the colloidal glass transition. Here, we propose a new system consisting of a colloidal suspension of non-deformable spherical silica nanoparticles, in which PNiPAM hydrogel spheres of ∼100-200μm size are suspended. These non-colloidal "mesogels" allow for controlling the sample volume effectively available to the silica nanoparticles and hence their φ, thanks to the T-induced change in mesogels' volume. Using optical microscopy, we first show that the mesogels retain their ability to change size with T when suspended in Ludox suspensions, similarly as in water. We then show that their size is independent of the sample thermal history such that a well-defined, reversible relationship between T and φ may be established. Finally, we use space-resolved dynamic light scattering to demonstrate that, upon varying T, our system exhibits a broad range of dynamical behaviors across the glass transition and beyond, comparable with those exhibited by a series of distinct silica nanoparticle suspensions of various φ.

A broad perspective to particle-laden fluid interfaces systems: from chemically homogeneous particles to active colloids

Particles adsorbed to fluid interfaces are ubiquitous in industry, nature or life. The wide range of properties arising from the assembly of particles at fluid interface has stimulated an intense research activity on shed light to the most fundamental physico-chemical aspects of these systems. These include the mechanisms driving the equilibration of the interfacial layers, trapping energy, specific inter-particle interactions and the response of the particle-laden interface to mechanical perturbations and flows. The understanding of the physico-chemistry of particle-laden interfaces becomes essential for taking advantage of the particle capacity to stabilize interfaces for the preparation of different dispersed systems (emulsions, foams or colloidosomes) and the fabrication of new reconfigurable interface-dominated devices. This review presents a detailed overview of the physico-chemical aspects that determine the behavior of particles trapped at fluid interfaces. This has been combined with some examples of real and potential applications of these systems in technological and industrial fields. It is expected that this information can provide a general perspective of the topic that can be exploited for researchers and technologist non-specialized in the study of particle-laden interfaces, or for experienced researcher seeking new questions to solve.

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.

Are cell jamming and unjamming essential in tissue development?

The last decade has seen a surge of evidence supporting the existence of the transition of the multicellular tissue from a collective material phase that is regarded as being jammed to a collective material phase that is regarded as being unjammed. The jammed phase is solid-like and effectively 'frozen', and therefore is associated with tissue homeostasis, rigidity, and mechanical stability. The unjammed phase, by contrast, is fluid-like and effectively 'melted', and therefore is associated with mechanical fluidity, plasticity and malleability that are required in dynamic multicellular processes that sculpt organ microstructure. Such multicellular sculpturing, for example, occurs during embryogenesis, growth and remodeling. Although unjamming and jamming events in the multicellular collective are reminiscent of those that occur in the inert granular collective, such as grain in a hopper that can flow or clog, the analogy is instructive but limited, and the implications for cell biology remain unclear. Here we ask, are the cellular jamming transition and its inverse --the unjamming transition-- mere epiphenomena? That is, are they dispensable downstream events that accompany but neither cause nor quench these core multicellular processes? Drawing from selected examples in developmental biology, here we suggest the hypothesis that, to the contrary, the graded departure from a jammed phase enables controlled degrees of malleability as might be required in developmental dynamics. We further suggest that the coordinated approach to a jammed phase progressively slows those dynamics and ultimately enables long-term mechanical stability as might be required in the mature homeostatic multicellular tissue.

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.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

Interparticle Adhesion Regulates the Surface Roughness of Growing Dense Three-Dimensional Active Particle Aggregates

Activity and self-generated motion are fundamental features observed in many living and nonliving systems. Given that interparticle adhesive forces can regulate particle dynamics, we investigate how interparticle adhesion strength controls the boundary growth and roughness of active particle aggregates. Using particle based simulations incorporating both activity (birth, death, and growth) and systematic physical interactions (elasticity and adhesion), we establish that interparticle adhesion strength (f^ad) controls the surface roughness of a densely packed three-dimensional(3D) active particle aggregate expanding into a highly viscous medium. We discover that the surface roughness of a 3D active particle aggregate increases in proportion to the interparticle adhesion strength (f^ad) and show that asymmetry in the radial and transverse active particle mean-squared displacement (MSD) suppresses 3D surface roughness at lower adhesion strengths. By analyzing the statistical properties of particle displacements at the aggregate periphery, we determine that the 3D surface roughness is driven by the movement of active particle toward the core at high interparticle adhesion strengths. Our results elucidate the physics controlling the expansion of adhesive 3D active particle collectives into a highly viscous medium, with implications into understanding stochastic interface growth in active matter systems characterized by self-generation of particles.

Tagged-particle motion of Percus-Yevick hard spheres from first principles

We develop a first-principles-based generalized mode-coupling theory (GMCT) for the tagged-particle motion of glassy systems. This theory establishes a hierarchy of coupled integro-differential equations for self-multi-point density correlation functions, which can formally be extended up to infinite order. We use our GMCT framework to calculate the self-nonergodicity parameters and the self-intermediate scattering function for the Percus-Yevick hard-sphere system based on the first few levels of the GMCT hierarchy. We also test the scaling laws in the α- and β-relaxation regimes near the glass-transition singularity. Furthermore, we study the mean-square displacement and the Stokes-Einstein relation in the supercooled regime. We find that qualitatively our GMCT results share many similarities with the well-established predictions from standard mode-coupling theory, but the quantitative results change, and typically improve, by increasing the GMCT closure level. However, we also demonstrate on general theoretical grounds that the current GMCT framework is unable to account for violation of the Stokes-Einstein relation, underlining the need for further improvements in the first-principles description of glassy dynamics.

How to study a persistent active glassy system

We explore glassy dynamics of dense assemblies of soft particles that are self-propelled by active forces. These forces have a fixed amplitude and a propulsion direction that varies on a timescaleτ_p, the persistence timescale. Numerical simulations of such active glasses are computationally challenging when the dynamics is governed by large persistence times. We describe in detail a recently proposed scheme that allows one to study directly the dynamics in the large persistence time limit, on timescales around and well above the persistence time. We discuss the idea behind the proposed scheme, which we call 'activity-driven dynamics', as well as its numerical implementation. We establish that our prescription faithfully reproduces all dynamical quantities in the appropriate limitτ_p→ ∞. We deploy the approach to explore in detail the statistics of Eshelby-like plastic events in the steady state dynamics of a dense and intermittent active glass.

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.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Multiple dynamic regimes in a coarsening foam

Intermittent dynamics driven by internal stress imbalances in disordered systems is a fascinating yet poorly understood phenomenon. Here, we study it for a coarsening foam. By exploiting differential dynamic microscopy and particle tracking we determine the dynamical characteristics of the foam at different ages in reciprocal and direct space, respectively. At all wavevectors q investigated, the intermediate scattering function exhibits a compressed exponential decay. However, the access to unprecedentedly small values of q highlights the existence of two distinct regimes for the q-dependence of the foam relaxation rate Γ(q). At high q, Γ(q) ∼ q consistent with directionally-persistent and intermittent bubble displacements. At low q, we find the surprising scaling Γ(q) ∼ q ^δ , with δ = 1.6 ± 0.2. The analysis of the bubble displacement distribution in real space reveals the existence of a displacement cut-off of the order of the bubble diameter. Introducing such cut-off length in an existing model, describing stress-driven dynamics in disordered systems, fully accounts for the observed behavior in direct and reciprocal space.

Multiple Types of Aging in Active Glasses

Recent experiments and simulations have revealed glassy features in, e.g., cytoplasm, living tissues and dense assemblies of self-propelled colloids. This leads to a fundamental question: how do these nonequilibrium (active) amorphous materials differ from conventional passive glasses, created by lowering temperature or increasing density? To address this we investigate the aging after a quench to an almost arrested state of a model active glass former, a Kob-Andersen glass in two dimensions. Each constituent particle is driven by a constant propulsion force whose direction diffuses over time. Using extensive molecular dynamics simulations we reveal rich aging behavior of this dense active matter system: short persistence times of the active forcing give effective thermal aging; in the opposite limit we find a two-step aging process with active athermal aging at short times and activity-driven aging at late times. We develop a dedicated simulation method that gives access to this longtime scaling regime for highly persistent active forces.

Universality in incompressible active fluid: Effect of nonlocal shear stress

Phase transitions in active fluids attracted significant attention within the last decades. Recent results show [L. Chen et al., New J. Phys. 17, 042002 (2015)10.1088/1367-2630/17/4/042002] that an order-disorder phase transition in incompressible active fluids belongs to a new universality class. In this work, we further investigate this type of phase transition and focus on the effect of long-range interactions. This is achieved by introducing a nonlocal shear stress into the hydrodynamic description, which leads to superdiffusion of the velocity field, and can be viewed as a result of the active particles performing Lévy walks. The universal properties in the critical region are derived by performing a perturbative renormalization group analysis of the corresponding response functional within the one-loop approximation. We show that the effect of nonlocal shear stress decreases the upper critical dimension of the model, and can lead to the irrelevance of the active fluid self-advection with the resulting model belonging to an unusual long-range Model A universality class not reported before, to our knowledge. Moreover, when the degree of nonlocality is sufficiently high all nonlinearities become irrelevant and the mean-field description is valid in any spatial dimension.

Emergence of Multiscale Dynamics in Colloidal Gels

To gain insight into the kinetics of colloidal gel evolution at low particle volume fractions ϕ, we utilize differential dynamic microscopy to investigate particle aggregation, geometric percolation, and the subsequent transition to nonergodic dynamics. We report the emergence of unexpectedly rich multiscale dynamics upon the onset of nonergodicity, which separates the wave vectors q into three different regimes. In the high-q domain, the gel exhibits ϕ-independent internal vibrations of fractal clusters. The intermediate-q domain is dominated by density fluctuations at the length scale of the clusters, as evidenced by the q independence of the relaxation time τ. In the low-q domain, the scaling of τ as q^{-3} suggests that the network appears homogeneous. The transitions between these three regimes introduce two characteristic length scales, distinct from the cluster size.

Phase separation in a two-dimensional binary colloidal mixture by quorum sensing activity

We present results from Langevin dynamics simulations of a glassy active-passive mixture of soft-repulsive binary colloidal disks. Activity on the smaller particles is applied according to the quorum sensing scheme, in which a smaller particle will be active for a persistence time if its local nearest neighbors are equal to or greater than a certain threshold value. We start with a passive glassy state of the system and apply activity to the smaller particles, which shows a nonmonotonous glassy character of the active particles with the persistence time of the active force, from its passive limit (zero activity). On the other hand, passive particles of the active-passive mixture phase separate at the intermediate persistence time of the active force, resulting in the hexatic-liquid and solid-liquid phases. Thus, our system shows three regimes as active glass, phase separation, and active liquid, as the persistence time increases from its smaller values. We show that the solidlike and hexatic phases consisting of passive large particles are stable due to the smaller momentum transfer from active to passive particles, compared to the higher persistence time where the positional and orientational ordering vanishes. Our model is relevant to active biological systems, where glassy dynamics is present, e.g., bacterial cytoplasm, biological tissues, dense quorum sensing bacteria, and synthetic smart amorphous glasses.