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

Physical Review Research

Hierarchical algebraic structure and recursive solution of the Bose-Hubbard model

Anonymous, 戴连蓉, J. P. Draayer

Published Sep 1, 2026• 0 citations
Publisher

Abstract

We uncover a previously unrecognized hierarchical algebraic structure in the Bose-Hubbard model that enables a systematic construction of its spectrum without relying on integrability. We develop a recursive framework in which, within each fixed particle-number sector, the eigenproblem of a p -site system is reduced to a closed set of algebraic relations determined entirely by the eigenpairs and single-particle transfer amplitudes of the ( p − 1 ) -site subsystem. This reduction yields two complementary sets of Bethe-type equations whose solutions exhaust the spectrum. The underlying structure originates from the locality of the Hamiltonian, which constrains the coupling between subsystems and gives rise to an emergent sparse, near-tridiagonal representation of the many-body problem. As an explicit demonstration, we construct the complete spectrum of the three-site Bose-Hubbard model from that of the two-site system for particle numbers n ≤ 3 , obtaining exact agreement with direct diagonalization. The resulting formulation substantially reduces computational cost and reveals a form of hierarchical solvability in a nonintegrable quantum many-body system. More broadly, our results suggest that recursive algebraic constructions may provide a general framework for analyzing lattice models beyond the conventional integrable paradigm.

DOI
10.1103/psg9-5rv3

Related papers

Exact Solution of Two-Site Bose–Hubbard Model with Generic Open Boundaries*

Zhirong Xin, Tao Yang, Kun Hao +1 more

Abstract The Bose–Hubbard model is a paradigm for the study of strongly correlated bosonic systems. We study the two-site Bose–Hubbard model with generic integrable open boundaries specified by the most general non-diagonal reflecting matrices. Besides the inhomogeneous parameters, the model itself has three free boundary parameters, which break the U(1)-symmetry, in other words, break the particle number conservation. The Hamiltonian H under these circumstances is constructed. With the help of the off-diagonal Bethe Ansatz method, we successfully obtain the corresponding Bethe Ansatz equations as well as the eigenvalues.

Communications in Theoretical Physics•Dec 1, 2015•1 citation

Scar subspaces stabilized by algebraic closure: Beyond equally-spaced spectra and exact solvability

Chihiro Matsui

We propose algebraic closure as a stabilizing mechanism for quantum many-body scar subspaces, showing that scar states need not rely on an embedded equally spaced spectrum or analytical solvability. We formulate Lie-algebraically closed invariant subspaces and construct local quantum many-body Hamiltonians that preserve an $\mathfrak{su}(N)$-invariant scar subspace but need not act diagonally within it. Such non-diagonal action mixes the conventional tower states while preserving the multidirectional lattice-like spectral structure, allowing nonthermal dynamics and persistent multifrequency oscillations to survive even when individual scar eigenstates become analytically intractable.

cond-mat.stat-mechCross-listedarXiv:2604.11015•Apr 14, 2026•0 citations

Exact Solution of the Bose-Hubbard Model with Unidirectional Hopping

Mingchen Zheng, Yi Qiao, Yupeng Wang +2 more

A one-dimensional Bose-Hubbard model with unidirectional hopping is shown to be exactly solvable. Applying the algebraic Bethe ansatz method, we prove the integrability of the model and derive the Bethe ansatz equations. The exact eigenvalue spectrum can be obtained by solving these equations. The distribution of Bethe roots reveals the presence of a superfluid-Mott insulator transition at the ground state, and the critical point is determined. By adjusting the boundary parameter, we demonstrate the existence of a non-Hermitian skin effect even in the presence of interaction, but it is completely suppressed for the Mott insulator state in the thermodynamical limit. Our result represents a new class of exactly solvable non-Hermitian many-body systems, which has no Hermitian correspondence and can be used as a benchmark for various numerical techniques developed for non-Hermitian many-body systems.

arXiv:2305.00439•Physical Review Letters•Feb 23, 2024•40 citations

Resolvent-Based Self-Consistent Framework with Hierarchical Correlation Expansion for Strongly Correlated Many-Body Systems

Zhi-qiang Huang, Qing-yu Cai

We develop a nonperturbative framework for generic nonintegrable many-body systems that reorganizes the expansion of diagonal Green's functions. Starting from exact projection identities and the spectral representation of the resolvent, we derive a recursive hierarchy for the self-energy in which cross-correlated propagation processes are systematically rewritten in terms of diagonal resolvents. Under a diagonal closure approximation, the hierarchy becomes formally closed yet remains systematically improvable. The framework combines two nonperturbative mechanisms. First, a Lanczos continued-fraction representation provides a recursive single-resolvent structure that naturally produces non-Lorentzian spectral features beyond self-consistent Born approximations. Second, an exact projected multi-resolvent hierarchy introduces nonlocal frequency couplings through products of resolvents and their Hilbert transforms. These contributions mix parity sectors under energy reflection and generate spectral skewness, which is absent in single-resolvent closures. To solve the resulting equations, we employ a hierarchy of Lorentzian, Gaussian, and Voigt-type ansätze together with an effective Faddeeva self-energy representation ensuring analyticity and causality. Spectral broadening, distribution tails, and higher-order fluctuations emerge from the interplay between continued-fraction recursion and multi-resolvent correlations. The framework requires no small expansion parameters or diagrammatic truncations, relying instead on ETH-type statistical assumptions appropriate for dense chaotic spectra. It provides a unified route from microscopic interactions to emergent spectral structure, revealing a progression from single-pole self-consistent dynamics to continued-fraction renormalization and finally to multi-resolvent interference effects.

quant-pharXiv:2604.00606•Apr 2, 2026•0 citations

Bethe states for the two-site Bose–Hubbard model: A binomial approach

Gilberto N. Santos Filho, Changrim Ahn, Angela Foerster +1 more

We calculate explicitly the Bethe vectors states by the algebraic Bethe ansatz method with the gl(2)-invariant R-matrix for the two-site Bose–Hubbard model. Using a binomial expansion of the n-th power of a sum of two operators we get and solve a recursion equation. We calculate the scalar product and the norm of the Bethe vectors states. The form factors of the imbalance current operator are also computed.

arXiv:1503.07885•Physics Letters B•May 5, 2015•5 citations

Exact solution of the Bose-Hubbard model on the Bethe lattice

Guilhem Semerjian, Marco Tarzia, Francesco Zamponi

The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.

arXiv:0904.3075•Physical Review B•Jul 31, 2009•42 citations
We uncover a previously unrecognized hierarchical algebraic structure in the Bose-Hubbard model that enables a systematic construction of its spectrum without relying on integrability. We develop a recursive framework in which, within each fixed particle-number sector, the eigenproblem of a p -site system is reduced to a closed set of algebraic relations determined entirely by the eigenpairs and single-particle transfer amplitudes of the ( p − 1 ) -site subsystem. This reduction yields two complementary sets of Bethe-type equations whose solutions exhaust the spectrum. The underlying structure originates from the locality of the Hamiltonian, which constrains the coupling between subsystems and gives rise to an emergent sparse, near-tridiagonal representation of the many-body problem. As an explicit demonstration, we construct the complete spectrum of the three-site Bose-Hubbard model from that of the two-site system for particle numbers n ≤ 3 , obtaining exact agreement with direct diagonalization. The resulting formulation substantially reduces computational cost and reveals a form of hierarchical solvability in a nonintegrable quantum many-body system. More broadly, our results suggest that recursive algebraic constructions may provide a general framework for analyzing lattice models beyond the conventional integrable paradigm.
Abstract The Bose–Hubbard model is a paradigm for the study of strongly correlated bosonic systems. We study the two-site Bose–Hubbard model with generic integrable open boundaries specified by the most general non-diagonal reflecting matrices. Besides the inhomogeneous parameters, the model itself has three free boundary parameters, which break the U(1)-symmetry, in other words, break the particle number conservation. The Hamiltonian H under these circumstances is constructed. With the help of the off-diagonal Bethe Ansatz method, we successfully obtain the corresponding Bethe Ansatz equations as well as the eigenvalues.
We propose algebraic closure as a stabilizing mechanism for quantum many-body scar subspaces, showing that scar states need not rely on an embedded equally spaced spectrum or analytical solvability. We formulate Lie-algebraically closed invariant subspaces and construct local quantum many-body Hamiltonians that preserve an su(N)\mathfrak{su}(N)su(N)-invariant scar subspace but need not act diagonally within it. Such non-diagonal action mixes the conventional tower states while preserving the multidirectional lattice-like spectral structure, allowing nonthermal dynamics and persistent multifrequency oscillations to survive even when individual scar eigenstates become analytically intractable.
A one-dimensional Bose-Hubbard model with unidirectional hopping is shown to be exactly solvable. Applying the algebraic Bethe ansatz method, we prove the integrability of the model and derive the Bethe ansatz equations. The exact eigenvalue spectrum can be obtained by solving these equations. The distribution of Bethe roots reveals the presence of a superfluid-Mott insulator transition at the ground state, and the critical point is determined. By adjusting the boundary parameter, we demonstrate the existence of a non-Hermitian skin effect even in the presence of interaction, but it is completely suppressed for the Mott insulator state in the thermodynamical limit. Our result represents a new class of exactly solvable non-Hermitian many-body systems, which has no Hermitian correspondence and can be used as a benchmark for various numerical techniques developed for non-Hermitian many-body systems.
We develop a nonperturbative framework for generic nonintegrable many-body systems that reorganizes the expansion of diagonal Green's functions. Starting from exact projection identities and the spectral representation of the resolvent, we derive a recursive hierarchy for the self-energy in which cross-correlated propagation processes are systematically rewritten in terms of diagonal resolvents. Under a diagonal closure approximation, the hierarchy becomes formally closed yet remains systematically improvable. The framework combines two nonperturbative mechanisms. First, a Lanczos continued-fraction representation provides a recursive single-resolvent structure that naturally produces non-Lorentzian spectral features beyond self-consistent Born approximations. Second, an exact projected multi-resolvent hierarchy introduces nonlocal frequency couplings through products of resolvents and their Hilbert transforms. These contributions mix parity sectors under energy reflection and generate spectral skewness, which is absent in single-resolvent closures. To solve the resulting equations, we employ a hierarchy of Lorentzian, Gaussian, and Voigt-type ansätze together with an effective Faddeeva self-energy representation ensuring analyticity and causality. Spectral broadening, distribution tails, and higher-order fluctuations emerge from the interplay between continued-fraction recursion and multi-resolvent correlations. The framework requires no small expansion parameters or diagrammatic truncations, relying instead on ETH-type statistical assumptions appropriate for dense chaotic spectra. It provides a unified route from microscopic interactions to emergent spectral structure, revealing a progression from single-pole self-consistent dynamics to continued-fraction renormalization and finally to multi-resolvent interference effects.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.