Quantum-Geometric Design of Lattice Generalized Landau Levels

Abstract: We design lattice models with tailored quantum geometry, including generalized Landau levels (LLs) satisfying the integrated trace condition and higher-Chern bands with ideal quantum geometry. Our models with $N=2$, $3$, and $4$ sublattices include a generalized Haldane model ($N=2$ honeycomb lattice model) with Gaussian-decaying hoppings realizable in twisted bilayer MoTe$_2$, and $N \geq 3$ models with exponentially decaying hoppings. Exact diagonalization reveals fractional Chern insulators in the generalized zeroth LL bands of all three models, a Moore-Read state in the generalized first LL band of the $N=4$ model, and various interaction-driven topological phases$\unicode{x2013}$including integer and fractional anomalous Hall crystals and a multicomponent Halperin state$\unicode{x2013}$in the ideal higher-Chern band of the $N=3$ model. Informed by quantum geometry, our work provides a pathway for lattice realizations of Landau-level and beyond-Landau-level physics.

• The paper develops a density-modulated basis and Gram–Schmidt method that creates ideal generalized Landau-level bands on lattices with 2–4 sublattices, including Chern numbers −1, −2, and −3.
• The construction produces analytically controlled hopping models, unifies Gaussian-decay bands with generalized Kapit–Mueller models, and quantitatively reproduces effective honeycomb hoppings in twisted bilayer MoTe₂.
• Exact diagonalization reveals Moore–Read, Halperin, Laughlin, integer anomalous Hall crystal, and fractional anomalous Hall crystal states, while showing that linear independence and emergent higher-band symmetries require lattice- and filling-specific checks.

Overview

This paper constructs explicit lattice Hamiltonians whose single-particle bands realize generalized Landau levels (gLLs) — lattice analogues of the $n$th Landau level (nLL) with ideal quantum geometry — on non-Bravais lattices with $N = 2, 3, 4$ sublattices. The central construction uses a "density-modulated basis" built from magnetic Bloch wave functions sampled at sublattice positions, followed by Gram–Schmidt orthogonalization. The resulting bands include a generalized 0LL band of Chern number $\mathcal{C} = -1$, an ideal higher-Chern band with $\mathcal{C} = -(N-1)$, and, for $N \geq 3$, isolated generalized 1LL bands supporting Moore–Read-type physics at half filling. The authors further demonstrate that the effective honeycomb model of twisted bilayer MoTe$_2$ (tMoTe$_2$) quantitatively matches their $N=2$ construction, and they use exact diagonalization (ED) to identify integer and fractional anomalous Hall crystal states and fractionalized phases in the engineered bands.

Magnetic Bloch wave functions and the density-modulated basis

The construction starts from the magnetic Bloch wave function $\Psi_{n,\mathbf{k}}(\mathbf{r})$ for the nLL on a torus, expressed via Haldane's modified Weierstrass sigma function (2607.08702). These wave functions satisfy magnetic translational symmetry, position–momentum duality, and quasi-periodicity in momentum space. The key object is the density-modulated basis state

$$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$

where $N = 2, 3, 4$0 is a site-dependent gauge factor and $N = 2, 3, 4$1 are sublattice positions. For an $N = 2, 3, 4$2-sublattice model, applying Gram–Schmidt orthogonalization to $N = 2, 3, 4$3 yields $N = 2, 3, 4$4 orthonormal states $N = 2, 3, 4$5 spanning a flat band with ideal quantum geometry; the remaining state $N = 2, 3, 4$6, fixed by completeness, forms an ideal higher-Chern band.

The authors are careful to note that this procedure requires the density-modulated states to be linearly independent — verified explicitly via positive-definite Gram matrices for their $N = 2, 3, 4$7 models — but that the condition is not automatic: on the kagome lattice, $N = 2, 3, 4$8 for all three sublattices, so $N = 2, 3, 4$9 vanishes and the construction fails for the 1LL. This is an important caveat: each lattice realization must be checked individually.

Ideal higher Chern band: proof of holomorphicity

The paper proves that $\mathcal{C} = -1$0 has ideal quantum geometry by showing its momentum-space coefficients are holomorphic in $\mathcal{C} = -1$1 up to a Gaussian normalization factor. The proof exploits two facts: (i) the Gram–Schmidt transformation matrix $\mathcal{C} = -1$2 is lower-triangular with diagonal entries $\mathcal{C} = -1$3 independent of $\mathcal{C} = -1$4, so $\mathcal{C} = -1$5 is $\mathcal{C} = -1$6-independent; and (ii) $\mathcal{C} = -1$7 can be expanded in $\mathcal{C} = -1$8-derivatives of $\mathcal{C} = -1$9, making the coefficients $\mathcal{C} = -(N-1)$0 holomorphic after stripping the Gaussian factor. Consequently,

$\mathcal{C} = -(N-1)$1

with $\mathcal{C} = -(N-1)$2 holomorphic, which satisfies the ideal trace condition $\mathcal{C} = -(N-1)$3 exactly. This establishes that the top band carries Chern number $\mathcal{C} = -(N-1)$4 with perfectly uniform quantum metric trace equal to the Berry curvature magnitude — the defining property of an ideal flatband suitable for exact fractionalized parent states.

Explicit hopping models

The hopping parameters are obtained analytically by Fourier transforming the Hamiltonian matrix elements. For the Gaussian-decay model ($\mathcal{C} = -(N-1)$5, honeycomb), where the generalized 0LL is a zero-energy band and all higher levels are degenerate, the hoppings take closed form:

$\mathcal{C} = -(N-1)$6

with $\mathcal{C} = -(N-1)$7 decaying as $\mathcal{C} = -(N-1)$8. A notable structural result is that this Gaussian-decay model is equivalent to the generalized Kapit–Mueller model on non-Bravais lattices [Dong2020Exact] up to a site-dependent gauge transformation — unifying two previously distinct exact-flatband constructions. Truncation to next-nearest neighbors recovers the Haldane model.

For $\mathcal{C} = -(N-1)$9 and $N \geq 3$0, the energy spectrum is chosen to lift the degeneracy and isolate the generalized 1LL band. The resulting hoppings decay exponentially rather than Gaussianly (shown on log plots), with dominant first- through fourth-neighbor terms. The models respect explicit point-group symmetries ($N \geq 3$1/$N \geq 3$2 for $N \geq 3$3; $N \geq 3$4/$N \geq 3$5 for $N \geq 3$6). Numerically, truncating hoppings beyond $N \geq 3$7 preserves the band structures semi-quantitatively and retains nearly ideal quantum geometry, indicating practical short-range realizations.

Connection to twisted bilayer MoTe$N \geq 3$8

A strong quantitative claim of the paper is that tMoTe$N \geq 3$9 at magic angle realizes the $_2$0 model. Constructing layer-polarized Wannier functions from the continuum model and projecting onto them yields an effective honeycomb tight-binding model. Comparing relative hoppings $_2$1 against the model's $_2$2:

Parameter set	$_2$3	$_2$4	$_2$5	$_2$6
Wu et al.	$_2$7	$_2$8	$_2$9	$_2$0
Reddy et al.	$_2$1	$_2$2	$_2$3	$_2$4
Wang et al.	$_2$5	$_2$6	$_2$7	$_2$8
$_2$9 model	$N=2$0	$N=2$1	$N=2$2	$N=2$3

The match is quantitative for $N=2$4 across three independent parameter sets (magic angles $N=2$5, $N=2$6, $N=2$7). This implies that the topmost moiré valence band of tMoTe$N=2$8 — known to have nearly ideal quantum geometry and near-vanishing bandwidth — is effectively a lattice realization of the generalized 0LL, providing microscopic justification for mapping fractional Chern insulator physics in tMoTe$N=2$9 onto Landau-level phenomenology.

Many-body phases from exact diagonalization

The ED results establish several correlated phases in the engineered bands:

• Generalized 1LL at $\Psi_{n,\mathbf{k}}(\mathbf{r})$0 ($\Psi_{n,\mathbf{k}}(\mathbf{r})$1): PES shows an entanglement gap with Moore–Read counting (13338 and 18571 below-gap levels for $\Psi_{n,\mathbf{k}}(\mathbf{r})$2 clusters), evidencing a non-Abelian paired state.
• $\Psi_{n,\mathbf{k}}(\mathbf{r})$3 band ($\Psi_{n,\mathbf{k}}(\mathbf{r})$4): Halperin (332)-type state at $\Psi_{n,\mathbf{k}}(\mathbf{r})$5 with $\Psi_{n,\mathbf{k}}(\mathbf{r})$6 GPP counting of 1360, and a Laughlin-type state at $\Psi_{n,\mathbf{k}}(\mathbf{r})$7 with $\Psi_{n,\mathbf{k}}(\mathbf{r})$8 counting of 2280.
• Integer anomalous Hall crystal (AHC): at $\Psi_{n,\mathbf{k}}(\mathbf{r})$9 in the ideal higher Chern band, constructed analytically as product states of filled ideal subbands with emergent SU(2) structure parameterized on a Bloch sphere; these variational states achieve >99% overlap weight with the ED ground-state manifold.
• $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$0 band ($$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$1): at $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$2, a threefold quasi-degenerate gapped state with many-body Chern number $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$3 (integer AHC); at $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$4, a sevenfold manifold with $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$5 (multicomponent Halperin); at $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$6, a 30-fold (84-fold for $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$7) degeneracy consistent with SU(3) × Laughlin product structure, i.e., a fractional AHC with $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$8.

The AHC construction itself is technically notable: the ideal subbands created by $$|e_{n,\mathbf{k}}\rangle = \sum_s \mathcal{B}(\boldsymbol{\tau}_s)\Psi_{n,\mathbf{k}}(\boldsymbol{\tau}_s)|\psi_{\mathbf{k},s}\rangle,$$9 remain anti-holomorphic in $N = 2, 3, 4$00 up to normalization, so the charge-ordered crystals inherit ideal quantum geometry — extending the ideal-flatband program from fluids to broken-symmetry crystalline states.

Limitations and open questions

Several caveats qualify the results. First, the linear independence of density-modulated basis states is assumption-dependent and fails on specific lattices (kagome for the 1LL), so the framework does not universally apply. Second, the emergence of SU(3) symmetry in the $N = 2, 3, 4$01 band appears filling-dependent: it is well developed at $N = 2, 3, 4$02 but evidently strongly broken at $N = 2, 3, 4$03, where only threefold rather than $N = 2, 3, 4$04-fold degeneracy is observed — the mechanism behind this asymmetry is not resolved. Third, the equivalence between the Gaussian-decay and Kapit–Mueller models holds for one flux quantum per unit cell, and the analytic hopping formulas apply strictly to the zero-energy-band case; the $N = 2, 3, 4$05 hoppings require numerical Fourier transforms. Finally, the tMoTe$N = 2, 3, 4$06 comparison is made at the single-particle level; whether interaction-driven phases in the material map one-to-one onto those found in the ideal models remains an open question.

Conclusion

This work provides a systematic, analytically controlled recipe for engineering lattice bands that reproduce generalized Landau levels with provably ideal quantum geometry, unifying the Gaussian-decay and Kapit–Mueller constructions and demonstrating quantitative realization in twisted bilayer MoTe$N = 2, 3, 4$07. The combination of exact band design with ED evidence for Moore–Read, Halperin, Laughlin, and (fractional) anomalous Hall crystal phases makes these models a controlled setting for studying zero-field fractional quantum Hall physics, while leaving open the symmetry-breaking mechanisms governing higher-Chern-band crystal states at partial fillings.