Finding defects in glasses through machine learning.

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

Two-step devitrification of ultrastable glasses.

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

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

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

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

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

Depletion of Two-Level Systems in Ultrastable Computer-Generated Glasses.

Amorphous solids exhibit quasiuniversal low temperature anomalies whose origin has been ascribed to localized tunneling defects. Using an advanced Monte Carlo procedure, we create in silico glasses spanning from hyperquenched to ultrastable glasses. Using a multidimensional path-finding protocol, we locate tunneling defects with energy splittings smaller than k_{B}T_{Q}, with T_{Q} the temperature below which quantum effects are relevant (T_{Q}≈1 K in most experiments). We find that as the stability of a glass increases, its energy landscape as well as the manner in which it is probed tend to deplete the density of tunneling defects, as observed in recent experiments. We explore the real-space nature of tunneling defects, and find that they are mostly localized to a few atoms, but are occasionally dramatically delocalized.

Nature of excitations and defects in structural glasses.

The nature of defects in amorphous materials, analogous to vacancies and dislocations in crystals, remains elusive. Here, we explore their nature in a three-dimensional microscopic model glass-former that describes granular, colloidal, atomic and molecular glasses by changing the temperature and density. We find that all glasses evolve in a very rough energy landscape, with a hierarchy of barrier sizes corresponding to both localized and delocalized excitations. Collective excitations dominate in the jamming regime relevant for granular and colloidal glasses. By moving gradually to larger densities describing atomic and molecular glasses, the system crosses over to a regime dominated by localized defects and relatively simpler landscapes. We quantify the energy and temperature scales associated to these defects and their evolution with density. Our results pave the way to a systematic study of low-temperature physics in a broad range of physical conditions and glassy materials.

Does the Adam-Gibbs relation hold in simulated supercooled liquids?

We perform stringent tests of thermodynamic theories of the glass transition over the experimentally relevant temperature regime for several simulated glass-formers. The swap Monte Carlo algorithm is used to estimate the configurational entropy and static point-to-set lengthscale, and careful extrapolations are used for the relaxation times. We first quantify the relation between configurational entropy and the point-to-set lengthscale in two and three dimensions. We then show that the Adam-Gibbs relation is generally violated in simulated models for the experimentally relevant time window. Collecting experimental data for several supercooled molecular liquids, we show that the same trends are observed experimentally. Deviations from the Adam-Gibbs relation remain compatible with random first order transition theory and may account for the reported discrepancies between Kauzmann and Vogel-Fulcher-Tammann temperatures. Alternatively, they may also indicate that even near Tg thermodynamics is not the only driving force for slow dynamics.

Rejuvenation and Memory Effects in a Structural Glass.

We show numerically that a three-dimensional model for structural glass displays aging, rejuvenation, and memory effects when subjected to a temperature cycle. These effects indicate that the free energy landscape of structural glasses may possess the complex hierarchical structure that characterizes materials such as spin and polymer glasses. We use the theoretical concept of marginal stability to interpret our results, and explain in which physical conditions a complex aging dynamics can emerge in dense supercooled liquids, paving the way for future experimental studies of complex aging dynamics in colloidal and granular glasses.

Configurational entropy of glass-forming liquids.

The configurational entropy is one of the most important thermodynamic quantities characterizing supercooled liquids approaching the glass transition. Despite decades of experimental, theoretical, and computational investigation, a widely accepted definition of the configurational entropy is missing, its quantitative characterization remains fraught with difficulties, misconceptions, and paradoxes, and its physical relevance is vividly debated. Motivated by recent computational progress, we offer a pedagogical perspective on the configurational entropy in glass-forming liquids. We first explain why the configurational entropy has become a key quantity to describe glassy materials, from early empirical observations to modern theoretical treatments. We explain why practical measurements necessarily require approximations that make its physical interpretation delicate. We then demonstrate that computer simulations have become an invaluable tool to obtain precise, nonambiguous, and experimentally relevant measurements of the configurational entropy. We describe a panel of available computational tools, offering for each method a critical discussion. This perspective should be useful to both experimentalists and theoreticians interested in glassy materials and complex systems.

Marginally stable phases in mean-field structural glasses.

A novel form of amorphous matter characterized by marginal stability was recently discovered in the mean-field theory of structural glasses. Using this approach, we provide complete phase diagrams delimiting the location of the marginally stable glass phase for a large variety of pair interactions and physical conditions, extensively exploring physical regimes relevant to granular matter, foams, emulsions, hard and soft colloids, and molecular glasses. We find that all types of glasses may become marginally stable, but the extent of the marginally stable phase highly depends on the preparation protocol. Our results suggest that marginal phases should be observable for colloidal and non-Brownian particles near jamming and for poorly annealed glasses. For well-annealed glasses, two distinct marginal phases are predicted. Our study unifies previous results on marginal stability in mean-field models and will be useful to guide numerical simulations and experiments aimed at detecting marginal stability in finite-dimensional amorphous materials.

Absence of Marginal Stability in a Structural Glass.

Marginally stable solids have peculiar physical properties that were first analyzed in the context of the jamming transition. We theoretically investigate the existence of marginal stability in a prototypical model for structural glass formers, combining analytical calculations in infinite dimensions to computer simulations in three dimensions. While mean-field theory predicts the existence of a Gardner phase transition towards a marginally stable glass phase at low temperatures, simulations show no hint of diverging time scales or length scales, but reveal instead the presence of sparse localized defects. Our results suggest that the Gardner transition is deeply affected by finite dimensional fluctuations, and raise issues about the relevance of marginal stability in structural glasses far away from jamming.

Cages and anomalous diffusion in vibrated dense granular media.

A vertically shaken granular medium hosts a blade rotating around a fixed vertical axis, which acts as a mesorheological probe. At high densities, independently of the shaking intensity, the blade's dynamics shows strong caging effects, marked by transient subdiffusion and a maximum in the velocity power density spectrum, at a resonant frequency ~10 Hz. Interpreting the data through a diffusing harmonic cage model allows us to retrieve the elastic constant of the granular medium and its collective diffusion coefficient. For high frequencies f, a tail ~1/f in the velocity power density spectrum reveals nontrivial correlations in the intracage microdynamics. At very long times (larger than 10 s), a superdiffusive behavior emerges, ballistic in the most extreme cases. Consistently, the distribution of slow velocity inversion times τ displays a power-law decay, likely due to persistent collective fluctuations of the host medium.

Measurements of the dielectric and viscoelastic constants in mixtures of 4,4'-n-octyl-cyanobiphenyl and biphenyl.

We performed measurements of the dielectric constants, splay elastic constant, and rotational viscosity of the nematic phase of mixtures of 4,4'-n-octyl-cyanobiphenyl (8CB) and biphenyl (BP). In contrast with previous results of DasGupta et al. [Phys. Rev. E 63, 041703 (2001); Phys. Lett. A 288, 323 (2001)], we do not find any anomaly of these constants when the smectic-A phase is approached at all concentrations of BP. These results are compatible with recent calorimetric measurements of Denolf et al. [Phys. Rev. Lett. 97, 107801 (2006); Phys. Rev. E 76, 051702 (2007)] and the absence of a tricritical point in the phase diagram. The origin of the anomalies observed by DasGupta et al. at large concentration of BP is also briefly discussed and a likely explanation in terms of biphenyl evaporation is proposed.