Search

In multilayer moir\'e heterostructures, the interference of multiple twist angles ubiquitously leads to tunable ultra-long-wavelength patterns known as supermoir\'e lattices. However, their impact on the system's many-body electronic phase diagram remains largely unexplored. We present local compressibility measurements revealing numerous incompressible states resulting from supermoir\'e-lattice-scale isospin symmetry breaking driven by strong interactions. By using the supermoir\'e lattice occupancy as a probe of isospin symmetry, we observe an unexpected doubling of the miniband filling near $\nu=-2$, possibly indicating a hidden phase transition or normal-state pairing proximal to the superconducting phase. Our work establishes supermoir\'e lattices as a tunable parameter for designing novel quantum phases and an effective tool for unraveling correlated phenomena in moir\'e materials.

We propose a minimal "three-patch model" for the anomalous Hall crystal (AHC), a topological electronic state that spontaneously breaks both time-reversal symmetry and continuous translation symmetry. The proposal for this state is inspired by the recently observed integer and fractional quantum Hall states in rhombohedral multilayer graphene at zero magnetic field. There, interaction effects appear to amplify the effects of a weak moir\'e potential, leading to the formation of stable, isolated Chern bands. It has been further shown that Chern bands are stabilized in mean field calculations even without a moir\'e potential, enabling a realization of the AHC state. Our model is built upon the dissection of the Brillouin zone into patches centered around high symmetry points. Within this model, the wavefunctions at high symmetry points fully determine the topology and energetics of the state. We extract two quantum geometrical phases of the non-interacting wavefunctions that control the stability of the topologically nontrivial AHC state. The model predicts that the AHC state wins over the topological trivial Wigner crystal in a wide range of parameters, and agrees very well with the results of full self-consistent Hartree-Fock calculations of the rhombohedral multilayer graphene Hamiltonian., Comment: 15+10 pages, 9+5 figures

The rise of moir\'{e} materials has led to experimental realizations of integer and fractional Chern insulators in small or vanishing magnetic fields. At the same time, a set of minimal conditions sufficient to guarantee a Abelian fractional state in a flat band were identified, namely "ideal" or "vortexable" quantum geometry. Such vortexable bands share essential features with the lowest Landau level, while excluding the need for more fine-tuned aspects such as flat Berry curvature. A natural and important generalization is to ask if such conditions can be extended to capture the quantum geometry of higher Landau levels, particularly the first (1LL), where non-Abelian states at $\nu = 1/2,2/5$ are known to be competitive. The possibility of realizing these states at zero magnetic field , and perhaps even more exotic ones, could become a reality if we could identify the essential structure of the 1LL in Chern bands. In this work, we introduce a precise definition of 1LL quantum geometry, along with a figure of merit that measures how closely a given band approaches the 1LL. We apply the definition to identify two models with 1LL structure -- a toy model of double bilayer twisted graphene and a more realistic model of strained Bernal graphene., Comment: 12 pages, 5 figures

Recent experiments on rhombohedral pentalayer graphene flakes with a substrate induced moir\'e potential have identified both Chern insulators and fractional Quantum Hall states in the absence of an applied magnetic field. Surprisingly, these states are observed in strong displacement fields where the effects of the moir\'e lattice are weak, and seem to be readily accessed without fine-tuning. To address these experimental puzzles we study an interacting model of electrons in this geometry, first within the self-consistent Hartree-Fock (SCHF) approximation. We find an isolated Chern band with Chern number $|C|=1$, that moreover is relatively flat and shows good quantum geometry. Exact diagonalization and density matrix renormalization group methods at fractional filling establish the presence of fractional quantum anomalous Hall (FQAH) states. The $|C|=1$ band in SCHF is remarkably robust to varying microscopic parameters, and is also found in the $N_L=4$ and $N_L=6$ layer systems. Remarkably, it appears stable even to switching off the moir\'e potential, pointing to spontaneous breaking of translation symmetry. We term this topological crystalline state the ``anomalous Hall crystal" (AHC), and argue that it constitutes a general mechanism for creating stable Chern bands in rhombohedral graphene. Our work elucidates the physics behind the recent rhombohedral pentalayer graphene observations, predicts the appearance of the same phase in other systems, and opens the door to studying the interplay between electronic topology and spontaneous translation symmetry breaking., Comment: 5+9 pages, 4+4 figures

The highly tunable nature of synthetic quantum materials -- both in the solid-state and cold atom contexts -- invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moir\'e materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce a minimal 2D lattice model that incorporates exactly these features: time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall (QSH) insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. At the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the QSH insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be $p$-wave -- reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice., Comment: 27 pages (with 8 figures) + 35 pages supplementary (with 14 figures)

The pursuit of exotic phases of matter outside of the extreme conditions of a quantizing magnetic field has been a long standing quest of solid state physics. Recent experiments have observed spontaneous valley polarization and fractional Chern insulators (FCIs) in zero magnetic field in twisted bilayers of MoTe$_2$, at partial filling of the topological valence band ($\nu=-2/3$ and $-3/5$). We study the topological valence band at $\textit{half}$ filling, using exact diagonalization and density matrix renormalization group calculations. We discover a composite Fermi liquid (CFL) phase even at zero magnetic field that covers a large portion of the phase diagram centered around twist angle ${\sim}3.6^\circ$. The CFL is a non-Fermi liquid phase that shows metallic behavior despite the absence of Landau quasiparticles. We discuss experimental implications including the competition between the CFL and a Fermi liquid, which can be tuned with a displacement field. The topological valence band has excellent quantum geometry over a wide range of twist angles and a small bandwidth that is, remarkably, reduced by interactions for a range of angles. These key properties are responsible for stabilizing the exotic quantum Hall phases we find at zero field. Finally, we present an optical signature involving "extinguished" optical responses as a means to identify Chern bands with ideal quantum geometry., Comment: 5+22 pages, 3+13 figures

When the two layers of a twisted moir\'e system are subject to different degrees of strain, the effect is amplified by the inverse twist angle, e.g., by a factor of 50 in magic angle twisted bilayer graphene (TBG). Samples of TBG typically have heterostrains of 0.1-0.7%, increasing the bandwidth of the "flat'' bands by as much as tenfold, placing TBG in an intermediate coupling regime. Here we study the phase diagram of TBG in the presence of heterostrain with unbiased, large-scale density matrix renormalization group calculations (bond dimension $\chi=24576$), including all spin and valley degrees of freedom. Working at filling $\nu = -3$, we find a strain of $0.05\%$ drives a transition from a quantized anomalous Hall insulator into an incommensurate-Kekul\'e spiral (IKS) phase. This peculiar order, proposed and studied at mean-field level by Kwan et al (PRX 11, 041063), breaks both valley conservation and translation symmetry $\hat{T}$, but preserves a modified translation symmetry $\hat{T}'$ with moir\'e-incommensurate phase modulation. Even higher strains drive the system to a fully symmetric metal., Comment: 4.5 pages and three figures, plus appendices

Fractional Chern insulators realize the remarkable physics of the fractional quantum Hall effect (FQHE) in crystalline systems with Chern bands. The lowest Landau level (LLL) is known to host the FQHE, but not all Chern bands are suitable for realizing fractional Chern insulators (FCI). Previous approaches to stabilizing FCIs focused on mimicking the LLL through momentum space criteria. Here instead we take a real-space perspective by introducing the notion of vortexability. Vortexable Chern bands admit a fixed operator that introduces vortices into any band wavefunction while keeping the state entirely within the same band. Vortexable bands admit trial wavefunctions for FCI states, akin to Laughlin states. In the absence of dispersion and for sufficiently short ranged interactions, these FCI states are the ground state -- independent of the distribution of Berry curvature. Vortexable bands are much more general than the LLL, and we showcase a recipe for constructing them. We exhibit diverse examples in graphene-based systems with or without magnetic field, and with any Chern number. A special class of vortexable bands is shown to be equivalent to the momentum space "trace condition" or "ideal band condition". In addition, we also identify a more general form of vortexability that goes beyond this criterion. We introduce a modified measure that quantifies deviations from general vortexability which can be applied to generic Chern bands to identify promising FCI platforms., Comment: 9 pages, 3 figures main text. 26 pages, 4 figures including supplement. Corrections to example E2 of V1

Motivated by recent experiments that observed low-frequency second-order optical responses in doped striped superconductors, here we investigate the nonlinear electrodynamics of systems exhibiting a charge density wave (CDW) order parameter. Due to the Bragg scattering off the CDW order, an incoming spatially homogeneous electric field in addition to zero momentum current generates Umklapp currents that are modulated in space at momenta of the reciprocal CDW lattice. In particular, here we predict and microscopically evaluate the Umklapp shift current, a finite momentum analog of the regular shift current which represents the second-order optical process that downconverts homogeneous AC electric field into low-frequency, zero momentum current. Specifically, we evaluate real-time response functions within mean-field theory via the Keldysh technique and use the Peierls substitution to compute observables at finite momenta in lattice models. We find that systems with certain lattice symmetries (such as inversion symmetry), where the regular shift current is disallowed, may give rise to the Umklapp one. We apply our framework to investigate lattice symmetries in layered materials with helical-like stripes and show that both types of shift currents provide insight into the nature of intertwined phases of matter. Finally, we discuss the relation of our findings to recent experiments in striped superconductors.

Fractional Chern insulators (FCIs) are lattice analogues of fractional quantum Hall states that may provide a new avenue toward manipulating non-abelian excitations. Early theoretical studies have predicted their existence in systems with energetically flat Chern bands and highlighted the critical role of a particular quantum band geometry. Thus far, however, FCI states have only been observed in Bernal-stacked bilayer graphene aligned with hexagonal boron nitride (BLG/hBN), in which a very large magnetic field is responsible for the existence of the Chern bands, precluding the realization of FCIs at zero field and limiting its potential for applications. By contrast, magic angle twisted bilayer graphene (MATBG) supports flat Chern bands at zero magnetic field, and therefore offers a promising route toward stabilizing zero-field FCIs. Here we report the observation of eight FCI states at low magnetic field in MATBG enabled by high-resolution local compressibility measurements. The first of these states emerge at 5 T, and their appearance is accompanied by the simultaneous disappearance of nearby topologically-trivial charge density wave states. Unlike the BLG/hBN platform, we demonstrate that the principal role of the weak magnetic field here is merely to redistribute the Berry curvature of the native Chern bands and thereby realize a quantum band geometry favorable for the emergence of FCIs. Our findings strongly suggest that FCIs may be realized at zero magnetic field and pave the way for the exploration and manipulation of anyonic excitations in moir\'e systems with native flat Chern bands.

We present transport measurements of bilayer graphene with 1.38{\deg} interlayer twist and apparent additional alignment to its hexagonal boron nitride cladding. As with other devices with twist angles substantially larger than the magic angle of 1.1{\deg}, we do not observe correlated insulating states or band reorganization. However, we do observe several highly unusual behaviors in magnetotransport. For a large range of densities around half filling of the moir\'e bands, magnetoresistance is large and quadratic. Over these same densities, the magnetoresistance minima corresponding to gaps between Landau levels split and bend as a function of density and field. We reproduce the same splitting and bending behavior in a simple tight-binding model of Hofstadter's butterfly on a square lattice with anisotropic hopping terms. These features appear to be a generic class of experimental manifestations of Hofstadter's butterfly and may provide insight into the emergent states of twisted bilayer graphene., Comment: 8 pages, 4 figures; updated supplemental material; added additional device measurements and updated toy model

The interplay between strong electron-electron interactions and band topology can lead to novel electronic states that spontaneously break symmetries. The discovery of flat bands in magic-angle twisted bilayer graphene (MATBG) with nontrivial topology has provided a unique platform in which to search for new symmetry-broken phases. Recent scanning tunneling microscopy and transport experiments have revealed a sequence of topological insulating phases in MATBG with Chern numbers $C=\pm 3, \, \pm 2, \, \pm 1$ near moir\'e band filling factors $\nu = \pm 1, \, \pm 2, \, \pm 3$, corresponding to a simple pattern of flavor-symmetry-breaking Chern insulators. Here, we report high-resolution local compressibility measurements of MATBG with a scanning single electron transistor that reveal a new sequence of incompressible states with unexpected Chern numbers observed down to zero magnetic field. We find that the Chern numbers for eight of the observed incompressible states are incompatible with the simple picture in which the $C= \pm 1$ bands are sequentially filled. We show that the emergence of these unusual incompressible phases can be understood as a consequence of broken translation symmetry that doubles the moir\'e unit cell and splits each $C=\pm 1$ band into a $C=\pm 1$ band and a $C=0$ band. Our findings significantly expand the known phase diagram of MATBG, and shed light onto the origin of the close competition between different correlated phases in the system.

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 ($\epsilon \sim 0.1 - 0.2 \%$) drive a zero-temperature phase transition between the symmetry-broken Kramers intervalley-coherent insulator and a nematic semi-metal. 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 opens up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moire materials, in future work., Comment: corrected a typo in arxiv author list

We perform high-field magnetization measurements on the triangular lattice antiferromagnet Fe$_{1/3}$NbS$_2$. We observe a plateau in the magnetization centered at approximately half the saturation magnetization over a wide range of temperature and magnetic field. From density functional theory calculations, we determine a likely set of magnetic exchange constants. Incorporating these constants into a minimal Hamiltonian model of our material, we find that the plateau and of the $Z_3$ symmetry breaking ground state both arise from interplane and intraplane antiferromagnetic interactions acting in competition. These findings are pertinent to the magneto-electric properties of Fe$_{1/3}$NbS$_2$, which allow electrical switching of antiferromagnetic textures at relatively low current densities., Comment: 5 pages, 4 figures, supplemental material on request

We present a new method for compressing matrix product operators (MPOs) which represent sums of local terms, such as Hamiltonians. Just as with area law states, such local operators may be fully specified with a small amount of information per site. Standard matrix product state (MPS) tools are ill-suited to this case, due to extensive Schmidt values that coexist with intensive ones, and Jordan blocks in the transfer matrix. We ameliorate these issues by introducing an "almost Schmidt decomposition" that respects locality. Our method is "$\varepsilon$-close" to the accuracy of MPS-based methods for finite MPOs, and extends seamlessly to the thermodynamic limit, where MPS techniques are inapplicable. In the framework of control theory, our method generalizes Kung's algorithm for model order reduction. Several examples are provided, including an all-MPO version of the operator recursion method (Lanczos algorithm) directly in the thermodynamic limit. All results are accompanied by practical algorithms, well-suited for the large MPOs that arise in DMRG for long-range or quasi-2D models., Comment: 31 pages, 3 figures, 6 practical algorithms; v2: corrected algorithm for triangular iMPO with unit cells; v3: added an Appendix and several corrections suggested by reviewers. v3 is essentially identical to the journal version

We present a hypothesis for the universal properties of operators evolving under Hamiltonian dynamics in many-body systems. The hypothesis states that successive Lanczos coefficients in the continued fraction expansion of the Green's functions grow linearly with rate $\alpha$ in generic systems, with an extra logarithmic correction in 1d. The rate $\alpha$ --- an experimental observable --- governs the exponential growth of operator complexity in a sense we make precise. This exponential growth even prevails beyond semiclassical or large-$N$ limits. Moreover, $\alpha$ upper bounds a large class of operator complexity measures, including the out-of-time-order correlator. As a result, we obtain a sharp bound on Lyapunov exponents $\lambda_L \leq 2 \alpha$, which complements and improves the known universal low-temperature bound $\lambda_L \leq 2 \pi T$. We illustrate our results in paradigmatic examples such as non-integrable spin chains, the Sachdev-Ye-Kitaev model, and classical models. Finally we use the hypothesis in conjunction with the recursion method to develop a technique for computing diffusion constants., Comment: 18+9 pages, 10 figures, 1 table; accepted version

We argue that symmetry-broken phases proximate in phase space to symmetry-protected topological phases can exhibit dynamical signatures of topological physics. This dynamical, symmetry-protected "topological" regime is characterized by anomalously long edge coherence times due to the topological decoration of quasiparticle excitations, even if the underlying zero-temperature ground state is in a non-topological, symmetry-broken state. The dramatic enhancement of coherence can even persist at infinite temperature due to prethermalization. We find exponentially long edge coherence times that are stable to symmetry-preserving perturbations, and not the result of integrability.

Nonlinear optical responses are a crucial probe of physical systems including periodic solids. In the absence of electron-electron interactions, they are calculable with standard perturbation theory starting from the band structure of Bloch electrons, but the resulting formulas are often large and unwieldy, involving many energy denominators from intermediate states. This work gives a Feynman diagram approach to calculating non-linear responses. This diagrammatic method is a systematic way to perform perturbation theory, which often offers shorter derivations and also provides a natural interpretation of nonlinear responses in terms of physical processes. Applying this method to second-order responses concisely reproduces formulas for the second-harmonic, shift current. We then apply this method to third-order responses and derive formulas for third-harmonic generation and self-focusing of light, which can be directly applied to tight-binding models. Third-order responses in the semiclasscial regime include a Berry curvature quadrupole term, whose importance is discussed including symmetry considerations and when the Berry curvature quadrupole becomes the leading contribution. The method is applied to compute third-order optical responses for a model Weyl semimetal, where we find a new topological contribution that diverges in a clean material, as well as resonances with a peculiar linear character., Comment: v2 is as-published. This corrects several small error, particularly on second-order responses

We introduce a systematic construction of a gapless symmetry protected topological phase in one dimension by "decorating" the domain walls of Luttinger liquids. The resulting strongly interacting phases provide a concrete example of a gapless symmetry protected topological (gSPT) phase with robust symmetry-protected edge modes. Using boundary conformal field theory arguments, we show that while the bulks of such gSPT phases are identical to conventional Luttinger liquids, their boundary critical behavior is controlled by a different, strongly-coupled renormalization group fixed point. Our results are checked against extensive density matrix renormalization group calculations., Comment: 10 pages, 3 figures

We introduce exactly solvable gapless quantum systems in $d$ dimensions that support symmetry protected topological (SPT) edge modes. Our construction leads to long-range entangled, critical points or phases that can be interpreted as critical condensates of domain walls "decorated" with dimension $(d-1)$ SPT systems. Using a combination of field theory and exact lattice results, we argue that such gapless SPT systems have symmetry-protected topological edge modes that can be either gapless or symmetry-broken, leading to unusual surface critical properties. Despite the absence of a bulk gap, these edge modes are robust against arbitrary symmetry-preserving local perturbations near the edges. In two dimensions, we construct wavefunctions that can also be interpreted as unusual quantum critical points with diffusive scaling in the bulk but ballistic edge dynamics., Comment: 16 pages, 8 figures. v2: published version

We investigate the entanglement properties of an infinite class of excited states in the quantum Lifshitz model (QLM). The presence of a conformal quantum critical point in the QLM makes it unusually tractable for a model above one spatial dimension, enabling the ground state entanglement entropy for an arbitrary domain to be expressed in terms of geometrical and topological quantities. Here we extend this result to excited states and find that the entanglement can be naturally written in terms of quantities which we dub "entanglement propagator amplitudes" (EPAs). EPAs are geometrical probabilities that we explicitly calculate and interpret. A comparison of lattice and continuum results demonstrates that EPAs are universal. This work shows that the QLM is an example of a 2+1d field theory where the universal behavior of excited-state entanglement may be computed analytically., Comment: Published version. Invited contribution to the special issue of J. Phys. A: "John Cardy's scale-invariant journey in low dimensions: a special issue for his 70th birthday"

Multi-loop scattering amplitudes in N=4 Yang-Mills theory possess cluster algebra structure. In order to develop a computational framework which exploits this connection, we show how to construct bases of Goncharov polylogarithm functions, at any weight, whose symbol alphabet consists of cluster coordinates on the $A_n$ cluster algebra. Using such a basis we present a new expression for the 2-loop 6-particle NMHV amplitude which makes some of its cluster structure manifest., Comment: 32 pages; v2: minor corrections and clarifications