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.
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.
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.
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.
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.
The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.