---
title: Generalized Wigner Crystals in Moiré TMDs
url: https://www.emergentmind.com/topics/generalized-wigner-crystals-gwcs
type: topic
---

# Generalized Wigner Crystals in Moiré TMDs

Generalized Wigner crystals (GWCs) are interaction-driven charge-ordered insulating states that form on a lattice at commensurate fractional filling when longer-range Coulomb repulsion dominates over kinetic delocalization. They are “Wigner-like” because further-range interactions drive charge localization, and “Mott-like” because lattice commensurability is essential to the ordering pattern. In moiré transition-metal dichalcogenide (TMD) systems, GWCs have become a central framework for interpreting incompressible states at fractional fillings, especially on triangular superlattices supporting \(\sqrt{3}\times\sqrt{3}\), stripe, and honeycomb charge order [2312.03828], [2009.13536].

## 1. Definition and conceptual scope

A conventional Wigner crystal is the low-density electron solid of a continuum system, where long-range Coulomb repulsion overwhelms kinetic energy and spontaneously selects a crystalline arrangement. A Mott insulator, by contrast, is a lattice insulator driven primarily by on-site repulsion, typically at integer filling. A GWC lies between these limits: the lattice is present from the outset, but the decisive interaction physics is not purely local. Fractional filling, lattice commensurability, and extended repulsion cooperate to produce a pinned, incompressible charge crystal rather than a uniform correlated metal or a purely on-site Mott state [2312.03828], [2009.13536].

In moiré materials this distinction is sharpened by the periodic moiré potential itself. The lattice not only constrains the allowed charge patterns but also pins them, so the insulating state is incompressible through a pinning gap even without disorder. In one survey of twisted TMD homobilayers and heterobilayers, the densities required to satisfy the Mott criterion \(n_{\rm Mott}^{1/2} a_\ast = 1\) were found to be nearly three orders of magnitude larger than the maximal density for GWC formation, indicating that typical operating densities are far from the Mott regime [2009.13536].

Although the modern literature is dominated by moiré TMDs, the notion is broader. Exact diagonalization in partially filled topological flat bands found low-density lattice Wigner crystals whose shape is largely controlled by boundary conditions and classical point-particle geometry rather than by band topology, which is a lattice-generalized form of Wigner crystallization rather than a purely continuum one [1712.06007]. In a broader generalized sense, finite one-dimensional few-electron nanotube crystals also realize interaction-dominated charge ordering with suppressed exchange and a spin-incoherent regime [1803.08523].

## 2. Moiré platforms and microscopic descriptions

The canonical GWC platforms are triangular moiré superlattices in TMD heterobilayers and twisted bilayers, including WSe\(_2\)/WS\(_2\), MoTe\(_2\)/WSe\(_2\), twisted WS\(_2\), and twisted MoSe\(_2\). In zero-twist or near-aligned heterobilayers, electrons or holes occupy localized Wannier orbitals on an emergent triangular lattice with \(D_3\) point-group symmetry; in twisted systems, narrow minibands strongly suppress kinetic energy, making Coulomb interactions dominant at special fractional fillings [2112.08624], [2210.15168], [2209.12830].

A common microscopic description is the extended Hubbard model on the triangular lattice,
\[
H = - \sum_{ij\sigma} t_{ij}\, c^\dagger_{i\sigma} c_{j\sigma}
+ U \sum_i n_{i\uparrow}n_{i\downarrow}
+ \sum_{i<j} V_{ij} n_i n_j,
\]
or, in nearest-neighbor form,
\[
H =  -t\sum_{\langle ij\rangle}\left(c_i^\dag c_j + H.c.\right)
+ U\sum_{i} n_{i,\uparrow} n_{i,\downarrow}
+ V_1 \sum_{\langle ij\rangle} n_i n_j.
\]
In moiré TMD applications, the \(\uparrow,\downarrow\) labels often represent valley flavors rather than literal spin, and the physically relevant interaction is frequently a dual-gate screened long-range Coulomb potential obtained by the method of images rather than a strictly truncated short-range form [2312.03828], [2604.19870].

Projected continuum models are also widely used. In that approach, one starts from a moiré single-particle Hamiltonian such as
\[
H_0=-\frac{\hbar^2}{2m^*}{\bm k}^2+\Delta({\bm r}),
\qquad
\Delta({\bm r})=2V_m\sum_{j=1,3,5}\cos({\bm b}_j\cdot{\bm r}+\psi),
\]
and projects Coulomb interactions into the topmost moiré band. Momentum-space exact diagonalization, Hartree–Fock, DMRG, finite-temperature Lanczos calculations, classical Monte Carlo, and overdamped classical dynamics have all been deployed to study complementary aspects of GWC formation, magnetism, melting, and transport [2210.15168], [2308.03020], [2406.15663].

A recurrent theoretical issue is the treatment of interaction range. Long-range Coulomb tails materially affect both the selected charge pattern and the relevant energy scales. In triangular moiré systems at \(n=1/3\) and \(2/3\), the full long-range interaction robustly selects the \(\sqrt{3}\times\sqrt{3}\) crystal for realistic gate separations, whereas truncations to \(V_1\), \(V_2\), or \(V_3\) can reshuffle the energetic hierarchy and stabilize competing dimer, stripe, trimer, or four-mer-like states [2409.13814].

## 3. Commensurate charge order and competing crystalline patterns

The most prominent GWC patterns in triangular moiré lattices recur across theory and experiment.

| Filling | Charge pattern | Representative context |
|---|---|---|
| \(1/4\) | \(2\times2\) triangular crystal | finite-temperature Lanczos study [2604.19870] |
| \(1/3\) | \(\sqrt{3}\times\sqrt{3}\) triangular crystal | STM, DMRG, ED [2106.10599], [2312.03828] |
| \(2/5\) | columnar dimer crystal | classical melting study [2112.08624] |
| \(1/2\) | stripe crystal | STM and transport studies [2106.10599], [2505.16681] |
| \(3/5\) | elongated honeycomb lattice | Hartree–Fock interpretation of transport [2505.16681] |
| \(2/3\) | honeycomb, or \((110)\), pattern | STM and DMRG [2106.10599], [2409.11202] |

At \(\nu=1/3\), the standard state is the \(\sqrt{3}\times\sqrt{3}\) GWC, in which one of the three triangular sublattices is preferentially occupied and the unit cell is enlarged by a factor of three. In momentum space the corresponding ordering wavevectors were given as
\[
\mathbf{k}_1=\left(\frac{2\pi}{3a}, \frac{2\pi}{\sqrt{3}a}\right),
\qquad
\mathbf{k}_2=\left(\frac{4\pi}{3a}, 0\right).
\]
At \(\nu=2/3\), the complementary filling pattern occupies two of the three sublattices and forms a honeycomb network; at \(\nu=1/2\), the favored state breaks \(C_3\) symmetry and forms stripes [2312.03828], [2106.10599].

These apparently simple commensurate patterns sit inside a more intricate energy landscape. Hartree–Fock calculations for hole-doped WSe\(_2\)/WS\(_2\) found many competing ordered and disordered states across \(\nu=\frac17,\frac14,\frac13,\frac25,\frac12,\frac35,\frac23,\frac34,\frac67,1\), separated by energies of a fraction of \(\sim 1\) meV/hole. Stability analysis of the Hartree–Fock Hessian showed that many of these are genuine metastable minima rather than numerical artifacts, while HCI, CCSD, and CCSD(T) found correlation energies that are small in absolute size but large enough to reorder nearly degenerate Hartree–Fock states [2308.03020].

Charge frustration is central to this competition. In finite-range \(V_1\text{-}V_2\text{-}V_3\) descriptions, the system can cluster particles to avoid short-range penalties even without attraction, generating dimer crystals, stripes, trimers, and related patterns. A particularly notable result is the “pinball” phase: a partially quantum melted GWC with a \(2\times2\) triangular charge pattern plus delocalized carriers, with coexisting solid-like and liquid-like features and no classical analog [2409.13814]. This suggests that the boundary between “crystal” and “metal” is not always binary in interaction-frustrated regimes.

## 4. Magnetic structure and internal degrees of freedom

Once the charge sector crystallizes, the residual spin or valley degrees of freedom generate nontrivial magnetic order. In the \(\nu=1/3\) triangular-lattice GWC, both DMRG and momentum-space exact diagonalization found antiferromagnetic \(120^\circ\) order on the occupied triangular sublattice. In DMRG on YC6 triangular cylinders, the GWC phase showed long-range antiferromagnetic order in the valley-spin channel, with spin structure factor peaks at the corners of the reduced Brillouin zone generated by the \(\sqrt{3}\times\sqrt{3}\) charge order; an appendix reconstruction of the six-peak pattern in \(S(\mathbf{k})\) was fully consistent with the classic triangular-lattice \(120^\circ\) antiferromagnet [2312.03828], [2210.15168].

At \(\nu=2/3\), where the occupied sites form a honeycomb network, antiferromagnetism again appears naturally, but with greater tunability. In projected continuum exact diagonalization, the effective nearest-neighbor exchange
\[
J_1(2/3)=\frac{4t_1^2}{U_0-U_1}-2X_1
\]
can change sign because direct exchange competes strongly with superexchange; reducing dielectric screening can therefore drive a crossover from antiferromagnetism to ferromagnetism, visible as a change in the dominant spin-structure-factor peaks from \({\bm\kappa}\) to \({\bm\gamma}\) [2210.15168].

Spin-valley locking modifies this picture further. In a single-band extended-Hubbard description of WSe\(_2\)/WS\(_2\) at \(n=2/3\), increasing nearest-neighbor repulsion \(V\) drives a metal-to-insulator transition concomitant with the emergence of the honeycomb \((110)\) GWC. Within that phase, the occupied honeycomb network exhibits out-of-plane antiferromagnetic order in \(\mathcal{S}^z_{ij}\), but growing in-plane correlations \(\mathcal{S}^{xy}_{ij}\) produce spin canting; the spin structure factor develops peaks at \(\mathbf{K}\) and \(\mathbf{K}'\) together with a growing \(\mathbf{\Gamma}\) peak, and this mixed character disappears when the spin-valley-locking phase \(\phi\) is set to zero [2409.11202].

A complementary semiclassical perspective comes from multiparticle ring exchange. For a two-dimensional triangular Wigner crystal in a commensurate triangular or honeycomb moiré potential, odd-particle exchanges favor ferromagnetism while even-particle exchanges favor antiferromagnetism. Increasing the moiré potential suppresses longer odd exchanges and promotes two-particle exchange, driving a transition to \(120^\circ\) Néel order. The calculated crossover occurs at \(\epsilon_c \approx 0.0033\) for the triangular potential and \(\epsilon_c \approx 0.26\) for the honeycomb potential, implying that triangular commensurate GWCs are far more susceptible to antiferromagnetic locking than honeycomb ones [2502.00323].

## 5. Quantum melting, thermal melting, depinning, and transport

One of the central questions in the field is how a GWC disappears as bandwidth, temperature, or carrier density is varied. In a DMRG study of the triangular-lattice extended Hubbard model at \(\nu=1/3\), a single direct first-order transition was found between a Fermi liquid and the \(\sqrt{3}\times\sqrt{3}\) GWC, with no evidence for an intermediate phase in the accessible parameter regime. Along the cut \(V_1/U=0.4\), the transition occurs around \(U/t \approx 11.8\): the quasiparticle residue \(Z\) extracted from the jump in \(n(\mathbf{k})\) drops abruptly as the charge-order structure factor \(N(\mathbf{k})\) simultaneously develops peaks at the \(K\) points [2312.03828].

Other approaches produce different melting phenomenology. Momentum-space exact diagonalization of projected continuum moiré models found a broader “melting regime” in which charge order survives while quantum fluctuations generate spin dynamics and excitonic particle-hole bound states; the charge gap then collapses with only a small jump compared with the atomic-limit gap, and the onset of strong quantum melting in finite-size calculations was placed roughly around \(W/V_m \sim 0.05\text{–}0.1\) [2210.15168]. Electrical transport in twisted bilayer MoSe\(_2\) reported continuous quantum melting transitions of GWCs at \(\nu = 2/5, 1/2, 3/5, 2/3, 8/9, 10/9,\) and \(4/3\), driven by doping density, magnetic field, and displacement field, with quantum critical scaling
\[
R_{xx}(T) = R_c(T)\, f\!\left(\frac{T}{T_0(\delta B)}\right),\qquad
T_0(\delta B) \sim |\delta B|^{z\nu},
\]
and exponents \(z\nu = 1.34 \pm 0.02\) at \(\nu=1/2\) and \(1.31 \pm 0.03\) at \(\nu=2/3\) [2505.16681]. This suggests that the apparent order of the melting transition is model- and tuning-dependent.

At finite temperature, the melting sequence can be even richer. Classical Monte Carlo for a strong-coupling triangular lattice gas found that the \(\nu=1/3\) GWC first gives way, upon slight doping, to an isotropic compressible hexagonal domain state formed by domain walls between the three registries of the \(1/3\) crystal; further density increase produces a type-II nematic associated with fragmented \(2/5\) columnar dimer order, while dilution from \(\nu=1/2\) produces a type-I nematic as a short-ranged version of the stripe crystal. The nematic order parameter transforms as the two-dimensional \(E\) representation of the moiré point group, with type-I and type-II distinguished by \(\langle\cos(6\theta)\rangle = +1\) and \(-1\), respectively [2112.08624].

Quantum corrections do not simply reduce transition temperatures. Finite-temperature Lanczos calculations on an extended triangular-lattice Hubbard model showed that classical estimates can miss experimental melting temperatures by as much as \(\sim 55\%\) at some fillings, and that the sign of the quantum shift depends on filling: at \(n=1/3\), increasing \(t\) lowers \(T_c\), but at \(n=1/2\) and \(n=1/4\), \(T_c\) increases with \(t\) by as much as \(30\!-\!40\%\) on the finite clusters studied. In the perturbative regime, \(\delta T_c \propto t^2\), with the sign determined by whether kinetic processes lower domain-wall energies or stabilize the ordered free energy by reducing the entropy jump [2604.19870].

Driven dynamics add another layer. Overdamped simulations of long-range repulsive charges on a hexagonal substrate found strongly pinned commensurate states at \(\nu=1/3\), \(1/2\), and \(2/3\). At \(\nu=1/3\) and \(1/2\), \(F_c = F_p\); at \(\nu=1/2\), hole doping produces a two-step depinning in which anti-kinks first slide along the stripes before the whole crystal depins, while at \(\nu=2/3\) the pinned honeycomb state depins discontinuously into a floating hexagonal lattice. Transport is strongly asymmetric on the two sides of commensuration because vacancies, interstitials, kinks, and anti-kinks have distinct mobilities [2406.15663]. STM imaging of density-driven melting in near-\(60^\circ\) twisted MoSe\(_2\) bilayers likewise revealed pronounced electron–hole asymmetry: hole-doped GWCs melt into interaction-driven disordered states with telegraph noise and short-range remnants of the parent order, whereas electron-doped GWCs melt into spatially uniform liquid-like compressible states because the reconstructed electron and hole pockets couple differently to the GWC gap [2512.16050].

## 6. Experimental probes, spectroscopy, and broader directions

The decisive experimental advance for GWCs in moiré TMDs was direct real-space imaging. A non-invasive STM technique using a monolayer graphene sensing layer above WSe\(_2\)/WS\(_2\), separated by an hBN spacer of about \(d = 5\,\text{nm}\), directly visualized different charge configurations at \(n=1/3\), \(1/2\), and \(2/3\) on a moiré lattice with \(L_M \approx 8\,\text{nm}\). The \(n=1/3\) and \(n=2/3\) states were observed as triangular and honeycomb lattices, respectively, while the \(n=1/2\) state formed stripes with spontaneously broken \(C_3\) symmetry. Bias-dependent expanding rings identified tip-induced discharge as the dominant imaging mechanism rather than simple topographic contrast [2106.10599].

Scanning single-electron charging (SSEC) spectroscopy then extended this program from static charge order to charge excitations. In twisted WS\(_2\) moiré heterostructures, a graphene sensing layer enabled single-electron charge resolution and \(\sim 1\) nm spatial resolution, allowing direct imaging of electron and hole quasiparticles across the GWC gap. Electron and hole wavefunctions were found to be complementary in real space, and thermodynamic gaps were extracted as
\[
\Delta_{n=1/3} = 52 \pm 8\ \text{meV},\qquad
\Delta_{n=2/3} = 47 \pm 8\ \text{meV},\qquad
\Delta_{n=1} = 107 \pm 16\ \text{meV},
\]
using a tip-gating efficiency \(\alpha = 0.16 \pm 0.02\) [2209.12830].

Transport has become comparably informative. In twisted bilayer MoSe\(_2\), \(R_{xx}\) peaks, activation gaps, and nonlinear \(dV/dI\) depinning features identified GWCs at \(\nu = 2/5, 1/2, 3/5, 2/3, 8/9, 10/9,\) and \(4/3\). In zero magnetic field, the maximum activation gaps were reported as \(2.2\) K at \(\nu=2/3\), \(4.6\) K at \(\nu=3/5\), and \(5.2\) K at \(\nu=1/2\), while the collapse of the low-bias \(dV/dI\) peak at large bias was interpreted as depinning or sliding of the crystal [2505.16681].

Excited-state correlations now extend the subject beyond static insulating order. First-principles GW-Bethe–Salpeter calculations for angle-aligned MoSe\(_2\)/MoS\(_2\) at hole fillings \(v_h=1/3\) and \(2/3\) found “Wigner-crystalline excitons” whose electron and hole densities both inherit the broken translational symmetry of the GWC ground state. The exciton binding energies were estimated as \(\sim 118\) meV at \(v_h=1/3\) and \(\sim 98\) meV at \(v_h=2/3\), compared with a kinetic scale of only about \(2\text{–}3\) meV, and the lowest-energy WCEs were found to be optically dark because the valence and conduction states reside in different valleys [2509.08211].

Several unresolved issues remain active. Long-range Coulomb tails, gate screening, and moiré geometry materially reshape the phase diagram, and effective short-range descriptions require careful renormalization rather than naive truncation [2409.13814]. Magnetic crossover scales are predicted to lie in the range \(T_m \sim 0.2\ \mathrm{K} - 1\ \mathrm{K}\), with particle–hole differences of a few hundred mK, motivating spin-sensitive probes such as NV-center magnetometry, spin-polarized STM, or nanoSQUID measurements [2409.13814]. More generally, the present body of work suggests that GWCs are best viewed not as isolated fractional insulators but as the organizing centers of a broader hierarchy of correlated phases—nematic fluids, domain-wall metals, pinball states, sliding crystals, and excitonic states—whose accessibility is controlled by commensurability, interaction range, screening, and bandwidth tuning [2112.08624], [2406.15663], [2509.08211].

Source: https://www.emergentmind.com/topics/generalized-wigner-crystals-gwcs