- Open Access
Efficient quantum thermal state preparation via local driving: Lindbladian simulation with provable guarantees
Phys. Rev. B 114, 014302 – Published 1 July, 2026
DOI: https://doi.org/10.1103/b6mc-w34s
Abstract
Preparing the thermal density matrix corresponding to a given Hamiltonian is a task of central interest across quantum many-body physics, and is particularly salient when attempting to study it with quantum computers. Although solved in principle by recent constructions of efficiently simulable Lindblad master equations—that provably have as a steady state [C.-F. Chen et al., Nature 646, 561 (2025)]—the implementation of these “exact Gibbs samplers” requires large-scale quantum computing resources and is hence challenging in practice on current or even near-term quantum devices. Here, we propose a scheme for approximately simulating an exact Gibbs sampler up to a rigorously bounded error that only requires the (repeated) implementation of three readily available ingredients: (a) analog simulation of ; (b) strictly local but time-dependent couplings to ancilla qubits; and (c) reset of the ancillas. We give rigorous guarantees on the difference between the fixed point reached by our protocol and the exact thermal state, which only depend on parameters of the protocol and its mixing time. The procedure is efficiently implementable on near-term devices if is local and the mixing time scales mildly with both system size and protocol parameters. While guaranteeing the latter for Hamiltonians of interest remains an important problem for future work, here we lay the groundwork for developing fully efficient thermal state preparation protocols on quantum simulators.
Physics Subject Headings (PhySH)
Article Text
References (77)
- S. Lloyd, Universal quantum simulators, Science 273, 1073 (1996).
- B. M. Terhal and D. P. DiVincenzo, Problem of equilibration and the computation of correlation functions on a quantum computer, Phys. Rev. A 61, 022301 (2000).
- A. Riera, C. Gogolin, and J. Eisert, Thermalization in nature and on a quantum computer, Phys. Rev. Lett. 108, 080402 (2012).
- Y. Ge, A. Molnár, and J. I. Cirac, Rapid adiabatic preparation of injective projected entangled pair states and Gibbs states, Phys. Rev. Lett. 116, 080503 (2016).
- E. Mozgunov and D. Lidar, Completely positive master equation for arbitrary driving and small level spacing, Quantum 4, 227 (2020).
- S. McArdle, T. Jones, S. Endo, Y. Li, S. C. Benjamin, and X. Yuan, Variational ansatz-based quantum simulation of imaginary time evolution, npj Quantum Inf. 5, 75 (2019).
- O. Shtanko and R. Movassagh, Preparing thermal states on noiseless and noisy programmable quantum processors, arXiv:2112.14688.
- A. Schuckert, A. Bohrdt, E. Crane, and M. Knap, Probing finite-temperature observables in quantum simulators of spin systems with short-time dynamics, Phys. Rev. B 107, L140410 (2023).
- X. Mi, A. A. Michailidis, S. Shabani, K. C. Miao, P. V. Klimov, J. Lloyd, E. Rosenberg, R. Acharya, I. Aleiner, T. I. Andersen, et al., Stable quantum-correlated many-body states through engineered dissipation, Science 383, 1332 (2024).
- J. Lloyd, A. A. Michailidis, X. Mi, V. Smelyanskiy, and D. A. Abanin, Quasiparticle cooling algorithms for quantum many-body state preparation, PRX Quantum 6, 010361 (2025).
- A. Deshpande, M. Hinsche, S. Najafi, K. Sharma, R. Sweke, and C. Zoufal, Dynamic parameterized quantum circuits: Expressive and barren-plateau free, arXiv:2411.05760.
- M. Consiglio, Variational quantum algorithms for Gibbs state preparation, in Numerical Computations: Theory and Algorithms (Springer Nature, Switzerland, 2025), pp. 56–70.
- D. Poulin and P. Wocjan, Sampling from the thermal quantum Gibbs state and evaluating partition functions with a quantum computer, Phys. Rev. Lett. 103, 220502 (2009).
- E. Bilgin and S. Boixo, Preparing thermal states of quantum systems by dimension reduction, Phys. Rev. Lett. 105, 170405 (2010).
- K. Temme, T. J. Osborne, K. G. Vollbrecht, D. Poulin, and F. Verstraete, Quantum Metropolis sampling, Nature (London) 471, 87 (2011).
- P. Rall, C. Wang, and P. Wocjan, Thermal state preparation via rounding promises, Quantum 7, 1132 (2023).
- M.-H. Yung and A. Aspuru-Guzik, A quantum–quantum metropolis algorithm, Proc. Natl. Acad. Sci. USA 109, 754 (2012).
- A. N. Chowdhury and R. D. Somma, Quantum algorithms for Gibbs sampling and hitting-time estimation, Quantum Inf. Comput. 17, 41 (2017).
- M. Motta, C. Sun, A. T. K. Tan, M. J. O'Rourke, E. Ye, A. J. Minnich, F. G. S. L. Brandão, and G. K.-L. Chan, Determining eigenstates and thermal states on a quantum computer using quantum imaginary time evolution, Nat. Phys. 16, 205 (2020).
- Z. Holmes, G. Muraleedharan, R. D. Somma, Y. Subasi, and B. Şahinoğlu, Quantum algorithms from fluctuation theorems: Thermal-state preparation, Quantum 6, 825 (2022).
- P. Wocjan and K. Temme, Szegedy walk unitaries for quantum maps, Commun. Math. Phys. 402, 3201 (2023).
- D. Zhang, J. L. Bosse, and T. Cubitt, Dissipative quantum Gibbs sampling, arXiv:2304.04526.
- J. Jiang and S. Irani, Quantum Metropolis sampling via weak measurement, arXiv:2406.16023.
- M. J. Kastoryano and F. G. S. L. Brandão, Quantum Gibbs samplers: The commuting case, Commun. Math. Phys. 344, 915 (2016).
- C.-F. Chen and F. G. S. L. Brandão, Fast thermalization from the eigenstate thermalization hypothesis, arXiv:2112.07646.
- C.-F. Chen, M. J. Kastoryano, F. G. S. L. Brandão, and A. Gilyén, Quantum thermal state preparation, arXiv:2303.18224.
- C.-F. Chen, M. J. Kastoryano, and A. Gilyén, An efficient and exact noncommutative quantum Gibbs sampler, arXiv:2311.09207.
- A. Gilyén, C.-F. Chen, J. F. Doriguello, and M. J. Kastoryano, Quantum generalizations of Glauber and Metropolis dynamics, arXiv:2405.20322.
- Z. Ding, B. Li, and L. Lin, Efficient quantum Gibbs samplers with Kubo–Martin–Schwinger detailed balance condition, Commun. Math. Phys. 406, 67 (2025).
- C.-F. Chen, M. Kastoryano, F. G. S. L. Brandão, and A. Gilyén, Efficient quantum thermal simulation, Nature 646, 561 (2025).
- E. B. Davies, Markovian master equations, Commun. Math. Phys. 39, 91 (1974).
- E. B. Davies, Markovian master equations. II, Math. Ann. 219, 147 (1976).
- J. Guo, O. Hart, C.-F. Chen, A. J. Friedman, and A. Lucas, Designing open quantum systems with known steady states: Davies generators and beyond, Quantum 9, 1612 (2025).
- G. S. Agarwal, Open quantum Markovian systems and the microreversibility, Z. Phys. A: Hadrons Nucl. 258, 409 (1973).
- R. Alicki, On the detailed balance condition for non-Hamiltonian systems, Rep. Math. Phys. 10, 249 (1976).
- A. Kossakowski, A. Frigerio, V. Gorini, and M. Verri, Quantum detailed balance and KMS condition, Commun. Math. Phys. 57, 97 (1977).
- F. Fagnola and V. Umanita, Generators of detailed balance quantum Markov semigroups, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 10, 335 (2007).
- G. E. Santoro, R. Martonak, E. Tosatti, and R. Car, Theory of quantum annealing of an Ising spin glass, Science 295, 2427 (2002).
- B. Altshuler, H. Krovi, and J. Roland, Anderson localization makes adiabatic quantum optimization fail, Proc. Natl. Acad. Sci. USA 107, 12446 (2010).
- V. Bapst, L. Foini, F. Krzakala, G. Semerjian, and F. Zamponi, The quantum adiabatic algorithm applied to random optimization problems: The quantum spin glass perspective, Phys. Rep. 523, 127 (2013).
- T. Rakovszky, B. Placke, N. P. Breuckmann, and V. Khemani, Bottlenecks in quantum channels and finite temperature phases of matter, arXiv:2412.09598.
- B. Placke, T. Rakovszky, N. P. Breuckmann, and V. Khemani, Topological quantum spin glass order and its realization in qLDPC codes, arXiv:2412.13248.
- D. Gamarnik, B. T. Kiani, and A. Zlokapa, Slow mixing of quantum Gibbs samplers, arXiv:2411.04300.
- E. R. Anschuetz, Efficient learning implies quantum glassiness, Commun. Math. Phys. 407, 30 (2026).
- H. Chen, B. Li, J. Lu, and L. Ying, A randomized method for simulating Lindblad equations and thermal state preparation, Quantum 9, 1917 (2025).
- E. Brunner, L. Coopmans, G. Matos, M. Rosenkranz, F. Sauvage, and Y. Kikuchi, Lindblad engineering for quantum Gibbs state preparation under the eigenstate thermalization hypothesis, Quantum 9, 1843 (2025).
- H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, Oxford, 2002).
- L. A. Correa and J. Glatthard, Potential renormalisation, Lamb shift and mean-force Gibbs state–To shift or not to shift? Quantum 9, 1868 (2025).
- As in the classical case, this is sufficient but not necessary.
- E. Schrödinger, Statistical Thermodynamics (Dover Publications, Mineola, NY, 1967).
- G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys. 48, 119 (1976).
- K. Sharma and M. C. Tran, Hamiltonian simulation in the interaction picture using the Magnus expansion, arXiv:2404.02966.
- M. J. Kastoryano and K. Temme, Quantum logarithmic Sobolev inequalities and rapid mixing, J. Math. Phys. 54, 052202 (2013).
- K. Temme, M. J. Kastoryano, M. B. Ruskai, M. M. Wolf, and F. Verstraete, The -divergence and mixing times of quantum Markov processes, J. Math. Phys. 51, 122201 (2010).
- M. Žnidarič, Relaxation times of dissipative many-body quantum systems, Phys. Rev. E 92, 042143 (2015).
- I. Bardet, A. Capel, L. Gao, A. Lucia, D. Pérez-García, and C. Rouzé, Rapid thermalization of spin chain commuting Hamiltonians, Phys. Rev. Lett. 130, 060401 (2023).
- J. Kochanowski, A. M. Alhambra, A. Capel, and C. Rouzé, Rapid thermalization of dissipative many-body dynamics of commuting Hamiltonians, Commun. Math. Phys. 406, 176 (2025).
- Y. Tong and Y. Zhan, Fast mixing of weakly interacting fermionic systems at any temperature, arXiv:2501.00443.
- Z. Ding, Z. Landau, B. Li, L. Lin, and R. Zhang, Polynomial-time preparation of low-temperature Gibbs states for 2D toric code, arXiv:2410.01206.
- Š. Šmíd, R. Meister, M. Berta, and R. Bondesan, Polynomial time quantum Gibbs sampling for Fermi-Hubbard model at any temperature, Nat. Commun. 16, 10736 (2025).
- C. Rouzé, D. S. França, and Á. M. Alhambra, Efficient thermalization and universal quantum computing with quantum Gibbs samplers, Nat. Phys. (2026).
- C. Rouzé, D. S. França, and Á. M. Alhambra, Optimal quantum algorithm for Gibbs state preparation, Phys. Rev. Lett. 136, 060601 (2026).
- J. Lloyd and D. A. Abanin, Quantum thermal state preparation for near-term quantum processors, arXiv:2506.21318.
- Z. Ding, Y. Zhan, J. Preskill, and L. Lin, End-to-end efficient quantum thermal and ground state preparation made simple, arXiv:2508.05703.
- H. Kim and D. A. Huse, Ballistic spreading of entanglement in a diffusive nonintegrable system, Phys. Rev. Lett. 111, 127205 (2013).
- J. Dormand and P. Prince, A family of embedded Runge-Kutta formulae, J. Comput. Appl. Math. 6, 19 (1980).
- A. J. Daley, I. Bloch, C. Kokail, S. Flannigan, N. Pearson, M. Troyer, and P. Zoller, Practical quantum advantage in quantum simulation, Nature (London) 607, 667 (2022).
- C. Ong, S. A. Parameswaran, B. Placke, and D. Hahn, Rigorous error bounds for dissipative thermal state preparation from weak system-bath coupling, arXiv:2605.03011.
- H. Chen, Z. Ding, and R. Zhang, Overcoming the Lamb shift in system-bath interaction models via KMS detailed balance: High-accuracy thermalization with time-bounded interactions, arXiv:2604.15616.
- S. R. White, Minimally entangled typical quantum states at finite temperature, Phys. Rev. Lett. 102, 190601 (2009).
- T. Kuwahara, T. Mori, and K. Saito, Floquet–magnus theory and generic transient dynamics in periodically driven many-body quantum systems, Ann. Phys. 367, 96 (2016).
- E. M. Stoudenmire and S. R. White, Minimally entangled typical thermal state algorithms, New J. Phys. 12, 055026 (2010).
- M. Rigol, T. Bryant, and R. R. P. Singh, Numerical linked-cluster approach to quantum lattice models, Phys. Rev. Lett. 97, 187202 (2006).
- This follows since , where we used the triangle inequality to rewrite the commutator in terms of a product and then applied Hölder's inequality to rewrite the numerator.
- The original argument in Ref. [26] is in terms of matrices, but we can view these as superoperators written in the doubled space where density matrices are vectorized, so we give the argument in terms of superoperators here for brevity.
- R. Bhatia, Variational principles for eigenvalues, in Matrix Analysis (Springer, New York, 1997), pp. 57–83.
- Note that for any matrix , which yields the identity in Eq. (C22) if indexing is understood to be modulo .