Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2508.00126 · Physical Review Research

Efficient and simple Gibbs state preparation of the two-dimensional toric code via duality to classical Ising chains

Pablo Páez-Velasco, Niclas Schilling, Samuel O. Scalet, Frank Verstraete, Ángela Capel

Published May 15, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce the notion of polynomial-depth duality transformations, which relates two sets of operator algebras through a conjugation by a poly-depth quantum circuit, and make use of this to construct efficient Gibbs samplers for a variety of interesting quantum Hamiltonians as they are poly-depth dual to classical Hamiltonians. This is, for example, the case for the two-dimensional toric code, which is demonstrated to be poly-depth dual to two decoupled classical Ising spin chains for any system size, and we give evidence that such dualities hold for a wide class of stabilizer Hamiltonians. Additionally, we extend the above notion of duality to Lindbladians in order to show that mixing times and other quantities such as the spectral gap or the modified logarithmic Sobolev inequality are preserved under duality.

DOI
10.1103/sbxq-sqvp

Related papers

Efficient quantum Gibbs sampling of stabilizer codes using hybrid computation

Ivan H. C. Shum, Angela Capel

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 $\sim 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.

quant-pharXiv:2511.10839•New Journal of Physics•Nov 17, 2025•0 citations

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

Zhiyan Ding, Zeph A. Landau, Bowen Li +2 more

We propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

Journal of Mathematical Physics•May 1, 2026•1 citation

Gibbs state preparation for commuting Hamiltonian: Mapping to classical Gibbs sampling

Yeongwoo Hwang, Jiaqing Jiang

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.

arXiv:2410.04909•Quantum•Sep 18, 2026•1 citation

Rapid mixing of quantum spin chains at any finite temperature

Leeseok Kim

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 $\mathcal{O}\left(\log (n/\varepsilon)\right)$ 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 $\mathcal{O}(\text{polylog}(n/\varepsilon))$ depth of quantum circuits to prepare these Gibbs states. We hope these techniques will facilitate the analysis and design of other Gibbs samplers.

quant-pharXiv:2610.01190•Oct 2, 2026•0 citations

Fast Mixing of Quantum Spin Chains at All Temperatures

Thiago Bergamaschi, Chi-Fang Chen

It is shown that every one-dimensional Hamiltonian with short-range interactions admits a quantum Gibbs sampler [CKG23] with a system-size independent spectral gap at all finite temperatures. Consequently, their Gibbs states can be prepared in polylogarithmic depth, and satisfy exponential clustering of correlations, generalizing [Ara69].

arXiv:2510.08533•ACM Symposium on Theory of Computing (STOC)•Jun 9, 2026•1 citation

Code Swendsen-Wang Dynamics

Dominik Hangleiter, Nathan Ju, Umesh V. Vazirani

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.

arXiv:2510.08446•Oct 9, 2025•0 citations
We introduce the notion of polynomial-depth duality transformations, which relates two sets of operator algebras through a conjugation by a poly-depth quantum circuit, and make use of this to construct efficient Gibbs samplers for a variety of interesting quantum Hamiltonians as they are poly-depth dual to classical Hamiltonians. This is, for example, the case for the two-dimensional toric code, which is demonstrated to be poly-depth dual to two decoupled classical Ising spin chains for any system size, and we give evidence that such dualities hold for a wide class of stabilizer Hamiltonians. Additionally, we extend the above notion of duality to Lindbladians in order to show that mixing times and other quantities such as the spectral gap or the modified logarithmic Sobolev inequality are preserved under duality.
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\sim L/2∼L/2 quantum circuit depth to prepare the Gibbs state of the rotated surface code Hamiltonian, and LLL quantum circuit depth to prepare the Gibbs state of the toric code Hamiltonian, being LLL 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 propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.
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 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/ε))\mathcal{O}\left(\log (n/\varepsilon)\right)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 111. We further provide O(polylog(n/ε))\mathcal{O}(\text{polylog}(n/\varepsilon))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.
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.