Clustering Theorem for Bose-Hubbard class Gibbs states

We establish the exponential clustering of correlation functions for the high-temperature Gibbs states of Bose-Hubbard type models. To overcome the technical difficulties arising from the unboundedness of bosonic operators, we develop the interaction-picture cluster-expansion technique. This method also allows us to systematically bound the moments of the local particle number. This result provides an analytical justification for the low-boson-density condition frequently assumed in the study of bosonic many-body systems. As direct mathematical consequences of the clustering property, we derive a uniform upper bound on the specific heat density and establish a bosonic thermal area law with improved temperature dependence.

💡 Research Summary

This paper addresses a long‑standing open problem in quantum statistical mechanics: establishing exponential clustering of correlations for lattice boson systems at high temperature. While exponential decay of correlations (clustering) is well‑understood for spin and fermionic models with finite‑dimensional local Hilbert spaces, bosonic models pose additional challenges because each lattice site carries an infinite‑dimensional Hilbert space and the fundamental creation and annihilation operators are unbounded. Consequently, standard cluster‑expansion techniques, which rely on uniform operator bounds, cannot be applied directly.

The authors consider a broad class of Bose‑Hubbard‑type Hamiltonians defined on a finite, simple, connected graph (V,E) with bounded degree. The Hamiltonian consists of nearest‑neighbour hopping, possible parametric squeezing terms, and an on‑site repulsive interaction quadratic in the particle number, together with a chemical potential. The parameters (hopping amplitudes, squeezing strengths, interaction strengths, and chemical potentials) are assumed to be uniformly bounded. By invoking the Kato‑Rellich theorem, they first prove that the Hamiltonian is essentially self‑adjoint on the dense subspace of finite‑particle vectors and that it is bounded from below by a quadratic function of the local number operators. This lower bound guarantees that the Gibbs operator e^{‑βH} is trace‑class for any inverse temperature β>0, so the Gibbs state ρ_β = e^{‑βH}/Tr(e^{‑βH}) is well defined.

The central methodological innovation is the development of an “interaction‑picture cluster expansion”. The Hamiltonian is split into a non‑interacting part I (the hopping and squeezing terms) and an on‑site potential W. In the interaction picture with respect to W, the Gibbs operator is written as
e^{‑βH}=e^{‑βW} 𝒯 exp\bigl(‑∫₀^{β} I(τ)dτ\bigr),
where I(τ)=e^{τW} I e^{‑τW} and 𝒯 denotes time ordering. This representation allows the authors to treat the unbounded hopping operators as time‑dependent bounded perturbations, because the on‑site quadratic potential dominates and controls the growth of particle numbers.

A key technical step is the derivation of a “low‑boson‑density inequality”. Using the interaction‑picture expansion together with combinatorial bounds on the number of connected edge subsets (controlled by a graph growth constant σ), the authors prove that for any site x and any integer s≥1, the s‑th moment of the local number operator satisfies
⟨n_x^{s}⟩_β ≤ C^{s} s! ,
where C depends only on β and the model parameters but not on the system size. This factorial bound replaces the ad‑hoc low‑density assumption that has been used in many previous rigorous works on bosons.

With this inequality in hand, the authors establish the main clustering theorem. For two observables A_X and B_Y supported on disjoint regions X and Y, the connected correlation function obeys
|⟨A_X B_Y⟩_β – ⟨A_X⟩_β⟨B_Y⟩_β| ≤ ‖A_X‖‖B_Y‖ K e^{‑α dist(X,Y)} ,
where the decay rate α is proportional to the inverse temperature (α≈c β^{‑1}) and the prefactor K depends only on the maximal degree of the graph, the growth constant σ, and the bounded interaction parameters. Importantly, the bound is uniform in the total number of sites |V|, demonstrating true thermodynamic‑limit clustering.

The paper then explores two immediate physical consequences. First, a uniform upper bound on the specific heat density is derived, leading to a “quasi‑Dulong‑Petit” law: at sufficiently high temperature the specific heat per site remains bounded by a constant independent of β. Second, a bosonic thermal area law is proved. For any region Λ⊂V, the von Neumann entropy of the reduced Gibbs state satisfies
S(ρ_Λ) ≤ K′ |∂Λ| β^{‑γ} ,
with γ>0 reflecting an improved temperature dependence compared with earlier results. This shows that thermal entanglement in bosonic lattice systems scales with the boundary area even when the local Hilbert space is infinite‑dimensional.

The authors discuss optimality of the temperature threshold, compare their results with earlier works based on Feynman‑Kac representations, coherent‑state functional integrals, or small‑field assumptions, and argue that their interaction‑picture technique can be adapted to more general settings: graphs with long‑range edges, inhomogeneous interaction patterns, or driven systems where particle number is not conserved.

In conclusion, the paper delivers the first rigorous proof of exponential clustering for high‑temperature Gibbs states of a wide class of Bose‑Hubbard models, simultaneously providing a mathematically solid justification for the low‑boson‑density condition and yielding new bounds on thermodynamic quantities such as specific heat and thermal entanglement. The interaction‑picture cluster expansion introduced here opens a promising avenue for tackling other unbounded‑operator problems in quantum many‑body physics.