ACS Publications. Most Trusted. Most Cited. Most Read
Terahertz Fingerprint of Monolayer Wigner Crystals
My Activity
  • Open Access
Letter

Terahertz Fingerprint of Monolayer Wigner Crystals
Click to copy article linkArticle link copied!

Open PDFSupporting Information (1)

Nano Letters

Cite this: Nano Lett. 2022, 22, 3, 1311–1315
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.nanolett.1c04620
Published January 20, 2022

Copyright © 2022 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Abstract

Click to copy section linkSection link copied!

The strong Coulomb interaction in monolayer semiconductors represents a unique opportunity for the realization of Wigner crystals without external magnetic fields. In this work, we predict that the formation of monolayer Wigner crystals can be detected by their terahertz response spectrum, which exhibits a characteristic sequence of internal optical transitions. We apply the density matrix formalism to derive the internal quantum structure and the optical conductivity of the Wigner crystal and to microscopically analyze the multipeak shape of the obtained terahertz spectrum. Moreover, we predict a characteristic shift of the peak position as a function of charge density for different atomically thin materials and show how our results can be generalized to an arbitrary two-dimensional system.

This publication is licensed under

CC-BY 4.0 .
  • cc licence
  • by licence
Copyright © 2022 The Authors. Published by American Chemical Society
At low temperatures, electrons arrange in a crystal lattice (1) in order to minimize their repulsive energy, which represents one of the most intriguing quantum phase transitions. Ever since their prediction, the experimental realization of these Wigner crystals (WCs) has remained challenging, considering that the prerequisite of a dominating Coulomb versus kinetic energy requires very low temperatures and charge densities. Therefore, the first realizations of WCs (2−4) have used strong magnetic fields and quasi two-dimensional (2D) quantum well systems in order to quench the kinetic energy via Landau quantization. The advances in the research on transition-metal dichalcogenide (TMD) monolayers (5−7) have delivered new opportunities to realize electron crystallization. Several recent studies (8−12) have reported insulating states at fractional fillings of flat moiré bands, (13−15) which arise in the superlattices of twisted homobilayer and heterobilayer systems. However, similar to the application of magnetic fields, these so-called generalized Wigner crystals heavily rely on the external moiré potential stabilizing the lattice structure and can therefore not be considered as intrinsic electronic phases.
In contrast, TMDs in their monolayer form already represent a unique platform to achieve electron crystallization, as they exhibit large effective masses and a strong Coulomb interaction leading to prominent exciton physics. (16,17) Smoleński et al. have recently demonstrated signatures of a WC formed in hBN-encapsulated MoSe2 monolayers (18) at electron densities of up to n = 3 × 1011 cm–2 without any external fields. As an indicator for the crystallization, the authors have used the emergence of an additional peak in the reflectance spectrum at visible frequencies. This peak stems from excitons with finite center-of-mass momentum that couple to the light cone via Bragg scattering at the WC. Other rather indirect confirmations of the WC, such as a decrease in the conductivity, also have been applied to verify other WC systems. (2−4,8,10−12) However, a direct, more reliable and conclusive experimental proof of the WC, such as a direct spectral signature of the internal structure of the crystal, could not be provided so far.
In this work, we propose to exploit terahertz radiation to directly probe the internal quantum transitions of the WC in TMD monolayers, as illustrated in Figure 1. We have developed a fully quantum mechanical model of the interaction of monolayer WC with low frequency light. Our model is based on the many-particle density matrix formalism, where the Hamiltonian of interacting electrons is transformed to the self-consistent eigenbasis of the WC on a Hartree–Fock level. Consequently, we analyze the linear optical response of the system resulting from transitions between different bands of the electron crystal.

Figure 1

Figure 1. Sketch of the 2D Wigner crystal with a honeycomb lattice and alternating spin polarization. The colored curves underneath the particles illustrate their wave functions. A low-energy light field can probe the internal transition of the Wigner crystal toward excited quantum states to provide a conclusive proof of the phase transition.

We predict the emergence of a series of terahertz resonances whose position is a function of electron density and its characteristic peak sequence can be explained by the nature of the involved excited quantum states. Consequently, the predicted low frequency response can be used as a direct probe of the Wigner crystallization in arbitrary 2D systems of interacting fermions. We apply our microscopic model to determine the material realistic response of hBN-encapsulated WSe2 monolayers without any open parameters and explain the characteristic response shape with the help of our microscopic model. We predict that the central energy of the internal quantum resonances shifts from ∼4 meV at 109 cm–2 to 24 meV at 8 × 1010 cm–2. Moreover, we demonstrate that an enhanced effective mass and reduced dielectric screening, such as in MoSe2 on SiO2, can push the WC resonances up to 100 meV. Finally, we also present a unitless study for ideal 2D systems, which can be rescaled to arbitrary fermionic systems.

Microscopic Mean Field Theory

Previous numerical studies on WCs based on quantum Monte Carlo methods (19) have demonstrated that the low-density crystal phase is characterized by small correlation energies and that the WCs properties are well-reproduced within a Hartree–Fock description of the Coulomb interaction. (20) To obtain access to the optical response of the WC, we start with the Hamilton operator of interacting electrons,
(1)
which can, e.g., be related to the free charge carriers within n(p)-doped TMD monolayers, such that aσk(†) creates/annihilates electrons (holes) in the lowest conduction band (highest valence band) with spin up at the K point (σ = ) or spin down at the −K point (σ = ) and momentum k. For low temperature and densities, the electron dispersion εk is determined by an effective mass m*. (21) The Coulomb interaction is treated in a mean field approximation (20,22) via the nonlocal Hartree–Fock potential,
(2)
depending on the Fourier transform of the Coulomb potential Vq, which includes the dielectric environment of the TMD monolayer through the Rytova–Keldysh potential. (23,24) Moreover, the mean field potential also is dependent on the charge distribution entering via the Fourier transform of the Wigner distribution wkσ(q) = ⟨aσk+qaσk⟩. The Hamiltonian in eq 1 is now diagonalized via the basis transformation Aλσk = ∑Guσkλ*(G)aσ,k+G, where we sum over the reciprocal lattice vectors G of the WC. Here, uσkλ(G) is the complete set of Bloch functions fulfilling the periodic Hartree–Fock equation (see the Supporting Information (SI)). The latter is solved iteratively until a self-consistent solution is found, starting with a randomized charge distribution. Hence, in the new basis describing WC electrons the Hamiltonian reads
(3)
where the self-consistent solution requires that the lowest band of the WC (λ = 0) is completely filled (one electron per WC unit cell) and the rest is empty, i.e., fλσk = ⟨AλσkAλσk⟩ = δλ0.
Having determined the energy spectrum of WC electrons Eσkλ and their wave functions uσkλ(G), we can transform the electron-light interaction into the Wigner basis and compute the dynamics of the WC using the Heisenberg’s equation of motion. (25,26) In particular, assuming a weak perturbation through the electromagnetic field, we can derive the linear response of the system, which, e.g., can be specified via the diagonal components of the linear optical conductivity tensor reading:
(4)
where we have introduced the monolayer thickness d, (27) the normalization area A, and a phenomenological dephasing constant Γ. The oscillator strength of transitions between different WC sub-bands with energy ΔEσkλ = EσkλEσk0 for ê-polarized light is given by the current matrix element Jσkλ=e0/m*Guσkλ*(G)G·êuσk0(G). All details of the derivation and the used parameters are given in the SI. The microscopic model described above is not specific to TMD monolayers and can be applied to arbitrary 2D electron system. In principle, a similar approach can be used to compute the optical response also in generalized WCs within twisted bilayers.

Internal Quantum Structure of WCs

Now, we apply the developed model to the exemplary system of p-doped, hBN-encapsulated WSe2 monolayers at a density of n ≈ 3 × 1010 cm–2 (rs = 50), i.e., deep within the crystalline phase that was predicted to melt at rs = 30–40. (28−30) The temperature is assumed low enough such that only the lowest band of the WC is occupied and the thermal occupation of the excited states can be neglected. Moreover, we focus on the spin-unpolarized (antiferromagnetic) phase, since we assume that the charge density is controlled via an electrostatic gate. As a result, the Wigner crystal is in a chemical equilibrium with an unpolarized charge reservoir. Apart from the WC geometry, the main conclusions qualitatively also hold for spin (valley)-polarized systems. (31)Figure 2 shows the computed band structure for electrons with σ = (on a cut through the WCs Brillouin zone) as well as the wave functions of the monolayer WC. The band structure is illustrated in Figure 2a, where the background has been color-coded based on the optical activity,
We find a completely flat ground state at ∼17 meV below the edge of the single particle dispersion. We also observe a series of flat excited states, where the first is separated from the ground state by a Wigner gap of ∼10 meV. With the increasing band index, the excited states become more dense and dispersed. In order to interpret this result, we investigate the underlying nature of these quantum states and their respective wave functions.

Figure 2

Figure 2. Internal quantum structure of the WC in p-doped, hBN-encapsulated WSe2 at rs = 50. (a) Band structure of the WC, where only the lowest band is fully occupied. The background color illustrates the oscillator strength for optical transitions from the ground state. (b) The surface plot shows the charge density, whereas the projection contains the spin polarization ( is red, is blue). Also shown are the wave function at k = 0 of the (c) first excited state and (d) the state with the largest oscillator.

The fully occupied ground state is, in fact, the WC itself, i.e., the corresponding wave functions are eigenstates of the Coulomb potential created by their own collective charge distribution. Here, the charge density (Figure 2b) exhibits strongly localized peaks ordered in a honeycomb lattice (32) with alternating spin polarization (red and blue projections). Moreover, the localized charge distributions exhibit a significant triangular warping, which is usually neglected in variational approaches.
In contrast to the self-consistent ground state, the excited states are unoccupied and therefore do not contribute to the nonlocal Hartree–Fock potential (see the SI). Figure 2c shows the wave function of the first excited state at k = 0, whereas the contours indicate the peaks of the (red) and charge density (blue) from Figure 2b. In order to minimize the repulsive Coulomb energy, the excited state is strongly localized in the pockets between the charge peaks. Consequently, the wave function overlap with the ground state is small, which is reflected by the weak optical activity of the first excited state. With the increasing band index, we first find several other localized states with increasing kinetic energy and growing spatial extent. (15) Consequently, the localized orbitals of higher-order bands begin to overlap at higher quantum indices and we find dispersed bands at high energies corresponding to scattering states. Figure 2d shows the excited state with the largest optical activity for transitions from the occupied ground state. As result of its delocalized character and the beneficial exchange interaction between equally spin-polarized electrons, the overlap with the ground state (red contours) is large, resulting in a strong current matrix element. The proposed detection mechanism via low-energy excitations has previously been applied to WCs formed in GaAs quantum wells in strong magnetic fields, where a single peak has been observed at microwave energies, (33,34) the so-called pinning frequency. (35) In contrast, our microscopic model predicts that zero-field monolayer Wigner crystals should exhibit an entire series of resonances at much larger frequencies.
With the knowledge of the internal WC quantum structure, we can now consider the dynamics of the system under external perturbations, e.g., the absorption of light.

Optical Fingerprint of WCs

Now, we evaluate the linear optical conductivity from eq 4 in the energy range of the WCs interband transition. Figure 3 illustrates the real part of the resulting spectrum for two exemplary line widths. For the narrower line width of Γ = 0.1 meV, the spectrum exhibits a distinguished multipeak structure with a characteristic sequence of peak amplitudes. The series of resonances starts at the energy of the Wigner gap (cf. Figure 2a), reflecting the transition to the first excited state, which has a weak oscillator strength, as discussed above. The subsequent series of peaks corresponds to the transition to higher-order excited states, whose oscillator strength increases until a maximum is reached. This is a result of the growing kinetic energy of the excited states with increasing band index. With the increasing orbital size of the excited states, which are mostly localized in the pockets between the WC electrons, their overlap with the ground-state electrons becomes enhanced. After a maximum is reached for the excited state illustrated in Figure 2d, the oscillator strength decreases again. Here, the final states momentum spectrum uσkλ(G) begins to move toward large G vectors, corresponding to scattering states with a large kinetic energy. Consequently, the momentum space overlap determining Jσkλ shrinks.

Figure 3

Figure 3. Terahertz fingerprint of the Wigner crystal. The real part of the optical conductivity exhibits a characteristic series of internal quantum transitions, whose amplitude sequence can be explained by the nature of the involved excited states. Note that the strongly broadened response signal has been scaled by a factor of 2, whereas the Drude signal of free electrons (black line) has been multiplied by 100.

The calculated conductivity, on the order of 10 Ω–1 cm–1, corresponds to a rather low overall absorbance of the monolayer in the range of α = dRe(σ)/(cϵ0) ≈ 10–4, which is due to the nanometer thickness of the sample. However, comparably weak low-frequency responses stemming from intraexcitonic transitions have already been used to measure time-resolved exciton dynamics in TMD monolayers. (36)
For the larger line width of  Γ = 1.5 meV, the spectrum melts to a single peak centered at ∼15 meV. For comparison, we have also included the Drude response that is expected for an electron gas at the same density. This shows that, even for broad resonances, the terahertz spectrum is undergoing a dramatic change during the formation of the WC. However, we expect a much smaller line width than Γ = 1.5 meV, since the crystallization is occurring at very low temperatures and the phase space for electron–phonon scattering is strongly suppressed within the flat bands of the WC. Moreover, in order to ensure the stability of the WC, ultraclean monolayer samples that are almost free of impurities and disorder are necessary. These samples will, at the same time, also reduce inhomogeneous broadening of WC resonances. Under these conditions, we predict a clear multipeak shape that will allow one to distinguish between different transitions the internal quantum states of the WC.

Density Dependence

Now, we vary the charge density and investigate how the spectral response shifts as the lattice parameter of the WC, , is changed. To this end, we consider the center of the spectral response defined as ωpeak = ∫ωσ(ω)dω/∫σ(ω)dω, where the integrals are only performed over positive frequencies. Figure 4 shows the center energy of the optical response as a function of density for different TMD monolayers. The filled dots represent the values at rs = 30, such that for higher densities quantum melting of the WC is expected. For all three illustrated systems we consider p-doping, since the valence band in TMDs has a significantly larger mass then the conduction band, which enables larger critical densities before quantum melting sets in. Figure 4 illustrates that the optical response shifts significantly to higher energy at enhanced charge densities, which can serve as a clear indicator that the observed peaks indeed stem from transitions of the WC. With increasing density, the charge density peaks are more closely spaced, which enhances the repulsive Coulomb energy for the excited states and therefore increases the transition energies. As a result, we also find that the transition energies for samples on SiO2 are significantly larger, since the Coulomb interaction is much less screened by the substrate. In addition to the spectral shift of the center energy, the amplitude of the optical response also increases with growing densities, since the density of optically active particles is enhanced.

Figure 4

Figure 4. Central energy of the optical response signal as a function of the hole density for different monolayers/substrates. The inset shows the result for an ideal 2D system, where the only system parameter is given by rs and the energy scale is given in units of effective Hartree.

We have additionally performed calculations for an ideal 2D system with Vq ∝ 1/q, which reproduces the results for TMD monolayers at small densities. For an ideal 2D system, the only relevant system parameter is given by with the effective Bohr radius being defined as aB = 4πϵ0ϵ2/(m*e02). Using this density scale, the energies can be given in terms of effective Hartree EH=2/(m*aB2), such that the results given in the inset of Figure 4 can be rescaled to arbitrary 2D system. The inset furthermore demonstrates that the overall width of the optical response ΔEEpeakEgap also grows with increasing density, since the distance between excited states becomes larger.
Finally, we want to note that the susceptibility derived in our work results from an oscillating polarization current as a response to the electromagnetic field. Since the oscillating electrons (holes) are forming the WC crystal itself, the optical resonances found in our work could be interpreted in a classical picture as the excitation of long-range phonon modes of the WC. (32) Moreover, we want to emphasize that the applied model is only valid for weak excitation conditions, i.e., assuming a negligible population of the excited states. A lower bound for the critical power density of the exciting laser can be obtained by estimating a lower bound for the lifetime of excited states. Assuming that the lifetime is shorter than τ = 1 ns, the weak excitation regime can be specified by a power density of In(E1E0)/τ, which for the systems studied in Figures 2 and 3, yields I ≪ 50 mW/cm2. For stronger excitation powers, we expect optical nonlinearities to occur that could even result in a destruction of the WC.
In conclusion, our study demonstrates that the terahertz response of monolayer Wigner crystals can be used as an unambiguous fingerprint, exhibiting a characteristic sequence of internal quantum transitions. The inherent spectral shift of resonance frequencies with the charge density further represents a strong benchmark for the crystallization. We have predicted the peak positions for different TMD monolayers, as well as other 2D fermion systems. The presented results will guide future experiments toward the detection of Wigner crystallization and the developed approach can be further exploited to theoretically study the interaction dynamics in pure as well as generalized Wigner crystals in twisted bilayers.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.nanolett.1c04620.

  • Transformation of Hartree–Fock Hamiltonian into Wigner basis; material-specific parameters and interaction matrix elements; derivation of the optical response (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
    • Notes
      The authors declare no competing financial interest.

    Acknowledgments

    Click to copy section linkSection link copied!

    We acknowledge support from Deutsche Forschungsgemeinschaft (DFG) via SFB 1083 (Project B9) and the European Unions Horizon 2020 Research and Innovation Program, under Grant Agreement No. 881603 (Graphene Flagship).

    References

    Click to copy section linkSection link copied!

    This article references 36 other publications.

    1. 1
      Wigner, E. On the interaction of electrons in metals. Phys. Rev. 1934, 46, 1002,  DOI: 10.1103/PhysRev.46.1002
    2. 2
      Andrei, E.; Deville, G.; Glattli, D.; Williams, F.; Paris, E.; Etienne, B. Observation of a magnetically induced Wigner solid. Phys. Rev. Lett. 1988, 60, 2765,  DOI: 10.1103/PhysRevLett.60.2765
    3. 3
      Goldman, V.; Santos, M.; Shayegan, M.; Cunningham, J. Evidence for two-dimentional quantum Wigner crystal. Phys. Rev. Lett. 1990, 65, 2189,  DOI: 10.1103/PhysRevLett.65.2189
    4. 4
      Williams, F.; Wright, P.; Clark, R.; Andrei, E.; Deville, G.; Glattli, D.; Probst, O.; Etienne, B.; Dorin, C.; Foxon, C. Conduction threshold and pinning frequency of magnetically induced Wigner solid. Phys. Rev. Lett. 1991, 66, 3285,  DOI: 10.1103/PhysRevLett.66.3285
    5. 5
      Wang, G.; Chernikov, A.; Glazov, M. M.; Heinz, T. F.; Marie, X.; Amand, T.; Urbaszek, B. Colloquium: Excitons in atomically thin transition metal dichalcogenides. Rev. Mod. Phys. 2018, 90, 021001,  DOI: 10.1103/RevModPhys.90.021001
    6. 6
      Mueller, T.; Malic, E. Exciton physics and device application of two-dimensional transition metal dichalcogenide semiconductors. npj 2D Mater. Appl. 2018, 2, 112,  DOI: 10.1038/s41699-018-0074-2
    7. 7
      Raja, A.; Waldecker, L.; Zipfel, J.; Cho, Y.; Brem, S.; Ziegler, J. D.; Kulig, M.; Taniguchi, T.; Watanabe, K.; Malic, E. Dielectric disorder in two-dimensional materials. Nat. Nanotechnol. 2019, 14, 832837,  DOI: 10.1038/s41565-019-0520-0
    8. 8
      Xu, Y.; Liu, S.; Rhodes, D. A.; Watanabe, K.; Taniguchi, T.; Hone, J.; Elser, V.; Mak, K. F.; Shan, J. Correlated insulating states at fractional fillings of moiré superlattices. Nature 2020, 587, 214218,  DOI: 10.1038/s41586-020-2868-6
    9. 9
      Wang, L.; Shih, E.-M.; Ghiotto, A.; Xian, L.; Rhodes, D. A.; Tan, C.; Claassen, M.; Kennes, D. M.; Bai, Y.; Kim, B. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 2020, 19, 861866,  DOI: 10.1038/s41563-020-0708-6
    10. 10
      Zhou, Y.; Sung, J.; Brutschea, E.; Esterlis, I.; Wang, Y.; Scuri, G.; Gelly, R. J.; Heo, H.; Taniguchi, T.; Watanabe, K. Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure. Nature 2021, 595, 4852,  DOI: 10.1038/s41586-021-03560-w
    11. 11
      Huang, X.; Wang, T.; Miao, S.; Wang, C.; Li, Z.; Lian, Z.; Taniguchi, T.; Watanabe, K.; Okamoto, S.; Xiao, D. Correlated insulating states at fractional fillings of the WS2/WSe2 moiré lattice. Nat. Phys. 2021, 17, 715719,  DOI: 10.1038/s41567-021-01171-w
    12. 12
      Miao, S.; Wang, T.; Huang, X.; Chen, D.; Lian, Z.; Wang, C.; Blei, M.; Taniguchi, T.; Watanabe, K.; Tongay, S. Strong interaction between interlayer excitons and correlated electrons in WSe2/WS2 moiré superlattice. Nat. Commun. 2021, 12, 16,  DOI: 10.1038/s41467-021-23732-6
    13. 13
      Yu, H.; Liu, G.-B.; Tang, J.; Xu, X.; Yao, W. Moiré excitons: From programmable quantum emitter arrays to spin-orbit–coupled artificial lattices. Sci. Adv. 2017, 3, e1701696  DOI: 10.1126/sciadv.1701696
    14. 14
      Tran, K.; Moody, G.; Wu, F.; Lu, X.; Choi, J.; Kim, K.; Rai, A.; Sanchez, D. A.; Quan, J.; Singh, A. Evidence for moiré excitons in van der Waals heterostructures. Nature 2019, 567, 7175,  DOI: 10.1038/s41586-019-0975-z
    15. 15
      Brem, S.; Linderälv, C.; Erhart, P.; Malic, E. Tunable phases of moiré excitons in van der Waals heterostructures. Nano Lett. 2020, 20, 85348540,  DOI: 10.1021/acs.nanolett.0c03019
    16. 16
      Ugeda, M. M.; Bradley, A. J.; Shi, S.-F.; da Jornada, F. H.; Zhang, Y.; Qiu, D. Y.; Ruan, W.; Mo, S.-K.; Hussain, Z.; Shen, Z.-X. Giant bandgap renormalization and excitonic effects in a monolayer transition metal dichalcogenide semiconductor. Nat. Mater. 2014, 13, 10911095,  DOI: 10.1038/nmat4061
    17. 17
      Chernikov, A.; Berkelbach, T. C.; Hill, H. M.; Rigosi, A.; Li, Y.; Aslan, O. B.; Reichman, D. R.; Hybertsen, M. S.; Heinz, T. F. Exciton binding energy and nonhydrogenic Rydberg series in monolayer WS2. Phys. Rev. Lett. 2014, 113, 076802,  DOI: 10.1103/PhysRevLett.113.076802
    18. 18
      Smolenski, T.; Dolgirev, P. E.; Kuhlenkamp, C.; Popert, A.; Shimazaki, Y.; Back, P.; Lu, X.; Kroner, M.; Watanabe, K.; Taniguchi, T. Signatures of Wigner crystal of electrons in a monolayer semiconductor. Nature 2021, 595, 5357,  DOI: 10.1038/s41586-021-03590-4
    19. 19
      Drummond, N.; Radnai, Z.; Trail, J.; Towler, M.; Needs, R. Diffusion quantum Monte Carlo study of three-dimensional Wigner crystals. Phys. Rev. B 2004, 69, 085116,  DOI: 10.1103/PhysRevB.69.085116
    20. 20
      Trail, J.; Towler, M.; Needs, R. Unrestricted Hartree-Fock theory of Wigner crystals. Phys. Rev. B 2003, 68, 045107,  DOI: 10.1103/PhysRevB.68.045107
    21. 21
      Kormányos, A.; Burkard, G.; Gmitra, M.; Fabian, J.; Zólyomi, V.; Drummond, N. D.; Fal’ko, V. k-p theory for two-dimensional transition metal dichalcogenide semiconductors. 2D Materials 2015, 2, 022001,  DOI: 10.1088/2053-1583/2/2/022001
    22. 22
      Pan, H.; Wu, F.; Das Sarma, S. Quantum phase diagram of a Moiré-Hubbard model. Phys. Rev. B 2020, 102, 201104,  DOI: 10.1103/PhysRevB.102.201104
    23. 23
      Rytova, N. Screened potential of a point charge in a thin film. Moscow Univ. Phys. Bull. 1967, 3, 30
    24. 24
      Brem, S.; Ekman, A.; Christiansen, D.; Katsch, F.; Selig, M.; Robert, C.; Marie, X.; Urbaszek, B.; Knorr, A.; Malic, E. Phonon-assisted photoluminescence from indirect excitons in monolayers of transition-metal dichalcogenides. Nano Lett. 2020, 20, 28492856,  DOI: 10.1021/acs.nanolett.0c00633
    25. 25
      Lindberg, M.; Koch, S. W. Effective Bloch equations for semiconductors. Phys. Rev. B 1988, 38, 3342,  DOI: 10.1103/PhysRevB.38.3342
    26. 26
      Kira, M.; Koch, S. Many-body correlations and excitonic effects in semiconductor spectroscopy. Progress Quantum Electron. 2006, 30, 155296,  DOI: 10.1016/j.pquantelec.2006.12.002
    27. 27
      Laturia, A.; Van de Put, M. L.; Vandenberghe, W. G. Dielectric properties of hexagonal boron nitride and transition metal dichalcogenides: from monolayer to bulk. npj 2D Mater. Appl. 2018, 2, 17,  DOI: 10.1038/s41699-018-0050-x
    28. 28
      Yoon, J.; Li, C.; Shahar, D.; Tsui, D.; Shayegan, M. Wigner crystallization and metal-insulator transition of two-dimensional holes in GaAs at B = 0. Phys. Rev. Lett. 1999, 82, 1744,  DOI: 10.1103/PhysRevLett.82.1744
    29. 29
      Drummond, N.; Needs, R. Phase diagram of the low-density two-dimensional homogeneous electron gas. Phys. Rev. Lett. 2009, 102, 126402,  DOI: 10.1103/PhysRevLett.102.126402
    30. 30
      Zarenia, M.; Neilson, D.; Partoens, B.; Peeters, F. Wigner crystallization in transition metal dichalcogenides: A new approach to correlation energy. Phys. Rev. B 2017, 95, 115438,  DOI: 10.1103/PhysRevB.95.115438
    31. 31
      Tanatar, B.; Ceperley, D. M. Ground state of the two-dimensional electron gas. Phys. Rev. B 1989, 39, 5005,  DOI: 10.1103/PhysRevB.39.5005
    32. 32
      Bonsall, L.; Maradudin, A. Some static and dynamical properties of a two-dimensional Wigner crystal. Phys. Rev. B 1977, 15, 1959,  DOI: 10.1103/PhysRevB.15.1959
    33. 33
      Li, C.-C.; Yoon, J.; Engel, L.; Shahar, D.; Tsui, D.; Shayegan, M. Microwave resonance and weak pinning in two-dimensional hole systems at high magnetic fields. Phys. Rev. B 2000, 61, 10905,  DOI: 10.1103/PhysRevB.61.10905
    34. 34
      Chen, Y.; Lewis, R.; Engel, L.; Tsui, D.; Ye, P.; Pfeiffer, L.; West, K. Microwave resonance of the 2D Wigner crystal around integer Landau fillings. Phys. Rev. Lett. 2003, 91, 016801,  DOI: 10.1103/PhysRevLett.91.016801
    35. 35
      Chitra, R.; Giamarchi, T.; Le Doussal, P. Pinned Wigner crystals. Phys. Rev. B 2001, 65, 035312,  DOI: 10.1103/PhysRevB.65.035312
    36. 36
      Pöllmann, C.; Steinleitner, P.; Leierseder, U.; Nagler, P.; Plechinger, G.; Porer, M.; Bratschitsch, R.; Schüller, C.; Korn, T.; Huber, R. Resonant internal quantum transitions and femtosecond radiative decay of excitons in monolayer WSe2. Nat. Mater. 2015, 14, 889893,  DOI: 10.1038/nmat4356

    Cited By

    Click to copy section linkSection link copied!

    This article is cited by 11 publications.

    1. Junho Choi, Jacob Embley, Daria D. Blach, Raül Perea-Causín, Daniel Erkensten, Dong Seob Kim, Long Yuan, Woo Young Yoon, Takashi Taniguchi, Kenji Watanabe, Keiji Ueno, Emanuel Tutuc, Samuel Brem, Ermin Malic, Xiaoqin Li, Libai Huang. Fermi Pressure and Coulomb Repulsion Driven Rapid Hot Plasma Expansion in a van der Waals Heterostructure. Nano Letters 2023, 23 (10) , 4399-4405. https://doi.org/10.1021/acs.nanolett.3c00678
    2. Daniel Erkensten, Samuel Brem, Raul Perea-Causin, Ermin Malic. Stability of Wigner crystals and Mott insulators in twisted moiré structures. Physical Review B 2024, 110 (15) https://doi.org/10.1103/PhysRevB.110.155132
    3. Roman Ya. Kezerashvili, Shalva M. Tsiklauri, Anastasia Spiridonova. Controllable trions in buckled two-dimensional materials. Physical Review B 2024, 110 (3) https://doi.org/10.1103/PhysRevB.110.035425
    4. Pavel E. Dolgirev, Ilya Esterlis, Alexander A. Zibrov, Mikhail D. Lukin, Thierry Giamarchi, Eugene Demler. Local Noise Spectroscopy of Wigner Crystals in Two-Dimensional Materials. Physical Review Letters 2024, 132 (24) https://doi.org/10.1103/PhysRevLett.132.246504
    5. Willy Knorr, Samuel Brem, Giuseppe Meneghini, Ermin Malic. Polaron-induced changes in moiré exciton propagation in twisted van der Waals heterostructures. Nanoscale 2024, 16 (18) , 8996-9003. https://doi.org/10.1039/D4NR00136B
    6. Hongyi Yu, Jiyong Zhou. Phonons of electronic crystals in two‐dimensional semiconductor moiré patterns. Natural Sciences 2023, 3 (3) https://doi.org/10.1002/ntls.20220065
    7. Huifang Ma, Yanzhao Yang, Heng Jing, Wanshun Jiang, Wenyue Guo, Hao Ren. Semi-Empirical model to retrieve finite temperature terahertz absorption spectra using Morse potential. Chinese Journal of Chemical Physics 2023, 36 (1) , 15-24. https://doi.org/10.1063/1674-0068/cjcp2202032
    8. Xi Feng, Xieyu Chen, Yongchang Lu, Qingwei Wang, Li Niu, Quan Xu, Xueqian Zhang, Jiaguang Han, Weili Zhang. Direct Emission of Focused Terahertz Vortex Beams Using Indium‐Tin‐Oxide‐Based Fresnel Zone Plates. Advanced Optical Materials 2023, 11 (1) https://doi.org/10.1002/adom.202201628
    9. Willy Knorr, Samuel Brem, Giuseppe Meneghini, Ermin Malic. Exciton transport in a moiré potential: From hopping to dispersive regime. Physical Review Materials 2022, 6 (12) https://doi.org/10.1103/PhysRevMaterials.6.124002
    10. Raul Perea-Causin, Daniel Erkensten, Jamie M. Fitzgerald, Joshua J. P. Thompson, Roberto Rosati, Samuel Brem, Ermin Malic. Exciton optics, dynamics, and transport in atomically thin semiconductors. APL Materials 2022, 10 (10) https://doi.org/10.1063/5.0107665
    11. Daniel Erkensten, Samuel Brem, Raül Perea-Causín, Ermin Malic. Microscopic origin of anomalous interlayer exciton transport in van der Waals heterostructures. Physical Review Materials 2022, 6 (9) https://doi.org/10.1103/PhysRevMaterials.6.094006

    Nano Letters

    Cite this: Nano Lett. 2022, 22, 3, 1311–1315
    Click to copy citationCitation copied!
    https://doi.org/10.1021/acs.nanolett.1c04620
    Published January 20, 2022

    Copyright © 2022 The Authors. Published by American Chemical Society. This publication is licensed under

    CC-BY 4.0 .

    Article Views

    3268

    Altmetric

    -

    Citations

    Learn about these metrics

    Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

    Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

    The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

    • Abstract

      Figure 1

      Figure 1. Sketch of the 2D Wigner crystal with a honeycomb lattice and alternating spin polarization. The colored curves underneath the particles illustrate their wave functions. A low-energy light field can probe the internal transition of the Wigner crystal toward excited quantum states to provide a conclusive proof of the phase transition.

      Figure 2

      Figure 2. Internal quantum structure of the WC in p-doped, hBN-encapsulated WSe2 at rs = 50. (a) Band structure of the WC, where only the lowest band is fully occupied. The background color illustrates the oscillator strength for optical transitions from the ground state. (b) The surface plot shows the charge density, whereas the projection contains the spin polarization ( is red, is blue). Also shown are the wave function at k = 0 of the (c) first excited state and (d) the state with the largest oscillator.

      Figure 3

      Figure 3. Terahertz fingerprint of the Wigner crystal. The real part of the optical conductivity exhibits a characteristic series of internal quantum transitions, whose amplitude sequence can be explained by the nature of the involved excited states. Note that the strongly broadened response signal has been scaled by a factor of 2, whereas the Drude signal of free electrons (black line) has been multiplied by 100.

      Figure 4

      Figure 4. Central energy of the optical response signal as a function of the hole density for different monolayers/substrates. The inset shows the result for an ideal 2D system, where the only system parameter is given by rs and the energy scale is given in units of effective Hartree.

    • References


      This article references 36 other publications.

      1. 1
        Wigner, E. On the interaction of electrons in metals. Phys. Rev. 1934, 46, 1002,  DOI: 10.1103/PhysRev.46.1002
      2. 2
        Andrei, E.; Deville, G.; Glattli, D.; Williams, F.; Paris, E.; Etienne, B. Observation of a magnetically induced Wigner solid. Phys. Rev. Lett. 1988, 60, 2765,  DOI: 10.1103/PhysRevLett.60.2765
      3. 3
        Goldman, V.; Santos, M.; Shayegan, M.; Cunningham, J. Evidence for two-dimentional quantum Wigner crystal. Phys. Rev. Lett. 1990, 65, 2189,  DOI: 10.1103/PhysRevLett.65.2189
      4. 4
        Williams, F.; Wright, P.; Clark, R.; Andrei, E.; Deville, G.; Glattli, D.; Probst, O.; Etienne, B.; Dorin, C.; Foxon, C. Conduction threshold and pinning frequency of magnetically induced Wigner solid. Phys. Rev. Lett. 1991, 66, 3285,  DOI: 10.1103/PhysRevLett.66.3285
      5. 5
        Wang, G.; Chernikov, A.; Glazov, M. M.; Heinz, T. F.; Marie, X.; Amand, T.; Urbaszek, B. Colloquium: Excitons in atomically thin transition metal dichalcogenides. Rev. Mod. Phys. 2018, 90, 021001,  DOI: 10.1103/RevModPhys.90.021001
      6. 6
        Mueller, T.; Malic, E. Exciton physics and device application of two-dimensional transition metal dichalcogenide semiconductors. npj 2D Mater. Appl. 2018, 2, 112,  DOI: 10.1038/s41699-018-0074-2
      7. 7
        Raja, A.; Waldecker, L.; Zipfel, J.; Cho, Y.; Brem, S.; Ziegler, J. D.; Kulig, M.; Taniguchi, T.; Watanabe, K.; Malic, E. Dielectric disorder in two-dimensional materials. Nat. Nanotechnol. 2019, 14, 832837,  DOI: 10.1038/s41565-019-0520-0
      8. 8
        Xu, Y.; Liu, S.; Rhodes, D. A.; Watanabe, K.; Taniguchi, T.; Hone, J.; Elser, V.; Mak, K. F.; Shan, J. Correlated insulating states at fractional fillings of moiré superlattices. Nature 2020, 587, 214218,  DOI: 10.1038/s41586-020-2868-6
      9. 9
        Wang, L.; Shih, E.-M.; Ghiotto, A.; Xian, L.; Rhodes, D. A.; Tan, C.; Claassen, M.; Kennes, D. M.; Bai, Y.; Kim, B. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 2020, 19, 861866,  DOI: 10.1038/s41563-020-0708-6
      10. 10
        Zhou, Y.; Sung, J.; Brutschea, E.; Esterlis, I.; Wang, Y.; Scuri, G.; Gelly, R. J.; Heo, H.; Taniguchi, T.; Watanabe, K. Bilayer Wigner crystals in a transition metal dichalcogenide heterostructure. Nature 2021, 595, 4852,  DOI: 10.1038/s41586-021-03560-w
      11. 11
        Huang, X.; Wang, T.; Miao, S.; Wang, C.; Li, Z.; Lian, Z.; Taniguchi, T.; Watanabe, K.; Okamoto, S.; Xiao, D. Correlated insulating states at fractional fillings of the WS2/WSe2 moiré lattice. Nat. Phys. 2021, 17, 715719,  DOI: 10.1038/s41567-021-01171-w
      12. 12
        Miao, S.; Wang, T.; Huang, X.; Chen, D.; Lian, Z.; Wang, C.; Blei, M.; Taniguchi, T.; Watanabe, K.; Tongay, S. Strong interaction between interlayer excitons and correlated electrons in WSe2/WS2 moiré superlattice. Nat. Commun. 2021, 12, 16,  DOI: 10.1038/s41467-021-23732-6
      13. 13
        Yu, H.; Liu, G.-B.; Tang, J.; Xu, X.; Yao, W. Moiré excitons: From programmable quantum emitter arrays to spin-orbit–coupled artificial lattices. Sci. Adv. 2017, 3, e1701696  DOI: 10.1126/sciadv.1701696
      14. 14
        Tran, K.; Moody, G.; Wu, F.; Lu, X.; Choi, J.; Kim, K.; Rai, A.; Sanchez, D. A.; Quan, J.; Singh, A. Evidence for moiré excitons in van der Waals heterostructures. Nature 2019, 567, 7175,  DOI: 10.1038/s41586-019-0975-z
      15. 15
        Brem, S.; Linderälv, C.; Erhart, P.; Malic, E. Tunable phases of moiré excitons in van der Waals heterostructures. Nano Lett. 2020, 20, 85348540,  DOI: 10.1021/acs.nanolett.0c03019
      16. 16
        Ugeda, M. M.; Bradley, A. J.; Shi, S.-F.; da Jornada, F. H.; Zhang, Y.; Qiu, D. Y.; Ruan, W.; Mo, S.-K.; Hussain, Z.; Shen, Z.-X. Giant bandgap renormalization and excitonic effects in a monolayer transition metal dichalcogenide semiconductor. Nat. Mater. 2014, 13, 10911095,  DOI: 10.1038/nmat4061
      17. 17
        Chernikov, A.; Berkelbach, T. C.; Hill, H. M.; Rigosi, A.; Li, Y.; Aslan, O. B.; Reichman, D. R.; Hybertsen, M. S.; Heinz, T. F. Exciton binding energy and nonhydrogenic Rydberg series in monolayer WS2. Phys. Rev. Lett. 2014, 113, 076802,  DOI: 10.1103/PhysRevLett.113.076802
      18. 18
        Smolenski, T.; Dolgirev, P. E.; Kuhlenkamp, C.; Popert, A.; Shimazaki, Y.; Back, P.; Lu, X.; Kroner, M.; Watanabe, K.; Taniguchi, T. Signatures of Wigner crystal of electrons in a monolayer semiconductor. Nature 2021, 595, 5357,  DOI: 10.1038/s41586-021-03590-4
      19. 19
        Drummond, N.; Radnai, Z.; Trail, J.; Towler, M.; Needs, R. Diffusion quantum Monte Carlo study of three-dimensional Wigner crystals. Phys. Rev. B 2004, 69, 085116,  DOI: 10.1103/PhysRevB.69.085116
      20. 20
        Trail, J.; Towler, M.; Needs, R. Unrestricted Hartree-Fock theory of Wigner crystals. Phys. Rev. B 2003, 68, 045107,  DOI: 10.1103/PhysRevB.68.045107
      21. 21
        Kormányos, A.; Burkard, G.; Gmitra, M.; Fabian, J.; Zólyomi, V.; Drummond, N. D.; Fal’ko, V. k-p theory for two-dimensional transition metal dichalcogenide semiconductors. 2D Materials 2015, 2, 022001,  DOI: 10.1088/2053-1583/2/2/022001
      22. 22
        Pan, H.; Wu, F.; Das Sarma, S. Quantum phase diagram of a Moiré-Hubbard model. Phys. Rev. B 2020, 102, 201104,  DOI: 10.1103/PhysRevB.102.201104
      23. 23
        Rytova, N. Screened potential of a point charge in a thin film. Moscow Univ. Phys. Bull. 1967, 3, 30
      24. 24
        Brem, S.; Ekman, A.; Christiansen, D.; Katsch, F.; Selig, M.; Robert, C.; Marie, X.; Urbaszek, B.; Knorr, A.; Malic, E. Phonon-assisted photoluminescence from indirect excitons in monolayers of transition-metal dichalcogenides. Nano Lett. 2020, 20, 28492856,  DOI: 10.1021/acs.nanolett.0c00633
      25. 25
        Lindberg, M.; Koch, S. W. Effective Bloch equations for semiconductors. Phys. Rev. B 1988, 38, 3342,  DOI: 10.1103/PhysRevB.38.3342
      26. 26
        Kira, M.; Koch, S. Many-body correlations and excitonic effects in semiconductor spectroscopy. Progress Quantum Electron. 2006, 30, 155296,  DOI: 10.1016/j.pquantelec.2006.12.002
      27. 27
        Laturia, A.; Van de Put, M. L.; Vandenberghe, W. G. Dielectric properties of hexagonal boron nitride and transition metal dichalcogenides: from monolayer to bulk. npj 2D Mater. Appl. 2018, 2, 17,  DOI: 10.1038/s41699-018-0050-x
      28. 28
        Yoon, J.; Li, C.; Shahar, D.; Tsui, D.; Shayegan, M. Wigner crystallization and metal-insulator transition of two-dimensional holes in GaAs at B = 0. Phys. Rev. Lett. 1999, 82, 1744,  DOI: 10.1103/PhysRevLett.82.1744
      29. 29
        Drummond, N.; Needs, R. Phase diagram of the low-density two-dimensional homogeneous electron gas. Phys. Rev. Lett. 2009, 102, 126402,  DOI: 10.1103/PhysRevLett.102.126402
      30. 30
        Zarenia, M.; Neilson, D.; Partoens, B.; Peeters, F. Wigner crystallization in transition metal dichalcogenides: A new approach to correlation energy. Phys. Rev. B 2017, 95, 115438,  DOI: 10.1103/PhysRevB.95.115438
      31. 31
        Tanatar, B.; Ceperley, D. M. Ground state of the two-dimensional electron gas. Phys. Rev. B 1989, 39, 5005,  DOI: 10.1103/PhysRevB.39.5005
      32. 32
        Bonsall, L.; Maradudin, A. Some static and dynamical properties of a two-dimensional Wigner crystal. Phys. Rev. B 1977, 15, 1959,  DOI: 10.1103/PhysRevB.15.1959
      33. 33
        Li, C.-C.; Yoon, J.; Engel, L.; Shahar, D.; Tsui, D.; Shayegan, M. Microwave resonance and weak pinning in two-dimensional hole systems at high magnetic fields. Phys. Rev. B 2000, 61, 10905,  DOI: 10.1103/PhysRevB.61.10905
      34. 34
        Chen, Y.; Lewis, R.; Engel, L.; Tsui, D.; Ye, P.; Pfeiffer, L.; West, K. Microwave resonance of the 2D Wigner crystal around integer Landau fillings. Phys. Rev. Lett. 2003, 91, 016801,  DOI: 10.1103/PhysRevLett.91.016801
      35. 35
        Chitra, R.; Giamarchi, T.; Le Doussal, P. Pinned Wigner crystals. Phys. Rev. B 2001, 65, 035312,  DOI: 10.1103/PhysRevB.65.035312
      36. 36
        Pöllmann, C.; Steinleitner, P.; Leierseder, U.; Nagler, P.; Plechinger, G.; Porer, M.; Bratschitsch, R.; Schüller, C.; Korn, T.; Huber, R. Resonant internal quantum transitions and femtosecond radiative decay of excitons in monolayer WSe2. Nat. Mater. 2015, 14, 889893,  DOI: 10.1038/nmat4356
    • Supporting Information

      Supporting Information


      The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.nanolett.1c04620.

      • Transformation of Hartree–Fock Hamiltonian into Wigner basis; material-specific parameters and interaction matrix elements; derivation of the optical response (PDF)


      Terms & Conditions

      Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.