A detailed understanding of the vibrational heat transfer mechanisms between solids is essential for the efficient thermal engineering and control of nanomaterials. We investigate the frequency dependence of anharmonic scattering and interfacial thermal conduction between two acoustically mismatched solids in planar contact by calculating the spectral decomposition of the heat current flowing through an interface between two materials. The calculations are based on analyzing the correlations of atomic vibrations using the data extracted from nonequilibrium molecular dynamics simulations. Inelastic effects arising from anharmonic interactions are shown to significantly facilitate heat transfer between two mass-mismatched face-centered-cubic lattices even at frequencies exceeding the cutoff frequency of the heavier material due to (i) enhanced dissipation of evanescent vibrational modes and (ii) frequency-doubling and frequency-halving three-phonon energy transfer processes at the interface. The results provide substantial insight into interfacial energy transfer mechanisms, especially at high temperatures, where inelastic effects become important and other computational methods are ineffective.

Owing to their long phonon mean free paths (MFPs) and high thermal conductivity, carbon nanotubes (CNTs) are ideal candidates for, e.g., removing heat from electronic devices. It is unknown, however, how the intrinsic phonon MFPs depend on vibrational frequency in nonequilibrium. We determine the spectrally resolved phonon MFPs in isotopically pure CNTs from the spectral phonon transmission function calculated using nonequilibrium molecular dynamics, fully accounting for the resistive phonon-phonon scattering processes through the anharmonic terms of the interatomic potential energy function. Our results show that the effective room temperature MFPs of low-frequency phonons ($f<0.5$ THz) exceed $10\mu \mathrm{m}$, while the MFP of high-frequency phonons ($f\gtrsim 20$ THz) is in the range 10–100 nm. Because the determined MFPs directly reflect the resistance to energy flow, they can be used to accurately predict the thermal conductivity for arbitrary tube lengths by calculating a single frequency integral. The presented results and methods are expected to significantly improve the understanding of nonequilibrium thermal transport in low-dimensional nanostructures.

Understanding the design rules to obtain materials that enable a tight control of phonon transport over a broad range of frequencies would aid major developments in thermoelectric energy harvesting, heat management in microelectronics, and information and communication technology. Using atomistic simulations we show that the metamaterials approach relying on localized resonances is very promising to engineer heat transport at the nanoscale. Combining designed resonant structures to alloying can lead to extremely low thermal conductivity in silicon nanowires. The hybridization between resonant phonons and propagating modes greatly reduces the group velocities and the phonon mean free paths in the low frequency acoustic range below 4 THz. Concurrently, alloy scattering hinders the propagation of high frequency thermal phonons. Our calculations establish a rationale between the size, shape, and period of the resonant structures, and the thermal conductivity of the nanowire, and demonstrate that this approach is even effective to block phonon transport in wavelengths much longer than the size and period of the surface resonant structures. A further consequence of using resonant structures is that they are not expected to scatter electrons, which is beneficial for thermoelectric applications.

The high thermal conductivity of graphene and few-layer graphene undergoes severe degradations through contact with the substrate. Here we show experimentally that the thermal management of a micro heater is substantially improved by introducing alternative heat-escaping channels into a graphene-based film bonded to functionalized graphene oxide through amino-silane molecules. Using a resistance temperature probe for in situ monitoring we demonstrate that the hotspot temperature was lowered by ∼28 °C for a chip operating at 1,300 W cm⁻². Thermal resistance probed by pulsed photothermal reflectance measurements demonstrated an improved thermal coupling due to functionalization on the graphene–graphene oxide interface. Three functionalization molecules manifest distinct interfacial thermal transport behaviour, corroborating our atomistic calculations in unveiling the role of molecular chain length and functional groups. Molecular dynamics simulations reveal that the functionalization constrains the cross-plane phonon scattering, which in turn enhances in-plane heat conduction of the bonded graphene film by recovering the long flexural phonon lifetime.

Thermal transport through liquid-solid interfaces plays an important role in many chemical and biological processes, and better understanding of liquid-solid energy transfer is expected to enable improving the efficiency of thermally driven applications. We determine the spectral distribution of thermal current at liquid-solid interfaces from nonequilibrium molecular dynamics, delivering a detailed picture of the contributions of different vibrational modes to liquid-solid energy transfer. Our results show that surface modes located at the Brillouin zone edge and polarized along the liquid-solid surface normal play a crucial role in liquid-solid energy transfer. Strong liquid-solid adhesion allows also for the coupling of in-plane polarized modes in the solid with the liquid, enhancing the heat-transfer rate and enabling efficient energy transfer up to the cutoff frequency of the solid. Our results provide fundamental understanding of the energy-transfer mechanisms in liquid-solid systems and enable detailed investigations of energy transfer between, e.g., water and organic molecules.