Fermionic Monte Carlo study of a realistic model of twisted bilayer graphene

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. Our results include (i) the emergence of an insulating Kramers inter-valley coherent ground state in competition with a correlated semi-metal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions which 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.