Dynamic Gardner cross-over in a simple glass.

The criticality of the jamming transition responsible for amorphous solidification has been theoretically linked to the marginal stability of a thermodynamic Gardner phase. While the critical exponents of jamming appear independent of the preparation history, the pertinence of Gardner physics far from equilibrium is an open question. To fill this gap, we numerically study the nonequilibrium dynamics of hard disks compressed toward the jamming transition using a broad variety of protocols. We show that dynamic signatures of Gardner physics can be disentangled from the aging relaxation dynamics. We thus define a generic dynamic Gardner cross-over regardless of the history. Our results show that the jamming transition is always accessed by exploring increasingly complex landscape, resulting in anomalous microscopic relaxation dynamics that remains to be understood theoretically.

Two-step devitrification of ultrastable glasses.

The discovery of ultrastable glasses raises novel challenges about glassy systems. Recent experiments studied the macroscopic devitrification of ultrastable glasses into liquids upon heating but lacked microscopic resolution. We use molecular dynamics simulations to analyze the kinetics of this transformation. In the most stable systems, devitrification occurs after a very large time, but the liquid emerges in two steps. At short times, we observe the rare nucleation and slow growth of isolated droplets containing a liquid maintained under pressure by the rigidity of the surrounding glass. At large times, pressure is released after the droplets coalesce into large domains, which accelerates devitrification. This two-step process produces pronounced deviations from the classical Avrami kinetics and explains the emergence of a giant lengthscale characterizing the devitrification of bulk ultrastable glasses. Our study elucidates the nonequilibrium kinetics of glasses following a large temperature jump, which differs from both equilibrium relaxation and aging dynamics, and will guide future experimental studies.

Direct Numerical Analysis of Dynamic Facilitation in Glass-Forming Liquids.

We propose a computational strategy to quantify the temperature evolution of the timescales and length scales over which dynamic facilitation affects the relaxation dynamics of glass-forming liquids at low temperatures, which requires no assumption about the nature of the dynamics. In two glass models, we find that dynamic facilitation depends strongly on temperature, leading to a subdiffusive spreading of relaxation events which we characterize using a temperature-dependent dynamic exponent. We also establish that this temperature evolution represents a major contribution to the increase of the structural relaxation time.

Collective Relaxation Dynamics in a Three-Dimensional Lattice Glass Model.

We numerically elucidate the microscopic mechanisms controlling the relaxation dynamics of a three-dimensional lattice glass model that has static properties compatible with the approach to a random first-order transition. At low temperatures, the relaxation is triggered by a small population of particles with low-energy barriers forming mobile clusters. These emerging quasiparticles act as facilitating defects responsible for the spatially heterogeneous dynamics of the system, whose characteristic length scales remain strongly coupled to thermodynamic fluctuations. We compare our findings both with existing theoretical models and atomistic simulations of glass formers.

Irreversible Monte Carlo Algorithms for Hard Disk Glasses: From Event-Chain to Collective Swaps.

Equilibrium sampling of the configuration space in disordered systems requires algorithms that bypass the glassy slowing down of the physical dynamics. Irreversible Monte Carlo algorithms breaking detailed balance successfully accelerate sampling in some systems. We first implement an irreversible event-chain Monte Carlo algorithm in a model of continuously polydisperse hard disks. The effect of collective translational moves marginally affects the dynamics and results in a modest speedup that decreases with density. We then propose an irreversible algorithm performing collective particle swaps which outperforms all known Monte Carlo algorithms. We show that these collective swaps can also be used to prepare very dense jammed packings of disks.

Emerging Mesoscale Flows and Chaotic Advection in Dense Active Matter.

We study two models of overdamped self-propelled disks in two dimensions, with and without aligning interactions. Both models support active mesoscale flows, leading to chaotic advection and transport over large length scales in their homogeneous dense fluid states, away from dynamical arrest. They form streams and vortices reminiscent of multiscale flow patterns in turbulence. We show that the characteristics of these flows do not depend on the specific details of the active fluids, and result from the competition between crowding effects and persistent propulsions. This observation suggests that dense active suspensions of self-propelled particles present a type of "active turbulence" distinct from collective flows reported in other types of active 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.

Predicting Dynamic Heterogeneity in Glass-Forming Liquids by Physics-Inspired Machine Learning.

We introduce GlassMLP, a machine learning framework using physics-inspired structural input to predict the long-time dynamics in deeply supercooled liquids. We apply this deep neural network to atomistic models in 2D and 3D. Its performance is better than the state of the art while being more parsimonious in terms of training data and fitting parameters. GlassMLP quantitatively predicts four-point dynamic correlations and the geometry of dynamic heterogeneity. Transferability across system sizes allows us to efficiently probe the temperature evolution of spatial dynamic correlations, revealing a profound change with temperature in the geometry of rearranging regions.

Two-Dimensional Crystals far from Equilibrium.

When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behavior with no equilibrium counterpart. We study a two-dimensional particle model initially proposed to describe driven non-Brownian suspensions undergoing nonequilibrium absorbing phase transitions. We show that when the transition occurs at large density, the dynamics produces long-range crystalline order. In the ordered phase, long-range translational order is observed because equipartition of energy is lacking, phonons are suppressed, and density fluctuations are hyperuniform. Our study offers an explicit microscopic model where nonequilibrium violations of the Mermin-Wagner theorem stabilize crystalline order in two dimensions.

Sampling Efficiency of Transverse Forces in Dense Liquids.

Sampling the Boltzmann distribution using forces that violate detailed balance can be faster than with the equilibrium evolution, but the acceleration depends on the nature of the nonequilibrium drive and the physical situation. Here, we study the efficiency of forces transverse to energy gradients in dense liquids through a combination of techniques: Brownian dynamics simulations, exact infinite-dimensional calculation, and a mode-coupling approximation. We find that the sampling speedup varies nonmonotonically with temperature, and decreases as the system becomes more glassy. We characterize the interplay between the distance to equilibrium and the efficiency of transverse forces by means of odd transport coefficients.

Intermittent relaxation and avalanches in extremely persistent active matter.

We use numerical simulations to study the dynamics of dense assemblies of self-propelled particles in the limit of extremely large, but finite, persistence times. In this limit, the system evolves intermittently between mechanical equilibria where active forces balance interparticle interactions. We develop an efficient numerical strategy allowing us to resolve the statistical properties of elastic and plastic relaxation events caused by activity-driven fluctuations. The system relaxes via a succession of scale-free elastic events and broadly distributed plastic events that both depend on the system size. Correlations between plastic events lead to emergent dynamic facilitation and heterogeneous relaxation dynamics. Our results show that dynamical behaviour in extremely persistent active systems is qualitatively similar to that of sheared amorphous solids, yet with some important differences.

Front propagation in ultrastable glasses is dynamically heterogeneous.

Upon heating, ultrastable glassy films transform into liquids via a propagating equilibration front, resembling the heterogeneous melting of crystals. A microscopic understanding of this robust phenomenology is, however, lacking because experimental resolution is limited. We simulate the heterogeneous transformation kinetics of ultrastable configurations prepared using the swap Monte Carlo algorithm, thus allowing a direct comparison with experiments. We resolve the liquid-glass interface both in space and in time as well as the underlying particle motion responsible for its propagation. We perform a detailed statistical analysis of the interface geometry and kinetics over a broad range of temperatures. We show that the dynamic heterogeneity of the bulk liquid is passed on to the front that propagates heterogeneously in space and intermittently in time. This observation allows us to relate the averaged front velocity to the equilibrium diffusion coefficient of the liquid. We suggest that an experimental characterization of the interface geometry during the heterogeneous devitrification of ultrastable glassy films could provide direct experimental access to the long-sought characteristic length scale of dynamic heterogeneity in bulk supercooled liquids.

Microscopic observation of two-level systems in a metallic glass model.

The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements.

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.

Static self-induced heterogeneity in glass-forming liquids: Overlap as a microscope.

We propose and numerically implement a local probe of the static self-induced heterogeneity characterizing glass-forming liquids. This method relies on the equilibrium statistics of the overlap between pairs of configurations measured in mesoscopic cavities with unconstrained boundaries. By systematically changing the location of the probed cavity, we directly detect spatial variations of the overlap fluctuations. We provide a detailed analysis of the statistics of a local estimate of the configurational entropy, and we infer an estimate of the surface tension between amorphous states, ingredients that are both at the basis of the random first-order transition theory of glass formation. Our results represent the first direct attempt to visualize and quantify the self-induced heterogeneity underpinning the thermodynamics of glass formation. They pave the way for the development of coarse-grained effective theories and for a direct assessment of the role of thermodynamics in the activated dynamics of deeply supercooled liquids.

Collective dynamics in a glass-former with Mari-Kurchan interactions.

We numerically study the equilibrium relaxation dynamics of a two-dimensional Mari-Kurchan glass model. The tree-like structure of particle interactions forbids both nontrivial structural motifs and the emergence of a complex free-energy landscape leading to a thermodynamic glass transition, while the finite-dimensional nature of the model prevents the existence of a mode-coupling singularity. Nevertheless, the equilibrium relaxation dynamics is shown to be in excellent agreement with simulations performed in conventional glass-formers. Averaged time correlation functions display a phenomenology typical of supercooled liquids, including the emergence of an excess signal in relaxation spectra at intermediate frequencies. We show that this evolution is accompanied by strong signatures of collective and heterogeneous dynamics that cannot be interpreted in terms of single particle hopping and emerge from dynamic facilitation. Our study demonstrates that an off-lattice interacting particle model with extremely simple structural correlations displays quantitatively realistic glassy dynamics.

Self-Induced Heterogeneity in Deeply Supercooled Liquids.

A theoretical treatment of deeply supercooled liquids is difficult because their properties emerge from spatial inhomogeneities that are self-induced, transient, and nanoscopic. I use computer simulations to analyze self-induced static and dynamic heterogeneity in equilibrium systems approaching the experimental glass transition. I characterize the broad sample-to-sample fluctuations of salient dynamic and thermodynamic properties in elementary mesoscopic systems. Findings regarding local lifetimes and distributions of dynamic heterogeneity are in excellent agreement with recent single molecule studies. Surprisingly broad thermodynamic fluctuations are also found, which correlate well with dynamic fluctuations, thus providing a local test of the thermodynamic origin of slow dynamics.

Excess wings and asymmetric relaxation spectra in a facilitated trap model.

In a recent computer study, we have shown that the combination of spatially heterogeneous dynamics and kinetic facilitation provides a microscopic explanation for the emergence of excess wings in deeply supercooled liquids. Motivated by these findings, we construct a minimal empirical model to describe this physics and introduce dynamic facilitation in the trap model, which was initially developed to capture the thermally activated dynamics of glassy systems. We fully characterize the relaxation dynamics of this facilitated trap model varying the functional form of energy distributions and the strength of dynamic facilitation, combining numerical results and analytic arguments. Dynamic facilitation generically accelerates the relaxation of the deepest traps, thus making relaxation spectra strongly asymmetric, with an apparent "excess" signal at high frequencies. For well-chosen values of the parameters, the obtained spectra mimic experimental results for organic liquids displaying an excess wing. Overall, our results identify the minimal physical ingredients needed to describe excess processes in the relaxation spectra of supercooled liquids.

Random critical point separates brittle and ductile yielding transitions in amorphous materials.

We combine an analytically solvable mean-field elasto-plastic model with molecular dynamics simulations of a generic glass former to demonstrate that, depending on their preparation protocol, amorphous materials can yield in two qualitatively distinct ways. We show that well-annealed systems yield in a discontinuous brittle way, as metallic and molecular glasses do. Yielding corresponds in this case to a first-order nonequilibrium phase transition. As the degree of annealing decreases, the first-order character becomes weaker and the transition terminates in a second-order critical point in the universality class of an Ising model in a random field. For even more poorly annealed systems, yielding becomes a smooth crossover, representative of the ductile rheological behavior generically observed in foams, emulsions, and colloidal glasses. Our results show that the variety of yielding behaviors found in amorphous materials does not necessarily result from the diversity of particle interactions or microscopic dynamics but is instead unified by carefully considering the role of the initial stability of the system.

Glassy Behavior of Sticky Spheres: What Lies beyond Experimental Timescales?

We use the swap Monte Carlo algorithm to analyze the glassy behavior of sticky spheres in equilibrium conditions at densities where conventional simulations and experiments fail to reach equilibrium, beyond predicted phase transitions and dynamic singularities. We demonstrate the existence of a unique ergodic region comprising all the distinct phases previously reported, except for a phase-separated region at strong adhesion. All structural and dynamic observables evolve gradually within this ergodic region, the physics evolving smoothly from well-known hard sphere glassy behavior at small adhesions and large densities, to a more complex glassy regime characterized by unusually broad distributions of relaxation timescales and length scales at large adhesions.