Abstract
Rhombohedral-stacked multilayer graphene hosts a pair of flat bands touching at zero energy, which should give rise to correlated electron phenomena that can be tuned further by an electric field. Moreover, when electron correlation breaks the isospin symmetry, the valley-dependent Berry phase at zero energy may give rise to topologically non-trivial states. Here we measure electron transport through hexagonal boron nitride-encapsulated pentalayer graphene down to 100 mK. We observed a correlated insulating state with resistance at the megaohm level or greater at charge density n = 0 and displacement field D = 0. Tight-binding calculations predict a metallic ground state under these conditions. By increasing D, we observed a Chern insulator state with C = −5 and two other states with C = −3 at a magnetic field of around 1 T. At high D and n, we observed isospin-polarized quarter- and half-metals. Hence, rhombohedral pentalayer graphene exhibits two different types of Fermi-surface instability, one driven by a pair of flat bands touching at zero energy, and one induced by the Stoner mechanism in a single flat band. Our results establish rhombohedral multilayer graphene as a suitable system for exploring intertwined electron correlation and topology phenomena in natural graphitic materials without the need for moiré superlattice engineering.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The data shown in the main figures are available from the Harvard Dataverse Repository at https://doi.org/10.7910/DVN/ISWXLA. The datasets generated during and/or analysed during this study are available from the corresponding author upon reasonable request.
References
Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).
Sharpe, A. L. et al. Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene. Science 365, 605–608 (2019).
Serlin, M. et al. Intrinsic quantized anomalous Hall effect in a moiré heterostructure. Science 367, 900–903 (2020).
Chen, G. et al. Evidence of a gate-tunable Mott insulator in a trilayer graphene moiré superlattice. Nat. Phys. 15, 237–241 (2019).
Chen, G. et al. Signatures of tunable superconductivity in a trilayer graphene moiré superlattice. Nature 572, 215–219 (2019).
Chen, G. et al. Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice. Nature 579, 56–61 (2020).
Regan, E. C. et al. Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. Nature 579, 359–363 (2020).
Tang, Y. et al. Simulation of Hubbard model physics in WSe2/WS2 moiré superlattices. Nature 579, 353–358 (2020).
Wang, L. et al. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861–866 (2020).
Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059–1064 (2019).
Cao, Y. et al. Tunable correlated states and spin-polarized phases in twisted bilayer–bilayer graphene. Nature 583, 215–220 (2020).
Shen, C. et al. Correlated states in twisted double bilayer graphene. Nat. Phys. 16, 520–525 (2020).
Liu, X. et al. Tunable spin-polarized correlated states in twisted double bilayer graphene. Nature 583, 221–225 (2020).
Polshyn, H. et al. Electrical switching of magnetic order in an orbital Chern insulator. Nature 588, 66–70 (2020).
Chen, G. et al. Tunable orbital ferromagnetism at noninteger filling of a moiré superlattice. Nano Lett. 22, 238–245 (2022).
Min, H. & MacDonald, A. H. Electronic structure of multilayer graphene. Prog. Theor. Phys. Suppl. 176, 227–252 (2008).
Zhang, F., Sahu, B., Min, H. & MacDonald, A. H. Band structure of ABC-stacked graphene trilayers. Phys. Rev. B 82, 035409 (2010).
Zhang, F., Jung, J., Fiete, G. A., Niu, Q. & MacDonald, A. H. Spontaneous quantum Hall states in chirally stacked few-layer graphene systems. Phys. Rev. Lett. 106, 156801 (2011).
Koshino, M. & McCann, E. Trigonal warping and Berry’s phase Nπ in ABC-stacked multilayer graphene. Phys. Rev. B 80, 165409 (2009).
Yang, N., Li, C., Tang, Y. & Yelgel, C. Electronic structure of ABC-stacked multilayer graphene and trigonal warping: a first principles calculation. J. Phys. Conf. Ser. 707, 012022 (2016).
Shi, Y. et al. Electronic phase separation in multilayer rhombohedral graphite. Nature 584, 210–214 (2020).
Zhou, H. et al. Half- and quarter-metals in rhombohedral trilayer graphene. Nature 598, 429–433 (2021).
Zhou, H., Xie, T., Taniguchi, T., Watanabe, K. & Young, A. F. Superconductivity in rhombohedral trilayer graphene. Nature 598, 434–438 (2021).
Pierucci, D. et al. Evidence for flat bands near the Fermi level in epitaxial rhombohedral multilayer graphene. ACS Nano 9, 5432–5439 (2015).
Kerelsky, A. et al. Moiréless correlations in ABCA graphene. Proc. Natl Acad. Sci. USA 118, e2017366118 (2021).
Weitz, R. T., Allen, M. T., Feldman, B. E., Martin, J. & Yacoby, A. Broken-symmetry states in doubly gated suspended bilayer graphene. Science 330, 812–816 (2010).
Freitag, F., Trbovic, J., Weiss, M. & Schönenberger, C. Spontaneously gapped ground state in suspended bilayer graphene. Phys. Rev. Lett. 108, 76602 (2012).
Bao, W. et al. Stacking-dependent band gap and quantum transport in trilayer graphene. Nat. Phys. 7, 948–952 (2011).
Velasco, J. et al. Transport spectroscopy of symmetry-broken insulating states in bilayer graphene. Nat. Nanotechnol. 7, 156–160 (2012).
Myhro, K. et al. Large tunable intrinsic gap in rhombohedral-stacked tetralayer graphene at half filling. 2D Mater. 5, 045013 (2018).
Pamuk, B., Baima, J., Mauri, F. & Calandra, M. Magnetic gap opening in rhombohedral-stacked multilayer graphene from first principles. Phys. Rev. B 95, 075422 (2017).
Jia, J., Gorbar, E. V. & Gusynin, V. P. Gap generation in ABC-stacked multilayer graphene: screening versus band flattening. Phys. Rev. B 88, 205428 (2013).
Zhou, H. et al. Isospin magnetism and spin-polarized superconductivity in Bernal bilayer graphene. Science 375, 774–778 (2022).
de la Barrera, S. C. et al. Cascade of isospin phase transitions in Bernal-stacked bilayer graphene at zero magnetic field. Nat. Phys. 18, 771–775 (2022).
Seiler, A. M. et al. Quantum cascade of correlated phases in trigonally warped bilayer graphene. Nature 608, 298–302 (2022).
Zhang, F., Min, H., Polini, M. & MacDonald, A. H. Spontaneous inversion symmetry breaking in graphene bilayers. Phys. Rev. B 81, 041402 (2010).
Jung, J., Zhang, F. & MacDonald, A. H. Lattice theory of pseudospin ferromagnetism in bilayer graphene: competing interaction-induced quantum Hall states. Phys. Rev. B 83, 115408 (2011).
Nandkishore, R. & Levitov, L. Quantum anomalous Hall state in bilayer graphene. Phys. Rev. B 82, 115124 (2010).
Vafek, O. & Yang, K. Many-body instability of Coulomb interacting bilayer graphene: renormalization group approach. Phys. Rev. B 81, 041401 (2010).
Lemonik, Y., Aleiner, I. & Fal’Ko, V. I. Competing nematic, antiferromagnetic, and spin-flux orders in the ground state of bilayer graphene. Phys. Rev. B 85, 245451 (2012).
Kharitonov, M. Antiferromagnetic state in bilayer graphene. Phys. Rev. B 86, 195435 (2012).
Xu, D. H. et al. Stacking order, interaction, and weak surface magnetism in layered graphene sheets. Phys. Rev. B 86, 201404 (2012).
Sun, K., Yao, H., Fradkin, E. & Kivelson, S. A. Topological insulators and nematic phases from spontaneous symmetry breaking in 2D Fermi systems with a quadratic band crossing. Phys. Rev. Lett. 103, 046811 (2009).
Streda, P. Theory of quantised Hall conductivity in two dimensions. J. Phys. C 15, L717 (1982).
Li, J., Tupikov, Y., Watanabe, K., Taniguchi, T. & Zhu, J. Effective Landau level diagram of bilayer graphene. Phys. Rev. Lett. 120, 47701 (2018).
McCann, E. & Fal’ko, V. I. Landau-level degeneracy and quantum Hall effect in a graphite bilayer. Phys. Rev. Lett. 96, 86805 (2006).
Slizovskiy, S., McCann, E., Koshino, M. & Fal’ko, V. I. Films of rhombohedral graphite as two-dimensional topological semimetals. Commun. Phys. 2, 164 (2019).
Kopnin, N. B., Ijäs, M., Harju, A. & Heikkilä, T. T. High-temperature surface superconductivity in rhombohedral graphite. Phys. Rev. B 87, 140503 (2013).
Ghazaryan, A., Holder, T., Berg, E. & Serbyn, M. Multilayer graphenes as a platform for interaction-driven physics and topological superconductivity. Phys. Rev. B 107, 104502 (2023).
Calvera, V., Kivelson, S. A. & Berg, E. Pseudo-spin order of Wigner crystals in multi-valley electron gases. Low Temp. Phys. 49, 679–700 (2023).
Ju, L. et al. Topological valley transport at bilayer graphene domain walls. Nature 520, 650–655 (2015).
Li, H. et al. Electrode-free anodic oxidation nanolithography of low-dimensional materials. Nano Lett. 18, 8011–8015 (2018).
Geisenhof, F. R. et al. Quantum anomalous Hall octet driven by orbital magnetism in bilayer graphene. Nature 598, 53–58 (2021).
Varlet, A. et al. Anomalous sequence of quantum Hall liquids revealing a tunable Lifshitz transition in bilayer graphene. Phys. Rev. Lett. 113, 116602 (2014).
Acknowledgements
We acknowledge helpful discussions with F. Zhang, T. Senthil, L. Levitov, L. Fu, Z. Dong and A. Patri. L.J. acknowledges support from a Sloan Fellowship. Work by Tonghang Han was supported by NSF grant number DMR- 2225925. The device fabrication of this work was supported by the STC Center for Integrated Quantum Materials, NSF grant number DMR-1231319. Device fabrication was carried out at the Harvard Center for Nanoscale Systems and MIT.Nano. Part of the device fabrication was supported by the USD(R&E) under contract no. FA8702-15-D-0001. K.W. and T.T. acknowledge support from the JSPS KAKENHI (grant numbers 20H00354, 21H05233 and 23H02052) and World Premier International Research Center Initiative (WPI), MEXT, Japan. H.P. acknowledges support by NSF grant number PHY-1506284 and AFOSR grant number FA9550-21-1-0216. A portion of this work was performed at the National High Magnetic Field Laboratory, which is supported by the National Science Foundation Cooperative Agreement No. DMR-2128556* and the State of Florida.
Author information
Authors and Affiliations
Contributions
L.J. supervised the project. Tonghang Han, Z.L., G.S., J.S. and J.W. performed the DC magneto-transport measurements. Tonghang Han and Tianyi Han fabricated the devices. K.W. and T.T. grew the hBN single crystals. H.P. contributed to the data analysis. All authors discussed the results and wrote the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Nanotechnology thanks Yuan Cao and the other, anonymous, reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Identification of rhombohedral pentalayer graphene.
a, AFM topography map of a pentalayer graphene sample on a SiO2/Si substrate. (The small region on the right corresponds to a graphene tetralayer.) b, Near-field infrared nanoscopy image of the same pentalayer graphene sample as in a, showing different contrast on the pentalayer region. The bright region corresponds to Bernal stacking and the darker region corresponds to rhombohedral stacking. c, Raman spectra taken at rhombohedral and Bernal stacking domains in the pentalayer graphene sample as in b. ω is the Raman shift.
Extended Data Fig. 2 Single particle band structure and density of states of rhombohedral multilayer graphene.
a–f, Tight-binding calculation of single-particle band structure and density of states (DOS) for rhombohedral multilayer graphene (layer number N = 2, 3, 4, 5, 7, 9). Due to the remote hopping terms, the band structure deviates from E ≈ kN at low energy. The rhombohedral pentalayer graphene has the flattest band among all layer numbers.
Extended Data Fig. 3 Correlated insulator at n = D = 0 in pentalayer rhombohedral graphene.
a, The calculated single-particle density of states (DOS) in Bernal bilayer, rhombohedral trilayer (time a factor of 15 and 3 for comparison), and pentalayer graphene at D = 0. The blue and orange shaded areas depict the DOS of the valence and conduction band. b–d, The n-D Rxx map of the bilayer, trilayer and pentalayer graphene measured at 2 K. Pentalayer graphene has a much larger DOS and band overlap at n = 0 compared to the bilayer and trilayer case and expects to be more conducting at D = n = 0. However, both bilayer and trilayer graphene are conducting at D = n = 0, while an insulating state appears at D = n = 0 in pentalayer graphene, indicating the non-single-particle origin of the insulating state. The off-diagonal line in the n-D map of pentalayer graphene is due to the big contact resistance when the bottom gate is zero.
Extended Data Fig. 4 Phase diagram for both electron and hole doping.
a, b, Color plots of four-probe resistance Rxx as a function of carrier density n and displacement field D for the hole doping side and electron doping side measured at B = 0 and a temperature of 100 mK. Colored dots label different phases including band insulator (BI), correlated insulator (CI), spin-polarized half metal (SPHM), isospin-polarized quarter metal (IPQM) and unpolarized metal (UP). c, d, Hall resistance Rxy and longitudinal resistance Rxx as a function of the out-of-plane magnetic field at the red dot in a. e, f, Rxy and Rxx at the yellow dot in a. The quarter metal c & d shows a clear anomalous Hall effect and magnetic hysteresis, indicating a net valley polarization. While the half metal e & f does not show anomalous Hall effect or magnetic hysteresis, indicating the absence of net valley polarization. Therefore, we conclude the half metal to be spin polarized but valley unpolarized.
Extended Data Fig. 5 Additional temperature dependence data of the correlated insulator.
a, Temperature dependence of the four-probe resistance Rxx measured at charge neutrality (n = 0). BI and CI stand for band insulator and correlated insulator. The correlated insulator develops below ~ 30 K. The two white arrows indicate the semimetal phase with an anomalous temperature dependence at the low-temperature region, discussed in Fig. 4c. b, Rxx versus D at even higher temperatures. As the temperature is increased, the resistive state at D = 0 disappears and evolves to a dip in Rxx. Circles trace the center of the bumps in Rxx.
Extended Data Fig. 6 The C = -3 state at D = 0.21 V/nm.
a, b, 2D color plot of Rxx and Rxy versus carrier density n and out-of-plane magnetic field B taken at D = 0.21 V/nm at a temperature of 100 mK. The dashed line indicates the C = -3 state. The C = -3 state is the only visible state on the hole doping side at low magnetic fields, in contrast to the electron side where all Landau levels appear at a similar magnetic field.
Extended Data Fig. 7 Tracing the C = -3 state towards zero magnetic field.
a, 2D color plot of dRxx/dB versus carrier density n and out-of-plane magnetic field B taken at D = 0.11 V/nm at a temperature of 100 mK. The dashed line indicates the C = -3 state. b, 2D color plot of dRxy/dn versus carrier density n and out-of-plane magnetic field B taken at D = 0.11 V/nm at a temperature of 100 mK. The dashed line indicates the C = -3 state. The C = -3 state is the only visible state on the hole doping side at low magnetic field and traces all the way to 0 T. On the electron side, a complete set of Landau levels is observed and disappears at around 1.5 T.
Extended Data Fig. 8 Chern insulators at the negative side of D.
a, 2D color plots of Rxx (upper panel) and Rxy (lower panel) at B = 1 T, revealing three Chern insulator states at the hole-doped side in the gap-closing range of D. The state with a Chern number C = -5 happens at D = -0.15 V/nm, while two states with C = -3 happen at D = -0.093 V/nm and D = -0.2 V/nm, as indicated by the dashed lines. b, c, d, 2D color plots of Rxx (upper panel) and Rxy (lower panel) versus n and B at D = -0.15 V/nm, -0.093 V/nm and -0.2 V/nm respectively. The dashed lines indicate the n-B relation of the Chern insulators as predicted by the Streda formula. These results are consistent with those of the positive D side.
Extended Data Fig. 9 Evolution of the band structure at charge-neutrality with D.
a, Rxx measured at n = 0 at 2 K. b, The band structure schematic of each isospin flavor at charge-neutrality under different D values. The colored dots correspond to the states in a, including LAF (layer antiferromagnet), SM (semimetal), QAH (quantum anomalous Hall) and LP (layer polarized state). The color of each band represents the valley Chern number (black: 5/2, red: -5/2), while the band with a grey color indicates the gap-closing case. c, The layer polarization of the LP, QAH and LAF state. The layer polarization of both conduction and valence band for each isospin flavor is shown, where green and orange color corresponds to the K’ and K valley. At D = 0, the system starts with the LAF state where charges are evenly distributed in the top and bottom layer in a spin-polarized manner (point A). As D increases, the gap of spin up (down) flavor expands (shrinks) due to its layer configuration (point B). It is important to note that the gap sizes of the two valleys with the same spin may differ. Consequently, one of the two gaps closes and reopens first (point C), leading to a situation where the gap of the K valley inverts while the gap of the K’ valley remains the same, resulting in the so-called QAH state (point D). Moreover, the layer polarization becomes partially polarized at this stage. If D continues to increase, the gap of the other valley eventually closes and reopens (point E), leading to the fully layer-polarized state (point F).
Extended Data Fig. 10 Temperature dependence of the SPHM, IPQM and UP states.
Temperature-dependent Rxx at D = 0.16 V/nm and n = 0 (red), D = 0 and n = -2.5*1012cm-2 (light grey, UP), D = 0.26 V/nm and n = -1.9*1012cm-2 (dark grey, SPHM), and D = 0.26 V/nm and n = -0.8*1012cm-2 (black, IPQM).
Supplementary information
Supplementary information
Supplementary Figs. 1–8 and Discussion (Sections I–VII).
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Han, T., Lu, Z., Scuri, G. et al. Correlated insulator and Chern insulators in pentalayer rhombohedral-stacked graphene. Nat. Nanotechnol. 19, 181–187 (2024). https://doi.org/10.1038/s41565-023-01520-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41565-023-01520-1
This article is cited by
-
Rhombohedral graphene goes correlated at four or five layers
Nature Nanotechnology (2024)
-
Layer-polarized ferromagnetism in rhombohedral multilayer graphene
Nature Communications (2024)
-
Spontaneous broken-symmetry insulator and metals in tetralayer rhombohedral graphene
Nature Nanotechnology (2024)
-
Fractional quantum anomalous Hall effect in multilayer graphene
Nature (2024)
-
Orbital multiferroicity in pentalayer rhombohedral graphene
Nature (2023)