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

The role of excitations in supercooled liquids: Density, geometry, and relaxation dynamics

Low-energy excitations play a key role in all condensed-matter systems, yet there is limited understanding of their nature in glasses, where they correspond to local rearrangements of groups of particles. Here, we introduce an algorithm to systematically uncover these excitations up to the activation energy scale relevant to structural relaxation. We use it in a model system to measure the density of states on a scale never achieved before, confirming that this quantity shifts to higher energy under cooling, precisely as the activation energy does. Second, we show that the excitations' energetic and spatial features allow one to predict with great accuracy the dynamic propensity, i.e., the location of future relaxation dynamics. Finally, we find that excitations have a primary field whose properties, including the displacement of the most mobile particle, scale as a power-law of their activation energy and are independent of temperature. Additionally, they exhibit an outer deformation field that depends on the material's stability and, therefore, on temperature. We build a scaling description of these findings. Overall, our analysis supports that excitations play a crucial role in regulating relaxation dynamics near the glass transition, effectively suppressing the transition to dynamical arrest predicted by mean-field theories while also being strongly influenced by it.

Characterizing the Slow Dynamics of the Swap Monte Carlo Algorithm

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

Local vs. cooperative: Unraveling glass transition mechanisms with SEER

Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

The nature of non-phononic excitations in disordered systems

The frequency scaling exponent of low-frequency excitations in microscopically small glasses, which do not allow for the existence of waves (phonons), has been in the focus of the recent literature. The density of states g(ω) of these modes obeys an ωs scaling, where the exponent s, ranging between 2 and 5, depends on the quenching protocol. The orgin of these findings remains controversal. Here we show, using heterogeneous-elasticity theory, that in a marginally-stable glass sample g(ω) follows a Debye-like scaling (s = 2), and the associated excitations (type-I) are of random-matrix type. Further, using a generalisation of the theory, we demonstrate that in more stable samples, other, (type-II) excitations prevail, which are non-irrotational oscillations, associated with local frozen-in stresses. The corresponding frequency scaling exponent s is governed by the statistics of small values of the stresses and, therefore, depends on the details of the interaction potential.

Low-frequency vibrational density of states of ordinary and ultra-stable glasses

A remarkable feature of disordered solids distinct from crystals is the violation of the Debye scaling law of the low-frequency vibrational density of states. Because the low-frequency vibration is responsible for many properties of solids, it is crucial to elucidate it for disordered solids. Numerous recent studies have suggested power-law scalings of the low-frequency vibrational density of states, but the scaling exponent is currently under intensive debate. Here, by classifying disordered solids into stable and unstable ones, we find two distinct and robust scaling exponents for non-phononic modes at low frequencies. Using the competition of these two scalings, we clarify the variation of the scaling exponent and hence reconcile the debate. Via the study of both ordinary and ultra-stable glasses, our work reveals a comprehensive picture of the low-frequency vibration of disordered solids and sheds light on the low-frequency vibrational features of ultra-stable glasses on approaching the ideal glass.

Structure of quantum supercooled liquids

Supercooled liquids show a drastic slowdown in the dynamics with decreasing temperature, while their structure remains similar to that of normal liquids. In this paper, the structural features in a quantum supercooled liquid are explored in terms of cages defined using the Voronoi polyhedra and characterized in terms of their volumes and geometries. The cage volume fluctuations are sensitive to the quantum effects, and decrease as the glass transition is approached by varying the quantumness. This is in contrast to the classical case where the volumes are insensitive to temperature variations as one approaches the transition. The cage geometry becomes more spherical upon increasing quantumness from zero, pushing the system closer to the glass transition. The cage geometry is found to be significantly correlated with asymmetry in the position uncertainty of the caged particle in the strongly quantum regime.

Depletion of two-level systems in highly stable glasses with different molecular ordering

Recent findings of structural glasses with extremely high kinetic and thermodynamic stability have attracted much attention. The question has been raised as to whether the well-known, low-temperature “glassy anomalies” (attributed to the presence of two-level systems [TLS] and the “boson peak”) persist or not in these ultrastable glasses of much lower configurational entropy. To resolve previous contradictory results, we study a particular type of ultrastable glass, TPD, which can be prepared by physical vapor deposition in a highly-stable state with different degrees of layering and molecular orientation, and also as a conventional glass and in crystalline state. After a thorough characterization of the different samples prepared, we have measured their specific heat down to 0.4 K. Whereas the conventional glass exhibits the typical glassy behaviour and the crystal the expected Debye cubic dependence at very low temperatures, a strong depletion of the TLS contribution is found in both kinds of ultrastable glass, regardless of their layering and molecular ordering.

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.

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.

Quantum uncertainty effects in the dynamics of supercooled liquids: A molecular dynamics study

Dynamics of density fluctuations in quantum supercooled liquids is analyzed using molecular dynamics simulations. In contrast to the classical case, the uncertainty in the particle position (delocalization of quantum particle in space) leads to significant differences in the dynamics of quantum liquids, both in the short- and long-time limits. The effect of uncertainty is found to be significant for length scales smaller than the uncertainty itself, and diminishes as the length scale grows. The dynamic heterogeneity of the system at short times is enhanced due to uncertainty. In the intermediate (β-relaxation) time regime, the heterogeneity tends to get suppressed due to quantum uncertainty. The probability distribution of particle displacements shows highly nonclassical behavior with double-peak structure at short timescales.

Low-frequency vibrational states in ideal glasses with random pinning

Glasses exhibit spatially localized vibrations in the low-frequency regime. These localized modes emerge below the boson peak frequency ω_{BP}, and their vibrational densities of state follow g(ω)∝ω^{4} (ω is frequency). Here, we attempt to address how the localized vibrations behave through the ideal glass transition. To do this, we employ a random pinning method, which enables us to study the thermodynamic glass transition. We find that the localized vibrations survive even in equilibrium glass states. Remarkably, the localized vibrations still maintain the properties of appearance below ω_{BP} and g(ω)∝ω^{4}. Our results provide important insight into the material properties of ideal glasses.

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

We numerically study the equilibrium relaxation dynamics of a two-dimensional Mari-Kurchan glass model. The tree-like structure of particle interactions forbids both 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.

Nuclear quantum effects on the dynamics and glass behavior of a monatomic liquid with two liquid states

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid-liquid phase transition and liquid-liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., Tg LDL(P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., Tg HDL(P) increases upon compression. (iii) NQE shift both Tg LDL(P) and Tg HDL(P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA-HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA-HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA-HDA transformation.

Ultrastable glasses: new perspectives for an old problem

Ultrastable glasses (mostly prepared from the vapor phase under optimized deposition conditions) represent a unique class of materials with low enthalpies and high kinetic stabilities. These highly stable and dense glasses show unique physicochemical properties, such as high thermal stability, improved mechanical properties or anomalous transitions into the supercooled liquid, offering unprecedented opportunities to understand many aspects of the glassy state. Their improved properties with respect to liquid-cooled glasses also open new prospects to their use in applications where liquid-cooled glasses failed or where not considered as usable materials. In this review article we summarize the state of the art of vapor-deposited (and other) ultrastable glasses with a focus on the mechanism of equilibration, the transformation to the liquid state and the low temperature properties. The review contains information on organic, metallic, polymeric and chalcogenide glasses and an updated list with relevant properties of all materials known today to form a stable glass.

Mean-field description for the architecture of low-energy excitations in glasses

In amorphous materials, groups of particles can rearrange locally into a new stable configuration. Such elementary excitations are key as they determine the response to external stresses, as well as to thermal and quantum fluctuations. Yet, understanding what controls their geometry remains a challenge. Here we build a scaling description of the geometry and energy of low-energy excitations in terms of the distance to an instability, as predicted, for instance, at the dynamical transition in mean-field approaches of supercooled liquids. We successfully test our predictions in ultrastable computer glasses, with a gapped spectrum and an ungapped (regular) spectrum. Overall, our approach explains why excitations become less extended, with a higher energy and displacement scale upon cooling.

Emergence of two-level systems in glass formers: a kinetic Monte Carlo study

Using a distinguishable-particle lattice model based on void-induced dynamics, we successfully reproduce the well-known linear relation between heat capacity and temperature at very low temperatures. The heat capacity is dominated by two-level systems formed due to the strong localization of voids to two neighboring sites, and can be exactly calculated in the limit of ultrastable glasses. Similar but weaker localization at higher temperatures accounts for glass transition. The result supports the conventional two-level tunneling picture by revealing how two-level systems emerge from random particle interactions, which also cause glass transition. Our approach provides a unified framework for relating microscopic dynamics of glasses at room and cryogenic temperatures.

Low-energy quasilocalized excitations in structural glasses

Glassy solids exhibit a wide variety of generic thermomechanical properties, ranging from universal anomalous specific heat at cryogenic temperatures to nonlinear plastic yielding and failure under external driving forces, which qualitatively differ from their crystalline counterparts. For a long time, it has been believed that many of these properties are intimately related to nonphononic, low-energy quasilocalized excitations (QLEs) in glasses. Indeed, recent computer simulations have conclusively revealed that the self-organization of glasses during vitrification upon cooling from a melt leads to the emergence of such QLEs. In this Perspective, we review developments over the past three decades toward understanding the emergence of QLEs in structural glasses and the degree of universality in their statistical and structural properties. We discuss the challenges and difficulties that hindered progress in achieving these goals and review the frameworks put forward to overcome them. We conclude with an outlook on future research directions and open questions.

Emergent Fractal Energy Landscape as the Origin of Stress-Accelerated Dynamics in Amorphous Solids

The ageing dynamics in a multiplicity of metastable glasses are investigated at various thermomechanical conditions. By using data analytics to deconvolute the integral effects of environmental factors (e.g., energy level, temperature, stress), and by directly scrutinizing the minimum energy pathways for local excitations, we demonstrate external shear would make the system's energy landscape surprisingly fractal and create an emergent low-barrier mode with highly tortuous pathways, leading to an accelerated relaxation. This finding marks a departure from the classic picture of shear-induced simple bias of energy landscape. The insights and implications of this study are also discussed.

Microscopic Structural Evolution during Ultrastable Metallic Glass Formation

By decreasing the rate of physical vapor deposition, ZrCuAl metallic glasses with improved stability and mechanical performances can be formed, while the microscopic structural mechanisms remain unclear. Here, with scanning transmission electron microscopy and high-energy synchrotron X-ray diffraction, we found that the metallic glass deposited at a higher rate exhibits a heterogeneous structure with compositional fluctuations at a distance of a few nanometers, which gradually disappear on decreasing the deposition rate; eventually, a homogeneous structure is developed approaching ultrastability. This microscopic structural evolution suggests the existence of the following two dynamical processes during ultrastable metallic glass formation: a faster diffusion process driven by the kinetic energy of the depositing atoms, which results in nanoscale compositional fluctuations, and a slower collective relaxation process that eliminates the compositional and structural heterogeneity, equilibrates the deposited atoms, and strengthens the local atomic connectivity.

Glass Dynamics Deep in the Energy Landscape

When a liquid is cooled, progress down the energy landscape is arrested near the glass transition temperature T_g. In principle, lower energy states can be accessed by waiting for further equilibration, but the rough energy landscape of glasses quickly leads to kinetics on geologically slow time scales below T_g. Over the past decade, progress has been made probing deeper into the energy landscape via several techniques. By looking at bulk and surface diffusion, using layered deposition that promotes equilibration, imaging glass surfaces with faster dynamics below T_g, and optically exciting glasses, experiments have moved into a regime of ultrastable, low energy glasses that was difficult to access in the past. At the same time, both simulations and energy landscape theory based on a random first order transition (RFOT) have tackled systems that include surfaces, optical excitation, and interfacial dynamics. Here we review some of the recent experimental work, and how energy landscape theory illuminates glassy dynamics well below the glass transition temperature by making direct connections between configurational entropy, energy landscape barriers, and the resulting dynamics.

Thermal origin of quasilocalized excitations in glasses

Key aspects of glasses are controlled by the presence of excitations in which a group of particles can rearrange. Surprisingly, recent observations indicate that their density is dramatically reduced and their size decreases as the temperature of the supercooled liquid is lowered. Some theories predict these excitations to cause a gap in the spectrum of quasilocalized modes of the Hessian that grows upon cooling, while others predict a pseudogap D_{L}(ω)∼ω^{α}. To unify these views and observations, we generate glassy configurations of controlled gap magnitude ω_{c} at temperature T=0, using so-called breathing particles, and study how such gapped states respond to thermal fluctuations. We find that (i) the gap always fills up at finite T with D_{L}(ω)≈A_{4}(T)ω^{4} and A_{4}∼exp(-E_{a}/T) at low T, (ii) E_{a} rapidly grows with ω_{c}, in reasonable agreement with a simple scaling prediction E_{a}∼ω_{c}^{4} and (iii) at larger ω_{c} excitations involve fewer particles, as we rationalize, and eventually become stringlike. We propose an interpretation of mean-field theories of the glass transition, in which the modes beyond the gap act as an excitation reservoir, from which a pseudogap distribution is populated with its magnitude rapidly decreasing at lower T. We discuss how this picture unifies the rarefaction as well as the decreasing size of excitations upon cooling, together with a stringlike relaxation occurring near the glass transition.

Chemical control of peptide material phase transitions

Progressive solute-rich polymer phase transitions provide pathways for achieving ordered supramolecular assemblies. Intrinsically disordered protein domains specifically regulate information in biological networks via conformational ordering. Here we consider a molecular tagging strategy to control ordering transitions in polymeric materials and provide a proof-of-principle minimal peptide phase network captured with a dynamic chemical network.

Substrate initiated assembly of a dynamic chemical network.

Intermittent rearrangements accompanying thermal fluctuations distinguish glasses from crystals

Glasses exhibit vibrational and thermal properties that are markedly different from those of crystals. While recent works have advanced our understanding of vibrational excitations in glasses in the harmonic approximation limit, efforts in understanding finite-temperature anharmonic processes have been limited. In crystals, phonon-phonon coupling provides an extremely efficient mechanism for anharmonic decay that is also important in glasses. By using extensive molecular dynamics simulation of model atomic systems, here we first describe, both numerically and analytically, the anharmonic couplings in the crystal and the glass by focusing on the temperature dependence of the associated decay rates. Next, we show that an additional anharmonic channel of different origin emerges in the amorphous case, which induces unconventional intermittent rearrangements of particles. We have found that thermal vibrations in glasses trigger transitions among numerous different local minima of the energy landscape, which, however, are located within the same wide (meta)basin. These processes generate motions that are different from both diffusive and out-of-equilibrium aging dynamics. We suggest that (i) the observed intermittent rearrangements accompanying thermal fluctuations are crucial features distinguishing glasses from crystals and (ii) they can be considered as relics of the liquid state that survive the complete dynamic arrest taking place at the glass transition temperature.

Stable glassy configurations of the Kob-Andersen model using swap Monte Carlo

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 generalizations 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 optimize them. We provide several physical measurements indicating that the proposed algorithms considerably enhance mechanical and thermodynamic stability in the KA model, including a transition toward brittle yielding behavior. Our results thus pave the way for future studies of stable glasses using the KA model.

Pinching a glass reveals key properties of its soft spots

It is now well established that glasses feature quasilocalized nonphononic excitations-coined "soft spots"-, which follow a universal [Formula: see text] density of states in the limit of low frequencies ω. All glass-specific properties, such as the dependence on the preparation protocol or composition, are encapsulated in the nonuniversal prefactor of the universal [Formula: see text] law. The prefactor, however, is a composite quantity that incorporates information both about the number of quasilocalized nonphononic excitations and their characteristic stiffness, in an apparently inseparable manner. We show that by pinching a glass-i.e., by probing its response to force dipoles-one can disentangle and independently extract these two fundamental pieces of physical information. This analysis reveals that the number of quasilocalized nonphononic excitations follows a Boltzmann-like law in terms of the parent temperature from which the glass is quenched. The latter, sometimes termed the fictive (or effective) temperature, plays important roles in nonequilibrium thermodynamic approaches to the relaxation, flow, and deformation of glasses. The analysis also shows that the characteristic stiffness of quasilocalized nonphononic excitations can be related to their characteristic size, a long sought-for length scale. These results show that important physical information, which is relevant for various key questions in glass physics, can be obtained through pinching a glass.