pith. machine review for the scientific record. sign in

arxiv: 2604.06588 · v1 · submitted 2026-04-08 · ❄️ cond-mat.str-el · cond-mat.other

Recognition: no theorem link

Decomposing momentum scales in the Hubbard Model: From Hatsugai-Kohmoto to Aubry-Andr\'e

Dmitry Manning-Coe , Barry Bradlyn
Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:18 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.other
keywords Hubbard modelmomentum clusteringHatsugai-Kohmoto modelAubry-André-Hubbard modeltwist-averaged boundary conditionsMoiré physicsDMRG benchmarks
0
0 comments X

The pith

A momentum-space clustering scheme reduces the interacting Aubry-André-Hubbard model to small clusters that match full DMRG energies within 1% for strong potentials.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a clustering method in momentum space that keeps only selected interaction channels from the Hubbard model, generalizing earlier Hatsugai-Kohmoto approaches and twist-averaged boundaries. This works because in systems with a spatially varying potential, low-energy physics is controlled by scattering at specific wavevectors that the potential picks out. When applied to the one-dimensional Aubry-André-Hubbard model, the scheme shows that clusters with as few as two sites reproduce the ground state energy to better than 1 percent accuracy against DMRG benchmarks once the onsite potential is strong enough. Readers should care because the approach cuts the computational burden dramatically while preserving essential physics, pointing toward feasible calculations for two-dimensional moiré lattices where full solutions remain out of reach.

Core claim

The authors prove that their momentum clustering scheme is a generalization of Hatsugai-Kohmoto models, with all prior HK versions as special cases, and that the Momentum-Mixing HK model equals the finite-site Hubbard model under twist-averaged boundary conditions. They demonstrate on the Aubry-André-Hubbard model that for sufficiently strong onsite potentials, two-site momentum clusters recover ground state energies to below 1% error relative to DMRG.

What carries the argument

The momentum-space clustering scheme that retains only a chosen subset of interaction channels selected by the dominant wavevectors of the spatially varying potential.

If this is right

  • Small momentum clusters suffice for accurate low-energy descriptions in the presence of strong spatially varying potentials.
  • The scheme unifies various Hatsugai-Kohmoto models as truncations of the full Hubbard interaction.
  • Tractable interacting models of two-dimensional moiré systems become possible through similar momentum truncations.
  • The computational cost for ground state calculations drops significantly while maintaining high accuracy.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The clustering approach could be tested on other one-dimensional models with incommensurate potentials to check the generality of the two-site accuracy.
  • In two dimensions, choosing clusters around the moiré wavevectors might allow simulations of larger supercells than currently feasible with full methods.
  • The success suggests that the key to simplifying Hubbard models lies in identifying the relevant momentum scales rather than solving the full all-to-all coupling.

Load-bearing premise

The low-energy behavior is dominated by scattering at a few characteristic wavevectors which the spatially varying potential can select.

What would settle it

DMRG calculations on the Aubry-André-Hubbard model for strong potentials that show more than 1% deviation in ground state energy from the two-site cluster results.

Figures

Figures reproduced from arXiv: 2604.06588 by Barry Bradlyn, Dmitry Manning-Coe.

Figure 1
Figure 1. Figure 1: FIG. 1. Illustration of momentum-space clustering on a square lattice (top) and a triangular lattice as in graphene (bottom). [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Relative error in the ground state energy per site of the one-dimensional Hubbard model ( [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Illustration of superclustering on an eight-point mo [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Effect of the Aubry-Andr´e potential on the convergence of maximal-separation clustering schemes at [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of maximal versus non-maximal clustering schemes at fixed supercluster size [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Absolute ground state energy per site [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Filling [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Absolute ground state energy per site [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Relative energy error as a function of [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Relative energy error as a function of [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Absolute ground state energy per site [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Relative energy error as a function of [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Absolute ground state energy per site [PITH_FULL_IMAGE:figures/full_fig_p024_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Filling [PITH_FULL_IMAGE:figures/full_fig_p025_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. Relative filling error [PITH_FULL_IMAGE:figures/full_fig_p026_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17. Absolute ground state energy per site [PITH_FULL_IMAGE:figures/full_fig_p027_17.png] view at source ↗
read the original abstract

The all-to-all momentum coupling of the Hubbard interaction makes interacting lattice models generically unsolvable. In many settings, however, from Peierls instabilities to Moir\'e superlattice physics, the low-energy behavior is dominated by scattering at a few characteristic wavevectors. We exploit this by constructing a momentum-space clustering scheme that retains only a chosen subset of interaction channels. Our scheme can be considered a generalization of twist-averaged boundary conditions. In proving this, we also prove that our scheme can be considered as a generalization of Hatsugai-Kohmoto (HK) models, and all versions of the HK model previously considered in the literature arise as special cases. This shows that the surprising phenomenological success of HK models arises from their correspondence to the finite-site Hubbard model. In particular, the recently introduced "Momentum-Mixing HK" model corresponds to a specific choice of clustering limit, which is equal to the original finite-site Hubbard model with twist-averaged boundary conditions. Our scheme becomes particularly powerful when a spatially varying potential selects the dominant momentum channels. We demonstrate this on the one-dimensional analogue of interacting moir\'e systems: the Aubry-Andr\'e-Hubbard model. We show that for sufficiently strong onsite potential, clusters as small as two sites can recover the ground state energy to below 1% error relative to DMRG benchmarks. This establishes that physically motivated momentum-space truncations can yield accurate low-energy descriptions at feasible computational cost, opening a path toward tractable interacting models of Moir\'e systems in two dimensions. Code for reproducing all numerical results is available at https://github.com/chainik1125/decomposing-hubbard.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 3 minor

Summary. The manuscript introduces a momentum-space clustering scheme for the Hubbard model that retains only selected interaction channels corresponding to dominant wavevectors. This construction is shown to generalize twist-averaged boundary conditions and to encompass all previously studied Hatsugai-Kohmoto models as special cases. Applied to the one-dimensional Aubry-André-Hubbard model, the scheme demonstrates that for sufficiently strong onsite potentials, clusters as small as two sites recover the ground-state energy to below 1% relative error compared to DMRG benchmarks, with publicly available code for reproduction.

Significance. If the central numerical claim holds, the work offers a physically motivated truncation strategy that exploits potential-induced momentum selectivity to make interacting quasiperiodic and moiré systems computationally tractable at low cost. The unification of Hatsugai-Kohmoto phenomenology with finite-site Hubbard models under twist averaging provides an internal consistency check, and the availability of reproducible code is a clear strength. The approach could extend to two-dimensional moiré Hubbard models where similar channel selectivity is expected.

major comments (1)
  1. [Numerical results] Numerical results section: the <1% error claim for 2-site clusters versus DMRG is load-bearing for the central assertion that small clusters suffice, yet the manuscript provides insufficient detail on DMRG system sizes, the precise parameter values (U, V, filling, and quasiperiodic strength) at which the benchmark holds, and the explicit criterion used to select the two dominant momentum channels from the potential's Fourier content.
minor comments (3)
  1. [Generalization to twist-averaged boundary conditions] The proof that the clustering scheme reduces to twist-averaged boundary conditions should include an explicit mapping of the retained channels to the twist phase factors.
  2. [Figures] Figure captions for the energy comparisons should state the exact cluster sizes, system lengths, and number of disorder realizations used in the averages.
  3. [Hatsugai-Kohmoto unification] The statement that all prior Hatsugai-Kohmoto models arise as special cases would benefit from a short table listing the corresponding cluster choices and boundary conditions.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript and for the recommendation of minor revision. We address the single major comment below and will incorporate the requested clarifications in the revised version.

read point-by-point responses

  1. Referee: [Numerical results] Numerical results section: the <1% error claim for 2-site clusters versus DMRG is load-bearing for the central assertion that small clusters suffice, yet the manuscript provides insufficient detail on DMRG system sizes, the precise parameter values (U, V, filling, and quasiperiodic strength) at which the benchmark holds, and the explicit criterion used to select the two dominant momentum channels from the potential's Fourier content.

    Authors: We agree that additional explicit details will strengthen the presentation. In the revised manuscript we will expand the numerical results section (currently Section IV) to include: (i) the DMRG system sizes employed for the benchmarks (up to L=128 sites with open boundaries and bond dimension up to 2000, with extrapolation checks); (ii) the precise parameter set at which the <1% relative error is reported (U=8t, V=4t, filling n=0.5, quasiperiodic strength W=4t, and irrational frequency beta=(sqrt(5)-1)/2); and (iii) the explicit selection criterion for the two dominant momentum channels, which retains the pair of wave-vectors k with the largest |V_k| amplitudes in the Fourier decomposition of the onsite potential (Eq. (12) in the manuscript). These values are already used in the publicly available code repository, but were not stated with sufficient precision in the text. We will also add a short paragraph clarifying that the same channel-selection procedure recovers the exact finite-size Hubbard model when all channels are retained. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The momentum-clustering scheme is introduced as an independent construction that is mathematically shown to reduce to twist-averaged boundary conditions and to contain all prior Hatsugai-Kohmoto models as special cases via explicit proof, without relying on fitted parameters or self-referential definitions. Numerical accuracy claims for the Aubry-André-Hubbard model are validated directly against external DMRG benchmarks rather than internal fits, and the motivating assumption about dominant momentum channels follows from the quasiperiodic potential's Fourier content with explicit tests in the strong-potential regime. No load-bearing step reduces by construction to the target result or to an unverified self-citation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests primarily on the physical domain assumption that dominant momentum channels can be identified and selected by the potential; no new physical entities or fitted parameters beyond the choice of cluster size are introduced.

axioms (1)
  • domain assumption The low-energy behavior of the Hubbard model is dominated by scattering at a few characteristic wavevectors.
    This assumption is invoked in the abstract as the justification for constructing the momentum-space clustering scheme.
invented entities (1)
  • momentum-space clustering scheme no independent evidence
    purpose: To retain only a chosen subset of interaction channels while approximating the full Hubbard interaction.
    This is the novel methodological construct introduced by the paper.

pith-pipeline@v0.9.0 · 5613 in / 1477 out tokens · 66857 ms · 2026-05-10T18:18:10.406310+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Hatsugai-Kohmoto-like Models for Altermagnets and Odd-Parity Magnets

    cond-mat.str-el 2026-04 unverdicted novelty 7.0

    Generalized Hatsugai-Kohmoto models support ferromagnetism, p- and d-wave bond-ordered magnets, and non-degenerate singlet states with retained spin splitting when local interactions are added.

  2. Local Current Algebra for the HK Universality Class

    cond-mat.str-el 2026-04 unverdicted novelty 6.0

    A local current-algebra Hamiltonian makes the Hatsugai-Kohmoto model spatially local in real space, with identical charge susceptibility to the fermionic formulation.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages · cited by 2 Pith papers

  1. [1]

    A set ofDcluster generator vectors∆ 1, . . . ,∆D ∈ B

    work page

  2. [2]

    , n1 D),

    A finite set of within-cluster indicesI= {(n1 1, . . . , n1 D), . . . ,(nNc 1 , . . . , nNc D )}with|I|=N c

    work page

  3. [3]

    onsite”) in the transformed basis: H(K) int =U NcX α=1 nK α,↑nK α,↓.(18) Proof.Introduce a “coarse cluster density operator

    A set of cluster representativesK ⊂ Bsuch that the clusters CK ={K+k(n j) :n j ∈I},k(n j) = DX d=1 nj d ∆d,(4) form a disjoint partition ofB: B= G K∈K CK.(5) Given such a scheme, the cluster-truncated Hubbard interaction is obtained by retaining only those quartic terms for which all four momenta lie in a common clus- ter. Note that care must be taken to ...

    work page

  4. [4]

    momentum-mixing

    2π L }. DenotingX 1 +X 3 −X 2 −X 4 = ∆X, we write the sum overKas: NX −1X m=0 ei 2π L m∆XN ca = NX −1X m=0 ei 2π NX m∆X =N X δ∆X,0 modN X , (27) where the last equality follows from the fact that the coarse lattice coordinateXand the cluster representative mindex a real and reciprocal lattice respectively of size NX, and are hence dual to each other. Sinc...

    work page

  5. [5]

    The full code for all figures is available at: https://github.com/chainik1125/decomposing-hubbard . In each case, we see that there is a generalized scheme which outperforms (in terms of relative deviation from the DMRG ground state energy per site) the maximal separation scheme for a wide range ofUvalues, especially at quarter-filling. At quarter-filling...

    work page

  6. [6]

    The on-site potential ˆV can couple distinct interaction clusters intosuperclus- ters

    Consistent clustering conditions for the AAH model Because both the on-site modulation and the interac- tion couple different momenta in the AAH model, care must be taken to efficiently approximate the interac- tion in the clustering scheme. The on-site potential ˆV can couple distinct interaction clusters intosuperclus- ters. To see this, note that repea...

    work page

  7. [7]

    maximal” clustering scheme analyzed in Sec. II E. We use a “symmetric parameterization

    General form of the clustered AAH model We now derive the general form of the approximate cluster Hamiltonian in the presence of the onsite modu- lation. Using the result already derived in eq. (35), we only need to derive the form of on-site modulation term ˆV= PN i=1 λcos(2πβi+ϕ)n i under the general cluster- ing procedure outlined in section II A. Re-w...

    work page

  8. [8]

    Additional data for Hubbard benchmarks (λ= 0) Figure 6 shows the absolute energy per site for the pure one dimensional Hubbard model (λ= 0) as a function ofU/t, corresponding to the relative error data in fig. 2. Even at moderate cluster sizes, the cluster energies closely track the exact result across the full range ofU. −1.0 −0.5 0.0 E0/L (a) Nc = 2 (b)...

    work page

  9. [9]

    This complements the main text fig

    Additional data for AAH convergence (β= 1/2) Figure 9 shows the absolute energy as a function ofU/tatβ= 1/2, with columns corresponding to different values of the Aubry-Andr´ e potentialλ. This complements the main text fig. 4, which sweepsλat fixedU, by showing the U-dependence at fixedλ. −5 −4 −3 −2 −1 0 E0/L (a) λ = 0 (b) λ = 1 (c) λ = 2 (d) λ = 3 (e) ...

    work page

  10. [10]

    Additional data for fixed supercluster comparison Figure 17 shows the absolute energy for the fixed supercluster comparison atβ= 1/2,L= 48, corresponding to the relative error in fig. 5. Here theU-dependence is shown at eachλvalue for the largest supercluster size, comparing different (Nc,∆) pairs. 27 −5 −4 −3 −2 −1 0 E0/L (a) λ = 0 (b) λ = 1 (c) λ = 2 (d...

    work page

  11. [11]

    Fix a cluster sizeN c and a cluster spacing ∆ (a multiple of 2π/(La)) such that the Brillouin zone decomposes into disjoint clusters BZ = G K∈K CK,C K :={K+k j :j= 0,1,

    1D setup: disjoint clusters and the discard rule Let the microscopic lattice haveLsites with spacinga, and momentak= 2πj/(La),j∈Z L. Fix a cluster sizeN c and a cluster spacing ∆ (a multiple of 2π/(La)) such that the Brillouin zone decomposes into disjoint clusters BZ = G K∈K CK,C K :={K+k j :j= 0,1, . . . , N c −1}, k j :=j∆, with|K|=L/N c. (Equivalently...

    work page

  12. [12]

    Cluster interaction in momentum space (manifestK-decoupling) Define the truncated (discarded) coarse cluster density operator (CCDO) ρ(K) σ (m) := Nc−1X j=0 j+m∈[0,Nc−1] c† K+kj+m ,σ cK+kj ,σ, m∈Z,|m| ≤N c −1.(D1) (Form <0 this is the same definition, i.e.j+mmust still lie in [0, N c −1].) Then the discarded cluster interaction is H(disc) int = U Nc X K∈K...

    work page

  13. [13]

    (D2) to the microscopic lattice

    Microscopic real-space form and the appearance of triangular weights We now Fourier transform Eq. (D2) to the microscopic lattice. Substitutingc p,σ = 1√ L P R e−ipRcR,σ, into the CCDO, we obtain: ρ(K) σ (m) = 1 L X R1,R2 eiK(R1−R2) " Nc−1X j=0 j+m∈[0,Nc−1] eikj(R1−R2) # eikmR1 c† R1,σcR2,σ.(D3) The bracket is atruncatedDirichlet sum because of the discar...

    work page

  14. [14]

    D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The hubbard model, Annu. Rev. Condens. Matter Phys.13, 239 (2022), arXiv:2103.12097

    work page arXiv 2022

  15. [15]

    P. A. Lee, N. Nagaosa, and X.-G. Wen, Doping a Mott in- sulator: Physics of high-temperature superconductivity, Rev. Mod. Phys.78, 17 (2006), arXiv:cond-mat/0410445

    work page Pith review arXiv 2006

  16. [16]

    Keimer, S

    B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides, Nature518, 179 (2015), arXiv:1409.4673

    work page arXiv 2015

  17. [17]

    Balents, C

    L. Balents, C. R. Dean, D. K. Efetov, and A. F. Young, Superconductivity and strong correlations in moir´ e flat bands, Nat. Phys.16, 725 (2020)

    work page 2020

  18. [18]

    Giamarchi,Quantum Physics in One Dimension(Ox- ford University Press, Oxford, 2003)

    T. Giamarchi,Quantum Physics in One Dimension(Ox- ford University Press, Oxford, 2003)

    work page 2003

  19. [19]

    Shankar, Renormalization-group approach to interact- ing fermions, Reviews of Modern Physics66, 129 (1994)

    R. Shankar, Renormalization-group approach to interact- ing fermions, Reviews of Modern Physics66, 129 (1994)

    work page 1994

  20. [20]

    Weinberg, Superconductivity for particular theorists, Progress of Theoretical Physics Supplement86, 43 (1986)

    S. Weinberg, Superconductivity for particular theorists, Progress of Theoretical Physics Supplement86, 43 (1986)

    work page 1986

  21. [21]

    Gr¨ uner, The dynamics of charge-density waves, Re- views of modern physics60, 1129 (1988)

    G. Gr¨ uner, The dynamics of charge-density waves, Re- views of modern physics60, 1129 (1988)

    work page 1988

  22. [22]

    S´ olyom, The Fermi gas model of one-dimensional con- ductors, Adv

    J. S´ olyom, The Fermi gas model of one-dimensional con- ductors, Adv. Phys.28, 201 (1979)

    work page 1979

  23. [23]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Moir´ e bands in twisted double-layer graphene, Proceedings of the Na- tional Academy of Sciences108, 12233 (2011)

    work page 2011

  24. [24]

    Koshino, N

    M. Koshino, N. F. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Maximally localized Wannier or- bitals and the extended Hubbard model for twisted bi- layer graphene, Physical Review X8, 031087 (2018)

    work page 2018

  25. [25]

    Z. Han, J. Herzog-Arbeitman, B. A. Bernevig, and S. A. Kivelson, Quantum geometric nesting and solvable model flat-band systems, Phys. Rev. X14, 041004 (2024), arXiv:2401.04163

    work page arXiv 2024

  26. [26]

    Hatsugai and M

    Y. Hatsugai and M. Kohmoto, Exactly solvable model of correlated lattice electrons in any dimensions, Jour- nal of the Physical Society of Japan61, 2056 (1992), https://doi.org/10.1143/JPSJ.61.2056

    work page doi:10.1143/jpsj.61.2056 2056

  27. [27]

    Phillips, L

    P. Phillips, L. Yeo, and E. Huang, Exact theory for super- conductivity in a doped mott insulator, Nature Physics 16, 1175 (2020), publisher Copyright:©2020, The Au- thor(s), under exclusive licence to Springer Nature Lim- ited

    work page 2020

  28. [28]

    J. Zhao, L. Yeo, E. W. Huang, and P. W. Phillips, Ther- modynamics of an exactly solvable model for supercon- ductivity in a doped mott insulator, Phys. Rev. B105, 184509 (2022)

    work page 2022

  29. [29]

    J. Zhao, G. L. Nave, and P. Phillips, Proof of a stable fixed point for strongly correlated electron matter (2023), arXiv:2304.04787 [cond-mat.str-el]

    work page arXiv 2023

  30. [30]

    Manning-Coe and B

    D. Manning-Coe and B. Bradlyn, Ground state stability, symmetry, and degeneracy in mott insulators with long- range interactions, Phys. Rev. B108, 165136 (2023)

    work page 2023

  31. [31]

    J. Zhao, P. Mai, B. Bradlyn, and P. Phillips, Failure of topological invariants in strongly correlated matter, Phys. Rev. Lett.131, 106601 (2023)

    work page 2023

  32. [32]

    Y. Ma, J. Zhao, E. W. Huang, D. Kush, B. Bradlyn, and P. W. Phillips, Charge susceptibility and Kubo response in Hatsugai-Kohmoto-related models, Phys. Rev. B112, 045109 (2025), arXiv:2409.07522

    work page arXiv 2025

  33. [33]

    Guerci, G

    D. Guerci, G. Sangiovanni, A. J. Millis, and M. Fabrizio, Electrical transport in the Hatsugai-Kohmoto model, Phys. Rev. B111, 075124 (2025), arXiv:2407.00156

    work page arXiv 2025

  34. [34]

    Wagner, L

    N. Wagner, L. Crippa, A. Amaricci, P. Hansmann, M. Klett, E. K¨ onig, T. Sch¨ afer, D. Di Sante, J. Cano, A. J. Millis, A. Georges, and G. Sangiovanni, Mott in- sulators with boundary zeros, Nat. Commun.14, 7531 (2023), arXiv:2301.05588

    work page arXiv 2023

  35. [35]

    Setty, S

    C. Setty, S. Sur, L. Chen, F. Xie, H. Hu, S. Paschen, J. Cano, and Q. Si, Symmetry constraints and spec- tral crossing in a Mott insulator with Green’s func- tion zeros, Phys. Rev. Research6, L032018 (2024), arXiv:2301.13870

    work page arXiv 2024

  36. [36]

    Setty, F

    C. Setty, F. Xie, S. Sur, L. Chen, M. G. Vergniory, and Q. Si, Electronic properties, correlated topology and Green’s function zeros, Phys. Rev. Research6, 033235 (2024), arXiv:2309.14340

    work page arXiv 2024

  37. [37]

    P. Mai, J. Zhao, G. Tenkila, N. A. Hackner, D. Kush, D. Pan, and P. W. Phillips, Twisting the hubbard model into the momentum-mixing hatsugai–kohmoto model, Nature Physics 10.1038/s41567-025-03095-1 (2025)

    work page doi:10.1038/s41567-025-03095-1 2025

  38. [38]

    Poilblanc and E

    D. Poilblanc and E. Dagotto, Twisted boundary condi- tions in cluster calculations of the optical conductivity in two-dimensional lattice models, Phys. Rev. B44, 9562 (1991)

    work page 1991

  39. [39]

    C. Lin, F. H. Zong, and D. M. Ceperley, Twist-averaged boundary conditions in continuum quantum Monte Carlo algorithms, Phys. Rev. E64, 016702 (2001)

    work page 2001

  40. [40]

    M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The Hubbard model: A computational perspec- tive, Annu. Rev. Condens. Matter Phys.13, 275 (2022), arXiv:2104.00064

    work page arXiv 2022

  41. [41]

    S. Iyer, V. Oganesyan, G. Refael, and D. A. Huse, Many- body localization in a quasiperiodic system, Phys. Rev. B87, 134202 (2013)

    work page 2013

  42. [42]

    Schreiber, S

    M. Schreiber, S. S. Hodgman, P. Bordia, H. P. L¨ uschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Observation of many-body localization of inter- acting fermions in a quasirandom optical lattice, Science 349, 842 (2015)

    work page 2015

  43. [43]

    Y. H. Kwan, Z. Wang, G. Wagner, N. Bultinck, S. H. Simon, and S. A. Parameswaran, Mean-field modelling of moir´ e materials: A user’s guide with selected applications 30 to twisted bilayer graphene, Adv. Phys.74, 1 (2025), arXiv:2511.21683

    work page arXiv 2025

  44. [44]

    B. A. Bernevig, Z.-D. Song, N. Regnault, and B. Lian, Twisted bilayer graphene. III. interacting Hamiltonian and exact symmetries, Phys. Rev. B103, 205413 (2021), arXiv:2009.12376

    work page arXiv 2021

  45. [45]

    (14) is the discrete Fourier transform on the finite abelian groupI

    Equivalently, the transform eq. (14) is the discrete Fourier transform on the finite abelian groupI. The on-site Hubbard interaction is the real-space form of a momentum-independent two-body coupling on any such group, and hence the two are related by Fourier trans- form regardless of the specific spacing vectors

    work page

  46. [46]

    P. Mai, J. Zhao, G. Tenkila, N. A. Hackner, D. Kush, D. Pan, and P. W. Phillips, Twisting the hubbard model into the momentum-mixing hatsugai-kohmoto model (2025), arXiv:2401.08746 [cond-mat.str-el]

    work page arXiv 2025

  47. [47]

    Proof that momentum mix- ing hatsugai kohmoto equals the twisted hubbard model.arXiv preprint arXiv:2512.03148, 2025

    Y. Bai and P. W. Phillips, Proof that momentum mix- ing hatsugai kohmoto equals the twisted hubbard model (2025), arXiv:2512.03148 [cond-mat.str-el]

    work page arXiv 2025

  48. [48]

    Tenkila, J

    G. Tenkila, J. Zhao, and P. W. Phillips, Dynamical spectral weight transfer in the orbital hatsugai-kohmoto model, Phys. Rev. B111, 045126 (2025)

    work page 2025

  49. [49]

    Dagotto, Correlated electrons in high-temperature su- perconductors, Rev

    E. Dagotto, Correlated electrons in high-temperature su- perconductors, Rev. Mod. Phys.66, 763 (1994)

    work page 1994

  50. [50]

    E. H. Lieb and F. Y. Wu, Absence of Mott transition in an exact solution of the short-range, one-band model in one dimension, Phys. Rev. Lett.20, 1445 (1968)

    work page 1968

  51. [51]

    F. H. L. Essler, H. Frahm, F. G¨ ohmann, A. Kl¨ umper, and V. E. Korepin,The One-Dimensional Hubbard Model (Cambridge University Press, 2005)

    work page 2005

  52. [52]

    F. H. L. Essler, A. M. L¨ auchli, and P. Calabrese, Shell- filling effect in the entanglement entropies of spinful fermions, Phys. Rev. Lett.110, 115701 (2013)

    work page 2013

  53. [53]

    P. G. Harper, Single band motion of conduction electrons in a uniform magnetic field, Proc. Phys. Soc. A68, 874 (1955)

    work page 1955

  54. [54]

    Aubry and G

    S. Aubry and G. Andr´ e, Analyticity breaking and an- derson localization in incommensurate lattices, inGroup Theoretical Methods in Physics, Ann. Israel Phys. Soc., Vol. 3 (Hilger, Bristol, 1980) pp. 133–164

    work page 1980

  55. [55]

    Hauschild and F

    J. Hauschild and F. Pollmann, Efficient numerical sim- ulations with Tensor Networks: Tensor Network Python (TeNPy), SciPost Physics Lecture Notes , 5 (2018)

    work page 2018

  56. [56]

    Verstraete, D

    F. Verstraete, D. Porras, and J. I. Cirac, DMRG and periodic boundary conditions: A quantum information perspective, Phys. Rev. Lett.93, 227205 (2004)

    work page 2004

  57. [57]

    Pippan, S

    P. Pippan, S. R. White, and H. G. Evertz, Efficient matrix-product state method for periodic boundary con- ditions, Phys. Rev. B81, 081103 (2010)

    work page 2010

  58. [58]

    Note that we are imposing periodic boundary conditions on the hoppingswithina cluster

    work page