Equivalence of the boson peak in glasses to the transverse acoustic van Hove singularity in crystals

Temperature dependence of spatial nanoheterogeneities of shear modulus in supercooled glycerol

The boson peak in the terahertz vibrational spectrum carries information about nano-heterogeneities in the shear modulus in glass formers. Its evolution upon heating or cooling in a supercooled liquid state may shed light on the temperature dependence of heterogeneities. For this purpose, an analysis of the light scattering spectra of supercooled glycerol in the spectral range of the boson peak and fast relaxation was carried out and the parameters of the boson peak in the temperature range 180-330 K were determined. The temperature dependent frequency of the boson peak was then expressed in terms of the mean-square amplitude of the shear modulus fluctuations. This was done using the heterogeneous elasticity theory in combination with the perturbation theory on small fluctuations and Ioffe-Regel criterion for transverse vibrations in glass formers. The contribution of structural relaxation effects to phonon damping becomes significant with increasing temperature. It is shown here that structural relaxation largely determines the temperature dependence of the mean-square fluctuations of the shear modulus at high temperatures. By solving the inverse problem, the temperature dependence of shear modulus fluctuations was obtained. It shows a rapid decrease above ∼250 K with a linear extrapolation going to zero at the so-called Arrhenius temperature TA = 350 K. Comparison with literature data on the Landau-Placzek ratio shows that they have a similar temperature dependence at T < TA, which is explained by the appearance of nanometer scale spatial heterogeneities below TA. This is confirmed by the temperature dependence of the amplitude of the boson peak.

Uncovering the phonon spectra and lattice dynamics of plastically deformable InSe van der Waals crystals

Stacking two-dimensional (2D) van der Waals (vdW) materials in a layered bulk structure provides an appealing platform for the emergence of exotic physical properties. As a vdW crystal with exceptional plasticity, InSe offers the opportunity to explore various effects arising from the coupling of its peculiar mechanical behaviors and other physical properties. Here, we employ neutron scattering techniques to investigate the correlations of plastic interlayer slip, lattice anharmonicity, and thermal transport in InSe crystals. Not only are the interlayer slip direction and magnitude well captured by shifts in the Bragg reflections, but we also observe a deviation from the expected Debye behaviour in the heat capacity and lattice thermal conductivity. Combining the experimental data with first-principles calculations, we tentatively attribute the observed evidence of strong phonon-phonon interactions to a combination of a large acoustic-optical frequency resonance and a nesting effect. These findings correlate the macroscopic plastic slip and the microscopic lattice dynamics, providing insights into the mechano-thermo coupling and modulation in 2D vdW materials.

Stringlet excitation model of the boson peak

The boson peak (BP), a low-energy excess in the vibrational density of states over the Debye contribution, is often identified as a characteristic of amorphous solid materials. Despite decades of efforts, its microscopic origin still remains a mystery. Recently, it has been proposed, and corroborated with simulations, that the BP might stem from intrinsic localized modes involving one-dimensional (1D) string-like excitations ("stringlets"). We build on a theory originally proposed by Lund that describes the localized modes as 1D vibrating strings, but we specify the stringlet size distribution to be exponential, as observed in simulations. We provide an analytical prediction for the BP frequency ωBP in the temperature regime well below the observed glass transition temperature Tg. The prediction involves no free parameters and accords quantitatively with prior simulation observations in 2D and 3D model glasses based on inverse power law potentials. The comparison of the string model to observations is more uncertain when compared to simulations of an Al-Sm metallic glass material at temperatures well above Tg. Nonetheless, our stringlet model of the BP naturally reproduces the softening of the BP frequency upon heating and offers an analytical explanation for the experimentally observed scaling with the shear modulus in the glass state and changes in this scaling in simulations of glass-forming liquids. Finally, the theoretical analysis highlights the existence of a strong damping for the stringlet modes above Tg, which leads to a large low-frequency contribution to the 3D vibrational density of states, observed in both experiments and simulations.

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.

Boson-peak vibrational modes in glasses feature hybridized phononic and quasilocalized excitations

A hallmark of structural glasses and other disordered solids is the emergence of excess low-frequency vibrations on top of the Debye spectrum DDebye(ω) of phonons (ω denotes the vibrational frequency), which exist in any solid whose Hamiltonian is translationally invariant. These excess vibrations-a signature of which is a THz peak in the reduced density of states D(ω)/DDebye(ω), known as the boson peak-have resisted a complete theoretical understanding for decades. Here, we provide direct numerical evidence that vibrations near the boson peak consist of hybridizations of phonons with many quasilocalized excitations; the latter have recently been shown to generically populate the low-frequency tail of the vibrational spectra of structural glasses quenched from a melt and of disordered crystals. Our results suggest that quasilocalized excitations exist up to and in the vicinity of the boson-peak frequency and, hence, constitute the fundamental building blocks of the excess vibrational modes in glasses.

Effect of the alkali vs iron ratio on glass transition temperature and vibrational properties of synthetic basalt-like glasses

Volcanic eruptions generate huge amounts of material with a wide range of compositions and therefore different physicochemical properties. We present a combined Raman and calorimetric study carried out on four synthetic basaltic glasses with different alkali vs iron ratio which spans the typical compositions of basalts on Earth. Differential scanning calorimetry shows that changes of this ratio modify the glass transition interval whereas Raman spectra allow to gain insight about the structure of the glass in the microscopic and macroscopic range. Indeed, our Raman analysis is extended from the high frequency region, characterized by the molecular peaks, to the very low frequency region where glasses exhibit the boson peak. Spectra show a variation of the non-bridging oxygens number that affects the medium range order of the glass and the network interconnections. In the considered substitution interval, the boson peak shape is conserved while its position shift upwards. This means that increasing the alkali vs iron content, the elastic medium hardens but it does not change nature. This study emphasizes the importance of considering the full-range spectra when analysing multicomponent or natural systems with small chemical variations.

Extreme ultraviolet transient gratings: A tool for nanoscale photoacoustics

Collective lattice dynamics determine essential aspects of condensed matter, such as elastic and thermal properties. These exhibit strong dependence on the length-scale, reflecting the marked wavevector dependence of lattice excitations. The extreme ultraviolet transient grating (EUV TG) approach has demonstrated the potential of accessing a wavevector range corresponding to the 10s of nm length-scale, representing a spatial scale of the highest relevance for fundamental physics and forefront technology, previously inaccessible by optical TG and other inelastic scattering methods. In this manuscript we report on the capabilities of this technique in the context of probing thermoelastic properties of matter, both in the bulk and at the surface, as well as discussing future developments and practical considerations.

Understanding the emergence of the boson peak in molecular glasses

A common feature of glasses is the "boson peak", observed as an excess in the heat capacity over the crystal or as an additional peak in the terahertz vibrational spectrum. The microscopic origins of this peak are not well understood; the emergence of locally ordered structures has been put forward as a possible candidate. Here, we show that depolarised Raman scattering in liquids consisting of highly symmetric molecules can be used to isolate the boson peak, allowing its detailed observation from the liquid into the glass. The boson peak in the vibrational spectrum matches the excess heat capacity. As the boson peak intensifies on cooling, wide-angle x-ray scattering shows the simultaneous appearance of a pre-peak due to molecular clusters consisting of circa 20 molecules. Atomistic molecular dynamics simulations indicate that these are caused by over-coordinated molecules. These findings represent an essential step toward our understanding of the physics of vitrification.

Upper bound of fragility from spatial fluctuations of shear modulus and boson peak in glasses

It is shown that the normalized rms fluctuation of the shear modulus on the medium-range order scale in glasses correlates with fragility: the higher fragility, the smaller the fluctuation amplitude. The latter is calculated within the heterogeneous elasticity theory using the data on the boson peak in glasses. On a smaller scale corresponding to cooperative structural relaxation, the normalized rms fluctuation of the infinite-frequency shear modulus was estimated using the data on the decoupling of viscosity and diffusion in supercooled liquids. These fluctuations are much smaller in amplitude, and, in contrast, they increase with increasing fragility. Extrapolation predicts intersection of both rms fluctuations and disappearing of the boson peak at the upper limit to fragility ≈180.

Role of Optical Phonons and Anharmonicity in the Appearance of the Heat Capacity Boson Peak-like Anomaly in Fully Ordered Molecular Crystals

We demonstrate that the heat capacity Boson peak (BP)-like anomaly appearing in fully ordered anharmonic molecular crystals emerges as a result of the strong interactions between propagating (acoustic) and low-energy quasi-localized (optical) phonons. In particular, we experimentally determine the low-temperature (<30 K) specific heat of the molecular crystal benzophenone and those of several of its fully ordered bromine derivatives. Subsequently, by means of theoretical first-principles methods based on density functional theory, we estimate the corresponding phonon dispersions and vibrational density of states. Our results reveal two possible mechanisms for the emergence of the BP-like anomaly: (i) acoustic-optic phonon avoided crossing, which gives rise to a pseudo-van Hove singularity in the acoustic phonon branches, and (ii) piling up of low-frequency optical phonons, which are quasi degenerate with longitudinal acoustic modes and lead to a surge in the vibrational density of states at low energies.

Investigation of the vibrational density of states of sodium carboxymethyl starch glass via terahertz time-domain spectroscopy

We investigated the vibrational density of states of sodium carboxymethyl starch (CM-starch) by terahertz (THz) time-domain spectroscopy. The CM-starch showed a broad peak at ∼3 THz. The structure of the peak was similar to those corresponding to glucose-based polymer glasses possessing hydrogen bonds. The boson peak (BP) appeared at 1.16 THz at the lowest temperature and disappeared because of the existence of excess wing at higher temperatures. However, based on our novel BP frequency determination method using the inflection point of the extinction coefficient, the BP frequency showed almost no dependence on temperature. Further, the chain length dependence of the BP frequency of the glucose-based glasses showed that the BP frequency of the polymer glass was slightly lower than that of the monomer glass. The power law behaviour of the absorption coefficient suggested the existence of fractons, and the fractal dimension was estimated to be 2.33.

Average density of states of amorphous Hamiltonians: role of phonon mediated coupling of nano-clusters

Based on a description of an amorphous solid as a collection of coupled nanosize molecular clusters referred as basic blocks, we analyse the statistical properties of its Hamiltonian. The information is then used to derive the ensemble averaged density of the vibrational states (non-phonon) which turns out to be a Gaussian in the bulk of the spectrum and an Airy function in the low frequency regime. A comparison with experimental data for six glasses confirms validity of our theoretical predictions.

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.

Unifying Description of the Vibrational Anomalies of Amorphous Materials

The vibrational density of states D(ω) of solids controls their thermal and transport properties. In crystals, the low-frequency modes are extended phonons distributed in frequency according to Debye's law, D(ω)∝ω^{2}. In amorphous solids, phonons are damped, and at low frequency D(ω) comprises extended modes in excess over Debye's prediction, leading to the so-called boson peak in D(ω)/ω^{2} at ω_{bp}, and quasilocalized ones. Here we show that boson peak and phonon attenuation in the Rayleigh scattering regime are related, as suggested by correlated fluctuating elasticity theory, and that amorphous materials can be described as homogeneous isotropic elastic media punctuated by quasilocalized modes acting as elastic heterogeneities. Our numerical results resolve the conflict between theoretical approaches attributing amorphous solids' vibrational anomalies to elastic disorder and localized defects.

Boson-peak-like anomaly caused by transverse phonon softening in strain glass

Strain glass is a glassy state with frozen ferroelastic/martensitic nanodomains in shape memory alloys, yet its nature remains unclear. Here, we report a glassy feature in strain glass that was thought to be only present in structural glasses. An abnormal hump is observed in strain glass around 10 K upon normalizing the specific heat by cubed temperature, similar to the boson peak in metallic glass. The simulation studies show that this boson-peak-like anomaly is caused by the phonon softening of the non-transforming matrix surrounding martensitic domains, which occurs in a transverse acoustic branch not associated with the martensitic transformation displacements. Therefore, this anomaly neither is a relic of van Hove singularity nor can be explained by other theories relying on structural disorder, while it verifies a recent theoretical model without any assumptions of disorder. This work might provide fresh insights in understanding the nature of glassy states and associated vibrational properties.

Heat capacity anomalies of the molecular crystal 1-fluoro-adamantane at low temperatures

Disorder-disorder phase transitions are rare in nature. Here, we present a comprehensive low-temperature experimental and theoretical study of the heat capacity and vibrational density of states of 1-fluoro-adamantane (C₁₀H₁₅F), an intriguing molecular crystal that presents a continuous disorder-disorder phase transition at T = 180 K and a low-temperature tetragonal phase that exhibits fractional fluorine occupancy. It is shown that fluorine occupancy disorder in the low-T phase of 1-fluoro-adamantane gives rise to the appearance of low-temperature glassy features in the corresponding specific heat (i.e., "boson peak" -BP-) and vibrational density of states. We identify the inflation of low-energy optical modes as the main responsible for the appearance of such glassy heat-capacity features and propose a straightforward correlation between the first localized optical mode and maximum BP temperature for disordered molecular crystals (either occupational or orientational). Thus, the present study provides new physical insights into the possible origins of the BP appearing in disordered materials and expands the set of molecular crystals in which "glassy-like" heat-capacity features have been observed.

Low temperature heat capacity of nanosize amorphous solids

Contrary to previous studies based on bulk materials, we analyze the behavior of the vibrational density of states (VDOS) and specific heat of an amorphous solid of nano-scales (∼3 nm) at low temperature. With Hamiltonian formulation based on the intermolecular dispersion forces, our analysis indicates a universal semi-circle form of the average VDOS in the bulk of the spectrum along with a super-exponentially increasing behavior in its edge. The latter in turn leads to a specific heat with a superlinear temperature (T) dependence belowT< 1Keven at nano-scales, and, surprisingly agreeing with the experiments although the latter are carried out at macroscopic scales. The omnipresence of dispersion forces at microscopic scales indicates the application of our results to other disordered materials too.

Understanding the scaling of boson peak through insensitivity of elastic heterogeneity to bending rigidity in polymer glasses

Amorphous materials exhibit peculiar mechanical and vibrational properties, including non-affine elastic responses and excess vibrational states, i.e., the so-called boson peak (BP). For polymer glasses, these properties are considered to be affected by the bending rigidity of the constituent polymer chains. In our recent work [Tomoshige,et al2019,Sci. Rep.919514], we have revealed simple relationships between the variations of vibrational properties and the global elastic properties: the response of the BP scales only with that of the global shear modulus. This observation suggests that the spatial heterogeneity of the local shear modulus distribution is insensitive to changes in the bending rigidity. Here, we demonstrate the insensitivity of elastic heterogeneity by directly measuring the local shear modulus distribution. We also study transverse sound wave propagation, which is also shown to scale only with the global shear modulus. Through these analyses, we conclude that the bending rigidity does not alter the spatial heterogeneity of the local shear modulus distribution, which yields vibrational and acoustic properties that are controlled solely by the global shear modulus of a polymer glass.

Open and Anisotropic Soft Regions in a Model Polymer Glass

The vibrational dynamics of a model polymer glass is studied by Molecular Dynamics simulations. The focus is on the “soft” monomers with high participation to the lower-frequency vibrational modes contributing to the thermodynamic anomalies of glasses. To better evidence their role, the threshold to qualify monomers as soft is made severe, allowing for the use of systems with limited size. A marked tendency of soft monomers to form quasi-local clusters involving up to 15 monomers is evidenced. Each chain contributes to a cluster up to about three monomers and a single cluster involves a monomer belonging to about 2–3 chains. Clusters with monomers belonging to a single chain are rare. The open and tenuous character of the clusters is revealed by their fractal dimension df<2. The inertia tensor of the soft clusters evidences their strong anisotropy in shape and remarkable linear correlation of the two largest eigenvalues. Owing to the limited size of the system, finite-size effects, as well as dependence of the results on the adopted polymer length, cannot be ruled out.

Low-Temperature Heat Capacity Anomalies in Ordered and Disordered Phases of Normal and Deuterated Thiophene

We measured the specific heat C_p of normal (C₄H₄S) and deuterated (C₄D₄S) thiophene in the temperature interval of 1 ≤ T, K ≤ 25. C₄H₄S exhibits a metastable phase II₂ and a stable phase V, both with frozen orientational disorder (OD), whereas C₄D₄S exhibits a metastable phase II₂, which is analogous to the OD phase II₂ of C₄H₄S and a fully ordered stable phase V. Our measurements demonstrate the existence of a large bump in the heat capacity of both stable and metastable C₄D₄S and C₄H₄S phases at temperatures of ∼10 K, which significantly departs from the expected Debye temperature behavior of C_p ≈ T³. This case study demonstrates that the identified low-temperature C_p anomaly, typically referred to as a "Boson-peak" in the context of glassy crystals, is not exclusive of disordered materials.

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.

Role of Anharmonic Interactions for Vibration Density of States in α-Cristobalite

The vibrational density of states (VDOS) of solids in the low-energy regime controls the thermal and transport properties of materials, such as heat capacity, heat conduction, free energy and entropy. In α-Cristobalite, the low-frequency part of vibration density of states (VDOS) has many common features with the Boson peak in silica glass of matched densities. Recent theoretical work reported that anharmonic phonon–phonon interactions were critical for the low-frequency part of VDOS in α-Cristobalite. Therefore, it is urgent to identify the role of different anharmonic interactions from first principles. In this paper, we focus on the main peak of the low-frequency part of VDOS in α-Cristobalite. Calculated by our own developed codes and first principles, we find that the quartic anharmonic interaction can increase the frequency of the peak, while the cubic anharmonic can reduce the frequency and change the shape of the peak. Meanwhile, the anharmonic interactions are critical for the temperature effect. Therefore, we calculated the temperature-dependent property of the peak. We find that the frequency of the peak is directly proportional to the temperature. The atomic displacement patterns of different temperatures also confirm the above conclusion. All our calculations converged well. Moreover, our basic results agree well with other published results. Finally, we highlight that our codes offer a general and reliable way to calculate the VDOS with temperature.

Physics of phonon-polaritons in amorphous materials

The nature of bosonic excitations in disordered materials has remained elusive due to the difficulties in defining key concepts such as quasi-particles in the presence of disorder. We report on an experimental observation of phonon-polaritons in glasses, including a prominent boson peak (BP), i.e., excess of THz modes over the Debye law. A theoretical framework based on the concept of diffusons is developed to describe the broadening linewidth of the polariton due to disorder-induced scattering. It is shown here for the first time that the BP frequency and the Ioffe-Regel (IR) crossover frequency of the polariton collapse onto the same power-law decay with the diffusivity of the bosonic excitation. This analysis dismisses the hypothesis of the BP being caused by a relic of the van Hove singularity. The presented framework establishes a new methodology to analyze bosonic excitations in amorphous media, well beyond the traditional case of acoustic phonons, and establishes the IR crossover as the fundamental physical mechanism behind the BP.

Low frequency vibrational density of state of highly permeable super glassy polynorbornenes - the Boson peak

Inelastic incoherent neutron time-of-flight scattering was employed to measure the low frequency density of states for a series of addition polynorbornenes with bulky side groups. The rigid main chain in combination with the bulky side groups give rise to a microporosity of these polymers in the solid state. The microporosity characterized by the BET surfaces area varies systematically in the considered series. Such materials have some possible application as active separation layer in gas separation membranes. All investigated materials show excess contributions to the Debye type density of states characteristic for glasses known as Boson peak. The maximum position of the Boson peak shifts to lower frequency values with increasing microporosity. Data for PIM-1 and Matrimid included for comparison are in good agreement to this dependency. This result supports the sound wave interpretation of the Boson peak.

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.

Detection of boson peak and fractal dynamics of disordered systems using terahertz spectroscopy

The boson peak is a largely unexplained excitation found universally in the terahertz vibrational spectra of disordered systems; the so-called fracton is a vibrational excitation associated with the self-similar structure of monomers in polymeric glasses. We demonstrate that such excitations can be detected using terahertz spectroscopy. In the case of fractal structures, we determine the infrared light-vibration coupling coefficient for the fracton region and show that information concerning the fractal and fracton dimensions appears in the exponent of the absorption coefficient. Finally, using terahertz time-domain spectroscopy and low-frequency Raman scattering, we experimentally observe these universal excitations in a protein (lysozyme) system that has an intrinsically disordered and fractal structure and argue that the system should be considered a single supramolecule. These findings are applicable to amorphous and fractal objects in general and will be valuable for understanding universal dynamics of disordered systems via terahertz light.

Probing the Nanostructure of Neutron-Irradiated Diamond Using Raman Spectroscopy

Disordering of crystal lattice induced by irradiation with fast neutrons and other high-energy particles is used for the deep modification of electrical and optical properties of diamonds via significant nanoscale restructuring and defects engineering. Raman spectroscopy was employed to investigate the nature of radiation damage below the critical graphitization level created when chemical vapor deposition and natural diamonds are irradiated by fast neutrons with fluencies from 1 × 10¹⁸ to 3 × 10²⁰ cm⁻² and annealed at the 100–1700 °C range. The significant changes in the diamond Raman spectra versus the neutron-irradiated conditions are associated with the formation of intrinsic irradiation-induced defects that do not completely destroy the crystalline feature but decrease the phonon coherence length as the neutron dose increases. It was shown that the Raman spectrum of radiation-damaged diamonds is determined by the phonon confinement effect and that the boson peak is present in the Raman spectra up to annealing at 800–1000 °C. Three groups of defect-induced bands (first group = 260, 495, and 730 cm⁻¹; second group = 230, 500, 530, 685, and 760 cm^–1; and third group = 335, 1390, 1415, and 1740 cm⁻¹) were observed in Raman spectra of fast-neutron-irradiated diamonds.

Observation of High-Frequency Transverse Phonons in Metallic Glasses

Using inelastic neutron scattering and molecular dynamics simulations on a model Zr-Cu-Al metallic glass, we show that transverse phonons persist well into the high-frequency regime, and can be detected at large momentum transfer. Furthermore, the apparent peak width of the transverse phonons was found to follow the static structure factor. The one-to-one correspondence, which was demonstrated for both Zr-Cu-Al metallic glass and a three-dimensional Lennard-Jones model glass, suggests a universal correlation between the phonon dynamics and the underlying disordered structure. This remarkable correlation, not found for longitudinal phonons, underscores the key role that transverse phonons hold for understanding the structure-dynamics relationship in disordered materials.

Molecular vibration and Boson peak analysis of glucose polymers and ester via terahertz spectroscopy

Complex permittivity spectra were obtained herein by performing broadband terahertz (THz) spectroscopy on cellulose, paramylon, and paramylon ester. Absorption peaks observed for cellulose and paramylon at approximately 3 THz are attributed to hydrogen bonds. In addition, a broad absorption peak around 2 THz was observed for all the polymers, demonstrating a general feature of polymer glasses derived from weak interatomic van der Waals forces. The boson peak was observed for cellulose and paramylon ester. The boson peak frequency for cellulose nearly equaled that for glassy glucose-a unit structure of the cellulose polymer. Additionally, the insensitivity of cellulose to the polymerization degree was consistent with recent results obtained via molecular dynamics simulations. In contrast, the boson peak frequency of paramylon ester was markedly smaller than that of cellulose. These results demonstrate the importance of hydrogen bonds as determinants of the boson peak frequency.

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.

Density scaling of structure and dynamics of an ionic liquid

The lines in the pressure–temperature phase diagram with constant conductivity are found to be lines where other dynamic variables as well as the molecular structure factor peak are constant, while charge ordering changes.

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

Structural and Thermal Properties in Formamidinium and Cs-Mixed Lead Halides

Formamidinium [FA, HC(NH₂)₂⁺] lead iodide and its cation mixture have attracted interest as potentials in applications for efficient solar cells superior to well-known methylammonium lead iodide. We investigated the crystal structure and thermodynamic properties of high-quality single crystals of FA_1-xCs_xPbI₃ for x = 0 and 0.1 through X-ray diffraction and heat capacity measurements. Both α-FA_0.9Cs_0.1PbI₃ as well as α-FAPbI₃ crystallize in a cubic Pm3̅m structure with orientationally disordered FA molecules confined in the nondistorted Pb-I framework. In FAPbI₃, we observed a second-order transition at 280 K and two first-order transitions at 141.2 and 130.2 K in between β- and γ-phases instead of the previously known single β-γ transition. After doping with 10% Cs, the multiple first-order transitions disappeared, leading to phase transitions emerging at 300 and 149 K with second-order character. We moreover observed low-energy localized modes for both compounds, which is presumably tied to anomalous thermal motion, rattling, of the FA molecule.

Frenkel line crossover of confined supercritical fluids

We investigate the temperature evolution of dynamics and structure of partially confined Lennard Jones (LJ) fluids in supercritical phase along an isobaric line in the P-T phase diagram using molecular dynamics simulations. We compare the Frenkel line (FL) crossover features of partially confined LJ fluids to that of the bulk LJ fluids in supercritical phase. Five different spacings have been chosen in this study and the FL crossover characteristics have been monitored for each of these spacings for temperatures ranging from 240 K to 1500 K keeping the pressure fixed at 5000 bar. We characterize the FL crossover using density of states (DoS) function and find that partially confined supercritical fluids (SCF) exhibit a progressive shift of FL crossover point to higher temperatures for smaller spacings. While the DoS perpendicular to the walls shows persistent oscillatory modes, the parallel component exhibits a smooth crossover from an oscillatory to non-oscillatory characteristics representative of FL crossover. We find that the vanishing of peaks in DoS parallel to the walls indicates that the SCF no longer supports shear mode excitations and could serve as an identifier of the FL crossover for confined systems just as is done for the bulk. Layer heights of density profiles, self-diffusivity and the peak heights of radial distribution function parallel to the walls also feature the FL crossover consistent with the DoS criteria. Surprisingly, self-diffusivity undergoes an Arrhenius to super-Arrhenius crossover at low temperatures for smaller spacings as a result of enhanced structural order evidenced via pair-excess entropy. This feature, typical of glass-forming liquids and binary supercooled liquids, is found to develop from the glass-like characteristic slowdown and strong caging in confined supercritical fluid, evidenced via mean squared displacement and velocity autocorrelation function respectively, over intermediate timescales.

Tracking the Connection between Disorder and Energy Landscape in Glasses Using Geologically Hyperaged Amber

Fossil amber offers the unique opportunity to investigate an amorphous material that has been exploring its energy landscape for more than 110 million years of natural aging. By applying different X-ray scattering methods to amber before and after annealing the sample to erase its thermal history, we identify a link between the potential energy landscape and the structural and vibrational properties of glasses. We find that hyperaging induces a depletion of the vibrational density of states in the terahertz region, also ruling the sound dispersion and attenuation properties of the corresponding acoustic waves. Critically, this is accompanied by a densification with structural implications different in nature from that caused by hydrostatic compression. Our results, rationalized within the framework of fluctuating elasticity theory, reveal how upon approaching the bottom of the potential energy landscape (9% decrease in the fictive temperature) the elastic matrix becomes increasingly less disordered (6%) and longer-range correlated (22%).

Structural Parameter of Orientational Order to Predict the Boson Vibrational Anomaly in Glasses

It has so far remained a major challenge to quantitatively predict the boson peak, a THz vibrational anomaly universal for glasses, from features in the amorphous structure. Using molecular dynamics simulations of a model Cu_{50}Zr_{50} glass, we decompose the boson peak to contributions from atoms residing in different types of Voronoi polyhedra. We then introduce a microscopic structural parameter to depict the "orientational order," using the vector pointing from the center atom to the farthest vertex of its Voronoi coordination polyhedron. This order parameter represents the most probable direction of transverse vibration at low frequencies. Its magnitude scales linearly with the boson peak intensity, and its spatial distribution accounts for the quasilocalized modes. This correlation is shown to be universal for different types of glasses.

Boson Peak Decouples from Elasticity in Glasses with Low Connectivity

We perform molecular-dynamics simulations of the vibrational and elastoplastic properties of polymeric glasses and crystals and the corresponding atomic systems. We evidence that the elastic scaling of the density of states in the low-frequency boson peak (BP) region is different in crystals and glasses. Also, we see that the BP of the polymeric glass is nearly coincident with the one of the atomic glasses, thus revealing that the former-unlike the elasticity-is controlled by nonbonding interactions only. Our results suggest that the interpretation of the BP in terms of the macroscopic elasticity, discussed in highly connected systems, does not hold for systems with low connectivity.

SciPhon: a data analysis software for nuclear resonant inelastic X-ray scattering with applications to Fe, Kr, Sn, Eu and Dy

The synchrotron radiation technique of nuclear resonant inelastic X-ray scattering (NRIXS), also known as nuclear resonance vibrational spectroscopy or nuclear inelastic scattering, provides a wealth of information on the vibrational properties of solids. It has found applications in studies of lattice dynamics and elasticity, superconductivity, heme biochemistry, seismology, isotope geochemistry and many other fields. It involves probing the vibrational modes of solids by using the nuclear resonance of Mössbauer isotopes such as ⁵⁷Fe, ⁸³Kr, ¹¹⁹Sn, ¹⁵¹Eu and ¹⁶¹Dy. After data reduction, it provides the partial phonon density of states of the Mössbauer isotope that is investigated, as well as many other derived quantities such as the mean force constant of the chemical bonds and the Debye velocity. The data reduction is, however, not straightforward and involves removal of the elastic peak, normalization and Fourier-Log transformation. Furthermore, some of the quantities derived are highly sensitive to details in the baseline correction. A software package and several novel procedures to streamline and hopefully improve the reduction of the NRIXS data generated at sector 3ID of the Advanced Photon Source have been developed. The graphical user interface software is named SciPhon and runs as a Mathematica package. It is easily portable to other platforms and can be easily adapted for reducing data generated at other beamlines. Several tests and comparisons are presented that demonstrate the usefulness of this software, whose results have already been used in several publications. Here, the SciPhon software is used to reduce Kr, Sn, Eu and Dy NRIXS data, and potential implications for interpreting natural isotopic variations in those systems are discussed.

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.

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.

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.

Boson peak, heterogeneity and intermediate-range order in binary SiO2-Al2O3 glasses

In binary aluminosilicate liquids and glasses, heterogeneity on intermediate length scale is a crucial factor for optical fiber performance, determining the lower limit of optical attenuation and Rayleigh scattering, but also clustering and precipitation of optically active dopants, for example, in the fabrication of high-power laser gain media. Here, we consider the low-frequency vibrational modes of such materials for assessing structural heterogeneity on molecular scale. We determine the vibrational density of states VDoS g(ω) using low-temperature heat capacity data. From correlation with low-frequency Raman spectroscopy, we obtain the Raman coupling coefficient. Both experiments allow for the extraction of the average dynamic correlation length as a function of alumina content. We find that this value decreases from about 3.9 nm to 3.3 nm when mildly increasing the alumina content from zero (vitreous silica) to 7 mol%. At the same time, the average inter-particle distance increases slightly due to the presence of oxygen tricluster species. In accordance with Loewensteinian dynamics, this proves that mild alumina doping increases structural homogeneity on molecular scale.

Boson peak dynamics of natural polymer starch investigated by terahertz time-domain spectroscopy and low-frequency Raman scattering

Terahertz time-domain spectroscopy and low-frequency Raman scattering were performed on the natural polymer starch to investigate the boson peak (BP) dynamics. In the infrared spectrum, the BP was observed at 0.99THz at the lowest temperature. Compared to the result from a previous study for vitreous glucose, both the frequency of the BP and absorption coefficient show lower values than those of the vitreous glucose. These behaviors originate from the longer correlation length of the medium-range order and lower concentration of hydroxyl groups in the starch. In the Raman spectrum, the BP was observed at 1.1THz at room temperature, although the BP was not observed around room temperature due to the excess wing of the fast relaxation modes in the infrared spectrum. The temperature dependence of ε″(ν) during the heating process and cooling process shows a hysteresis below 230K. During the heating process, kinks were observed at 140K and 230K. These kinks are attributed to the β-relaxation and the β_wet-relaxation, respectively.

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.

Evidence of a one-dimensional thermodynamic phase diagram for simple glass-formers

Glass formers show motional processes over an extremely broad range of timescales, covering more than ten orders of magnitude, meaning that a full understanding of the glass transition needs to comprise this tremendous range in timescales. Here we report simultaneous dielectric and neutron spectroscopy investigations of three glass-forming liquids, probing in a single experiment the full range of dynamics. For two van der Waals liquids, we locate in the pressure-temperature phase diagram lines of identical dynamics of the molecules on both second and picosecond timescales. This confirms predictions of the isomorph theory and effectively reduces the phase diagram from two to one dimension. The implication is that dynamics on widely different timescales are governed by the same underlying mechanisms.

Anomalies in the low frequency vibrational density of states for a polymer with intrinsic microporosity – the Boson peak of PIM-1

For the first time the low frequency vibrational density of states is reported for a polymer with intrinsic microporosity.