Abstract 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.
Abstract 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.
Convergence properties of the variational cluster approach with respect to the variational parameter space, cluster size, and boundary conditions of the reference system are investigated and discussed for bosonic many-body systems. Specifically, the variational cluster approach is applied to the one-dimensional Bose-Hubbard model, which exhibits a quantum phase transition from Mott to superfluid phase. In order to benchmark the variational cluster approach, results for the phase boundary delimiting the first Mott lobe are compared with essentially exact density matrix renormalization group data. Furthermore, static quantities, such as the ground state energy and the one-particle density matrix are compared with high-order strong coupling perturbation theory results. For reference systems with open boundary conditions the variational parameter space is extended by an additional variational parameter which allows for a more uniform particle density on the reference system and thus drastically improves the results. It turns out that the variational cluster approach yields accurate results with relatively low-computational effort for both the phase boundary as well as the static properties of the one-dimensional Bose-Hubbard model, even at the tip of the first Mott lobe where correlation effects are most pronounced.
The entanglement area law is a universal principle that characterizes quantum many-body phases and underpins tensor network algorithms. Traditionally, its validity has been limited to systems with short-range interactions and bounded local energy. Achieving a complete generalization that removes both of these constraints has been a longstanding goal in quantum many-body theory, especially for interacting boson systems where unbounded energy presents intrinsic difficulties. In this work, we rigorously prove the area law for one-dimensional interacting boson systems with long-range interactions, covering broad models including the Bose-Hubbard and ϕ4 classes. Furthermore, we establish an efficiency guarantee for Matrix-Product-State approximations of the ground states, offering a practical route to numerical simulation. One of our main technical contributions is a general method for Hilbert space dimension reduction, whose applicability extends to arbitrary spatial dimensions. These results address two major challenges simultaneously and provide important foundations for simulating long-range cold atomic systems. Entanglement area law has been established in systems with short range interactions and bounded local energy. Here the authors extend the proof to 1D bosonic systems with long-range interactions and unbounded local energy and outline implications for matrix product state approximations of the ground states.
While recent advances have established efficient quantum algorithms for preparing Gibbs states of finite-dimensional systems, comparable complexity results for bosonic and other infinite-dimensional models remain unexplored. We introduce the first general rigorous Gibbs sampling framework for bosonic many-body systems, showing that physically relevant bosonic models admit gapped dissipative generators, enabling efficient preparation of thermal states, provided that the spectral gap scales favourably with the number of modes. Although our results hold for broad classes of models, we illustrate them using Bose-Hubbard Hamiltonians, both within and beyond the mean-field regime. In both cases, we show that the associated dissipative generators maintain a positive spectral gap, thereby implying exponential convergence to the thermal state. For the Gibbs sampler corresponding to the full Bose-Hubbard model, we obtain a spectral-gap lower bound Ce−cn, where n is the number of modes. We apply our results to provide a Gibbs-state preparation algorithm on qubit hardware, with runtime polynomial in the number of modes and the inverse spectral gap, and thereby obtain a quantum algorithm to compute thermal properties of the model. This provides the first mathematically controlled route to Gibbs sampling in infinite-dimensional systems, with implications for quantum simulation, thermalization, and many-body complexity, where quantum advantages may arise.
We consider the many-body quantum Gibbs state for the Bose–Hubbard model on a finite graph at positive temperature. We scale the interaction with the inverse temperature, corresponding to a mean-field limit where the temperature is of the order of the average particle number. For this model it is known that the many-body Gibbs state converges, as temperature goes to infinity, to the Gibbs measure of a discrete nonlinear Schrödinger equation, i.e., a Gibbs measure defined in terms of a one-body theory. In this article we extend these results by proving an expansion to any order of the many-body Gibbs state with inverse temperature as a small parameter. The coefficients in the expansion can be calculated as vacuum expectation values using a recursive formula, and we compute the first two coefficients explicitly.
In the context of C* dynamical systems, we consider a locally compact group G acting by *-automorphisms on a C* algebra U of observables, and assume a state of U that satisfies the clustering property with respect to a net of group elements of G. That is, the two-point connected correlation function vanishes in the limit on the net, when one observable is translated under the group action. Then we show that all higher order connected correlation functions (Ursell functions, or classical cumulants) and all free correlation functions (free cumulants, from free probability) vanish at the same rate in that limit. Additionally, we show that mean clustering, also called ergodicity, extends to higher order correlations. We then apply those results to equilibrium states of quantum spin lattice models. Under certain assumptions on the range of the interaction, high temperature Gibbs states are known to be exponentially clustering w.r.t. space translations. Combined with the Lieb-Robinson bound, one obtains exponential clustering for space-time translations outside the Lieb-Robinson light-cone. Therefore, by our present results, all the higher order connected and free correlation functions will vanish exponentially under space-time translations outside the Lieb-Robinson light cone, in high temperature Gibbs states. Another consequence is that their long-time averaging over a space-time ray vanishes for almost every ray velocity.
A physical system is said to satisfy a thermal area law if the mutual information between two adjacent regions in the Gibbs state is controlled by the area of their boundary. Lattice bosons have recently gained significant interest because they can be precisely tuned in experiments and bosonic codes can be employed in quantum error correction to circumvent classical no-go theorems. However, the proofs of many basic information-theoretic inequalities such as the thermal area law break down for bosons because their interactions are unbounded. Here, we rigorously derive a thermal area law for a class of bosonic Hamiltonians in any dimension which includes the paradigmatic Bose-Hubbard model. The main idea to go beyond bounded interactions is to introduce a quasi-free reference state with artificially decreased chemical potential by means of a double Peierls-Bogoliubov estimate.