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:1703.02979 · Physical Review Letters

Entanglement Entropy of Eigenstates of Quadratic Fermionic Hamiltonians

Lev Vidmar, Lucas Hackl, Eugenio Bianchi A, Marcos Rigol

Published Jul 11, 2017• 134 citations
arXivHTMLPDF

Abstract

In a seminal paper [D. N. Page, Phys. Rev. Lett. 71, 1291 (1993)], Page proved that the average entanglement entropy of subsystems of random pure states is ${S}_{\text{ave}}\ensuremath{\simeq}\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}(1/2){\mathcal{D}}_{\mathrm{A}}^{2}/\mathcal{D}$ for $1\ensuremath{\ll}{\mathcal{D}}_{\mathrm{A}}\ensuremath{\le}\sqrt{\mathcal{D}}$, where ${\mathcal{D}}_{\mathrm{A}}$ and $\mathcal{D}$ are the Hilbert space dimensions of the subsystem and the system, respectively. Hence, typical pure states are (nearly) maximally entangled. We develop tools to compute the average entanglement entropy $⟨S⟩$ of all eigenstates of quadratic fermionic Hamiltonians. In particular, we derive exact bounds for the most general translationally invariant models $\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}(\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}{)}^{2}/\mathrm{ln}\mathcal{D}\ensuremath{\le}⟨S⟩\ensuremath{\le}\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}[1/(2\mathrm{ln}2)](\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}{)}^{2}/\mathrm{ln}\mathcal{D}$. Consequently, we prove that (i) if the subsystem size is a finite fraction of the system size, then $⟨S⟩<\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}$ in the thermodynamic limit; i.e., the average over eigenstates of the Hamiltonian departs from the result for typical pure states, and (ii) in the limit in which the subsystem size is a vanishing fraction of the system size, the average entanglement entropy is maximal; i.e., typical eigenstates of such Hamiltonians exhibit eigenstate thermalization.

DOI
10.1103/physrevlett.119.020601

Related papers

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

Page curve for fermionic Gaussian states

Eugenio Bianchi A, Lucas Hackl, Mario Kieburg

In a seminal paper, Page found the exact formula for the average entanglement entropy for a pure random state. We consider the analogous problem for the ensemble of pure fermionic Gaussian states, which plays a crucial role in the context of random free Hamiltonians. Using recent results from random matrix theory, we show that the average entanglement entropy of pure random fermionic Gaussian states in a subsystem of ${N}_{A}$ out of $N$ degrees of freedom is given by ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{\mathrm{G}}=(N\ensuremath{-}\frac{1}{2})\mathrm{\ensuremath{\Psi}}(2N)+(\frac{1}{4}\ensuremath{-}{N}_{A})\mathrm{\ensuremath{\Psi}}(N)+(\frac{1}{2}+{N}_{A}\ensuremath{-}N)\mathrm{\ensuremath{\Psi}}(2N\ensuremath{-}2{N}_{A})\ensuremath{-}\frac{1}{4}\mathrm{\ensuremath{\Psi}}(N\ensuremath{-}{N}_{A})\ensuremath{-}{N}_{A}$, where $\mathrm{\ensuremath{\Psi}}$ is the digamma function. Its asymptotic behavior in the thermodynamic limit is given by ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{\mathrm{G}}=N(log2\ensuremath{-}1)f+N(f\ensuremath{-}1)log(1\ensuremath{-}f)+\frac{1}{2}f+\frac{1}{4}log(1\ensuremath{-}f)\phantom{\rule{0.16em}{0ex}}+\phantom{\rule{0.16em}{0ex}}O(1/N)$, where $f={N}_{A}/N\ensuremath{\le}1/2$. Remarkably, its leading order agrees with the average over eigenstates of random quadratic Hamiltonians with number conservation, as found by \L{}yd\ifmmode \dot{z}\else \.{z}\fi{}ba, Rigol, and Vidmar. Finally, we compute the variance in the thermodynamic limit, given by the constant ${lim}_{N\ensuremath{\rightarrow}\ensuremath{\infty}}{(\mathrm{\ensuremath{\Delta}}{S}_{A})}_{\mathrm{G}}^{2}=\frac{1}{2}[f+{f}^{2}+log(1\ensuremath{-}f)]$.

arXiv:2103.05416•Physical review. B./Physical review. B•Jun 29, 2021•63 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

Average eigenstate entanglement entropy of the XY chain in a transverse field and its universality for translationally invariant quadratic fermionic models

Lucas Hackl, Lev Vidmar, Marcos Rigol +1 more

We recently showed [Phys. Rev. Lett. 121, 220602 (2018)] that the average bipartite entanglement entropy of all energy eigenstates of the quantum Ising chain exhibits a universal (for translationally invariant quadratic fermionic models) leading term that scales linearly with the subsystem's volume, while in the thermodynamic limit the first subleading correction does not vanish at the critical field (it only depends on the ratio $f$ between the volume of the subsystem and volume of the system) and vanishes otherwise. Here we show, analytically for bounds and numerically for averages, that the same remains true for the spin-1/2 XY chain in a transverse magnetic field. We then tighten the bounds for the coefficient of the universal volume-law term, which is a concave function of $f$. We develop a systematic approach to compute upper and lower bounds, and we provide explicit analytic expressions for up to the fourth-order bounds.

arXiv:1812.08757•Physical review. B./Physical review. B•Feb 11, 2019•64 citations

Eigenstate Entanglement Entropy in Random Quadratic Hamiltonians

Patrycja Łydżba, Marcos Rigol, Lev Vidmar

The eigenstate entanglement entropy is a powerful tool to distinguish integrable from generic quantum-chaotic models. In integrable models, the average eigenstate entanglement entropy (over all Hamiltonian eigenstates) has a volume-law coefficient that generally depends on the subsystem fraction. In contrast, it is maximal (subsystem fraction independent) in quantum-chaotic models. Using random matrix theory for quadratic Hamiltonians, we obtain a closed-form expression for the average eigenstate entanglement entropy as a function of the subsystem fraction. We test it against numerical results for the quadratic Sachdev-Ye-Kitaev model and show that it describes the results for the power-law random banded matrix model (in the delocalized regime). We show that localization in quasimomentum space produces (small) deviations from our analytic predictions.

arXiv:2006.11302•Physical Review Letters•Oct 27, 2020•82 citations

Eigenstate entanglement entropy in Bose-Hubbard models

G. Medoš, L. Vidmar

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.

cond-mat.stat-mechCross-listedarXiv:2604.09180•Apr 13, 2026•0 citations
In a seminal paper [D. N. Page, Phys. Rev. Lett. 71, 1291 (1993)], Page proved that the average entanglement entropy of subsystems of random pure states is Save\ensuremath≃lnDA\ensuremath−(1/2)DA2/D{S}_{\text{ave}}\ensuremath{\simeq}\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}(1/2){\mathcal{D}}_{\mathrm{A}}^{2}/\mathcal{D}Save​\ensuremath≃lnDA​\ensuremath−(1/2)DA2​/D for 1\ensuremath≪DA\ensuremath≤D1\ensuremath{\ll}{\mathcal{D}}_{\mathrm{A}}\ensuremath{\le}\sqrt{\mathcal{D}}1\ensuremath≪DA​\ensuremath≤D​, where DA{\mathcal{D}}_{\mathrm{A}}DA​ and D\mathcal{D}D are the Hilbert space dimensions of the subsystem and the system, respectively. Hence, typical pure states are (nearly) maximally entangled. We develop tools to compute the average entanglement entropy ⟨S⟩⟨S⟩⟨S⟩ of all eigenstates of quadratic fermionic Hamiltonians. In particular, we derive exact bounds for the most general translationally invariant models lnDA\ensuremath−(lnDA)2/lnD\ensuremath≤⟨S⟩\ensuremath≤lnDA\ensuremath−[1/(2ln2)](lnDA)2/lnD\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}(\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}{)}^{2}/\mathrm{ln}\mathcal{D}\ensuremath{\le}⟨S⟩\ensuremath{\le}\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}\ensuremath{-}[1/(2\mathrm{ln}2)](\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}{)}^{2}/\mathrm{ln}\mathcal{D}lnDA​\ensuremath−(lnDA​)2/lnD\ensuremath≤⟨S⟩\ensuremath≤lnDA​\ensuremath−[1/(2ln2)](lnDA​)2/lnD. Consequently, we prove that (i) if the subsystem size is a finite fraction of the system size, then ⟨S⟩<lnDA⟨S⟩<\mathrm{ln}{\mathcal{D}}_{\mathrm{A}}⟨S⟩<lnDA​ in the thermodynamic limit; i.e., the average over eigenstates of the Hamiltonian departs from the result for typical pure states, and (ii) in the limit in which the subsystem size is a vanishing fraction of the system size, the average entanglement entropy is maximal; i.e., typical eigenstates of such Hamiltonians exhibit eigenstate thermalization.
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.
In a seminal paper, Page found the exact formula for the average entanglement entropy for a pure random state. We consider the analogous problem for the ensemble of pure fermionic Gaussian states, which plays a crucial role in the context of random free Hamiltonians. Using recent results from random matrix theory, we show that the average entanglement entropy of pure random fermionic Gaussian states in a subsystem of NA{N}_{A}NA​ out of NNN degrees of freedom is given by ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{\mathrm{G}}=(N\ensuremath{-}\frac{1}{2})\mathrm{\ensuremath{\Psi}}(2N)+(\frac{1}{4}\ensuremath{-}{N}_{A})\mathrm{\ensuremath{\Psi}}(N)+(\frac{1}{2}+{N}_{A}\ensuremath{-}N)\mathrm{\ensuremath{\Psi}}(2N\ensuremath{-}2{N}_{A})\ensuremath{-}\frac{1}{4}\mathrm{\ensuremath{\Psi}}(N\ensuremath{-}{N}_{A})\ensuremath{-}{N}_{A}$, where $\mathrm{\ensuremath{\Psi}}$ is the digamma function. Its asymptotic behavior in the thermodynamic limit is given by ${\ensuremath{\langle}{S}_{A}\ensuremath{\rangle}}_{\mathrm{G}}=N(log2\ensuremath{-}1)f+N(f\ensuremath{-}1)log(1\ensuremath{-}f)+\frac{1}{2}f+\frac{1}{4}log(1\ensuremath{-}f)\phantom{\rule{0.16em}{0ex}}+\phantom{\rule{0.16em}{0ex}}O(1/N)$, where $f={N}_{A}/N\ensuremath{\le}1/2$. Remarkably, its leading order agrees with the average over eigenstates of random quadratic Hamiltonians with number conservation, as found by \L{}yd\ifmmode \dot{z}\else \.{z}\fi{}ba, Rigol, and Vidmar. Finally, we compute the variance in the thermodynamic limit, given by the constant ${lim}_{N\ensuremath{\rightarrow}\ensuremath{\infty}}{(\mathrm{\ensuremath{\Delta}}{S}_{A})}_{\mathrm{G}}^{2}=\frac{1}{2}[f+{f}^{2}+log(1\ensuremath{-}f)]$.
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.
We recently showed [Phys. Rev. Lett. 121, 220602 (2018)] that the average bipartite entanglement entropy of all energy eigenstates of the quantum Ising chain exhibits a universal (for translationally invariant quadratic fermionic models) leading term that scales linearly with the subsystem's volume, while in the thermodynamic limit the first subleading correction does not vanish at the critical field (it only depends on the ratio fff between the volume of the subsystem and volume of the system) and vanishes otherwise. Here we show, analytically for bounds and numerically for averages, that the same remains true for the spin-1/2 XY chain in a transverse magnetic field. We then tighten the bounds for the coefficient of the universal volume-law term, which is a concave function of fff. We develop a systematic approach to compute upper and lower bounds, and we provide explicit analytic expressions for up to the fourth-order bounds.
The eigenstate entanglement entropy is a powerful tool to distinguish integrable from generic quantum-chaotic models. In integrable models, the average eigenstate entanglement entropy (over all Hamiltonian eigenstates) has a volume-law coefficient that generally depends on the subsystem fraction. In contrast, it is maximal (subsystem fraction independent) in quantum-chaotic models. Using random matrix theory for quadratic Hamiltonians, we obtain a closed-form expression for the average eigenstate entanglement entropy as a function of the subsystem fraction. We test it against numerical results for the quadratic Sachdev-Ye-Kitaev model and show that it describes the results for the power-law random banded matrix model (in the delocalized regime). We show that localization in quasimomentum space produces (small) deviations from our analytic predictions.
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.