Anharmonic damping of terahertz acoustic waves in a network glass and its effect on the density of vibrational states

Properties of stable ensembles of Euclidean random matrices

We study the spectrum of a system of coupled disordered harmonic oscillators in the thermodynamic limit. This Euclidean random matrix ensemble has been suggested as a model for the low temperature vibrational properties of glass. Exact numerical diagonalization is performed in three and two spatial dimensions, which is accompanied by a detailed finite size analysis. It reveals a low-frequency regime of sound waves that are damped by Rayleigh scattering. At large frequencies localized modes exist. In between, the central peak in the vibrational density of states is well described by Wigner's semicircle law for not too large disorder, as is expected for simple random matrix systems. We compare our results with predictions from two recent self-consistent field theories.

Vibrational Phenomena in Glasses at Low Temperatures Captured by Field Theory of Disordered Harmonic Oscillators

We investigate the vibrational properties of topologically disordered materials by analytically studying particles that harmonically oscillate around random positions. Exploiting classical field theory in the thermodynamic limit at T=0, we build up a self-consistent model by analyzing the Hessian utilizing Euclidean random matrix theory. In accordance with earlier findings [T. S. Grigera et al.J. Stat. Mech. (2011) P02015.JSMTC61742-546810.1088/1742-5468/2011/02/P02015], we take nonplanar diagrams into account to correctly address multiple local scattering events. By doing so, we end up with a first principles theory that can predict the main anomalies of athermal disordered materials, including the boson peak, sound softening, and Rayleigh damping of sound. In the vibrational density of states, the sound modes lead to Debye's law for small frequencies. Additionally, an excess appears in the density of states starting as ω^{4} in the low frequency limit, which is attributed to (quasi-) localized modes.

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.

Sound attenuation in two-dimensional glasses at finite temperatures

The thermal conductivity of glasses exhibits an unusual temperature dependence compared to their crystalline counterparts. Sound attenuation due to disorder in glasses was proposed to be important in rationalizing this special behavior. Simulation studies suggest that in the harmonic approximation, the sound attenuation follows Rayleigh scattering scaling at small wave vector in both two-dimensional (2D) and 3D glasses. The influence of the anharmonicity on sound attenuation has very recently been investigated numerically, but only in 3D glasses. Hence, it remains unknown in simulations how sound attenuation changes with the wave vector in 2D glasses when the anharmonicity comes into play. Here, we address this issue by performing computer simulations in low-temperature 2D glasses over a large range of glass stabilities. We find that the way the anharmonicity affects sound attenuation in 2D glasses is the same as that in 3D, thus revealing that numerically the influence of the anharmonicity on sound attenuation does not rely on the spatial dimension.

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.

Sound damping in soft particle packings: the interplay between configurational disorder and inelasticity

We numerically investigate sound damping in disordered two-dimensional soft particle packings. We simulate evolution of standing waves of particle displacements and analyze time correlation functions of particle velocities and power spectra. We control the strength of inelastic interactions between the particles in contact to show how the inelasticity affects anomalous sound characteristics of disordered systems: Increasing the strength of inelastic interactions, we find that (i) sound softening vanishes and (ii) attenuation coefficients exhibit a transition from the Rayleigh law to quadratic growth. We also report (iii) how the Ioffe-Regel limit frequencies depend on the strength of inelasticity as useful information for experiments and applications of the sound in disordered media. Our findings suggest that sound damping in soft particle packings is determined by the interplay between elastic heterogeneities and inelasticity.

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.

Analytical prediction of logarithmic Rayleigh scattering in amorphous solids from tensorial heterogeneous elasticity with power-law disorder

The damping or attenuation coefficient of sound waves in solids due to impurities scales with the wavevector to the fourth power, also known as Rayleigh scattering. In amorphous solids, Rayleigh scattering may be enhanced by a logarithmic factor although computer simulations offer conflicting conclusions regarding this enhancement and its microscopic origin. We present a tensorial replica field-theoretic derivation based on heterogeneous or fluctuating elasticity (HE), which shows that long-range (power-law) spatial correlations of the elastic constants, is the origin of the logarithmic enhancement to Rayleigh scattering of phonons in amorphous solids. We also consider the case of zero spatial fluctuations in the elastic constants, and of power-law decaying fluctuations in the internal stresses. Also in this case the logarithmic enhancement to the Rayleigh scattering law can be derived from the proposed tensorial HE framework.

Sound attenuation in finite-temperature stable glasses

The temperature dependence of the thermal conductivity of amorphous solids is markedly different from that of their crystalline counterparts, but exhibits universal behaviour. Sound attenuation is believed to be related to this universal behaviour. Recent computer simulations demonstrated that in the harmonic approximation sound attenuation Γ obeys quartic, Rayleigh scattering scaling for small wavevectors k and quadratic scaling for wavevectors above the Ioffe-Regel limit. However, simulations and experiments do not provide a clear picture of what to expect at finite temperatures where anharmonic effects become relevant. Here we study sound attenuation at finite temperatures for model glasses of various stability, from unstable glasses that exhibit rapid aging to glasses whose stability is equal to those created in laboratory experiments. We find several scaling laws depending on the temperature and stability of the glass. First, we find the large wavevector quadratic scaling to be unchanged at all temperatures. Second, we find that at small wavevectors Γ∼k^1.5 for an aging glass, but Γ∼k² when the glass does not age on the timescale of the calculation. For our most stable glass, we find that Γ∼k² at small wavevectors, then a crossover to Rayleigh scattering scaling Γ∼k⁴, followed by another crossover to the quadratic scaling at large wavevectors. Our computational observation of this quadratic behavior reconciles simulation, theory and experiment, and will advance the understanding of the temperature dependence of thermal conductivity of glasses.

Relaxation and vibrational properties in metal alloys and other disordered systems

The relaxation dynamics and the vibrational spectra of amorphous solids, such as metal alloys, have been intensely investigated as well separated topics in the past. The aim of this review is to summarize recent results in both these areas in an attempt to establish, or unveil, deeper connections between the two phenomena of relaxation and vibration. Theoretical progress in the area of slow relaxation dynamics of liquid and glassy systems and in the area of vibrational spectra of glasses and liquids is reviewed. After laying down a generic modelling framework to connect vibration and relaxation, the physics of metal alloys is considered where the emergence of power-law exponents has been identified both in the vibrational density of states (VDOS) as well as in density correlations. Also, theoretical frameworks which connect the VDOS to the relaxation behaviour and mechanical viscoelastic response in metallic glasses are reviewed. The same generic interpretative framework is then applied to the case of molecular glass formers where the emergence of stretched-exponential relaxation in dielectric relaxation can be put in quantitative relation with the VDOS by means of memory-function approaches. Further connections between relaxation and vibration are provided by the study of phonon linewidths in liquids and glasses, where a natural starting point is given by hydrodynamic theories. Finally, an agenda of outstanding issues including the appearance of compressed exponential relaxation in the intermediate scattering function of experimental and simulated systems (metal alloys, colloidal gels, jammed packings) is presented in light of available (or yet to be developed) mathematical models, and compared to non-exponential behaviour measured with macroscopic means such as mechanical spectroscopy/rheology.

Sound attenuation in stable glasses

Understanding the difference between the universal low-temperature properties of amorphous and crystalline solids requires an explanation for the stronger damping of long-wavelength phonons in amorphous solids. A longstanding sound attenuation scenario, resulting from a combination of experiments, theories, and simulations, leads to a quartic scaling of sound attenuation with the wavevector, which is commonly attributed to the Rayleigh scattering of sound. Modern computer simulations offer conflicting conclusions regarding the validity of this picture. We simulate glasses with an unprecedentedly broad range of stabilities to perform the first microscopic analysis of sound damping in model glass formers across a range of experimentally relevant preparation protocols. We present convincing evidence that quartic scaling is recovered for small wavevectors irrespective of the glass's stability. With increasing stability, the wavevector where the quartic scaling begins increases by approximately a factor of three and the sound attenuation decreases by over an order of magnitude. Our results uncover an intimate connection between glass stability and sound damping.

Fluctuating Elasticity Fails to Capture Anomalous Sound Scattering in Amorphous Solids

The fluctuating elasticity (FE) model, introduced phenomenologically and developed by Schirmacher [J. Non-Cryst. Solids 357, 518 (2011)JNCSBJ0022-309310.1016/j.jnoncrysol.2010.07.052], is today the only theoretical framework available to analyze low-temperature elastic acoustic scattering in glasses. Its existing formulations, which neglect the tensorial nature of elasticity and exclude long-range disorder correlations, predict that the acoustic damping coefficients obey the standard Rayleigh scaling law: Γ∼k^{d+1}, with k the acoustic wave vector, in dimension d. However, recent numerical data, supported by the analysis of existing experimental results, show that Γ does not obey this scaling law but Γ∼-k^{d+1}lnk. Here we analyze in detail how a fully tensorial FE model can be constructed as a long wavelength approximation of the elastic response of the discrete, atomistic, problem. We show that, although it incorporates all long-range correlations, it fails to capture the observed damping in two respects: (i) it misses the anomalous scaling, and predicts the standard Rayleigh law; (ii) it grossly underestimates the amplitude of scattering by about 2 orders of magnitude. This brings clear evidence that the small scale nonaffine displacement fields, although not simply reducible to local defects, play a crucial role in acoustic wave scattering and hence cannot be ignored.

Stochastic simulation of anharmonic dissipation. II. Harmonic bath potentials with quadratic couplings

The workhorse simulating the dissipative dynamics is mainly based on the harmonic bath potentials together with linear system-bath couplings, but a realistic bath always assumes anharmonicity. In this work, we extend the linear dissipation model to include quadratic couplings and suggest a stochastic simulation scheme for the anharmonic dissipation. We show that the non-Gaussian noises induced by the anharmonic bath can be rigorously constructed, and the resulting stochastic Liouville equation has the same form as that for the linear dissipation model. As a preliminary application, we use this stochastic method to investigate the vibration-induced symmetry breaking in two-level electronic systems and find that the characteristic function of the non-Gaussian noises determines the absorption and fluorescence spectra.

Rotational sound in disordered granular materials

We numerically study the evolution of elastic standing waves in disordered disk systems with a focus on the dispersion relations of rotational sound. As on a lattice, the rotational mode exhibits an optical-like dispersion relation in the high frequency regime, representing a shoulder in the vibrational density of states and fast oscillations of the autocorrelations of rotational velocities. If tangential stiffness between the disks is large enough, a lattice-based model perfectly describes the dispersion relation of the rotational mode. If it is comparable to or smaller than the normal stiffness, the model fails for short wavelengths. However, the dispersion relation then follows the model prediction for the transverse mode, implying that the fast oscillations of disks' rotations switch to acousticlike behavior. We evidence such a transition from rotational to transverse modes by analyzing their respective participation of different degrees of freedom to the eigenvectors.

Propagative and diffusive regimes of acoustic damping in bulk amorphous material

In amorphous solids, a non-negligible part of thermal conductivity results from phonon scattering on the structural disorder. The conversion of acoustic energy into thermal energy is often measured by the dynamical dtructure factor (DSF) thanks to inelastic neutron or x-ray scattering. The DSF is used to quantify the dispersion relation of phonons, together with their damping. However, the connection of the dynamical structure factor with dynamical attenuation of wave packets in glasses is still a matter of debate. We focus here on the analysis of wave-packet propagation in numerical models of amorphous silicon. We show that the damped harmonic oscillator model fits of the dynamical structure factors give a good estimate of the wave packets mean free path, only below the Ioffe-Regel frequency. Above the Ioffe-Regel frequency and below the mobility edge, a pure diffusive regime without a definite mean free path is observed. The high-frequency mobility edge is characteristic of a transition to localized vibrations. Below the Ioffe-Regel frequency, a mixed regime is evidenced at intermediate frequencies, with a coexistence of propagative and diffusive wave fronts. The transition between these different regimes is analyzed in detail and reveals a complex dynamics for energy transport, thus raising the question of the correct modeling of thermal transport in amorphous materials.

Damping of vibrational excitations in glasses at terahertz frequency: The case of 3-methylpentane

We report a compared analysis of inelastic X ray scattering (IXS) and of low frequency Raman data of glassy 3-methylpentane. The IXS spectra have been analysed allowing for the existence of two distinct excitations at each scattering wavevector obtaining a consistent interpretation of the spectra. In particular, this procedure allows us to interpret the linewidth of the modes in terms of a simple model which relates them to the width of the first sharp diffraction peak in the static structure factor. In this model, the width of the modes arises from the blurring of the dispersion curves which increases approaching the boundary of the first pseudo-Brillouin zone. The position of the boson peak contribution to the density of vibrational states derived from the Raman scattering measurements is in agreement with the interpretation of the two excitations in terms of a longitudinal mode and a transverse mode, the latter being a result of the mixed character of the transverse modes away from the center of the pseudo-Brillouin zone.

Anomalous phonon scattering and elastic correlations in amorphous solids

A major issue in materials science is why glasses present low-temperature thermal and vibrational properties that sharply differ from those of crystals. In particular, long-wavelength phonons are considerably more damped in glasses, yet it remains unclear how structural disorder at atomic scales affects such a macroscopic phenomenon. A plausible explanation is that phonons are scattered by local elastic heterogeneities that are essentially uncorrelated in space, a scenario known as Rayleigh scattering, which predicts that the damping of acoustic phonons scales with wavenumber k as k^d+1 (in dimension d). Here we demonstrate that phonon damping scales instead as - k^d+1 ln k, with this logarithmic enhancement originating from long-range spatial correlations of elastic disorder caused by similar stress correlations. Our work suggests that the presence of long-range spatial correlations of local stress and elasticity may well be the crucial feature that distinguishes amorphous solids from crystals.

THz Acoustic Spectroscopy by using Double Quantum Wells and Ultrafast Optical Spectroscopy

GaN is a pivotal material for acoustic transducers and acoustic spectroscopy in the THz regime, but its THz phonon properties have not been experimentally and comprehensively studied. In this report, we demonstrate how to use double quantum wells as a THz acoustic transducer for measuring generated acoustic phonons and deriving a broadband acoustic spectrum with continuous frequencies. We experimentally investigated the sub-THz frequency dependence of acoustic attenuation (i.e., phonon mean-free paths) in GaN, in addition to its physical origins such as anharmonic scattering, defect scattering, and boundary scattering. A new upper limit of attenuation caused by anharmonic scattering, which is lower than previously reported values, was obtained. Our results should be noteworthy for THz acoustic spectroscopy and for gaining a fundamental understanding of heat conduction.