No Arabic abstract
We study magic angle graphene in the presence of both strain and particle-hole symmetry breaking due to non-local inter-layer tunneling. We perform a self-consistent Hartree-Fock study that incorporates these effects alongside realistic interaction and substrate potentials, and explore a comprehensive set of competing orders including those that break translational symmetry at arbitrary wavevectors. We find that at all non-zero integer fillings very small strains, comparable to those measured in scanning tunneling experiments, stabilize a fundamentally new type of time-reversal symmetric and spatially non-uniform order. This order, which we dub the incommensurate Kekule spiral (IKS) order, spontaneously breaks both the emergent valley-charge conservation and moire translation symmetries, but preserves a modified translation symmetry $hat{T}$ -- which simultaneously shifts the spatial coordinates and rotates the $U(1)$ angle which characterizes the spontaneous inter-valley coherence. We discuss the phenomenological and microscopic properties of this order. We argue that our findings are consistent with all experimental observations reported so far, suggesting a unified explanation of the global phase diagram in terms of the IKS order.
Using exact diagonalization, we study the projected Hamiltonian with Coulomb interaction in the 8 flat bands of first magic angle twisted bilayer graphene. Employing the U(4) (U(4)$times$U(4)) symmetries in the nonchiral (chiral) flat band limit, we reduced the Hilbert space to an extent which allows for study around $ u=pm 3,pm2,pm1$ fillings. In the first chiral limit $w_0/w_1=0$ where $w_0$ ($w_1$) is the $AA$ ($AB$) stacking hopping, we find that the ground-states at these fillings are extremely well-described by Slater determinants in a so-called Chern basis, and the exactly solvable charge $pm1$ excitations found in [arXiv:2009.14200] are the lowest charge excitations up to system sizes $8times8$ (for restricted Hilbert space) in the chiral-flat limit. We also find that the Flat Metric Condition (FMC) used in [arXiv:2009.11301,2009.11872,2009.12376,2009.13530,2009.14200] for obtaining a series of exact ground-states and excitations holds in a large parameter space. For $ u=-3$, the ground state is the spin and valley polarized Chern insulator with $ u_C=pm1$ at $w_0/w_1lesssim0.9$ (0.3) with (without) FMC. At $ u=-2$, we can only numerically access the valley polarized sector, and we find a spin ferromagnetic phase when $w_0/w_1gtrsim0.5t$ where $tin[0,1]$ is the factor of rescaling of the actual TBG bandwidth, and a spin singlet phase otherwise, confirming the perturbative calculation [arXiv:2009.13530]. The analytic FMC ground state is, however, predicted in the intervalley coherent sector which we cannot access [arXiv:2009.13530]. For $ u=-3$ with/without FMC, when $w_0/w_1$ is large, the finite-size gap $Delta$ to the neutral excitations vanishes, leading to phase transitions. Further analysis of the ground state momentum sectors at $ u=-3$ suggests a competition among (nematic) metal, momentum $M_M$ ($pi$) stripe and $K_M$-CDW orders at large $w_0/w_1$.
The dominance of Coulomb interactions over kinetic energy of electrons in narrow, non-trivial moir{e} bands of magic-angle twisted bilayer graphene (TBG) gives rise to a variety of correlated phases such as correlated insulators, superconductivity, orbital ferromagnetism, Chern insulators and nematicity. Most of these phases occur at or near an integer number of carriers per moir{e} unit cell. Experimental demonstration of ordered states at fractional moir{e} band-fillings at zero applied magnetic field $B$, is a challenging pursuit. In this letter, we report the observation of states at half-integer band-fillings of $ u = 0.5$ and $3.5$ at $Bapprox 0$ in a TBG proximitized by a layer of tungsten diselenide (WSe$_2$). The magnetotransport data enables us to deduce features in the underlying band structure consistent with a spontaneously broken translational symmetry supercell with twice the area of the original TBG moir{e} cell. A series of Lifshitz transitions due to the changes in the topology of the Fermi surface implies the evolution of van Hove singularities (VHS) of the diverging density of states at a discrete set of partial fillings of flat bands. Further, we observe reset of charge carriers at $ u = 2, 3$. In addition to magnetotransport, we employ thermoelectricity as a tool to probe the system at $B=0$. Band structure calculations for a TBG moir{e} pattern, together with a commensurate density wave potential and spin-orbit coupling (SOC) terms, allow to obtain degeneracy-lifted, zone-folded moir{e} bands with spin-valley isospin ordering anisotropy that describe the states at half-integer fillings observed experimentally. Our results suggest the emergence of a spin-charge density wave ground state in TBG in the zero $B-$ field limit.
Graphene-based moir{e} systems have attracted considerable interest in recent years as they display a remarkable variety of correlated phenomena. Besides insulating and superconducting phases in the vicinity of integer fillings of the moir{e} unit cell, there is growing evidence for electronic nematic order both in twisted bilayer graphene and twisted double-bilayer graphene (tDBG), as signaled by the spontaneous breaking of the threefold rotational symmetry of the moir{e} superlattices. Here, we combine symmetry-based analysis with a microscopic continuum model to investigate the structure of the nematic phase of tDBG and its experimental manifestations. First, we perform a detailed comparison between the theoretically calculated local density of states and recent scanning tunneling microscopy data [arXiv:2009.11645] to resolve the internal structure of the nematic order parameter in terms of the layer, sublattice, spin, and valley degrees of freedom. We find strong evidence that the dominant contribution to the nematic order parameter comes from states at the moir{e} scale rather than at the microscopic scale of the individual graphene layers, which demonstrates the key role played by the moire degrees of freedom and confirms the correlated nature of the nematic phase in tDBG. Secondly, our analysis reveals an unprecedented tunability of the orientation of the nematic director in tDBG by an externally applied electric field, allowing the director to rotate away from high-symmetry crystalline directions. We compute the expected fingerprints of this rotation in both STM and transport experiments, providing feasible ways to probe it. Rooted in the strong sensitivity of the flat bands of tDBG to the displacement field, this effect opens an interesting route to the electrostatic control of electronic nematicity in moir{e} systems.
Focusing on the twist angle for the minimal commensurate structure, we perform nonperturbative calculations of electron dynamics in the twisted bilayer graphene (TBG) under intense laser fields. We show that the TBG exhibits enriched high-harmonic generation that cannot occur in monolayer or conventional bilayers. We elucidate the mechanism of these nonlinear responses by analyzing dynamical symmetries, momentum-resolved dynamics, and roles of interlayer coupling. Our results imply nonlinear optotwistronics, or controlling optical properties of layered materials by artificial twists.
Recently twisted bilayer graphene (t-BLG) emerges as a new strongly correlated physical platform near a magic twist angle, which hosts many exciting phenomena such as the Mott-like insulating phases, unconventional superconducting behavior and emergent ferromagnetism. Besides the apparent significance of band flatness, band topology may be another critical element in determining strongly correlated twistronics yet receives much less attention. Here we report compelling evidence for nontrivial noninteracting band topology of t-BLG moire Dirac bands through a systematic nonlocal transport study, in conjunction with an examination rooted in $K$-theory. The moire band topology of t-BLG manifests itself as two pronounced nonlocal responses in the electron and hole superlattice gaps. We further show that the nonlocal responses are robust to the interlayer electric field, twist angle, and edge termination, exhibiting a universal scaling law. While an unusual symmetry of t-BLG trivializes Berry curvature, we elucidate that two $Z_2$ invariants characterize the topology of the moire Dirac bands, validating the topological edge origin of the observed nonlocal responses. Our findings not only provide a new perspective for understanding the emerging strongly correlated phenomena in twisted van der Waals heterostructures, but also suggest a potential strategy to achieve topologically nontrivial metamaterials from topologically trivial quantum materials based on twist engineering.