Theory for the density of interacting quasilocalized modes in amorphous solids

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.

Enhanced vibrational stability in glass droplets

We show through simulations of amorphous solids prepared in open-boundary conditions that they possess significantly fewer low-frequency vibrational modes compared to their periodic boundary counterparts. Specifically, using measurements of the vibrational density of states, we find that the D ( ω ) ∼ ω 4 law changes to D ( ω ) ∼ ω δ with δ ≈ 5 in two dimensions and δ ≈ 4.5 in three dimensions. Crucially, this enhanced stability is achieved when utilizing slow annealing protocols to generate solid configurations. We perform an anharmonic analysis of the minima corresponding to the lowest frequency modes in such open-boundary systems and discuss their correlation with the density of states. A study of various system sizes further reveals that small systems display a higher degree of localization in vibrations. Lastly, we confine open-boundary solids in order to introduce macroscopic stresses in the system, which are absent in the unconfined system and find that the D ( ω ) ∼ ω 4 behavior is recovered.

Density of states below the first sound mode in 3D glasses

Glasses feature universally low-frequency excess vibrational modes beyond Debye prediction, which could help rationalize, e.g., the glasses’ unusual temperature dependence of thermal properties compared to crystalline solids. The way the density of states of these low-frequency excess modes D( ω) depends on the frequency ω has been debated for decades. Recent simulation studies of 3D glasses suggest that D( ω) scales universally with ω4 in a low-frequency regime below the first sound mode. However, no simulation study has ever probed as low frequencies as possible to test directly whether this quartic law could work all the way to extremely low frequencies. Here, we calculated D( ω) below the first sound mode in 3D glasses over a wide range of frequencies. We find D( ω) scales with ω β with β < 4 at very low frequencies examined, while the ω4 law works only in a limited intermediate-frequency regime in some glasses. Moreover, our further analysis suggests our observation does not depend on glass models or glass stabilities examined. The ω4 law of D( ω) below the first sound mode is dominant in current simulation studies of 3D glasses, and our direct observation of the breakdown of the quartic law at very low frequencies thus leaves an open but important question that may attract more future numerical and theoretical studies.

Low-Frequency Excess Vibrational Modes in Two-Dimensional Glasses

Glasses possess more low-frequency vibrational modes than predicted by Debye theory. These excess modes are crucial for the understanding of the low temperature thermal and mechanical properties of glasses, which differ from those of crystalline solids. Recent simulational studies suggest that the density of the excess modes scales with their frequency ω as ω^{4} in two and higher dimensions. Here, we present extensive numerical studies of two-dimensional model glass formers over a large range of glass stabilities. We find that the density of the excess modes follows D_{exc}(ω)∼ω^{2} up to around the boson peak, regardless of the glass stability. The stability dependence of the overall scale of D_{exc}(ω) correlates with the stability dependence of low-frequency sound attenuation. However, we also find that, in small systems, where the first sound mode is pushed to higher frequencies, at frequencies below the first sound mode, there are excess modes with a system size independent density of states that scales as ω^{3}.

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.

Ioffe-Regel criterion and viscoelastic properties of amorphous solids

We show that viscoelastic effects play a crucial role in the damping of vibrational modes in harmonic amorphous solids. The relaxation of a given plane elastic wave is described by a memory function of a semi-infinite one-dimensional mass-spring chain. The initial vibrational energy spreads from the first site of the chain to infinity. In the beginning of the chain, there is a barrier, which significantly reduces the decay of vibrational energy below the Ioffe-Regel frequency. To obtain the parameters of the chain, we present a numerically stable method, based on the Chebyshev expansion of the local vibrational density of states.

Novel elastic instability of amorphous solids in finite spatial dimensions

Recently, progress has been made in the understanding of anomalous vibrational excitations in amorphous solids. In the lowest-frequency region, the vibrational spectrum follows a non-Debye quartic law, which persists up to zero frequency without any frequency gap. This gapless vibrational density of states (vDOS) suggests that glasses are on the verge of instability. This feature of marginal stability is now highlighted as a key concept in the theories of glasses. In particular, the elasticity theory based on marginal stability predicts the gapless vDOS. However, this theory yields a quadratic law and not the quartic law. To address this inconsistency, we presented a new type of instability, which is different from the conventional one, and proposed that amorphous solids are marginally stable considering the new instability in the preceding study [M. Shimada, H. Mizuno and A. Ikeda, Soft Matter, 2020, 16, 7279]. In this study, we further extend and detail the results for these instabilities. By analyzing various examples of disorder, we demonstrate that real glasses in finite spatial dimensions can be marginally stable by the proposed novel instability.

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.

Properties of the density of shear transformations in driven amorphous solids

The strain load Δγthat triggers consecutive avalanches is a key observable in the slow deformation of amorphous solids. Its temporally averaged value ⟨Δγ⟩ displays a non-trivial system-size dependence that constitutes one of the distinguishing features of the yielding transition. Details of this dependence are not yet fully understood. We address this problem by means of theoretical analysis and simulations of elastoplastic models for amorphous solids. An accurate determination of the size dependence of ⟨Δγ⟩ leads to a precise evaluation of the steady-state distribution of local distances to instabilityx. We find that the usually assumed formP(x) ∼x^θ(withθbeing the so-called pseudo-gap exponent) is not accurate at lowxand that in generalP(x) tends to a system-size-dependentfinitelimit asx→ 0. We work out the consequences of this finite-size dependence standing on exact results for random-walks and disclosing an alternative interpretation of the mechanical noise felt by a reference site. We test our predictions in two- and three-dimensional elastoplastic models, showing the crucial influence of the saturation ofP(x) at smallxon the size dependence of ⟨Δγ⟩ and related scalings.

Continuum constitutive laws to describe acoustic attenuation in glasses

Nowadays metamaterials are at the focus of an intense research as promising for thermal and acoustic engineering. However, the computational cost associated to the large system size required for correctly simulating them imposes the use of finite-elements simulations, developing continuum models, able to grasp the physics at play without entering in the atomistic details. Still, a correct description should be able to reproduce not only the extrinsic scattering sources on waves propagation, as introduced by the metamaterial microstructure, but also the intrinsic wave attenuation of the material itself. This becomes dramatically important when the metamaterial is made out of a glass, which is intrinsically highly dissipative and with a wave attenuation strongly dependent on frequency. Here we propose a continuum mechanical model for a viscoelastic medium, able to bridge atomic and macroscopic scale in amorphous materials and describe phonon attenuation due to atomistic mechanisms, characterized by a defined frequency dependence. This represents a first decisive step for investigating the effect of a complex nano- or microstructure on acoustic attenuation, while including the atomistic contribution as well.

Vibrational spectrum derived from local mechanical response in disordered solids

The low-frequency vibrations of glasses are markedly different from those of crystals. These vibrations have recently been categorized into two types: spatially extended vibrations, whose vibrational density of states (vDOS) follows a non-Debye quadratic law, and quasilocalized vibrations (QLVs), whose vDOS follows a quartic law. The former are explained by elasticity theory with quenched disorder and microscopic replica theory as being a consequence of elastic instability, but the origin of the latter is still debated. Here, we show that the latter can also be directly derived from elasticity theory with quenched disorder. We find another elastic instability that the theory encompasses but that has been overlooked so far, namely, the instability of the system against a local dipolar force. This instability gives rise to an additional contribution to the vDOS, and the spatial structure and energetics of the mode originating from this instability are consistent with those of the QLVs. Finally, we construct a model in which the additional contribution to the vDOS follows a quartic law.

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.

Renormalization of elastic quadrupoles in amorphous solids

Plasticity in amorphous solids is mediated by localized quadrupolar instabilities, but the mechanism by which an amorphous solid eventually fails or melts is debated. In this work we argue that these phenomena can be investigated in the model problem of an elastic continuum with quadrupolar defects, at finite temperature. This problem is posed and the collective behavior of the defects is analytically investigated. Using both renormalization group and field-theoretic techniques, it is found that the model has a yielding/melting transition of spinodal type.

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.

Nonlinear quasilocalized excitations in glasses: True representatives of soft spots

Structural glasses formed by quenching a melt possess a population of soft quasilocalized excitations-often called "soft spots"-that are believed to play a key role in various thermodynamic, transport, and mechanical phenomena. Under a narrow set of circumstances, quasilocalized excitations assume the form of vibrational (normal) modes, that are readily obtained by a harmonic analysis of the multidimensional potential energy. In general, however, direct access to the population of quasilocalized modes via harmonic analysis is hindered by hybridizations with other low-energy excitations, e.g., phonons. In this series of papers we reintroduce and investigate the statistical-mechanical properties of a class of low-energy quasilocalized modes-coined here nonlinear quasilocalized excitations (NQEs)-that are defined via an anharmonic (nonlinear) analysis of the potential-energy landscape of a glass, and do not hybridize with other low-energy excitations. In this paper, we review the theoretical framework that embeds a micromechanical definition of NQEs. We demonstrate how harmonic quasilocalized modes hybridize with other soft excitations, whereas NQEs properly represent soft spots without hybridization. We show that NQEs' energies converge to the energies of the softest, nonhybridized harmonic quasilocalized modes, cementing their status as true representatives of soft spots in structural glasses. Finally, we perform a statistical analysis of the mechanical properties of NQEs, which results in a prediction for the distribution of potential-energy barriers that surround typical inherent states of structural glasses, as well as a prediction for the distribution of local strain thresholds to plastic instability.

Universal non-mean-field scaling in the density of states of amorphous solids

Amorphous solids have excess soft modes in addition to the phonon modes described by the Debye theory. Recent numerical results show that if the phonon modes are carefully removed, the density of state of the excess soft modes exhibit universal quartic scaling, independent of the interaction potential, preparation protocol, and spatial dimensions. We hereby provide a theoretical framework to describe this universal scaling behavior. For this purpose, we extend the mean-field theory to include the effects of finite-dimensional fluctuation. Based on a semiphenomenological argument, we show that mean-field quadratic scaling is replaced by the quartic scaling in finite dimensions. Furthermore, we apply our formalism to explain the pressure and protocol dependence of the excess soft modes.