Finite-size effects in the nonphononic density of states in computer glasses

Testing the heterogeneous-elasticity theory for low-energy excitations in structural glasses

Understanding the statistical mechanics of low-energy excitations in structural glasses has been the focus of extensive research efforts in the past decades due to their key roles in determining the low-temperature mechanical and transport properties of these intrinsically nonequilibrium materials. While it is established that glasses feature low-energy nonphononic excitations that follow a non-Debye vibrational density of states, we currently lack a well-founded theory of these fundamental objects and their vibrational spectra. A recent theory-that builds on the so-called heterogeneous-elasticity theory (HET) and its extensions-provides explicit predictions for the scaling of the low-frequency tail of the nonphononic spectrum of glasses, the localization properties of the vibrational modes that populate this tail, and its connections to glass formation history and to the form of the distribution of weak microscopic (interatomic) stresses. Here, we employ computer models of structural glasses to quantitatively test these predictions. Our findings do not support the HET's predictions regarding the nature and statistics of low-energy excitations in glasses, highlighting the need for additional theoretical developments.

Enumerating low-frequency nonphononic vibrations in computer glasses

In addition to Goldstone phonons that generically emerge in the low-frequency vibrational spectrum of any solid, crystalline or glassy, structural glasses also feature other low-frequency vibrational modes. The nature and statistical properties of these modes-often termed "excess modes"-have been the subject of decades-long investigation. Studying them, even using well-controlled computer glasses, has proven challenging due to strong spatial hybridization effects between phononic and nonphononic excitations, which hinder quantitative analyses of the nonphononic contribution DG(ω) to the total spectrum D(ω), per frequency ω. Here, using recent advances indicating that DG(ω)=D(ω)-DD(ω), where DD(ω) is Debye's spectrum of phonons, we present a simple and straightforward scheme to enumerate nonphononic modes in computer glasses. Our analysis establishes that nonphononic modes in computer glasses indeed make an additive contribution to the total spectrum, including in the presence of strong hybridizations. Moreover, it cleanly reveals the universal DG(ω)∼ω4 tail of the nonphononic spectrum, and opens the way for related analyses of experimental spectra of glasses.

Elasticity of self-organized frustrated disordered spring networks

There have been some interesting recent advances in understanding the notion of mechanical disorder in structural glasses and the statistical mechanics of these systems' low-energy excitations. Here we contribute to these advances by studying a minimal model for structural glasses' elasticity in which the degree of mechanical disorder-as characterized by recently introduced dimensionless quantifiers-is readily tunable over a very large range. We comprehensively investigate a number of scaling laws observed for various macro, meso and microscopic elastic properties, and rationalize them using scaling arguments. Interestingly, we demonstrate that the model features the universal quartic glassy vibrational density of states as seen in many atomistic and molecular models of structural glasses formed by cooling a melt. The emergence of this universal glassy spectrum highlights the role of self-organization (toward mechanical equilibrium) in its formation, and elucidates why models featuring structural frustration alone do not feature the same universal glassy spectrum. Finally, we discuss relations to existing work in the context of strain stiffening of elastic networks and of low-energy excitations in structural glasses, in addition to future research directions.

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.

Detecting low-energy quasilocalized excitations in computer glasses

Soft, quasilocalized excitations (QLEs) are known to generically emerge in a broad class of disordered solids and to govern many facets of the physics of glasses, from wave attenuation to plastic instabilities. In view of this key role of QLEs, shedding light upon several open questions in glass physics depends on the availability of computational tools that allow one to study QLEs' statistical mechanics. The latter is a formidable task since harmonic analyses are typically contaminated by hybridizations of QLEs with phononic excitations at low frequencies, obscuring a clear picture of QLEs' abundance, typical frequencies, and other important micromechanical properties. Here we present an efficient algorithm to detect the field of quasilocalized excitations in structural computer glasses. The algorithm introduced takes a computer-glass sample as input and outputs a library of QLEs embedded in that sample. We demonstrate the power of the algorithm by reporting the spectrum of glassy excitations in two-dimensional computer glasses featuring a huge range of mechanical stability, which is inaccessible using conventional harmonic analyses due to phonon hybridizations. Future applications are discussed.

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.

Non-phononic density of states of two-dimensional glasses revealed by random pinning

The vibrational density of states of glasses is considerably different from that of crystals. In particular, there exist spatially localized vibrational modes in glasses. The density of states of these non-phononic modes has been observed to follow g(ω) ∝ ω4, where ω is the frequency. However, in two-dimensional systems, the abundance of phonons makes it difficult to accurately determine this non-phononic density of states because they are strongly coupled to non-phononic modes and yield strong system-size and preparation-protocol dependencies. In this article, we utilize the random pinning method to suppress phonons and disentangle their coupling with non-phononic modes and successfully calculate their density of states as g(ω) ∝ ω4. We also study their localization properties and confirm that low-frequency non-phononic modes in pinned systems are truly localized without far-field contributions. We finally discuss the excess density of states over the Debye value that results from the hybridization of phonons and non-phononic modes.

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.

Disordered Crystals Reveal Soft Quasilocalized Glassy Excitations

Structural glasses formed by quenching a melt are known to host a population of low-energy quasilocalized (nonphononic) excitations whose frequencies ω follow a universal ∼ω^{4} distribution as ω→0, independently of the glass formation history, the interparticle interaction potential, or spatial dimension. Here, we show that the universal quartic law of nonphononic excitations also holds in disordered crystals featuring finite long-range order, which is absent in their glassy counterparts. We thus establish that the degree of universality of the quartic law extends beyond structural glasses quenched from a melt. We further find that disordered crystals, whose level of disorder can be continuously controlled, host many more quasilocalized excitations than expected based on their degree of mechanical disorder-quantified by the relative fluctuations of the shear modulus-as compared to structural glasses featuring a similar degree of mechanical disorder. Our results are related to glasslike anomalies experimentally observed in disordered crystals. More broadly, they constitute an important step toward tracing the essential ingredients necessary for the emergence of universal nonphononic excitations in disordered solids.

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.

Erratum: Low-Frequency Excess Vibrational Modes in Two-Dimensional Glasses [Phys. Rev. Lett. 127, 248001 (2021)]

This corrects the article DOI: 10.1103/PhysRevLett.127.248001.

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.

Bond-space operator disentangles quasilocalized and phononic modes in structural glasses

The origin of several emergent mechanical and dynamical properties of structural glasses is often attributed to populations of localized structural instabilities, coined quasilocalized modes (QLMs). Under a restricted set of circumstances, glassy QLMs can be revealed by analyzing computer glasses' vibrational spectra in the harmonic approximation. However, this analysis has limitations due to system-size effects and hybridization processes with low-energy phononic excitations (plane waves) that are omnipresent in elastic solids. Here we overcome these limitations by exploring the spectrum of a linear operator defined on the space of particle interactions (bonds) in a disordered material. We find that this bond-force-response operator offers a different interpretation of QLMs in glasses and cleanly recovers some of their important statistical and structural features. The analysis presented here reveals the dependence of the number density (per frequency) and spatial extent of QLMs on material preparation protocol (annealing). Finally, we discuss future research directions and possible extensions of this work.

Unified quantifier of mechanical disorder in solids

Mechanical disorder in solids, which is generated by a broad range of physical processes and controls various material properties, appears in a wide variety of forms. Defining unified and measurable dimensionless quantifiers, allowing quantitative comparison of mechanical disorder across widely different physical systems, is therefore an important goal. Two such coarse-grained dimensionless quantifiers (among others) appear in the literature: one is related to the spectral broadening of discrete phononic bands in finite-size systems (accessible through computer simulations) and the other is related to the spatial fluctuations of the shear modulus in macroscopically large systems. The latter has been recently shown to determine the amplitude of wave attenuation rates in the low-frequency limit (accessible through laboratory experiments). Here, using two alternative and complementary theoretical approaches linked to the vibrational spectra of solids, we derive a basic scaling relation between the two dimensionless quantifiers. This scaling relation, which is supported by simulational data, shows that the two apparently distinct quantifiers are in fact intrinsically related, giving rise to a unified quantifier of mechanical disorder in solids. We further discuss the obtained results in the context of the unjamming transition taking place in soft sphere packings at low confining pressures, in addition to their implications for our understanding of the low-frequency vibrational spectra of disordered solids in general, and in particular those of glassy systems.

Mechanical disorder of sticky-sphere glasses. I. Effect of attractive interactions

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions-referred to as "glass stickiness"-can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, negative stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

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.

Universality of the Nonphononic Vibrational Spectrum across Different Classes of Computer Glasses

It has been recently established that the low-frequency spectrum of simple computer glass models is populated by soft, quasilocalized nonphononic vibrational modes whose frequencies ω follow a gapless, universal distribution D(ω)∼ω^{4}. While this universal nonphononic spectrum has been shown to be robust to varying the glass history and spatial dimension, it has so far only been observed in simple computer glasses featuring radially symmetric, pairwise interaction potentials. Consequently, the relevance of the universality of nonphononic spectra seen in simple computer glasses to realistic laboratory glasses remains unclear. Here, we demonstrate the emergence of the universal ω^{4} nonphononic spectrum in a broad variety of realistic computer glass models, ranging from tetrahedral network glasses with three-body interactions, through molecular glasses and glassy polymers, to bulk metallic glasses. Taken together with previous observations, our results indicate that the low-frequency nonphononic vibrational spectrum of any glassy solid quenched from a melt features the universal ω^{4} law, independently of the nature of its microscopic interactions.

Statistical mechanics of local force dipole responses in computer glasses

Soft quasilocalized modes (QLMs) are universally featured by structural glasses quenched from a melt, and are involved in several glassy anomalies such as the low-temperature scaling of their thermal conductivity and specific heat, and sound attenuation at intermediate frequencies. In computer glasses, QLMs may assume the form of harmonic vibrational modes under a narrow set of circumstances; however, direct access to their full distribution over frequency is hindered by hybridizations of QLMs with other low-frequency modes (e.g., phonons). Previous studies to overcome this issue have demonstrated that the response of a glass to local force dipoles serves as a good proxy for its QLMs; we, therefore, study here the statistical-mechanical properties of these responses in computer glasses, over a large range of glass stabilities and in various spatial dimensions, with the goal of revealing properties of the yet-inaccessible full distribution of QLMs' frequencies. We find that as opposed to the spatial-dimension-independent universal distribution of QLMs' frequencies ω (and, consequently, also of their stiffness κ = ω²), the distribution of stiffnesses associated with responses to local force dipoles features a (weak) dependence on spatial dimension. We rationalize this dependence by introducing a lattice model that incorporates both the real-space profiles of QLMs-associated with dimension-dependent long-range elastic fields-and the universal statistical properties of their frequencies. Based on our findings, we propose a conjecture about the form of the full distribution of QLMs' frequencies and its protocol-dependence. Finally, we discuss possible connections of our findings to basic aspects of glass formation and deformation.

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.