Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Unconventional ferroelectricity in moiré heterostructures

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

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Quadratic, chiral fermions on a moiré potential in a BLG/BN moiré superlattice.
Fig. 2: Hysteretic transport behaviour for device H4.
Fig. 3: Measuring the out-of-plane electric polarization and a possible microscopic picture based on interlayer charge transfer.
Fig. 4: Robustness of the ferroelectric switching in the graphene/BN moiré system.

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

  1. Keimer, B. & Moore, J. E. The physics of quantum materials. Nat. Phys. 13, 1045–1055 (2017).

    CAS  Google Scholar 

  2. Tokura, Y., Kawasaki, M. & Nagaosa, N. Emergent functions of quantum materials. Nat. Phys. 13, 1056–1068 (2017).

    CAS  Google Scholar 

  3. 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).

    ADS  CAS  Google Scholar 

  4. Armitage, N. P., Mele, E. J. & Vishwanath, A. Weyl and Dirac semimetals in three-dimensional solids. Rev. Mod. Phys. 90, 015001 (2018).

    ADS  MathSciNet  CAS  Google Scholar 

  5. 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).

    ADS  Google Scholar 

  6. Li, G. et al. Observation of van Hove singularities in twisted graphene layers. Nat. Phys. 6, 109–113 (2010).

    Google Scholar 

  7. Bistritzer, R. & MacDonald, A. H. Moiré bands in twisted double-layer graphene. Proc. Natl Acad. Sci. USA 108, 12233–12237 (2011).

    ADS  CAS  PubMed  Google Scholar 

  8. Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).

    ADS  CAS  PubMed  Google Scholar 

  9. Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018).

    ADS  CAS  PubMed  Google Scholar 

  10. Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059–1064 (2019).

    ADS  CAS  PubMed  Google Scholar 

  11. Sharpe, A. L. et al. Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene. Science 365, 605–608 (2019).

    ADS  CAS  PubMed  Google Scholar 

  12. Serlin, M. et al. Intrinsic quantized anomalous Hall effect in a moiré heterostructure. Science 367, 900–903 (2020).

    ADS  CAS  PubMed  Google Scholar 

  13. Chen, G. et al. Evidence of a gate-tunable Mott insulator in a trilayer graphene moiré superlattice. Nat. Phys. 15, 237–241 (2019).

    CAS  Google Scholar 

  14. Chen, G. et al. Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice. Nature 579, 56–61 (2020); correction 581, E3 (2020).

    ADS  CAS  PubMed  Google Scholar 

  15. Burg, G. W. et al. Correlated insulating states in twisted double bilayer graphene. Phys. Rev. Lett. 123, 197702 (2019).

    ADS  CAS  PubMed  Google Scholar 

  16. Liu, X. et al. Spin-polarized correlated insulator and superconductor in twisted double bilayer graphene. Nature 583, 221–225 (2020).

    ADS  CAS  PubMed  Google Scholar 

  17. 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).

    ADS  CAS  PubMed  Google Scholar 

  18. Shen, C. et al. Correlated states in twisted double bilayer graphene. Nat. Phys. 16, 520–525 (2020).

    CAS  Google Scholar 

  19. Wang, L. et al. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861–866 (2020).

    CAS  PubMed  Google Scholar 

  20. Regan, E. C. et al. Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. Nature 579, 359–363 (2020).

    ADS  CAS  PubMed  Google Scholar 

  21. Tang, Y. et al. Simulation of Hubbard model physics in WSe2/WS2 moiré superlattices. Nature 579, 353–358 (2020).

    ADS  CAS  PubMed  Google Scholar 

  22. Nandkishore, R. & Levitov, L. Dynamical screening and excitonic instability in bilayer graphene. Phys. Rev. Lett. 104, 156803 (2010).

    ADS  PubMed  Google Scholar 

  23. 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).

    ADS  CAS  Google Scholar 

  24. Fu, L. Parity-breaking phases of spin–orbit-coupled metals with gyrotropic, ferroelectric, and multipolar orders. Phys. Rev. Lett. 115, 026401 (2015).

    ADS  PubMed  Google Scholar 

  25. Fernandes, R. M. & Venderbos, J. W. Nematicity with a twist: rotational symmetry breaking in a moiré superlattice. Sci. Adv. 6, eaba8834 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  26. 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).

    ADS  PubMed  Google Scholar 

  27. Mishra, A. & Lee, S. Topological multiferroic phases in the extended Kane–Mele–Hubbard model in the Hofstadter regime. Phys. Rev. B 98, 235124 (2018).

    ADS  CAS  Google Scholar 

  28. Cao, Y. et al. Nematicity and competing orders in superconducting magic-angle graphene. Preprint at https://arxiv.org/abs/2004.04148 (2020).

  29. Jiang, Y. et al. Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene. Nature 573, 91–95 (2019).

    ADS  CAS  PubMed  Google Scholar 

  30. 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).

    CAS  Google Scholar 

  31. 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).

    ADS  CAS  Google Scholar 

  32. McCann, E. & Koshino, M. The electronic properties of bilayer graphene. Rep. Prog. Phys. 76, 056503 (2013).

    ADS  PubMed  Google Scholar 

  33. 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).

    CAS  Google Scholar 

  34. Ju, L. et al. Topological valley transport at bilayer graphene domain walls. Nature 520, 650–655 (2015).

    ADS  CAS  PubMed  Google Scholar 

  35. Sui, M. et al. Gate-tunable topological valley transport in bilayer graphene. Nat. Phys. 11, 1027–1031 (2015).

    CAS  Google Scholar 

  36. 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).

    CAS  Google Scholar 

  37. Ju, L. et al. Tunable excitons in bilayer graphene. Science 358, 907–910 (2017).

    ADS  CAS  PubMed  Google Scholar 

  38. Maher, P. et al. Evidence for a spin phase transition at charge neutrality in bilayer graphene. Nat. Phys. 9, 154–158 (2013).

    CAS  Google Scholar 

  39. Hunt, B. et al. Direct measurement of discrete valley and orbital quantum numbers in bilayer graphene. Nat. Commun. 8, 948 (2017).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  40. 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).

    ADS  CAS  PubMed  Google Scholar 

  41. Bao, W. et al. Evidence for a spontaneous gapped state in ultraclean bilayer graphene. Proc. Natl Acad. Sci. USA 109, 10802–10805 (2012).

    ADS  CAS  PubMed  Google Scholar 

  42. Freitag, F., Trbovic, J., Weiss, M. & Schönenberger, C. Spontaneously gapped ground state in suspended bilayer graphene. Phys. Rev. Lett. 108, 076602 (2012).

    ADS  CAS  PubMed  Google Scholar 

  43. 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).

    ADS  CAS  PubMed  Google Scholar 

  44. Fei, Z. et al. Ferroelectric switching of a two-dimensional metal. Nature 560, 336–339 (2018).

    ADS  CAS  PubMed  Google Scholar 

  45. 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).

  46. 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).

    ADS  CAS  Google Scholar 

  47. Zhang, Y. et al. Direct observation of a widely tunable bandgap in bilayer graphene. Nature 459, 820–823 (2009).

    ADS  CAS  PubMed  Google Scholar 

  48. Young, A. F. & Levitov, L. S. Capacitance of graphene bilayer as a probe of layer-specific properties. Phys. Rev. B 84, 085441 (2011).

    ADS  Google Scholar 

  49. 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).

    ADS  CAS  PubMed  Google Scholar 

  50. Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012).

    CAS  Google Scholar 

  51. Dean, C. R. et al. Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013).

    ADS  CAS  PubMed  Google Scholar 

  52. Ponomarenko, L. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013).

    ADS  CAS  PubMed  Google Scholar 

  53. Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013).

    ADS  CAS  PubMed  Google Scholar 

  54. Finney, N. R. et al. Tunable crystal symmetry in graphene–boron nitride heterostructures with coexisting moiré superlattices. Nat. Nanotechnol. 14, 1029–1034 (2019).

    ADS  CAS  PubMed  Google Scholar 

  55. Novoselov, K. S. et al. Unconventional quantum Hall effect and Berry’s phase of 2π in bilayer graphene. Nat. Phys. 2, 177–180 (2006).

    Google Scholar 

  56. Craciun, M. et al. Trilayer graphene is a semimetal with a gate-tunable band overlap. Nat. Nanotechnol. 4, 383–388 (2009).

    ADS  CAS  PubMed  Google Scholar 

  57. Jhang, S. H. et al. Stacking-order dependent transport properties of trilayer graphene. Phys. Rev. B 84, 161408 (2011).

    ADS  Google Scholar 

  58. Wang, H., Wu, Y., Cong, C., Shang, J. & Yu, T. Hysteresis of electronic transport in graphene transistors. ACS Nano 4, 7221–7228 (2010).

    CAS  PubMed  Google Scholar 

  59. McGilly, L. et al. Visualization of moiré superlattices. Nat. Nanotechnol. 15, 580–584 (2020).

    ADS  CAS  PubMed  Google Scholar 

Download references

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

Authors

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

Correspondence to Qiong Ma or Pablo Jarillo-Herrero.

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.

Extended Data Table 1 Device parameters and characteristics for devices M1, N0, T1 and H1–H4

Supplementary information

Supplementary Information

This file contains Supplementary Notes and Discussions and includes Supplementary Figures 1–33.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-020-2970-9

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing