Beyond-mean-field theory of the correlated insulator states in magic-angle twisted bilayer graphene: fluctuation corrections, parameter and screening control, metastability, and the Chern-number criterion

Twisted bilayer graphene near the first magic angle hosts narrow moiré flat bands that develop interaction-induced insulating states at moiré fillings $n=0$ and $n=\pm p/4$ ($p=1,2,3$). The Hartree-Fock study of Xie and MacDonald (Phys. Rev. Lett. 124, 097601 (2020)) established a mean-field account of these states in the full graphene $\pi$-band continuum model with long-range Coulomb interactions, and identified a set of coupled open issues: fluctuations beyond Hartree-Fock, sensitivity to continuum-model parameters and sample-dependent screening, metastability and dynamics of flavor-polarized solutions, the microscopic criteria separating gapped states with nonzero Chern number from those with vanishing Chern number, together with the relation of both to $C_2\mathcal{T}$ breaking, and the lack of a quantitative comparison with experiment. This problem asks for a controlled, quantitative account of the correlated insulating states that resolves these issues jointly: fluctuation-corrected energetics, phase boundaries and collective modes across filling, twist angle, screening, and strain; the classification of when insulating states carry quantized anomalous Hall responses; the metastability and transition kinetics of competing states; and theory-grade comparison with transport and thermodynamic measurements. No fluctuation-corrected, parameter-resolved theory that quantitatively reproduces the observed insulating phases and their dynamics is known in the literature to date, so the problem remains open.

Background

Twisted bilayer graphene (TBG) is two graphene sheets rotated by a small angle $\theta$, which produces a long-wavelength moiré pattern — a periodic superlattice that folds the Dirac cones of graphene into narrow 'flat' bands when $\theta$ is close to the first magic angle $\theta \approx 1.05^\circ$ (the angle at which the Bistritzer-MacDonald continuum model finds the bands become dispersionless, enhancing interaction effects). Experiments starting in 2018 found insulating states at integer fractions of the moiré band filling, $n=\pm p/4$ with integer $p$, together with superconductivity, and in some devices a quantized anomalous Hall effect (a Hall resistance quantized to $h/e^2$ at zero magnetic field, produced by spontaneous magnetic order in a band with nonzero Chern number, the integer topological invariant that counts the winding of the band wavefunctions). The two valley degrees of freedom $K, K'$ of graphene are approximately conserved, and electrons carry an approximately $SU(4)$ spin-valley 'flavor' symmetry; insulating states arise by 'flavor polarization' — occupying some flavors and emptying others so that exchange energy is gained faster than the residual band width costs kinetic energy. In a single valley the only exact symmetries left are the combined two-fold rotation and time reversal $C_2\mathcal{T}$, the three-fold rotation $C_3$, and a mirror. $C_2\mathcal{T}$ protects Dirac points in the single-valley band structure, so a state that breaks $C_2\mathcal{T}$ can open a gap that carries Berry curvature; whether the resulting valley-polarized insulator has nonzero total Chern number then decides whether it exhibits a quantized anomalous Hall effect.

Xie and MacDonald (Phys. Rev. Lett. 124, 097601 (2020)) treated these questions with self-consistent Hartree-Fock (SCHF) calculations — a static mean-field approximation in which each electron moves in the average charge (Hartree) and exchange (Fock) potential of the others — performed in the full $\pi$-band Hilbert space of the continuum model rather than in a truncated flat-band model, with long-range Coulomb interactions characterized by an effective dielectric constant $\epsilon$ that models screening by the encapsulating substrate and gates. Their findings frame the problem: gaps open at neutrality over a wide range of twist angles, sometimes without breaking $C_2\mathcal{T}$; because the flat bands hybridize with remote (higher- and lower-energy) bands unless $\epsilon \gtrsim 20$, flat-band isolation is not guaranteed and $C_2\mathcal{T}$-breaking is not necessary for gap formation at moderate interaction strength; flavor-polarized insulators at $n=\pm p/4$ appear only near the magic angle and are in general only metastable (self-consistently gapped) mean-field solutions; and, since a valley-polarized gapped state may fail to show a quantized anomalous Hall effect either because $C_2\mathcal{T}$ is unbroken or because its Chern number happens to vanish, the microscopic criteria for anomalous Hall behavior at quarter fillings were left undetermined. Experiments (Lu et al., Nature 574, 653 (2019)) indeed showed device-dependent behavior in which some quarter-filled insulators are Chern insulators and some are not. The band topology was found to be sensitive not only to interactions but to the non-interacting model parameters — the ratio $T_{AA}/T_{AB}$ of intrasublattice to intersublattice interlayer tunneling, which encodes lattice relaxation, corrugation, and heterostrain (uniaxial strain that breaks the moiré $C_3$ symmetry).

Problem Statement

Provide a controlled, quantitative theory of the correlated insulating states of magic-angle twisted bilayer graphene that goes beyond the static Hartree-Fock description of Xie and MacDonald (Phys. Rev. Lett. 124, 097601 (2020)) and resolves, within one framework, the issues that study left open:

1. Determine how thermal and quantum fluctuations of the collective fields modify each of the stability, phase boundaries, ordering temperatures, and collective-mode spectra of the insulating states at charge neutrality and at fillings $n=\pm p/4$ with integer $p=1,2,3$, in particular the relative stability of $C_2\mathcal{T}$-breaking versus $C_2\mathcal{T}$-preserving gapped states.
2. Establish how band topology and gap formation depend on the continuum-model parameters (the tunneling ratio $T_{AA}/T_{AB}$, strain and corrugation) and on sample-dependent screening (gate distance and substrate dielectric constant), so that theory can be mapped quantitatively onto specific samples.
3. Characterize the metastability and dynamical accessibility of the flavor-polarized insulating solutions and the kinetics of transitions between competing broken-symmetry states, including domain formation and hysteresis.
4. Classify, as a function of the continuum-model and screening parameters identified in item 2, the microscopic criteria that decide for a gapped state at a given filling whether its total Chern number is nonzero (in which case it exhibits a quantized anomalous Hall effect) or vanishing; treat the breaking of $C_2\mathcal{T}$ and the vanishing of the Chern number as two independent conditions, so that the classification covers gapped states with vanishing Chern number both with and without broken $C_2\mathcal{T}$.
5. Achieve a quantitative comparison with experiment by incorporating screening by proximate gates and substrates, finite-temperature effects, and fluctuation corrections to the mean-field energetics.

The five objectives above are coupled by a single scientific objective: a predictive, fluctuation-corrected account of which insulating states occur in magic-angle twisted bilayer graphene under which controllable conditions, sufficient to reproduce and explain the measured gaps, phase boundaries, switching, and anomalous Hall response of real devices.

References

1. Ming Xie and A. H. MacDonald, "Nature of the Correlated Insulator States in Twisted Bilayer Graphene", Phys. Rev. Lett. 124, 097601 (2020), DOI 10.1103/PhysRevLett.124.097601, https://doi.org/10.1103/PhysRevLett.124.097601 (arXiv:1812.04213, https://arxiv.org/abs/1812.04213)
2. Xiaobo Lu, Petr Stepanov, Wei Yang, Ming Xie, Mohammed Ali Aamir, Ipsita Das, Carles Urgell, Kenji Watanabe, Takashi Taniguchi, Guangyu Zhang, Adrian Bachtold, Allan H. MacDonald, and Dmitri K. Efetov, "Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene", Nature 574, 653 (2019), DOI 10.1038/s41586-019-1695-0, https://doi.org/10.1038/s41586-019-1695-0
3. M. Serlin, C. L. Tschirhart, H. Polshyn, Y. Zhang, J. Zhu, K. Watanabe, T. Taniguchi, L. Balents, and A. F. Young, "Intrinsic quantized anomalous Hall effect in a moiré heterostructure", Science 367, 900 (2020), DOI 10.1126/science.aay5533, https://doi.org/10.1126/science.aay5533
4. Nick Bultinck, Shubhayu Chatterjee, and Michael P. Zaletel, "Mechanism for Anomalous Hall Ferromagnetism in Twisted Bilayer Graphene", Phys. Rev. Lett. 124, 166601 (2020), DOI 10.1103/PhysRevLett.124.166601, https://doi.org/10.1103/PhysRevLett.124.166601
5. Nick Bultinck, Eslam Khalaf, Shang Liu, Shubhayu Chatterjee, Ashvin Vishwanath, and Michael P. Zaletel, "Ground State and Hidden Symmetry of Magic-Angle Graphene at Even Integer Filling", Phys. Rev. X 10, 031034 (2020), DOI 10.1103/PhysRevX.10.031034, https://doi.org/10.1103/PhysRevX.10.031034
6. Shubhayu Chatterjee, Nick Bultinck, and Michael P. Zaletel, "Symmetry breaking and skyrmionic transport in twisted bilayer graphene", Phys. Rev. B 101, 165141 (2020), DOI 10.1103/PhysRevB.101.165141, https://doi.org/10.1103/PhysRevB.101.165141
7. H. Polshyn, J. Zhu, M. A. Kumar, Y. Zhang, F. Yang, C. L. Tschirhart, M. Serlin, K. Watanabe, T. Taniguchi, A. H. MacDonald, and A. F. Young, "Electrical switching of magnetic order in an orbital Chern insulator", Nature 588, 66 (2020), DOI 10.1038/s41586-020-2963-8, https://doi.org/10.1038/s41586-020-2963-8
8. Daniel E. Parker, Tomohiro Soejima, Johannes Hauschild, Michael P. Zaletel, and Nick Bultinck, "Strain-Induced Quantum Phase Transitions in Magic-Angle Graphene", Phys. Rev. Lett. 127, 027601 (2021), DOI 10.1103/PhysRevLett.127.027601, https://doi.org/10.1103/PhysRevLett.127.027601
9. Jie Wang, Yunqin Zheng, Andrew J. Millis, and Jennifer Cano, "Chiral approximation to twisted bilayer graphene: Exact intravalley inversion symmetry, nodal structure, and implications for higher magic angles", Phys. Rev. Research 3, 023155 (2021), DOI 10.1103/PhysRevResearch.3.023155, https://doi.org/10.1103/PhysRevResearch.3.023155
10. Tomohiro Soejima, Daniel E. Parker, Nick Bultinck, Johannes Hauschild, and Michael P. Zaletel, "Efficient simulation of moiré materials using the density matrix renormalization group", Phys. Rev. B 102, 205111 (2020), DOI 10.1103/PhysRevB.102.205111, https://doi.org/10.1103/PhysRevB.102.205111
11. Johannes S. Hofmann, Eslam Khalaf, Ashvin Vishwanath, Erez Berg, and Jong Yeon Lee, "Fermionic Monte Carlo Study of a Realistic Model of Twisted Bilayer Graphene", Phys. Rev. X 12, 011061 (2022), DOI 10.1103/PhysRevX.12.011061, https://doi.org/10.1103/PhysRevX.12.011061
12. Gautam Rai, Lorenzo Crippa, Dumitru Călugăru, Haoyu Hu, Francesca Paoletti, Luca de' Medici, Antoine Georges, B. Andrei Bernevig, Roser Valentí, Giorgio Sangiovanni, and Tim O. Wehling, "Dynamical Correlations and Order in Magic-Angle Twisted Bilayer Graphene", Phys. Rev. X 14, 031045 (2024), DOI 10.1103/PhysRevX.14.031045, https://doi.org/10.1103/PhysRevX.14.031045
13. Ajesh Kumar, Ming Xie, and Allan H. MacDonald, "Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene", Phys. Rev. B 104, 035119 (2021), DOI 10.1103/PhysRevB.104.035119, https://doi.org/10.1103/PhysRevB.104.035119
14. Yves H. Kwan, Yichen Hu, Steven H. Simon, and S. A. Parameswaran, "Exciton Band Topology in Spontaneous Quantum Anomalous Hall Insulators: Applications to Twisted Bilayer Graphene", Phys. Rev. Lett. 126, 137601 (2021), DOI 10.1103/PhysRevLett.126.137601, https://doi.org/10.1103/PhysRevLett.126.137601
15. Bogdan Andrei Bernevig, Biao Lian, Aditya Cowsik, Fang Xie, Nicolas Regnault, and Zhi-Da Song, "Twisted bilayer graphene. V. Exact analytic many-body excitations in Coulomb Hamiltonians: Charge gap, Goldstone modes, and absence of Cooper pairing", Phys. Rev. B 103, 205415 (2021), DOI 10.1103/PhysRevB.103.205415, https://doi.org/10.1103/PhysRevB.103.205415
16. Zachary A. H. Goodwin, Valerio Vitale, Xia Liang, Arash A. Mostofi, and Johannes Lischner, "Hartree theory calculations of quasiparticle properties in twisted bilayer graphene", Electronic Structure 2, 034001 (2020), DOI 10.1088/2516-1075/ab9f94, https://doi.org/10.1088/2516-1075/ab9f94