Abstract Twisted bilayber graphene (TBG) has revealed fascinating quantum many-body states, and unstrained chargeneutral TBG at 1.08° is understood, both experimentally and theoretically , as a correlated insulator due to the interplay of long-range Coulomb interactions and the quantum metrics of flat bands. However, the fate of the ground state as one continuously tunes the twisted angle is still largely unknown. Here, theoretically, by employing a newly developed momentum-space continuous-field quantum Monte Carlo method that can fully take into account the long-range Coulomb interactions and quantum metrics of flat bands with system sizes that were not possible before , we show that charge-neutral TBG realizes an angle-tuned quantum phase transition with Gross-Neveu criticality at Θc ~ 1.20(1)°, from the analyses of single-particle spectra, free energy and the order parameter of inter-valley coherence. This shows the exciting possibility that by tuning away from 1.08°, pristine TBG can mimic a fundamental phase transition of Dirac fermions.
The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.
Two-dimensional van der Waals materials offer a highly tunable platform for engineering electronic band structures and interactions. By employing techniques such as twisting, gating, or applying pressure, these systems enable precise control over the electronic excitation spectrum. In moir\'e bilayer graphene, the tunability facilitates the transition from a symmetric Dirac semimetal phase through a quantum critical point into an interaction-induced long-range ordered phase with a finite band gap. At charge neutrality, the ordered state proposed to emerge from twist-angle tuning is the Kramers intervalley-coherent insulator. In this case, the transition falls into the quantum universality class of the relativistic Gross-Neveu-XY model in
2+1 dimensions. Here, we refine estimates for the critical exponents characterizing this universality class using an expansion around the lower critical space-time dimension of two. We compute the order-parameter anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\varphi}}$ and the correlation-length exponent $\ensuremath{\nu}$ at one-loop order, and the fermion anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\psi}}$ at two-loop order. Combining these results with previous findings from the expansion around the upper critical dimension, we obtain improved estimates for the universal exponents in
2+1 dimensions via Pad\'e interpolation. For
Nf=4 four-component Dirac fermions, relevant to moir\'e bilayer graphene, we estimate $1/\ensuremath{\nu}=0.916(5), {\ensuremath{\eta}}_{\ensuremath{\varphi}}=0.926(13)$, and ${\ensuremath{\eta}}_{\ensuremath{\psi}}=0.0404(13)$. For
Nf=2, potentially relevant to recent tetralayer
WSe2 experiments, the Gross-Neveu-XY fixed point may be unstable owing to a fixed-point collision at
Nf,c, with ${N}_{\text{f,c}}=1+\sqrt{2}+\mathcal{O}(\ensuremath{\epsilon})$ in the expansion around the lower critical dimension.
We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG). The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation when either the two valleys or the two spin degrees of freedom are considered. By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point (CNP). Then, with the exact Green’s functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings.
Determining the nature of charge carriers is a fundamental goal in the study of strongly correlated electron systems. Here, we employ the continuous-field momentum-space quantum Monte Carlo method to reveal exotic "Dirac trion" excitations in the finite-temperature normal state of twisted bilayer graphene. While the ground state is a symmetry-breaking insulator with gapped (
∼ 20 meV) electron-like excitations, we show that a small temperature (
∼ 3 meV), well below the interaction scale, drives the system into a strongly fluctuating symmetric normal state. We demonstrate that this normal state hosts gapless excitations consisting of three-particle bound states, two electrons and one hole, that are exactly orthogonal to the higher-energy electrons at the zero-momentum gapless point. These Dirac trions have the remarkable property of being arbitrarily light despite being composed of heavy constituents, and their spectra can be easily tuned by varying the twist angle and interlayer hopping strength. Our unbiased quantum many-body computation sheds light on the Dirac trions in a projected correlated flat-band setting and opens the door for further investigation of many-body excitations in strongly correlated topological bands beyond Landau levels.
Abstract One of the most remarkable theoretical findings in magic angle twisted bilayer graphene (TBG) is the emergence of ferromagnetic Slater determinants as exact ground states for the interacting Hamiltonian at the chiral limit. This discovery provides an explanation for the correlated insulating phase which has been experimentally observed at half filling. This work is the first mathematical study of interacting models in magic angle graphene systems. These include not only TBG but also TBG-like systems featuring four flat bands per valley, and twisted trilayer graphene systems with equal twist angles. We identify symmetries of the chiral limit of the Bistritzer-MacDonald Hamiltonian that are responsible for characterizing the Hartree-Fock ground states as zero energy many-body ground states. Furthermore, for a general class of Hamiltonian, we establish criteria that the ferromagnetic Slater determinants are the unique ground states within the class of uniformly half-filled, translation invariant Slater determinants. We then demonstrate that these criteria can be explicitly verified for TBG and TBG-like systems at the chiral limit, using properties of Jacobi-
θ θ and Weierstrass-
℘ ℘ functions.
Abstract Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of ±1. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit.