We investigate the effect of uniaxial heterostrain on the interacting phase diagram of magic-angle twisted bilayer graphene. Using both self-consistent Hartree-Fock and density-matrix renormalization group calculations, we find that small strain values (ε∼0.1%-0.2%) drive a zero-temperature phase transition between the symmetry-broken "Kramers intervalley-coherent" insulator and a nematic semimetal. The critical strain lies within the range of experimentally observed strain values, and we therefore predict that strain is at least partly responsible for the sample-dependent experimental observations.

We present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

A growing body of experimental work suggests that magic angle twisted bilayer graphene exhibits a "cascade" of spontaneous symmetry-breaking transitions, sparking interest in the potential relationship between symmetry breaking and superconductivity. However, it has proven difficult to find experimental probes which can unambiguously identify the nature of the symmetry breaking. Here, we show how atomically resolved scanning tunneling microscopy can be used as a fingerprint of symmetry-breaking order. By analyzing the pattern of sublattice polarization and "Kekulé" distortions in small magnetic fields, order parameters for each of the most competitive symmetry-breaking states can be identified. In particular, we show that the "Kramers intervalley coherent state," which theoretical work predicts to be the ground state at even integer fillings, shows a Kekulé distortion which emerges only in a magnetic field.

We investigate twisted bilayer graphene (TBG) at filling ν=-3 in the presence of realistic heterostrain. Strain amplifies the band dispersion and drives the system beyond the strong-coupling regime of previous theoretical studies. We use DMRG to conduct an unbiased, large-scale numerical calculations that include all spin and valley degrees of freedom, up to bond dimension χ=24576. We establish a global phase diagram that unifies a number of theoretical and experimental results. Near zero strain we find an intervalley-coherent quantized anomalous Hall (QAH-IVC) state, a competitive strong-coupling order that evaded past numerical studies. A tiny strain around 0.05% drives a transition into an incommensurate Kekulé spiral (IKS) phase, supporting the mean-field prediction in [Kwan, Phys. Rev. X 11, 041063 (2021)2160-330810.1103/PhysRevX.11.041063]. Even higher strains above 0.2% favor a flavor-symmetric metallic order, which may explain metals found at ν=-3 in many experiments.