No Arabic abstract
Electron-electron interactions are intrinsically long ranged, but many models of strongly interacting electrons only take short-ranged interactions into account. Here, we present results of atomistic calculations including both long-ranged and short-ranged electron-electron interactions for the magnetic phase diagram of twisted bilayer graphene and demonstrate that qualitatively different results are obtained when long-ranged interactions are neglected. In particular, we use Hartree theory augmented with Hubbard interactions and calculate the interacting spin susceptibility at a range of doping levels and twist angles near the magic angle to identify the dominant magnetic instabilities. At the magic angle, mostly anti-ferromagnetic order is found, while ferromagnetism dominates at other twist angles. Moreover, long-ranged interactions significantly increase the twist angle window in which strong correlation phenomena can be expected. These findings are in good agreement with available experimental data.
Twisted bilayer graphene (tBLG) has recently emerged as a platform for hosting correlated phenomena, owing to the exceptionally flat band dispersion that results near interlayer twist angle $thetaapprox1.1^circ$. At low temperature a variety of phases are observed that appear to be driven by electron interactions including insulating states, superconductivity, and magnetism. Electrical transport in the high temperature regime has received less attention but is also highly anomalous, exhibiting gigantic resistance enhancement and non-monotonic temperature dependence. Here we report on the evolution of the scattering mechanisms in tBLG over a wide range of temperature and for twist angle varying from 0.75$^circ$ - 2$^circ$. We find that the resistivity, $rho$, exhibits three distinct phenomenological regimes as a function of temperature, $T$. At low $T$ the response is dominated by correlation and disorder physics; at high $T$ by thermal activation to higher moire subbands; and at intermediate temperatures $rho$ varies linearly with $T$. The $T$-linear response is much larger than in monolayer graphenefor all measured twist angles, and increases by more than three orders of magnitude for $theta$ near the flat-band condition. Our results point to the dominant role of electron-phonon scattering in twisted layer systems, with possible implications for the origin of the observed superconductivity.
We derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4)$times$U(4) symmetry, and we find that the exact ground states at integer filling $-4le ule 4$ relative to charge neutrality are Chern insulators of Chern numbers $ u_C=4-| u|,2-| u|,cdots,| u|-4$, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $ u=0,pm2$ has exact ground states of Chern number $0$, while $ u=pm1,pm3$ has perturbative ground states of Chern number $ u_C=pm1$, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $ u_C=0$ ($0<| u_C|<4-| u|$) at integer fillings $ u$ are fully (partially) intervalley coherent, while the insulator states with Chern number $| u_C|=4-| u|$ are valley polarized. However, for $0<| u_C|le4-| u|$, the fully intervalley coherent states are highly competitive (0.005meV/electron higher). At nonzero magnetic field $|B|>0$, a first-order phase transition for $ u=pm1,pm2$ from Chern number $ u_C=text{sgn}( u B)(2-| u|)$ to $ u_C=text{sgn}( u B)(4-| u|)$ is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.
In the vicinity of the magic angle in twisted bilayer graphene (TBG), the two low-energy van Hove singularities (VHSs) become exceedingly narrow1-10 and many exotic correlated states, such as superconductivity, ferromagnetism, and topological phases, are observed11-16. Heterostrain, which is almost unavoidable in the TBG, can modify its single-particle band structure and lead to novel properties of the TBG that have never been considered so far. Here, we show that heterostrain in a TBG near the magic angle generates a new zero-energy flat band between the two VHSs. Doping the TBG to partially fill the zero-energy flat band, we observe a correlation-induced gap of about 10 meV that splits the flat band. By applying perpendicular magnetic fields, a large and linear response of the gap to magnetic fields is observed, attributing to the emergence of large orbital magnetic moments in the TBG when valley degeneracy of the flat band is lifted by electron-electron interactions. The orbital magnetic moment per moire supercell is measured as about 15 uB in the TBG.
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 explore in detail the electronic phases of a system consisting of three non-colinear arrays of coupled quantum wires, each rotated 120 degrees with respect to the next. A perturbative renormalization-group analysis reveals that multiple correlated states can be stabilized: a $s$-wave or $d pm id$ superconductor, a charge density wave insulator, a two-dimensional Fermi liquid, and a 2D Luttinger liquid (also known as smectic metal or sliding Luttinger liquid). The model provides an effective description of electronic interactions in small-angle twisted bilayer graphene and we discuss its implications in relation to the recent observation of correlated and superconducting groundstates near commensurate densities in magic-angle twisted samples, as well as the ``strange metal behavior at finite temperatures as a natural outcome of the 2D Luttinger liquid phase.