Investigating Thermal Conductivity in Graphene Moiré Superlattices Through Approach-to-Equilibrium Molecular Dynamics

We investigate the thermal conductivity of graphene Moiré superlattices formed by twisting bilayer graphene (TBG) at small angles, employing approach-to-equilibrium molecular dynamics and lattice dynamics calculations based on the Boltzmann Transport Equation. Our simulations reveal a non-monotonic dependence of the thermal conductivity on the twisting angle, with a local minimum near the first magic angle (θ∼1.1∘). This behavior is attributed to the evolution of local stacking configurations—AA, AB, and saddle-point (SP)—across the Moiré superlattice, which strongly affect phonon transport. A detailed analysis of phonon mean free paths and mode-resolved thermal conductivity shows that AA stacking suppresses thermal transport while, AB and SP stackings exhibit enhanced conductivity owing to more efficient low-frequency phonon transport. Furthermore, we establish a direct correlation between the thermal conductivity of twisted structures and the relative abundance of stacking domains within the Moiré supercell. Our results demonstrate that even very small changes in twisting angle (<2∘) can lead to thermal conductivity variations of over 30%, emphasizing the high tunability of thermal transport in TBG.