Generalized Wigner Crystals in Moiré TMDs

• Generalized Wigner Crystals (GWCs) are interaction-driven charge-ordered insulators formed at commensurate fractional fillings, where longer-range Coulomb repulsion overcomes kinetic energy.
• They exhibit diverse charge patterns—such as √3×√3, stripe, and honeycomb orders—whose selection depends on lattice geometry, filling fraction, and interaction range.
• Advanced models and techniques (e.g., extended Hubbard models, exact diagonalization, and STM imaging) reveal complex quantum melting, magnetic ordering, and depinning dynamics in GWCs.

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 (Zhou et al., 2023, Padhi et al., 2020).

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 (Zhou et al., 2023, Padhi et al., 2020).

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 (Padhi et al., 2020).

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 (Jaworowski et al., 2017). 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 (Shapir et al., 2018).

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 (Matty et al., 2021, Morales-Durán et al., 2022, Li et al., 2022).

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,

$n_{\rm Mott}^{1/2} a_\ast = 1$0

In moiré TMD applications, the $n_{\rm Mott}^{1/2} a_\ast = 1$1 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 (Zhou et al., 2023, Kumar et al., 21 Apr 2026).

Projected continuum models are also widely used. In that approach, one starts from a moiré single-particle Hamiltonian such as

$n_{\rm Mott}^{1/2} a_\ast = 1$2

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 (Morales-Durán et al., 2022, Ung et al., 2023, Reichhardt et al., 2024).

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_{\rm Mott}^{1/2} a_\ast = 1$3 and $n_{\rm Mott}^{1/2} a_\ast = 1$4, the full long-range interaction robustly selects the $n_{\rm Mott}^{1/2} a_\ast = 1$5 crystal for realistic gate separations, whereas truncations to $n_{\rm Mott}^{1/2} a_\ast = 1$6, $n_{\rm Mott}^{1/2} a_\ast = 1$7, or $n_{\rm Mott}^{1/2} a_\ast = 1$8 can reshuffle the energetic hierarchy and stabilize competing dimer, stripe, trimer, or four-mer-like states (Kumar et al., 2024).

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
$n_{\rm Mott}^{1/2} a_\ast = 1$9	$_2$0 triangular crystal	finite-temperature Lanczos study (Kumar et al., 21 Apr 2026)
$_2$1	$_2$2 triangular crystal	STM, DMRG, ED (Li et al., 2021, Zhou et al., 2023)
$_2$3	columnar dimer crystal	classical melting study (Matty et al., 2021)
$_2$4	stripe crystal	STM and transport studies (Li et al., 2021, Zong et al., 22 May 2025)
$_2$5	elongated honeycomb lattice	Hartree–Fock interpretation of transport (Zong et al., 22 May 2025)
$_2$6	honeycomb, or $_2$7, pattern	STM and DMRG (Li et al., 2021, Biborski et al., 2024)

At $_2$8, the standard state is the $_2$9 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

$_2$0

At $_2$1, the complementary filling pattern occupies two of the three sublattices and forms a honeycomb network; at $_2$2, the favored state breaks $_2$3 symmetry and forms stripes (Zhou et al., 2023, Li et al., 2021).

These apparently simple commensurate patterns sit inside a more intricate energy landscape. Hartree–Fock calculations for hole-doped WSe$_2$4/WS$_2$5 found many competing ordered and disordered states across $_2$6, separated by energies of a fraction of $_2$7 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 (Ung et al., 2023).

Charge frustration is central to this competition. In finite-range $_2$8 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$9 triangular charge pattern plus delocalized carriers, with coexisting solid-like and liquid-like features and no classical analog (Kumar et al., 2024). 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 $_2$0 triangular-lattice GWC, both DMRG and momentum-space exact diagonalization found antiferromagnetic $_2$1 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 $_2$2 charge order; an appendix reconstruction of the six-peak pattern in $_2$3 was fully consistent with the classic triangular-lattice $_2$4 antiferromagnet (Zhou et al., 2023, Morales-Durán et al., 2022).

At $_2$5, 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

$_2$6

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 $_2$7 to $_2$8 (Morales-Durán et al., 2022).

Spin-valley locking modifies this picture further. In a single-band extended-Hubbard description of WSe$_2$9/WS$_2$0 at $_2$1, increasing nearest-neighbor repulsion $_2$2 drives a metal-to-insulator transition concomitant with the emergence of the honeycomb $_2$3 GWC. Within that phase, the occupied honeycomb network exhibits out-of-plane antiferromagnetic order in $_2$4, but growing in-plane correlations $_2$5 produce spin canting; the spin structure factor develops peaks at $_2$6 and $_2$7 together with a growing $_2$8 peak, and this mixed character disappears when the spin-valley-locking phase $_2$9 is set to zero (Biborski et al., 2024).

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 $_2$0 Néel order. The calculated crossover occurs at $_2$1 for the triangular potential and $_2$2 for the honeycomb potential, implying that triangular commensurate GWCs are far more susceptible to antiferromagnetic locking than honeycomb ones (Esterlis et al., 1 Feb 2025).

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 $_2$3, a single direct first-order transition was found between a Fermi liquid and the $_2$4 GWC, with no evidence for an intermediate phase in the accessible parameter regime. Along the cut $_2$5, the transition occurs around $_2$6: the quasiparticle residue $_2$7 extracted from the jump in $_2$8 drops abruptly as the charge-order structure factor $_2$9 simultaneously develops peaks at the $_2$0 points (Zhou et al., 2023).

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 $_2$1 (Morales-Durán et al., 2022). Electrical transport in twisted bilayer MoSe$_2$2 reported continuous quantum melting transitions of GWCs at $_2$3 and $_2$4, driven by doping density, magnetic field, and displacement field, with quantum critical scaling

$_2$5

and exponents $_2$6 at $_2$7 and $_2$8 at $_2$9 (Zong et al., 22 May 2025). 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 $D_3$0 GWC first gives way, upon slight doping, to an isotropic compressible hexagonal domain state formed by domain walls between the three registries of the $D_3$1 crystal; further density increase produces a type-II nematic associated with fragmented $D_3$2 columnar dimer order, while dilution from $D_3$3 produces a type-I nematic as a short-ranged version of the stripe crystal. The nematic order parameter transforms as the two-dimensional $D_3$4 representation of the moiré point group, with type-I and type-II distinguished by $D_3$5 and $D_3$6, respectively (Matty et al., 2021).

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 $D_3$7 at some fillings, and that the sign of the quantum shift depends on filling: at $D_3$8, increasing $D_3$9 lowers $$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,$$0, but at $$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,$$1 and $$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,$$2, $$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,$$3 increases with $$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,$$4 by as much as $$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,$$5 on the finite clusters studied. In the perturbative regime, $$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,$$6, with the sign determined by whether kinetic processes lower domain-wall energies or stabilize the ordered free energy by reducing the entropy jump (Kumar et al., 21 Apr 2026).

Driven dynamics add another layer. Overdamped simulations of long-range repulsive charges on a hexagonal substrate found strongly pinned commensurate states at $$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,$$7, $$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,$$8, and $$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,$$9. At $n_{\rm Mott}^{1/2} a_\ast = 1$00 and $n_{\rm Mott}^{1/2} a_\ast = 1$01, $n_{\rm Mott}^{1/2} a_\ast = 1$02; at $n_{\rm Mott}^{1/2} a_\ast = 1$03, hole doping produces a two-step depinning in which anti-kinks first slide along the stripes before the whole crystal depins, while at $n_{\rm Mott}^{1/2} a_\ast = 1$04 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 (Reichhardt et al., 2024). STM imaging of density-driven melting in near-$n_{\rm Mott}^{1/2} a_\ast = 1$05 twisted MoSe$n_{\rm Mott}^{1/2} a_\ast = 1$06 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 (Berger et al., 18 Dec 2025).

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$n_{\rm Mott}^{1/2} a_\ast = 1$07/WS$n_{\rm Mott}^{1/2} a_\ast = 1$08, separated by an hBN spacer of about $n_{\rm Mott}^{1/2} a_\ast = 1$09, directly visualized different charge configurations at $n_{\rm Mott}^{1/2} a_\ast = 1$10, $n_{\rm Mott}^{1/2} a_\ast = 1$11, and $n_{\rm Mott}^{1/2} a_\ast = 1$12 on a moiré lattice with $n_{\rm Mott}^{1/2} a_\ast = 1$13. The $n_{\rm Mott}^{1/2} a_\ast = 1$14 and $n_{\rm Mott}^{1/2} a_\ast = 1$15 states were observed as triangular and honeycomb lattices, respectively, while the $n_{\rm Mott}^{1/2} a_\ast = 1$16 state formed stripes with spontaneously broken $n_{\rm Mott}^{1/2} a_\ast = 1$17 symmetry. Bias-dependent expanding rings identified tip-induced discharge as the dominant imaging mechanism rather than simple topographic contrast (Li et al., 2021).

Scanning single-electron charging (SSEC) spectroscopy then extended this program from static charge order to charge excitations. In twisted WS$n_{\rm Mott}^{1/2} a_\ast = 1$18 moiré heterostructures, a graphene sensing layer enabled single-electron charge resolution and $n_{\rm Mott}^{1/2} a_\ast = 1$19 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

$n_{\rm Mott}^{1/2} a_\ast = 1$20

using a tip-gating efficiency $n_{\rm Mott}^{1/2} a_\ast = 1$21 (Li et al., 2022).

Transport has become comparably informative. In twisted bilayer MoSe$n_{\rm Mott}^{1/2} a_\ast = 1$22, $n_{\rm Mott}^{1/2} a_\ast = 1$23 peaks, activation gaps, and nonlinear $n_{\rm Mott}^{1/2} a_\ast = 1$24 depinning features identified GWCs at $n_{\rm Mott}^{1/2} a_\ast = 1$25 and $n_{\rm Mott}^{1/2} a_\ast = 1$26. In zero magnetic field, the maximum activation gaps were reported as $n_{\rm Mott}^{1/2} a_\ast = 1$27 K at $n_{\rm Mott}^{1/2} a_\ast = 1$28, $n_{\rm Mott}^{1/2} a_\ast = 1$29 K at $n_{\rm Mott}^{1/2} a_\ast = 1$30, and $n_{\rm Mott}^{1/2} a_\ast = 1$31 K at $n_{\rm Mott}^{1/2} a_\ast = 1$32, while the collapse of the low-bias $n_{\rm Mott}^{1/2} a_\ast = 1$33 peak at large bias was interpreted as depinning or sliding of the crystal (Zong et al., 22 May 2025).

Excited-state correlations now extend the subject beyond static insulating order. First-principles GW-Bethe–Salpeter calculations for angle-aligned MoSe$n_{\rm Mott}^{1/2} a_\ast = 1$34/MoS$n_{\rm Mott}^{1/2} a_\ast = 1$35 at hole fillings $n_{\rm Mott}^{1/2} a_\ast = 1$36 and $n_{\rm Mott}^{1/2} a_\ast = 1$37 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 $n_{\rm Mott}^{1/2} a_\ast = 1$38 meV at $n_{\rm Mott}^{1/2} a_\ast = 1$39 and $n_{\rm Mott}^{1/2} a_\ast = 1$40 meV at $n_{\rm Mott}^{1/2} a_\ast = 1$41, compared with a kinetic scale of only about $n_{\rm Mott}^{1/2} a_\ast = 1$42 meV, and the lowest-energy WCEs were found to be optically dark because the valence and conduction states reside in different valleys (You et al., 10 Sep 2025).

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 (Kumar et al., 2024). Magnetic crossover scales are predicted to lie in the range $n_{\rm Mott}^{1/2} a_\ast = 1$43, with particle–hole differences of a few hundred mK, motivating spin-sensitive probes such as NV-center magnetometry, spin-polarized STM, or nanoSQUID measurements (Kumar et al., 2024). 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 (Matty et al., 2021, Reichhardt et al., 2024, You et al., 10 Sep 2025).