Pushing the glass transition towards random close packing using self-propelled hard spheres

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.

Active Particles Knead Three-Dimensional Gels into Porous Structures

Colloidal gels are prime examples of functional materials exhibiting disordered, amorphous, yet metastable forms. They maintain stability through short-range attractive forces and their material properties are tunable by external forces. Combining persistent homology analyses and simulations of three-dimensional colloidal gels doped with active particles, we reveal novel dynamically evolving structures of colloidal gels. Specifically, we show that the local injection of energy by active dopants can lead to highly porous, yet compact gel structures that can significantly affect the transport of active particles within the modified colloidal gel. We further show how passive interfaces play a topologically significant role in interacting with active particles in 3D. The results open the door to an unexplored prospect of forming a wide variety of compact but highly heterogeneous and percolated porous media through active doping of 3D passive matter, with diverse implications in designing new functional materials for active ground remediation.

Non-equilibrium dynamic hyperuniform states

Disordered hyperuniform structures are an exotic state of matter having suppressed density fluctuations at large length-scale similar to perfect crystals and quasicrystals but without any long range orientational order. In the past decade, an increasing number of non-equilibrium systems were found to have dynamic hyperuniform states, which have emerged as a new research direction coupling both non-equilibrium physics and hyperuniformity. Here we review the recent progress in understanding dynamic hyperuniform states found in various non-equilibrium systems, including the critical hyperuniformity in absorbing phase transitions, non-equilibrium hyperuniform fluids and the hyperuniform structures in phase separating systems via spinodal decomposition.

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.

Structural fluctuations in active glasses

The glassy dynamics of dense active matter have recently become a topic of interest due to their importance in biological processes such as wound healing and tissue development. However, while the liquid-state properties of dense active matter have been studied in relation to the glass transition of active matter, the solid-state properties of active glasses have yet to be understood. In this work, we study the structural fluctuations in the active glasses composed of self-propelled particles. We develop a formalism to describe the solid-state properties of active glasses in the harmonic approximation limit and use it to analyze the displacement fields in the active glasses. Our findings reveal that the dynamics of high-frequency normal modes become quasi-static with respect to the active forces, and consequently, excitations of these modes are significantly suppressed. This leads to a violation of the equipartition law, suppression of particle displacements, and the apparent collective motion of active glasses. Overall, our results provide a fundamental understanding of the solid-state properties of active glasses.

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.

A mean-field theory for predicting single polymer collapse induced by neutral crowders

Macromolecular crowding can induce the collapse of a single long polymer into a globular form due to depletion forces of entropic nature. This phenomenon has been shown to play a significant role in compacting the genome within the bacterium Escherichia coli into a well-defined region of the cell known as the nucleoid. Motivated by the biological significance of this process, numerous theoretical and computational studies have searched for the primary determinants of the behavior of polymer-crowder phases. However, our understanding of this process remains incomplete and there is debate on a quantitatively unified description. In particular, different simulation studies with explicit crowders have proposed different order parameters as potential predictors for the collapse transition. In this work, we present a comprehensive analysis of published simulation data obtained from different sources. Based on the common behavior we find in this data, we develop a unified phenomenological model that we show to be predictive. Finally, to further validate the accuracy of the model, we conduct new simulations on polymers of various sizes, and investigate the role of jamming of the crowders.

Survival-Associated Cellular Response Maintained in Pancreatic Ductal Adenocarcinoma (PDAC) Switched Between Soft and Stiff 3D Microgel Culture

Pancreatic ductal adenocarcinoma (PDAC) accounts for about 90% of all pancreatic cancer cases. Five-year survival rates have remained below 12% since the 1970s, in part due to the difficulty in detection prior to metastasis (migration and invasion into neighboring organs and glands). Mechanical memory is a concept that has emerged over the past decade that may provide a path toward understanding how invading PDAC cells "remember" the mechanical properties of their diseased ("stiff", elastic modulus, E ≈ 10 kPa) microenvironment even while invading a healthy ("soft", E ≈ 1 kPa) microenvironment. Here, we investigated the role of mechanical priming by culturing a dilute suspension of PDAC (FG) cells within a 3D, rheologically tunable microgel platform from hydrogels with tunable mechanical properties. We conducted a suite of acute (short-term) priming studies where we cultured PDAC cells in either a soft (E ≈ 1 kPa) or stiff (E ≈ 10 kPa) environment for 6 h, then removed and placed them into a new soft or stiff 3D environment for another 18 h. Following these steps, we conducted RNA-seq analyses to quantify gene expression. Initial priming in the 3D culture showed persistent gene expression for the duration of the study, regardless of the subsequent environments (stiff or soft). Stiff 3D culture was associated with the downregulation of tumor suppressors (LATS1, BCAR3, CDKN2C), as well as the upregulation of cancer-associated genes (RAC3). Immunofluorescence staining (BCAR3, RAC3) further supported the persistence of this cellular response, with BCAR3 upregulated in soft culture and RAC3 upregulated in stiff-primed culture. Stiff-primed genes were stratified against patient data found in The Cancer Genome Atlas (TCGA). Upregulated genes in stiff-primed 3D culture were associated with decreased survival in patient data, suggesting a link between patient survival and mechanical priming.

Dynamics of self-propelled particles in vibrated dense granular media

We study a system consisting of a few self-propelled particles (SPPs) placed among a crowd of densely packed granular particles that are vertically vibrated in a two-dimensional circular confinement. Our experiments reveal two important findings. First, an SPP exhibits a fractal renewal process within the dense granular medium, which induces a superdiffusive behavior whose diffusion exponent increases with its aspect ratio. Second, the SPPs eventually reach the boundary and form a moving cluster, which transitions from the moving state to the static state as the number of SPPs is increased. These results suggest a simple and effective method of modulating the fluidity and directionality of granular systems via controlling the shape and the number of SPPs.

Stochastic motion under nonlinear friction representing shear thinning

We study stochastic motion under a nonlinear frictional force that levels off with increasing velocity. Specifically, our frictional force is of the so-called Coulomb-tanh type. At small speed, it increases approximately linearly with velocity, while at large speed, it approaches a constant magnitude, similarly to solid (dry, Coulomb) friction. In one spatial dimension, a formal analogy between the associated Fokker-Planck equation and the Schrödinger equation for a quantum mechanical oscillator in a nonharmonic Pöschl-Teller potential is revealed. Then, the stationary velocity statistics can be treated analytically. From such analytical considerations, we determine associated diffusion coefficients, which we confirm by agent-based simulations. Moreover, from such simulations and from numerically solving the associated Fokker-Planck equation, we find that the spatial distribution function, starting from an initial Gaussian shape, develops tails that appear exponential at intermediate timescales. At small magnitudes of stochastic driving, the velocity distribution resembles the case of linear friction, while at large magnitudes, it rather approaches the case of solid (dry, Coulomb) friction. The same is true for diffusion coefficients. In a certain sense thus interpolating between the two extreme cases of linear friction and solid (dry, Coulomb) friction, our approach should be useful to describe several cases of practical relevance. For instance, a reduced increase in friction with increasing relative speed is typical of shear-thinning behavior. Therefore, driven motion in shear-thinning environments is one specific example to which our description may be applied.

Different glassy characteristics are related to either caging or dynamical heterogeneity

Despite the enormous theoretical and application interests, a fundamental understanding of the glassy dynamics remains elusive. The static properties of glassy and ordinary liquids are similar, but their dynamics are dramatically different. What leads to this difference is the central puzzle of the field. Even the primary defining glassy characteristics, their implications, and if they are related to a single mechanism remain unclear. This lack of clarity is a severe hindrance to theoretical progress. Here, we combine analytical arguments and simulations of various systems in different dimensions and address these questions. Our results suggest that the myriad of glassy features are manifestations of two distinct mechanisms. Particle caging controls the mean, and coexisting slow- and fast-moving regions govern the distribution of particle displacements. All the other glassy characteristics are manifestations of these two mechanisms; thus, the Fickian yet non-Gaussian nature of glassy liquids is not surprising. We discover a crossover, from stretched exponential to a power law, in the behavior of the overlap function. This crossover is prominent in simulation data and forms the basis of our analyses. Our results have crucial implications on how the glassy dynamics data are analyzed, challenge some recent suggestions on the mechanisms governing glassy dynamics, and impose strict constraints that a correct theory of glasses must have.

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.

Motility-induced shear thickening in dense colloidal suspensions

Phase transitions and collective dynamics of active colloidal suspensions are fascinating topics in soft matter physics, particularly for out-of-equilibrium systems, which can lead to rich rheological behaviours in the presence of steady shear flow. Here the role of self-propulsion in the rheological response of a dense colloidal suspension is investigated by using particle-resolved Brownian dynamics simulations. First, the combined effect of activity and shear in the solid on the disordering transition of the suspension is analyzed. While both self-propulsion and shear destroy order and melt the system if critical values are exceeded, self-propulsion largely lowers the stress barrier needed to be overcome during the transition. We further explore the rheological response of the active sheared system once a steady state is reached. While passive suspensions show a solid-like behaviour, turning on particle motility fluidises the system. At low self-propulsion, the active suspension behaves in the steady state as a shear-thinning fluid. Increasing the self-propulsion changes the behaviour of the liquid from shear-thinning to shear-thickening. We attribute this to clustering in the sheared suspensions induced by motility. This new phenomenon of motility-induced shear thickening (MIST) can be used to tailor the rheological response of colloidal suspensions.

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.

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

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

Enhanced short time peak in four-point dynamic susceptibility in dense active glass-forming liquids

Active glassy systems are simple model systems that imitate complex biological processes. Sometimes, it becomes crucial to estimate the amount of activity present in such biological systems, such as predicting the progression rate of the cancer cells or the healing time of the wound, etc. In this work, we study a model active glassy system to quantify the degree of activity from the collective, long-wavelength fluctuations in the system. These long-wavelength fluctuations present themselves as an additional peak in the four-point dynamic susceptibility (χ₄(t)) apart from the usual peak at structural relaxation time. We then show how the degree of the activity at such a small timescale can be obtained by measuring the variation in χ₄(t) due to changing activity. A Detailed finite size analysis of the peak height of χ₄(t) suggests the existence of an intrinsic dynamic length scale that grows with increasing activity. Finally, we show that this peak height is a unique function of effective activity across all system sizes, serving as a possible parameter for characterizing the degree of activity in a system.

Unjamming and emergent nonreciprocity in active ploughing through a compressible viscoelastic fluid

A dilute suspension of active Brownian particles in a dense compressible viscoelastic fluid, forms a natural setting to study the emergence of nonreciprocity during a dynamical phase transition. At these densities, the transport of active particles is strongly influenced by the passive medium and shows a dynamical jamming transition as a function of activity and medium density. In the process, the compressible medium is actively churned up - for low activity, the active particle gets self-trapped in a cavity of its own making, while for large activity, the active particle ploughs through the medium, either accompanied by a moving anisotropic wake, or leaving a porous trail. A hydrodynamic approach makes it evident that the active particle generates a long-range density wake which breaks fore-aft symmetry, consistent with the simulations. Accounting for the back-reaction of the compressible medium leads to (i) dynamical jamming of the active particle, and (ii) a dynamical non-reciprocal attraction between two active particles moving along the same direction, with the trailing particle catching up with the leading one in finite time. We emphasize that these nonreciprocal effects appear only when the active particles are moving and so manifest in the vicinity of the jamming-unjamming transition.

Statistics for an object actively driven by spontaneous symmetry breaking into reversible directions

Propulsion of otherwise passive objects is achieved by mechanisms of active driving. We concentrate on cases in which the direction of active drive is subject to spontaneous symmetry breaking. In our case, this direction will be maintained until a large enough impulse by an additional stochastic force reverses it. Examples may be provided by self-propelled droplets, gliding bacteria stochastically reversing their propulsion direction, or nonpolar vibrated hoppers. The magnitude of active forcing is regarded as constant, and we include the effect of inertial contributions. Interestingly, this situation can formally be mapped to stochastic motion under (dry, solid) Coulomb friction, however, with a negative friction parameter. Diffusion coefficients are calculated by formal mapping to the situation of a quantum-mechanical harmonic oscillator exposed to an additional repulsive delta-potential. Results comprise a ditched or double-peaked velocity distribution and spatial statistics showing outward propagating maxima when starting from initially concentrated arrangements.

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.

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.

Dynamical anomalies and structural features of Active Brownian Particles characterised by two repulsive length scales

In this work we study a two-dimensional system composed by Active Brownian Particles (ABPs) interacting via a repulsive potential with two-length-scales, a soft shell and a hard-core. Depending on the ratio between the strength of the soft shell barrier and the activity, we find two regimes: If this ratio is much larger or smaller than 1, the observed behavior is comparable with ABPs interacting via a single length-scale potential. If this ratio is similar to 1, the two length-scales are relevant for both structure and dynamical properties. On the structural side, when the system exhibits a motility induced phase separation, the dense phase is characterised by new and more complex structures compared with the hexatic phase observed in single length-scale systems.On the dynamical side, as far as we are aware, this is the first representation of an anomalous dynamics in active particles.

Activity-driven phase separation and ordering kinetics of passive particles

The steady state and phase ordering kinetics in a pure active Brownian particle system are studied in recent years. In binary mixture of active and passive Brownian particles passive particles are used as probe to understand the properties of active medium. In our present study, we study the mixture of passive and active Brownian particles. Here, we aim to understand the steady state and kinetics of small passive particles in the mixture. In our system, the passive particles are small in size and large in number, whereas ABPs are large in size and small in number. The system is studied on a two-dimensional substrate using overdamped Langevin dynamic simulation. The steady state and kinetics of passive particles are studied for various size and activity of active particles. Passive particles are purely athermal in nature and have dynamics only due to bigger ABPs. For small size ratio and activity, the passive particles remain homogeneous in the system, whereas on increasing size ratio and activity they form periodic hexagonal close pack (HCP) spanning clusters in the system. We have also studied the kinetics of growing passive particle clusters. The mass of the largest cluster shows a much slower growth kinetics in contrast to conserved growth kinetics in ABP system. Our study provides an understanding of steady state and kinetics of passive particles in the presence of bigger active particles. The mixture can be thought of as effect of big microorganism moving in passive medium.

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.

Tracer dynamics in crowded active-particle suspensions

We discuss the dynamics of active Brownian particles (ABPs) in crowded environments through the mean-squared displacement (MSD) of active and passive tracer particles in both active and passive host systems. Exact equations for the MSD are derived using a projection operator technique, extending to dense systems the known solution for a single ABP. The interaction of the tracer particle with the host particles gives rise to strong memory effects. Evaluating these approximately in the framework of a recently developed mode-coupling theory for active Brownian particles (ABP-MCT), we discuss the various dynamical regimes that emerge: While self-propelled motion gives rise to super-diffusive MSD, at high densities, this competes with an interaction-induced sub-diffusive regime. The predictions of the theory are shown to be in good agreement with results obtained from an event-driven Brownian dynamics (ED-BD) simulation scheme for the dynamics of two-dimensional active Brownian hard disks.

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.

Effect of polydispersity on the dynamics of active Brownian particles

We numerically study the dynamics and the phases of self-propelled disk-shaped particles of different sizes with soft repulsive potential in two dimensions. Size diversity is introduced by the polydispersity index (PDI) ε, which is the width of the uniform distribution of the particle's radius. The self-propulsion speed of the particles controls the activity v. We observe enhanced dynamics for large size diversity among the particles. We calculate the effective diffusion coefficient D_{eff} in the steady state. The system exhibits four distinct phases, jammed phase with small D_{eff} for small activity and liquid phase with enhanced D_{eff} for large activity. The number fluctuation is larger and smaller than the equilibrium limit in the liquid and jammed phases, respectively. Further, the jammed phase is of two types: solid jammed and liquid jammed for small and large PDI. Whereas the liquid phase is called motility induced phase separation (MIPS) liquid for small PDI and for large PDI, we find enhanced diffusivity and call it the pure liquid phase. The system is studied for three packing densities ϕ, and the response of the system for polydispersity is the same for all ϕ's. Our study can help understand the behavior of cells of various sizes in a tissue, artificial self-driven granular particles, or living organisms of different sizes in a dense environment.

Phase Diagram of Active Brownian Spheres: Crystallization and the Metastability of Motility-Induced Phase Separation

Motility-induced phase separation (MIPS), the phenomenon in which purely repulsive active particles undergo a liquid-gas phase separation, is among the simplest and most widely studied examples of a nonequilibrium phase transition. Here, we show that states of MIPS coexistence are in fact only metastable for three-dimensional active Brownian particles over a very broad range of conditions, decaying at long times through an ordering transition we call active crystallization. At an activity just above the MIPS critical point, the liquid-gas binodal is superseded by the crystal-fluid coexistence curve, with solid, liquid, and gas all coexisting at the triple point where the two curves intersect. Nucleating an active crystal from a disordered fluid, however, requires a rare fluctuation that exhibits the nearly close-packed density of the solid phase. The corresponding barrier to crystallization is surmountable on a feasible timescale only at high activity, and only at fluid densities near maximal packing. The glassiness expected for such dense liquids at equilibrium is strongly mitigated by active forces, so that the lifetime of liquid-gas coexistence declines steadily with increasing activity, manifesting in simulations as a facile spontaneous crystallization at extremely high activity.

Characterization of MIPS in a suspension of repulsive active Brownian particles through dynamical features

We study a two-dimensional system composed by Active Brownian Particles (ABPs), focusing on the onset of Motility Induced Phase Separation (MIPS), by means of molecular dynamics simulations. For a pure hard-disk system with no translational diffusion, the phase diagram would be completely determined by their density and Péclet number. In our model, two additional effects are present: translational noise and the overlap of particles; we study the effects of both in the phase space. As we show, the second effect can be mitigated if we use, instead of the standard Weeks-Chandler-Andersen potential, a stiffer potential: the pseudo-hard sphere potential. Moreover, in determining the boundary of our phase space, we explore different approaches to detect MIPS and conclude that observing dynamical features, via the non-Gaussian parameter, is more efficient than observing structural ones, such as through the local density distribution function. We also demonstrate that the Vogel-Fulcher equation successfully reproduces the decay of the diffusion as a function of density, with the exception of very high densities. Thus, in this regard, the ABP system behaves similar to a fragile glass.

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.

Phase Transitions of Oppositely Charged Colloidal Particles Driven by Alternating Current Electric Field

We study systems containing oppositely charged colloidal particles under applied alternating current electric fields (AC fields) using overdamped Langevin dynamics simulations in three dimensions. We obtain jammed bands perpendicular to the field direction under intermediate frequencies and lanes parallel with the field under low frequencies. These structures also depend upon the particle charges. The pathway for generating jammed bands follows a stepwise mechanism, and intermediate bands are observed during lane formation in some systems. We investigate the component of the pressure tensors in the direction parallel to the field and observe that the jammed to lane transition occurs at a critical value for this pressure. We also find that the stable steady states appear to satisfy the principle of maximum entropy production. Our results may help to improve the understand of the underlying mechanisms for these types of dynamic phase transitions and the subsequent cooperative assemblies of colloidal particles under such non-equilibrium conditions.

Autonomously Probing Viscoelasticity in Disordered Suspensions

Recent experiments show a strong rotational diffusion enhancement for self-propelled microrheological probes in colloidal glasses. Here, we provide microscopic understanding using simulations with a frictional probe-medium coupling that converts active translation into rotation. Diffusive enhancement emerges from the medium's disordered structure and peaks at a second-order transition in the number of contacts. Our results reproduce the salient features of the colloidal glass experiment and support an effective description that is applicable to a broader class of viscoelastic suspensions.

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.

Universal reshaping of arrested colloidal gels via active doping

Colloids that interact via a short-range attraction serve as the primary building blocks for a broad range of self-assembled materials. However, one of the well-known drawbacks to this strategy is that these building blocks rapidly and readily condense into a metastable colloidal gel. Using computer simulations, we illustrate how the addition of a small fraction of purely repulsive self-propelled colloids, a technique referred to as active doping, can prevent the formation of this metastable gel state and drive the system toward its thermodynamically favored crystalline target structure. The simplicity and robust nature of this strategy offers a systematic and generic pathway to improving the self-assembly of a large number of complex colloidal structures. We discuss in detail the process by which this feat is accomplished and provide quantitative metrics for exploiting it to modulate the self-assembly. We provide evidence for the generic nature of this approach by demonstrating that it remains robust under a number of different anisotropic short-ranged pair interactions in both two and three dimensions. In addition, we report on a novel microphase in mixtures of passive and active colloids. For a broad range of self-propelling velocities, it is possible to stabilize a suspension of fairly monodisperse finite-size crystallites. Surprisingly, this microphase is also insensitive to the underlying pair interaction between building blocks. The active stabilization of these moderately sized monodisperse clusters is quite remarkable and should be of great utility in the design of hierarchical self-assembly strategies. This work further bolsters the notion that active forces can play a pivotal role in directing colloidal self-assembly.

Transport of self-propelled particles across a porous medium: trapping, clogging, and the Matthew effect

We study the transport of self-propelled particles from one free chamber to another across two stripe-like areas of dense porous medium. The medium is mimicked by arrays of obstacles. We find that active motion could greatly speed up the transport of particles. However, more and more particles become trapped in the obstacle arrays with the enhancement of activity. At high persistence (low rotational diffusion rate) and moderate particle concentration, we observe the Matthew effect in the aggregation of particles in the two obstacle arrays. This effect is weakened by introduction of randomness or deformability into the obstacle arrays. Moreover, the dependence on deformability shows the characteristics of first-order phase transition. In rare situations, the system could be stuck in a dynamic unstable state, e.g. the particles alternatively gather more in one of the two obstacle arrays, exhibiting oscillation of particle number between the arrays. Our results reveal new features in the transport of active objects in a complex medium and have implications for manipulating their collective assembly.

Anisotropic Exclusion Effect between Photocatalytic Ag/AgCl Janus Particles and Passive Beads in a Dense Colloidal Matrix

Synthetic nano- and micromotors interact with each other and their surroundings in a complex manner. Here, we report on the anisotropy of active-passive particle interaction in a soft matter system containing an immobile yet photochemical Ag/AgCl-based Janus particle embedded in a dense matrix of passive beads in pure water. The asymmetry in the chemical gradient around the Janus particle, triggered upon visible light illumination, distorts the isotropy of the surrounding electric potential and results in the repulsion of adjacent passive beads to a certain distance away from the Janus particle. This exclusion effect is found to be anisotropic with larger distances to passive beads in front of the Ag/AgCl cap of the Janus particle. We provide insight into this phenomenon by performing the angular analysis of the radii of exclusion and tracking their time evolution at the level of a single bead. Our study provides a novel fundamental insight into the collective behavior of a complex mixture of active and passive particles and is relevant for various application scenarios, e.g., particle transport at micro- and nanoscale and local chemical sensing.

Active colloidal molecules assembled via selective and directional bonds

The assembly of active and self-propelled particles is an emerging strategy to create dynamic materials otherwise impossible. However, control of the complex particle interactions remains challenging. Here, we show that various dynamic interactions of active patchy particles can be orchestrated by tuning the particle size, shape, composition, etc. This capability is manifested in establishing dynamic colloidal bonds that are highly selective and directional, which greatly expands the spectrum of colloidal structures and dynamics by assembly. For example, we demonstrate the formation of colloidal molecules with tunable bond angles and orientations. They exhibit controllable propulsion, steering, reconfiguration as well as other dynamic behaviors that collectively reflect the bond properties. The working principle is further extended to the co-assembly of synthetic particles with biological entities including living cells, giving rise to hybrid colloidal molecules of various types, for example, a colloidal carrousel structure. Our strategy should enable active systems to perform sophisticated tasks in future such as selective cell treatment.

The assembly of active and self-propelled particles is an emerging strategy to create dynamic materials otherwise impossible. Here, the authors show the assembly of active colloidal molecules with a wide spectrum of new structures and dynamics, conferred to them by highly selective and directional interactions.

Extreme active matter at high densities

We study the remarkable behaviour of dense active matter comprising self-propelled particles at large Péclet numbers, over a range of persistence times, from τ_p → 0, when the active fluid undergoes a slowing down of density relaxations leading to a glass transition as the active propulsion force f reduces, to τ_p → ∞, when as f reduces, the fluid jams at a critical point, with stresses along force-chains. For intermediate τ_p, a decrease in f drives the fluid through an intermittent phase before dynamical arrest at low f. This intermittency is a consequence of periods of jamming followed by bursts of plastic yielding associated with Eshelby deformations. On the other hand, an increase in f leads to an increase in the burst frequency; the correlated plastic events result in large scale vorticity and turbulence. Dense extreme active matter brings together the physics of glass, jamming, plasticity and turbulence, in a new state of driven classical matter.

While active matter exhibits unusual dynamics at low density, high density behavior has not been explored. Mandal et al. show that extreme dense active matter, shows a rich spectrum of behaviour from intermittent plastic bursts and turbulence, to glassy states and jamming in the limit of infinite persistence time.

Tuning the Glass Transition: Enhanced Crystallization of the Laves Phases in Nearly Hard Spheres

Colloidal crystals with a diamond and pyrochlore structure display wide photonic band gaps at low refractive index contrasts. However, these low-coordinated and open structures are notoriously difficult to self-assemble from colloids interacting with simple pair interactions. To circumvent these problems, one can self-assemble both structures in a closely packed MgCu₂ Laves phase from a binary mixture of colloidal spheres and then selectively remove one of the sublattices. Although Laves phases have been proven to be stable in a binary hard-sphere system, they have never been observed to spontaneously crystallize in such a fluid mixture in simulations nor in experiments of micron-sized hard spheres due to slow dynamics. Here we demonstrate, using computer simulations, that softness in the interparticle potential suppresses the degree of 5-fold symmetry in the binary fluid phase and enhances crystallization of Laves phases in nearly hard spheres.

Activity-modulated phase transition in a two-dimensional mixture of active and passive colloids

We study a two-dimensional binary mixture of active and passive colloids as an idealized model of an hybrid aggregate of living cells and inert particles. We perform molecular dynamics simulations of this system using two different thermostats, and we systematically investigate the effect of varying these two effective temperatures on the system behavior, as characterized by its density, structure and thermoelastic properties. Our results indicate that the presence of active colloids shifts the mixture towards the liquid state and renders it more deformable. Such system softening and melting effects due to the addition of active particles are larger than expected from a linear combination of temperatures of the active and passive components. This heightened effect becomes more pronounced as the effective temperature difference between the two components becomes larger. The binary mixture remains homogeneous for moderate colloidal activity, but segregation arises for large effective temperature difference. Our results provide insights to guide future experimental hybrid aggregate studies with promising biomedical applications.

Direct-imaging of light-driven colloidal Janus particles in weightlessness

We describe a highly integrated automated experiment module that allows us to investigate the active Brownian motion of light-driven colloidal Janus-particle suspensions. The module RAMSES (RAndom Motion of SElf-propelled particles in Space) is designed for the sounding rocket platform MAPHEUS (MAterialPHysikalische Experimente Unter Schwerelosigkeit). It allows us to perform experiments under weightlessness conditions in order to avoid sedimentation of the Janus particles and thus to study the spatially three-dimensional dynamics in the suspension. The module implements a newly developed strong homogeneous light source to excite self-propulsion in the Janus particles. The light source is realized through an array of high-power light-emitting diodes and replaces the conventional laser source, thus reducing heat dissipation and spatial extension of the experiment setup. The rocket module contains ten independent sample cells in order to ease the systematic study of the effect of control parameters such as light intensity or particle concentration and size in a single sounding-rocket flight. For each sample cell, transmitted light intensities are stored for postflight analysis in terms of differential dynamical microscopy.

Active Hard Spheres in Infinitely Many Dimensions

Few equilibrium-even less so nonequilibrium-statistical-mechanical models with continuous degrees of freedom can be solved exactly. Classical hard spheres in infinitely many space dimensions are a notable exception. We show that, even without resorting to a Boltzmann distribution, dimensionality is a powerful organizing device for exploring the stationary properties of active hard spheres evolving far from equilibrium. In infinite dimensions, we exactly compute the stationary state properties that govern and characterize the collective behavior of active hard spheres: the structure factor and the equation of state for the pressure. In turn, this allows us to account for motility-induced phase separation. Finally, we determine the crowding density at which the effective propulsion of a particle vanishes.