Moiré materials have emerged as a rich platform for exploring strong correlation effects in low dimensions, with twisted bilayer graphene (TBG) as a paradigmatic example. To distill the essential ingredients driving moiré-induced phases, a simplified one-dimensional analog -- a two-leg ladder with spatially modulated interleg hopping and a uniform magnetic flux -- was recently introduced. This model, which we refer to as the moiré ladder, features a nearly flat lowest-energy band in a suitable parameter regime, capturing the band-flattening mechanism of TBG. We investigate the ground-state phase diagram of the moiré ladder using a combination of bosonization and density matrix renormalization group (DMRG) techniques, and systematically disentangle the respective roles of the flux and the hopping modulation. At half filling, previous numerical work identified a metal-insulator transition at finite interaction strength and an unexpected ferromagnetic ground state. Revisiting this, we show that the metal-insulator transition can be understood perturbatively within bosonization, governed by the number of Fermi points. In contrast, the ferromagnetic correlations are nonperturbative and require both flux and spatial modulation -- neither alone is sufficient. We extend our analysis to other fillings: one-quarter, three-quarters, slightly above half filling (half filling plus two electrons), and slightly below half filling (half filling minus two electrons). At moderate interactions, we observe ferromagnetism below half filling and antiferromagnetism above; at stronger interactions, ferromagnetism dominates across all studied fillings. Crucially, the analysis demonstrates that periodic interleg hopping alone does not engender new correlated phases; the magnetic flux is essential for the observed unconventional behavior.
We study the single- and many-particle properties of a two-leg ladder model threaded by a flux with the legs coupled by a spatially varying term. Although a priori unrelated to twisted bilayer graphene (TBG), the model is found to have striking similarities: a quasi-flat low-energy band emerges with characteristics similar to that of magic angle TBG. We study the effect of interparticle interaction in our model using the density matrix renormalization group and find that when the band is quasi-flat, the ground state is a ferromagnetic Mott insulator. As the band becomes more dispersive, the system undergoes a ferromagnetic to antiferromagnetic transition. We discuss how our model is relevant not only to magic-angle physics in TBG, but also in the larger context of one-dimensional correlations and magnetism.
Abstract One of the most remarkable theoretical findings in magic angle twisted bilayer graphene (TBG) is the emergence of ferromagnetic Slater determinants as exact ground states for the interacting Hamiltonian at the chiral limit. This discovery provides an explanation for the correlated insulating phase which has been experimentally observed at half filling. This work is the first mathematical study of interacting models in magic angle graphene systems. These include not only TBG but also TBG-like systems featuring four flat bands per valley, and twisted trilayer graphene systems with equal twist angles. We identify symmetries of the chiral limit of the Bistritzer-MacDonald Hamiltonian that are responsible for characterizing the Hartree-Fock ground states as zero energy many-body ground states. Furthermore, for a general class of Hamiltonian, we establish criteria that the ferromagnetic Slater determinants are the unique ground states within the class of uniformly half-filled, translation invariant Slater determinants. We then demonstrate that these criteria can be explicitly verified for TBG and TBG-like systems at the chiral limit, using properties of Jacobi-
θ θ and Weierstrass-
℘ ℘ functions.
Abstract Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of ±1. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit.
Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of ±1. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit.
Mott physics in topological bands has recently attracted considerable attention, particularly in the context of twisted bilayer graphene (TBG). However, the essential ingredients for stabilizing this physics remain unclear. Here, we demonstrate a quantum Mott semimetal phase as the ground state within a one-dimensional spinful Hubbard model featuring only one orbital per unit cell, protected by inversion and particle-hole symmetries. We start from a two-orbital model where a localized
f orbital on the A sublattice hybridizes with a delocalized
c orbital on the B sublattice. Projecting the
f-orbital Hubbard
U onto the active flat band yields a lattice model with Wannier orbitals centered on the B sublattice. Similar to TBG, a momentum-space scale
k∗ emerges, setting the interaction range in the projected model to
1/k∗. While the ground state is ferromagnetic with only the Hubbard
U, introducing an inter-site antiferromagnetic spin coupling
J stabilizes a Mott semimetal
∗ phase with a central charge
c=3. Using exact diagonalization (ED) and density matrix renormalization group (DMRG) methods, we show that this phase hosts a spinful Dirac fermion coexisting with a neutral spin mode -- analogous to the fractionalized Fermi liquid (FL
∗) phase in higher dimensions. Furthermore, breaking particle-hole (PH) symmetry via dispersion transforms the Mott semimetal into a Mott insulator, which is separated from a distinct Mott insulating phase by a continuous transition with a polarization jump of
1/2. Our work provides the first unbiased evidence of a Mott semimetal ground state and demonstrates that this 1D model captures some essential aspects of TBG physics, despite lacking a Wannier obstruction.
We study the dimensional crossover from one to two dimensions in half-filled and lightly doped, weakly interacting N-leg Hubbard ladders. In this case, the Hubbard ladders are equivalent to an N-band model. Using renormalization-group techniques, we find that the half-filled ladders exhibit (in the spin sector) an odd-even effect only below a crossover energy ${E}_{c}\ensuremath{\propto}\mathrm{exp}[\ensuremath{-}\ensuremath{\alpha}\mathrm{exp}(\ensuremath{\gamma}N)]$ $(\ensuremath{\alpha}\ensuremath{\ll}1$ and $\ensuremath{\gamma}\ensuremath{\sim}1$ depend on the interaction strength and on the hopping matrix elements): Below
Ec, the dominant interactions take place within band pairs $(j,N+1\ensuremath{-}j)$ [and within the band
(N+1)/2 for N odd], such that even-leg ladders are an insulating spin liquid, while odd-leg ladders have one gapless spinon mode. In contrast, above the energy scale
Ec, all bands are interacting with each other and the system is a two-dimensional- (2D)-like (insulating) antiferromagnet; we obtain an analytical expression for the Hamiltonian which is similar to the 2D Heisenberg antiferromagnet. Bosonization techniques show that in the charge sector the Mott insulator is as well below and above
Ec of the same type as in
N/2 half-filled two-leg ladders. Doping away from half-filling, we find that the effect of an increasing doping $\ensuremath{\delta}$ is very similar to decreasing the number of legs N: In both cases interactions between unpaired bands are suppressed and thereby the antiferromagnetic correlations reduced. The resulting band pairs form then insulating spin liquids and when doped, there is a spin gap, but phase coherence exists only within the band pairs. At higher doping levels ${\ensuremath{\delta}}_{c}={\ensuremath{\delta}}_{c}(N),$ phase coherence between all band pairs sets in and the system becomes a 2D-like d-wave superconductor $[{\ensuremath{\delta}}_{c}(N)\ensuremath{\rightarrow}0$ for $\stackrel{\ensuremath{\rightarrow}}{N}\ensuremath{\infty}].$