## Strongly correlated Hofstadter subbands in minimally twisted bilayer graphene

Cheng Shen 1 ,a † , Yifei Guan 1 † , Davide Pizzirani , Zekang Zhou , Punam Barman , Kenji 2 1 1

Watanabe , Takashi Taniguchi , Steffen Wiedmann , Oleg V. Yazyev , and Mitali Banerjee 3 4 2 1 1 ⋆

1 Institute of Physics, Ecole Polytechnique F´d´rale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland e e

2 High Field Magnet Laboratory (HFML-EMFL), Radboud University,

Toernooiveld 7, 6525 ED Nijmegen, The Netherlands

3 Research Center for Functional Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan 4 International Center for Materials Nanoarchitectonics,

National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan

a

Current address:

School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China

† these authors contribute equally to the work and

⋆ corresponding author: mitali.banerjee@epfl.ch

Moir´ e superlattice in twisted bilayer graphene has been proven to be a versatile platform for exploring exotic quantum phases. Extensive investigations have been invoked focusing on the zeromagnetic-field phase diagram at the magic twist angle around θ = 1 1 , which has been indicated to . ° be an exclusive regime for exhibiting flat band with the interplay of strong electronic correlation and untrivial topology in the experiment so far. In contrast, electronic bands in non-magic-angle twisted bilayer graphene host dominant electronic kinetic energy compared to Coulomb interaction. By quenching the kinetic energy and enhancing Coulomb exchange interactions by means of an applied perpendicular magnetic field, here we unveil gapped flat Hofstadter subbands at large magnetic flux that yield correlated insulating states in minimally twisted bilayer graphene at θ = 0 41 . . ° These states appear with isospin symmetry breaking due to strong Coulomb interactions. Our work provides a platform for studying the phase transition of the strongly correlated Hofstadter spectrum.

## INTRODUCTION

Electrons in solids undergo a phase transition from conducting to a localized state when on-site Coulomb interaction dominates over the kinetic hopping energy. Such a phase transition, typically described by the MottHubbard model, requires a flat electronic band with a quenched kinetic energy. Recent progress in moir´ sue perlattice has revealed twisted bilayer graphene (TBG) as a platform system, in which a flat electronic band can arise from periodic interlayer coupling, and a comparable Coulomb interaction energy is approachable by tuning moir´ e wavelength through the tunable twist angle [1, 2]. Diverse quantum states with strong electronic correlation have been revealed in this system including correlated insulator [1], superconductivity [2], orbital magnet [3-6] and fractional Chern insulator [7]. Manipulating electronic flat bands in moir´ e superlattice provides the alluring potential to approach exotic quantum phases in condensed matters. To date, most work on TBG has focused on the so-called 'magic angle' of approximately 1 1 , . ° where flat bands are intrinsically emerging. With a fairly more dispersive band, the non-magic angle TBG typically requires external tuning to enable strong electronic correlation [8, 9].

Subjected to an out-of-plane magnetic field, the magnetic length l B provides the two-dimensional electron gas with tunable strength of Coulomb interactions. For electrons in a periodic lattice, the interplay between magnetic cyclotron energy and the periodic lattice potential gives rise to a non-interacting fractal Hofstadter butterfly spectrum with a recursive magnetic subband structure. The Coulomb interaction is intuitively ex- pected to affect such Hofstadter butterfly structure, yet not fully conclusive [10-13]. Adding an interacting Hubbard energy to Hofstadter Hamiltonian, the socalled Hofstadter-Hubbard model, is constructed to understand the correlated Hofstadter electrons. Experimentally, previous investigations have exhibited a noninteracting fractal spectrum by magnetoresistance oscillations in the graphene/hexagonal boron nitride (hBN) moir´ e heterostructure [14-16], and now extended to a strongly correlated regime in twisted graphene superlattice. For the latter case, the inherited strong Coulomb interaction from zero magnetic field results in the isospinsymmetry-broken Hofstadter subband ferromagnetism and translation-symmetry broken Chern insulator phases under small magnetic flux [17, 18]. Albeit the Hofstadter butterfly picture is generic to all twist angles, studies of correlated Hofstadter subbands in small twist angle TBG samples and the impact of Coulomb interaction associated with magnetic field are scarce. Here, we show the phase transition of Hofstadter subbands in a 0 41 . ° TBG at high magnetic flux. The magnetic-field-induced strong Coulomb interactions, ideal isolation, and flattening of Hofstadter subbands produce strong electronic correlation to induce various isospin-symmetry-broken correlated insulating states at partial fillings of moir´ unit e cell. Our results establish Hofstadter subbands in minimally twisted TBG as a platform for studying electronic correlations that open new opportunities as compared to the magic-angle TBG.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-e9c05cd5-8177-4a57-95a6-103ecab1c162.jpg)

R xx

FIG. 1. Hofstadter butterfly and band gap in small-angle twisted bilayer graphene. a , Schematics of moir´ pattern e and device structure. b , Dual gate dependence of four-probe resistance. c , Quantum oscillations in magnetic field characterized by longitudinal magneto resistance R xx (top panel) and Hall resistance R xy (middle panel). The bottom panel illustrates band insulators (red) and orbital Landau levels (black). The right y-axis shows the corresponding magnetic flux. d , Non-interacting numeric simulation of Hofstadter subband structure for 0 41 . ° TBG. Single-particle band gaps are illustrated with the Chern number C = 0. The Hofstadter subbands are marked by blue shading. e , Metal-insulator phase transition mediated by the perpendicular magnetic field. The lincuts are acquired from ( d ) (top panel) at moir´ unit cell filling e ν =-4 (blue), ν =0 (red), and ν =4 (yellow).

## GAPPED HOFSTADTER SUBBANDS

We use a standard dry van der Waals transfer technique to fabricate our stacking, which is composed of dual graphite gates, twisted graphene bilayers, and hexagonal boron nitride dielectric layers. We produced our TBG stackings with a targeting twist angle less than 1 , ° which is prone to reduction during the assembling process. In this study, we focus on TBG with 0 41 . ° twist angle, which is close to the predicted 'second magic angle' (SMA) (depending on the parameters of lattice reconstruction and interlayer hopping, SMA vary in the interval of 0 2 . ° -0 5 )[19, 20]. . ° We carried out four-terminal conductance measurements at a fridge temperature of around 8mK unless specified. SMA is speculated to host electronic correlation physics resembling that at the 1 1 . ° magic angle. However, our results show quite distinct behaviors at these two magic angles. At 1 1 . ° twist angle, TBG exhibits single-particle band gaps at one full filling of the moir´ unit cell (corresponding to four electrons or e holes in one moir´ unit cell) and thus flat bands isolated e from remote dispersive bands. Unlike 1 1 . ° TBG, our 0 41 . °

system features various metallic resistance peaks over a large doping regime in four-terminal transport measurements (Fig.1b). These peaks correspond to multiple full fillings of the moir´ unit cell but without gap openings. e Our non-interacting numeric simulations show multiple energy bands, despite a small energy scale of the order of 10 meV, overlapping with each other (Fig.S1). As a result, the effective bandwidth is enlarged to the extent that the correlated insulator phase at partial fillings is no longer favored. The dual gate measurements show a behavior independent of the displacement field, consistent with the strong interlayer coupling in TBG with a small twist angle(Fig.1b).

Our device exhibits a clear Hofstadter spectrum featuring intercepted cyclotron Landau levels from different fillings at rational magnetic flux Φ = Φ 0 /q , where q is integers and Φ 0 = h/e is the quantum magnetic flux. The well-determined magnetic flux helps us to reconfirm the twist angle via magnetic field B 0 = 4T for one quantum magnetic flux. Distinct from hetero-bilayer superlattices like graphene/hBN moir´ e structure, where moir´ wavee length is limited within a specific range due to the lat-

FIG. 2. Correlated insulators and phase transition in Hofstadter subbands. a , Longitudinal magneto resistance R xx as a function of charge filling ν and magnetic field B ⊥ . b , Schematic illustration of states in ( a ). Correlated insulators at ν =-3, -2 and -1, band insulators at ν =-4, 0 and 4, and orbital Landau levels are represented by blue, red, and black shadings, respectively. c , Linecuts of R xx and R xy as a function of charge filling ν . d , Emergent half-filling correlated insulator on electron side in a high magnetic field.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-9742db19-2e4a-4d93-b20f-9aa2153217da.jpg)

tice mismatch, twisted homo-bilayer can host extremely large moir´ wavelength at small twist angles. e Therefore, twisted homo-bilayer has the advantage of reaching high magnetic flux per moir´ e unit cell with lower magnetic fields that are easily available in the laboratory. Fig.1c show clear signs of two flux quanta where Landau levels prominently cross with each other. The resetting of Landau levels at rational magnetic quantum flux, shown by approximately zero Hall resistance R xy , originates from a renormalized zero effective magnetic field.

splitting, and Coulomb interaction when subjected to the perpendicular magnetic field. Like in graphene, the fieldinduced gap at CNP is dominated by the Coulomb interaction, featuring continuous increment in a magnetic field and hosting complicated ground states of charge density wave or intervalley coherence[24, 25]. This may interpret the inconsistency of field-dependence of the CNP gap between our theoretical calculation and experimental results.

Strikingly, metallic states of various full fillings at zero field develop into insulating states at the finite magnetic field, as characterized by quite high longitudinal resistance R xx and sign change of Hall resistance R xy over doping (Fig.1d, see also Fig.S3). Respectively, charge neutrality point (CNP) shows immediate gap opening in the field, while full fillings of ν = n/n 0 = 4 8 require , magnetic flux Φ = Φ 0 (here n is charge density and n 0 indicates charge density when moir´ unit cell is filled with e one electron). In Fig.1d, our non-interacting calculations of the Hofstadter spectrum through the Lanczos recursion method [21-23] predicts such single-particle band gaps in the magnetic field with Chern number C = 0 (see more calculations in Fig.S2), however, with a contrasting behavior of nearly closing when Φ &gt; 3Φ 0 at ν = 0 compared to experimental results. The low-energy band, considered roughly as renormalized mini Dirac cones, is associated with orbital cyclotron energy, spin/orbital Zeeman

We also note these gapped states are modulated by the crossing of the electron and hole cyclotron Landau levels at rational flux, compatible with our zero-field overlapped band structure and thus implying a candidate ground state of excitonic insulator at finite magnetic field. We note previous studies proposed a fragile topology derived from identical helicity of Dirac bands in the moir´ e Brillouin zone to interpret Landau level crossing with band gaps at CNP and full fillings [26, 27]. The exact mechanism underlying this crossing is beyond the scope of this work.

## CORRELATED INSULATORS IN HOFSTADTER SUBBANDS

Upon approaching Φ = 3Φ 0 we found that the quantum oscillations corresponding to orbital Landau levels embedded between band insulators at ν = 0 , ± 4 disap-

FIG. 3. Spin texture of correlated insulator at ν = 2 . a , Temperature dependence of longitudinal resistance R xx as a function of charge filling ν at perpendicular magnetic field B ⊥ = 11 7T. . b , In-plane magnetic field dependence of correlated insulator resistance at 1.5K and thermal-activation gap ∆ at ν = 2. The left panel shows a rotated sample measurement scheme. The total magnetic field B T is projected to in-plane B ∥ and out-of-plane B ⊥ component. B ⊥ is fixed at 12T. The inset figures show the spin texture transition from an antiferromagnetic to a canted antiferromagnetic state.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-d17f64a0-42bf-40d9-9ebd-8a615d04c529.jpg)

pear. When the moir´ e unit cell is partially filled with holes, these quantum oscillations give way to resistive states, which show no slopes of dn/dB ⊥ in the fan diagram (Fig.2a). Specifically, we observe these resistive states at half filling ( ν = -2, two holes per moir´ unit e cell) and quarter fillings ( ν = -1 , -3, one and three holes per moir´ e unit cell, respectively). These states emerge with a hierarchy of required magnetic fields. Resistive states also appear at fractional but non-fixed ν , which might be associated with changed singularities (like van Hove singularities) in Hofstadter subbands modulated by magnetic field. On the electron side, no state is found at partial fillings below B ⊥ = 14T, indicating a considerable particle-hole symmetry breaking in the electronic band structure, which is quite common in twisted and non-twisted crystalline multilayer graphene [28-32]. Increasing the field to Φ &gt; 4Φ , 0 we find a half-filling state at the electron side eventually appears (Fig.2d) (the high-resistance background is due to the resistive state of ν = 0 at high magnetic field which yields a bad contact).

To further elucidate these states at partial fillings, the Hall resistance R xy exhibiting sign change at integer filling ν = 2 as shown in Fig.2c, points to Fermi surface topology transition with gap opening. We note the Hall resistance R xy is obviously shifted by colossal longitu- dinal R xx component since the misalignment of opposite Hall bars is inevitable. Such an insulating state is further evidenced by the electronic thermal activation behavior of increasing resistance with decreasing temperature at ν = 2 (Fig.3a).

The states at partial fillings ν = -1 , ± 2 , -3 of Hofstadter subbands unambiguously indicate SU(4) spinvalley symmetry breaking induced by strong electronic correlation. To understand the underlying mechanism for strong electronic correlation in Hofstadter subbands, we consider the typical Stoner ferromagnetism model involving competing Coulomb exchange interaction U and density of state at Fermi level D E ( F ). In the Stoner picture, when D E ( F ) · U &gt; 1 is satisfied, the isospin flavor symmetry will be spontaneously broken, forming spin or valley-polarized correlated states. The magnetic field leads to high U and D E ( F ) to satisfy Stoner criteria, as we will address in the following.

An out-of-plane magnetic field B ⊥ substantially reconstructs the zero-field low-energy moir´ e bands. In our case, band reconstruction includes: (1) Hofstadter subbands isolation by gap openings at ν = 0 , ± 4; (2) reduced bandwidth of Hofstadter subbands. The crucial role of Hofstadter subband isolation for correlated insulators is reflected by the hierarchy that partial fillings ν = -1 , ± 2 , -3 require preformed gap states at ν = 0 , ± 4. The bandwidth of Hofstadter subbands is directly related to orbital Landau level gaps inside itself. In Fig.2a, the disappearance of quantum oscillation (illustrated as black lines in Fig.2b) in high field signifies reduced orbital Landau level gaps, as Landau level broadening is exclusively induced by the fixed disorder and temperature effect here. Thus, this phase transition directly points to a related flattening of Hofstadter subbands where high D E ( F ) is available. In Fig.S1b, the calculated Hofstadter subbands at 3Φ 0 show a narrow bandwidth smaller than 10meV. We roughly estimate the on-site Coulomb interaction energy in 0.41 ° TBG at the zero-field limit compared to magic angle [33] by only considering the change of moir´ wavelength. It yields one site Coulomb interaction energy of the same order as the calculated bandwidths, implying a correlated regime for Hofstadter electrons. Despite the real bandwidth or density of the state is hard to precisely acquire theoretically (due to the strong lattice relaxation effect for small angle TBG), we also emphasize the impact of Coulomb exchange interaction U between different isospins. Since U is proportional to e /l 2 B where l B = √ ℏ /eB ⊥ is the magnetic length, the Stoner criteria can also be satisfied with a much enhanced U in a high magnetic field. Observing correlated insulators of partial fillings in higher magnetic flux is consistent with the above analysis. In contrast, in magic-angle TBG, isospin symmetry breaking for Hofstadter subbands is realized at a low magnetic field due to inherited strong electronic correlation intrinsic to zero fields [24]. Hofstadter subbands flattening and enhanced Coulomb exchange interaction by magnetic flux, therefore providing an avenue for realizing strong correlations

FIG. 4. Isospin-symmetry-broken correlated insulators with nonzero Chern number. a , Mapping plot of longitudinal magnetoresistance R xx and Hall resistance R xy as a function of perpendicular magnetic field B ⊥ and charge filling ν . b , Linecuts of R xx and R xy acquired from ( a ) with respect to renormalized charge filling s ⋆ . Here s ⋆ = s and s ⋆ = s -4 for B ⊥ = 13 11T , and B ⊥ = 7T, respectively. s follows the Diophantine equation ν = C Φ Φ / 0 + . s The R xx linecuts are consecutively shifted upward by 1k Ω . The dashed line in the bottom figure represents the Hall resistance of 4 h/e 2 . c , Schematic illustration of states in ( a ). Band gaps with Chern number C = 0, C = 4, and correlated gaps are illustrated by red, black, and blue shadings, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-fa697a27-2607-4771-b50c-a628566da30e.jpg)

.

in a system that doesn't inherently host such a scenario.

## SPIN FLAVOR ORDERING AT HALF FILLING OF HOFSTADTER SUBBANDS

The isospin symmetry breaking leads to intricate isospin flavor ordering. Subjected to the competing spin Zeeman effect and Coulomb exchange interaction in a perpendicular magnetic field, the half-filling correlated insulator can be either a spin-unpolarized state due to an antiferromagnetic exchange interaction or a spinpolarized state due to the spin alignment by a magnetic field. We apply a tilted magnetic field B T , whose out-ofplane projection B ⊥ is fixed and in-plane B ∥ component is varied. As the in-plane magnetic field can only alter the electron's spin, not its orbit, the response to the in-plane component of the tilted magnetic field can thus help to unveil the spin flavor ordering of the ground state. Fig.3b shows the four-probe resistance and thermal-activation gap (fitted with Arrhenius law) at half filling. They exhibit non-monotonic behavior with respect to the inplane field, showing a minimum at B ∥ around 2 T to 4 T. The decreasing of the correlated gap by a small inplane field reveals an in-plane spin-unpolarized ground state with antiferromagnetic (AF) ordering at B ⊥ = B T .

The larger in-plane field alters the spin inclined to inplane ferromagnetic ordering, forming a canted antiferromagnetic (CAF) state with respect to the direction of the tilted field (as shown in Fig.3b). The canted angle θ is determined by the ratio of Zeeman energy gµ B B T and the antiferromagnetic exchange interaction, which depends only on B ⊥ (here, g is bare gyromagnetic ratio, µ B is Bohr magneton).

## CORRELATED GAPS WITH NON-ZERO CHERN NUMBER

Spanning from ν = 0 to ν = -4, all the insulator states exhibit zero Chern number as a characteristic of zero dn/dB in the phase diagram shown in Fig.2a according to the Streda formula C = h e · dn dB . In contrast to a common expectation that magnetic field endows Hofstadter subband with nonzero Chern numbers[17, 34], here the Hofstadter subband, as embedded between gaps with zero Chern numbers at ν = 0 and ν = -4, has a net Chern number ∆ C = C ν ( = 0) -C ν ( = -4) = 0. For correlated gap at partial fillings ν , its Chern number is C ν ( ) = C ν ( = -4) + ν ∆ C/ 4, which turns out to be always zero. In Fig.4, we depict symmetry-broken correlated states that correspond to filling the moir´ unit cell e with one, two, or three electrons in a regime of magnetic flux Φ &gt; 2Φ 0 and two electrons for Φ &lt; 2Φ . 0 These states show more pronounced symmetry breaking in higher magnetic flux, which is compatible with the above picture of band flattening and enhanced Coulomb interaction by a magnetic field. However, they all have

- [1] Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, and P. JarilloHerrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 556 , 80 (2018).
- [2] Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, Unconventional superconductivity in magic-angle graphene superlattices, Nature 556 , 43 (2018).
- [3] A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watanabe, T. Taniguchi, M. A. Kastner, and D. Goldhaber-Gordon, Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene, Science 365 , 605 (2019).
- [4] X. Lu, P. Stepanov, W. Yang, M. Xie, M. A. Aamir, I. Das, C. Urgell, K. Watanabe, T. Taniguchi, G. Zhang, A. Bachtold, A. H. MacDonald, and D. K. Efetov, Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene, Nature 574 , 653 (2019).
- [5] M. Serlin, C. L. Tschirhart, H. Polshyn, Y. Zhang, J. Zhu, K. Watanabe, T. Taniguchi, L. Balents, and A. F. Young, Intrinsic quantized anomalous Hall effect in a moir´ hete erostructure, Science 367 , 900 (2020).

a Chern number C = 4 with edge states as characterized by the zero longitudinal magnetoresistance R xx and Hall resistance R xy plateau of approximately 4 h/e 2 (Fig.4b). Described by Diophantine equation ν = C Φ Φ + , these / 0 s states are characterized by ( C, s ) = (4 , 1) (4 , , 2)(4 3) and , (4 6), respectively. , Having the same origin as correlated insulators in the first valence Hofstadter subbands, these correlated states appear by splitting higher-energy flat conduction Hofstadter subbands via the isospin symmetry breaking. The net Chern number of these Hofstadter subbands is ∆ C = 0 as well, as determined by the same Chern number C = 4 for the two adjacent orbital Landau level gaps in the Hofstadter spectrum. Splitting these Hofstadter subbands yields correlated gaps with Chern number C = 4.

̸

To conclude, the flat Hofstadter subbands in smallangle twisted bilayer graphene that are universal in various magnetic flux and charge fillings provide a new pathway for investigating strongly correlated electronic phases. The phase transition associated with isospin symmetry breaking observed in our work lays the foundation to study correlated Hofstadter spectrum and achieve intriguing fractional Chern insulating states in Hofstadter subbands with nonzero net Chern number∆ C = 0 if possible [7, 34]. The intricate ground states and isospin flavor ordering for correlated insulators at partial fillings and even band insulators in Hofstadter subbands appeal to more theoretical efforts. Further experimental investigations through applying higher magnetic flux, local compressibility measurements[7], etc., will help to shed light on these exotic quantum phases.

- [6] C. L. Tschirhart, M. Serlin, H. Polshyn, A. Shragai, Z. Xia, J. Zhu, Y. Zhang, K. Watanabe, T. Taniguchi, M. E. Huber, and A. F. Young, Imaging orbital ferromagnetism in a moir´ Chern insulator, Science e 372 , 1323 (2021).
- [7] Y. Xie, A. T. Pierce, J. M. Park, D. E. Parker, E. Khalaf, P. Ledwith, Y. Cao, S. H. Lee, S. Chen, P. R. Forrester, K. Watanabe, T. Taniguchi, A. Vishwanath, P. JarilloHerrero, and A. Yacoby, Fractional Chern insulators in magic-angle twisted bilayer graphene, Nature 600 , 439 (2021).
- [8] M. Yankowitz, S. Chen, H. Polshyn, Y. Zhang, K. Watanabe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, Tuning superconductivity in twisted bilayer graphene, Science 363 , 1059 (2019).
- [9] H. S. Arora, R. Polski, Y. Zhang, A. Thomson, Y. Choi, H. Kim, Z. Lin, I. Z. Wilson, X. Xu, J.-H. Chu, K. Watanabe, T. Taniguchi, J. Alicea, and S. Nadj-Perge, Superconductivity in metallic twisted bilayer graphene stabilized by WSe2, Nature 583 , 379 (2020).
- [10] K. Czajka, A. Gorczyca, M. M. Ma´ ska, and M. Mierzejewski, Hofstadter butterfly for a finite correlated system, Physical Review B 74 , 125116 (2006).
- [11] A. Barelli, J. Bellissard, P. Jacquod, and D. L. Shep-

- elyansky, Double Butterfly Spectrum for Two Interacting Particles in the Harper Model, Physical Review Letters 77 , 4752 (1996).
- [12] H. Doh and S.-H. S. Salk, Effects of electron correlations on the Hofstadter spectrum, Physical Review B 57 , 1312 (1998).
- [13] J. K. Ding, W. O. Wang, B. Moritz, Y. Schattner, E. W. Huang, and T. P. Devereaux, Thermodynamics of correlated electrons in a magnetic field, Communications Physics 5 , 204 (2022).
- [14] B. Hunt, J. D. Sanchez-Yamagishi, A. F. Young, M. Yankowitz, B. J. LeRoy, K. Watanabe, T. Taniguchi, P. Moon, M. Koshino, and P. Jarillo-Herrero, Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure, Science 340 , 1427 (2013).
- [15] L. A. Ponomarenko, R. V. Gorbachev, G. L. Yu, D. C. Elias, R. Jalil, A. A. Patel, A. Mishchenko, A. S. Mayorov, C. R. Woods, J. R. Wallbank, M. MuchaKruczynski, B. A. Piot, M. Potemski, I. V. Grigorieva, K. S. Novoselov, F. Guinea, V. I. Fal'ko, and A. K. Geim, Cloning of Dirac fermions in graphene superlattices, Nature 497 , 594 (2013).
- [16] C. R. Dean, L. Wang, P. Maher, C. Forsythe, F. Ghahari, Y. Gao, J. Katoch, M. Ishigami, P. Moon, M. Koshino, T. Taniguchi, K. Watanabe, K. L. Shepard, J. Hone, and P. Kim, Hofstadter's butterfly and the fractal quantum Hall effect in moir´ superlattices, Nature e 497 , 598 (2013).
- [17] Y. Saito, J. Ge, L. Rademaker, K. Watanabe, T. Taniguchi, D. A. Abanin, and A. F. Young, Hofstadter subband ferromagnetism and symmetry-broken Chern insulators in twisted bilayer graphene, Nature Physics 17 , 478 (2021).
- [18] J. Yu, B. A. Foutty, Z. Han, M. E. Barber, Y. Schattner, K. Watanabe, T. Taniguchi, P. Phillips, Z.-X. Shen, S. A. Kivelson, and B. E. Feldman, Correlated Hofstadter spectrum and flavour phase diagram in magic-angle twisted bilayer graphene, Nature Physics 18 , 825 (2022).
- [19] Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, All Magic Angles in Twisted Bilayer Graphene are Topological, Physical Review Letters 123 , 36401 (2019).
- [20] G. Tarnopolsky, A. J. Kruchkov, and A. Vishwanath, Origin of Magic Angles in Twisted Bilayer Graphene, Physical Review Letters 122 , 106405 (2019).
- [21] Q. Wu, J. Liu, Y. Guan, and O. V. Yazyev, Landau Levels as a Probe for Band Topology in Graphene Moir \ 'e Superlattices, Physical Review Letters 126 , 56401 (2021).
- [22] Y. Guan, O. V. Yazyev, and A. Kruchkov, Reentrant magic-angle phenomena in twisted bilayer graphene in integer magnetic fluxes, Physical Review B 106 , L121115 (2022).
- [23] Q. Wu, S. Zhang, H.-F. Song, M. Troyer, and A. A. Soluyanov, WannierTools: An open-source software package for novel topological materials, Computer Physics Communications 224 , 405 (2018).
- [24] A. F. Young, C. R. Dean, L. Wang, H. Ren, P. CaddenZimansky, K. Watanabe, T. Taniguchi, J. Hone, K. L. Shepard, and P. Kim, Spin and valley quantum Hall ferromagnetism in graphene, Nature Physics 8 , 550 (2012).
- [25] X. Liu, G. Farahi, C.-L. Chiu, Z. Papic, K. Watanabe, T. Taniguchi, M. P. Zaletel, and A. Yazdani, Visualizing broken symmetry and topological defects in a quantum Hall ferromagnet, Science 375 , 321 (2022).
- [26] B. Lian, F. Xie, and B. A. Bernevig, Landau level of
- fragile topology, Physical Review B 102 , 041402 (2020).
- [27] X. Lu, B. Lian, G. Chaudhary, B. A. Piot, G. Romagnoli, K. Watanabe, T. Taniguchi, M. Poggio, A. H. MacDonald, B. A. Bernevig, and D. K. Efetov, Multiple flat bands and topological Hofstadter butterfly in twisted bilayer graphene close to the second magic angle, Proceedings of the National Academy of Sciences 118 , e2100006118 (2021).
- [28] G. Chen, L. Jiang, S. Wu, B. Lyu, H. Li, B. L. Chittari, K. Watanabe, T. Taniguchi, Z. Shi, J. Jung, Y. Zhang, and F. Wang, Evidence of a gate-tunable Mott insulator in a trilayer graphene moir´ superlattice, Nature Physics e 15 , 237 (2019).
- [29] C. Shen, Y. Chu, Q. Wu, N. Li, S. Wang, Y. Zhao, J. Tang, J. Liu, J. Tian, and K. Watanabe, Correlated states in twisted double bilayer graphene, Nature Physics 16 , 520 (2020).
- [30] X. Liu, Z. Hao, E. Khalaf, J. Y. Lee, Y. Ronen, H. Yoo, D. H. Najafabadi, K. Watanabe, T. Taniguchi, and A. Vishwanath, Tunable spin-polarized correlated states in twisted double bilayer graphene, Nature 583 , 221 (2020).
- [31] Y. Cao, D. Rodan-Legrain, O. Rubies-Bigorda, J. M. Park, K. Watanabe, T. Taniguchi, and P. Jarillo-Herrero, Tunable correlated states and spin-polarized phases in twisted bilayer-bilayer graphene, Nature 583 , 215 (2020).
- [32] S. Chen, M. He, Y.-H. Zhang, V. Hsieh, Z. Fei, K. Watanabe, T. Taniguchi, D. H. Cobden, X. Xu, C. R. Dean, and M. Yankowitz, Electrically tunable correlated and topological states in twisted monolayer-bilayer graphene, Nature Physics 17 , 374 (2021).
- [33] D. Wong, K. P. Nuckolls, M. Oh, B. Lian, Y. Xie, S. Jeon, K. Watanabe, T. Taniguchi, B. A. Bernevig, and A. Yazdani, Cascade of electronic transitions in magic-angle twisted bilayer graphene, Nature 582 , 198 (2020).
- [34] E. M. Spanton, A. A. Zibrov, H. Zhou, T. Taniguchi, K. Watanabe, M. P. Zaletel, and A. F. Young, Observation of fractional Chern insulators in a van der Waals heterostructure, Science 360 , 62 (2018).
- [35] F. Gargiulo and O. V. Yazyev, Structural and electronic transformation in low-angle twisted bilayer graphene, 2D Materials 5 , 015019 (2018).
- [36] W. Zhu, T. Low, V. Perebeinos, A. A. Bol, Y. Zhu, H. Yan, J. Tersoff, and P. Avouris, Structure and Electronic Transport in Graphene Wrinkles, Nano Letters 12 , 3431 (2012).
- [37] J. C. Slater and G. F. Koster, Simplified LCAO Method for the Periodic Potential Problem, Physical Review 94 , 1498 (1954).
- [38] N. N. T. Nam and M. Koshino, Lattice relaxation and energy band modulation in twisted bilayer graphene, Physical Review B 96 , 075311 (2017).
- [39] N. Nakatsuji, T. Kawakami, and M. Koshino, Multiscale lattice relaxation in general twisted trilayer graphenes, Physical Review X 13 , 041007 (2023).
- [40] S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of Computational Physics 117 , 1 (1995).
- [41] J. H. Los and A. Fasolino, Intrinsic long-range bond-order potential for carbon: Performance in monte carlo simulations of graphitization, Physical Review B 68 , 024107 (2003).
- [42] A. N. Kolmogorov and V. H. Crespi, Registry-dependent

interlayer potential for graphitic systems, Physical Review B 71 , 235415 (2005).

## ACKNOWLEDGMENTS

We thank Ady Stern for the fruitful discussions. C.S. acknowledges the cleanroom facilities at Cmi (EPFL). We acknowledge the support of the HFML, member of the European Magnetic Field Laboratory (EMFL). M.B. acknowledges the support of SNSF Eccellenza grant No. PCEGP2 194528, and support from the QuantERA II Programme that has received funding from the European Union's Horizon 2020 research and innovation program under Grant Agreement No 101017733. Y.G. and O.V.Y. acknowledge support from the Swiss National Science Foundation (grant No. 204254). Computations were performed at the Swiss National Supercomputing Centre (CSCS) under project No. s1146 and the facilities of the Scientific IT and Application Support Center of EPFL. K.W. and T.T. acknowledge support from the JSPS KAKENHI (Grant Numbers 20H00354 and 23H02052) and World Premier International Research Center Initiative (WPI), MEXT, Japan.

## Author Contributions

C.S. fabricated devices and performed measurements ( &lt; 14T) and data analysis. Y.G. and O.V.Y. carried out numeric simulations. C.S.performed the high-field measurements with the help from D.P., Z.Z., P.B. and S.W.. T.T. and K.W. provided hBN crystals. M.B. supervised this project. C.S. wrote the paper with input from Y.G., O.V.Y. and M.B..

## Competing interests

The authors declare no competing interests.

## Data availability

The data supporting the findings of this study are available from the corresponding author upon reasonable request.

## SUPPLEMENTARY INFORMATION

## A. Device fabrication and transport measurement

We made stackings using the van der Waals assembly technique. Graphite, hexagonal boron nitride, and graphene flakes were exfoliated onto SiO 2 (285nm)/Si++ substrate and picked up at 100 C with a propylene carbonate (PC) ◦ film which is placed on a polydimethyl siloxane (PDMS) stamp. We twisted graphene bilayers with a targeting twist angle less than 1 ◦ . Stackings are fabricated with a Hall bar configuration via e-beam lithography and SF 6 /O 2 plasma etching and contacted with Cr/Au metal. Measurements are done with lock-in amplifiers in Bluefors dilution fridge at a base temperature of 8mK unless specified.

The twist angle is identified through the magnetic flux quanta Φ . 0 Φ 0 and the corresponding magnetic field B 0 are acquired from the Hofstadter spectrum. The twist angle is therefore obtained via θ = √ √ 3 B a / 0 2 (2Φ ) (here 0 a = 0 246 . nm is graphene lattice constant).

## B. Band-structure calculations

## 1. Tight-binding Hamiltonian

The electronic structure calculations are performed with the tight-binding Hamiltonian accounting for the atomistic details of twisted bilayer graphene in free-electron formalism,

$$H = \sum _ { i., j } t _ { i j } c _ { i } ^ { \dagger } c _ { j }$$

where c i ( c † i ) are the annihilation (creation) operators of the p z -orbital electron states on the i -th carbon atom. The hopping integral t ij between two p z orbitals is expressed as a sum of σ -type and π -type couplings following the Slater-Koster formalism [35-37]

$$t _ { i j } = V _ { 0 } ^ { \pi } \exp ( - \frac { r - a _ { 0 } } { r _ { 0 } } ) \sin ^ { 2 } \theta + V _ { 0 } ^ { \sigma } e x p ( - \frac { r - d _ { 0 } } { r _ { 0 } } ) \cos ^ { 2 } \theta,$$

where r = | r ij | and θ is the angle between r ij and e z , θ = ∠ ( r ij , (0 , 0 1)). , The numerical parameters have been adapted from Ref.[38]; hence V π 0 = -2 7 eV is the hopping integral between nearest-neighbor (NN) atoms in the same . layer; V σ 0 = 0 48 eV is the hopping integral between the aligned atoms in different layers. . Here, a 0 = 0 142 nm and . d 0 = 0 334 nm correspond to the nearest neighbor distance and the interlayer distance; . r 0 = 0 148 . a (we use a = √ 3 a 0 as the lattice constant of graphene) is the hopping strength decay length. The cutoff in real space is set to 9.84 ˚ , A outside which the coupling vanishes. It is known that the relaxation effect has a significant effect on the band structures of TBG [38, 39]. We use the classical force fields package LAMMPS [40] to relax the atomic structure. The classical force field contains the bond-order potential [41] as well as the modified version of the Kolmogorov-Crespi registrydependent potential for describing the interlayer interactions [42]. The Kolmogorov-Crespi long range potential is fitted to DFT results [35].

## 2. Chern numbers in Hofstadter gaps

The magnetotransport properties of the system are obtained by Hofstadter spectrum calculations through the Lanczos recursion methods. We calculate the chiral edge states in the Hofstadter gaps to probe the quantized Hall conductance in the topological gaps. Defining the external magnetic field B = B e z , finite width is taken along the y -direction, while the x -direction is periodic. For example, in Figure S2 we present the band structure plot in such a configuration for Φ = 3Φ . 0 The localization of the eigenstates on the edge is represented with the color-coding. The number of edge states on each edge corresponds to the Chern number in the gap.

## C. Twist-angle homogeneity

To identify the twist-angle homogeneity, we measure different pairs of electrodes. Fig.S4 shows the four-probe resistance. It's obvious that the charge density corresponding to full fillings of moir´ unit cell doesn't vary much e among different regions, indicating high twist-angle homogeneity with variation less than 0.007 ° .

## D. More data for correlated insulators

Fig.S5 shows extended data acquired from the device with twist angle θ = 0 37 . ° (we examine the twist angle via the magnetic flux quanta at 3.6T). All the features including the band gaps and Hofstadter butterfly resemble to the device of θ = 0 41 . ° discussed in the main text. We observed the resistive state at ν = 2 when B ⊥ &gt; 10T.

FIG. S1. Non-interacting band structures of TBG in the moir´ Brillouin zone at 0.41 twist angle. e ° a , Band structure at Φ = 0. b , Band structure at Φ = 3Φ . 0

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-bdae1c50-d4d9-47f1-a480-c925bdde9427.jpg)

.

FIG. S2. Example of the edge state band structure plot at Φ = 3Φ . The color of the points encodes the localization 0 of the states on one of the edges. The dashed box highlights flat Hofstadter subbands.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-eea354ed-3a3e-40ca-ad31-baeda1ef0c31.jpg)

.

FIG. S3. Extended data for band gaps . Longitudinal magneto-resistance R xx ( a ) and Hall resistance R xy ( b ) as a function of charge filling ν at B ⊥ =2T, 5T and 9T. R xx and R xy for different perpendicular magnetic fields are sequentially shifted by 25k Ω and 40k Ω , respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-2093e98c-142f-4018-a694-f979b3898d23.jpg)

.

FIG. S4. Twist angle homogeneity . Four-probe resistance R xx data acquired from different pairs of electrodes at zero field. The inset figures show the optical and atomic force microscopic images of the device.The error bars for two figures are 2 µm (left) and 1 µm (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-c7612114-7ace-47c7-a51f-484796cc20e9.jpg)

.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-c61134ae-cdb6-4304-b809-eaaf0a36f2a7.jpg)

FIG. S5. Extended data acquired from another device of θ = 0 37 . . ° a , Dual-gate mapping of four-probe resistance R xx . b , Longitudinal magneto-resistance as a function of back-gate voltage and perpendicular magnetic field. Top-gate voltage is fixed at 0V. The dash lines mark the magnetic flux quanta Φ = Φ 0 at 3.4T and Φ = 2Φ . 0 The arrow points to the emergent half-filling correlated insulator. Both figures are obtained at 240mK.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-d6495656-4d71-4cec-9c14-39ab05bded47.jpg)

.

FIG. S6. Arrhenius fitting of thermal activation gap at ν = 2. The original temperature dependence of four-probe resistance is depicted in the form of logarithm as a function of 1 /T . The device was rotated by 0 , ° 10 , ° 20 , ° 30 and 40 , ° ° corresponding to B ∥ =0T, 2.12T, 4.37T, 6.93T and 10.07T, respectively when B ⊥ is fixed at 12T.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-961a4a4a-e49b-401b-87c0-c107d927d31a.jpg)

.

FIG. S7. Non-interacting tight binding calculation of higher-energy Hofsadter subbands . The Hofstadter subbands are depicted by blue shadings, adjacent to the C = 4 Chern gaps.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00868v1-figure-526c1983-47c3-499b-b33e-81569acf2dba.jpg)

.
Cheng Shen Yifei Guan Davide Pizzirani Zekang Zhou Punam Barman Kenji Watanabe Takashi Taniguchi Steffen Wiedmann Oleg V. Yazyev Mitali Banerjee
Strongly correlated Hofstadter subbands in minimally twisted bilayer graphene
Moir\'e superlattice in twisted bilayer graphene has been proven to be a
versatile platform for exploring exotic quantum phases. Extensive
investigations have been invoked focusing on the zero-magnetic-field phase
diagram at the magic twist angle around $\theta=1.1\degree$, which has been
indicated to be an exclusive regime for exhibiting flat band with the interplay
of strong electronic correlation and untrivial topology in the experiment so
far. In contrast, electronic bands in non-magic-angle twisted bilayer graphene
host dominant electronic kinetic energy compared to Coulomb interaction. By
quenching the kinetic energy and enhancing Coulomb exchange interactions by
means of an applied perpendicular magnetic field, here we unveil gapped flat
Hofstadter subbands at large magnetic flux that yield correlated insulating
states in minimally twisted bilayer graphene at $\theta=0.41\degree$. These
states appear with isospin symmetry breaking due to strong Coulomb
interactions. Our work provides a platform for studying the phase transition of
the strongly correlated Hofstadter spectrum.