Anomalous vibrational properties in the continuum limit of glasses

Enhanced collective vibrations in granular materials

Granular materials are defined as collections of macroscopic dissipative particles. Although these systems are ubiquitous in our lives, the nature and the causes of their non-trivial collective dynamics still remain elusive and have attracted significant interest in non-equilibrium physics. Here, we focus on the vibrational dynamics of granular materials. While the vibrational dynamics of random packings have been examined concerning the jamming transition, previous research has overlooked the role of contact dissipations. We conducted numerical and analytical investigations into the vibrational dynamics of random packings influenced by the normal dissipative force, which is the simplest model for contact dissipations. Our findings reveal that the kinetic energy per mode diverges in the low-frequency range, following the scaling law with the frequency ωl, indicating that low-frequency modes experience strong excitation and that the equipartition of energy is violated. Additionally, the spatial structure factor of the velocity field displays the scaling law Sv(q) ∝ q-2 with the wavenumber q, which signifies that the velocity field has an infinitely long range. We demonstrate that these phenomena arise from the effects of weaker damping on softer modes, where the particle displacements parallel to the contacts are minimal in the low-frequency modes, rendering normal dissipation ineffective at damping these modes.

Jamming transition and normal modes of polydispersed soft particle packing

The jamming transition of soft particles characterized by narrow size distributions has been well studied by physicists. However, polydispersed systems are more relevant to engineering, and the influence of polydispersity on jamming phenomena is still unexplored. Here, we numerically investigate jamming transitions of polydispersed soft particles in two dimensions. We find that polydispersity strongly influences contact forces, local coordination, and the jamming transition density. In contrast, the critical scaling of pressure and elastic moduli is not affected by the particle size distribution. Consistent with this observation, we find that the vibrational density of states is also insensitive to the polydispersity. Our results suggest that, regardless of particle size distributions, both mechanical and vibrational properties of soft particle packings near jamming are governed by the distance to jamming.

Structural fluctuations in active glasses

The glassy dynamics of dense active matter have recently become a topic of interest due to their importance in biological processes such as wound healing and tissue development. However, while the liquid-state properties of dense active matter have been studied in relation to the glass transition of active matter, the solid-state properties of active glasses have yet to be understood. In this work, we study the structural fluctuations in the active glasses composed of self-propelled particles. We develop a formalism to describe the solid-state properties of active glasses in the harmonic approximation limit and use it to analyze the displacement fields in the active glasses. Our findings reveal that the dynamics of high-frequency normal modes become quasi-static with respect to the active forces, and consequently, excitations of these modes are significantly suppressed. This leads to a violation of the equipartition law, suppression of particle displacements, and the apparent collective motion of active glasses. Overall, our results provide a fundamental understanding of the solid-state properties of active glasses.

Vibrational lifetimes and viscoelastic properties of ultrastable glasses

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

The nature of non-phononic excitations in disordered systems

The frequency scaling exponent of low-frequency excitations in microscopically small glasses, which do not allow for the existence of waves (phonons), has been in the focus of the recent literature. The density of states g(ω) of these modes obeys an ωs scaling, where the exponent s, ranging between 2 and 5, depends on the quenching protocol. The orgin of these findings remains controversal. Here we show, using heterogeneous-elasticity theory, that in a marginally-stable glass sample g(ω) follows a Debye-like scaling (s = 2), and the associated excitations (type-I) are of random-matrix type. Further, using a generalisation of the theory, we demonstrate that in more stable samples, other, (type-II) excitations prevail, which are non-irrotational oscillations, associated with local frozen-in stresses. The corresponding frequency scaling exponent s is governed by the statistics of small values of the stresses and, therefore, depends on the details of the interaction potential.

Instantaneous normal modes of glass-forming liquids during the athermal relaxation process of the steepest descent algorithm

Understanding glass formation by quenching remains a challenge in soft condensed matter physics. Recent numerical studies on steepest descent dynamics, which is one of the simplest models of quenching, revealed that quenched liquids undergo slow relaxation with a power law towards mechanical equilibrium and that the late stage of this process is governed by local rearrangements of particles. These advances motivate the detailed study of instantaneous normal modes during the relaxation process because the glassy dynamics is considered to be governed by stationary points of the potential energy landscape. Here, we performed a normal mode analysis of configurations during the steepest descent dynamics and found that the dynamics is driven by almost flat directions of the potential energy landscape at long times. These directions correspond to localized modes and we characterized them in terms of their statistics and structure using methods developed in the study of local minima of the potential energy landscape.

Development of broad modulus profile upon polymer-polymer interface formation between immiscible glassy-rubbery domains

Interfaces of glassy materials such as thin films, blends, and composites create strong unidirectional gradients to the local heterogeneous dynamics that can be used to elucidate the length scales and mechanisms associated with the dynamic heterogeneity of glasses. We focus on bilayer films of two different polymers with very different glass transition temperatures ([Formula: see text]) where previous work has demonstrated a long-range (∼200 nm) profile in local [Formula: see text] is established between immiscible glassy and rubbery polymer domains when the polymer-polymer interface is formed to equilibrium. Here, we demonstrate that an equally long-ranged gradient in local modulus [Formula: see text] is established when the polymer-polymer interface ([Formula: see text]5 nm) is formed between domains of glassy polystyrene (PS) and rubbery poly(butadiene) (PB), consistent with previous reports of a broad [Formula: see text] profile in this system. A continuum physics model for the shear wave propagation caused by a quartz crystal microbalance across a PB/PS bilayer film is used to measure the viscoelastic properties of the bilayer during the evolution of the PB/PS interface showing the development of a broad gradient in local modulus [Formula: see text] spanning [Formula: see text]180 nm between the glassy and rubbery domains of PS and PB. We suggest these broad profiles in [Formula: see text] and [Formula: see text] arise from a coupling of the spectrum of vibrational modes across the polymer-polymer interface as a result of acoustic impedance matching of sound waves with [Formula: see text] nm during interface broadening that can then trigger density fluctuations in the neighboring domain.

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.

Shear-induced criticality in glasses shares qualitative similarities with the Gardner phase

Although glass phases are ubiquitously found in various soft matter systems, we are still far from a complete understanding of them. The concept of marginal stability predicted by infinite-dimensional mean-field theories is drawing attention as a candidate for a universal and distinguishing unique feature of glasses. While among theoretical predictions, the non-Debye scaling has indeed been observed universally over various classes of glasses, the Gardner phase is found only in a limited portion of them. In this work, we numerically demonstrate that plastic events observed in two-dimensional Lennard-Jones glasses under quasistatic shear exhibit statistical properties that are qualitatively consistent with the picture of an infinitely hierarchical energy landscape associated with the Gardner phase.

Johari-Goldstein β relaxation in glassy dynamics originates from two-scale energy landscape

Supercooled liquids undergo complicated structural relaxation processes, which have been a long-standing problem in both experimental and theoretical aspects of condensed matter physics. In particular, past experiments widely observed for many types of molecular liquids that relaxation dynamics separated into two distinct processes at low temperatures. One of the possible interpretations is that this separation originates from the two-scale hierarchical topography of the potential energy landscape; however, it has never been verified. Molecular dynamics simulations are a promising approach to tackle this issue, but we must overcome laborious difficulties. First, we must handle a model of molecular liquids that is computationally demanding compared to simple spherical models, which have been intensively studied but show only a slower process: α relaxation. Second, we must reach a sufficiently low-temperature regime where the two processes become well-separated. Here, we handle an asymmetric dimer system that exhibits a faster process: Johari-Goldstein β relaxation. Then, we employ the parallel tempering method to access the low-temperature regime. These laborious efforts enable us to investigate the potential energy landscape in detail and unveil the first direct evidence of the topographic hierarchy that induces the β relaxation. We also successfully characterize the microscopic motions of particles during each relaxation process. Finally, we study the correlation between low-frequency modes and two relaxation processes. Our results establish a fundamental and comprehensive understanding of experimentally observed relaxation dynamics in supercooled liquids.

BOTAN: BOnd TArgeting Network for prediction of slow glassy dynamics by machine learning relative motion

Recent developments in machine learning have enabled accurate predictions of the dynamics of slow structural relaxation in glass-forming systems. However, existing machine learning models for these tasks are mostly designed such that they learn a single dynamic quantity and relate it to the structural features of glassy liquids. In this study, we propose a graph neural network model, "BOnd TArgeting Network," that learns relative motion between neighboring pairs of particles, in addition to the self-motion of particles. By relating the structural features to these two different dynamical variables, the model autonomously acquires the ability to discern how the self motion of particles undergoing slow relaxation is affected by different dynamical processes, strain fluctuations and particle rearrangements, and thus can predict with high precision how slow structural relaxation develops in space and time.

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

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

Low-frequency vibrational states in ideal glasses with random pinning

Glasses exhibit spatially localized vibrations in the low-frequency regime. These localized modes emerge below the boson peak frequency ω_{BP}, and their vibrational densities of state follow g(ω)∝ω^{4} (ω is frequency). Here, we attempt to address how the localized vibrations behave through the ideal glass transition. To do this, we employ a random pinning method, which enables us to study the thermodynamic glass transition. We find that the localized vibrations survive even in equilibrium glass states. Remarkably, the localized vibrations still maintain the properties of appearance below ω_{BP} and g(ω)∝ω^{4}. Our results provide important insight into the material properties of ideal glasses.

Density of states below the first sound mode in 3D glasses

Glasses feature universally low-frequency excess vibrational modes beyond Debye prediction, which could help rationalize, e.g., the glasses’ unusual temperature dependence of thermal properties compared to crystalline solids. The way the density of states of these low-frequency excess modes D( ω) depends on the frequency ω has been debated for decades. Recent simulation studies of 3D glasses suggest that D( ω) scales universally with ω4 in a low-frequency regime below the first sound mode. However, no simulation study has ever probed as low frequencies as possible to test directly whether this quartic law could work all the way to extremely low frequencies. Here, we calculated D( ω) below the first sound mode in 3D glasses over a wide range of frequencies. We find D( ω) scales with ω β with β < 4 at very low frequencies examined, while the ω4 law works only in a limited intermediate-frequency regime in some glasses. Moreover, our further analysis suggests our observation does not depend on glass models or glass stabilities examined. The ω4 law of D( ω) below the first sound mode is dominant in current simulation studies of 3D glasses, and our direct observation of the breakdown of the quartic law at very low frequencies thus leaves an open but important question that may attract more future numerical and theoretical studies.

Phonon transport properties of particulate physical gels

Particulate physical gels are sparse, low-density amorphous materials in which clusters of glasses are connected to form a heterogeneous network structure. This structure is characterized by two length scales, ξ s and ξ G: ξ s measures the length of heterogeneities in the network structure and ξ G is the size of glassy clusters. Accordingly, the vibrational states (eigenmodes) of such a material also exhibit a multiscale nature with two characteristic frequencies, [Formula: see text] and ω G, which are associated with ξ s and ξ G, respectively: (i) phonon-like vibrations in the homogeneous medium at [Formula: see text], (ii) phonon-like vibrations in the heterogeneous medium at [Formula: see text], and (iii) disordered vibrations in the glassy clusters at ω > ω G. Here, we demonstrate that the multiscale characteristics seen in the static structures and vibrational states also extend to the phonon transport properties. Phonon transport exhibits two distinct crossovers at frequencies ω* and ω G (or at wavenumbers of [Formula: see text] and [Formula: see text]). In particular, both transverse and longitudinal phonons cross over between Rayleigh scattering at [Formula: see text] and diffusive damping at [Formula: see text]. Remarkably, the Ioffe–Regel limit is located at the very low frequency of ω*. Thus, phonon transport is localized above ω*, even where phonon-like vibrational states persist. This markedly strong scattering behavior is caused by the sparse, porous structure of the gel.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

Low-Frequency Excess Vibrational Modes in Two-Dimensional Glasses

Glasses possess more low-frequency vibrational modes than predicted by Debye theory. These excess modes are crucial for the understanding of the low temperature thermal and mechanical properties of glasses, which differ from those of crystalline solids. Recent simulational studies suggest that the density of the excess modes scales with their frequency ω as ω^{4} in two and higher dimensions. Here, we present extensive numerical studies of two-dimensional model glass formers over a large range of glass stabilities. We find that the density of the excess modes follows D_{exc}(ω)∼ω^{2} up to around the boson peak, regardless of the glass stability. The stability dependence of the overall scale of D_{exc}(ω) correlates with the stability dependence of low-frequency sound attenuation. However, we also find that, in small systems, where the first sound mode is pushed to higher frequencies, at frequencies below the first sound mode, there are excess modes with a system size independent density of states that scales as ω^{3}.

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.

Novel elastic instability of amorphous solids in finite spatial dimensions

Recently, progress has been made in the understanding of anomalous vibrational excitations in amorphous solids. In the lowest-frequency region, the vibrational spectrum follows a non-Debye quartic law, which persists up to zero frequency without any frequency gap. This gapless vibrational density of states (vDOS) suggests that glasses are on the verge of instability. This feature of marginal stability is now highlighted as a key concept in the theories of glasses. In particular, the elasticity theory based on marginal stability predicts the gapless vDOS. However, this theory yields a quadratic law and not the quartic law. To address this inconsistency, we presented a new type of instability, which is different from the conventional one, and proposed that amorphous solids are marginally stable considering the new instability in the preceding study [M. Shimada, H. Mizuno and A. Ikeda, Soft Matter, 2020, 16, 7279]. In this study, we further extend and detail the results for these instabilities. By analyzing various examples of disorder, we demonstrate that real glasses in finite spatial dimensions can be marginally stable by the proposed novel instability.

Intermittent rearrangements accompanying thermal fluctuations distinguish glasses from crystals

Glasses exhibit vibrational and thermal properties that are markedly different from those of crystals. While recent works have advanced our understanding of vibrational excitations in glasses in the harmonic approximation limit, efforts in understanding finite-temperature anharmonic processes have been limited. In crystals, phonon-phonon coupling provides an extremely efficient mechanism for anharmonic decay that is also important in glasses. By using extensive molecular dynamics simulation of model atomic systems, here we first describe, both numerically and analytically, the anharmonic couplings in the crystal and the glass by focusing on the temperature dependence of the associated decay rates. Next, we show that an additional anharmonic channel of different origin emerges in the amorphous case, which induces unconventional intermittent rearrangements of particles. We have found that thermal vibrations in glasses trigger transitions among numerous different local minima of the energy landscape, which, however, are located within the same wide (meta)basin. These processes generate motions that are different from both diffusive and out-of-equilibrium aging dynamics. We suggest that (i) the observed intermittent rearrangements accompanying thermal fluctuations are crucial features distinguishing glasses from crystals and (ii) they can be considered as relics of the liquid state that survive the complete dynamic arrest taking place at the glass transition temperature.

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

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

Vibrational spectrum derived from local mechanical response in disordered solids

The low-frequency vibrations of glasses are markedly different from those of crystals. These vibrations have recently been categorized into two types: spatially extended vibrations, whose vibrational density of states (vDOS) follows a non-Debye quadratic law, and quasilocalized vibrations (QLVs), whose vDOS follows a quartic law. The former are explained by elasticity theory with quenched disorder and microscopic replica theory as being a consequence of elastic instability, but the origin of the latter is still debated. Here, we show that the latter can also be directly derived from elasticity theory with quenched disorder. We find another elastic instability that the theory encompasses but that has been overlooked so far, namely, the instability of the system against a local dipolar force. This instability gives rise to an additional contribution to the vDOS, and the spatial structure and energetics of the mode originating from this instability are consistent with those of the QLVs. Finally, we construct a model in which the additional contribution to the vDOS follows a quartic law.

Statistical mechanics of local force dipole responses in computer glasses

Soft quasilocalized modes (QLMs) are universally featured by structural glasses quenched from a melt, and are involved in several glassy anomalies such as the low-temperature scaling of their thermal conductivity and specific heat, and sound attenuation at intermediate frequencies. In computer glasses, QLMs may assume the form of harmonic vibrational modes under a narrow set of circumstances; however, direct access to their full distribution over frequency is hindered by hybridizations of QLMs with other low-frequency modes (e.g., phonons). Previous studies to overcome this issue have demonstrated that the response of a glass to local force dipoles serves as a good proxy for its QLMs; we, therefore, study here the statistical-mechanical properties of these responses in computer glasses, over a large range of glass stabilities and in various spatial dimensions, with the goal of revealing properties of the yet-inaccessible full distribution of QLMs' frequencies. We find that as opposed to the spatial-dimension-independent universal distribution of QLMs' frequencies ω (and, consequently, also of their stiffness κ = ω²), the distribution of stiffnesses associated with responses to local force dipoles features a (weak) dependence on spatial dimension. We rationalize this dependence by introducing a lattice model that incorporates both the real-space profiles of QLMs-associated with dimension-dependent long-range elastic fields-and the universal statistical properties of their frequencies. Based on our findings, we propose a conjecture about the form of the full distribution of QLMs' frequencies and its protocol-dependence. Finally, we discuss possible connections of our findings to basic aspects of glass formation and deformation.

Low-frequency vibrations of jammed packings in large spatial dimensions

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

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.

Low-frequency vibrational modes of stable glasses

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

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.