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

Odd/Even or Half ? Entanglement Anomaly in the Bose-Hubbard model

Humberto E. López-Hernández, Santiago F. Caballero-Benitez

cond-mat.quant-gasphysics.atom-phquant-phSubmitted Sep 17, 2026• Announced Sep 18, 2026• 0 citations
arXivHTMLPDF

Abstract

The area law relates the bipartite entanglement entropy of a quantum many-body ground state to the size of the boundary between the subsystems, but the geometry of this boundary is rarely discussed. We inspect this in the 1D Bose-Hubbard model at fixed density by comparing four spatial bipartitions of the periodic lattice: first half, second half, even sites, and odd sites; sharing the same number of sites but differing in how the boundary is arranged. We find analytical limits with perturbation theory: in the Mott insulator the contiguous cut obeys the area lay while the alternating cut obeys a volume law $S\propto N_s$, in this sense an anomaly, and for the superfluid both cuts colapse to the binomial saturation due to delocalization of the state. We formulate these limits as a statement about the many-body problem using a generalized slave-boson approach based on mean-field with quantum fluctuations while verifying with Exact Diagonalization (ED) for small lattice sizes and Densitiy Matrix Renormalization Group (DMRG) simulations for $N_s\gg 1$. The slave-boson Gaussian ground state allows to compute the entanglement entropy from a reduced correlation matrix for any desired bipartition consistent with ED and DMRG results. Using slave bosons the computational cost is set by the local cutoff $n_{\max}$ rather than the Hilbert space dimension, so we can reach lattice sizes far beyond ED. Our method is capable of establishing the partition-dependent scaling laws as a many-body feature, not only a finite-size effect, in great agreement with the ED for $N_s\in[4,10]$ and DMRG for larger lattice sizes.

Comments
7 pages, 2 tables, 3 figures and supplemental material with 9 pages, 3 tables and 4 Figures

Related papers

Entanglement Entropy across the Superfluid-Insulator Transition: A Signature of Bosonic Criticality

Irénée Frérot, Tommaso Roscilde

We study the entanglement entropy and entanglement spectrum of the paradigmatic Bose-Hubbard model, describing strongly correlated bosons on a lattice. The use of a controlled approximation-the slave-boson approach-allows us to study entanglement in all regimes of the model (and, most importantly, across its superfluid-Mott-insulator transition) at a minimal cost. We find that the area-law scaling of entanglement-verified in all the phases-exhibits a sharp singularity at the transition. The singularity is greatly enhanced when the transition is crossed at fixed, integer filling, due to a richer entanglement spectrum containing an additional gapless mode, which descends from the amplitude (Higgs) mode of the global excitation spectrum-while this mode remains gapped at the generic (commensurate-incommensurate) transition with variable filling. Hence, the entanglement properties contain a unique signature of the two different forms of bosonic criticality exhibited by the Bose-Hubbard model.

arXiv:1512.00805•Physical Review Letters•May 11, 2016•70 citations

Entanglement structure of the Hubbard model in momentum space

G. Ehlers, Jenő Sólyom, Örs Legeza +1 more

We study the properties of the ground states of the one- and two-dimensional Hubbard models at half-filling and moderate doping using entanglement-based measures, which we calculate numerically using the momentum-space density matrix renormalization group (DMRG). In particular, we investigate quantities such as the single-site entropy and two-site mutual information of single-particle momentum states as well as the behavior of the bipartite subsystem entropy for partitions in momentum space. The distribution of these quantities in momentum space gives insight into the fundamental nature of the ground state, insight that can be used to make contact with weak-coupling-based analytic approaches and to optimize numerical methods, the momentum-space DMRG in particular. We study the site and subsystem entropies as a function of interaction strength $U$ and system size. In both the one- and two-dimensional cases, we find that the subsystem entropy scales proportionally to ${U}^{2}$ for weak $U$ and proportionally to volume. Nevertheless, the optimized momentum-space DMRG can provide variationally accurate results for the two-dimensional Hubbard model at weak coupling for moderate system sizes.

arXiv:1508.07477•Physical Review B•Dec 8, 2015•34 citations

Entanglement in one-dimensional bosonic chains

C.A. Avila, Juan José Mendoza-Arenas, Roberto Franco +1 more

"We studied the one-dimensional Bose-Hubbard model, by means of the Density Matrix Renormalization Group (DMRG). The U parameter, which contains the Coulombian repulsion information, is fixed to U = 1 and it is performed a swept over the hopping parameter t. The entanglement is investigated by implementing the one-site, two-site and block von Neumann entropy looking for which of these measures give account of the quantum phase transition from Mott insulator to superfluid. The one-site and two-site entropy do not give information of the phase transition of the system. In contrast, the block entropy shows, in the Mott insulator phase a behavior that saturates rapidly to a finite value and in the superfluid phase shows a logarithmic behavior that agrees with the mean field theory. From the derivative of the block entropy average as a function of t we found the phase transition takes place at tc =0.15."

Revista Mexicana de Física•Jan 1, 2012•0 citations

A path integral ground state Monte Carlo algorithm for entanglement of lattice bosons

Emanuel Casiano-Diaz, C. M. Herdman, Adrian Del Maestro

A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.

arXiv:2207.11301•SciPost Physics•Mar 29, 2023•6 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

Even-odd entanglement in boson and spin systems

R. Rossignoli, N. Canosa, J. M. Matera

We examine the entanglement entropy of the even half of a translationally invariant finite chain or lattice in its ground state. This entropy measures the entanglement between the even and odd halves (each forming a ``comb'' of $n/2$ sites) and can be expected to be extensive for short-range couplings away from criticality. We first consider bosonic systems with quadratic couplings, where analytic expressions for arbitrary dimensions can be provided. The bosonic treatment is then applied to finite spin chains and arrays by means of the random-phase approximation. Results for first-neighbor anisotropic $\mathit{XY}$ couplings indicate that, while at strong magnetic fields this entropy is strictly extensive, at weak fields important deviations arise, stemming from parity-breaking effects and the presence of a factorizing field (in the vicinity of which it becomes size-independent and identical to the entropy of a contiguous half). Exact numerical results for small spin $s$ chains are shown to be in agreement with the bosonic random-phase approximation prediction.

arXiv:1104.4005•Physical Review A•Apr 22, 2011•6 citations
The area law relates the bipartite entanglement entropy of a quantum many-body ground state to the size of the boundary between the subsystems, but the geometry of this boundary is rarely discussed. We inspect this in the 1D Bose-Hubbard model at fixed density by comparing four spatial bipartitions of the periodic lattice: first half, second half, even sites, and odd sites; sharing the same number of sites but differing in how the boundary is arranged. We find analytical limits with perturbation theory: in the Mott insulator the contiguous cut obeys the area lay while the alternating cut obeys a volume law S∝NsS\propto N_sS∝Ns​, in this sense an anomaly, and for the superfluid both cuts colapse to the binomial saturation due to delocalization of the state. We formulate these limits as a statement about the many-body problem using a generalized slave-boson approach based on mean-field with quantum fluctuations while verifying with Exact Diagonalization (ED) for small lattice sizes and Densitiy Matrix Renormalization Group (DMRG) simulations for Ns≫1N_s\gg 1Ns​≫1. The slave-boson Gaussian ground state allows to compute the entanglement entropy from a reduced correlation matrix for any desired bipartition consistent with ED and DMRG results. Using slave bosons the computational cost is set by the local cutoff nmax⁡n_{\max}nmax​ rather than the Hilbert space dimension, so we can reach lattice sizes far beyond ED. Our method is capable of establishing the partition-dependent scaling laws as a many-body feature, not only a finite-size effect, in great agreement with the ED for Ns∈[4,10]N_s\in[4,10]Ns​∈[4,10] and DMRG for larger lattice sizes.
We study the entanglement entropy and entanglement spectrum of the paradigmatic Bose-Hubbard model, describing strongly correlated bosons on a lattice. The use of a controlled approximation-the slave-boson approach-allows us to study entanglement in all regimes of the model (and, most importantly, across its superfluid-Mott-insulator transition) at a minimal cost. We find that the area-law scaling of entanglement-verified in all the phases-exhibits a sharp singularity at the transition. The singularity is greatly enhanced when the transition is crossed at fixed, integer filling, due to a richer entanglement spectrum containing an additional gapless mode, which descends from the amplitude (Higgs) mode of the global excitation spectrum-while this mode remains gapped at the generic (commensurate-incommensurate) transition with variable filling. Hence, the entanglement properties contain a unique signature of the two different forms of bosonic criticality exhibited by the Bose-Hubbard model.
We study the properties of the ground states of the one- and two-dimensional Hubbard models at half-filling and moderate doping using entanglement-based measures, which we calculate numerically using the momentum-space density matrix renormalization group (DMRG). In particular, we investigate quantities such as the single-site entropy and two-site mutual information of single-particle momentum states as well as the behavior of the bipartite subsystem entropy for partitions in momentum space. The distribution of these quantities in momentum space gives insight into the fundamental nature of the ground state, insight that can be used to make contact with weak-coupling-based analytic approaches and to optimize numerical methods, the momentum-space DMRG in particular. We study the site and subsystem entropies as a function of interaction strength UUU and system size. In both the one- and two-dimensional cases, we find that the subsystem entropy scales proportionally to U2{U}^{2}U2 for weak UUU and proportionally to volume. Nevertheless, the optimized momentum-space DMRG can provide variationally accurate results for the two-dimensional Hubbard model at weak coupling for moderate system sizes.
"We studied the one-dimensional Bose-Hubbard model, by means of the Density Matrix Renormalization Group (DMRG). The U parameter, which contains the Coulombian repulsion information, is fixed to U = 1 and it is performed a swept over the hopping parameter t. The entanglement is investigated by implementing the one-site, two-site and block von Neumann entropy looking for which of these measures give account of the quantum phase transition from Mott insulator to superfluid. The one-site and two-site entropy do not give information of the phase transition of the system. In contrast, the block entropy shows, in the Mott insulator phase a behavior that saturates rapidly to a finite value and in the superfluid phase shows a logarithmic behavior that agrees with the mean field theory. From the derivative of the block entropy average as a function of t we found the phase transition takes place at tc =0.15."
A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.
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.
We examine the entanglement entropy of the even half of a translationally invariant finite chain or lattice in its ground state. This entropy measures the entanglement between the even and odd halves (each forming a ``comb'' of n/2n/2n/2 sites) and can be expected to be extensive for short-range couplings away from criticality. We first consider bosonic systems with quadratic couplings, where analytic expressions for arbitrary dimensions can be provided. The bosonic treatment is then applied to finite spin chains and arrays by means of the random-phase approximation. Results for first-neighbor anisotropic XY\mathit{XY}XY couplings indicate that, while at strong magnetic fields this entropy is strictly extensive, at weak fields important deviations arise, stemming from parity-breaking effects and the presence of a factorizing field (in the vicinity of which it becomes size-independent and identical to the entropy of a contiguous half). Exact numerical results for small spin sss chains are shown to be in agreement with the bosonic random-phase approximation prediction.