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:2604.09180

Eigenstate entanglement entropy in Bose-Hubbard models

G. Medoš, L. Vidmar

cond-mat.stat-mechcond-mat.quant-gascond-mat.str-elquant-phSubmitted Apr 10, 2026• Announced Apr 13, 2026• 0 citations
arXivHTMLPDF

Abstract

While the eigenstate entanglement entropy has been extensively studied for fermionic systems, much less is known about bosonic systems. Here, we study the entanglement entropy of mid-spectrum eigenstates of Bose-Hubbard models, focusing on weakly disordered models with and without particle-number conservation, and contrasting them with the translationally-invariant model. We analyze the volume-law and O(1) contributions to the entanglement entropy via the averages over mid-spectrum eigenstates and the corresponding distributions. We derive the volume-law coefficient of the entanglement entropy by generalizing the mean-field approach from [Phys. Rev. Lett. 119, 220603 (2017)] to many-body systems with a tunable local bosonic cutoff, which agrees with previous analytical and numerical results from [Phys. Rev. B 110, 235154 (2024)]. We show that the volume-law contribution to the entanglement entropy does not change upon breaking translational invariance via on-site disorder. We then numerically study the role of the subleading O(1) contribution to the entanglement entropy. We find that, in the particle-number conserving case, it exhibits a nontrivial dependence on the particle-number density and the local bosonic cutoff, while without particle-number conservation, results suggest the emergence of a universal O(1) contribution beyond the random pure state predictions.

Related papers

Volume-Law Entanglement Entropy of Typical Pure Quantum States

Eugenio Bianchi A, Hackl, L, Kieburg, M +2 more

The entanglement entropy of subsystems of typical eigenstates of quantum many-body Hamiltonians has recently been conjectured to be a diagnostic of quantum chaos and integrability. In quantum chaotic systems it has been found to behave as in typical pure states, while in integrable systems it has been found to behave as in typical pure Gaussian states. In this tutorial, we provide a pedagogical introduction to known results about the entanglement entropy of subsystems of typical pure states and of typical pure Gaussian states. They both exhibit a leading term that scales with the volume of the subsystem, when smaller than one half of the volume of the system, but the prefactor of the volume law is fundamentally different. It is constant (and maximal) for typical pure states, and it depends on the ratio between the volume of the subsystem and of the entire system for typical pure Gaussian states. Since particle-number conservation plays an important role in many physical Hamiltonians, we discuss its effect on typical pure states and on typical pure Gaussian states. We prove that, while the behavior of the leading volume-law terms does not change qualitatively, the nature of the subleading terms can change. In particular, subleading corrections can appear that depend on the square root of the volume of the subsystem. We unveil the origin of those corrections. Finally, we discuss the connection between the entanglement entropy of typical pure states and analytical results obtained in the context of random matrix theory, as well as numerical results obtained for physical Hamiltonians.

•Jul 1, 2022•123 citations

Entanglement in bipartite systems with symmetry: coupled chaotic kicked Bose-Hubbard systems

J. Himmelsbach, M. F. I. Kieler, A. Bäcker

We investigate the average eigenstate entanglement in a bipartite many-body system which exhibits a breaking of two local conservation laws into a global conserved quantity. Such a setting is realized by the particle number conservation in Bose-Hubbard systems. We devise a corresponding random matrix model which captures the universal features of this symmetry breaking and allows for applying powerful random matrix methods. By combining the concept of symmetry resolved entanglement with perturbation theory for quantum chaotic systems we obtain a universal entanglement transition depending on a single fundamental parameter. Furthermore, the symmetry resolved entanglement allows for separating the genuine entanglement from the part which originates from the conserved quantity. This latter contribution is quantified by the number entropy. For this it is shown that the symmetry breaking generates a localization of the eigenstates due to a banded structure of the time-evolution operator. By extrapolating the results beyond the perturbative regime we obtain an analytic description of the full transition.

quant-pharXiv:2607.29165•Aug 3, 2026•0 citations

Volume Law and Quantum Criticality in the Entanglement Entropy of Excited Eigenstates of the Quantum Ising Model

Lev Vidmar, Lucas Hackl, Eugenio Bianchi A +1 more

Much has been learned about universal properties of entanglement entropies in ground states of quantum many-body lattice systems. Here we unveil universal properties of the average bipartite entanglement entropy of eigenstates of the paradigmatic quantum Ising model in one dimension. The leading term exhibits a volume-law scaling that we argue is universal for translationally invariant quadratic models. The subleading term is constant at the critical field for the quantum phase transition and vanishes otherwise (in the thermodynamic limit); i.e., the critical field can be identified from subleading corrections to the average (over all eigenstates) entanglement entropy.

arXiv:1808.08963•Physical Review Letters•Nov 28, 2018•92 citations

Entanglement Entropy of Typical Eigenstates of Translationally Invariant Quadratic Hamiltonians

Lev Vidmar, Lucas Hackl, Eugenio Bianchi A +1 more

In a seminal paper [Phys. Rev. Lett. 71, 1291 (1993)], Page conjectured that the entanglement entropy of typical pure states is $S_{\rm typ}\simeq\ln{\cal D}_{\rm A} - (1/2) {\cal D}_{\rm A}^2/{\cal D}$, where ${\cal D}_{\rm A}$ and ${\cal D}$ are the Hilbert space dimensions of the subsystem and the system, respectively, and typicality is with respect to the Haar measure. Typical pure states are hence (nearly) maximally entangled. We prove that $S^+_1=\ln{\cal D}_{\rm A} - [1/(2\ln2)] (\ln{\cal D}_{\rm A})^2/\ln{\cal D}$ is an upper bound for the average entanglement entropy of all many-body fermionic eigenstates of a class of translationally invariant quadratic Hamiltonians in one dimension, and that the variance vanishes with increasing system size. Consequently, whenever the size of the subsystem is not a vanishing fraction of the size of the system, typical eigenstates of such Hamiltonians have a volume-law entanglement but do not become maximally entangled as ${\cal D}\rightarrow\infty$, i.e., they are qualitatively different from typical pure states. We also prove that, in the limit in which the subsystem size is a vanishing fraction of the system size, the average entanglement entropy is maximal.

•Mar 8, 2017•0 citations

Average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians

Maksymilian Kliczkowski, Rafał Świętek, Lev Vidmar +1 more

To which degree the average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians agrees with that of random pure states is a question that has attracted considerable attention in the recent years. While there is substantial evidence that the leading (volume-law) terms are identical, which and how subleading terms differ between them is less clear. Here we carry out state-of-the-art full exact diagonalization calculations of clean spin-1/2 XYZ and XXZ chains with integrability breaking terms to address this question in the absence and presence of U(1) symmetry, respectively. We first introduce the notion of maximally chaotic regime, for the chain sizes amenable to full exact diagonalization calculations, as the regime in Hamiltonian parameters in which the level spacing ratio, the distribution of eigenstate coefficients, and the entanglement entropy are closest to the random matrix theory predictions. In this regime, we carry out a finite-size scaling analysis of the subleading terms of the average entanglement entropy of midspectrum eigenstates when different fractions ν of the spectrum are included in the average. We find indications that, for ν→0, the magnitude of the negative O(1) correction is only slightly greater than the one predicted for random pure states. For finite ν, following a phenomenological approach, we derive a simple expression that describes the numerically observed ν dependence of the O(1) deviation from the prediction for random pure states.

arXiv:2303.13577•Physical review. E•Jun 20, 2023•22 citations

Entanglement Entropy of Eigenstates of Quantum Chaotic Hamiltonians

Lev Vidmar, Marcos Rigol

In quantum statistical mechanics, it is of fundamental interest to understand how close the bipartite entanglement entropy of eigenstates of quantum chaotic Hamiltonians is to maximal. For random pure states in the Hilbert space, the average entanglement entropy is known to be nearly maximal, with a deviation that is, at most, a constant. Here we prove that, in a system that is away from half filling and divided in two equal halves, an upper bound for the average entanglement entropy of random pure states with a fixed particle number and normally distributed real coefficients exhibits a deviation from the maximal value that grows with the square root of the volume of the system. Exact numerical results for highly excited eigenstates of a particle number conserving quantum chaotic model indicate that the bound is saturated with increasing system size.

arXiv:1708.08453•Physical Review Letters•Nov 29, 2017•175 citations
While the eigenstate entanglement entropy has been extensively studied for fermionic systems, much less is known about bosonic systems. Here, we study the entanglement entropy of mid-spectrum eigenstates of Bose-Hubbard models, focusing on weakly disordered models with and without particle-number conservation, and contrasting them with the translationally-invariant model. We analyze the volume-law and O(1) contributions to the entanglement entropy via the averages over mid-spectrum eigenstates and the corresponding distributions. We derive the volume-law coefficient of the entanglement entropy by generalizing the mean-field approach from [Phys. Rev. Lett. 119, 220603 (2017)] to many-body systems with a tunable local bosonic cutoff, which agrees with previous analytical and numerical results from [Phys. Rev. B 110, 235154 (2024)]. We show that the volume-law contribution to the entanglement entropy does not change upon breaking translational invariance via on-site disorder. We then numerically study the role of the subleading O(1) contribution to the entanglement entropy. We find that, in the particle-number conserving case, it exhibits a nontrivial dependence on the particle-number density and the local bosonic cutoff, while without particle-number conservation, results suggest the emergence of a universal O(1) contribution beyond the random pure state predictions.
The entanglement entropy of subsystems of typical eigenstates of quantum many-body Hamiltonians has recently been conjectured to be a diagnostic of quantum chaos and integrability. In quantum chaotic systems it has been found to behave as in typical pure states, while in integrable systems it has been found to behave as in typical pure Gaussian states. In this tutorial, we provide a pedagogical introduction to known results about the entanglement entropy of subsystems of typical pure states and of typical pure Gaussian states. They both exhibit a leading term that scales with the volume of the subsystem, when smaller than one half of the volume of the system, but the prefactor of the volume law is fundamentally different. It is constant (and maximal) for typical pure states, and it depends on the ratio between the volume of the subsystem and of the entire system for typical pure Gaussian states. Since particle-number conservation plays an important role in many physical Hamiltonians, we discuss its effect on typical pure states and on typical pure Gaussian states. We prove that, while the behavior of the leading volume-law terms does not change qualitatively, the nature of the subleading terms can change. In particular, subleading corrections can appear that depend on the square root of the volume of the subsystem. We unveil the origin of those corrections. Finally, we discuss the connection between the entanglement entropy of typical pure states and analytical results obtained in the context of random matrix theory, as well as numerical results obtained for physical Hamiltonians.
We investigate the average eigenstate entanglement in a bipartite many-body system which exhibits a breaking of two local conservation laws into a global conserved quantity. Such a setting is realized by the particle number conservation in Bose-Hubbard systems. We devise a corresponding random matrix model which captures the universal features of this symmetry breaking and allows for applying powerful random matrix methods. By combining the concept of symmetry resolved entanglement with perturbation theory for quantum chaotic systems we obtain a universal entanglement transition depending on a single fundamental parameter. Furthermore, the symmetry resolved entanglement allows for separating the genuine entanglement from the part which originates from the conserved quantity. This latter contribution is quantified by the number entropy. For this it is shown that the symmetry breaking generates a localization of the eigenstates due to a banded structure of the time-evolution operator. By extrapolating the results beyond the perturbative regime we obtain an analytic description of the full transition.
Much has been learned about universal properties of entanglement entropies in ground states of quantum many-body lattice systems. Here we unveil universal properties of the average bipartite entanglement entropy of eigenstates of the paradigmatic quantum Ising model in one dimension. The leading term exhibits a volume-law scaling that we argue is universal for translationally invariant quadratic models. The subleading term is constant at the critical field for the quantum phase transition and vanishes otherwise (in the thermodynamic limit); i.e., the critical field can be identified from subleading corrections to the average (over all eigenstates) entanglement entropy.
In a seminal paper [Phys. Rev. Lett. 71, 1291 (1993)], Page conjectured that the entanglement entropy of typical pure states is Styp≃ln⁡DA−(1/2)DA2/DS_{\rm typ}\simeq\ln{\cal D}_{\rm A} - (1/2) {\cal D}_{\rm A}^2/{\cal D}Styp​≃lnDA​−(1/2)DA2​/D, where DA{\cal D}_{\rm A}DA​ and D{\cal D}D are the Hilbert space dimensions of the subsystem and the system, respectively, and typicality is with respect to the Haar measure. Typical pure states are hence (nearly) maximally entangled. We prove that S1+=ln⁡DA−[1/(2ln⁡2)](ln⁡DA)2/ln⁡DS^+_1=\ln{\cal D}_{\rm A} - [1/(2\ln2)] (\ln{\cal D}_{\rm A})^2/\ln{\cal D}S1+​=lnDA​−[1/(2ln2)](lnDA​)2/lnD is an upper bound for the average entanglement entropy of all many-body fermionic eigenstates of a class of translationally invariant quadratic Hamiltonians in one dimension, and that the variance vanishes with increasing system size. Consequently, whenever the size of the subsystem is not a vanishing fraction of the size of the system, typical eigenstates of such Hamiltonians have a volume-law entanglement but do not become maximally entangled as D→∞{\cal D}\rightarrow\inftyD→∞, i.e., they are qualitatively different from typical pure states. We also prove that, in the limit in which the subsystem size is a vanishing fraction of the system size, the average entanglement entropy is maximal.
To which degree the average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians agrees with that of random pure states is a question that has attracted considerable attention in the recent years. While there is substantial evidence that the leading (volume-law) terms are identical, which and how subleading terms differ between them is less clear. Here we carry out state-of-the-art full exact diagonalization calculations of clean spin-1/2 XYZ and XXZ chains with integrability breaking terms to address this question in the absence and presence of U(1) symmetry, respectively. We first introduce the notion of maximally chaotic regime, for the chain sizes amenable to full exact diagonalization calculations, as the regime in Hamiltonian parameters in which the level spacing ratio, the distribution of eigenstate coefficients, and the entanglement entropy are closest to the random matrix theory predictions. In this regime, we carry out a finite-size scaling analysis of the subleading terms of the average entanglement entropy of midspectrum eigenstates when different fractions ν of the spectrum are included in the average. We find indications that, for ν→0, the magnitude of the negative O(1) correction is only slightly greater than the one predicted for random pure states. For finite ν, following a phenomenological approach, we derive a simple expression that describes the numerically observed ν dependence of the O(1) deviation from the prediction for random pure states.
In quantum statistical mechanics, it is of fundamental interest to understand how close the bipartite entanglement entropy of eigenstates of quantum chaotic Hamiltonians is to maximal. For random pure states in the Hilbert space, the average entanglement entropy is known to be nearly maximal, with a deviation that is, at most, a constant. Here we prove that, in a system that is away from half filling and divided in two equal halves, an upper bound for the average entanglement entropy of random pure states with a fixed particle number and normally distributed real coefficients exhibits a deviation from the maximal value that grows with the square root of the volume of the system. Exact numerical results for highly excited eigenstates of a particle number conserving quantum chaotic model indicate that the bound is saturated with increasing system size.