Soft Modes, Localization, and Two-Level Systems in Spin Glasses

Vibrational lifetimes and viscoelastic properties of ultrastable glasses

Amorphous solids are viscoelastic. They dissipate energy when deformed at finite rate and finite temperature. We here use analytic theory and molecular simulations to demonstrate that linear viscoelastic dissipation can be directly related to the static and dynamic properties of the fundamental vibrational excitations of an amorphous system. We study ultrastable glasses that do not age, i.e., that remain in stable minima of the potential energy surface at finite temperature. Our simulations show four types of vibrational modes, which differ in spatial localization, similarity to plane waves and vibrational lifetimes. At frequencies below the Boson peak, the viscoelastic response can be split into contributions from plane-wave and quasilocalized modes. We derive a parameter-free expression for the viscoelastic storage and loss moduli for both of these modes. Our results show that the dynamics of microscopic dissipation, in particular the lifetimes of the modes, determine the viscoelastic response only at high frequency. Quasilocalized modes dominate the linear viscoelastic response at intermediate frequencies below the Boson peak.

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.

Study of the de Almeida-Thouless line in the one-dimensional diluted power-law XY spin glass

We study the de Almeida-Thouless (AT) line in the one-dimensional power-law diluted XY spin-glass model, in which the probability that two spins separated by a distance r interact with each other, decays as 1/r^{2σ}. Tuning the exponent σ is equivalent to changing the space dimension of a short-range model. We develop a heat bath algorithm to equilibrate XY spins; using this in conjunction with the standard parallel tempering and overrelaxation sweeps, we carry out large-scale Monte Carlo simulations. For σ=0.6, which is in the mean-field regime above six dimensions-it is similar to being in 10 dimensions-we find clear evidence for an AT line. For σ=0.75 and σ=0.85, which are in the non-mean-field regime and similar to four and three dimensions, respectively, our data is like that found in previous studies of the Ising and Heisenberg spin glasses when reducing the temperature at fixed field. For σ=0.75, there is evidence from finite-size-scaling studies for an AT transition but for σ=0.85, the evidence for a transition is nonexistent. We have also studied these systems at fixed temperature varying the field and discovered that at both σ=0.75 and at σ=0.85 there is evidence of an AT transition! Confusingly, the correlation length and spin-glass susceptibility as a function of the field are both entirely consistent with the predictions of the droplet picture and hence the nonexistence of an AT line. In the usual finite-size critical point scaling studies used to provide evidence for an AT transition, there is seemingly good evidence for an AT line at σ=0.75 for small values of the system size N, which is strengthening as N is increased, but for N>2048 the trend changes and the evidence then weakens as N is further increased. We have also studied with fewer bond realizations the system at σ=0.70, which is the analog of a system with short-range interactions just below six dimensions, and found that it is similar in its behavior to the system at σ=0.75 but with larger finite-size corrections. The evidence from our simulations points to the complete absence of the AT line in dimensions outside the mean-field region and to the correctness of the droplet picture. Previous simulations which suggested there was an AT line can be attributed to the consequences of studying systems which are just too small. The collapse of our data to the droplet scaling form is poor for σ=0.75 and to some extent also for σ=0.85, when the correlation length becomes of the order of the length of the system, due to the existence of excitations which only cost a free energy of O(1), just as envisaged in the TNT picture of the ordered state of spin glasses. However, for the case of σ=0.85 we can provide evidence that for larger system sizes, droplet scaling will prevail even when the correlation length is comparable to the system size.

Gas-liquid phase separation at zero temperature: mechanical interpretation and implications for gelation

The relationship between glasses and gels has been intensely debated for decades; however, the transition between these two phases remains elusive. To investigate a gel formation process in the zero-temperature limit and its relation to the glass phase, we conducted numerical experiments on athermal quasistatic decompression. During decompression, the system experiences a cavitation event similar to phase separation and this is a gelation process at zero temperature. A normal mode analysis revealed that the phase separation is signaled by the vanishing of the lowest eigenenergy, similar to plastic events of glasses under shear. One primary difference from the shear-induced plasticity is that the vanishing mode experiences a qualitative change in its spatial energy distribution at the phase separation point. These findings enable us to define the glass-gel phase boundary based on mechanics.

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.

Universal density of low-frequency states in silica glass at finite temperatures

The theoretical understanding of the low-frequency modes in amorphous solids at finite temperature is still incomplete. The study of the relevant modes is obscured by the dressing of interparticle forces by collision-induced momentum transfer that is unavoidable at finite temperatures. Recently, it was proposed that low-frequency modes of vibrations around the thermally averaged configurations deserve special attention. In simple model glasses with bare binary interactions, these included quasilocalized modes whose density of states appears to be universal, depending on the frequencies as D(ω)∼ω^{4}, in agreement with the similar law that is obtained with bare forces at zero temperature. In this paper, we report investigations of a model of silica glass at finite temperature; here the bare forces include binary and ternary interactions. Nevertheless, we can establish the validity of the universal law of the density of quasilocalized modes also in this richer and more realistic model 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.

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

We combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field-nonaffinity-which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding "positive sound dispersion" phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

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.

Debye Brownian oscillator and Debye-type noise: A series solution versus Monte Carlo simulation

For the Debye Brownian oscillator, we present a series solution to the generalized Langevin equation describing the motion of a particle. The external potential is considered to be a harmonic potential and the spectral density of driven noise is a hard cutoff at high finite frequencies. The results are in agreement with both numerical calculations and Monte Carlo simulations. We demonstrate abnormal weak ergodic breaking; specifically, the long-time average of the observable vanishes but the corresponding ensemble average continues to oscillate with time. This Debye Brownian oscillator does not arrive at an equilibrium state and undergoes underdamped-like motion for any model parameter. Nevertheless, ergodic behavior and equilibrium can be recovered concurrently using a strong bound potential. We give an understanding of the behavior as being the consequence of discrete breather modes in the lattices similar to the formation of an additional periodic signal. Furthermore, we compare the results calculated by cutting off separately the spectral density and the correlation function of colored noise.

Spin glass behavior and magnetic boson peak in a structural glass of a magnetic ionic liquid

Glassy magnetic behavior has been observed in a wide range of crystalline magnetic materials called spin glass. Here, we report spin glass behavior in a structural glass of a magnetic ionic liquid, C4mimFeCl4. Magnetization measurements demonstrate that an antiferromagnetic ordering occurs at TN = 2.3 K in the crystalline state, while a spin glass transition occurs at TSG = 0.4 K in the structural glass state. In addition, localized magnetic excitations were found in the spin glass state by inelastic neutron scattering, in contrast to spin-wave excitations in the ordered phase of the crystalline sample. The localized excitation was scaled by the Bose population factor below TSG and gradually disappeared above TSG. This feature is highly reminiscent of boson peaks commonly observed in structural glasses. We suggest the "magnetic" boson peak to be one of the inherent dynamics of a spin glass state.

Universal Density of Low-Frequency States in Amorphous Solids at Finite Temperatures

It has been established that the low frequency quasilocalized modes of amorphous solids at zero temperature exhibit universal density of states, depending on the frequencies as D(ω)∼ω^{4}. It remains an open question whether this universal law extends to finite temperatures. In this Letter we show that well quenched model glasses at temperatures as high as T_{g}/3 possess the same universal density of states. The only condition required is that average particle positions stabilize before thermal diffusion destroys the cage structure of the material. The universal density of quasilocalized low frequency modes refers then to vibrations around the thermally averaged configuration of the material.

A structural approach to vibrational properties ranging from crystals to disordered systems

Many scientists generally attribute the vibrational anomalies of disordered solids to the structural disorder, which, however, is still under intense debate. Here we conduct simulations on two-dimensional packings with a finite temperature, whose structure is tuned from a crystalline configuration to an amorphous one, then the amorphous from very dense state to a relatively loose state. By measuring the vibrational density of states and the reduced density of states, we clearly observe the evolution of the boson peak with the change of the disorder and volume fractions. Meanwhile, to understand the structural origin of this anomaly, we identify the soft regimes of all systems with a novel machine-learning method, where the "softness", a local structural quantity, is defined. Interestingly, we find a strong monotonic relationship between the shape of the boson peak and the softness as well as its spatial heterogeneity, suggesting that the softness of a system may be a new structural approach to the anomalous vibrational properties of amorphous solids.

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.

Universal Low-Frequency Vibrational Modes in Silica Glasses

It was recently shown that different simple models of glass formers with binary interactions define a universality class in terms of the density of states of their quasilocalized low-frequency modes. Explicitly, once the hybridization with standard Debye (extended) modes is avoided, a number of such models exhibit a universal density of states, depending on the mode frequencies as D(ω)∼ω^{4}. It is unknown, however, how wide this universality class is, and whether it also pertains to more realistic models of glass formers. To address this issue we present analysis of the quasilocalized modes in silica, a network glass that has both binary and ternary interactions. We conclude that in three dimensions silica exhibits the very same frequency dependence at low frequencies, suggesting that this universal form is a generic consequence of amorphous glassiness.

Effective traplike activated dynamics in a continuous landscape

We use a simple model to extend network models for activated dynamics to a continuous landscape with a well-defined notion of distance and a direct connection to many-body systems. The model consists of a tracer in a high-dimensional funnel landscape with no disorder. We find a nonequilibrium low-temperature phase with aging dynamics that is effectively equivalent to that of models with built-in disorder, such as the trap model, step model and random energy model. Finally, we compare entropy with energy-driven activation, and we remark that the former is more robust to the choice of the dynamics since it does not depend on whether one uses local or global updates.

Vibrational Properties of Hard and Soft Spheres Are Unified at Jamming

The unconventional thermal properties of jammed amorphous solids are directly related to their density of vibrational states. While the vibrational spectrum of jammed soft sphere solids has been fully described, the vibrational spectrum of hard spheres, a model glass former often related to physical colloidal glasses, is still unknown due to the difficulty of treating the nonanalytic interaction potential. We bypass this difficulty using the recently described effective interaction potential for the free energy of thermal hard spheres. By minimizing this effective free energy, we mimic the rapid compression of hard spheres and produce typical configurations of the thermal system. We measure the resulting vibrational spectrum and characterize its evolution toward the jamming point where configurations of hard and soft spheres are trivially unified. For densities approaching jamming from below, we observe low-frequency modes which agree with those found in numerical simulations of jammed soft spheres. Our measurements of the vibrational structure demonstrate that the jamming universality extends away from jamming: hard sphere thermal systems below jamming exhibit the same vibrational spectra as thermal and athermal soft sphere systems above the transition.

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.

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.

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.

Fast generation of ultrastable computer glasses by minimization of an augmented potential energy

We present a model and protocol that enable the generation of extremely stable computer glasses at minimal computational cost. The protocol consists of an instantaneous quench in an augmented potential energy landscape, with particle radii as additional degrees of freedom. We demonstrate how our glasses' mechanical stability, which is readily tunable in our approach, is reflected in both microscopic and macroscopic observables. Our observations indicate that the stability of our computer glasses is at least comparable to that of computer glasses generated by the celebrated Swap Monte Carlo algorithm. Strikingly, some key properties support even qualitatively enhanced stability in our scheme: the density of quasilocalized excitations displays a gap in our most stable computer glasses, whose magnitude scales with the polydispersity of the particles. We explain this observation, which is consistent with the lack of plasticity we observe at small stress. It also suggests that these glasses are depleted from two-level systems, similarly to experimental vapor-deposited ultrastable glasses.

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.

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.

Continuum limit of the vibrational properties of amorphous solids

The low-frequency vibrational and low-temperature thermal properties of amorphous solids are markedly different from those of crystalline solids. This situation is counterintuitive because all solid materials are expected to behave as a homogeneous elastic body in the continuum limit, in which vibrational modes are phonons that follow the Debye law. A number of phenomenological explanations for this situation have been proposed, which assume elastic heterogeneities, soft localized vibrations, and so on. Microscopic mean-field theories have recently been developed to predict the universal non-Debye scaling law. Considering these theoretical arguments, it is absolutely necessary to directly observe the nature of the low-frequency vibrations of amorphous solids and determine the laws that such vibrations obey. Herein, we perform an extremely large-scale vibrational mode analysis of a model amorphous solid. We find that the scaling law predicted by the mean-field theory is violated at low frequency, and in the continuum limit, the vibrational modes converge to a mixture of phonon modes that follow the Debye law and soft localized modes that follow another universal non-Debye scaling law.

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

Computational studies of supercooled liquids often focus on various analyses of their "underlying inherent states"-the glassy configurations at zero temperature obtained by an infinitely fast (instantaneous) quench from equilibrium supercooled states. Similar protocols are also regularly employed in investigations of the unjamming transition at which the rigidity of decompressed soft-sphere packings is lost. Here we investigate the statistics and localization properties of low-frequency vibrational modes of glassy configurations obtained by such instantaneous quenches. We show that the density of vibrational modes grows as ω^{β} with β depending on the parent temperature T_{0} from which the glassy configurations were instantaneously quenched. For quenches from high temperature liquid states we find β≈3, whereas β appears to approach the previously observed value β=4 as T_{0} approaches the glass transition temperature. We discuss the consistency of our findings with the theoretical framework of the soft potential model, and contrast them with similar measurements performed on configurations obtained by continuous quenches at finite cooling rates. Our results suggest that any physical quench at rates sufficiently slower than the inverse vibrational time scale-including all physically realistic quenching rates of molecular or atomistic glasses-would result in a glass whose density of vibrational modes is universally characterized by β=4.

Metastable minima of the Heisenberg spin glass in a random magnetic field

We have studied zero-temperature metastable minima in classical m-vector component spin glasses in the presence of m-component random fields for two models, the Sherrington-Kirkpatrick (SK) model and the Viana-Bray (VB) model. For the SK model we have calculated analytically its complexity (the log of the number of minima) for both the annealed case where one averages the number of minima before taking the log and the quenched case where one averages the complexity itself, both for fields above and below the de Almeida-Thouless (AT) field, which is finite for m>2. We have done numerical quenches starting from a random initial state (infinite temperature state) by putting spins parallel to their local fields until there is no further decrease of the energy and found that in zero field it always produces minima that have zero overlap with each other. For the m=2 and m=3 cases in the SK model the final energy reached in the quench is very close to the energy E_{c} at which the overlap of the states would acquire replica symmetry-breaking features. These minima have marginal stability and will have long-range correlations between them. In the SK limit we have analytically studied the density of states ρ(λ) of the Hessian matrix in the annealed approximation. Despite the fact that in the presence of a random field there are no continuous symmetries, the spectrum extends down to zero with the usual sqrt[λ] form for the density of states for fields below the AT field. However, when the random field is larger than the AT field, there is a gap in the spectrum, which closes up as the AT field is approached. The VB model behaves differently and seems rather similar to studies of the three-dimensional Heisenberg spin glass in a random vector field.

Universal Non-Debye Scaling in the Density of States of Amorphous Solids

At the jamming transition, amorphous packings are known to display anomalous vibrational modes with a density of states (DOS) that remains constant at low frequency. The scaling of the DOS at higher packing fractions remains, however, unclear. One might expect to find a simple Debye scaling, but recent results from effective medium theory and the exact solution of mean-field models both predict an anomalous, non-Debye scaling. Being mean-field in nature, however, these solutions are only strictly valid in the limit of infinite spatial dimension, and it is unclear what value they have for finite-dimensional systems. Here, we study packings of soft spheres in dimensions 3 through 7 and find, away from jamming, a universal non-Debye scaling of the DOS that is consistent with the mean-field predictions. We also consider how the soft mode participation ratio evolves as dimension increases.

Statistics and Properties of Low-Frequency Vibrational Modes in Structural Glasses

Low-frequency vibrational modes play a central role in determining various basic properties of glasses, yet their statistical and mechanical properties are not fully understood. Using extensive numerical simulations of several model glasses in three dimensions, we show that in systems of linear size L sufficiently smaller than a crossover size L_{D}, the low-frequency tail of the density of states follows D(ω)∼ω^{4} up to the vicinity of the lowest Goldstone mode frequency. We find that the sample-to-sample statistics of the minimal vibrational frequency in systems of size L<L_{D} is Weibullian, with scaling exponents in excellent agreement with the ω^{4} law. We further show that the lowest-frequency modes are spatially quasilocalized and that their localization and associated quartic anharmonicity are largely frequency independent. The effect of preparation protocols on the low-frequency modes is elucidated, and a number of glassy length scales are briefly discussed.