Effect of instantaneous and continuous quenches on the density of vibrational modes in model glasses

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.

Self-driven configurational dynamics in frustrated spring-mass systems

Various physical systems relax mechanical frustration through configurational rearrangements. We examine such rearrangements via Hamiltonian dynamics of simple internally stressed harmonic four-mass systems. We demonstrate theoretically and numerically how mechanical frustration controls the underlying potential energy landscape. Then, we examine the harmonic four-mass systems' Hamiltonian dynamics and relate the onset of chaotic motion to self-driven rearrangements. We show such configurational dynamics may occur without strong precursors, rendering such dynamics seemingly spontaneous.

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.

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}.

Internally Stressed and Positionally Disordered Minimal Complexes Yield Glasslike Nonphononic Excitations

Glasses, unlike their crystalline counterparts, exhibit low-frequency nonphononic excitations whose frequencies ω follow a universal D(ω)∼ω^{4} density of states. The process of glass formation generates positional disorder intertwined with mechanical frustration, posing fundamental challenges in understanding the origins of glassy nonphononic excitations. Here we suggest that minimal complexes-mechanically frustrated and positionally disordered local structures-embody the minimal physical ingredients needed to generate glasslike excitations. We investigate the individual effects of mechanical frustration and positional disorder on the vibrational spectrum of isolated minimal complexes, and demonstrate that ensembles of marginally stable minimal complexes yield D(ω)∼ω^{4}. Furthermore, glasslike excitations emerge by embedding a single minimal complex within a perfect lattice. Consequently, minimal complexes offer a conceptual framework to understand glasslike excitations from first principles, as well as a practical computational method for introducing them into solids.

Mechanical disorder of sticky-sphere glasses. II. Thermomechanical inannealability

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history dependence and what it is possibly affected by are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles-referred to here as glass "stickiness." We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses; the latter are therefore deemed "thermomechanically inannealable." The implications of our findings and future research directions are discussed.

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.

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.

Vibrational density of states of amorphous solids with long-ranged power-law-correlated disorder in elasticity

A theory of vibrational excitations based on power-law spatial correlations in the elastic constants (or equivalently in the internal stress) is derived, in order to determine the vibrational density of states D([Formula: see text]) of disordered solids. The results provide the first prediction of a boson peak in amorphous materials where spatial correlations in the internal stresses (or elastic constants) are of power-law form, as is often the case in experimental systems, leading to a logarithmic enhancement of (Rayleigh) phonon attenuation. A logarithmic correction of the form [Formula: see text] is predicted to occur in the plot of the reduced excess DOS for frequencies around the boson peak in 3D. Moreover, the theory provides scaling laws of the density of states in the low-frequency region, including a [Formula: see text] regime in 3D, and provides information about how the boson peak intensity depends on the strength of power-law decay of fluctuations in elastic constants or internal stress. Analytical expressions are also derived for the dynamic structure factor for longitudinal excitations, which include a logarithmic correction factor, and numerical calculations are presented supporting the assumptions used in the theory.

Key role of retardation and non-locality in sound propagation in amorphous solids as evidenced by a projection formalism

We investigate acoustic propagation in amorphous solids by constructing a projection formalism based on separating atomic vibrations into two, "phonon" (P) and "non-phonon" (NP), subspaces corresponding to large and small wavelengths. For a pairwise interaction model, we show the existence of a "natural" separation lengthscale, determined by structural disorder, for which the isolated P subspace presents the acoustic properties of a nearly homogenous (Debye-like) elastic continuum, while the NP one encapsulates all small scale non-affinity effects. The NP eigenstates then play the role of dynamical scatterers for the phonons. However, at variance with a conjecture of defect theories, their spectra present a finite low frequency gap, which turns out to lie around the Boson peak frequency, and only a small fraction of them are highly localized. We then show that small scale disorder effects can be rigorously reduced to the existence, in the Navier-like wave equation of the continuum, of a generalized elasticity tensor, which is not only retarded, since scatterers are dynamical, but also non-local. The full neglect of both retardation and non-locality suffices to account for most of the corrections to Born macroscopic moduli. However, these two features are responsible for sound speed dispersion and have quite a significant effect on the magnitude of sound attenuation. Although it remains open how they impact the asymptotic, large wavelength scaling of sound damping, our findings rule out the possibility of representing an amorphous solid by an inhomogeneous elastic continuum with the standard (i.e., local and static) elastic moduli.

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.

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.

Low-frequency vibrations of jammed packings in large spatial dimensions

Amorphous packings prepared in the vicinity of the jamming transition play a central role in theoretical studies of the vibrational spectrum of glasses. Two mean-field theories predict that the vibrational density of states g(ω) obeys a characteristic power law, g(ω)∼ω^{2}, called the non-Debye scaling in the low-frequency region. Numerical studies have, however, reported that this scaling breaks down at low frequencies, due to finite-dimensional effects. In this study, we prepare amorphous packings of up to 128000 particles in spatial dimensions from d=3 to d=9 to characterize the range of validity of the non-Debye scaling. Our numerical results suggest that the non-Debye scaling is obeyed down to a frequency that gradually decreases as d increases, and possibly vanishes for large d, in agreement with mean-field predictions. We also show that the prestress is an efficient control parameter to quantitatively compare packings across different spatial dimensions.

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.

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

The universal form of the density of nonphononic, quasilocalized vibrational modes of frequency ω in structural glasses, D(ω), was predicted theoretically decades ago, but only recently revealed in numerical simulations. In particular, it has been recently established that, in generic computer glasses, D(ω) increases from zero frequency as ω^{4}, independent of spatial dimension and of microscopic details. However, it has been shown [Lerner and Bouchbinder, Phys. Rev. E 96, 020104(R) (2017)2470-004510.1103/PhysRevE.96.020104] that the preparation protocol employed to create glassy samples may affect the form of their resulting D(ω): glassy samples rapidly quenched from high-temperature liquid states were shown to feature D(ω)∼ω^{β} with β<4, presumably limiting the degree of universality of the ω^{4} law. Here we show that exponents β<4 are seen only in small glassy samples quenched from high-temperature liquid states-whose sizes are comparable to or smaller than the size of the disordered core of soft quasilocalized vibrations-while larger glassy samples made with the same protocol feature the universal ω^{4} law. Our results demonstrate that observations of β<4 in the nonphononic density of states stem from finite-size effects, and we thus conclude that the ω^{4} law should be featured by any sufficiently large glass quenched from a melt.

Energy transport in glasses

The temperature dependence of the thermal conductivity is linked to the nature of the energy transport at a frequency ω, which is quantified by thermal diffusivity d(ω). Here we study d(ω) for a poorly annealed glass and a highly stable glass prepared using the swap Monte Carlo algorithm. To calculate d(ω), we excite wave packets and find that the energy moves diffusively for high frequencies up to a maximum frequency, beyond which the energy stays localized. At intermediate frequencies, we find a linear increase of the square of the width of the wave packet with time, which allows for a robust calculation of d(ω), but the wave packet is no longer well described by a Gaussian as for high frequencies. In this intermediate regime, there is a transition from a nearly frequency independent thermal diffusivity at high frequencies to d(ω) ∼ ω-4 at low frequencies. For low frequencies the sound waves are responsible for energy transport and the energy moves ballistically. The low frequency behavior can be predicted using sound attenuation coefficients.

Relation between Heterogeneous Frozen Regions in Supercooled Liquids and Non-Debye Spectrum in the Corresponding Glasses

Recent numerical studies on glassy systems provide evidence for a population of non-Goldstone modes (NGMs) in the low-frequency spectrum of the vibrational density of states D(ω). Similarly to Goldstone modes (GMs), i.e., phonons in solids, NGMs are soft low-energy excitations. However, differently from GMs, NGMs are localized excitations. Here we first show that the parental temperature T^{*} modifies the GM/NGM ratio in D(ω). In particular, the phonon attenuation is reflected in a parental temperature dependency of the exponent s(T^{*}) in the low-frequency power law D(ω)∼ω^{s(T^{*})}, with 2≤s(T^{*})≤4. Second, by comparing s(T^{*}) with s(p), i.e., the same quantity obtained by pinning a p particle fraction, we suggest that s(T^{*}) reflects the presence of dynamical heterogeneous regions of size ξ^{3}∝p. Finally, we provide an estimate of ξ as a function of T^{*}, finding a mild power law divergence, ξ∼(T^{*}-T_{d})^{-α/3}, with T_{d} the dynamical crossover temperature and α falling in the range α∈[0.8,1.0].

Impact of jamming criticality on low-temperature anomalies in structural glasses

We present a mechanism for the anomalous behavior of the specific heat in low-temperature amorphous solids. The analytic solution of a mean-field model belonging to the same universality class as high-dimensional glasses, the spherical perceptron, suggests that there exists a cross-over temperature above which the specific heat scales linearly with temperature, while below it, a cubic scaling is displayed. This relies on two crucial features of the phase diagram: (i) the marginal stability of the free-energy landscape, which induces a gapless phase responsible for the emergence of a power-law scaling; and (ii) the vicinity of the classical jamming critical point, as the cross-over temperature gets lowered when approaching it. This scenario arises from a direct study of the thermodynamics of the system in the quantum regime, where we show that, contrary to crystals, the Debye approximation does not hold.

Anisotropic structural predictor in glassy materials

There is growing evidence that relaxation in glassy materials, both spontaneous and externally driven, is mediated by localized soft spots. Recent progress made it possible to identify the soft spots inside glassy structures and to quantify their degree of softness. These softness measures, however, are typically scalars, not taking into account the tensorial, anisotropic nature of soft spots, which implies orientation-dependent coupling to external deformation. Here, we derive from first principles the linear response coupling between the local heat capacity of glasses, previously shown to provide a measure of glassy softness, and external deformation in different directions. We first show that this linear response quantity follows an anomalous, fat-tailed distribution related to the universal ω^{4} density of states of quasilocalized, nonphononic excitations in glasses. We then construct a structural predictor as the product of the local heat capacity and its linear response to external deformation, and show that it offers an enhanced predictability of plastic rearrangements under deformation in different directions, compared to the purely scalar predictor.

Revealing Inherent Structural Characteristics of Jammed Particulate Packings

We look for inherent structural characteristics hidden behind amorphous solid formation by using zero-temperature jammed packings of frictionless particles as models. Differently from previous geometrical approaches, we introduce a microscopic mechanical or vibrational order parameter Ψ, which characterizes the susceptibility of particle motion to infinitesimal thermal excitation. We show that (i) the distribution of Ψ has a power-law tail toward high Ψ and (ii) the spatial organization of Ψ is characterized by a nontrivial scale-free correlation. Both findings (i) and (ii) are regarded as a real-space manifestation of marginal stability due to critical self-organization of jammed packings toward mechanical equilibrium. Furthermore, we find that the power-law exponent of the Ψ distribution tail shows a critical-like scaling behavior toward the unjamming transition, which unveils an intriguing interplay between jamming criticality and marginal stability. Our microscopic order parameter provides new structural insights into the marginal stability and instability of jammed packings and may shed light on the important common structural feature of amorphous solids.

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.

Theory for the density of interacting quasilocalized modes in amorphous solids

Quasilocalized modes appear in the vibrational spectrum of amorphous solids at low frequency. Though never formalized, these modes are believed to have a close relationship with other important local excitations, including shear transformations and two-level systems. We provide a theory for their frequency density, D_{L}(ω)∼ω^{α}, that establishes this link for systems at zero temperature under quasistatic loading. It predicts two regimes depending on the density of shear transformations P(x)∼x^{θ} (with x the additional stress needed to trigger a shear transformation). If θ>1/4, then α=4 and a finite fraction of quasilocalized modes form shear transformations, whose amplitudes vanish at low frequencies. If θ<1/4, then α=3+4θ and all quasilocalized modes form shear transformations with a finite amplitude at vanishing frequencies. We confirm our predictions numerically.

Low-frequency vibrational modes of stable glasses

Unusual features of the vibrational density of states D(ω) of glasses allow one to rationalize their peculiar low-temperature properties. Simulational studies of D(ω) have been restricted to studying poorly annealed glasses that may not be relevant to experiments. Here we report on D(ω) of zero-temperature glasses with kinetic stabilities ranging from poorly annealed to ultrastable glasses. For all preparations, the low-frequency part of D(ω) splits between extended and quasi-localized modes. Extended modes exhibit a boson peak crossing over to Debye behavior (Dex(ω) ~ ω2) at low-frequency, with a strong correlation between the two regimes. Quasi-localized modes obey Dloc(ω) ~ ω4, irrespective of the stability. The prefactor of this quartic law decreases with increasing stability, and the corresponding modes become more localized and sparser. Our work is the first numerical observation of quasi-localized modes in a regime relevant to experiments, and it establishes a direct connection between glasses' stability and their soft vibrational modes.

Probing the non-Debye low-frequency excitations in glasses through random pinning

We investigate the properties of the low-frequency spectrum in the density of states [Formula: see text] of a 3D model glass former. To magnify the non-Debye sector of the spectrum, we introduce a random pinning field that freezes a finite particle fraction to break the translational invariance and shifts all of the vibrational frequencies of the extended modes toward higher frequencies. We show that non-Debye soft localized modes progressively emerge as the fraction p of pinned particles increases. Moreover, the low-frequency tail of [Formula: see text] goes to zero as a power law [Formula: see text], with [Formula: see text] and [Formula: see text] above a threshold fraction [Formula: see text].

Universal Nonphononic Density of States in 2D, 3D, and 4D Glasses

It is now well established that structural glasses possess disorder- and frustration-induced soft quasilocalized excitations, which play key roles in various glassy phenomena. Recent work has established that in model glass formers in three dimensions, these nonphononic soft excitations may assume the form of quasilocalized, harmonic vibrational modes whose frequency follows a universal density of states D(ω)∼ω^{4}, independently of microscopic details, and for a broad range of glass preparation protocols. Here, we further establish the universality of the nonphononic density of vibrational modes by direct measurements in model structural glasses in two dimensions and four dimensions. We also investigate their degree of localization, which is generally weaker in lower spatial dimensions, giving rise to a pronounced system-size dependence of the nonphononic density of states in two dimensions, but not in higher dimensions. Finally, we identify a fundamental glassy frequency scale ω_{c} above which the universal ω^{4} law breaks down.

Mean-field model for the density of states of jammed soft spheres

We propose a class of mean-field models for the isostatic transition of systems of soft spheres, in which the contact network is modeled as a random graph and each contact is associated to d degrees of freedom. We study such models in the hypostatic, isostatic, and hyperstatic regimes. The density of states is evaluated by both the cavity method and exact diagonalization of the dynamical matrix. We show that the model correctly reproduces the main features of the density of states of real packings and, moreover, it predicts the presence of localized modes near the lower band edge. Finally, the behavior of the density of states D(ω)∼ω^{α} for ω→0 in the hyperstatic regime is studied. We find that the model predicts a nontrivial dependence of α on the details of the coordination distribution.

Escape time, relaxation, and sticky states of a softened Henon-Heiles model: Low-frequency vibrational mode effects and glass relaxation

Here we study the relaxation of a chain consisting of three masses joined by nonlinear springs and periodic conditions when the stiffness is weakened. This system, when expressed in their normal coordinates, yields a softened Henon-Heiles system. By reducing the stiffness of one low-frequency vibrational mode, a faster relaxation is enabled. This is due to a reduction of the energy barrier heights along the softened normal mode as well as for a widening of the opening channels of the energy landscape in configurational space. The relaxation is for the most part exponential, and can be explained by a simple flux equation. Yet, for some initial conditions the relaxation follows as a power law, and in many cases there is a regime change from exponential to power-law decay. We pinpoint the initial conditions for the power-law decay, finding two regions of sticky states. For such states, quasiperiodic orbits are found since almost for all components of the initial momentum orientation, the system is trapped inside two pockets of configurational space. The softened Henon-Heiles model presented here is intended as the simplest model in order to understand the interplay of rigidity, nonlinear interactions and relaxation for nonequilibrium systems such as glass-forming melts or soft matter. Our softened system can be applied to model β relaxation in glasses and suggest that local reorientational jumps can have an exponential and a nonexponential contribution for relaxation, the latter due to asymmetric molecules sticking in cages for certain orientations.

Frustration-induced internal stresses are responsible for quasilocalized modes in structural glasses

It has been recently shown [E. Lerner, G. Düring, and E. Bouchbinder, Phys. Rev. Lett. 117, 035501 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.035501] that the nonphononic vibrational modes of structural glasses at low frequencies ω are quasilocalized and follow a universal density of states D(ω)∼ω^{4}. Here we show that the gapless nature of the observed density of states depends on the existence of internal stresses that generically emerge in glasses due to frustration, thus elucidating a basic element underlying this universal behavior. Similarly to jammed particulate packings, low-frequency modes in structural glasses emerge from a balance between a local elasticity term and an internal stress term in the dynamical matrix, where the difference between them is orders of magnitude smaller than their typical magnitude. By artificially reducing the magnitude of internal stresses in a computer glass former in three dimensions, we show that a gap is formed in the density of states below which no vibrational modes exist, thus demonstrating the crucial importance of internal stresses. Finally, we show that while better annealing the glass upon cooling from the liquid state significantly reduces its internal stresses, the self-organizational processes during cooling render the gapless D(ω)∼ω^{4} density of state unaffected.

Anomalous vibrational properties in the continuum limit of glasses

The low-temperature thermal properties of glasses are anomalous with respect to those of crystals. These thermal anomalies indicate that the low-frequency vibrational properties of glasses differ from those of crystals. Recent studies revealed that, in the simplest model of glasses, i.e., the harmonic potential system, phonon modes coexist with soft localized modes in the low-frequency (continuum) limit. However, the nature of low-frequency vibrational modes of more realistic models is still controversial. In the present work, we study the Lennard-Jones (LJ) system using large-scale molecular-dynamics (MD) simulation and establish that the vibrational property of the LJ glass converges to coexistence of the phonon modes and the soft localized modes in the continuum limit as in the case of the harmonic potential system. Importantly, we find that the low-frequency vibrations are rather sensitive to the numerical scheme of potential truncation, which is usually implemented in the MD simulation, and this is the reason why contradictory arguments have been reported by previous works. We also discuss the physical origin of this sensitiveness by means of a linear stability analysis.

Instabilities of Jammed Packings of Frictionless Spheres Under Load

We consider the contribution to the density of vibrational states and the distribution of energy barrier heights of incipient instabilities in a glass modeled by a jammed packing of spheres. On approaching an instability, the frequency of a normal mode and the height of the energy barrier to cross into a new ground state both vanish. These instabilities produce a contribution to the density of vibrational states that scales as ω^{3} at low frequencies ω, and a contribution to the distribution of energy barriers ΔH that scales as ΔH^{-1/3} at low barrier heights.

Erratum: Statistics and Properties of Low-Frequency Vibrational Modes in Structural Glasses [Phys. Rev. Lett. 117, 035501 (2016)]

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