Correlated Insulators in Moiré Systems
A correlated insulator in a moiré system is a charge-incompressible or charge-localized phase at a commensurate moiré filling whose insulating behavior requires electron interactions beyond a calibrated independent-particle miniband description. This definition deliberately separates an observation from its interpretation. A resistance peak at an integer or fractional filling is evidence for inhibited transport; it is not yet proof of a thermodynamic gap, a Mott state, a Wigner crystal, or any particular order.
At one carrier per moiré cell, repulsion can suppress local charge fluctuations while leaving spin or valley degrees of freedom active. At one-third filling, longer-range repulsion can instead select a repeating pattern of occupied cells. These examples motivate the central question: does insulation arise from local charge constraints, spatial order, a reconstructed band structure, or a combination?
Required background. Moiré Superlattices supplies the unit-cell and filling scales. Flat Bands supplies the distinction between suppressed bandwidth and interaction-driven order. Mott Insulators supplies charge-incompressibility criteria and the tests that distinguish a Mott mechanism from other causes of insulating transport.
Interaction-driven insulation in moiré bands
Section titled “Interaction-driven insulation in moiré bands”We will first translate gate density into carriers per cell, then compare local constraints with enlarged-period charge order. Thermodynamic, transport and imaging measurements answer different parts of this question. The general charge-gap definition is developed in Mott Insulators; here we use it to interpret moiré devices.
Filling-Controlled Phases
Section titled “Filling-Controlled Phases”Filling is a calibrated, signed quantity
Section titled “Filling is a calibrated, signed quantity”For a device with local carrier density , reference density , and moiré-cell area , define
The sign, reference, and active-manifold capacity must accompany the number. In one convention is charge neutrality; in another it is an empty conduction miniband or a filled valence miniband. A four-flavor manifold may be called full at , whereas a normalized convention calls the same point filling one. “Half filling” is therefore ambiguous unless the number of active one-particle states per cell is stated.
For a triangular moiré lattice of period , the cell area is .
The inferred density of one carrier per cell is only as accurate as the local twist, strain, reconstruction, and electrostatic calibration. A small spatial variation of produces a doubled relative variation in and hence shifts the density assigned to a fixed .
Integer filling does not identify a mechanism
Section titled “Integer filling does not identify a mechanism”At integer , at least three conceptually different gaps are possible.
- A band gap appears because an integer number of one-particle subbands is filled.
- A Mott-like gap appears because changing the local integer occupation costs interaction energy even though the noninteracting active band would be partially filled.
- A flavor-ordering gap appears when interactions polarize or coherently mix internal states and thereby reconstruct the bands.
These are limiting descriptions, not always sharply separated phases. A flavor-polarized state can also have strong local charge correlations; a lattice relaxation can produce a one-body gap that interactions enhance; and a Chern insulator can arise after flavor polarization of a topological band.
Fractional filling gives a stronger hint of translation-breaking or topological order, but it is still not decisive. A larger structural supercell, single-particle miniband splitting, disorder-induced puddles, or a fractional Chern phase can all produce nontrivial fractions. The correct denominator is the capacity of the experimentally established primitive cell, not an assumed ideal cell.
The tuning space is multidimensional
Section titled “The tuning space is multidimensional”Gate voltages can vary density, displacement field, dielectric screening, and layer polarization at once. Twist angle changes the bandwidth, wave functions, topology, and relaxation. Pressure changes tunneling and structure. Magnetic field couples to spin, orbital magnetization, and Chern number. A useful phase diagram therefore reports the experimentally controlled coordinates rather than compressing them into a single nominal ratio such as .
Mott-Like and Wigner-Like Pictures
Section titled “Mott-Like and Wigner-Like Pictures”A minimal extended lattice model
Section titled “A minimal extended lattice model”When isolated moiré bands admit sufficiently localized orbitals, a useful starting point is
Assume so the hopping is Hermitian. Here label moiré orbitals or sites, label spin, valley, layer, or orbital flavors, and . The local term can include intra-orbital repulsion, interorbital repulsion, Hund-like exchange, and flavor anisotropy. Nonlocal exchange, density-assisted hopping, pair hopping, and coupling to remote bands may be important. In this partition : onsite interactions are assigned to , while contains intersite density interactions. For an unscreened long-range interaction, include the neutralizing background; equivalently use density deviations with the associated one-body terms absorbed into the potential. The Coulomb matrix depends on Wannier extent, gates, dielectric interfaces, and dynamical screening. Wu et al. (2018) construct a generalized triangular-lattice Hubbard description for a specified class of TMD moiré bands; this is a model-building example, not an identification of every moiré insulator.
This Hamiltonian is a hypothesis whose basis and truncation must be validated. Topological or fragile-topological active bands may obstruct a symmetry-preserving, exponentially localized, few-orbital representation. Even without an obstruction, broad moiré Wannier functions make neighboring interaction matrix elements comparable enough that an onsite-only model can fail.
Mott-like localization
Section titled “Mott-like localization”In the simplest site-centered picture, each effective site carries a commensurate integer occupation. Charge motion creates an energetically costly empty–double-occupation pair, while spin or valley degrees of freedom remain at lower energy. The charge density need not enlarge the moiré translation cell.
A Mott-like interpretation becomes compelling when the noninteracting active bands would be metallic, charge addition costs a finite energy, and local moments persist above the ordering temperature. These observations connect the gap to a local charge constraint. A gap that appears only when an order parameter develops may instead be largely caused by band reconstruction. The two effects can coexist, so temperature and bandwidth tuning are more informative than the label alone.
The canonical many-body charge-gap definition and its distinction from spectral, optical, and transport gaps are developed in Mott Insulators.
Generalized Wigner crystallization
Section titled “Generalized Wigner crystallization”A continuum Wigner crystal minimizes long-range Coulomb energy at low density by spontaneously forming a lattice. A generalized Wigner crystal in a moiré system instead is pinned to the pre-existing moiré lattice and forms a larger commensurate charge pattern. On a triangular lattice at , for example, occupying one of three sublattices produces a pattern with one carrier per enlarged cell.
In a classical occupancy limit,
Padhi, Chitra, and Phillips (2021) analyze Wigner crystallization in an underlying moiré potential, including the role of lattice pinning in the charge gap. In a lattice description, longer-range terms select among stripes, bubbles, triangular crystals, and other commensurate arrangements. Quantum hopping melts or dresses these patterns. Because the moiré potential already breaks continuous translation symmetry, “Wigner crystal” does not mean the same symmetry breaking as in a translationally invariant electron gas; the experimentally testable statement is enlarged-period charge order driven substantially by nonlocal repulsion.
A charge structure factor can diagnose that enlargement:
Subtracting the spatially averaged filling retains a static charge modulation. Subtracting each local mean instead would define a connected fluctuation correlator and remove the elastic Bragg contribution in a selected symmetry-broken state.
Top: occupy one of the triangular lattice’s three sublattices. The dashed charge cell has three times the primitive area. Bottom: for a perfect static pattern on periodic lattices, . Independent random occupations with probability instead give an ensemble mean at nonzero momentum. A peak that grows with sample size distinguishes coherent long-range order from local density fluctuations. These are calculated ideal limits, not measured data.
Bragg scaling of at a wavevector not equivalent to a primitive moiré reciprocal vector supports long-range charge order. A gap at fractional filling without such evidence should be called a fractional correlated insulator until its order or topology is established.
The boundary can be a crossover
Section titled “The boundary can be a crossover”“Mott” and “Wigner” emphasize different pieces of the same extended-interaction problem. At one carrier per moiré site, onsite avoidance may dominate while nonlocal repulsion reshapes the charge distribution inside each cell. At fractional filling, a charge crystal can still have strong local suppression of double occupancy. Broad orbitals may shift charge continuously between site-centered and interstitial regions.
Accordingly, labels such as Mott–Wigner state can be useful if the measured regime genuinely combines local constraints and nonlocal charge order. The name alone does not establish either contribution. The practical discriminator is whether an enlarged charge unit cell is required and observed, together with how the gap evolves when screening, bandwidth, and filling are tuned independently.
Spin and Valley Order
Section titled “Spin and Valley Order”Charge localization leaves internal degrees of freedom that can order at scales far below the charge gap. For a spin–valley spinor , convenient local operators are
acts on spin and on valley. Valley polarization has . Intervalley coherence instead uses an off-diagonal order parameter, for example
which breaks an approximate valley symmetry when intervalley scattering is weak. Layer polarization and combined spin–valley order require analogous operators. The symmetry actually available is material dependent: graphene often has approximate spin and valley flavor symmetries, while spin–orbit coupling in TMD valence bands can lock spin to valley and make the effective exchange strongly anisotropic.
For one spin-degenerate orbital per site, real nearest-neighbor hopping, one spin-1/2 fermion per site, and a stable singly occupied regime with , the simplest virtual-hopping estimate is
Here multiplies in an SU(2)-symmetric spin Hamiltonian. The intermediate empty–doubly occupied pair loses the original bond energy and gains , giving the denominator; contributions from other neighbors cancel only in a uniform environment. This is not a generic moiré exchange formula, and a multi-flavor permutation Hamiltonian uses a different normalization. Triangular geometry frustrates antiferromagnetism; complex hopping phases, multiple orbitals, direct exchange, Hund coupling, and longer-range virtual processes generate anisotropic, ring-exchange, and chiral terms. If is not the dominant excitation energy, even the denominator changes.
The hierarchy of transition scales is informative. A charge gap can open at a temperature much higher than spin or valley ordering, leaving an intermediate local-moment regime. Conversely, an insulating gap that appears only with a flavor-order parameter may be primarily Slater-like. Magnetic-field dependence alone is insufficient because Zeeman coupling, orbital moments, Chern magnetization, and field-dependent screening can all shift a gap.
Experimental Signatures
Section titled “Experimental Signatures”Establish the moiré cell and filling first
Section titled “Establish the moiré cell and filling first”Microscopy or diffraction determines the local moiré period and reconstruction; electrostatics determines the density reference and gate response. Landau fans and spectroscopy help determine the active-band capacity. These measurements should sample the region that carries the measured response, since local twist and strain can make different parts of a device reach a given filling at different gate voltages. A feature repeating every nominal suggests commensurability, but does not determine whether the electronic cell has enlarged.
Thermodynamic incompressibility is the central charge test
Section titled “Thermodynamic incompressibility is the central charge test”At zero temperature, a homogeneous charge-gapped state has fixed density over a finite interval of chemical potential. On a plot with density as the horizontal axis, this becomes a jump in chemical potential. Finite temperature and spatial variations round both descriptions. Local single-electron transistors, penetration-field capacitance, gate capacitance, and calibrated optical sensors can access related quantities. For this carrier-density convention, use the inverse density response
with a singular or strongly enhanced response at a zero-temperature charge gap. Real devices add geometric capacitance, electrostatic image charges, finite temperature, inhomogeneous broadening, and negative electronic compressibility in nearby metallic regions. These contributions must be modeled before converting a measured capacitance dip or peak into an energy gap.
In TMD moiré heterostructures, excitons in the active or a nearby sensing layer can respond sharply to electronic screening. Such optical signatures have revealed integer and fractional incompressible states even where direct electrical contacts are difficult. The exciton is a transducer: its energy or oscillator strength must be calibrated against charge response, and it does not by itself identify the electron ordering pattern.
Transport locates phases but rarely names them
Section titled “Transport locates phases but rarely names them”A thermally activated channel is often fitted as
The factor of two assumes symmetric activation of mobile particles and holes; variable-range hopping, percolation, edge conduction, contact resistance, and a temperature-dependent prefactor change the interpretation. is a transport scale, not automatically the many-body charge gap.
The most reliable transport claims reproduce filling locations across contacts and thermal cycles, map the response over displacement field and magnetic field, and compare it with a thermodynamic probe in the same density range. A resistance maximum that follows a high-resistance device region or disappears after changing contacts is weak bulk evidence.
Determine the broken symmetry
Section titled “Determine the broken symmetry”Charge order can be tested by scanning tunneling microscopy, local compressibility, microwave impedance microscopy, diffraction-sensitive methods, or spatially resolved optical response. The required observation is a reproducible modulation with a wavevector and phase relation not already imposed by the primitive structure.
Spin and valley order require probes coupled to those degrees of freedom: magnetization, susceptibility, circular dichroism, Kerr rotation, spin-resolved spectroscopy, field-angle dependence, or collective excitations. Hysteresis can support ferromagnetism but can also arise from domains and slow electrostatics. An order assignment should identify its conjugate field, symmetry, domain structure, and transition or crossover scale.
For a gapped two-dimensional phase, the Středa relation connects the density change at fixed chemical potential to the Hall response. With and signed chosen consistently with the carrier and Hall conventions,
Here is integer for an integer Chern phase; fractional Hall phases can have a rational Hall coefficient instead. A resistance-peak trajectory need not hold chemical potential fixed, so its slope is not automatically this derivative. Together with a robust Hall plateau and vanishing longitudinal response, an integer slope supports Chern number . A sloped line in a density–field map by itself is not quantization; contact mixing, parallel channels, and incomplete bulk insulation must be excluded. Fractional Hall topology requires correspondingly stronger evidence, including fractional quantization and many-body gap systematics.
Measurements are most informative in combination. A chemical-potential jump establishes a charge gap; spatial imaging distinguishes uniform filling from charge order; magnetic and optical probes constrain flavor order; Hall measurements test topology. A model should explain these observations with one physically consistent parameter set and their uncertainties, rather than requiring unrelated adjustments for each measurement.
Relation to Hubbard Models
Section titled “Relation to Hubbard Models”Moiré tunability makes Hubbard language tempting
Section titled “Moiré tunability makes Hubbard language tempting”Moiré periods enlarge electronic orbitals and reduce kinetic scales. Gates can alter dielectric screening, displacement fields can change orbital character, and twist can tune hopping. This makes the ratio of interaction to bandwidth unusually controllable. It does not make every moiré material a one-band Hubbard model.
The active orbitals must reproduce the relevant dispersion and wavefunctions. Their spatial extent sets the interaction matrix elements, while the filling convention fixes the model occupancy. Comparing a calculated optical spectrum or capacitance with a measurement also requires the corresponding observable and electrostatic environment, not only the many-body energy levels.
The hierarchy is often better described as an extended multi-flavor Hubbard family. A triangular TMD miniband may motivate a one-orbital triangular model in a restricted regime. Twisted graphene generally requires valley, sublattice, topology, remote-band, and interaction-form-factor care. At fractional filling with sizable , a lattice-gas or charge-order description may be more transparent than an onsite Hubbard reduction.
Parameters are linked to the device
Section titled “Parameters are linked to the device”With the displayed convention , Wannier matrix elements schematically obey
The minus in matches . Here contains the entire projected one-body matrix, including diagonal onsite energies through ; if an onsite one-body term is instead written separately, remove those matrix elements from so they are counted once. in this partition contains local interactions. Here is the screened interaction appropriate to the chosen low-energy partition, and only in the simplest density–density estimate. Metallic gates change the long-distance tail; layer-polarized wave functions change both hopping and interactions; and screening by remote bands can be frequency dependent. Quoting from is a scale estimate, not a complete downfolding.
A qualitative sequence of insulating fillings is a useful first test. Predicting gaps, ordering patterns and their response to screening with the same parameters is substantially more restrictive. Hubbard Physics in Materials develops this comparison between a device and its effective model.
Representative Experimental Lessons
Section titled “Representative Experimental Lessons”Twisted graphene
Section titled “Twisted graphene”The 2018 observation by Cao et al. of a half-filled insulating state near the first magic angle established that narrow graphene moiré bands can host interaction-driven phases. Subsequent compressibility, transport, and local-probe studies revealed cascades of flavor transitions, substantial device dependence, and close competition among correlated insulators, Chern states, superconductivity, and compressible metals. The early “Mott-like” analogy was productive, but a literal single-orbital Mott assignment is not universal: topology, flavor polarization, remote bands, strain, and screening all matter.
Semiconductor moiré systems
Section titled “Semiconductor moiré systems”In aligned WSe/WS, Tang et al. (2020) probed charge and magnetic response near a Mott-like filling, while Regan et al. (2020) optically detected insulating states at one hole per cell and at one-third and two-thirds filling. Related heterostructures revealed further incompressible states through optical sensing, capacitance, and transport. Integer states can realize Mott-like local constraints; fractions such as and motivate generalized Wigner or charge-order pictures because long-range repulsion can favor enlarged cells. The real-space imaging of H. Li et al. (2021) resolves charge patterns at these fractions and at one-half filling, adding spatial evidence to earlier gap-based assignments. Later observations of stripe-like patterns, tunable charge gaps, and continuous metal–insulator transitions demonstrate both the power and the nonuniversality of a single minimal model.
TMD homobilayers add twist-tunable bandwidth, displacement-field control, and strong spin–orbit-coupled valley structure. They are especially useful for testing triangular-lattice Hubbard ideas, but the relevant orbital center, valley content, and topology can change across twist and field. “TMD moiré” is a platform family, not one Hamiltonian.
Distinguishing the main possibilities
Section titled “Distinguishing the main possibilities”| Observation | What it establishes and what remains open |
|---|---|
| Integer-filling gap | Charge addition is costly; band completion, local constraints and flavor reconstruction remain possible. |
| Fractional-filling gap | The filling is special; charge order, topology or an enlarged structural cell must be distinguished. |
| Enlarged-period charge pattern | Translation symmetry is reduced; interaction and structural contributions still need comparison. |
| Activated transport | Mobile carriers face an energy scale; this need not equal the thermodynamic gap. |
Several mechanisms may coexist. A charge crystal can suppress double occupancy, and a flavor-polarized state can carry a Chern number. Moiré Topology develops the Hall tests, while Moiré Superconductivity examines neighboring paired phases.
Exercises
Section titled “Exercises”Exercise 1: converting density to filling
Section titled “Exercise 1: converting density to filling”A triangular moiré lattice has . Find its unit-cell area and the carrier density corresponding to one carrier per cell. Express the density in .
Solution
The area is
Therefore
If the local period were larger, the cell area would be about larger and the inferred one-per-cell density about smaller.
Exercise 2: commensurability versus mechanism
Section titled “Exercise 2: commensurability versus mechanism”At a device shows a sharp resistance peak. A continuum calculation predicts a single-particle gap at the same density, and the measured gap changes only weakly when a nearby screening gate is moved. What is the strongest justified claim?
Solution
The data establish a commensurate insulating transport feature consistent with a band gap. They do not require a correlated insulator because the calibrated one-body model already provides a candidate explanation and screening produces little response. A bulk thermodynamic gap measurement would strengthen the insulating assignment, but interaction necessity would still require a discrepancy with the one-body prediction or another interaction-specific signature.
Exercise 3: a fractional charge pattern
Section titled “Exercise 3: a fractional charge pattern”On a triangular moiré lattice, construct the simplest charge pattern that avoids occupied nearest neighbors. How many primitive sites lie in its charge unit cell, and what new ordering wavevector should a structure-factor measurement find?
Solution
Color the triangular lattice with three sublattices and occupy one color. The occupied sites form a triangular lattice with spacing times the original spacing. The charge cell contains three primitive moiré sites and one carrier.
The order is conventionally called . Its Bragg peaks occur at the inequivalent points of the primitive triangular-lattice Brillouin zone, modulo reciprocal lattice vectors. For a periodic sample with sites and dimensions divisible by three, the Fourier amplitude of this perfect pattern at has magnitude . Therefore . Subtracting each site’s local mean would instead remove the static pattern and miss this peak. A fractional gap alone does not demonstrate a crystal.
Exercise 4: exchange scale
Section titled “Exercise 4: exchange scale”Use the simplest one-band estimate for and . Find in millielectronvolts and kelvin. Take .
Solution
The leading estimate gives
Thus . This estimate is meaningful only if the assumed local charge excitation dominates and other orbitals, direct exchange, complex hopping, and higher-order processes are small.
Exercise 5: three gaps
Section titled “Exercise 5: three gaps”A capacitance measurement gives addition and removal edges separated by . Transport fitted over a restricted range gives , while an optical sensor shifts at the same filling. Are the three observations inconsistent?
Solution
No. The separation estimates a thermodynamic charge gap after electrostatic corrections. The smaller transport scale can reflect mobility edges, disorder-assisted hopping, or in-gap conduction. The optical shift supports a change in screening or compressibility but requires its own calibration and need not equal either energy. Agreement in filling and systematic temperature dependence is valuable even when the extracted scales differ.
Exercise 6: combining charge and Hall measurements
Section titled “Exercise 6: combining charge and Hall measurements”A fractional-filling state has a chemical-potential jump, no detected charge modulation, a Hall resistivity magnitude , longitudinal resistivity , and a density–field slope close to . Assume an isotropic homogeneous sheet. Compute and explain what the observations establish and what remains open.
Solution
The chemical-potential jump establishes bulk incompressibility at the measured resolution. Inverting the resistivity tensor gives
Thus a Hall resistivity of does not by itself imply a quantized Hall conductivity. The field slope is consistent with an integer Chern response only if it measures the fixed-chemical-potential derivative. The Hall value is suggestive. Because the longitudinal resistance is nonzero and exact Hall quantization is not established, the data do not yet prove a Chern insulator. The absence of detected charge modulation bounds order only at the probe’s resolution and wavevector coverage. Improved longitudinal insulation, metrological Hall accuracy, nonlocal or edge tests, and more complete imaging would discriminate a Chern phase from an inhomogeneous charge-ordered or parallel-conduction state.
References
Section titled “References”- Y. Cao et al., “Correlated Insulator Behaviour at Half-Filling in Magic-Angle Graphene Superlattices,” Nature 556, 80–84 (2018), doi:10.1038/nature26154.
- H. Li et al., “Imaging Two-Dimensional Generalized Wigner Crystals,” Nature 597, 650–654 (2021), doi:10.1038/s41586-021-03874-9.
- B. Padhi, R. Chitra, and P. W. Phillips, “Generalized Wigner Crystallization in Moiré Materials,” Physical Review B 103, 125146 (2021), doi:10.1103/PhysRevB.103.125146.
- E. C. Regan et al., “Mott and Generalized Wigner Crystal States in WSe/WS Moiré Superlattices,” Nature 579, 359–363 (2020), doi:10.1038/s41586-020-2092-4.
- Y. Tang et al., “Simulation of Hubbard Model Physics in WSe/WS Moiré Superlattices,” Nature 579, 353–358 (2020), doi:10.1038/s41586-020-2085-3.
- F. Wu, T. Lovorn, E. Tutuc, and A. H. MacDonald, “Hubbard Model Physics in Transition Metal Dichalcogenide Moiré Bands,” Physical Review Letters 121, 026402 (2018), doi:10.1103/PhysRevLett.121.026402.