Abstract
The discovery of correlated states and superconductivity in magic-angle twisted bilayer graphene (MATBG) established a new platform to explore interaction-driven and topological phenomena. However, despite multitudes of correlated phases observed in moiré systems, robust superconductivity appears the least common, found only in MATBG and more recently in magic-angle twisted trilayer graphene. Here we report the experimental realization of superconducting magic-angle twisted four-layer and five-layer graphene, hence establishing alternating twist magic-angle multilayer graphene as a robust family of moiré superconductors. This finding suggests that the flat bands shared by the members play a central role in the superconductivity. Our measurements in parallel magnetic fields, in particular the investigation of Pauli limit violation and spontaneous rotational symmetry breaking, reveal a clear distinction between the N = 2 and N > 2-layer structures, consistent with the difference between their orbital responses to magnetic fields. Our results expand the emergent family of moiré superconductors, providing new insight with potential implications for design of new superconducting materials platforms.
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 that support the current study are available from Harvard Dataverse50.
Code availability
The codes that support the current study are available from Harvard Dataverse50.
References
Andrei, E. Y. et al. The marvels of moiré materials. Nat. Rev. Mater. 6, 201–206 (2021).
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).
Chen, G. et al. Evidence of a gate-tunable Mott insulator in a trilayer graphene moiré superlattice. Nat. Phys. 15, 237 (2019).
Burg, G. W. et al. Correlated insulating states in twisted double bilayer graphene. Phys. Rev. Lett. 123, 197702 (2019).
Shen, C. et al. Correlated states in twisted double bilayer graphene. Nat. Phys. 16, 520–525 (2020).
Cao, Y. et al. Tunable correlated states and spin-polarized phases in twisted bilayer-bilayer graphene. Nature 583, 215–220 (2020).
Liu, X. et al. Tunable spin-polarized correlated states in twisted double bilayer graphene. Nature 583, 221–225 (2020).
He, M. et al. Symmetry breaking in twisted double bilayer graphene. Nat. Phys. 17, 26–30 (2020).
Polshyn, H. et al. Electrical switching of magnetic order in an orbital Chern insulator. Nature 588, 66–70 (2020).
Xu, S. et al. Tunable van Hove singularities and correlated states in twisted monolayer-bilayer graphene. Nat. Phys. 17, 619–626 (2021).
Chen, S. et al. Electrically tunable correlated and topological states in twisted monolayer–bilayer graphene. Nat. Phys. 17, 374–380 (2021).
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).
Xu, Y. et al. Correlated insulating states at fractional fillings of moiré superlattices. Nature 587, 214–218 (2020).
Jin, C. et al. Stripe phases in WSe2/WS2 moiré superlattices. Nat. Mater. 20, 940–944 (2021).
Li, T. et al. Continuous Mott transition in semiconductor moiré superlattices. Nature 597, 350–354 (2021).
Li, T. et al. Quantum anomalous Hall effect from intertwined moiré bands. Nature 600, 641–646 (2021).
Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059–1064 (2019).
Lu, X. et al. Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene. Nature 574, 653–657 (2019).
Chen, G. et al. Signatures of tunable superconductivity in a trilayer graphene moiré superlattice. Nature 572, 215–219 (2019).
Zhang, X. et al. Correlated insulating states and transport signature of superconductivity in twisted trilayer graphene superlattices. Phys. Rev. Lett. 127, 166802 (2021).
Park, J. M., Cao, Y., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Tunable strongly coupled superconductivity in magic-angle twisted trilayer graphene. Nature 590, 249–255 (2021).
Hao, Z. et al. Electric field-tunable superconductivity in alternating-twist magic-angle trilayer graphene. Science 371, 1133–1138 (2021).
Cao, Y., Park, J. M., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Pauli-limit violation and re-entrant superconductivity in moiré graphene. Nature 595, 526–531 (2021).
Suárez Morell, E., Correa, J. D., Vargas, P., Pacheco, M. & Barticevic, Z. Flat bands in slightly twisted bilayer graphene: tight-binding calculations. Phys. Rev. B 82, 121407 (2010).
Bistritzer, R. & MacDonald, A. H. Moiré bands in twisted double-layer graphene. Proc. Natl Acad. Sci. USA 108, 12233–12237 (2011).
Lopes dos Santos, J. M. B., Peres, N. M. R. & Castro Neto, A. H. Continuum model of the twisted graphene bilayer. Phys. Rev. B 86, 155449 (2012).
Xie, Y. et al. Spectroscopic signatures of many-body correlations in magic-angle twisted bilayer graphene. Nature 572, 101–105 (2019).
Wong, D. et al. Cascade of electronic transitions in magic-angle twisted bilayer graphene. Nature 582, 198–202 (2020).
Zondiner, U. et al. Cascade of phase transitions and Dirac revivals in magic-angle graphene. Nature 582, 203–208 (2020).
Park, J. M., Cao, Y., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Flavour Hund’s coupling, Chern gaps and charge diffusivity in moiré graphene. Nature 592, 43–48 (2021).
Oh, M. et al. Evidence for unconventional superconductivity in twisted bilayer graphene. Nature 600, 240–245 (2021).
Pierce, A. T. et al. Unconventional sequence of correlated Chern insulators in magic-angle twisted bilayer graphene. Nat. Phys. 17, 1210–1215 (2021).
Xie, Y. et al. Fractional Chern insulators in magic-angle twisted bilayer graphene. Nature 600, 439–443 (2021).
Khalaf, E., Kruchkov, A. J., Tarnopolsky, G. & Vishwanath, A. Magic angle hierarchy in twisted graphene multilayers. Phys. Rev. B 100, 085109 (2019).
Nam, N. N. T. & Koshino, M. Lattice relaxation and energy band modulation in twisted bilayer graphene. Phys. Rev. B 96, 075311 (2017).
Turkel, S. et al. Orderly disorder in magic-angle twisted trilayer graphene. Science 376, 193–199 (2022).
Rodan-Legrain, D. et al. Highly tunable junctions and non-local Josephson effect in magic-angle graphene tunnelling devices. Nat. Nanotechnol. 116, 769–775 (2021).
de Vries, F. K. et al. Gate-defined Josephson junctions in magic-angle twisted bilayer graphene. Nat. Nanotechnol. 16, 760–763 (2021).
Cao, Y. et al. Strange metal in magic-angle graphene with near Planckian dissipation. Phys. Rev. Lett. 124, 076801 (2020).
Polshyn, H. et al. Large linear-in-temperature resistivity in twisted bilayer graphene. Nat. Phys. 15, 1011–1016 (2019).
Dean, C. R. et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 5, 722–726 (2010).
Qin, W. & MacDonald, A. H. In-plane critical magnetic fields in magic-angle twisted trilayer graphene. Phys. Rev. Lett.127, 097001 (2021).
Cao, Y. et al. Nematicity and competing orders in superconducting magic-angle graphene. Science 372, 264–271 (2021).
Ledwith, P. J. et al. TB or not TB? Contrasting properties of twisted bilayer graphene and the alternating twist η-layer structures (η=3, 4, 5, …). Preprint at https://arxiv.org/abs/2111.11060 (2021).
Lake, E. & Senthil, T. Reentrant superconductivity through a quantum Lifshitz transition in twisted trilayer graphene. Phys. Rev. B 104, 174505 (2021).
Moshe, M., Kogan, V. G. & Mints, R. G. Edge-type Josephson junctions in narrow thin-film strips. Phys. Rev. B 78, 020510 (2008).
Park, J. et al. Data for robust superconductivity in magic-angle multilayer graphene family. Harvard Dataverse https://doi.org/10.7910/DVN/XZV9ON (2022).
Acknowledgements
We thank A. Vishwanath, E. Khalaf and P. Ledwith for fruitful discussions. This work was primarily supported by the US Department of Energy (DOE), Office of Basic Energy Sciences (BES), Division of Materials Sciences and Engineering under award DE-SC0001819 (J.M.P. and S.S.). Help with transport measurements and data analysis were supported by the National Science Foundation (DMR-1809802), and the STC Center for Integrated Quantum Materials (NSF grant no. DMR-1231319) (Y.C.). Help with device fabrication was supported by the Air Force Office of Scientific Research (AFOSR) 2DMAGIC MURI FA9550-19-1-0390 (L.X.). P.J.-H. acknowledges support from the Gordon and Betty Moore Foundation’s EPiQS Initiative through grant no. GBMF9463. P.J.-H. acknowledges partial support by the Fundación Ramon Areces and the CIFAR Quantum Materials programme. The development of new nanofabrication and characterization techniques enabling this work was supported by the US DOE Office of Science, BES, under award DE-SC0019300. K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by MEXT, Japan (grant no. JPMXP0112101001), JSPS KAKENHI (grant no. JP20H00354), CREST(grant no. JPMJCR15F3) and JST. This work made use of the Materials Research Science and Engineering Center Shared Experimental Facilities supported by the National Science Foundation (grant no. DMR-0819762) and of Harvard’s Center for Nanoscale Systems, supported by the NSF (grant no. ECS-0335765).
Author information
Authors and Affiliations
Contributions
J.M.P. and Y.C. fabricated the samples and performed transport measurements, with help from L.-Q.X. and S.S. K.W. and T.T. provided hBN samples. J.M.P. and Y.C. performed numerical simulations. J.M.P., Y.C. and P.J-H. performed data analysis, discussed the results and wrote the manuscript with input from all co-authors.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Materials thanks Chun Ning Lau, Aaron Sharpe and the other, anonymous, reviewer(s) for their contribution to the peer review.
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 Calculated band structure of MATNG in the presence of electric displacement field \(\frac{D}{{\varepsilon }_{0}}=0.5\) V/nm.
The calculation uses the screened electrostatic potential calculated as illustrated in Extended Data Fig. 6c. For all N>2 MATNG structures, the electric displacement field hybridizes the flat bands with other dispersive bands.
Extended Data Fig. 2 ν-D phase diagrams of MAT4G and MAT5G devices.
Devices 4A-C shown in (a-c) are MAT4G devices, while Devices 5A-E shown in (d-i) are MAT5G devices. Device 5B-1 and 5B-2 are two different contacts from the same device that shares the same top gate and back gate. All measurements were performed at T ~ 200 mK. In all MAT5G devices, we find the superconductivity region to be either approaching or exceeding ν = 4, which is indicated by blue dashed lines in (d-i).
Extended Data Fig. 3 Landau fan diagrams of MAT4G and MAT5G devices.
All measurements were performed at T ~ 200 mK.
Extended Data Fig. 4 Superconducting properties in other MAT4G and MAT5G devices.
(a-c) Temperature-dependent Vxx − I curves in device 4A, 4C and 5D. (d-f) ν − T phase diagrams in device 4A, 4C and 5D. (g-i) TBKT and Ginzburg-Landau coherent length ξGL in device 4C, 5D and 5E. The error bars in (g-i) correspond to ξGL extracted using p=30% (point), 35% (lower bound), and 25% (upper bound).
Extended Data Fig. 5 ν-D map of BKT transition temperature in device 4B, 5A, 5B-1, and 5C.
(a-d) TBKT as a function of ν and D in device 4B (a), 5A (b), 5B-1 (c), and 5C (d).
Extended Data Fig. 6 Modeling layer-inhomogeneous screening in MATNG.
(a) Schematic of N graphene layers being gated by a top gate (TG) and bottom gate (BG). The electric displacement field between the graphene layers is reduced compared to the field outside the stack due to screening. (b) By assuming a finite density of states \({{{\mathcal{D}}}}\) on each graphene layer, we model the N-layer system as a capacitor network, where \({{{{\rm{c}}}}}_{{{{\rm{q}}}}}={{{{\rm{e}}}}}^{2}{{{\mathcal{D}}}}\) is the quantum capacitance and cg = ε0/d is the geometric capacitance. (c) Calculated electrostatic potential vi on each layer for MATBG, MATTG, MAT4G, MAT5G and MAT6G (N = 2, …, 6), assuming bandwidth W=20 meV, and external displacement field \(\frac{{{{\rm{D}}}}}{{\varepsilon }_{0}}=0.5V/nm\). The twist angle is the same as in Fig. 4, except for 6L which is 1.99°. Here the midplane (the plane of the \(\frac{{{{\rm{N}}}}+1}{2}\)-th layer when N is odd and the midplane between \(\frac{{{{\rm{N}}}}}{2}\)-th layer and (\(\frac{{{{\rm{N}}}}}{2}+1\))-th layer when N is even) is set to be zero both in layer position and in electrostatic potential. For comparison, the electrostatic potential without screening eDdx/ε0 is shown as the dashed line, where \(-\frac{({{{\rm{N}}}}-1)}{2} < {{{\rm{x}}}} < \frac{({{{\rm{N}}}}-1)}{2}\) is the layer position (horizontal axis) and d=0.34 nm is the interlayer distance).
Extended Data Fig. 7 Typical optical microscope images for MAT4G and MAT5G devices.
The illustration shows the split-gate geometry for the Josephson junction. Data from all contacts including the ones not shown in the figures exhibit robust superconductivity.
Extended Data Fig. 8 In-plane magnetic field response in other MAT4G and MAT5G devices.
(a-d) Absence of rotational symmetry breaking in other MAT4G and MAT5G devices. The filling factors for (a-d) are ν = − 3.41, − 1.95, − 3.50, − 2.33, respectively, and the electric displacement fields are \(\frac{{{{\rm{D}}}}}{{\varepsilon }_{0}}=0.29V/nm,0.29V/nm,0V/nm,0V/nm\), respectively. (d-e) Pauli limit violation in Device 4C measured at ν = 2.29 and \(\frac{{{{\rm{D}}}}}{{\varepsilon }_{0}}=0V/nm\) (d) and Device 5E measured at ν = 2.59 and \(\frac{{{{\rm{D}}}}}{{\varepsilon }_{0}}=0V/nm\) (e).
Extended Data Fig. 9 D-independent Pauli limit violation in MAT4G.
(a-b) shows the D-B∥and D-T map of Rxx in device 4B. We find the critical in-plane field and critical temperature follow similar trend with D, indicating that the PVR, which is proportional to the ratio between them, is largely independent of D. (c-d) Same trend is observed in device 4C.
Rights and permissions
About this article
Cite this article
Park, J.M., Cao, Y., Xia, LQ. et al. Robust superconductivity in magic-angle multilayer graphene family. Nat. Mater. 21, 877–883 (2022). https://doi.org/10.1038/s41563-022-01287-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41563-022-01287-1
This article is cited by
-
Engineering interlayer hybridization in van der Waals bilayers
Nature Reviews Materials (2024)
-
Electric field tunable bandgap in twisted double trilayer graphene
npj 2D Materials and Applications (2024)
-
Twisted-layer boron nitride ceramic with high deformability and strength
Nature (2024)
-
Interlayer exciton dynamics of transition metal dichalcogenide heterostructures under electric fields
Nano Research (2024)
-
Mixed-dimensional moiré systems of twisted graphitic thin films
Nature (2023)