Theory of low-energy Raman scattering in glasses

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

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

Microscopic analysis of sound attenuation in low-temperature amorphous solids reveals quantitative importance of non-affine effects

Sound attenuation in low-temperature amorphous solids originates from their disordered structure. However, its detailed mechanism is still being debated. Here, we analyze sound attenuation starting directly from the microscopic equations of motion. We derive an exact expression for the zero-temperature sound damping coefficient. We verify that the sound damping coefficients calculated from our expression agree very well with results from independent simulations of sound attenuation. Small wavevector analysis of our expression shows that sound attenuation is primarily determined by the non-affine displacements’ contribution to the sound wave propagation coefficient coming from the frequency shell of the sound wave. Our expression involves only quantities that pertain to solids’ static configurations. It can be used to evaluate the low-temperature sound damping coefficients without directly simulating sound attenuation.

Probing Small-Angle Molecular Motions with EPR Spectroscopy: Dynamical Transition and Molecular Packing in Disordered Solids

Disordered molecular solids present a rather broad class of substances of different origin—amorphous polymers, materials for photonics and optoelectronics, amorphous pharmaceutics, simple molecular glass formers, and others. Frozen biological media in many respects also may be referred to this class. Theoretical description of dynamics and structure of disordered solids still does not exist, and only some phenomenological models can be developed to explain results of particular experiments. Among different experimental approaches, electron paramagnetic resonance (EPR) applied to spin probes and labels also can deliver useful information. EPR allows probing small-angle orientational molecular motions (molecular librations), which intrinsically are inherent to all molecular solids. EPR is employed in its conventional continuous wave (CW) and pulsed—electron spin echo (ESE)—versions. CW EPR spectra are sensitive to dynamical librations of molecules while ESE probes stochastic molecular librations. In this review, different manifestations of small-angle motions in EPR of spin probes and labels are discussed. It is shown that CW-EPR-detected dynamical librations provide information on dynamical transition in these media, similar to that explored with neutron scattering, and ESE-detected stochastic librations allow elucidating some features of nanoscale molecular packing. The possible EPR applications are analyzed for gel-phase lipid bilayers, for biological membranes interacting with proteins, peptides and cryoprotectants, for supercooled ionic liquids (ILs) and supercooled deep eutectic solvents (DESs), for globular proteins and intrinsically disordered proteins (IDPs), and for some other molecular solids.

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.

Phonon, plasmon and electronic properties of surfaces and interfaces of periodic W/Si and Si/W multilayers

The phonon and plasmon excitations and electronic properties of interfaces of periodic W/Si and Si/W multilayer structures were investigated. The Boson band originated from quasilocal surface acoustic phonons for ultrathin Si layers, excited by Raman scattering. In confined Si layers, a small fraction of crystalline Si nanoclusters were embedded within a large volume fraction of amorphous Si (a-Si) nanoclusters. The size of the a-Si nanoclusters was smaller for the thinner Si layer in the periodic layers. The plasmon energy in the Si layer was blueshifted with a decrease in the thickness of this layer. This was explained by the size-dependent quantization of plasmon shift. The valence band spectra comprised a substantial fine structure, which is associated with the interaction of valence orbitals of the W and Si atoms at the interface boundaries. For thinner Si layers, the binding interaction of W5d and Si3p states leads to the splitting of the density of states near the Fermi level in the energy range of 1.5-5 eV. However, the energy splitting with two maxima was observed at 0.7 and 2.4 eV for thicker layers. Thus, the results of X-ray photoelectron spectroscopy have indicated that the interface of W/Si multilayers consists of metal-enriched tungsten silicide. Both the atomic structure and the elemental composition of the silicide were modified with a change in the thickness of the Si layers. This novel investigation could be essential for designing nanomirrors with higher reflectivity.

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.

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.

Aging effects on thermal conductivity of glass-forming liquids

Thermal conductivity of a model glass-forming system in the liquid and glass states is studied using extensive numerical simulations. We show that near the glass transition temperature, where the structural relaxation time becomes very long, the measured thermal conductivity decreases with increasing age. Second, the thermal conductivity of the disordered solid obtained at low temperatures is found to depend on the cooling rate with which it was prepared. For the cooling rates accessible in simulations, lower cooling rates lead to lower thermal conductivity. Our analysis links this decrease of the thermal conductivity with increased exploration of lower-energy inherent structures of the underlying potential energy landscape. Further, we show that the lowering of conductivity for lower-energy inherent structures is related to the high-frequency harmonic modes associated with the inherent structure being less extended. Possible effects of considering relatively small systems and fast cooling rates in the simulations are discussed.

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Boson peak and Ioffe-Regel criterion in amorphous siliconlike materials: The effect of bond directionality

The vibrational properties of model amorphous materials are studied by combining complete analysis of the vibration modes, dynamical structure factor, and energy diffusivity with exact diagonalization of the dynamical matrix and the kernel polynomial method, which allows a study of very large system sizes. Different materials are studied that differ only by the bending rigidity of the interactions in a Stillinger-Weber modelization used to describe amorphous silicon. The local bending rigidity can thus be used as a control parameter, to tune the sound velocity together with local bonds directionality. It is shown that for all the systems studied, the upper limit of the Boson peak corresponds to the Ioffe-Regel criterion for transverse waves, as well as to a minimum of the diffusivity. The Boson peak is followed by a diffusivity's increase supported by longitudinal phonons. The Ioffe-Regel criterion for transverse waves corresponds to a common characteristic mean-free path of 5-7 Å (which is slightly bigger for longitudinal phonons), while the fine structure of the vibrational density of states is shown to be sensitive to the local bending rigidity.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Effect of glassy modes on electron spin-lattice relaxation in solid ethanol

Electron spin-lattice relaxation (SLR) of TEMPO radical was measured in the crystalline and glassy states of deuterated ethanol in the temperature range 5-80K using X-band electron paramagnetic resonance (EPR). The measured SLR rates are higher in the glassy than in crystalline state and the excess SLR rate in glassy state is much lower than in ethanol. This result suggests that extra modes in glassy state, i.e. glassy modes, produce the excess SLR rate via the electron-nuclear dipolar (END) interaction between the electron spin of radical and the matrix protons or deuterons. Using the soft-potential model and assuming the END interaction between the electron spin and the matrix protons, the contributions to SLR rate of various mechanisms of glassy modes were theoretically analyzed. The evaluations of SLR rates in glassy ethanol indicate two main mechanisms of glassy modes: thermally activated relaxation of double-well systems and phonon-induced relaxation of quasi-harmonic local modes. The SLR rates induced by these mechanisms correlate well with the experimental data.

Acoustic dissipation and density of states in liquid, supercooled, and glassy glycerol

Combined Brillouin spectra collected at visible, ultraviolet, and x-ray frequencies are used to reconstruct the imaginary part of the acoustic compliance J'' over a wide frequency range between 0.5 GHz and 5 THz. For liquid, supercooled, and glassy glycerol, J'' is found to be linearly dependent on the tagged-particle susceptibility measured by incoherent neutron scattering up to ≃1  THz, giving evidence of a clear relation between acoustic power dissipation and density of states. A simple but general formalism is presented to quantitatively explain this relation, thus clarifying the connection between the quasielastic component observed in neutron scattering experiments and the fast relaxation dynamics probed by Brillouin scattering.

The role of nucleation in vitrification of supercooled liquids

Low-frequency Raman and differential scanning calorimetry (DSC) investigations were carried out during a structural transformation of supercooled liquid salol (phenyl salicylate) in a wide temperature range. DSC experiments indicate that in the supercooled liquid salol at temperature ~40 K above the glass transition temperature metastable nuclei start to form. During subsequent cooling the nuclei become an important element of the glass structure, and thereby are considered as a measure of the intermediate range order in this glass. It was shown that the crystalline structure of the metastable nuclei differ from that of the stable nuclei. Low-frequency Raman spectra of the glassy salol show a broad band in the spectral range from 14.5 to 17.2 cm(-1); the so called 'Boson peak', which can be interpreted in terms of its relationship to the formation of structured clusters, with typical sizes in the nanometer range (critical radii).

Influence of pressure on quasielastic scattering in glasses: relationship to the boson peak

We report unexpectedly strong variations in the quasielastic scattering (QES) intensity in glasses under pressure. Analysis of the data reveals strong correlations between pressure-induced changes in the QES intensity and the intensity of the boson peak. This observation emphasizes a direct relationship between these two components of the fast dynamics. In addition, we observe changes of the QES spectral shape that can be interpreted as pressure-induced variations in the underlying energy landscape.

Raman scattering and the low-frequency vibrational spectrum of glasses

We present a theory of low-frequency Raman scattering in glasses, based on the concept that light couples to the elastic strains via spatially fluctuating elasto-optic (Pockels) constants. We show that the Raman intensity is not proportional to the vibrational density of states (as was widely believed), but to a convolution of Pockels constant correlation functions with the dynamic strain susceptibilities of the glass. Using the dynamic susceptibilities of a system with fluctuating elastic constants we are able for the first time to describe the Raman intensity and the anomalous vibration spectrum of a glass on the same footing. Good agreement between the theory and experiment for the Raman spectrum, the density of states, and the specific heat is demonstrated at the example of glassy As(2)S(3).

Acoustic attenuation in glasses and its relation with the boson peak

A theory for the vibrational dynamics in disordered solids [W. Schirmacher, Europhys. Lett. 73, 892 (2006), based on the random spatial variation of the shear modulus, has been applied to determine the wave vector (k) dependence of the Brillouin peak position (Omega(k)) and width (Gamma(k)), as well as the density of vibrational states [g(omega)], in disordered systems. As a result, we give a firm theoretical ground to the ubiquitous k2 dependence of Gamma(k) observed in glasses. Moreover, we derive a quantitative relation between the excess of the density of states (the boson peak) and Gamma(k), two quantities that were not considered related before. The successful comparison of this relation with the outcome of experiments and numerical simulations gives further support to the theory.