-
Monte Carlo simulations of glass-forming liquids beyond Metropolis
Authors:
Ludovic Berthier,
Federico Ghimenti Frédéric van Wijland
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 28 June, 2024;
originally announced June 2024.
-
Irreversible Boltzmann samplers in dense liquids: weak-coupling approximation and mode-coupling theory
Authors:
Federico Ghimenti,
Ludovic Berthier,
Grzegorz Szamel,
Frédéric van Wijland
Abstract:
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 quanti…
▽ More
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.
△ Less
Submitted 23 April, 2024;
originally announced April 2024.
-
Normalizing flows as an enhanced sampling method for atomistic supercooled liquids
Authors:
Gerhard Jung,
Giulio Biroli,
Ludovic Berthier
Abstract:
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 d…
▽ More
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.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Transverse forces and glassy liquids in infinite dimensions
Authors:
Federico Ghimenti,
Ludovic Berthier,
Grzegorz Szamel,
Frédéric van Wijland
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 17 June, 2024; v1 submitted 16 February, 2024;
originally announced February 2024.
-
Irreversible Monte Carlo algorithms for hard disk glasses: from event-chain to collective swaps
Authors:
Federico Ghimenti,
Ludovic Berthier,
Frédéric van Wijland
Abstract:
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 tra…
▽ More
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.
△ Less
Submitted 11 July, 2024; v1 submitted 9 February, 2024;
originally announced February 2024.
-
Yielding and plasticity in amorphous solids
Authors:
Ludovic Berthier,
Giulio Biroli,
M. Lisa Manning,
Francesco Zamponi
Abstract:
The physics of disordered media, from metallic glasses to colloidal suspensions, granular matter and biological tissues, offers difficult challenges because it often occurs far from equilibrium, in materials lacking symmetries and evolving through complex energy landscapes. Here, we review recent theoretical efforts to provide microscopic insights into the mechanical properties of amorphous media…
▽ More
The physics of disordered media, from metallic glasses to colloidal suspensions, granular matter and biological tissues, offers difficult challenges because it often occurs far from equilibrium, in materials lacking symmetries and evolving through complex energy landscapes. Here, we review recent theoretical efforts to provide microscopic insights into the mechanical properties of amorphous media using approaches from statistical mechanics as unifying frameworks. We cover both the initial regime corresponding to small deformations, and the yielding transition marking a change between elastic response and plastic flow. We discuss the specific features arising for systems evolving near a jamming transition, and extend our discussion to recent studies of the rheology of dense biological and active materials.
△ Less
Submitted 17 January, 2024;
originally announced January 2024.
-
Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions
Authors:
Thibaut Divoux,
Elisabeth Agoritsas,
Stefano Aime,
Catherine Barentin,
Jean-Louis Barrat,
Roberto Benzi,
Ludovic Berthier,
Dapeng Bi,
Giulio Biroli,
Daniel Bonn,
Philippe Bourrianne,
Mehdi Bouzid,
Emanuela Del Gado,
Hélène Delanoë-Ayari,
Kasra Farain,
Suzanne Fielding,
Matthias Fuchs,
Jasper van der Gucht,
Silke Henkes,
Maziyar Jalaal,
Yogesh M. Joshi,
Anaël Lemaître,
Robert L. Leheny,
Sébastien Manneville,
Kirsten Martens
, et al. (15 additional authors not shown)
Abstract:
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 worksh…
▽ More
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.
△ Less
Submitted 21 December, 2023;
originally announced December 2023.
-
Roadmap on machine learning glassy liquids
Authors:
Gerhard Jung,
Rinske M. Alkemade,
Victor Bapst,
Daniele Coslovich,
Laura Filion,
François P. Landes,
Andrea Liu,
Francesco Saverio Pezzicoli,
Hayato Shiba,
Giovanni Volpe,
Francesco Zamponi,
Ludovic Berthier,
Giulio Biroli
Abstract:
Unraveling the connections between microscopic structure, emergent physical properties, and slow dynamics has long been a challenge in the field of the glass transition. The absence of clear visible structural order in amorphous configurations complicates the identification of the key features related to structural relaxation and transport properties. The difficulty in sampling equilibrated config…
▽ More
Unraveling the connections between microscopic structure, emergent physical properties, and slow dynamics has long been a challenge in the field of the glass transition. The absence of clear visible structural order in amorphous configurations complicates the identification of the key features related to structural relaxation and transport properties. The difficulty in sampling equilibrated configurations at low temperatures hampers thorough numerical and theoretical investigations. This roadmap 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 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 glassy liquids. To foster a collaborative community effort, the article introduces the "GlassBench" dataset, providing simulation data and benchmarks for both two-dimensional and three-dimensional glass-formers. Emphasizing the importance of benchmarks, we identify critical metrics for comparing 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 understanding of glassy liquids.
△ Less
Submitted 23 November, 2023;
originally announced November 2023.
-
Dynamic heterogeneity at the experimental glass transition predicted by transferable machine learning
Authors:
Gerhard Jung,
Giulio Biroli,
Ludovic Berthier
Abstract:
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 transiti…
▽ More
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.
△ Less
Submitted 31 October, 2023;
originally announced October 2023.
-
Direct numerical analysis of dynamic facilitation in glass-forming liquids
Authors:
Cecilia Herrero,
Ludovic Berthier
Abstract:
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 sprea…
▽ More
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.
△ Less
Submitted 24 June, 2024; v1 submitted 25 October, 2023;
originally announced October 2023.
-
Collective relaxation dynamics in a three-dimensional lattice glass model
Authors:
Yoshihiko Nishikawa,
Ludovic Berthier
Abstract:
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 facili…
▽ More
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.
△ Less
Submitted 7 February, 2024; v1 submitted 16 July, 2023;
originally announced July 2023.
-
Sampling efficiency of transverse forces in dense liquids
Authors:
Federico Ghimenti,
Ludovic Berthier,
Grzegorz Szamel,
Frédéric van Wijland
Abstract:
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-…
▽ More
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.
△ Less
Submitted 14 December, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
Emerging mesoscale flows and chaotic advection in dense active matter
Authors:
Yann-Edwin Keta,
Juliane Klamser,
Robert L. Jack,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Front propagation in ultrastable glasses is dynamically heterogeneous
Authors:
Cecilia Herrero,
Mark D. Ediger,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 April, 2023;
originally announced April 2023.
-
Two-dimensional crystals far from equilibrium
Authors:
Leonardo Galliano,
Michael E. Cates,
Ludovic Berthier
Abstract:
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 cr…
▽ More
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.
△ Less
Submitted 24 May, 2023; v1 submitted 22 February, 2023;
originally announced February 2023.
-
Intermittent relaxation and avalanches in extremely persistent active matter
Authors:
Yann-Edwin Keta,
Rituparno Mandal,
Peter Sollich,
Robert L. Jack,
Ludovic Berthier
Abstract:
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 elasti…
▽ More
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.
△ Less
Submitted 5 June, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Finding defects in glasses through machine learning
Authors:
Simone Ciarella,
Dmytro Khomenko,
Ludovic Berthier,
Felix C. Mocanu,
David R. Reichman,
Camille Scalliet,
Francesco Zamponi
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 May, 2023; v1 submitted 11 December, 2022;
originally announced December 2022.
-
Predicting dynamic heterogeneity in glass-forming liquids by physics-inspired machine learning
Authors:
Gerhard Jung,
Giulio Biroli,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 28 September, 2023; v1 submitted 29 October, 2022;
originally announced October 2022.
-
Are supercooled liquids Fickian yet non Gaussian?
Authors:
Ludovic Berthier,
Elijah Flenner,
Grzegorz Szamel
Abstract:
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.…
▽ More
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.
△ Less
Submitted 4 October, 2023; v1 submitted 13 October, 2022;
originally announced October 2022.
-
Two-step devitrification of ultrastable glasses
Authors:
Cecilia Herrero,
Camille Scalliet,
M. D. Ediger,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 January, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Dynamic Gardner crossover in a simple structural glass
Authors:
Qinyi Liao,
Ludovic Berthier,
Hai-Jun Zhou,
Ning Xu
Abstract:
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 dynami…
▽ More
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.
△ Less
Submitted 11 June, 2023; v1 submitted 22 September, 2022;
originally announced September 2022.
-
Microscopic observation of two-level systems in a metallic glass model
Authors:
Felix C. Mocanu,
Ludovic Berthier,
Simone Ciarella,
Dmytro Khomenko,
David R. Reichman,
Camille Scalliet,
Francesco Zamponi
Abstract:
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 configurati…
▽ More
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.
△ Less
Submitted 4 January, 2023; v1 submitted 20 September, 2022;
originally announced September 2022.
-
Modern computational studies of the glass transition
Authors:
Ludovic Berthier,
David R. Reichman
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 January, 2023; v1 submitted 3 August, 2022;
originally announced August 2022.
-
Is glass a state of matter?
Authors:
Benjamin Guiselin,
Gilles Tarjus,
Ludovic Berthier
Abstract:
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 essentiall…
▽ More
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.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Thirty milliseconds in the life of a supercooled liquid
Authors:
Camille Scalliet,
Benjamin Guiselin,
Ludovic Berthier
Abstract:
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 parti…
▽ More
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.
△ Less
Submitted 15 December, 2022; v1 submitted 1 July, 2022;
originally announced July 2022.
-
Computer simulations of the glass transition and glassy materials
Authors:
Jean-Louis Barrat,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 October, 2022; v1 submitted 2 June, 2022;
originally announced June 2022.
-
Collective dynamics in a glass-former with Mari-Kurchan interactions
Authors:
Yoshihiko Nishikawa,
Atsushi Ikeda,
Ludovic Berthier
Abstract:
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. Ne…
▽ More
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.
△ Less
Submitted 23 June, 2022; v1 submitted 11 April, 2022;
originally announced April 2022.
-
Violation of the fluctuation-dissipation theorem and effective temperatures in spin ice
Authors:
Valentin Raban,
Ludovic Berthier,
Peter C. W. Holdsworth
Abstract:
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, to…
▽ More
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.
△ Less
Submitted 12 April, 2022; v1 submitted 14 February, 2022;
originally announced February 2022.
-
Static self-induced heterogeneity in glass-forming liquids: Overlap as a microscope
Authors:
Benjamin Guiselin,
Gilles Tarjus,
Ludovic Berthier
Abstract:
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 overl…
▽ More
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.
△ Less
Submitted 6 May, 2022; v1 submitted 25 January, 2022;
originally announced January 2022.
-
Disordered collective motion in dense assemblies of persistent particles
Authors:
Yann-Edwin Keta,
Robert L. Jack,
Ludovic Berthier
Abstract:
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, char…
▽ More
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
△ Less
Submitted 23 June, 2022; v1 submitted 13 January, 2022;
originally announced January 2022.
-
Excess wings and asymmetric relaxation spectra in a facilitated trap model
Authors:
Camille Scalliet,
Benjamin Guiselin,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 August, 2021; v1 submitted 3 June, 2021;
originally announced June 2021.
-
Relaxation dynamics in the energy landscape of glass-forming liquids
Authors:
Yoshihiko Nishikawa,
Misaki Ozawa,
Atsushi Ikeda,
Pinaki Chaudhuri,
Ludovic Berthier
Abstract:
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 findin…
▽ More
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.
△ Less
Submitted 6 April, 2022; v1 submitted 3 June, 2021;
originally announced June 2021.
-
Statistical mechanics of coupled supercooled liquids in finite dimensions
Authors:
Benjamin Guiselin,
Ludovic Berthier,
Gilles Tarjus
Abstract:
We study the statistical mechanics of supercooled liquids when the system evolves at a temperature $T$ with a field $ε$ 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…
▽ More
We study the statistical mechanics of supercooled liquids when the system evolves at a temperature $T$ with a field $ε$ 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 $(ε,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 $ε$ 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.
△ Less
Submitted 26 January, 2022; v1 submitted 19 May, 2021;
originally announced May 2021.
-
Microscopic origin of excess wings in relaxation spectra of supercooled liquids
Authors:
Benjamin Guiselin,
Camille Scalliet,
Ludovic Berthier
Abstract:
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 computati…
▽ More
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.
△ Less
Submitted 6 May, 2022; v1 submitted 2 March, 2021;
originally announced March 2021.
-
Rare events and disorder control the brittle yielding of well-annealed amorphous solids
Authors:
Misaki Ozawa,
Ludovic Berthier,
Giulio Biroli,
Gilles Tarjus
Abstract:
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 fa…
▽ More
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.
△ Less
Submitted 11 July, 2022; v1 submitted 10 February, 2021;
originally announced February 2021.
-
Self-induced heterogeneity in deeply supercooled liquids
Authors:
Ludovic Berthier
Abstract:
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 dy…
▽ More
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.
△ Less
Submitted 23 August, 2021; v1 submitted 23 October, 2020;
originally announced October 2020.
-
Glassy behaviour of sticky spheres: What lies beyond experimental timescales?
Authors:
Christopher J. Fullerton,
Ludovic Berthier
Abstract:
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-separa…
▽ More
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.
△ Less
Submitted 18 December, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Relaxation dynamics of non-Brownian spheres below jamming
Authors:
Yoshihiko Nishikawa,
Atsushi Ikeda,
Ludovic Berthier
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 9 February, 2021; v1 submitted 18 July, 2020;
originally announced July 2020.
-
On the overlap between configurations in glassy liquids
Authors:
Benjamin Guiselin,
Gilles Tarjus,
Ludovic Berthier
Abstract:
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/σ$ of the inter-particle distance, associated with how precisely similar two configurations must be for be…
▽ More
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/σ$ 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/σ$, 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/σ$. 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/σ$ 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.
△ Less
Submitted 26 November, 2020; v1 submitted 15 July, 2020;
originally announced July 2020.
-
Stable glassy configurations of the Kob-Andersen model using swap Monte Carlo
Authors:
Anshul D. S. Parmar,
Benjamin Guiselin,
Ludovic Berthier
Abstract:
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 stabilit…
▽ More
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.
△ Less
Submitted 6 October, 2020; v1 submitted 18 June, 2020;
originally announced June 2020.
-
Glasses and aging: A Statistical Mechanics Perspective
Authors:
Francesco Arceri,
François P. Landes,
Ludovic Berthier,
Giulio Biroli
Abstract:
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 re…
▽ More
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.
△ Less
Submitted 8 October, 2020; v1 submitted 17 June, 2020;
originally announced June 2020.
-
How to "measure" a structural relaxation time that is too long to be measured?
Authors:
Ludovic Berthier,
Mark D. Ediger
Abstract:
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 c…
▽ More
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.
△ Less
Submitted 22 July, 2020; v1 submitted 13 May, 2020;
originally announced May 2020.
-
Random-field Ising model criticality in a glass-forming liquid
Authors:
Benjamin Guiselin,
Ludovic Berthier,
Gilles Tarjus
Abstract:
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…
▽ More
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.
△ Less
Submitted 23 October, 2020; v1 submitted 22 April, 2020;
originally announced April 2020.
-
Predicting plasticity in disordered solids from structural indicators
Authors:
D. Richard,
M. Ozawa,
S. Patinet,
E. Stanifer,
B. Shang,
S. A. Ridout,
B. Xu,
G. Zhang,
P. K. Morse,
J. -L. Barrat,
L. Berthier,
M. L. Falk,
P. Guan,
A. J. Liu,
K. Martens,
S. Sastry,
D. Vandembroucq,
E. Lerner,
M. L. Manning
Abstract:
Amorphous solids lack long-range order. Therefore identifying structural defects -- akin to dislocations in crystalline solids -- that carry plastic flow in these systems remains a daunting challenge. By comparing many different structural indicators in computational models of glasses, under a variety of conditions we carefully assess which of these indicators are able to robustly identify the str…
▽ More
Amorphous solids lack long-range order. Therefore identifying structural defects -- akin to dislocations in crystalline solids -- that carry plastic flow in these systems remains a daunting challenge. By comparing many different structural indicators in computational models of glasses, under a variety of conditions we carefully assess which of these indicators are able to robustly identify the structural defects responsible for plastic flow in amorphous solids. We further demonstrate that the density of defects changes as a function of material preparation and strain in a manner that is highly correlated with the macroscopic material response. Our work represents an important step towards predicting how and when an amorphous solid will fail from its microscopic structure.
△ Less
Submitted 23 November, 2020; v1 submitted 25 March, 2020;
originally announced March 2020.
-
Ultrastable metallic glasses in silico
Authors:
Anshul D. S. Parmar,
Misaki Ozawa,
Ludovic Berthier
Abstract:
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 met…
▽ More
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.
△ Less
Submitted 8 October, 2020; v1 submitted 4 February, 2020;
originally announced February 2020.
-
Analogies between growing dense active matter and soft driven glasses
Authors:
Elsen Tjhung,
Ludovic Berthier
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 December, 2020; v1 submitted 3 February, 2020;
originally announced February 2020.
-
Finite-dimensional vestige of spinodal criticality above the dynamical glass transition
Authors:
Ludovic Berthier,
Patrick Charbonneau,
Joyjit Kundu
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 December, 2019;
originally announced December 2019.
-
Brittle yielding of amorphous solids at finite shear rates
Authors:
Murari Singh,
Misaki Ozawa,
Ludovic Berthier
Abstract:
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. W…
▽ More
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.
△ Less
Submitted 13 December, 2019;
originally announced December 2019.
-
Role of fluctuations in the yielding transition of two-dimensional glasses
Authors:
Misaki Ozawa,
Ludovic Berthier,
Giulio Biroli,
Gilles Tarjus
Abstract:
We numerically study yielding in two-dimensional glasses which are generated with a very wide range of stabilities by swap Monte-Carlo simulations and then slowly deformed at zero temperature. We provide strong numerical evidence that stable glasses yield via a nonequilibrium discontinuous transition in the thermodynamic limit. A critical point separates this brittle yielding from the ductile one…
▽ More
We numerically study yielding in two-dimensional glasses which are generated with a very wide range of stabilities by swap Monte-Carlo simulations and then slowly deformed at zero temperature. We provide strong numerical evidence that stable glasses yield via a nonequilibrium discontinuous transition in the thermodynamic limit. A critical point separates this brittle yielding from the ductile one observed in less stable glasses. We find that two-dimensional glasses yield similarly to their three-dimensional counterparts but display larger sample-to-sample disorder-induced fluctuations, stronger finite-size effects, and rougher spatial wandering of the observed shear bands. These findings strongly constrain effective theories of yielding.
△ Less
Submitted 22 May, 2020; v1 submitted 12 December, 2019;
originally announced December 2019.
-
Glass stability changes the nature of yielding under oscillatory shear
Authors:
Wei-Ting Yeh,
Misaki Ozawa,
Kunimasa Miyazaki,
Takeshi Kawasaki,
Ludovic Berthier
Abstract:
We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our resul…
▽ More
We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our results disentangle the relative role of mechanical and thermal annealing on the mechanical properties of amorphous solids, which is relevant for various experimental situations from the rheology of soft materials to fatigue failure in metallic glasses.
△ Less
Submitted 6 June, 2020; v1 submitted 29 November, 2019;
originally announced November 2019.