No Arabic abstract
We calculate the interactions between the Wannier functions of the 8-orbital model for twisted bilayer graphene (TBG). In this model, two orbitals per valley centered at the AA regions, the AA-p orbitals, account for the most part of the spectral weight of the flats bands. Exchange and assisted-hopping terms between these orbitals are found to be small. Therefore, the low energy properties of TBG will be determined by the density-density interactions. These interactions decay with the distance much faster than in the two orbital model, following a 1/r law in the absence of gates. The magnitude of the largest interaction in the model, the onsite term between the flat band orbitals, is controlled by the size of the AA regions and is estimated to be ~ 40 meV. To screen this interaction, the metallic gates have to be placed at a distance smaller than 5 nm. For larger distances only the long-range part of the interaction is substantially screened. The model reproduces the band deformation induced by doping found in other approaches within the Hartree approximation. Such deformation reveals the presence of other orbitals in the flat bands and is sensitive to the inclusion of the interactions involving them.
We present a simple model that we believe captures the key aspects of the competition between superconducting and insulating states in twisted bilayer graphene. Within this model, the superconducting phase is primary, and arises at generic fillings, but is interrupted by the insulator at commensurate fillings. Importantly, the insulator forms because of electron-electron interactions, but the model is agnostic as to the superconducting pairing mechanism, which need not originate with electron-electron interactions. The model is composed of a collection of crossed one-dimensional quantum wires whose intersections form a superlattice. At each superlattice point, we place a locally superconducting puddle which can exchange Cooper pairs with the quantum wires. We analyze this model assuming weak wire-puddle and wire-wire couplings. We show that for a range of repulsive intrawire interactions, the system is superconducting at `generic incommensurate fillings, with the superconductivity being `interrupted by an insulating phase at commensurate fillings. We further show that the gapped insulating states at commensurate fillings give way to gapless states upon application of external Zeeman fields. These features are consistent with experimental observations in magic-angle twisted bilayer graphenes despite the distinct microscopic details. We further study the full phase diagram of this model and discover that it contains several distinct correlated insulating states, which we characterize herein.
The recently observed superconductivity in twisted bilayer graphene emerges from insulating states believed to arise from electronic correlations. While there have been many proposals to explain the insulating behaviour, the commensurability at which these states appear suggests that they are Mott insulators. Here we focus on the insulating states with $pm 2$ electrons or holes with respect to the charge neutrality point. We show that the theoretical expectations for the Mott insulating states are not compatible with the experimentally observed dependence on temperature and magnetic field if, as frequently assumed, only the correlations between electrons on the same site are included. We argue that the inclusion of non-local (inter-site) correlations in the treatment of the Hubbard model can bring the predictions for the magnetic and temperature dependencies of the Mott transition to an agreement with experiments and have consequences for the critical interactions, the size of the gap, and possible pseudogap physics. The importance of the inter-site correlations to explain the experimental observations indicates that the observed insulating gap is not the one between the Hubbard bands and that antiferromagnetic-like correlations play a key role in the Mott transition.
We introduce and analyze a model that sheds light on the interplay between correlated insulating states, superconductivity, and flavor-symmetry breaking in magic angle twisted bilayer graphene. Using a variational mean-field theory, we determine the normal-state phase diagram of our model as a function of the band filling. The model features robust insulators at even integer fillings, occasional weaker insulators at odd integer fillings, and a pattern of flavor-symmetry breaking at non-integer fillings. Adding a phonon-mediated inter-valley retarded attractive interaction, we obtain strong-coupling superconducting domes, whose structure is in qualitative agreement with experiments. Our model elucidates how the intricate form of the interactions and the particle-hole asymmetry of the electronic structure determine the phase diagram. It also explains how subtle differences between devices may lead to the different behaviors observed experimentally. A similar model can be applied with minor modifications to other moir{e} systems, such as twisted trilayer graphene.
Twisted bilayer graphene exhibits a panoply of many-body phenomena that are intimately tied to the appearance of narrow and well isolated electronic bands near magic-angle. The microscopic ingredients that are responsible for the complex experimental phenomenology include electron-electron (phonon) interactions and non-trivial Bloch wavefunctions associated with the narrow bands. Inspired by recent experiments, we focus here on an interplay of two independent interaction-induced phenomena on superconductivity. We analyze the combined effects of Coulomb interaction driven band-flattening and phonon-mediated attraction due to the exchange of multiple electron-phonon umklapp processes, as a function of filling and twist angle. The former leads to a filling-dependent enhancement of the renormalized density of states, which contributes to a robust increase in the tendency towards pairing in a range of angles near magic-angle. In addition, the minimal spatial extent associated with the Wannier functions develops a non-trivial enhancement as a result of these many-body renormalizations, which can further contribute towards stabilizing the superconducting state over a wider range of fillings and twist-angles.
Moire systems displaying flat bands have emerged as novel platforms to study correlated electron phenomena. Insulating and superconducting states appear upon doping magic angle twisted bilayer graphene (TBG), and there is evidence of correlation induced effects at the charge neutrality point (CNP) which could originate from spontaneous symmetry breaking. Our theoretical calculations show how optical conductivity measurements can distinguish different symmetry breaking states, and reveal the nature of the correlated states. In the specific case of nematic order, which breaks the discrete rotational symmetry of the lattice, we find that the Dirac cones are displaced, not only in momentum space but also in energy, inducing finite Drude weight at the CNP. We also show that the sign of the dc conductivity anisotropy induced by a nematic order depends on the degree of lattice relaxation, the doping and the nature of the symmetry breaking.