We use unrestricted Hartree-Fock, density matrix renormalization group, and variational projected entangled pair state calculations to investigate the ground state phase diagram of the triangular lattice Hubbard model at "half doping" relative to single occupancy, i.e. at a filling of $(1\pm \frac{1}{2})$ electrons per site. The electron-doped case has a nested Fermi surface in the non-interacting limit, and hence a weak-coupling instability towards density-wave orders whose wavevectors are determined by Fermi surface nesting conditions. We find that at moderate to strong interaction strengths other spatially-modulated orders arise, with wavevectors distinct from the nesting vectors. In particular, we identify a series closely-competing itinerant long-wavelength magnetically ordered states, yielding to uniform ferromagnetic order at the largest interaction strengths. For half-hole doping and a similar range of interaction strengths, our data indicate that magnetic orders are most likely absent.

We construct an interacting Wannier model for both AA-stacked and AB-stacked twisted SnSe2, revealing a rich landscape of correlated quantum phases. For the AA-stacked case, the system is effectively described by a three-orbital triangular lattice model, where each orbital corresponds to a valley and exhibits an approximate one-dimensional hopping structure due to a new momentum-space non-symmorphic symmetry. By exploring the interacting phase diagram using a combination of theoretical methods, including Hartree-Fock mean-field theory and exact solutions of the spin model in certain limits, we identify several exotic quantum phases. These include a dimerized phase with finite residual entropy, valence bond solids, and quantum paramagnetism. In the AB-stacked case, the system realizes an interacting kagome lattice model, where the Wannier orbitals associated with the three valleys form three sublattices. In the strong coupling regime, we use cluster mean-field methods to demonstrate the emergence of a classical spin liquid phase due to the frustrated lattice structure. The high tunability of the moiré system, which allows control over both the filling and interaction strength (via twist angle), renders twisted SnSe2 a versatile platform for realizing a wide range of exotic correlated quantum phases.

We introduce and study new interacting topological states that arise in time-reversal symmetric bands with an underlying obstruction to forming localized states. If the U(1) valley symmetry linked to independent charge conservation in each time-reversal sector is spontaneously broken, the corresponding `excitonic' order parameter is forced to form a topologically non-trivial texture across the Brillouin zone. We show that the resulting phase, which we dub a textured exciton insulator, cannot be given a local-moment description due to a form of delicate topology. Using toy models of bands with Chern or Euler obstructions to localization we construct explicit examples of the Chern or Euler texture insulators (CTIs or ETIs) they support, and demonstrate that these are generically competitive ground states at intermediate coupling. We construct field theories that capture the response properties of these new states. Finally, we identify the incommensurate Kekulé spiral phase observed in magic-angle bi- and trilayer graphene as a concrete realization of an ETI.

We show that Hubbard models with nearest-neighbor hopping and a nearest-neighbor hardcore constraint exhibit `maximal' Hilbert space fragmentation in many lattices of arbitrary dimension $d$. Focusing on the $d=1$ rhombus chain and the $d=2$ Lieb lattice, we demonstrate that the fragmentation is strong for all fillings in the thermodynamic limit, and explicitly construct all emergent integrals of motion, which include an extensive set of higher-form symmetries. Blockades consisting of frozen particles partition the system in real space, leading to anomalous dynamics. Our results are potentially relevant to optical lattices of dipolar and Rydberg-dressed atoms.

We explore the phase diagrams of moir\'e materials in search of a new class of intervalley-coherent correlated insulating state: the Chern texture insulator (CTI). This phase of matter, proposed in a companion paper, breaks valley $U(1)$ symmetry in a nontrivial fashion wherein the valley order parameter is forced to texture in momentum space as a consequence of band topology. Using detailed Hartree-Fock studies, we establish that the CTI emerges as an energetically competitive intermediate-coupling ground state in several moir\'e systems which lack a twofold rotation symmetry that forbids the single-particle topology essential to the formation of the CTI valley texture.

We study the nonunitary dynamics of a class of quantum circuits based on stochastically measuring check operators of subsystem quantum error-correcting codes, such as the Bacon-Shor code and its various generalizations. Our focus is on how properties of the underlying code are imprinted onto the measurement-only dynamics. We find that in a large class of codes with nonlocal stabilizer generators, at late times there is generically a nonlocal contribution to the subsystem entanglement entropy that scales with the subsystem size. The nonlocal stabilizer generators can also induce slow dynamics, since depending on the rate of competing measurements the associated degrees of freedom can take exponentially long (in system size) to purify (i.e., disentangle from the environment when starting from a mixed state) and to scramble (i.e., become entangled with the rest of the system when starting from a product state). Concretely, we consider circuits in which the nonlocal stabilizer generators of the underlying subsystem code take the form of subsystem symmetries. We present a systematic study of the phase diagrams and relevant time scales in two and three spatial dimensions for both Calderbank-Shor-Steane (CSS) and non-CSS codes, focusing in particular on the link between slow measurement-only dynamics and the geometry of the subsystem symmetry. A key finding of our paper is that slowly purifying or scrambling degrees of freedom appear to emerge only in codes for which subsystem symmetries are nonlocally generated, a strict subset of those whose symmetries are simply nonlocal. We comment on the link between our results on subsystem codes and the phenomenon of Hilbert-space fragmentation in light of their shared algebraic structure.

Graphene moir\'e systems are ideal environments for investigating complex phase diagrams and gaining fundamental insights into the mechanisms underlying exotic states of matter, as they permit controlled manipulation of electronic properties. Magic-angle twisted trilayer graphene (MATTG) has emerged as a key platform to explore moir\'e superconductivity, owing to the robustness of its superconducting order and the displacement-field tunability of its energy bands. Recent measurements strongly suggest that superconductivity in MATTG is unconventional. Here, we report the first direct observation of double-dome superconductivity in MATTG. The temperature, magnetic field, and bias current dependence of the superconductivity of doped holes collectively show that it is significantly suppressed near moir\'e filling $\nu^* = -2.6$, leading to a double dome in the phase diagram within a finite window of the displacement field. The temperature dependence of the normal-state resistance and the $I-V$ curves straddling $\nu^*$ are suggestive of a phase transition and the potentially distinct nature of superconductivity in the two domes. Hartree-Fock calculations incorporating mild strain yield an incommensurate Kekul\'e spiral state whose effective spin polarization peaks in the regime where superconductivity is suppressed in experiments. This allows us to draw conclusions about the normal state as well as the unconventional nature of the superconducting order parameter.

It is a matter of current debate whether the gate-tunable superconductivity in twisted bilayer graphene is phonon mediated or arises from electron-electron interactions. The recent observation of the strong coupling of electrons to so-called K-phonon modes in angle-resolved photoemission spectroscopy experiments has resuscitated early proposals that K phonons drive superconductivity. We show that the bandwidth-enhancing effect of interactions drastically weakens both the intrinsic susceptibility towards pairing as well as the screening of Coulomb repulsion that is essential for the phonon attraction to dominate at low temperature. This rules out purely K-phonon-mediated superconductivity with the observed transition temperature of -1 K. We conclude that the unflattening of bands by Coulomb interactions challenges any purely phonon-driven pairing mechanism, and must be addressed by a successful theory of superconductivity in moir & eacute; graphene.

We study the phase diagram of magic-angle twisted symmetric trilayer graphene in the presence of uniaxial heterostrain and interlayer displacement field. For experimentally reasonable strain, our mean-field analysis finds robust Kekulé spiral order whose doping-dependent ordering vector is incommensurate with the moiré superlattice, consistent with recent scanning tunneling microscopy experiments, and paralleling the behavior of closely related twisted bilayer graphene (TBG) systems. Strikingly, we identify a possibility absent in TBG: the existence of Kekulé spiral order even at zero strain for experimentally realistic values of the interlayer potential in a trilayer. Our studies also reveal a complex pattern of charge transfer between weakly and strongly dispersive bands in strained trilayer samples as the density is tuned by electrostatic gating, that can be understood intuitively in terms of the “cascades” in the compressibility of magic-angle TBG. Published by the American Physical Society 2024

Recent scanning tunneling microscopy experiments [K.P. Nuckolls et al., arXiv:2303.00024] have revealed the ubiquity of Kekul\'e charge-density wave order in magic-angle twisted bilayer graphene. Most samples are moderately strained and show `incommensurate Kekul\'e spiral' (IKS) order involving a graphene-scale charge density distortion uniaxially modulated on the scale of the moir\'e superlattice, in accord with theoretical predictions. However, ultra-low strain samples instead show graphene-scale Kekul\'e charge order that is uniform on the moir\'e scale. This order, especially prominent near filling factor $\nu=-2$, is unanticipated by theory which predicts a time-reversal breaking Kekul\'e current order at low strain. We show that including the coupling of moir\'e electrons to graphene-scale optical zone-corner (ZC) phonons stabilizes a uniform Kekul\'e charge ordered state at $|\nu|=2$ with a quantized topological (spin or anomalous Hall) response. Our work clarifies how this phonon-driven selection of electronic order emerges in the strong-coupling regime of moir\'e graphene.

We show that the presence of anyons in the excitation spectrum of a two-dimensional system can be inferred from nonlinear spectroscopic quantities. In particular, we consider pump-probe spectroscopy, where a sample is irradiated by two light pulses with an adjustable time delay between them. The relevant response coefficient exhibits a universal form that originates from the statistical phase acquired when anyons created by the first pulse braid around those created by the second. This behavior is shown to be qualitatively unchanged by nonuniversal physics including nonstatistical interactions and small nonzero temperatures. In magnetic systems, the signal of interest can be measured using currently available terahertz-domain probes, highlighting the potential usefulness of nonlinear spectroscopic techniques in the search for quantum spin liquids.

We study a model of spinless fermions on the honeycomb lattice with nearest-neighbor exclusion and extended repulsive interactions that exhibits `lattice supersymmetry' [P. Fendley, K. Schoutens, and J. de Boer, Phys. Rev. Lett. 90, 120402 (2003)]. Using a combination of exact diagonalization of large ($N\leq56$ site) systems, mean-field numerics, and symmetry analysis, we establish a rich phase structure as a function of fermion density, that includes non-Fermi liquid behavior, resonating charge stripes, domain-wall and bubble physics, and identify a finite range of fillings with extensive ground state degeneracy and both gapped and gapless spectra. We comment on the stability of our results to relaxing the stringent requirements for supersymmetry, and on their possible broader relevance to systems of strongly-correlated electrons with extended repulsive interactions.

We study the behavior of linear and nonlinear spectroscopic quantities in two-dimensional topologically ordered systems, which host anyonic excitations exhibiting fractional statistics. We highlight the role that braiding phases between anyons have on the dynamics of such quasiparticles, which as we show dictates the behavior of both linear response coefficients at finite temperatures, as well as nonlinear pump-probe response coefficients. These quantities, which act as probes of temporal correlations in the system, are shown to obey distinctive universal forms at sufficiently long timescales. As well as providing an experimentally measurable fingerprint of anyonic statistics, the universal behavior that we find also demonstrates anomalously fast thermal relaxation: correlation functions decay as a "squished exponential" C(t) similar to exp(-[t/tau]3/2) at long times. We attribute this unusual asymptotic form to the nonlocal nature of interactions between anyons, which allows relaxation to occur much faster than in systems with quasiparticles interacting via local, nonstatistical interactions. While our results apply to any Abelian or non-Abelian topological phase in two-dimensions, we discuss in particular the implications for candidate quantum spin liquid materials, wherein the relevant quantities can be measured using pre-existing time-resolved terahertz-domain spectroscopic techniques.

Owing to the negligible spin-orbit coupling in graphene, there is an additional near-exact SU(2) spin symmetry present throughout, that we will not comment on explicitly.

We introduce an SU(M)-symmetric disordered bipartite spin model with unusual characteristics. Although superficially similar to the Sachdev-Ye model, it has several markedly different properties for M>2. In particular, it has a large non-trivial nullspace whose dimension grows exponentially with system size. The states in this nullspace are frustration-free, and are ground states when the interactions are ferromagnetic. The exponential growth of the nullspace leads to Hilbert-space fragmentation and a violation of the eigenstate thermalization hypothesis. We demonstrate that the commutant algebra responsible for this fragmentation is a non-trivial subalgebra of the Read-Saleur commutant algebra of certain nearest-neighbour models such as the spin-1 biquadratic spin chain. We also discuss the low-energy behaviour of correlations for the disordered version of this model in the limit of a large number of spins and large M, using techniques similar to those applied to the SY model. We conclude by generalizing the Shiraishi-Mori embedding formalism to non-local models, and apply it to turn some of our nullspace states into quantum many-body scars.

A striking series of experiments have observed superconductivity in Bernal-stacked bilayer graphene (BBG) when the energy bands are flattened by applying an electrical displacement field. Intriguingly, superconductivity manifests only at non-zero magnetic fields, or when spin-orbit coupling is induced in BBG by coupling to a substrate. We present detailed functional renormalization group and random-phase approximation calculations that provide a unified explanation for the superconducting mechanism in both cases. Both calculations yield a purely electronic $p$-wave instability of the Kohn-Luttinger (KL) type. The latter can be enhanced either by magnetic fields or Ising spin-orbit coupling, naturally explaining the behaviour seen in experiments.

We study close-packed dimers on the quasiperiodic Ammann-Beenker (AB) graph, that was recently shown to have the unusual feature that hard-core dimer constraints are exactly reproduced at successive discrete length scales. This observation led to a conjecture that it would be possible to construct an exact real-space decimation scheme where each iteration preserves both the quasiperiodic tiling structure and the constraint. Here, we confirm this conjecture by explicitly constructing the corresponding renormalization group transformation and show, using large-scale Monte Carlo simulations, that the dimer distributions flow to a fixed point with non-zero dimer potentials. We use the fixed-point Hamiltonian to demonstrate the existence of slowly decaying dimer correlations. We thus identify a remarkable example of a classical statistical mechanical model whose properties are controlled by the fixed point of an exact renormalization group procedure exhibiting discrete scale invariance but lacking translational and continuous rotational symmetries.

We demonstrate that nonlinear response functions in many-body systems carry a sharp signature of interactions between gapped low-energy quasiparticles. Such interactions are challenging to deduce from linear response measurements. The signature takes the form of a divergent-in-time contribution to the response-linear in time in the case when quasiparticles propagate ballistically-that is absent for free bosonic excitations. We give a physically transparent semiclassical picture of this singular behavior. While the semiclassical picture applies to a broad class of systems we benchmark it in two simple models: in the Ising chain using a form-factor expansion, and in a nonintegrable model-the spin-1 Affleck-Kennedy-Lieb-Tasaki chain-using time-dependent density matrix renormalization group simulations. We comment on extensions of these results to finite temperatures.

The flat electronic bands in magic-angle twisted bilayer graphene (MATBG) host a variety of correlated insulating ground states, many of which are predicted to support charged excitations with topologically non-trivial spin and/or valley skyrmion textures. However, it has remained challenging to experimentally address their ground state order and excitations, both because some of the proposed states do not couple directly to experimental probes, and because they are highly sensitive to spatial inhomogeneities in real samples. Here, using a scanning single-electron transistor, we observe thermodynamic gaps at even integer moiré filling factors at low magnetic fields. We find evidence of a field-tuned crossover from charged spin skyrmions to bare particle-like excitations, suggesting that the underlying ground state belongs to the manifold of strong-coupling insulators. From the spatial dependence of these states and the chemical potential variation within the flat bands, we infer a link between the stability of the correlated ground states and local twist angle and strain. Our work advances the microscopic understanding of the correlated insulators in MATBG and their unconventional excitations.

We present microscopic models of spin ladders which exhibit continuous critical surfaces whose properties and existence, unusually, cannot be inferred from those of the flanking phases. These models exhibit either `multiversality' -- the presence of different universality classes over finite regions of a critical surface separating two distinct phases -- or its close cousin, `unnecessary criticality'-- the presence of a stable critical surface within a single, possibly trivial, phase. We elucidate these properties using Abelian bosonization and density-matrix renormalization-group simulations, and attempt to distill the key ingredients required to generalize these considerations.