Chinese Physics Letters, 2021, Vol. 38, No. 4, Article code 047301Express Letter Emergence of Chern Insulating States in Non-Magic Angle Twisted Bilayer Graphene Cheng Shen (沈成)1,2, Jianghua Ying (应江华)1,2, Le Liu (刘乐)1,2, Jianpeng Liu (刘健鹏)3,4, Na Li (李娜)1,2,6, Shuopei Wang (王硕培)1,2,6, Jian Tang (汤建)1,2, Yanchong Zhao (赵岩翀)1,2, Yanbang Chu (褚衍邦)1,2, Kenji Watanabe7, Takashi Taniguchi8, Rong Yang (杨蓉)1,5,6, Dongxia Shi (时东霞)1,2,5, Fanming Qu (屈凡明)1,2,6, Li Lu (吕力)1,2,6, Wei Yang (杨威)1,2,6*, and Guangyu Zhang (张广宇)1,2,5,6* Affiliations 1Beijing National Laboratory for Condensed Matter Physics; Key Laboratory for Nanoscale Physics and Devices, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 2School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China 3School of Physical Sciences and Technology, ShanghaiTech University, Shanghai 200031, China 4ShanghaiTech Laboratory for Topological Physics, ShanghaiTech University, Shanghai 200031, China 5Beijing Key Laboratory for Nanomaterials and Nanodevices, Beijing 100190, China 6Songshan-Lake Materials Laboratory, Dongguan 523808, China 7Research Center for Functional Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan 8International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan Received 11 February 2021; accepted 9 March 2021; published online 17 March 2021 *Corresponding authors. Email: wei.yang@iphy.ac.cn; gyzhang@iphy.ac.cn
Citation Text: Chen C, Ying J H, Liu L, Liu J P, and Li N et al. 2021 Chin. Phys. Lett. 38 047301    Abstract Twisting two layers into a magic angle (MA) of $\sim$$1.1^{\circ} is found essential to create low energy flat bands and the resulting correlated insulating, superconducting, and magnetic phases in twisted bilayer graphene (TBG). While most of previous works focus on revealing these emergent states in MA-TBG, a study of the twist angle dependence, which helps to map an evolution of these phases, is yet less explored. Here, we report a magneto-transport study on one non-magic angle TBG device, whose twist angle \theta changes from 1.25^{\circ} at one end to 1.43^{\circ} at the other. For \theta =1.25^{\circ} we observe an emergence of topological insulating states at hole side with a sequence of Chern number \left| C \right|=4-\left| v \right|, where v is the number of electrons (holes) in moiré unite cell. When \theta >1.25^{\circ}, the Chern insulator from flat band disappears and evolves into fractal Hofstadter butterfly quantum Hall insulator where magnetic flux in one moiré unite cell matters. Our observations will stimulate further theoretical and experimental investigations on the relationship between electron interactions and non-trivial band topology. DOI:10.1088/0256-307X/38/4/047301 © 2021 Chinese Physics Society Article Text Twist angle plays an important role in twisted bilayer graphene (TBG) and other two-dimensional twisted moiré systems. By twisting two layers into a magic angle (MA) of \sim$$1.1^{\circ}$ in TBG, moiré bands are generated with a narrow band width of $\sim$10 meV, so-called flat band where it favors electron interactions over kinetic energy,[1–9] contributing to the realizations of correlated insulating,[2] superconducting,[3–5] and orbital magnetic phases.[5,10–12] As the twist angle is increased further to non-magic angle (NMA), the moiré band width is increased very fast and the bands are no longer flat, where the electrons tend to lose correlation.[1,13,14] Recently, aside from the observation of non-zero Chern bands in TBG,[5,11,12] twisted multilayer graphene[15,16] and ABC-stacked trilayer graphene[17] where inversion symmetry ($C_{2}$) is broken, the strong electron interactions in MA-TBG are found able to break time reversal symmetry ($T$), untie the degeneracy between the flat bands, and thus reveal a sequence of non-trivial topological insulating states with Chern number $|C|=4-|v|$,[18–21] where $v=n/n_{0}$ is the number of electrons (or holes) filled in one moiré unit cell, $n$ is the carrier density, and $n_{0}$ is carrier density of one electron per moiré unit cell area. Then, intuitively one may ask how strong an electron interaction is needed to produce these Chern bands, and if it survives even for NMA-TBG or not. The question is interesting and important, and yet to be explored. In this Letter, we try to bridge the gap by studying magneto transport properties of an NMA-TBG device with different twist angles. A perpendicular magnetic field  ($B$)  would have impacts on  TBG  in  two different ways. First, the magnetic field tends to recombine the flat bands of TBG into a series of fractal Landau levels (LLs), i.e., the so-called Hofstadter butterfly spectra, in which the recurring fractal bands are dependent on the number of magnetic fluxes in each moiré primitive cell. Second, aside from the formation of LLs, the magnetic field also induces splitting between the flat bands with opposite Chern numbers $+$/$-1$ due to the orbital Zeeman effect.[22]  Such orbital Zeeman splitting originates from the intrinsic orbital magnetism of the topological flat bands of TBG, and can be dramatically enhanced by Coulomb interactions, which would give rise to a series of time-reversal broken moiré Chern bands without the necessity of forming LLs. Note that quantum Hall insulator is in principle a Chern insulator, where the LL filling factor corresponds to the Chern number. To avoid confusion, we use “Chern” specifically for moiré Chern bands in this Letter. These moiré Chern bands under (weak) perpendicular magnetic field $B$ are manifested by quantized Hall conductance $\sigma_{xy}=Ce^{2}/h$, where $e$ is the electron charge and $h$ is Planck's constant. The Chern number could also be obtained from a linear slop of the quantized conductance in Landau fan diagram by $C=(h/e)(dn/dB)$, where $n$ is the carrier density. Although the ratio between the characteristic Coulomb interaction strength $U$ and the moiré bandwidth $W$ is significantly reduced for NMA TBG, it is still of fundamental interest to reveal the interplay and competition between the fractal LL spectra and the topological moiré bands driven by (interaction-enhanced) orbital Zeeman splittings. We shed light on such an intriguing problem by measuring Landau fan diagrams in NMA TBG samples with different twist angles.
Fig. 1. Hofstadter butterfly spectra in TBG with $\theta =1.38^{\circ}$. (a) Schematic side view of TBG device structure. (b) Two-terminal conductance as a function of carrier density $n$ acquired among different Hall bar pairs and at temperature $T=1.7$ K. The scale bar in the inset is 4 µm. (c) Schematics of twist angle distribution. (d) Transverse Hall resistance $R_{xy}$ mapping plot versus electron or hole filling of moiré unit cell and perpendicular magnetic fields. (e) Schematic illustration of Wannier diagram in panel (d). (f) Line cuts at varied fields from panel (d), showing suppression of quantum Hall ferromagnetism near CNP.
A schematic illustration of the device is shown in Fig. 1(a). The twisted graphene layers are encapsulated by hBN, and they are not aligned with hBN substrate in order to hold the $C_{2}$ sublattice symmetry (Fig. S1 in Supplementary Information). The non-magic angle nature is revealed in the transfer curves at different positions as shown in Fig. 1(b), from which the extracted twist angle changes from 1.25$^{\circ}$ at one end to 1.43$^{\circ}$ at the other in Fig. 1(c), and correlated insulating states are absent in our device. The twist angles are further quantified and validated in the magneto-transport, regardless of a moderate inhomogeneity of $\pm 0.03^{\circ}$ for each device presented in the main text. The details for assignment of twist angle and related discussions of twist inhomogeneity are presented in the first section of Supplementary Information. We start with magneto transport studies at a twist angle $\theta \sim 1.38^{\circ}$, where fractal Hofstadter butterfly spectra dominate. The measurements are performed at 20 mK, and a color mapping of Hall resistance $R_{xy}$ as a function of carrier density and magnetic field is shown in Fig. 1(d). In the TBG moiré system, a single particle picture of fractal LL spectra can be derived from the interplay of periodic interlayer moiré potential and magnetic fields.[23] They are featured by the recurring embryo LLs which emanate from magnetic flux per moiré unit cell $\mathit{\Phi}=\mathit{\Phi}_{0}(p/q)$, where $p$ and $q$ are co-prime integers, $\mathit{\Phi}_{0}=h/e$ is the quantum magnetic flux. These embryo LLs are well understood as replica minibands at effective magnetic field $B_{\rm eff}=\pm \vert B-B_{p/q}\vert$, with spin-valley symmetry breaking driven by electron interactions.[24–28] Notably, as depicted in Fig. 1(e), one could identify LLs from charge neutral point (CNP) at $v = 0$ and that from full filling at $v = 4$, and the fractal LLs appear only in the regimes between two rational fluxes at high magnetic field where LLs from $v=0$ and from $v=4$ intersect, consistent with the Hofstadter butterfly spectra in the previous well-studied graphene/hBN superlattice heterostructure.[24–28] Additionally, quantum Hall ferromagnets (QHFMs), well developed here with LL filling factors $v_{\rm LL}=\pm 1$, $\pm 2$, $\pm 3$, 5, 6, 7 near charge neutral point (CNP), are formed by lifting the 8-fold spin, valley, and layer degeneracies. Close to the Hofstadter gaps at $\mathit{\Phi} =\mathit{\Phi} _{0}/q$, the system experiences a reverse Stoner transition due to superlattice-modulated bandwidth broadening, which is manifested by the suppression of QHFM. Instead of being quantized plateaus, the Hall resistances develop into a series of peaks and valleys with positive and negative values respectively, and the alternating Hall resistance sign change [see the red curve of Fig. 1(f)] is a direct evidence of recurring fractal Hofstadter Butterfly bands.
Fig. 2. Symmetry-broken Chern insulators in TBG with $\theta =1.25^{\circ}$. (a) Landau fan diagram of TBG acquired at temperature $T=20$ mK. (b) Schematic illustration of Wannier diagram. Grey and dark blue lines represent integer quantum Hall (IQH) insulators originated from CNP and gap edge of fulfilling, respectively. Light blue line shows symmetry-broken Chern insulators (SCIs) interrupted by fractal Hofstadter gaps at rational magnetic flux. Orange shades represent band gap at CNP and fulfillings. (c) Schematic of Chern number texture of flat band under $C_{2}$ symmetry or time-reversal symmetry breaking. The dashed lines indicate Fermi level location at $v=-1$ filling.
Next, we discuss the transport measurements at $\theta = \sim$$1.25^{\circ}, an angle closer to magic angle of \sim$$1.1^{\circ}$, and we observed a phases diagram beyond Hofstadter butterfly spectra. Similar to the device with $\theta = \sim$$1.38^{\circ}, correlated insulating states are absent at zero magnetic field, and LLs from v = 0 and from v = 4 at 20 mK are clearly shown in Fig. 2(a). The major difference lies in the emergence of longitudinal resistance R_{xx} minima, as depicted by the red dashed lines in Fig. 2(b). These states are characterized by a series of Chern number (or LL filling factor) of -1, -2, -3, which could be traced to integer moiré band fillings of v = -3, -2, -1 at B = 0 T, respectively. For simplicity, we name these states by (C, v), that is, (-1, -3), (-2, -2), and (-3, -1). These states are not originated from the fractal bands in Hofstadter Butterfly spectra. Take the state (-2, -2), for example, the onset magnetic field is very low, at B < 2 T (\sim0.05\mathit{\Phi}_{0}), which suggests that magnetic flux plays a trivial role [Fig. 2(a) and Fig. S6]. In addition, the state threads through several rational fluxes in a fan diagram, which is also in contrast to the expected restriction to one rational flux for fractal minibands. Similar threading behaviors are found in the states (-1, -3) and (-3, -1), which emerge at B=8 T and 12 T respectively. It is also noted that the state (-2, -2) is robust, and it shows a thermal activation gap comparable to LL (v_{\rm LL} =-2) originated from CNP (Fig. S4). The twist angle inhomogeneity here precludes well-developed Hall plateaus for these states (Fig. S5). We interpret these states as the emergence of topologically nontrivial moiré Chern insulators, which have also been observed in MA-TBG in recent related works.[18–21,29,30] As discussed in previous theoretical works, the low-energy flat bands in TBG can be interpreted as 8 valley-spin degenerate Chern bands with opposite Chern numbers +/-1, which carry opposite orbital magnetization and exhibit opposite orbital g factor.[22] As a result, once a (weak) perpendicular magnetic field is applied, time-reversal symmetry would be broken due to the orbital Zeeman effect, which splits the 8-fold degenerate flat bands into two sets of 4-fold degenerate Chern bands with C=+1 and C=-1 respectively, equivalent to zeroth pseudo LLs.[31] If an integer number (v) of the 4-fold degenerate Chern bands are filled (emptified) on the electron (hole) side, a gap would be opened up between the occupied and unoccupied bands due to exchange interactions, leading to an interaction-driven Chern insulator with Chern number |C|=4-|v|. Therefore, one would expect to see a sequence of Chern insulator states (-1, -3), (-2, -2), (-3, -1) as the filling factor v decreases from -1 to -3, which is clearly marked in Fig. 2(b). Such a picture is similar to the QHFM for cyclotron LLs in graphene.[32,33] In a finite magnetic field, Coulomb interactions tend to first break either one of the spin or valley symmetry forming twofold degenerate LLs, and eventually lift all fourfold degeneracy. Coming back to the flat bands of TBG, the above picture implies that a symmetry-broken Chern insulator (SCI) (-2, -2) should appear first with the increase of magnetic field, followed by the (-1, -3) or (-3, -1) SCI states. This argument is consistent with our observed hierarchy of SCIs, and explains why a tiny field (\sim 2 T) is required to establish the (-2, -2) state. It is worthwhile to note that although an approximate particle–hole symmetry is present in the continuum model of TBG, it is absent in a more realistic band structure of TBG including the effects of atomic corrugations,[34] in which the conduction flat band has a wider bandwidth than the valence band, which suppresses the emergence of the SCIs on the electron side. Our results demonstrate that the emergence of non-trivial band topology is not restricted to the magic angle regime. Aside from the R_{xy} data acquired at T=20 mK for \theta =1.38^{\circ}, we also measured R_{xx} responses for various \theta of 1.25^{\circ}, 1.38^{\circ}, and 1.43^{\circ} at an elevated temperature T=1.7 K, with B=0–9 T. For band filling -4 < v < 0, SCI (-2, -2) and (-1, -3) are observed in the \theta =1.25^{\circ} TBG device, but not the \theta =1.38^{\circ} and 1.43^{\circ} devices (see Fig. 1, Fig. 3, Fig. S7, and Fig. S8). Experimentally, it is surprising to find the SCI in the 1.25^{\circ} device for two reasons, one is the non-magic twist angle, and the other is a moderate twist angle inhomogeneity. Usually, the inhomogeneity would destroy the correlated behavior and impede the observation of SCI. However, seeing is believing, the observation of SCIs in turn suggests that the relatively large twist angle does not play a dominant role here, and an explanation of microscopic mechanism behind as well as what role played by strain is beyond our work. And instead, we tend to offer a tentative yet quantitative analysis of electron interaction and kinetic energy. We found that our observation of Chern bands at 1.25^{\circ} is very close to the crossover regime where Coulomb interaction and moiré band bandwidth are comparable (see Fig. S9). Since the moiré bandwidth W increases almost linearly with the twist angle (when \theta is larger than the MA), the disappearance of the SCI states at larger twist angles can be explained by the reduced interaction effects due to the enhanced kinetic energy. Note that the above analysis electron interaction is kind of rough, and it strongly depends on the dielectric constant and also on how the electron separation is treated. The critical twist angle below which such SCIs can emerge requires more detailed exploration in the future. As discussed above, the SCI states observed in the 1.25^{\circ} device is interpreted as an interaction-driven symmetry breaking state triggered by a tiny onset magnetic field. The gap in the SCI state is generated by electrons' Coulomb interactions, while the topological nature of the gap is triggered by a tiny B field through the orbital Zeeman coupling. Thus, the resulted symmetry breaking state exhibits an orbital ferromagnetic order with nonzero Chern numbers. At larger twist angles, the Coulomb interaction may not be strong enough to overcome the moiré kinetic energy, then the symmetry-breaking scenario sketched above cannot happen and the system would still preserve valley and spin symmetries (despite the non-interacting spin and orbital Zeeman splittings induced by the small B fields, which are negligible compared to the bandwidth without considering interaction effects). However, another possible scenario is that the absence of SCI states in 1.38^{\circ} and 1.43^{\circ} TBG is accounted by the intrinsic trivial topology of the low-energy bands at larger twist angles such that the system remains topologically trivial even if the interactions are strong enough to drive the system into a symmetry-breaking state at larger twist angles. Although the intervalley couplings at the single particle are still negligible (\sim0.05–0.1 meV) at 1.38^{\circ} and 1.43^{\circ}, it could be dramatically enhanced by electron–phonon coupling,[35] thus we cannot rule out this possibility. A method to distinguish these two scenarios is by further elevating magnetic fields to enhance exchange Coulomb energy which is proportional to \sqrt B or reducing device disorders, to fully break spin and valley isospin symmetry. If the nontrivial topology is preserved at larger twist angles, one would expect to see the emergence of SCIs with fully lifted valley-spin degeneracy at stronger magnetic fields. The specific mechanism for the interplay between electrons' interactions and the nontrivial band topology would be further elucidated by checking if the SCI states would emerge or not at high-enough magnetic fields. Lastly, we show the interplay and competition between these Chern bands, cyclotron LLs, and also Hofstadter minibands in magnetic fields. In fact, fractal Hofstadter Butterfly spectra are greatly suppressed with \theta \sim 1.25^{\circ} in Fig. 2(a). Usually a relatively smaller twist angle gives a longer moiré wavelength and a stronger moiré potential, which should help to develop the fractal bands. Instead, the color mapping in Fig. 2(a) yields a Landau fan diagram with almost no trace of fractal bands fanning out at a rational filling of quantum magnetic flux, which is beyond the framework of single-particle Hofstadter Butterfly picture as demonstrated in Figs. 1(d) and 1(e) with \theta=\sim$$1.38^{\circ}$. This argument is further supported in Figs. 3(d) and 3(e) at elevated temperature of $T = 1.7$ K. The fractal Hofstadter butterfly survives as the Brown–Zak oscillations[36] in Fig. 3(e) for $\theta \sim 1.38^{\circ}$, while they are absent in Fig. 3(d) for $\theta \sim 1.25^{\circ}$. The suppression of conventional fractal Hofstadter bands suggests a competition, which is sensitive to electron interactions between LL quantization effect and non-trivial topological effect at zero magnetic field from moiré Chern bands.
Fig. 3. Landau fan diagram in TBG with respect to twist angle $\theta$. (a), (b), (c) Density of states (DOS) versus band energy for varied twist angles $\theta =1.25^{\circ}$, 1.38$^{\circ}$, and 1.43$^{\circ}$, respectively, according to the Bistritzer–MacDonald model[1] calculations with a hopping energy of 110.7 meV. Orange shades represent single-particle band gap at $v=\pm 4$ fillings. (d), (e), (f) Landau fan diagrams for varied twist angles in one device. Dashed lines in panels (a) and (b) show $R_{xx}$ minima of SCI and Brown–Zak oscillations, respectively.
Here, we give a simple and self-consistent qualitative explanation in the framework of a competition between orbital Zeeman effect and cyclotron quantization. In one way, the orbital Zeeman effects from the nontrivial topology of the flat bands in TBG are characterized by effective orbital $g$-factor,[37–39] which may be enhanced by electron–electron interactions. In another way, the cyclotron quantization is proportional to Fermi velocity (or inversely proportional to effective mass). At a magic angle, the moiré band is ultra-flat with the Fermi velocity greatly suppressed, i.e., orbital effects dominate over cyclotron quantization, and thus moiré Chern bands are dominating in the Landau fan diagram with a greatly suppressed LL quantization and eventually suppressed fractal Hofstadter butterfly spectra. When the twist angle increases to a non-magic angle, the Fermi velocity tends to increase and thus leads to a situation where both orbital effects and cyclotron quantization matter. As a result, moiré Chern bands and cyclotron LLs coexist and show intermediate interactions, as shown by instances that the trajectory of SCI ($-2$, $-2$) beads between rational fluxes $p\mathit{\Phi}_{0}/10$ and is crossed by Landau fan emanating from CNP [Fig. 2(b)]. At some points, e.g., by increasing the twist angle or by tuning doping level, the orbital effects fail and cyclotron quantization will win, eventually moiré Chern bands give way to Hofstadter butterfly fractal gaps, as manifested by prevailed Brown–Zak oscillations in electron branch of 1.25$^{\circ}$ NMA-TBG and also both electron and hole branches of 1.38$^{\circ}$ NMA-TBG. For a 1.43$^{\circ}$ device, the disappearance of moiré Chern bands is strongly related to the reduced $U/W$ for an increased moiré bandwidth, and the disappearance of Brown–Zak oscillations is due to a smaller moiré period and weaker moiré potential, in line with fan diagram of a 1.8$^{\circ}$ TBG device in Ref. [13] where two graphene layers are weakly coupled. In summary, our results demonstrate non-trivial topology for low-energy flat bands in non-magic-angle TBG with broken $T$ symmetry. The stabilization of Chern insulators requires both a tiny magnetic field and strong electron interactions stemmed from flat bands, implying a related $T$ symmetry breaking mechanism. Our studies also point out a crucial role of electron interactions in shaping Landau level phases in TBG. These discoveries would help us unveil the mystery of electron correlation effects in band topology and even correlated insulators and superconductivities. Acknowledgements.—We appreciate the helpful discussion with Q. Wu and Y. Guan. The authors thank the finical supports from the National Key R&D program (Grant No. 2020YFA0309604), the National Natural Science Foundation of China (Grant Nos. 61888102, 11834017, and 12074413), the Strategic Priority Research Program of CAS (Grant Nos. XDB30000000 and XDB33000000), the Key-Area Research and Development Program of Guangdong Province (Grant No. 2020B0101340001), and the Research Program of Beijing Academy of Quantum Information Sciences (Grant No. Y18G11). J. L. acknowledges the start-up grant of ShanghaiTech University and the National Key R&D Program (Grant No. 2020YFA0309601). K.W. and T.T. acknowledge supports from the Elemental Strategy Initiative conducted by the MEXT, Japan (Grant No. JPMXP0112101001), JSPS KAKENHI (Grant No. JP20H00354), and the CREST (JPMJCR15F3), JST. Author Contributions.—C.S., W.Y., and G.Z. conceived the project. C.S. fabricated the devices and performed the transport measurements above 1.5 K. C.S., J.Y., F.Q., and L.L. performed the transport measurements in dilution refrigerator. L.L. provided continuum model calculations. K.W. and T.T. provided hexagonal boron nitride crystals. C.S., W.Y., and G.Z. analyzed the data. C.S., J.L., W.Y., and G.Z. wrote the paper. All authors discussed and commented on this work.
References Moire bands in twisted double-layer grapheneCorrelated insulator behaviour at half-filling in magic-angle graphene superlatticesUnconventional superconductivity in magic-angle graphene superlatticesTuning superconductivity in twisted bilayer grapheneSuperconductors, orbital magnets and correlated states in magic-angle bilayer grapheneElectronic correlations in twisted bilayer graphene near the magic angleMaximized electron interactions at the magic angle in twisted bilayer grapheneSpectroscopic signatures of many-body correlations in magic-angle twisted bilayer grapheneCharge order and broken rotational symmetry in magic-angle twisted bilayer grapheneNearly flat Chern bands in moiré superlatticesEmergent ferromagnetism near three-quarters filling in twisted bilayer grapheneIntrinsic quantized anomalous Hall effect in a moiré heterostructureSuperlattice-Induced Insulating States and Valley-Protected Orbits in Twisted Bilayer GrapheneTunable moiré bands and strong correlations in small-twist-angle bilayer grapheneNonvolatile switching of magnetic order by electric fields in an orbital Chern insulatorElectrically tunable correlated and topological states in twisted monolayer-bilayer grapheneTunable correlated Chern insulator and ferromagnetism in a moiré superlatticeChern Insulators and Topological Flat-bands in Magic-angle Twisted Bilayer GrapheneStrongly Correlated Chern Insulators in Magic-Angle Twisted Bilayer GrapheneHofstadter subband ferromagnetism and symmetry broken Chern insulators in twisted bilayer grapheneSymmetry broken Chern insulators and magic series of Rashba-like Landau level crossings in magic angle bilayer grapheneQuantum Valley Hall Effect, Orbital Magnetism, and Anomalous Hall Effect in Twisted Multilayer Graphene SystemsMoiré butterflies in twisted bilayer grapheneHofstadter’s butterfly and the fractal quantum Hall effect in moiré superlatticesCloning of Dirac fermions in graphene superlatticesMassive Dirac Fermions and Hofstadter Butterfly in a van der Waals HeterostructureHierarchy of Hofstadter states and replica quantum Hall ferromagnetism in graphene superlatticesHofstadter Butterfly and Many-Body Effects in Epitaxial Graphene SuperlatticeTracing out Correlated Chern Insulators in Magic Angle Twisted Bilayer GrapheneFlavour Hund's Coupling, Correlated Chern Gaps, and Diffusivity in Moiré Flat BandsPseudo Landau level representation of twisted bilayer graphene: Band topology and implications on the correlated insulating phaseQuantum Hall Ferromagnetism in GrapheneSpin and valley quantum Hall ferromagnetism in grapheneCrucial role of atomic corrugation on the flat bands and energy gaps of twisted bilayer graphene at the magic angle $\theta \sim 1.{08}^{\circ }$Valley Jahn-Teller Effect in Twisted Bilayer GrapheneHigh-temperature quantum oscillations caused by recurring Bloch states in graphene superlatticesTheory of correlated insulating behaviour and spin-triplet superconductivity in twisted double bilayer grapheneFirst principle calculation of the effective Zeeman's couplings in topological materialsChiral orbital current and anomalous magnetic moment in gapped graphene
 [1] Bistritzer R and MacDonald A H 2011 Proc. Natl. Acad. Sci. USA 108 12233 [2] Cao Y, Fatemi V, Demir A, Fang S, Tomarken S L, Luo J Y, Sanchez-Yamagishi J D, Watanabe K, Taniguchi T, Kaxiras E, Ashoori R C and Jarillo-Herrero P 2018 Nature 556 80 [3] Cao Y, Fatemi V, Fang S, Watanabe K, Taniguchi T, Kaxiras E and Jarillo-Herrero P 2018 Nature 556 43 [4] Yankowitz M, Chen S, Polshyn H, Zhang Y, Watanabe K, Taniguchi T, Graf D, Young A F and Dean C R 2019 Science 363 1059 [5] Lu X, Stepanov P, Yang W, Xie M, Aamir M A, Das I, Urgell C, Watanabe K, Taniguchi T, Zhang G, Bachtold A, MacDonald A H and Efetov D K 2019 Nature 574 653 [6] Choi Y, Kemmer J, Peng Y, Thomson A, Arora H, Polski R, Zhang Y, Ren H, Alicea J, Refael G, Oppen F V, Watanabe K, Taniguchi T and Nadj-Perge S 2019 Nat. Phys. 15 1174 [7] Kerelsky A, McGilly L J, Kennes D M, Xian L, Yankowitz M, Chen S, Watanabe K, Taniguchi T, Hone J, Dean C, Rubio A and Pasupathy A N 2019 Nature 572 95 [8] Xie Y, Lian B, Jäck B, Liu X, Chiu C L, Watanabe K, Taniguchi T, Bernevig B A and Yazdani A 2019 Nature 572 101 [9] Jiang Y, Lai X, Watanabe K, Taniguchi T, Haule K, Mao J and Andrei E Y 2019 Nature 573 91 [10] Zhang Y H, Mao D, Cao Y, Jarillo-Herrero P and Senthil T 2019 Phys. Rev. B 99 075127 [11] Sharpe A L, Fox E J, Barnard A W, Finney J, Watanabe K, Taniguchi T, Kastner M A and Goldhaber-Gordon D 2019 Science 365 605 [12] Serlin M, Tschirhart C L, Polshyn H, Zhang Y, Zhu J, Watanabe K, Taniguchi T, Balents L and Young A F 2020 Science 367 900 [13] Cao Y, Luo J Y, Fatemi V, Fang S, Sanchez-Yamagishi J D, Watanabe K, Taniguchi T, Kaxiras E and Jarillo-Herrero P 2016 Phys. Rev. Lett. 117 116804 [14] Kim K, DaSilvab A, Huangc S, Fallahazada B, Larentisa S, Taniguchid T, Watanabe K, LeRoyc B J, MacDonald A H and Tutuc E 2017 Proc. Natl. Acad. Sci. USA 114 3364 [15] Polshyn H, Zhu J, Kumar M A, Zhang Y, Yang F, Tschirhart C L, Serlin M, Watanabe K, Taniguchi T, MacDonald A H and Young A F 2020 arXiv:2004.11353v1 [16] Chen S, He M, Zhang Y H, Hsieh V, Fei Z, Watanabe K, Taniguchi T, Cobden D H, Xu X, Dean C R and Yankowitz M 2020 arXiv:2004.11340v1 [17] Chen G, Sharpe A L, Fox E J, Zhang Y H, Wang S, Jiang L, Lyu B, Li H, Watanabe K, Taniguchi T, Shi Z, Senthil T, Goldhaber-Gordon D, Zhang Y and Wang F 2020 Nature 579 56 [18] Wu S, Zhang Z, Watanabe K, Taniguchi T and Andrei E Y 2020 arXiv:2007.03735 [19] Nuckolls K P, Oh M, Wong D, Lian B, Watanabe K, Taniguchi T, Bernevig B A and Yazdani A 2020 arXiv:2007.03810 [20] Saito Y, Ge J, Rademaker L, Watanabe K, Taniguchi T, Abanin D A and Young A F 2020 arXiv:2007.06115 [21] Das I, Lu X, Herzog-Arbeitman J, Song Z D, Watanabe K, Taniguchi T, Bernevig B A and Efetov D K 2020 arXiv:2007.13390 [22] Liu J, Ma Z, Gao J and Dai X 2019 Phys. Rev. X 9 031021 [23] Bistritzer R and MacDonald A H 2011 Phys. Rev. B 84 035440 [24] Dean C R, Wang L, Maher P, Forsythe C, Ghahari F, Gao Y, Katoch J, Ishigami M, Moon P, Koshino M, Taniguchi T, Watanabe K, Shepard K L, Hone J and Kim P 2013 Nature 497 598 [25] Ponomarenko L A, Gorbachev R V, Yu G L, Elias D C, Jalil R, Patel A A, Mishchenko A, Mayorov A S, Woods C R, Wallbank J R, Mucha-Kruczynski M, Piot B A, Potemski M, Grigorieva I V, Novoselov K S, Guinea F, Fal'ko V I and Geim A K 2013 Nature 497 594 [26] Hunt B, Sanchez-Yamagishi J D, Young A F, Yankowitz M, LeRoy B J, Watanabe K, Taniguchi T, Moon P, Koshino M, Jarillo-Herrero P and Ashoori R C 2013 Science 340 1427 [27] Yu G L, Gorbachev R V, Tu J S, Kretinin A V, Cao Y, Jalil R, Withers F, Ponomarenko L A, Piot B A, Potemski M, Elias D C, Chen X, Watanabe K, Taniguchi T, Grigorieva I V, Novoselov K S, Fal'ko V I, Geim A K and Mishchenko A 2014 Nat. Phys. 10 525 [28] Yang W, Lu X, Chen G, Wu S, Xie G, Cheng M, Wang D, Yang R, Shi D, Watanabe K, Taniguchi T, Voisin C, Plaçais B, Zhang Y and Zhang G 2016 Nano Lett. 16 2387 [29] Choi Y, Kim H, Peng Y, Thomson A, Lewandowski C, Polski R, Zhang Y, Arora H S, Watanabe K, Taniguchi T, Alicea J and Nadj-Perge S 2020 arXiv:2008.11746 [30] Park J M, Cao Y, Watanabe K, Taniguchi T and Jarillo-Herrero P 2020 arXiv:2008.12296 [31] Liu J, Liu J and Dai X 2019 Phys. Rev. B 99 155415 [32] Nomura K and MacDonald A H 2006 Phys. Rev. Lett. 96 256602 [33] Young A F, Dean C R, Wang L, Ren H, Cadden-Zimansky P, Watanabe K, Taniguchi T, Hone J, Shepard K L and Kim P 2012 Nat. Phys. 8 550 [34] Lucignano P, Cataudella D A V, Ninno D and Cantele G 2019 Phys. Rev. B 99 195419 [35] Angeli M, Tosatti E and Fabrizio M 2019 Phys. Rev. X 9 041010 [36] Kumar R K, Chen X, Auton G H, Mishchenko A, Bandurin D A, Morozov S V, Cao Y, Khestanova E, Shalom M B, Kretinin A V, Novoselov K S, Eaves L, Grigorieva I V, Ponomarenko L A, Fal'ko V I and Geim A K 2017 Science 357 181 [37] Lee J Y, Khalaf E, Liu S, Liu X, Hao Z, Kim P and Vishwanath A 2019 Nat. Commun. 10 5333 [38] Song Z D, Sun S, Xu Y F, Nie S M, Weng H M, Fang Z and Dai X 2015 arXiv:1512.05084 [39] Koshino M 2011 Phys. Rev. B 84 125427