Motivated by recent experiments probing the wavefunctions of magic-angle twisted bilayer graphene (tBLG), we perform large-scale first-principles calculations of tBLG with full atomic relaxation across a wide range of twist angles down to 0. 99^. Focusing on the magic angle, we compute wavefunctions of the low energy bands, resolving atomic-scale details and moiré-scale patterns that form triangular, honeycomb, and Kagome lattices. By tuning the interlayer interactions, we illustrate the formation of the flat bands from isolated monolayers and the emergence of the band inversion and fragile topology at a sufficiently large interaction strength. We identify strong indicators of a new topological phase transition with increasing interlayer interaction strength, achievable with external pressure or a decrease in the twist angle. When this transition occurs, the upper and lower flat bands exchange their wavefunction character, which may explain why superconductivity appears with electron doping below the magic angle. Our study demonstrates the feasibility of using first-principles wavefunctions to help interpret experimental signatures of topological and correlated phases in tBLG.
Zhu et al. (Fri,) studied this question.