Atomistic Simulations of Thermal Transport Across Interfaces

Jingjing Shi, Purdue University

Abstract

The rapid advance in modern electronics and photonics is pushing device design to the micro- and nano-scale, and the resulting high power density imposes immense challenges to thermal management. Promising materials like carbon nanotubes (CNTs) and graphene offer high thermal conductivity in the axial (or in-plane) directions, but their thermal transport in the radial (or cross-plane) directions are poor, limiting their applications. Hierarchical structures like pillared graphene, which is composed of many CNT-graphene junctions, have been proposed. However, thermal interfacial resistance is a critical issue for thermal management of these systems. In this work, we have systematically explored thermal transport across interfaces, particularly in pillared graphene and silicon/heavy-silicon. First, by recognizing that thermal resistance of the 3D pillared graphene architecture primarily comes from CNT-graphene junctions, a simple network model of thermal transport in pillared graphene structure is developed. Using non-equilibrium molecular dynamics (NEMD), the resistance across an individual CNT-graphene junction with sp2 covalent bonds is found to be around 6 × 10−11 m2K/W, which is significantly lower than typical values reported for planar interfaces between dissimilar materials. Interestingly, when the CNT pillar length is small, the interfacial resistance of the sp2 covalent junction is found to decrease as the CNT pillar length decreases, suggesting the presence of coherence effects. The junction resistance Rj is eventually used in the network model to estimate the effective thermal conductivity, and the results agree well with direct MD simulation data, demonstrating the effectiveness of our model. Then we identify three different mechanisms which can lead to thermal resistances across the pillared graphene junction: the material mismatch (phonon propagates from CNT to graphene), the non-planar junction (the phonon propagation direction must change), and defects (there are six heptagons at each junction). The NEMD results show that three mechanisms lead to similar resistance at the CNT-graphene junction, each at around 2.5 × 10−11 m2K/W. Further, we have predicted the transmission function of individual phonon mode using the wave packet method at CNT-graphene junction. Intriguing phonon polarization conversion behavior is observed for most incident phonon modes. It is found that the polarization conversion dominates the transmission and is more significant at larger phonon wavelength. We attribute such unique phonon polarization conversion behavior to the dimensional mismatch across CNT-graphene interface. It is found that the transmission functions at the junction cannot be predicted by the conventional acoustic mismatch models due to the existence of dimensional mismatch. Further analysis shows that, the dimensionally mismatched interface, on one hand tends to reduce the transmission and conductance due to defects and the change of phonon propagation direction at the interface, while on the other hand tends to enhance the transmission and conductance due to the new phonon transport channel introduced by polarization conversion. Finally, we address that many recent experiments have shown that the measured thermal boundary conductances (TBCs) significantly exceed those calculated using the Landauer approach. We identify that a key assumption that an interface is a local equilibrium system (different modes of phonons on each side of the interfaces are at the emitted phonon temperature Te), is generally invalid and can contribute to the discrepancy.

Degree

Ph.D.

Advisors

Fisher, Purdue University.

Subject Area

Acoustics|Atomic physics|Mathematics|Physics|Thermodynamics

COinS