Search

The ultraslow dynamics of glass-formers has been explained by two views considered as mutually exclusive: one invokes locally hindered mobility, the other rests on the complexity of the configuration space. Here we demonstrate that the evolution responds strongly to the details of the dynamics by changing the speed of time-flow: it has time-reparametrization softness. This finding reconciles both views: while local constraints reparametrize the flow of time, the global landscape determines relationships between different correlations at the same times. We show that modern algorithms developed to accelerate the relaxation to equilibrium act by changing the time reparametrization. Their success thus relies on their ability to exploit reparametrization softness. We conjecture that these results extend beyond the realm of glasses to the optimization of more general constraint satisfaction problems and to broader classes of algorithms.

The swap Monte Carlo algorithm introduces non-physical dynamic rules to accelerate the exploration of the configuration space of supercooled liquids. Its success raises deep questions regarding the nature and physical origin of the slow dynamics of dense liquids, and how it is affected by swap moves. We provide a detailed analysis of the slow dynamics generated by the swap Monte Carlo algorithm at very low temperatures in two glass-forming models. We find that the slowing down of the swap dynamics is qualitatively distinct from its local Monte Carlo counterpart, with considerably suppressed dynamic heterogeneity both at single-particle and collective levels. Our results suggest that local kinetic constraints are drastically reduced by swap moves, leading to nearly Gaussian and diffusive dynamics and weakly growing dynamic correlation lengthscales. The comparison between static and dynamic fluctuations shows that swap Monte Carlo is a nearly optimal local equilibrium algorithm, suggesting that further progress should necessarily involve collective or driven algorithms., Comment: 15 pages, 13 figures

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-temperatures properties in supercooled liquids., Comment: 12 pages, 13 figures. Submitted to "Special Topic on Monte Carlo Methods, 70 years after Metropolis et al. (1953)" in J. Chem. Phys

Exerting a nonequilibrium drive on an otherwise equilibrium Langevin process brings the dynamics out of equilibrium but can also speedup the approach to the Boltzmann steady-state. Transverse forces are a minimal framework to achieve dynamical acceleration of the Boltzmann sampling. We consider a simple liquid in three space dimensions subjected to additional transverse pairwise forces, and quantify the extent to which transverse forces accelerate the dynamics. We first explore the dynamics of a tracer in a weak coupling regime describing high temperatures. The resulting acceleration is correlated with a monotonous increase of the magnitude of odd transport coefficients (mobility and diffusivity) with the amplitude of the transverse drive. We then develop a nonequilibrium version of the mode-coupling theory able to capture the effect of transverse forces, and more generally of forces created by additional degrees of freedom. Based on an analysis of transport coefficients, both odd and longitudinal, both for the collective modes and for a tracer particle, we find a systematic acceleration of the dynamics. Quantitatively, the gain, which is guaranteed throughout the ergodic phase, turns out to be a decreasing function of temperature beyond a temperature crossover, in particular as the glass transition is approached. Our theoretical results are in good agreement with available numerical results.

Normalizing flows can transform a simple prior probability distribution into a more complex target distribution. Here, we evaluate the ability and efficiency of generative machine learning methods to sample the Boltzmann distribution of an atomistic model for glass-forming liquids. This is a notoriously difficult task, as it amounts to ergodically exploring the complex free energy landscape of a disordered and frustrated many-body system. We optimize a normalizing flow model to successfully transform high-temperature configurations of a dense liquid into low-temperature ones, near the glass transition. We perform a detailed comparative analysis with established enhanced sampling techniques developed in the physics literature to assess and rank the performance of normalizing flows against state-of-the-art algorithms. We demonstrate that machine learning methods are very promising, showing a large speedup over conventional molecular dynamics. Normalizing flows show performances comparable to parallel tempering and population annealing, while still falling far behind the swap Monte Carlo algorithm. Our study highlights the potential of generative machine learning models in scientific computing for complex systems, but also points to some of its current limitations and the need for further improvement.

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.

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 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.

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., Comment: 21 pages, 7 figures, perspective on the workshop 'Yield stress and fluidization in brittle and ductile amorphous systems' held from January 9 to 13, 2023 at the Lorentz Center in Leiden

Unraveling the connections between microscopic structure, emergent physical properties, and slow dynamics has long been a challenge when studying the glass transition. The absence of clear visible structural order in amorphous configurations complicates the identification of the key physical mechanisms underpinning slow dynamics. The difficulty in sampling equilibrated configurations at low temperatures hampers thorough numerical and theoretical investigations. This perspective article explores the potential of machine learning (ML) techniques to face these challenges, building on the algorithms that have revolutionized computer vision and image recognition. We present recent successful ML applications, as well as many open problems for the future, such as transferability and interpretability of ML approaches. We highlight new ideas and directions in which ML could provide breakthroughs to better understand the fundamental mechanisms at play in glass-forming liquids. To foster a collaborative community effort, this article also introduces the ``GlassBench" dataset, providing simulation data and benchmarks for both two-dimensional and three-dimensional glass-formers. We propose critical metrics to compare the performance of emerging ML methodologies, in line with benchmarking practices in image and text recognition. The goal of this roadmap is to provide guidelines for the development of ML techniques in systems displaying slow dynamics, while inspiring new directions to improve our theoretical understanding of glassy liquids.

We develop a transferable machine learning model which predicts structural relaxation from amorphous supercooled liquid structures. The trained networks are able to predict dynamic heterogeneity across a broad range of temperatures and time scales with excellent accuracy and transferability. We use the network transferability to predict dynamic heterogeneity down to the experimental glass transition temperature, $T_g$, where structural relaxation cannot be analyzed using molecular dynamics simulations. The results indicate that the strength, the geometry and the characteristic length scale of the dynamic heterogeneity evolve much more slowly near $T_g$ compared to their evolution at higher temperatures. Our results show that machine learning techniques can provide physical insights on the nature of the glass transition that cannot be gained using conventional simulation techniques.

We propose a computational strategy to quantify the temperature evolution of the timescales and lengthscales over which dynamic facilitation affects the relaxation dynamics of glass-forming liquids at low temperatures, that 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.

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 lengthscales remain strongly coupled to thermodynamic fluctuations. We compare our findings both with existing theoretical models and atomistic simulations of glass-formers., Comment: 7 pages, 5 figures + supplemental material

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 non-monotonically 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.

We study two models of overdamped self-propelled disks in two dimensions, with and without aligning interactions. Active mesoscale flows leading to chaotic advection emerge in both models in the homogeneous dense fluid away from dynamical arrest, forming 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. Our results suggest that dense active fluids present a type of `active turbulence' distinct from collective flows reported in other types of active systems., Comment: 6 pages, 4 figs

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 direct comparison with experiments. We resolve the liquid-glass interface both in space and 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 which 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 characterisation of the interface geometry during the heterogeneous devitrification of ultrastable glassy films would provide direct experimental access to the long-sought characteristic lengthscale of dynamic heterogeneity in bulk supercooled liquids.

When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behaviour 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., Comment: 6 pages, 4 figures. V2 accepted for publication at Phys. Rev. Lett

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., Comment: Authors' accepted version for publication in Soft Matter

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. Because of 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.

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.

Comment on `Fickian Non-Gaussian Diffusion in Glass-Forming Liquids', by Rusciano et al., Phys. Rev. Lett. 128, 168001 (2022). In a recent Letter, Rusciano et al. examined the statistics of individual particles displacements in two-dimensional glass-formers and concluded that the corresponding probability distribution is non-Gaussian in a time regime where the mean-squared displacement is Fickian. Here, we clarify that the multiple length scales and time scales reported in this work have either been characterized before, or are not well-defined. This leads us to dispute the conclusions that glass-formers display Fickian non-Gaussian behaviour and that this analogy fruitfully addresses the central questions regarding the nature of dynamic heterogeneity in these systems., Comment: 2 pages, no figure

The discovery of ultrastable glasses has raised 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 analyse 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 characterising 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.

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 towards 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 crossover regardless of the history. Our results show that the jamming transition is always accessed by exploring increasingly complex landscape, resulting in the anomalous microscopic relaxation dynamics that remains to be understood theoretically.

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., Comment: 17 pages, 12 figures

The physics of the glass transition and amorphous materials continues to attract the attention of a wide research community after decades of effort. Supercooled liquids and glasses have been studied numerically since the advent of molecular dynamics and Monte Carlo simulations in the last century. Computer studies have greatly enhanced both experimental discoveries and theoretical developments and constitute an active and continually expanding research field. Our goal in this review is to provide a modern perspective on this area. We describe the need to go beyond canonical methods to attack a problem that is notoriously difficult in terms of time scales, length scales, and physical observables. We first summarise recent algorithmic developments to achieve enhanced sampling and faster equilibration using replica exchange methods, cluster and swap Monte Carlo algorithms, and other techniques. We then review some major recent advances afforded by these novel tools regarding the statistical mechanical description of the liquid-to-glass transition as well as the mechanical, vibrational and thermal properties of the glassy solid. We finally describe some important challenges for future research., Comment: 17 pages, 5 figures, to be published in Nature Reviews Physics

Glass is everywhere. We use and are surrounded by glass objects which make tangible the reality of glass as a distinct state of matter. Yet, glass as we know it is usually obtained by cooling a liquid sufficiently rapidly below its melting point to avoid crystallisation. The viscosity of this supercooled liquid increases by many orders of magnitude upon cooling, until the liquid becomes essentially arrested on experimental timescales below the ``glass transition'' temperature. From a structural viewpoint, the obtained glass still very much resembles the disordered liquid, but from a mechanical viewpoint, it is as rigid as an ordered crystal. Does glass qualify as a separate state of matter? We provide a pedagogical perspective on this question using basic statistical mechanical concepts. We recall the definitions of states of matter and of phase transitions between them. We review recent theoretical results suggesting why and how an ``ideal glass'' can indeed be defined as a separate equilibrium state of matter. We discuss recent success of computer simulations trying to analyse this glass state. We close with some experimental perspectives., Comment: 10 pages, 4 figures. Pedagogical article based on talk given at 16th International Conference on the Physics of Non-Crystalline Solids, Canterbury (2022)

We combine the swap Monte Carlo algorithm to long multi-CPU molecular dynamics simulations to analyse the equilibrium relaxation dynamics of model supercooled liquids over a time window covering ten orders of magnitude for temperatures down to the experimental glass transition temperature $T_g$. The analysis of \rev{several} time correlation functions coupled to spatio-temporal resolution of particle motion allow us to elucidate the nature of the equilibrium dynamics in deeply supercooled liquids. We find that structural relaxation starts at early times in rare localised regions characterised by a waiting time distribution that develops a power law near $T_g$. At longer times, relaxation events accumulate with increasing probability in these regions as $T_g$ is approached. This accumulation leads to a power-law growth of the linear extension of relaxed domains with time with a large, temperature-dependent dynamic exponent. Past the average relaxation time, unrelaxed domains slowly shrink with time due to relaxation events happening at their boundaries. Our results provide a complete microscopic description of the particle motion responsible for key experimental signatures of glassy dynamics, from the shape and temperature evolution of relaxation spectra to the core features of dynamic heterogeneity. They also provide a microscopic basis to understand the emergence of dynamic facilitation in deeply supercooled liquids and allow us to critically reassess theoretical descriptions of the glass transition., Comment: 22 pages, 21 figures

We provide an overview of the different types of computational techniques developed over the years to study supercooled liquids, glassy materials and the physics of the glass transition. We organise these numerical strategies into four broad families. For each of them, we describe the general ideas without discussing any technical details. We summarise the type of questions which can be addressed by any given approach and outline the main results which have been obtained. Finally we describe two important directions for future computational studies of glassy systems., Comment: 16 pages, 1 figure. Submitted to the French Academy of Sciences to celebrate International Year of Glass 2022

We numerically study the equilibrium relaxation dynamics of a two-dimensional Mari-Kurchan glass model. The tree-like structure of particle interactions forbids both non-trivial 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 which 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., Comment: 12 pages, 7 figures

We present numerical tests of the fluctuation-dissipation theorem (FDT) in the dumbbell model of spin ice with parameters suitable for dysprosium titanate. The tests are made for local spin variables, magnetic monopole density, and energy. We are able to achieve local equilibrium in which the FDT is satisfied down to $T=0.4$ K below which the system completely freezes. Non-equilibrium dynamics, together with violation of the FDT, are nonetheless observed following a thermal quench into the non-contractable monopole pair regime. Despite FDT violation, an approximate linear response regime allows for the identification of effective non-equilibrium temperatures which are different for each variable. The spin variable appears hotter than the heat reservoir, the monopole concentration responds with a lower effective temperature while the energy has a negative effective temperature. Results are discussed in the context of the monopole picture of spin ice and compared to the structure of FDT violations in other glassy materials. Prospectives for future experiments are reviewed., Comment: 15 pages; 9 figs

We propose and numerically implement a local probe of the static self-induced heterogeneity characterizing glass-forming liquids. The 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., Comment: 13 pages, 5 figures

We explore the emergence of nonequilibrium collective motion in disordered non-thermal 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, Comment: Authors' accepted version for publication in Physical Review Letters

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 characterise 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 relaxation spectra of supercooled liquids., Comment: 11 pages, 9 figures

We numerically study the zero-temperature relaxation dynamics of several glass-forming models to their inherent structures, following quenches from equilibrium configurations sampled across a wide range of initial temperatures. In a mean-field Mari-Kurchan model, we find that relaxation changes from a power-law to an exponential decay below a well-defined temperature, consistent with recent findings in mean-field $p$-spin models. By contrast, for finite-dimensional systems, the relaxation is always algebraic, with a non-trivial universal exponent at high temperatures crossing over to a harmonic value at low temperatures. We demonstrate that this apparent evolution is controlled by a temperature-dependent population of localised glassy excitations. Our work unifies several recent lines of studies aiming at a detailed characterisation of the complex potential energy landscape of glass-formers, and challenges both mean-field and real space descriptions of glasses., Comment: 17 pages, 10 figs. V2 accepted for publication at Phys Rev X

We study the statistical mechanics of supercooled liquids when the system evolves at a temperature $T$ with a field $\epsilon$ linearly coupled to its overlap with a reference configuration of the same liquid sampled at a temperature $T_0$. We use mean-field theory to fully characterize the influence of the reference temperature $T_0$, and we mainly study the case of a fixed, low-$T_0$ value in computer simulations. We numerically investigate the extended phase diagram in the $(\epsilon,T)$ plane of model glass-forming liquids in spatial dimensions $d=2$ and $d=3$, relying on umbrella sampling and reweighting techniques. For both $2d$ and $3d$ cases, a similar phenomenology with nontrivial thermodynamic fluctuations of the overlap is observed at low temperatures, but a detailed finite-size analysis reveals qualitatively distinct behaviors. We establish the existence of a first-order transition line for nonzero $\epsilon$ ending in a critical point in the universality class of the random-field Ising model (RFIM) in $d=3$. In $d=2$ instead, no phase transition is found in large enough systems at least down to temperatures below the extrapolated calorimetric glass transition temperature $T_g$. Our results confirm that glass-forming liquid samples of limited size display the thermodynamic fluctuations expected for finite systems undergoing a random first-order transition. They also support the relevance of the physics of the RFIM for supercooled liquids, which may then explain the qualitative difference between $2d$ and $3d$ glass-formers., Comment: Main text of 14 pages, 11 figures + two appendices of 14 pages, 5 figures. Accepted for publication in SciPost Physics

Glass formation is encountered in diverse materials. Experiments have revealed that dynamic relaxation spectra of supercooled liquids generically become asymmetric near the glass transition temperature, $T_g$, where an extended power law emerges at high frequencies. The microscopic origin of this "wing" remains unknown, and was so far inaccessible to simulations. Here, we develop a novel computational approach and study the equilibrium dynamics of model supercooled liquids near $T_g$. We demonstrate the emergence of a power law wing in numerical spectra, which originates from relaxation at rare, localised regions over broadly-distributed timescales. We rationalise the asymmetric shape of relaxation spectra by constructing an empirical model associating heterogeneous activated dynamics with dynamic facilitation, which are the two minimal physical ingredients revealed by our simulations. Our work offers a glimpse of the molecular motion responsible for glass formation at relevant experimental conditions., Comment: 5 pages, 4 figures, supplementary materials

We use atomistic computer simulations to provide a microscopic description of the brittle failure of amorphous materials, and we assess the role of rare events and quenched disorder. We argue that brittle yielding originates at rare soft regions, similarly to Griffiths effects in disordered systems. We numerically demonstrate how localized plastic events in such soft regions trigger macroscopic failure via the propagation of a shear band. This physical picture, which no longer holds in poorly annealed ductile materials, allows us to discuss the role of finite size effects in brittle yielding and reinforces the similarities between yielding and other disorder-controlled nonequilibrium phase transitions., Comment: 12 pages, 10 figures

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 analyse self-induced static and dynamic heterogeneity in equilibrium systems approaching the experimental glass transition. I characterise 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 dynamics fluctuations, thus providing a local test of the thermodynamic origin of slow dynamics., Comment: 6 pages, 4 figs

We use the swap Monte Carlo algorithm to analyse the glassy behaviour 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 behaviour at small adhesions and large densities, to a more complex glassy regime characterised by unusually-broad distributions of relaxation timescales and lengthscales at large adhesions., Comment: 6 pages, 4 figures

We numerically study the relaxation dynamics and associated criticality of non-Brownian frictionless spheres below jamming in spatial dimensions $d=2$, $3$, $4$, and $8$, and in the mean-field Mari-Kurchan model. We discover non-trivial finite-size and volume fraction dependences of the relaxation time associated to the relaxation of unjammed packings. In particular, the relaxation time is shown to diverge logarithmically with system size at any density below jamming, and no critical exponent can characterise its behaviour approaching jamming. In mean-field, the relaxation time is instead well-defined: it diverges at jamming with a critical exponent that we determine numerically and differs from an earlier mean-field prediction. We rationalise the finite $d$ logarithmic divergence using an extreme-value statistics argument in which the relaxation time is dominated by the most connected region of the system. The same argument shows that the earlier proposition that relaxation dynamics and shear viscosity are directly related breaks down in large systems. The shear viscosity of non-Brownian packings is well-defined in all $d$ in the thermodynamic limit, but large finite-size effects plague its measurement close to jamming., Comment: 11 pages, 8 figures

The overlap, or similarity, between liquid configurations is at the core of the mean-field description of the glass transition, and remains a useful concept when studying three-dimensional glass-forming liquids. In liquids, however, the overlap involves a tolerance, typically of a fraction $a/\sigma$ of the inter-particle distance, associated with how precisely similar two configurations must be for belonging to the same physically relevant "state". Here, we systematically investigate the dependence of the overlap fluctuations and of the resulting phase diagram when the tolerance is varied over a large range. We show that while the location of the dynamical and thermodynamic glass transitions (if present) is independent of $a/\sigma$, that of the critical point associated with a transition between a low- and a high-overlap phases in the presence of an applied source nontrivially depends on the value of $a/\sigma$. We rationalize our findings by using liquid-state theory and the hypernetted chain (HNC) approximation for correlation functions. In addition, we confirm the theoretical trends by studying a three-dimensional glass-former by computer simulations. We show in particular that \rev{a range of $a/\sigma$ below what is commonly considered maximizes the temperature of the critical point, pushing it up in a liquid region where viscosity is low and computer investigations are easier due to a significantly faster equilibration., Comment: 19 pages, 15 figures Accepted for publication in JCP

The swap Monte Carlo algorithm allows the preparation of highly stable glassy configurations for a number of glass-formers, but is inefficient for some models, such as the much studied binary Kob-Andersen (KA) mixture. We have recently developed generalisations to the KA model where swap can be very effective. Here, we show that these models can in turn be used to considerably enhance the stability of glassy configurations in the original KA model at no computational cost. We successfully develop several numerical strategies both in and out of equilibrium to achieve this goal and show how to optimise them. We provide several physical measurements indicating that the proposed algorithms considerably enhance mechanical and thermodynamic stability in the KA model, including a transition towards brittle yielding behaviour. Our results thus pave the way for future studies of stable glasses using the KA model., Comment: 7 pages, 5 figs

We review the field of the glass transition, glassy dynamics and aging from a statistical mechanics perspective. We give a brief introduction to the subject and explain the main phenomenology encountered in glassy systems, with a particular emphasis on spatially heterogeneous dynamics. We review the main theoretical approaches currently available to account for these glassy phenomena, including recent developments regarding mean-field theory of liquids and glasses, novel computational tools, and connections to the jamming transition. Finally, the physics of aging and off-equilibrium dynamics exhibited by glassy materials is discussed., Comment: 50 pages, 24 figs. This is an updated version of a chapter initially written in 2009 for the Encyclopedia of Complexity and Systems Science (Springer). arXiv admin note: text overlap with arXiv:1011.2578

It has recently become possible to prepare ultrastable glassy materials characterised by structural relaxation times which vastly exceed the duration of any feasible experiment. Similarly, new algorithms have led to the production of ultrastable computer glasses. Is it possible to obtain a reliable estimate of a structural relaxation time that is too long to be measured? We review, organise, and critically discuss various methods to estimate very long relaxation times. We also perform computer simulations of three dimensional ultrastable hard spheres glasses to test and quantitatively compare some of these methods for a single model system. The various estimation methods disagree significantly and it is not yet clear how to accurately estimate extremely long relaxation times., Comment: 17 pages, 10 figures; version accepted for publication at J. Chem. Phys

We use computer simulations to investigate the extended phase diagram of a supercooled liquid linearly coupled to a quenched reference configuration. An extensive finite-size scaling analysis demonstrates the existence of a random-field Ising model (RFIM) critical point and of a first-order transition line, in agreement with recent field-theoretical approaches. The dynamics in the vicinity of this critical point resembles the peculiar activated scaling of RFIM-like systems, and the overlap autocorrelation displays a logarithmic stretching. Our study demonstrates RFIM criticality in the thermodynamic limit for a three-dimensional supercooled liquids at equilibrium., Comment: 9 pages, 6 figures To be published in PRE

We develop a generic strategy and simple numerical models for multi-component metallic glasses for which the swap Monte Carlo algorithm can produce highly stable equilibrium configurations equivalent to experimental systems cooled more than $10^7$ times slower than in conventional simulations. This paves the way for a deeper understanding of thermodynamic, dynamic, and mechanical properties of metallic glasses. As first applications, we considerably extend configurational entropy measurements down to the experimental glass temperature, and demonstrate a qualitative change of the mechanical response of metallic glasses of increasing stability towards brittleness., Comment: 11 pages, 9 figs

We develop a minimal model to describe growing dense active matter such as biological tissues, bacterial colonies and biofilms, that are driven by a competition between particle division and steric repulsion. We provide a detailed numerical analysis of collective and single particle dynamics. We show that the microscopic dynamics can be understood as the superposition of an affine radial component due to the global growth, and of a more complex non-affine component which displays features typical of driven soft glassy materials, such as aging, compressed exponential decay of time correlation functions, and a crossover from superdiffusive behaviour at short scales to subdiffusive behaviour at larger scales. This analogy emerges because particle division at the microscale leads to a global expansion which then plays a role analogous to shear flow in soft driven glasses. We conclude that growing dense active matter and sheared dense suspensions can generically be described by the same underlying physics., Comment: 14 pages, 6 figures

Finite-dimensional signatures of spinodal criticality are notoriously difficult to come by. The dynamical transition of glass-forming liquids, first described by mode-coupling theory, is a spinodal instability preempted by thermally activated processes that also limit how close the instability can be approached. We combine numerical tools to directly observe vestiges of the spinodal criticality in finite-dimensional glass formers. We use the swap Monte Carlo algorithm to efficiently thermalise configurations beyond the mode-coupling crossover, and analyze their dynamics using a scheme to screen out activated processes, in spatial dimensions ranging from $d=3$ to $d=9$. We observe a strong softening of the mean-field square-root singularity in $d=3$ that is progressively restored as $d$ increases above $d=8$, in surprisingly good agreement with perturbation theory., Comment: 6 pages, 6 figures

Amorphous solids display a ductile to brittle transition as the kinetic stability of the quiescent glass is increased, which leads to a material failure controlled by the sudden emergence of a macroscopic shear band in quasi-static protocols. We numerically study how finite deformation rates influence ductile and brittle yielding behaviors using model glasses in two and three spatial dimensions. We find that a finite shear rate systematically enhances the stress overshoot of poorly-annealed systems, without necessarily producing shear bands. For well-annealed systems, the non-equilibrium discontinuous yielding transition is smeared out by finite shear rates and it is accompanied by the emergence of multiple shear bands that have been also reported in metallic glass experiments. We show that the typical size of the bands and the distance between them increases algebraically with the inverse shear rate. We provide a dynamic scaling argument for the corresponding lengthscale, based on the competition between the deformation rate and the propagation time of the shear bands., Comment: 12 pages, 8 figures