Abstract
Quantum particles on a lattice with competing long-range interactions are ubiquitous in physics; transition metal oxides1,2, layered molecular crystals3 and trapped-ion arrays4 are a few examples. In the strongly interacting regime, these systems often show a rich variety of quantum many-body ground states that challenge theory2. The emergence of transition metal dichalcogenide moiré superlattices provides a highly controllable platform in which to study long-range electronic correlations5,6,7,8,9,10,11,12. Here we report an observation of nearly two dozen correlated insulating states at fractional fillings of tungsten diselenide/tungsten disulfide moiré superlattices. This finding is enabled by a new optical sensing technique that is based on the sensitivity to the dielectric environment of the exciton excited states in a single-layer semiconductor of tungsten diselenide. The cascade of insulating states shows an energy ordering that is nearly symmetric about a filling factor of half a particle per superlattice site. We propose a series of charge-ordered states at commensurate filling fractions that range from generalized Wigner crystals7 to charge density waves. Our study lays the groundwork for using moiré superlattices to simulate a wealth of quantum many-body problems that are described by the two-dimensional extended Hubbard model3,13,14 or spin models with long-range chargeâcharge and exchange interactions15,16.
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
Imada, M., Fujimori, A. & Tokura, Y. Metalâinsulator transitions. Rev. Mod. Phys. 70, 1039â1263 (1998).
Dagotto, E. Complexity in strongly correlated electronic systems. Science 309, 257â262 (2005).
Hotta, C. Theories on frustrated electrons in two-dimensional organic solids. Crystals 2, 1155â1200 (2012).
Lahaye, T., Menotti, C., Santos, L., Lewenstein, M. & Pfau, T. The physics of dipolar bosonic quantum gases. Rep. Prog. Phys. 72, 126401 (2009).
Zhu, Q., Tu, M. W. Y., Tong, Q. & Yao, W. Gate tuning from exciton superfluid to quantum anomalous Hall in van der Waals heterobilayer. Sci. Adv. 5, eaau6120 (2019).
Tang, Y. et al. Simulation of Hubbard model physics in WSe2/WS2 moiré superlattices. Nature 579, 353â358 (2020).
Regan, E. C. et al. Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. Nature 579, 359â363 (2020).
Wu, F., Lovorn, T., Tutuc, E. & Macdonald, A. H. Hubbard model physics in transition metal dichalcogenide moiré bands. Phys. Rev. Lett. 121, 026402 (2018).
Wu, F., Lovorn, T., Tutuc, E., Martin, I. & Macdonald, A. H. Topological insulators in twisted transition metal dichalcogenide homobilayers. Phys. Rev. Lett. 122, 086402 (2019).
Shimazaki, Y. et al. Strongly correlated electrons and hybrid excitons in a moiré heterostructure. Nature 580, 472â477 (2020).
Wang, L. et al. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861â866 (2020).
Zhang, Z. et al. Flat bands in twisted bilayer transition metal dichalcogenides. Nat. Phys. https://doi.org/10.1038/s41567-020-0958-x (2020).
Pietig, R., Bulla, R. & Blawid, S. Reentrant charge order transition in the extended Hubbard model. Phys. Rev. Lett. 82, 4046â4049 (1999).
Tocchio, L. F., Gros, C., Zhang, X. F. & Eggert, S. Phase diagram of the triangular extended Hubbard model. Phys. Rev. Lett. 113, 246405 (2014).
McKenzie, R. H., Merino, J., Marston, J. B. & Sushkov, O. P. Charge ordering and antiferromagnetic exchange in layered molecular crystals of the θ type. Phys. Rev. B 64, 085109 (2001).
Porras, D. & Cirac, J. I. Quantum manipulation of trapped ions in two dimensional coulomb crystals. Phys. Rev. Lett. 96, 250501 (2006).
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).
Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059â1064 (2019).
Serlin, M. et al. Intrinsic quantized anomalous Hall effect in a moiré heterostructure. Science 367, 900â903 (2020).
Lu, X. et al. Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene. Nature 574, 653â657 (2019).
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); correction 581, E3 (2020).
Padhi, B., Setty, C. & Phillips, P. W. Doped twisted bilayer graphene near magic angles: proximity to Wigner crystallization, not Mott insulation. Nano Lett. 18, 6175â6180 (2018).
Wang, G. et al. Excitons in atomically thin transition metal dichalcogenides. Rev. Mod. Phys. 90, 021001 (2018).
Raja, A. et al. Coulomb engineering of the bandgap and excitons in two-dimensional materials. Nat. Commun. 8, 15251 (2017).
Sidler, M. et al. Fermi polaron-polaritons in charge-tunable atomically thin semiconductors. Nat. Phys. 13, 255â261 (2017).
Gao, S., Liang, Y., Spataru, C. D. & Yang, L. Dynamical excitonic effects in doped two-dimensional semiconductors. Nano Lett. 16, 5568â5573 (2016).
Jin, C. et al. Observation of moiré excitons in WSe2/WS2 heterostructure superlattices. Nature 567, 76â80 (2019); correction 569, E7 (2019).
Tran, K. et al. Evidence for moiré excitons in van der Waals heterostructures. Nature 567, 71â75 (2019).
Alexeev, E. M. et al. Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures. Nature 567, 81â86 (2019); correction 572, E8 (2019).
Seyler, K. L. et al. Signatures of moiré-trapped valley excitons in MoSe2/WSe2 heterobilayers. Nature 567, 66â70 (2019).
Zhang, Y., Yuan, N. F. Q. & Fu, L. Moiré quantum chemistry: charge transfer in transition metal dichalcogenide superlattices. Preprint at https://arxiv.org/abs/1910.14061 (2019).
Slagle, K. & Fu, L. Charge transfer excitations, pair density waves, and superconductivity in moiré materials. Preprint at https://arxiv.org/abs/2003.13690 (2020).
Martin, I. & Batista, C. D. Itinerant electron-driven chiral magnetic ordering and spontaneous quantum hall effect in triangular lattice models. Phys. Rev. Lett. 101, 156402 (2008).
Nandkishore, R., Chern, G. W. & Chubukov, A. V. Itinerant half-metal spin-density-wave state on the hexagonal lattice. Phys. Rev. Lett. 108, 227204 (2012).
Nandkishore, R., Thomale, R. & Chubukov, A. V. Superconductivity from weak repulsion in hexagonal lattice systems. Phys. Rev. B 89, 144501 (2014).
Wang, L. et al. One-dimensional electrical contact to a two-dimensional material. Science 342, 614â617 (2013).
Chernikov, A. et al. Exciton binding energy and nonhydrogenic Rydberg series in monolayer WS2. Phys. Rev. Lett. 113, 076802 (2014).
Stier, A. V. et al. Magnetooptics of exciton Rydberg states in a monolayer semiconductor. Phys. Rev. Lett. 120, 057405 (2018).
He, K. et al. Tightly bound excitons in monolayer WSe2. Phys. Rev. Lett. 113, 026803 (2014).
Efimkin, D. K. & MacDonald, A. H. Many-body theory of trion absorption features in two-dimensional semiconductors. Phys. Rev. B 95, 035417 (2017).
Movva, H. C. P. et al. Density-dependent quantum Hall states and Zeeman splitting in monolayer and bilayer WSe2. Phys. Rev. Lett. 118, 247701 (2017).
Acknowledgements
We thank L. Fu for fruitful discussions. Research was primarily supported by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award number DE-SC0019481 (optical spectroscopy and growth of WSe2 crystals). Device fabrication was supported by US Office of Naval Research under award number N00014-18-1-2368. The growth of the hBN crystals was supported by the Elemental Strategy Initiative of MEXT, Japan and CREST (JPMJCR15F3), JST. K.F.M. also acknowledges support from David and Lucille Packard Fellowship.
Author information
Authors and Affiliations
Contributions
Y.X. fabricated the devices, performed the measurements and analysed the data. V.E. performed theoretical calculations. S.L., D.A.R. and J.H. grew the bulk WSe2 crystals, and K.W. and T.T. grew the bulk hBN crystals. Y.X., K.F.M. and J.S. designed the scientific objectives and oversaw the project. Y.X., V.E., K.F.M. and J.S. co-wrote the manuscript. All authors discussed the results and commented on the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature thanks Xiaobo Lu, Fengcheng Wu and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
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 Additional doping-dependence data and analysis for the control device.
a, Gate-dependent reflection contrast (ÎR/R0) spectrum near the WSe2 (sensor) 2s transition energy (same data as in Fig. 1d, middle). The WS2 sample is electron-doped above Vg = 1.5 V (white dashed line). b, Representative linecuts of a. Curves are shifted vertically for clarity. The red and blue dashed curves highlight the redshifted and blueshifted branches, respectively. c, Spectral weight of the two branches as a function of Vg (calculated after removing the background taken at 1 V and 4 V for the redshifted and blueshifted branches, respectively). The blueshifted branch has been multiplied by a factor of 0.3. The redshifted feature emerges while the blueshifted feature quickly goes to zero at Vg â 1.5 V.
Extended Data Fig. 2 Assignment of the filling factor of the insulating states.
For each insulating state, the position (filled squares) and the FWHM (horizontal bars) in gate voltage are determined through a Lorentzian fit of the 2s exciton resonance energy as a function of gate voltage. The red lines are linear fits using four established states v = 2, v = 1, v = 2/3 and v = 1/3 on the electron and hole sides independently. Both slopes correspond to 0.25 filling per volt. The filling factor of the other states is assigned using the slope as described in Methods.
Extended Data Fig. 3 Results from different regions of the device.
aâe, Gate-dependent reflection contrast spectrum measured at different regions of the device P1 (a), P2 (b), P3 (c), P4 (d) and P5 (e) at 1.6 K. All share the same x axis, given at the bottom. The correlated states that can be identified are labelled on the top axes.
Extended Data Fig. 4 Results from a different device.
a, Gate-dependent reflection contrast spectrum of a different sample as descried in the Methods at 1.6 K. The filling factor is shown on the right axis. A fixed back gate voltage of 8 V is applied to dope the contact region to form good electrical contact between the TMD moiré superlattice and the contact electrode. The top gate voltage Vg is swept to tune the electron doping density. Below Vg = â9 V (dashed white line), the sample is charge neutral. b, As in a, focusing on the 2s transition in the sensor. The identified insulating states are marked with their corresponding filling factors.
Extended Data Fig. 5 Analysis of the â1/3 and â2/3 states.
a, b, Two horizontal line cuts of Fig. 2a at 1.8403 eV (a) and 1.8383 eV (b). These energies correspond to the 2s exciton peak energy for the v = â1/3 and v = â2/3 states, respectively. They appear as peaks in âÎR/R0. Results at different temperatures (in ascending order from bottom to top, 1.6 K, 4 K, 6 K, 10 K, 13 K, 15 K, 17 K, 19 K, 22 K, 26 K, 29 K, 32 K, 35 K, 38 K, 40 K and 50 K) are vertically displaced for clarity. The red curves are Lorentzian fits to the peaks for the corresponding states. c, Reflection contrast spectrum (âÎR/R0) for Vg = â 3.3 V (v = â2/3) at 1.6 K. The red area underneath the 2s peak is integrated to obtain the 2s spectral weight. d, The 2s spectral weight as a function of temperature for state v = â1/3 (black symbols) and v = â2/3 (red symbols). The v = â1/3 state has a slightly higher TC than the v = â2/3 state. The lines are guides to the eye.
Extended Data Fig. 6 First derivative of data shown in Fig. 2a with respect to energy.
This is used to evaluate the gate voltage and width of the less well-developed insulating states such as those with |v|Â >Â 1. The black arrows highlight the enhanced features for the 1/7 and 6/7 states.
Extended Data Fig. 7 Optical response of the WSe2/WS2 moiré superlattice.
a, b, Gate-dependent reflection contrast spectrum in regions of the device without (a) and with (b) the WSe2 sensor at 1.6Â K. The Fermi level is inside the WSe2/WS2 bandgap between the dashed lines.
Extended Data Fig. 8 Contour plots of additional gate-dependent reflection contrast spectrum at higher temperatures.
From top to bottom, T = 67 K, 80 K and 150 K. All share the same x axis, given at bottom.
Extended Data Fig. 9 Transition temperature to the charge-ordered state (simulation).
Transition temperatures are determined from the peak in the temperature dependence of the heat capacity Cp for fillings 1/7, 1/4, 1/3, 2/5 and 1/2. Temperature is given in units of the energy e2/(4Ïεε0a) and d/a is fixed to be 10.
Supplementary information
Rights and permissions
About this article
Cite this article
Xu, Y., Liu, S., Rhodes, D.A. et al. Correlated insulating states at fractional fillings of moiré superlattices. Nature 587, 214â218 (2020). https://doi.org/10.1038/s41586-020-2868-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-020-2868-6
This article is cited by
-
Dipolar many-body complexes and their interactions in stacked 2D heterobilayers
Nature Reviews Physics (2024)
-
Ultrafast atomic-scale scanning tunnelling spectroscopy of a single vacancy in a monolayer crystal
Nature Photonics (2024)
-
Observation of spin polarons in a frustrated moiré Hubbard system
Nature Physics (2024)
-
Excitonic Mott insulator in a Bose-Fermi-Hubbard system of moiré WS2/WSe2 heterobilayer
Nature Communications (2024)
-
Remote imprinting of moiré lattices
Nature Materials (2024)