Observation of Reentrant Correlated Insulators and Interaction-Driven Fermi-Surface Reconstructions at One Magnetic Flux Quantum per Moiré Unit Cell in Magic-Angle Twisted Bilayer Graphene

Emergent phases in graphene flat bands

Electronic correlations in two-dimensional materials play a crucial role in stabilising emergent phases of matter. The realisation of correlation-driven phenomena in graphene has remained a longstanding goal, primarily due to the absence of strong electron-electron interactions within its low-energy bands. In this context, magic-angle twisted bilayer graphene has recently emerged as a novel platform featuring correlated phases favoured by the low-energy flat bands of the underlying moiré superlattice. Notably, the observation of correlated insulators and superconductivity, and the interplay between these phases have garnered significant attention. A wealth of correlated phases with unprecedented tunability was discovered subsequently, including orbital ferromagnetism, Chern insulators, strange metallicity, density waves, and nematicity. However, a comprehensive understanding of these closely competing phases remains elusive. The ability to controllably twist and stack multiple graphene layers has enabled the creation of a whole new family of moiré superlattices with myriad properties. Here, we review the progress and development achieved so far, encompassing the rich phase diagrams offered by these graphene-based moiré systems. Additionally, we discuss multiple phases recently observed in non-moiré multilayer graphene systems. Finally, we outline future opportunities and challenges for the exploration of hidden phases in this new generation of moiré materials.

Topological heavy fermions in magnetic field

The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization of THFM to include the magnetic field B is desired, but currently missing. Here we provide a systematic derivation of the THFM in B and solve the resulting model to obtain the interacting Hofstadter spectra for single particle charged excitations. While naive minimal substitution within THFM fails to correctly account for the total number of magnetic subbands within the narrow band i.e., its total Chern number, our method-based on projecting the light and heavy fermions onto the irreducible representations of the magnetic translation group- reproduces the correct total Chern number. Analytical results presented here offer an intuitive understanding of the nature of the (strongly interacting) Hofstadter bands.

Double and Quadruple Flat Bands Tuned by Alternative Magnetic Fluxes in Twisted Bilayer Graphene

Twisted bilayer graphene (TBG) can host the moiré energy flat bands with twofold degeneracy serving as a fruitful playground for strong correlations and topological phases. However, the number of degeneracy is not limited to two. Introducing a spatially alternative magnetic field, we report that the induced magnetic phase becomes an additional controllable parameter and leads to an undiscovered generation of fourfold degenerate flat bands. This emergence stems from the band inversion at the Γ point near the Fermi level with a variation of both twisted angle and magnetic phase. We present the conditions for the emergence of multifold degenerate flat bands, which are associated with the eigenvalue degeneracy of a Birman-Schwinger operator. Using holomorphic functions, which explain the origin of the double flat bands in the conventional TBG, we can generate analytical wave functions in the magnetic TBG to show absolute flatness with fourfold degeneracy. Moreover, we identify an orbital-related intervalley coherent state as the many-body ground state at charge neutrality. In contrast, the conventional TBG has only two moiré energy flat bands, and the highly degenerate flat bands with additional orbital channels in this magnetic platform might bring richer correlation physics.

Emergent Correlated Phases in Rhombohedral Trilayer Graphene Induced by Proximity Spin-Orbit and Exchange Coupling

The impact of proximity-induced spin-orbit and exchange coupling on the correlated phase diagram of rhombohedral trilayer graphene (RTG) is investigated theoretically. By employing ab initio-fitted effective models of RTG encapsulated by transition metal dichalcogenides (spin-orbit proximity effect) and ferromagnetic Cr_{2}Ge_{2}Te_{6} (exchange proximity effect), we incorporate the Coulomb interactions within the random-phase approximation to explore potential correlated phases at different displacement fields and doping. We find a rich spectrum of spin-valley resolved Stoner and intervalley coherence instabilities induced by the spin-orbit proximity effects, such as the emergence of a spin-valley-coherent phase due to the presence of valley-Zeeman coupling. Similarly, proximity exchange removes the phase degeneracies by biasing the spin direction, enabling a magnetocorrelation effect-strong sensitivity of the correlated phases to the relative magnetization orientations (parallel or antiparallel) of the encapsulating ferromagnetic layers.

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings ν=0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at ν=0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023).PRLTAO0031-900710.1103/PhysRevLett.131.026502, Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).PRLTAO0031-900710.1103/PhysRevLett.129.047601] and find the ordered states to have lower energy at ν=0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at ν=0,±1,±2 down to ∼2  K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling ν=±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Finally, we study the effect of a C_{3z}-rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C_{3z} symmetry) can be stabilized at sufficiently large strain at ν=0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund's Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

We apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.047601], to obtain its Kondo lattice limit. In this limit localized f electrons on a triangular lattice interact with topological conduction c electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0,±1,±2 are controlled by the heavy f electrons, while ν=±3 is at the border of a phase transition between two f-electron fillings. For ν=0,±1,±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the f moments in the full model and analytically prove the SU(4) Hund's rule for the ground state which maintains that two f electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite "bond-operators" method. We then compute the excitation spectrum of the f moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the f moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

Hofstadter Topology with Real Space Invariants and Reentrant Projective Symmetries

Adding magnetic flux to a band structure breaks Bloch's theorem by realizing a projective representation of the translation group. The resulting Hofstadter spectrum encodes the nonperturbative response of the bands to flux. Depending on their topology, adding flux can enforce a bulk gap closing (a Hofstadter semimetal) or boundary state pumping (a Hofstadter topological insulator). In this Letter, we present a real space classification of these Hofstadter phases. We give topological indices in terms of symmetry-protected real space invariants, which reveal the bulk and boundary responses of fragile topological states to flux. In fact, we find that the flux periodicity in tight-binding models causes the symmetries which are broken by the magnetic field to reenter at strong flux where they form projective point group representations. We completely classify the reentrant projective point groups and find that the Schur multipliers which define them are Arahanov-Bohm phases calculated along the bonds of the crystal. We find that a nontrivial Schur multiplier is enough to predict and protect the Hofstadter response with only zero-flux topology.

Unconventional self-similar Hofstadter superconductivity from repulsive interactions

Fractal Hofstadter bands have become widely accessible with the advent of moiré superlattices, opening the door to studies of the effect of interactions in these systems. In this work we employ a renormalization group (RG) analysis to demonstrate that the combination of repulsive interactions with the presence of a tunable manifold of Van Hove singularities provides a new mechanism for driving unconventional superconductivity in Hofstadter bands. Specifically, the number of Van Hove singularities at the Fermi energy can be controlled by varying the flux per unit cell and the electronic filling, leading to instabilities toward nodal superconductivity and chiral topological superconductivity with Chern number [Formula: see text]. The latter is characterized by a self-similar fixed trajectory of the RG flow and an emerging self-similarity symmetry of the order parameter. Our results establish Hofstadter quantum materials such as moiré heterostructures as promising platforms for realizing novel reentrant Hofstadter superconductors.

Moiré Phonons in Magic-Angle Twisted Bilayer Graphene

Magic-angle twisted bilayer graphene (TBG) has attracted significant interest recently due to the discoveries of diverse correlated and topological states. In this work, we study the phonon properties in magic-angle TBG based on many-body classical potential and interatomic forces generated by a deep neural network trained with data from ab initio calculations. We have discovered a number of soft modes which can exhibit dipolar, quadrupolar, and octupolar vibrational patterns in real space, as well as some time-reversal breaking chiral phonon modes. We have further studied the phonon effects on the electronic structures by freezing certain soft phonon modes. We find that if a soft quadrupolar phonon mode is assumed to be frozen, the system would exhibit a charge order which is perfectly consistent with recent experiments. Moreover, once some low-frequency C_2z-breaking modes get frozen, the Dirac points at the charge neutrality point would be gapped out, which provides an alternative perspective to the origin of correlated insulator state at charge neutrality point.

Reentrant Correlated Insulators in Twisted Bilayer Graphene at 25 T (2π Flux)

Twisted bilayer graphene (TBG) is remarkable for its topological flat bands, which drive strongly interacting physics at integer fillings, and its simple theoretical description facilitated by the Bistritzer-MacDonald Hamiltonian, a continuum model coupling two Dirac fermions. Because of the large moiré unit cell, TBG offers the unprecedented opportunity to observe reentrant Hofstadter phases in laboratory-strength magnetic fields near 25 T. This Letter is devoted to magic angle TBG at 2π flux where the magnetic translation group commutes. We use a newly developed gauge-invariant formalism to determine the exact single-particle band structure and topology. We find that the characteristic TBG flat bands reemerge at 2π flux, but, due to the magnetic field breaking C_{2z}T, they split and acquire Chern number ±1. We show that reentrant correlated insulating states appear at 2π flux driven by the Coulomb interaction at integer fillings, and we predict the characteristic Landau fans from their excitation spectrum.