Abstract
The constituent particles of matter can arrange themselves in various ways, giving rise to emergent phenomena that can be surprisingly rich and often cannot be understood by studying only the individual constituents. Discovering and understanding the emergence of such phenomena in quantum materialsâespecially those in which multiple degrees of freedom or energy scales are delicately balancedâis of fundamental interest to condensed-matter research1,2. Here we report on the surprising observation of emergent ferroelectricity in graphene-based moiré heterostructures. Ferroelectric materials show electrically switchable electric dipoles, which are usually formed by spatial separation between the average centres of positive and negative charge within the unit cell. On this basis, it is difficult to imagine grapheneâa material composed of only carbon atomsâexhibiting ferroelectricity3. However, in this work we realize switchable ferroelectricity in Bernal-stacked bilayer graphene sandwiched between two hexagonal boron nitride layers. By introducing a moiré superlattice potential (via aligning bilayer graphene with the top and/or bottom boron nitride crystals), we observe prominent and robust hysteretic behaviour of the graphene resistance with an externally applied out-of-plane displacement field. Our systematic transport measurements reveal a rich and striking response as a function of displacement field and electron filling, and beyond the framework of conventional ferroelectrics. We further directly probe the ferroelectric polarization through a non-local monolayer graphene sensor. Our results suggest an unconventional, odd-parity electronic ordering in the bilayer graphene/boron nitride moiré system. This emergent moiré ferroelectricity may enable ultrafast, programmable and atomically thin carbon-based memory devices.
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 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- 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 plots within this paper and other findings of this study are available from the corresponding authors upon reasonable request.
References
Keimer, B. & Moore, J. E. The physics of quantum materials. Nat. Phys. 13, 1045â1055 (2017).
Tokura, Y., Kawasaki, M. & Nagaosa, N. Emergent functions of quantum materials. Nat. Phys. 13, 1056â1068 (2017).
Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic properties of graphene. Rev. Mod. Phys. 81, 109â162 (2009).
Armitage, N. P., Mele, E. J. & Vishwanath, A. Weyl and Dirac semimetals in three-dimensional solids. Rev. Mod. Phys. 90, 015001 (2018).
Suárez Morell, E., Correa, J., Vargas, P., Pacheco, M. & Barticevic, Z. Flat bands in slightly twisted bilayer graphene: tight-binding calculations. Phys. Rev. B 82, 121407 (2010).
Li, G. et al. Observation of van Hove singularities in twisted graphene layers. Nat. Phys. 6, 109â113 (2010).
Bistritzer, R. & MacDonald, A. H. Moiré bands in twisted double-layer graphene. Proc. Natl Acad. Sci. USA 108, 12233â12237 (2011).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43â50 (2018).
Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80â84 (2018).
Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059â1064 (2019).
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. Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice. Nature 579, 56â61 (2020); correction 581, E3 (2020).
Burg, G. W. et al. Correlated insulating states in twisted double bilayer graphene. Phys. Rev. Lett. 123, 197702 (2019).
Liu, X. et al. Spin-polarized correlated insulator and superconductor in twisted double bilayer graphene. Nature 583, 221â225 (2020).
Cao, Y. et al. Tunable correlated states and spin-polarized phases in twisted bilayerâbilayer graphene. Nature 583, 215â220 (2020); correction 583, 215â220 (2020).
Shen, C. et al. Correlated states in twisted double bilayer graphene. Nat. Phys. 16, 520â525 (2020).
Wang, L. et al. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861â866 (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).
Nandkishore, R. & Levitov, L. Dynamical screening and excitonic instability in bilayer graphene. Phys. Rev. Lett. 104, 156803 (2010).
Fradkin, E., Kivelson, S. A., Lawler, M. J., Eisenstein, J. P. & Mackenzie, A. P. Nematic Fermi fluids in condensed matter physics. Annu. Rev. Condens. Matter Phys. 1, 153â178 (2010).
Fu, L. Parity-breaking phases of spinâorbit-coupled metals with gyrotropic, ferroelectric, and multipolar orders. Phys. Rev. Lett. 115, 026401 (2015).
Fernandes, R. M. & Venderbos, J. W. Nematicity with a twist: rotational symmetry breaking in a moiré superlattice. Sci. Adv. 6, eaba8834 (2020).
Kozii, V. & Fu, L. Odd-parity superconductivity in the vicinity of inversion symmetry breaking in spinâorbit-coupled systems. Phys. Rev. Lett. 115, 207002 (2015).
Mishra, A. & Lee, S. Topological multiferroic phases in the extended KaneâMeleâHubbard model in the Hofstadter regime. Phys. Rev. B 98, 235124 (2018).
Cao, Y. et al. Nematicity and competing orders in superconducting magic-angle graphene. Preprint at https://arxiv.org/abs/2004.04148 (2020).
Jiang, Y. et al. Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene. Nature 573, 91â95 (2019).
Choi, Y. et al. Imaging electronic correlations in twisted bilayer graphene near the magic angle. Nat. Phys. 15, 1174â1180 (2019); correction 15, 1205 (2019).
Zhang, Y.-H., Mao, D., Cao, Y., Jarillo-Herrero, P. & Senthil, T. Nearly flat Chern bands in moiré superlattices. Phys. Rev. B 99, 075127 (2019).
McCann, E. & Koshino, M. The electronic properties of bilayer graphene. Rep. Prog. Phys. 76, 056503 (2013).
Li, J., Martin, I., Büttiker, M. & Morpurgo, A. F. Topological origin of subgap conductance in insulating bilayer graphene. Nat. Phys. 7, 38â42 (2011).
Ju, L. et al. Topological valley transport at bilayer graphene domain walls. Nature 520, 650â655 (2015).
Sui, M. et al. Gate-tunable topological valley transport in bilayer graphene. Nat. Phys. 11, 1027â1031 (2015).
Shimazaki, Y. et al. Generation and detection of pure valley current by electrically induced Berry curvature in bilayer graphene. Nat. Phys. 11, 1032â1036 (2015).
Ju, L. et al. Tunable excitons in bilayer graphene. Science 358, 907â910 (2017).
Maher, P. et al. Evidence for a spin phase transition at charge neutrality in bilayer graphene. Nat. Phys. 9, 154â158 (2013).
Hunt, B. et al. Direct measurement of discrete valley and orbital quantum numbers in bilayer graphene. Nat. Commun. 8, 948 (2017).
Weitz, R. T., Allen, M., Feldman, B., Martin, J. & Yacoby, A. Broken-symmetry states in doubly gated suspended bilayer graphene. Science 330, 812â816 (2010).
Bao, W. et al. Evidence for a spontaneous gapped state in ultraclean bilayer graphene. Proc. Natl Acad. Sci. USA 109, 10802â10805 (2012).
Freitag, F., Trbovic, J., Weiss, M. & Schönenberger, C. Spontaneously gapped ground state in suspended bilayer graphene. Phys. Rev. Lett. 108, 076602 (2012).
Nam, Y., Ki, D.-K., Soler-Delgado, D. & Morpurgo, A. F. A family of finite-temperature electronic phase transitions in graphene multilayers. Science 362, 324â328 (2018).
Fei, Z. et al. Ferroelectric switching of a two-dimensional metal. Nature 560, 336â339 (2018).
Zhang, Y., Yuan, N. F. & Fu, L. Moiré quantum chemistry: charge transfer in transition metal dichalcogenide superlattices. Preprint at https://arxiv.org/abs/1910.14061 (2019).
Katayama, Y., Tsui, D., Manoharan, H., Parihar, S. & Shayegan, M. Charge transfer at double-layer to single-layer transition in double-quantum-well systems. Phys. Rev. B 52, 14817â14824 (1995).
Zhang, Y. et al. Direct observation of a widely tunable bandgap in bilayer graphene. Nature 459, 820â823 (2009).
Young, A. F. & Levitov, L. S. Capacitance of graphene bilayer as a probe of layer-specific properties. Phys. Rev. B 84, 085441 (2011).
Li, Y. et al. Probing symmetry properties of few-layer MoS2 and h-BN by optical second-harmonic generation. Nano Lett. 13, 3329â3333 (2013).
Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382â386 (2012).
Dean, C. R. et al. Hofstadterâs butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598â602 (2013).
Ponomarenko, L. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594â597 (2013).
Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427â1430 (2013).
Finney, N. R. et al. Tunable crystal symmetry in grapheneâboron nitride heterostructures with coexisting moiré superlattices. Nat. Nanotechnol. 14, 1029â1034 (2019).
Novoselov, K. S. et al. Unconventional quantum Hall effect and Berryâs phase of 2Ï in bilayer graphene. Nat. Phys. 2, 177â180 (2006).
Craciun, M. et al. Trilayer graphene is a semimetal with a gate-tunable band overlap. Nat. Nanotechnol. 4, 383â388 (2009).
Jhang, S. H. et al. Stacking-order dependent transport properties of trilayer graphene. Phys. Rev. B 84, 161408 (2011).
Wang, H., Wu, Y., Cong, C., Shang, J. & Yu, T. Hysteresis of electronic transport in graphene transistors. ACS Nano 4, 7221â7228 (2010).
McGilly, L. et al. Visualization of moiré superlattices. Nat. Nanotechnol. 15, 580â584 (2020).
Acknowledgements
We thank D. Bandurin, V. Fatemi, L. Levitov, Y. Lin, J. Mundy, R. Ramesh, J. Sanchez-Yamagishi, H. Shen, J. Song, S. Todadri, A. Vishwanath and N. Yuan for discussions; and T. Dinh for initial efforts on this project. Work in the P.J.-H. group was supported by the US DOE, BES Office, Division of Materials Sciences and Engineering under award DE-SC0001819 (device fabrication and transport measurements), the Center for the Advancement of Topological Semimetals, an Energy Frontier Research Center funded by the US Department of Energy Office of Science, through the Ames Laboratory under contract DE-AC02-07CH11358 (data analysis), and the Gordon and Betty Moore Foundationâs EPiQS Initiative through grant GBMF9643 to P.J.-H. The development of new nanofabrication and characterization techniques enabling this work has been supported by the US DOE Office of Science, BES, under award DE-SC0019300. Partial support for measurement and characterization training was through AFOSR grant FA9550-16-1-0382. This work made use of the Materials Research Science and Engineering Center Shared Experimental Facilities supported by the National Science Foundation (NSF) (grant number DMR-0819762). N.G. and S.-Y.X. acknowledge support from DOE, BES DMSE (data taking and analysis), and National Science Foundation under grant number DMR-1809815 (manuscript writing). K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan, grant number JPMXP0112101001, JSPS KAKENHI grant numbers JP20H00354 and the CREST(JPMJCR15F3), JST. R.A. (capacitance measurements), Z.B., Y.Z. and L.F. (theory) acknowledge support from NSF Science and Technology Center for Integrated Quantum Materials grant DMR-1231319. M.-H.L. was supported by Taiwan Ministry of Science and Technology (MOST) under grant numbers 109-2112-M-006-020-MY3 and 108-2638-M-006-002-MY2. N.M. and J.K. acknowledge the support by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES) under award DE-SC0020042.
Author information
Authors and Affiliations
Contributions
Q.M. and S.-Y.X. conceived the idea and experiment. Z.Z. fabricated devices, performed transport measurements and analysed data under the supervision of Q.M. and P.J.-H. S.d.l.B. performed capacitance measurements with the help of Z.Z. under the supervision of R.A. Z.B., Y.Z. and N.K. performed theoretical modelling and bandstructure calculations under the supervision of S.-Y.X., Q.M., N.G. and L.F. M.-H.L. performed the simulation of the transport behaviours for the hysteretic devices. N.M. performed second-harmonic generation measurements of BN flakes under the supervision of J.K. and W.A.T. K.W. and T.T. grew the bulk BN single crystals. All authors discussed the results and wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature thanks Kayoung Lee and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Resistance data summary for devices M1, N0, T1, H1 and H4.
Device schematic (a, d, g, j, n), dual-gate resistance map (b, e, h, k, o), nâD map (forward and backward) (l, p), and resistance as a function of externally applied displacement field at zero doping (c, f, i, m, q) for representative devices M1 (aâc), N0 (dâf), T1 (gâi), H1 (jâm) and H4 (nâq). The line traces along the electric-field direction are marked by black dashed lines. The superlattice resistance peaks are marked by cyan dashed lines. Note that the horizontal resistance line in h stems from a region of the sample that is only controlled by the top gate. Note that line traces in q are taken at next = 0 from the nâD map (p), whereas line traces in Fig. 1f in the main text are taken at VTG = 0 from the dual-gate map (Supplementary Fig. 25e, f), hence the difference in resistance magnitude.
Extended Data Fig. 2 Shifting of the hysteretic behaviour in device H4.
aâc, Forward (a) and backward (b) scans of the four-probe longitudinal resistance as a function of VTG and VBG and their difference (c). dâf, gâi, Same measurements as in aâc, except that VBG is swept between 0 V and 50 V (dâf) and between â40 V and 0 V (gâi). We present a phenomenological model to simulate the resistance maps in Supplementary Section V.3.2.
Extended Data Fig. 3 Hysteretic transport behaviour for device H2.
aâd, The four-probe resistance as a function of VBG (range, â10 V to 10 V) and VTG with ranges from 0 V to 5 V (a),from 5 V to â5 V (b), from â5 V to 10 V (c), and from 10 V to â5 V (d). The scan sequences are specified in the insets.
Extended Data Fig. 4 Shifting of the hysteretic behaviour in device H2.
aâc, Four-probe longitudinal resistance as a function of VTG and VBG for the forward (a) and backward (b) scans and their difference (c). dâi, Same measurements as in aâc, except that VTG is measured within 0Â V to 5 V (dâf) and within â5 V to 0Â V (gâi).
Extended Data Fig. 5 Hysteretic signature in Hall measurements for device H2.
a, b, The resistance measured while sweeping the externally applied displacement field Dext in the forward (a) and backward (b) direction at each fixed carrier density next. The carrier density scan direction is from the negative to positive values. câe, Carrier density extracted from Hall measurements along the lines L1 (c), L2 (d) and L3 (e) denoted in a. Red and blue curves were taken during the forward and backward scan of Dext, respectively.
Extended Data Fig. 6 Probing the out-of-plane electrical polarization using the top MLG sensor in device H4.
a, Experimental configuration for measurements of the conductance of the BLG. b, Measured conductance of the BLG as a function of VBG (top gate is grounded). The red and blue curves correspond to the forward and backward VBG scans, respectively. The vertical dashed lines denote the VBG values that correspond to the charge-neutrality point of the BLG for forward and backward scans. c, Experimental configuration for the measurements of the conductance of the tp-MLG with gate voltage VBG applied to bottom metal gate and the BLG grounded. d, Measured conductance of the tp-MLG as a function of VBG with the experimental configuration in c. The gate voltages are the same as in b (Methods). The red and blue dots denote the conductance of the tp-MLG when the BLG is charge neutral. e, Experimental configuration for the measurements of the conductance of the tp-MLG with gate voltage VBL applied to the BLG. f, Measurements of the conductance of the tp-MLG as a function of VBL with the experimental configuration in e. The conductance of the monolayer at the red and blue dots in d can be inversely mapped to two different VBL values, which corresponds to the difference of electrostatic potentials on the top surface of BLG induced by ferroelectric switching (Methods). The in-plane bias voltage VSD was kept below 1Â mV for all the measurements.
Extended Data Fig. 7 Hysteresis signature in the bulk electronic compressibility of device H2.
a, b, Bottom capacitance Cb between the bottom gate and BLG as a function of the externally applied field, Dext as the fast-scan axis, and gate-defined carrier density, next as the slow-scan axis. The white arrows indicate the sweep direction of Dext in each panel. Deviations of the capacitance from the geometric value reflect modulations in the electronic compressibility, ân/âμ, from the total area of BLG overlapping the bottom gate. Data were collected by sweeping the displacement field at each fixed carrier density, as in Fig. 2h, i. Dark features indicate regions of incompressibility resulting from the opening of a gap in the BLG. The gapless point, a compressible state with high Cb, is achieved at a finite Dext that depends on the sweep direction. c, Forward and backward traces from a and b at a fixed next. d, Resistance traces at the same density showing resistance peaks corresponding to the incompressible features in c. e, Circuit schematic of the bottom gate capacitance measurement, including a two-stage cryogenic amplifier (enclosed in dashed box). Capacitance is measured by applying a small a.c. excitation voltage to the bottom gate, δVBG, while also applying a nearly 180° out-of-phase signal, δVref, to a reference capacitor, Cref to null the voltage at the bridge balance point, (B). Deviations in the balanced signal caused by variations in compressibility are amplified by two high electron-mobility transistors and measured at the drain of the second stage, δVout. Carrier density next and external field Dext are controlled by top- and bottom-gate d.c. voltages VTG and VBG, in the same way as in the transport measurements.
Extended Data Fig. 8 Independence of the hysteretic behaviour on the sweep rate.
aâc, Forward (red) and backward (blue) sweeps of the bottom-gate capacitance, Cb, from device H2 at fixed carrier density next with sweep rates of 2.2 mV nmâ1 sâ1 (a), 4.5 mV nmâ1 sâ1 (b), and 9.8 mV nmâ1 sâ1 (c). Sweep rates shown in each panel denote the rate at which the externally applied displacement field Dext/ε0 was ramped in the BN dielectric layers. No noticeable variation was observed in the capacitance features for the large range of sweep rates.
Supplementary information
Supplementary Information
This file contains Supplementary Notes and Discussions and includes Supplementary Figures 1â33.
Rights and permissions
About this article
Cite this article
Zheng, Z., Ma, Q., Bi, Z. et al. Unconventional ferroelectricity in moiré heterostructures. Nature 588, 71â76 (2020). https://doi.org/10.1038/s41586-020-2970-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-020-2970-9
This article is cited by
-
Engineering band structures of two-dimensional materials with remote moiré ferroelectricity
Nature Communications (2024)
-
Robust flat bands in twisted trilayer graphene moiré quasicrystals
Nature Communications (2024)
-
Atomic-level polarization reversal in sliding ferroelectric semiconductors
Nature Communications (2024)
-
Multiresistance states in ferro- and antiferroelectric trilayer boron nitride
Nature Communications (2024)
-
A 2D ferroelectric vortex pattern in twisted BaTiO3 freestanding layers
Nature (2024)