Recent advances in quantum Gibbs sampling leave open the central question of rapid mixing near and below phase transitions. This challenge is especially relevant for code Hamiltonians whose Gibbs states capture phenomena such as the thermal stability of quantum topological order. In this work, we formulate a new Markov chain, Code Swendsen-Wang dynamics, which uses global updates to prepare the Gibbs states of arbitrary code Hamiltonians. We establish Code Swendsen-Wang dynamics as the right generalization of Swendsen-Wang dynamics for the Ising model to quantum and classical code Hamiltonians: it mixes rapidly for all previously known code Hamiltonians with efficient Gibbs samplers, resolves the central open case of the 4D toric code, and meets fundamental barriers exactly at first-order phase transitions.
Recent advances in quantum Gibbs sampling leave open the central question of rapid mixing near and below phase transitions. This challenge is especially relevant for code Hamiltonians whose Gibbs states capture phenomena such as the thermal stability of quantum topological order. In this work, we formulate a new Markov chain, Code Swendsen-Wang dynamics, which uses global updates to prepare the Gibbs states of arbitrary code Hamiltonians. We establish Code Swendsen-Wang dynamics as the right generalization of Swendsen-Wang dynamics for the Ising model to quantum and classical code Hamiltonians: it mixes rapidly for all previously known code Hamiltonians with efficient Gibbs samplers, resolves the central open case of the 4D toric code, and meets fundamental barriers exactly at first-order phase transitions.
Gibbs state preparation, or Gibbs sampling, is a key computational technique extensively used in physics, statistics, and other scientific fields. Recent efforts for designing fast mixing Gibbs samplers for quantum Hamiltonians have largely focused on commuting local Hamiltonians (CLHs), a non-trivial subclass of Hamiltonians which include highly entangled systems such as the Toric code and quantum double model. Most previous Gibbs samplers relied on simulating the Davies generator, which is a Lindbladian associated with the thermalization process in nature.Instead of using the Davies generator, we design a different Gibbs sampler for various CLHs by giving a reduction to classical Hamiltonians, in the sense that one can efficiently prepare the Gibbs state for some CLH H on a quantum computer as long as one can efficiently do classical Gibbs sampling for the corresponding classical Hamiltonian H ( c ) . We demonstrate that our Gibbs sampler is able to replicate state-of-the-art results as well as prepare the Gibbs state in regimes which were previously unknown, such as the low temperature region, as long as there exists fast mixing Gibbs samplers for the corresponding classical Hamiltonians. Our reductions are as follows.– If H is a 2-local qudit CLH, then H ( c ) is a 2-local qudit classical Hamiltonian.– If H is a 4-local qubit CLH on 2D lattice and there are no classical qubits, then H ( c ) is a 2-local qudit classical Hamiltonian on a planar graph. As an example, our algorithm can prepare the Gibbs state for the (defected) Toric code at any non-zero temperature in O ( n 2 p o l y ( l o g n ) ) time.– If H is a 4-local qubit CLH on 2D lattice and there are classical qubits, assuming that quantum terms are uniformly correctable, then H ( c ) is a constant-local classical Hamiltonian.
We present hybrid Gibbs sampling algorithms for the stabilizer code Hamiltonians of the rotated surface code and the toric code with only local quantum algorithms, using ∼L/2 quantum circuit depth to prepare the Gibbs state of the rotated surface code Hamiltonian, and L quantum circuit depth to prepare the Gibbs state of the toric code Hamiltonian, being L the side of the side of the square lattice. We further show that if we allow for non-local gates, the Gibbs state of the periodic 1D Ising model can be prepared in logarithmic depth and linearly many simultaneous measurements.
We prove that a quasi-local quantum Gibbs sampler rapidly mixes to the Gibbs state of any 1D finite-range Hamiltonian at all finite temperatures in O(log(n/ε)) time. Each update traces out a contiguous block of constant length and reconstructs it using the Petz recovery map, following early proposals for quantum heat-bath Gibbs samplers by Kastoryano and Brandão~\cite{kastoryano2016quantum}. Our proof directly bounds the distance between the updated state and the target Gibbs state using quantum Wasserstein distance of order 1. We further provide O(polylog(n/ε)) depth of quantum circuits to prepare these Gibbs states. We hope these techniques will facilitate the analysis and design of other Gibbs samplers.
I study the mixing time of quantum Gibbs samplers (Davies / Chen– Kastoryano–Gilyén generators) at low temperature. I prove a locality lower bound: any local, KMS–detailed-balanced sampler whose target has an order- parameter bottleneck has an Arrhenius-small spectral gap, by a quantum Cheeger argument, and hence mixing time eΩ(β∆F). The bound is a quantum Cheeger argument and holds for every KMS–detailed-balanced sampler. The coherence cannot lift an order-parameter bottleneck, thereby the obstruction survives quantization. It is not, however, intrinsic to the target. I give a 3-local, detailed-balanced Lindbladian on an enlarged spin+bond register whose spin marginal realizes the Fortuin–Kasteleyn cluster update and which, for commuting ferromagnets, has an inverse-polynomial gap at every temperature and prepares σ in poly(n)log(1/ε) gates. This gives an unconditional separation from the local-reversible class in which locality (not coherence) is the operative resource. As a conjecture, I isolate the genuinely quantum question about cluster updates for noncommuting targets, where no nonnegative bond representation exists and the cluster must be grown coherently.
It is of great interest to understand the thermalization of open quantum many-body systems, and how quantum computers are able to efficiently simulate that process. A recently introduced dissipative evolution, inspired by existing models of open system thermalization, has been shown to be efficiently implementable on a quantum computer. Here, we prove that, at high enough temperatures, this evolution reaches the Gibbs state in time scaling logarithmically with system size. The result holds for Hamiltonians that satisfy the Lieb-Robinson bound, such as local Hamiltonians on a lattice, and includes long-range systems. To the best of our knowledge, these are the first results rigorously establishing the rapid mixing property of high-temperature quantum Gibbs samplers, which is known to give the fastest possible speed for thermalization in the many-body setting. We then apply our result to the problem of estimating partition functions at high temperature, showing an improved performance over previous classical and quantum algorithms.
The problem of simulating the thermal behavior of quantum systems remains a central open challenge in quantum computing. Unlike well-established quantum algorithms for unitary dynamics, algorithms for preparing thermal states—crucial for probing equilibrium behavior—became available only recently with breakthrough algorithms based on the simulation of well-designed dissipative processes, a quantum-analog to Markov chain Monte Carlo algorithms. We show a way to implement these algorithms avoiding expensive block encoding and relying only on dense local circuits, akin to Hamiltonian simulation. Specifically, our method leverages spatial truncation and Trotterization of exact quasilocal dissipative processes. We rigorously prove that the approximations we use have little effect on rapid mixing at high temperatures and allow convergence to the thermal state with small bounded error. Moreover, we accompany our analytical results with numerical simulations that show that this method, unlike previously thought, is within the reach of current generation of quantum hardware. These results provide the first provably efficient quantum thermalization protocol implementable on near-term quantum devices, offering a concrete path toward practical simulation of equilibrium quantum phenomena.