Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2503.10760 · Physical Review A

Solving the SU( N ) orbital Hatsugai-Kohmoto model

Nico A. Hackner, Peizhi Mai, Philip Phillips

Published May 29, 2025• 1 citation
arXivHTMLPDF

Abstract

We show that the physics of the SU($N$) Hubbard model can be realistically simulated with the recently developed orbital Hatsugai-Kohmoto model. In this approach, the momentum mixing absent from the band Hatsugai-Kohmoto model is included by grouping $n$-Hubbard atoms into a cluster. We take advantage of the rapid convergence of this scheme ($1/n^2$) and show that for $n\approx O(10)$, quantitative agreement with the determinantal quantum Monte Carlo method arises for the double occupancy as well as qualitative agreement for the filling and compressibility across the Mott transition. Additional features of this work are that we can obtain low-temperature physics, which could be qualitatively different from its high-temperature counterpart, and dynamical quantities without resorting to analytical continuation, thereby establishing the orbital Hatsugai-Kohmoto model as a useful alternative perspective for studying strong correlations in SU($N$) systems.

DOI
10.1103/1qzc-6c6n

Related papers

Dynamical spectral weight transfer in the orbital Hatsugai-Kohmoto model

Gaurav Tenkila, Jinchao Zhao, Philip Phillips

The shift in the spectral weight distribution under doping is termed the spectral weight transfer (SWT) and is a key indicator of strongly correlated physics. The static SWT ($t=0$) for doping fraction $x$ scales as $x$ for semiconductors and as $2x$ for Mott insulators. The SWT receives corrections in $t/U$ due to the noncommutativity of the kinetic and interaction terms in the Hamiltonian, termed the dynamic SWT. We compute explicitly the low-energy spectral weight (LESW) in the lower Hubbard band using the orbital Hatsugai-Kohmoto (OHK) model. We show that the dynamical mixing enters the LESW with a universal slope independent of the number of orbitals. This independence on orbital number corroborates the rapid convergence of the OHK model to the Hubbard model. As a result, OHK is an effective simulator of Hubbard physics.

Physical review. B./Physical review. B•Jan 10, 2025•7 citations

Superconductivity in a Two-Orbital Hatsugai-Kohmoto Model at Half Filling

Nico Hahn, R. Matthias Geilhufe

Multiband superconductivity gives rise to a rich landscape of possible pairing states. Here we study superconductivity emerging from a half-filled normal state in a two-orbital extension of the Hatsugai-Kohmoto model, an exactly solvable model of correlated electrons with momentum-local interactions, which provides a minimal framework to explore the interplay of strong correlations, orbital structure and pairing symmetry. For a system with point-group symmetry $\rm D_{4h}$, we classify the symmetry-allowed superconducting gap structures, taking into account spin, orbital and momentum degrees of freedom. We further compute the critical temperature and the superconducting order parameter for selected pairing channels as functions of interaction and pairing strength within a mean-field treatment. Our results provide a systematic framework for analyzing superconductivity in the orbital Hatsugai-Kohmoto model and extend symmetry-based approaches to correlated multiband settings.

arXiv:2605.13259•May 13, 2026•0 citations

Decomposing momentum scales in the Hubbard Model: From Hatsugai-Kohmoto to Aubry-André

Dmitry Manning-Coe, Barry Bradlyn

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é 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é systems: the Aubry-André-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é systems in two dimensions. Code for reproducing all numerical results is available at https://github.com/chainik1125/decomposing-hubbard.

arXiv:2604.06588•Apr 8, 2026•0 citations

The SU(N) Holstein Model

Chunhan Feng, Linh Pham, G. G. Batrouni +1 more

From the condensed matter physics perspective, the most natural single orbital tight-binding Hamiltonians, and hence the most widely studied, contain two fermionic species, corresponding to spin up and spin down electrons. In cold atom systems, however, SU(N) symmetry, in which $N > 2$ fermionic species reside within a single band, also occurs. In order to understand such experiments, the SU(N) Hubbard model has been increasingly studied. Here we present determinant Quantum Monte Carlo simulations of the SU(N) {\it Holstein} Hamiltonian, in which $N$ fermionic species couple to a single local phonon mode. We show that at half filling it has an insulating charge density wave phase (CDW) at low temperatures, in which empty sites alternate with sites with $N$ particles. We determine the $N=3$ CDW phase diagram in the temperature, $T$, versus electron-phonon coupling, $α$, plane at fixed phonon frequency $ω_0$ and half-filling $ρ=1.5$. The critical temperature $T_c$ for $N=3$ can be as high as twice the maximum attainable for $N=2$. We also obtain the $N$ dependence of $T_c$ for a representative, fixed, $ω_0$ and $α$.

arXiv:2606.25240•Jun 23, 2026•0 citations

Real-space analysis of Hatsugai-Kohmoto interaction

Jan Skolimowski

The Hatsugai-Kohmoto interaction model has gained a lot of attention in recent years, due to the fact it is exactly solvable in momentum space in any dimension while capturing some key features of the Mott phase. Here a one-dimensional lattice model with this interaction is approached from the real-space perspective, to explore how breaking the translation invariance of a lattice affects the intuition built by studying the exact solution in $k$ space. The ground state properties of chains with periodic and open boundary conditions are calculated and compared with both the exact solution in momentum space as well as with analogous solutions of the Hubbard model. The results show that introducing hard edges enhances the ferromagnetic correlations and the system undergoes a magnetic transition before reaching the strong coupling limit. Understanding the impact of hard edges is a crucial step toward answering the looming question of the existence of edge states and other topological phenomena in systems with this type of interaction.

arXiv:2402.08557•Physical review. B./Physical review. B•Apr 16, 2024•9 citations

Dynamical Spectral Weight Transfer in the Orbital HK Model

Gaurav Tenkila, Jinchao Zhao, Philip Phillips

We compute explicitly the low-energy spectral weight (LESW) in the lower Hubbard band using the orbital Hatsugai-Kohmoto (OHK) model. We show that the dynamical mixing enters the LESW with a universal slope independent of the number of orbitals. This independence on orbital number corroborates the rapid convergence of the OHK model to the Hubbard model. As a result, OHK is an effective simulator of Hubbard physics.

arXiv:2406.17846•Jun 25, 2024•0 citations
We show that the physics of the SU(NNN) Hubbard model can be realistically simulated with the recently developed orbital Hatsugai-Kohmoto model. In this approach, the momentum mixing absent from the band Hatsugai-Kohmoto model is included by grouping nnn-Hubbard atoms into a cluster. We take advantage of the rapid convergence of this scheme (1/n21/n^21/n2) and show that for n≈O(10)n\approx O(10)n≈O(10), quantitative agreement with the determinantal quantum Monte Carlo method arises for the double occupancy as well as qualitative agreement for the filling and compressibility across the Mott transition. Additional features of this work are that we can obtain low-temperature physics, which could be qualitatively different from its high-temperature counterpart, and dynamical quantities without resorting to analytical continuation, thereby establishing the orbital Hatsugai-Kohmoto model as a useful alternative perspective for studying strong correlations in SU(NNN) systems.
The shift in the spectral weight distribution under doping is termed the spectral weight transfer (SWT) and is a key indicator of strongly correlated physics. The static SWT (t=0t=0t=0) for doping fraction xxx scales as xxx for semiconductors and as 2x2x2x for Mott insulators. The SWT receives corrections in t/Ut/Ut/U due to the noncommutativity of the kinetic and interaction terms in the Hamiltonian, termed the dynamic SWT. We compute explicitly the low-energy spectral weight (LESW) in the lower Hubbard band using the orbital Hatsugai-Kohmoto (OHK) model. We show that the dynamical mixing enters the LESW with a universal slope independent of the number of orbitals. This independence on orbital number corroborates the rapid convergence of the OHK model to the Hubbard model. As a result, OHK is an effective simulator of Hubbard physics.
Multiband superconductivity gives rise to a rich landscape of possible pairing states. Here we study superconductivity emerging from a half-filled normal state in a two-orbital extension of the Hatsugai-Kohmoto model, an exactly solvable model of correlated electrons with momentum-local interactions, which provides a minimal framework to explore the interplay of strong correlations, orbital structure and pairing symmetry. For a system with point-group symmetry D4h\rm D_{4h}D4h​, we classify the symmetry-allowed superconducting gap structures, taking into account spin, orbital and momentum degrees of freedom. We further compute the critical temperature and the superconducting order parameter for selected pairing channels as functions of interaction and pairing strength within a mean-field treatment. Our results provide a systematic framework for analyzing superconductivity in the orbital Hatsugai-Kohmoto model and extend symmetry-based approaches to correlated multiband settings.
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é 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é systems: the Aubry-André-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é systems in two dimensions. Code for reproducing all numerical results is available at https://github.com/chainik1125/decomposing-hubbard.
From the condensed matter physics perspective, the most natural single orbital tight-binding Hamiltonians, and hence the most widely studied, contain two fermionic species, corresponding to spin up and spin down electrons. In cold atom systems, however, SU(N) symmetry, in which N>2N > 2N>2 fermionic species reside within a single band, also occurs. In order to understand such experiments, the SU(N) Hubbard model has been increasingly studied. Here we present determinant Quantum Monte Carlo simulations of the SU(N) {\it Holstein} Hamiltonian, in which NNN fermionic species couple to a single local phonon mode. We show that at half filling it has an insulating charge density wave phase (CDW) at low temperatures, in which empty sites alternate with sites with NNN particles. We determine the N=3N=3N=3 CDW phase diagram in the temperature, TTT, versus electron-phonon coupling, ααα, plane at fixed phonon frequency ω0ω_0ω0​ and half-filling ρ=1.5ρ=1.5ρ=1.5. The critical temperature TcT_cTc​ for N=3N=3N=3 can be as high as twice the maximum attainable for N=2N=2N=2. We also obtain the NNN dependence of TcT_cTc​ for a representative, fixed, ω0ω_0ω0​ and ααα.
The Hatsugai-Kohmoto interaction model has gained a lot of attention in recent years, due to the fact it is exactly solvable in momentum space in any dimension while capturing some key features of the Mott phase. Here a one-dimensional lattice model with this interaction is approached from the real-space perspective, to explore how breaking the translation invariance of a lattice affects the intuition built by studying the exact solution in kkk space. The ground state properties of chains with periodic and open boundary conditions are calculated and compared with both the exact solution in momentum space as well as with analogous solutions of the Hubbard model. The results show that introducing hard edges enhances the ferromagnetic correlations and the system undergoes a magnetic transition before reaching the strong coupling limit. Understanding the impact of hard edges is a crucial step toward answering the looming question of the existence of edge states and other topological phenomena in systems with this type of interaction.