Boson Peak Decouples from Elasticity in Glasses with Low Connectivity

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.

Systematic Coarse-graining of Epoxy Resins with Machine Learning-Informed Energy Renormalization

A persistent challenge in predictive molecular modeling of thermoset polymers is to capture the effects of chemical composition and degree of crosslinking (DC) on dynamical and mechanical properties with high computational efficiency. We established a new coarse-graining (CG) approach that combines the energy renormalization method with Gaussian process surrogate models of the molecular dynamics simulations. This allows a machine-learning informed functional calibration of DC-dependent CG force field parameters. Taking versatile epoxy resins consisting of Bisphenol A diglycidyl ether combined with curing agent of either 4,4-Diaminodicyclohexylmethane or polyoxypropylene diamines, we demonstrated excellent agreement between all-atom and CG predictions for density, Debye-Waller factor, Young's modulus and yield stress at any DC. We further introduce a surrogate model enabled simplification of the functional forms of 14 non-bonded calibration parameters by quantifying the uncertainty of a candidate set of high-dimensional/flexible calibration functions. The framework established provides an efficient methodology for chemistry-specific, large-scale investigations of the dynamics and mechanics of epoxy resins.

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.

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.

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.

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.

Universal Origin of Boson Peak Vibrational Anomalies in Ordered Crystals and in Amorphous Materials

The vibrational spectra of solids, both ordered and amorphous, in the low-energy regime, control the thermal and transport properties of materials, from heat capacity to heat conduction, electron-phonon couplings, conventional superconductivity, etc. The old Debye model of vibrational spectra at low energy gives the vibrational density of states (VDOS) as proportional to the frequency squared, but in many materials the spectrum departs from this law which results in a peak upon normalizing the VDOS by frequency squared, which is known as the "boson peak." A description of the VDOS of solids (both crystals and glasses) is presented starting from first principles. Without using any assumptions whatsoever of disorder in the material, it is shown that the boson peak in the VDOS of both ordered crystals and glasses arises naturally from the competition between elastic mode propagation and diffusive damping. The theory explains the recent experimental observations of boson peak in perfectly ordered crystals, which cannot be explained based on previous theoretical frameworks. The theory also explains, for the first time, how the vibrational spectrum changes with the atomic density of the solid, and explains recent experimental observations of this effect.