Search

Phonon quasi-particles have been monumental in microscopically understanding thermodynamics and transport properties in condensed matter for decades. Phonons have one-to-one correspondence with harmonic eigenstates and their energies are often described by simple independent harmonic oscillator models. Higher order terms in the potential energy lead to interactions among them, resulting in finite lifetimes and frequency shifts, even in perfect crystals. However, increasing evidence including constant volume heat capacity different from the expected Dulong-Petit law suggests the need for re-evaluation of phonons having harmonic energies. In this work, we explicitly examine inter-mode dependence of phonon energies of a prototypical crystal, silicon, through energy covariance calculations and demonstrate the concerted nature of phonon energies even at 300 K, questioning independent harmonic oscillator assumptions commonly used for phonon energy descriptions of thermodynamics and transport., Comment: Same content as the previous version but in a letter form

Thermal transport properties of amorphous carbon has attracted increasing attention due to its extreme thermal properties: It has been reported to have among the highest thermal conductivity for bulk amorphous solids up to $\sim$ 37 Wm\textsuperscript{-1}K\textsuperscript{-1}, comparable to crystalline sapphire ($\alpha$-Al\textsubscript{2}O\textsubscript{3}). Further, large density dependence in thermal conductivity demonstrates a potential for largely tunable thermal conductivity. However, mechanism behind the high thermal conductivity and its large density dependence remains elusive due to many variables at play. In this work, we perform large-scale ($\sim$ 10\textsuperscript{5} atoms) molecular dynamics simulations utilizing a machine learning potential based on neural networks. Through spectral decomposition of thermal conductivity which enables a quantum correction to classical heat capacity, we find that propagating vibrational excitations govern thermal transport in amorphous carbon ($\sim$ 100 \% of thermal conductivity) in sharp contrast to the conventional wisdom that diffusive vibrational excitations dominate thermal transport in amorphous solids. Instead, this remarkable behavior resembles thermal transport in simple crystals. Moreover, our temperature dependent spectral diffusivity and velocity current correlation analyses reveal that the density dependent thermal conductivity originates from anharmonicity sensitive propagating excitations. Our work suggests a novel insight and design principle into developing mechanically hard, thermally conductive amorphous solids.

Decomposition of atomic motion into individual normal modes has led to remarkable success in microscopically understanding thermal properties and thermodynamics in simple solids. We start this chapter with an example of decomposing atomic motion of a simple monatomic linear chain crystal into normal modes followed by a more general, classical normal mode formalism. Different classifications of normal modes such as phonons, propagons, diffusons, and locons are introduced. Finally, heat capacity and thermal conductivity predictions from the normal mode formalism are demonstrated.

Developing microscopic understanding of the thermal properties of liquids is challenging due to their strong dynamic disorder, which prevents characterization of the atomic degrees of freedom. There have been significant research interests in the past few decades to extend the normal mode analysis for solids to instantaneous structures of liquids. However, the nature of normal modes that arise from these unstable structures is still elusive. In this work, we explore the instantaneous eigenmodes of dynamical matrices of various Lennard-Jones argon liquid/gas systems at high temperatures and show that the normal modes can be interpreted as an interpolation of T \to \infty (gas) and T = 0 (solid) mode descriptions. We find that normal modes become increasingly collisional and translational, recovering atomistic gas-like behavior rather than vibrational with increase in temperature, suggesting that normal modes in liquids may be described by both solid-like and gas-like modes.

Understanding thermodynamics in liquids at the atomic level is challenging because of strong atomic interactions and lack of symmetry. Recent prior theoretical works have focused on describing heat capacity of liquids in terms of phonon-like excitations but often rely on fitting parameters and ad hoc assumptions. In this work, we perform microscopic analysis on instantaneous normal modes and velocity autocorrelations on molecular dynamics simulations of single element systems over wide ranges of temperature (up to $10^8$ K) and pressure (up to 1 TPa). Our results demonstrate that heat capacity of liquids can be described by a combination of both solid-like and gas-like degrees of freedom, leading to a unified framework to describe heat capacity of all three phases of matter: solid, liquid, and gas.

Metallic glasses exhibit fast mechanical relaxations at temperatures well below the glass transition, one of which shows little variation with temperature known as nearly constant loss (NCL). Despite the important implications of this phenomenon to in aging and deformation, the origin of the relaxation is unclear. Through molecular dynamics simulations of a model metallic glass, Cu_64.5Zr_35.5, we implement dynamic mechanical analysis with system stress decomposed into atomic-level stresses to identify the group of atoms responsible for NCL. This work demonstrates that NCL relaxation is due to fully transient groups of atoms that become normal over picosecond timescales. They are spatially distributed throughout the glass and have no outstanding features, rather than defect-like as previously reported.

The attenuation rate of vibrational excitations in various metallic liquids and glasses has been reported to change from the quadratic dependence on wavevector at low wavevectors to the linear dependence at high wavevectors. However, the origin of this behavior is not clear. Here, the analysis of this phenomenon through molecular dynamics is presented for prototypical metallic liquids, Cu56Zr44 and Fe. It is shown that the crossover wavevector is strongly correlated with the structural coherence length characterizing coarse-grained density correlations. We suggest that the linear dependence is caused by scattering of vibrational excitations by structural activation processes with low activation energies which are distinctively observed in metallic systems.

Normal mode decomposition of atomic vibrations has been used to provide microscopic under-standing of thermal transport in amorphous solids for decades. In normal mode methods, it is naturally assumed that atoms vibrate around their equilibrium positions and that individual normal modes are the fundamental vibrational excitations transporting heat. With the abundance of predictions from normal mode methods and experimental measurements now available, we care-fully analyze these calculations in amorphous silicon, a model amorphous solid. We find a number of discrepancies, suggesting that treating individual normal modes as fundamental heat carriers may not be accurate in amorphous solids. Further, our classical and ab-initio molecular dynamics simulations of amorphous silicon demonstrate a large degree of atomic diffusion, especially at high temperatures, leading to the conclusion that thermal transport in amorphous solids could be better described starting from the perspective of liquid dynamics rather than from crystalline solids

The heat-carrying acoustic excitations of amorphous silicon are of interest because their mean free paths may approach micron scales at room temperature. Despite extensive investigation, the origin of the weak acoustic damping in the heat-carrying frequencies remains a topic of debate. Here, we report measurements of the thermal conductivity mean free path accumulation function in amorphous silicon thin films from 60 - 315 K using transient grating spectroscopy. With additional picosecond acoustics measurements and considering the known frequency-dependencies of damping mechanisms in glasses, we reconstruct the mean free paths from $\sim 0.1-3$ THz. The mean free paths are independent of temperature and exhibit a Rayleigh scattering trend over most of this frequency range. The observed trend is inconsistent with the predictions of numerical studies based on normal mode analysis but agrees with diverse measurements on other glasses. The micron-scale MFPs in amorphous Si arise from the absence of anharmonic or two-level system damping in the sub-THz frequencies, leading to heat-carrying acoustic excitations with room-temperature damping comparable to that of other glasses at cryogenic temperatures.

Thermal atomic vibrations in amorphous solids can be distinguished by whether they propagate as elastic waves or do not propagate due to lack of atomic periodicity. In a-Si, prior works concluded that non-propagating waves are the dominant contributors to heat transport, while propagating waves are restricted to frequencies less than a few THz and are scattered by anharmonicity. Here, we present a lattice and molecular dynamics analysis of vibrations in a-Si that supports a qualitatively different picture in which propagating elastic waves dominate the thermal conduction and are scattered by elastic fluctuations rather than anharmonicity. We explicitly demonstrate the propagating nature of vibration with frequency approaching 10 THz using a triggered wave computational experiment. Our work suggests that most heat is carried by propagating elastic waves in a-Si and demonstrates a route to achieve extreme thermal properties in amorphous materials by manipulating elastic fluctuations.

Pure amorphous solids are traditionally considered to set the lower bound of thermal conductivity due to their disordered atomic structure that impedes vibrational energy transport. However, the lower limits for thermal conductivity in heterogeneous amorphous solids and the physical mechanisms underlying these limits remain unclear. Here, we use equilibrium molecular dynamics to show that an amorphous SiGe nanocomposite can possess thermal conductivity substantially lower than those of the amorphous Si and Ge constituents. Normal mode analysis indicates that the presence of the Ge inclusion localizes vibrational modes with frequency above the Ge cutoff in the Si host, drastically reducing their ability to transport heat. This observation suggests a general route to achieve exceptionally low thermal conductivity in fully dense solids by restricting the vibrational density of states available for transport in heterogeneous amorphous nanocomposites.