Simon Becker, Cambyse Rouzé, Robert Salzmann
4 min
Abstract
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. 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. Our argument in the multi-mode case is based on a finite-rank reduction of the dissipative dynamics, which allows us to control the generator via compact perturbations and deduce the discreteness of the spectrum and the stability of the gap. We apply our results to provide efficient preparation of the corresponding Gibbs state on qubit hardware, and by that a quantum algorithm to compute thermal properties of the associated 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.
Alex: Right—like controlled tweaks to known-good systems preserve quick settling. Does this lead to actual computations, say for properties like energy or correlations?
Sam: Yes, it enables rigorous runtimes for preparing those states on qubit-based quantum hardware. For instance, to estimate free energy—which measures the system's stability at temperature, like total cost balancing energy and entropy—they connect a simple quadratic hopping Hamiltonian to the full model along a path, computing averages of the interaction term in states along the way. The paper shows this works with polynomial resources in system size and precision, a notable efficiency over classical limits.
Alex: So the perturbations unlock practical thermal properties without truncation errors. That's a meaningful step for bosonic simulations.
Sam: Precisely. It lays a foundation for quantum advantages in many-body thermal tasks, like optical lattices, though quantifying gap scaling remains open. The approach—reference models plus finite-rank stability—offers a robust path forward.
Alex: That's a grounded step forward for handling these infinite-particle challenges. Thanks for breaking it down, Sam. And that's our look at quantum Gibbs sampling for Bose-Hubbard thermal states. Thanks for listening to ResearchPod.