Search

We study magic angle graphene in the presence of both strain and particle-hole symmetry breaking due to nonlocal interlayer 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 wave vectors. We find that at all nonzero integer fillings very small strains, comparable to those measured in scanning tunneling experiments, stabilize a fundamentally new type of time-reversal-symmetric and spatially nonuniform order. This order, which we dub the “incommensurate Kekulé spiral” (IKS) order, spontaneously breaks both the emergent valley-charge conservation and moiré translation symmetries but preserves a modified translation symmetry T[over ^]^{′}—which simultaneously shifts the spatial coordinates and rotates the U(1) angle which characterizes the spontaneous intervalley 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.

Recent experiments on rhombohedral pentalayer graphene with a substrate-induced moiré potential have identified both Chern insulators and fractional quantum Hall states at zero magnetic field. Surprisingly, these states are observed in strong displacement fields where the effects of the moiré lattice are weak, and seem to be readily accessed without fine-tuning. To address these experimental puzzles, we study a model of interacting electrons in this geometry. Within self-consistent Hartree-Fock (SCHF) calculations, we find an isolated Chern band with small bandwidth and good quantum geometry. Exact diagonalization and density-matrix renormalization group calculations both confirm the band hosts fractional quantum Hall states without a magnetic field. Remarkably, the Chern band is stable at a wide range of angles, at four through six rhombohedral layers, at varying rhombohedral hopping parameters, and-most strikingly-survives in SCHF calculations when the moiré potential vanishes. In this limit, the state spontaneously breaks time-reversal and translation symmetry simultaneously, giving a topological crystalline state that we term the "anomalous Hall crystal." We argue this is a general mechanism to create stable Chern bands in rhombohedral multilayer graphene, opening the door to studying the interplay between electronic topology, fractionalization, and spontaneous translation symmetry breaking.

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_{-} is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_{-} to gapping out the edge. In this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν=1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1)_{2} and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_{-}, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

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 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.

Motivated by conjectures in holography relating the entanglement of purification and reflected entropy to the entanglement wedge cross section, we introduce two related non-negative measures of tripartite entanglement g and h. We prove structure theorems which show that states with nonzero g or h have nontrivial tripartite entanglement. We then establish that in one dimension these tripartite entanglement measures are universal quantities that depend only on the emergent low-energy theory. For a gapped system, we argue that either g≠0 and h=0 or g=h=0, depending on whether the ground state has long-range order. For a critical system, we develop a numerical algorithm for computing g and h from a lattice model. We compute g and h for various CFTs and show that h depends only on the central charge whereas g depends on the whole operator content.