We construct lower bounds to the spectral gap of a family of Lindblad generators known as Davies maps. These maps describe the thermalization of quantum systems weakly coupled to a heat bath. The steady state of these systems is given by the Gibbs distribution with respect to the system Hamiltonian. The bounds can be evaluated explicitly, when the eigenbasis and the spectrum of the Hamiltonian is known. A crucial assumption is that the spectrum of the Hamiltonian is non-degenerate. Furthermore, we provide a counterexample to the conjecture, that the convergence rate is always determined by the gap of the associated Pauli master equation. We conclude that the full dynamics of the Lindblad generator has to be considered. Finally, we present several physical example systems for which the bound to the spectral gap is evaluated.
We construct lower bounds to the spectral gap of a family of Lindblad generators known as Davies maps. These maps describe the thermalization of quantum systems weakly coupled to a heat bath. The steady state of these systems is given by the Gibbs distribution with respect to the system Hamiltonian. The bounds can be evaluated explicitly, when the eigenbasis and the spectrum of the Hamiltonian is known. A crucial assumption is that the spectrum of the Hamiltonian is non-degenerate. Furthermore, we provide a counterexample to the conjecture, that the convergence rate is always determined by the gap of the associated Pauli master equation. We conclude that the full dynamics of the Lindblad generator has to be considered. Finally, we present several physical example systems for which the bound to the spectral gap is evaluated.
Many physical phenomena, including thermalization in open quantum systems and quantum Gibbs sampling, are modeled by Lindbladians approximating a system weakly coupled to a bath. Understanding the convergence speed of these Lindbladians to their steady states is crucial for bounding algorithmic runtimes and thermalization timescales. We study two such families of processes: one characterizing a repeated-interaction Gibbs sampling algorithm, and another modeling open many-body quantum thermalization. We prove that both converge in polynomial time for several non-commuting systems, including high-temperature local lattices, weakly interacting fermions, and 1D spin chains. These results demonstrate that simple dissipative quantum algorithms can prepare complex Gibbs states and that Lindblad dynamics accurately capture thermal relaxation. Our proofs rely on a novel technical result that extrapolates spectral gap lower bounds from quasi-local Lindbladians to the non-local generators governing these dynamics.
Despite extensive study, our understanding of quantum Markov chains remains far less complete than that of their classical counterparts. [Temme'13] observed that the Davies Lindbladian, a well-studied model of quantum Markov dynamics, contains an embedded classical Markov generator, raising the natural question of how the convergence properties of the quantum and classical dynamics are related. While [Temme'13] showed that the spectral gap of the Davies Lindbladian can be much smaller than that of the embedded classical generator for certain highly structured Hamiltonians, we show that if the spectrum of the Hamiltonian does not contain long arithmetic progressions, then the two spectral gaps must be comparable. As a consequence, we prove that for a large class of Hamiltonians, including those obtained by perturbing a fixed Hamiltonian with a generic external field, the quantum spectral gap remains within a constant factor of the classical spectral gap. Our result aligns with physical intuition and enables the application of classical Markov chain techniques to the quantum setting. The proof is based on showing that any ``off-diagonal'' eigenvector of the Davies generator can be used to construct an observable which commutes with the Hamiltonian and has a Lindbladian Rayleigh quotient which can be upper bounded in terms of that of the original eigenvector's Lindbladian Rayleigh quotient. Thus, a spectral gap for such observables implies a spectral gap for the full Davies generator.
The mean-field Heisenberg ferromagnet is a quantum spin model on the complete graph with isotropic spin-1/2 interactions. This non-commuting Hamiltonian is permutation and SU(2) invariant, and its Gibbs states undergo an SU(2) symmetry breaking phase transition at inverse temperature β=2. We consider the associated Davies generator, a canonical model of open-system thermalization, and prove tight asymptotic estimates for its spectral gap at all noncritical temperatures. For fixed β<2, the gap as a function of number of qubits n is Θ(1), while for fixed β>2 the gap is Θ(n−1). The matching upper bound of the spectral gap is witnessed by the total magnetization order parameter, suggesting that the low-temperature (β>2) slowdown is associated with broken continuous symmetry. Two key ingredients in our approach are a comparison argument, which introduces auxiliary generators to bound dissipation on nontrivial representations of the symmetry groups SU(2) and Sn, and a decomposition of the space of observables into spherical tensor operators to reveal a form of monotonicity.
Preparing Gibbs states through dissipative dynamics requires controlling convergence for the chosen Hamiltonian H and coupling operators. However, efficient implementation requires convergence guarantees from the spectral gap of the generator of the dynamics, which remains a crucial challenge. Recent results address the challenge by restricting the dynamics with structural assumptions about the Hamiltonian and its interactions. For the Chen-Kastoryano-Gilyen (CKG) construction and arbitrary finite-dimensional Hamiltonians, we prove an explicit comparison that converts a spectral gap estimate for the infinite-temperature generator into a lower bound on the spectral gap of the CKG quantum Gibbs sampler at every finite positive inverse temperature β. It allows a single spectral gap estimate at infinite temperature to certify convergence across Hamiltonians and temperatures. Such an estimate is often known from a classical spectral gap estimate, such as the classical random walk spectral gap. For Hodge Laplacians of simplicial complexes as Hamiltonians, we construct coupling operators whose generator at infinite temperature has the same spectral gap as a classical simplicial down-up walk. A positive spectral gap of the walk and bounded spectral width W=λmax(H)−λmin(H) give a generator spectral gap independent of the number of simplices at every fixed temperature. We then establish conditions under which the sampler approximately prepares zero-energy (harmonic) states in polynomial evolution time. This provides a conditional convergence guarantee for the state preparation in the thermal approach to quantum topological data analysis.
Starting from an arbitrary full-rank state of a lattice quantum spin system, we define a "canonical purified Hamiltonian" and characterize its spectral gap in terms of a spatial mixing condition (or correlation decay) of the state. When the state considered is a Gibbs state of a local, commuting Hamiltonian at positive temperature, we show that the spectral gap of the canonical purified Hamiltonian provides a lower bound to the spectral gap of a large class of reversible generators of quantum Markov semigroup, including local and ergodic Davies generators. As an application of our construction, we show that the mixing condition is always satisfied for any finite-range 1D model, as well as by Kitaev's quantum double models.
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.