Monte Carlo simulations of glass-forming liquids beyond Metropolis.

Monte Carlo simulations are widely employed to measure the physical properties of glass-forming liquids in thermal equilibrium. Combined with local Monte Carlo moves, the Metropolis algorithm can also be used to simulate the relaxation dynamics, thus offering an efficient alternative to molecular dynamics. Monte Carlo simulations are, however, more versatile because carefully designed Monte Carlo algorithms can more efficiently sample the rugged free energy landscape characteristic of glassy systems. After a brief overview of Monte Carlo studies of glass-formers, we define and implement a series of Monte Carlo algorithms in a three-dimensional model of polydisperse hard spheres. We show that the standard local Metropolis algorithm is the slowest and that implementing collective moves or breaking detailed balance enhances the efficiency of the Monte Carlo sampling. We use time correlation functions to provide a microscopic interpretation of these observations. Seventy years after its invention, the Monte Carlo method remains the most efficient and versatile tool to compute low-temperature properties in supercooled liquids.

Transverse forces and glassy liquids in infinite dimensions.

We explore the dynamics of a simple liquid whose particles, in addition to standard potential-based interactions, are also subjected to transverse forces preserving the Boltzmann distribution. We derive the effective dynamics of one and two tracer particles in the infinite-dimensional limit. We determine the amount of acceleration of the dynamics caused by the transverse forces, in particular in the vicinity of the glass transition. We analyze the emergence and evolution of odd transport phenomena induced by the transverse forces.

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.

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.

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.

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.

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.

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.

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.

Finding defects in glasses through machine learning.

Structural defects control the kinetic, thermodynamic and mechanical properties of glasses. For instance, rare quantum tunneling two-level systems (TLS) govern the physics of glasses at very low temperature. Due to their extremely low density, it is very hard to directly identify them in computer simulations. We introduce a machine learning approach to efficiently explore the potential energy landscape of glass models and identify desired classes of defects. We focus in particular on TLS and we design an algorithm that is able to rapidly predict the quantum splitting between any two amorphous configurations produced by classical simulations. This in turn allows us to shift the computational effort towards the collection and identification of a larger number of TLS, rather than the useless characterization of non-tunneling defects which are much more abundant. Finally, we interpret our machine learning model to understand how TLS are identified and characterized, thus giving direct physical insight into their microscopic nature.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.