Variational Thermal State Preparation on Digital Quantum Processors Assisted by Matrix Product States
Pith reviewed 2026-05-18 03:59 UTC · model grok-4.3
The pith
A variational framework uses classical matrix product states to guide quantum circuits toward accurate finite-temperature Gibbs states.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By combining scalable classical matrix product state computation of the free energy with variational optimization of a hardware-efficient ansatz on quantum hardware, the method produces high-fidelity approximations to Gibbs states for one-dimensional lattices of 30 sites and two-dimensional lattices of 6 by 6 sites, with a working demonstration on a 156-qubit superconducting processor.
What carries the argument
Classical matrix product state evaluation of the Helmholtz free energy that guides variational optimization of a hardware-efficient quantum ansatz.
If this is right
- Observables including energy density, susceptibility, specific heat, and two-point correlations match exact analytic results in one dimension and quantum Monte Carlo results in two dimensions.
- The approach scales to systems requiring up to 44 qubits in numerical simulations.
- Error-mitigated execution on a 156-qubit IBM Heron processor yields usable approximations to the Gibbs state of a 30-site transverse-field Ising model.
- The framework directly supports thermal sampling tasks in quantum simulation and optimization.
Where Pith is reading between the lines
- Off-loading free-energy evaluation to classical matrix product states can lower the quantum circuit depth needed for thermal-state preparation.
- The same hybrid loop could be applied to other variational problems where classical pre-computation supplies a reliable cost function.
- Success on present-day hardware suggests the method may serve as a testbed for studying how classical assistance affects convergence in variational quantum algorithms.
Load-bearing premise
The variational optimization of the hardware-efficient ansatz, guided by the classical MPS free-energy value, converges to a quantum state whose measured observables agree with the true Gibbs state within the reported error bars.
What would settle it
A measurement on the 30-site transverse-field Ising model that shows energy or susceptibility after mitigation deviating from the MPS-predicted value by more than the claimed 50 percent error reduction.
Figures
read the original abstract
The preparation of quantum Gibbs states at finite temperatures is a cornerstone of quantum computation, enabling applications in quantum simulation of many-body systems, machine learning via quantum Boltzmann machines, and optimization through thermal sampling techniques. In this work, we introduce a variational framework that leverages matrix product states for the efficient classical evaluation of the Helmholtz free energy, combining scalable entanglement entropy computation with a hardware efficient ansatz to accurately approximate thermal states in one- and two-dimensional systems. We conduct extensive benchmarking on key observables, including energy density, susceptibility, specific heat, and two-point correlations, comparing against exact analytical results for 1D systems and quantum Monte Carlo simulations for 2D lattices across various temperatures and ansatz configurations. Our large-scale numerical simulations demonstrate the capability to prepare high-quality Gibbs states for 1D lattice models with up to 30 sites and 2D systems with up to 6x6 sites, using up to 44 qubits. Finally, we demonstrate the framework's practical viability on a 156-qubit IBM Heron processor by preparing the approximate Gibbs state of a 30-site transverse-field Ising model. Leveraging a combination of error mitigation techniques, we reduce the relative errors in energy and susceptibility measurements by over 50% compared to unmitigated results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a hybrid variational framework for preparing approximate quantum Gibbs states on digital quantum hardware. It employs classical matrix product states (MPS) to evaluate the Helmholtz free energy, guiding the optimization of parameters in a hardware-efficient ansatz. The approach is benchmarked via comparisons of energy density, susceptibility, specific heat, and correlation functions against exact 1D analytics and 2D quantum Monte Carlo results for lattices up to 30 sites (1D) and 6x6 (2D), using up to 44 qubits. A hardware demonstration prepares an approximate Gibbs state for a 30-site transverse-field Ising model on a 156-qubit IBM Heron processor, with error mitigation yielding over 50% reduction in relative errors for energy and susceptibility.
Significance. If the central claim holds, the work offers a practical route to thermal state preparation on near-term quantum devices by leveraging classical MPS for scalable free-energy guidance. This hybrid strategy could advance quantum simulation of finite-temperature many-body physics and related applications. Strengths include the scale of the numerical benchmarks (1D/2D systems with direct comparisons to exact/QMC data) and the real-hardware demonstration with mitigation; these provide concrete evidence of feasibility beyond small-system proofs of principle.
major comments (2)
- [Numerical Results / Benchmarking] The benchmarking sections (numerical simulations for 1D and 2D models) report close agreement between optimized observables and exact/QMC references, but provide limited analysis of ansatz expressivity or optimization convergence. Without diagnostics such as entanglement entropy matching the thermal value, fidelity estimates where feasible, or checks for multiple minima in the free-energy landscape, it remains possible that the hardware-efficient ansatz reaches states reproducing low-order observables within error bars while deviating from the true Gibbs state in higher moments or structure, particularly at lower temperatures. This directly affects the validity of the 'high-quality Gibbs states' claim.
- [Hardware Implementation] In the hardware demonstration section for the 30-site TFIM on the IBM Heron processor, the >50% error reduction via mitigation is presented as evidence of practical viability. However, the manuscript does not quantify how mitigation techniques (including any post-selection) alter the prepared state's proximity to the target thermal ensemble versus simply improving expectation values of the measured observables. This is load-bearing for interpreting the hardware results as successful Gibbs-state preparation rather than mitigated sampling of a variational state.
minor comments (2)
- [Abstract and Numerical Results] Clarify the precise qubit count and circuit depth for the 2D 6x6 simulations versus the 1D cases to avoid ambiguity in the 'up to 44 qubits' statement.
- [Method] The description of the MPS-assisted free-energy evaluation would benefit from an explicit equation or pseudocode showing how the entanglement entropy computation is combined with the variational parameters.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive comments, which have helped us improve the clarity and rigor of the manuscript. We address each major comment in detail below and have made revisions where appropriate to strengthen the presentation of our results.
read point-by-point responses
-
Referee: [Numerical Results / Benchmarking] The benchmarking sections (numerical simulations for 1D and 2D models) report close agreement between optimized observables and exact/QMC references, but provide limited analysis of ansatz expressivity or optimization convergence. Without diagnostics such as entanglement entropy matching the thermal value, fidelity estimates where feasible, or checks for multiple minima in the free-energy landscape, it remains possible that the hardware-efficient ansatz reaches states reproducing low-order observables within error bars while deviating from the true Gibbs state in higher moments or structure, particularly at lower temperatures. This directly affects the validity of the 'high-quality Gibbs states' claim.
Authors: We agree that additional diagnostics would strengthen the evidence for the quality of the approximated Gibbs states. In the revised manuscript, we have added a new subsection discussing the entanglement entropy of the optimized variational states for the 1D models, demonstrating close agreement with the exact thermal entanglement entropy obtained from the transfer-matrix method. We also include convergence plots of the free-energy optimization for representative temperatures and system sizes, showing consistent convergence across multiple random initial parameter sets with no evidence of trapping in significantly suboptimal minima. While direct fidelity estimation remains intractable for the largest systems considered, the combination of agreement on energy, susceptibility, specific heat, and correlation functions with the entropy diagnostic supports that the states capture the essential structure of the thermal ensemble beyond low-order observables. These additions are now reflected in the updated numerical results and discussion sections. revision: yes
-
Referee: [Hardware Implementation] In the hardware demonstration section for the 30-site TFIM on the IBM Heron processor, the >50% error reduction via mitigation is presented as evidence of practical viability. However, the manuscript does not quantify how mitigation techniques (including any post-selection) alter the prepared state's proximity to the target thermal ensemble versus simply improving expectation values of the measured observables. This is load-bearing for interpreting the hardware results as successful Gibbs-state preparation rather than mitigated sampling of a variational state.
Authors: We acknowledge the importance of this distinction. Because full quantum state tomography or direct fidelity computation to the target Gibbs state is infeasible on current hardware for a 30-site system, we cannot provide a quantitative metric of the prepared state's distance to the ideal thermal ensemble. The variational optimization itself is performed classically via MPS evaluation of the free energy, and the hardware step serves to prepare the resulting ansatz and extract observables. In the revised manuscript, we have clarified the role of error mitigation: it is applied post-preparation to reduce bias and variance in the measured expectation values, thereby enabling a more reliable comparison to the classically optimized targets. We have added an explicit discussion of this limitation, noting that hardware noise affects the fidelity of state preparation but that the mitigation improves the accuracy of the extracted observables without altering the underlying variational guarantee in the ideal case. This revised framing better contextualizes the hardware results as a demonstration of practical measurement rather than a direct claim of enhanced state fidelity. revision: partial
Circularity Check
No significant circularity; claims rest on external benchmarks and hardware execution
full rationale
The paper optimizes a hardware-efficient ansatz via classical MPS evaluation of the Helmholtz free energy, then validates the resulting quantum state by direct comparison of observables (energy density, susceptibility, specific heat, correlations) to exact 1D solutions and independent QMC simulations for 2D lattices up to 6x6. The 30-site TFIM demonstration on the 156-qubit IBM Heron processor further uses real hardware execution with error mitigation. These external references and the separation between classical MPS free-energy computation and quantum variational preparation mean the central claims do not reduce to self-definition, fitted-input renaming, or self-citation chains. The derivation chain is self-contained against independent benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- ansatz variational parameters
- MPS bond dimension
axioms (1)
- domain assumption Matrix product states provide an efficient and accurate classical representation for computing the Helmholtz free energy and entanglement entropy of the thermal states considered.
Reference graph
Works this paper leans on
-
[1]
Thermal energy density We first focus on estimating the energy density, defined as ε= ⟨H⟩ β N = Tr(ρβH) N ,(7) whereNis the number of spins in the system andρ β is the Gibbs state at inverse temperatureβ. In our work- flow, the Gibbs state is approximated by the variationally prepared mixed stateρ(θ ∗ β), whereθ ∗ β are the optimal circuit parameters forβ...
-
[2]
The factor ofN 2 in the denominators of Eqs
Magnetic susceptibility and specific heat Next, we present the thermal estimates of two other important observables, the magnetic susceptibility and specific heat [27], χ= β N 2 ⟨M 2 tot⟩β − ⟨Mtot⟩2 β ,(9) cv = β2 N 2 ⟨H 2⟩β − ⟨H⟩ 2 β ,(10) whereM tot =PN i=1 σz i is the total magnetization oper- ator, withσ z i being the Pauli-Zoperator acting on the i-t...
-
[3]
Two-point correlation functions As alluded to in the previous subsection, the locality of an operator can significantly impact the accuracy of its thermal expectation value estimated with the varia- tionally prepared Gibbs state. We have seen that the susceptibility estimates for the 1D TFIM are noticeably less accurate than the energy estimates across a ...
-
[4]
I. M. Georgescu, S. Ashhab, and F. Nori, Quantum sim- ulation, Rev. Mod. Phys.86, 153 (2014)
work page 2014
-
[5]
´A. M. Alhambra, Quantum Many-Body Systems in Ther- mal Equilibrium, PRX Quantum4, 040201 (2023)
work page 2023
-
[6]
M. H. Amin, E. Andriyash, J. Rolfe, B. Kulchytskyy, and R. Melko, Quantum Boltzmann Machine, Phys. Rev. X 8, 021050 (2018)
work page 2018
- [7]
-
[8]
S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, Optimiza- tion by Simulated Annealing, Science220, 671 (1983)
work page 1983
-
[9]
F. G. Brandao and K. M. Svore, Quantum Speed-Ups for Solving Semidefinite Programs, in2017 IEEE 58th Annual Symposium on Foundations of Computer Science (FOCS)(IEEE, Berkeley, CA, 2017) pp. 415–426
work page 2017
- [10]
-
[11]
R. Oliveira and B. M. Terhal, The complexity of quan- tum spin systems on a two-dimensional square lattice, Quantum Info. Comput.8, 900 (2008)
work page 2008
-
[12]
C.-F. Chen, M. J. Kastoryano, F. G. S. L. Brand˜ ao, and A. Gily´ en, Quantum Thermal State Preparation (2023), arXiv:2303.18224 [math-ph, physics:quant-ph]
work page internal anchor Pith review arXiv 2023
-
[13]
D. Poulin and P. Wocjan, Sampling from the Thermal Quantum Gibbs State and Evaluating Partition Func- tions with a Quantum Computer, Phys. Rev. Lett.103, 220502 (2009)
work page 2009
-
[14]
A. N. Chowdhury and R. D. Somma, Quantum algo- rithms for Gibbs sampling and hitting-time estimation (2016), arXiv:1603.02940 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2016
- [15]
-
[16]
C.-F. Chen and F. G. S. L. Brand˜ ao, Fast Thermalization from the Eigenstate Thermalization Hypothesis (2023), arXiv:2112.07646 [quant-ph]
-
[17]
P. Wocjan and K. Temme, Szegedy Walk Unitaries for Quantum Maps, Commun. Math. Phys.402, 3201 (2023)
work page 2023
-
[18]
P. Rall, C. Wang, and P. Wocjan, Thermal State Prepa- ration via Rounding Promises, Quantum7, 1132 (2023)
work page 2023
-
[19]
L. Lin, Dissipative Preparation of Many-Body Quantum States: Towards Practical Quantum Advantage (2025), arXiv:2505.21308 [quant-ph]
-
[20]
Y. S. Teo, H. Zhu, B.-G. Englert, J. ˇReh´ aˇ cek, and Z. Hradil, Quantum-State Reconstruction by Maximiz- ing Likelihood and Entropy, Phys. Rev. Lett.107, 020404 (2011)
work page 2011
-
[21]
K. Banaszek, M. Cramer, and D. Gross, Focus on quan- tum tomography, New J. Phys.15, 125020 (2013)
work page 2013
- [22]
- [23]
-
[24]
D. Zhu, S. Johri, N. M. Linke, K. A. Landsman, C. Huerta Alderete, N. H. Nguyen, A. Y. Matsuura, T. H. Hsieh, and C. Monroe, Generation of thermofield double states and critical ground states with a quantum com- 15 puter, Proc. Natl. Acad. Sci. U.S.A.117, 25402 (2020)
work page 2020
-
[25]
R. Sagastizabal, S. P. Premaratne, B. A. Klaver, M. A. Rol, V. Negˆ ırneac, M. S. Moreira, X. Zou, S. Johri, N. Muthusubramanian, M. Beekman, C. Zachariadis, V. P. Ostroukh, N. Haider, A. Bruno, A. Y. Mat- suura, and L. DiCarlo, Variational preparation of finite- temperature states on a quantum computer, npj Quan- tum Inf7, 130 (2021)
work page 2021
- [26]
-
[27]
Y. Wang, G. Li, and X. Wang, Variational Quantum Gibbs State Preparation with a Truncated Taylor Series, Phys. Rev. Applied16, 054035 (2021)
work page 2021
-
[28]
M. Consiglio, J. Settino, A. Giordano, C. Mastroianni, F. Plastina, S. Lorenzo, S. Maniscalco, J. Goold, and T. J. G. Apollaro, Variational Gibbs state preparation on noisy intermediate-scale quantum devices, Phys. Rev. A110, 012445 (2024)
work page 2024
-
[29]
B. Sambasivam, K. Sherbert, K. Shirali, N. J. Mayhall, E. Barnes, and S. E. Economou, TEPID-ADAPT: Adap- tive variational method for simultaneous preparation of low-temperature Gibbs and low-lying eigenstates (2025), arXiv:2503.14490 [quant-ph]
-
[30]
X. Wang, X. Feng, T. Hartung, K. Jansen, and P. Stor- nati, Critical behavior of the Ising model by preparing the thermal state on a quantum computer, Phys. Rev. A 108, 022612 (2023)
work page 2023
-
[31]
Y. Ilin and I. Arad, Dissipative Variational Quantum Algorithms for Gibbs State Preparation, IEEE Transac- tions on Quantum Engineering6, 1 (2025)
work page 2025
-
[32]
A. Francis, D. Zhu, C. Huerta Alderete, S. Johri, X. Xiao, J. K. Freericks, C. Monroe, N. M. Linke, and A. F. Kemper, Many-body thermodynamics on quantum com- puters via partition function zeros, Science Advances7, eabf2447 (2021)
work page 2021
-
[33]
Adiabatic preparation of thermal states and entropy-noise relation on noisy quantum computers
E. Granet and H. Dreyer, Adiabatic preparation of ther- mal states and entropy-noise relation on noisy quantum computers (2025), arXiv:2509.05206 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[34]
R. Haghshenas, J. Gray, A. C. Potter, and G. K.-L. Chan, Variational Power of Quantum Circuit Tensor Networks, Phys. Rev. X12, 011047 (2022)
work page 2022
- [35]
-
[36]
N. Robertson, A. Akhriev, J. Vala, and S. Zhuk, Approx- imate Quantum Compiling for Quantum Simulation: A Tensor Network Based Approach, ACM Trans. Quantum Comput.6, 1 (2025)
work page 2025
-
[37]
B. Jaderberg, G. Pennington, K. V. Marshall, L. W. Anderson, A. Agarwal, L. P. Lindoy, I. Rungger, S. Mensa, and J. Crain, Variational preparation of nor- mal matrix product states on quantum computers (2025), arXiv:2503.09683 [quant-ph]
-
[38]
R. Or´ us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)
work page 2014
-
[39]
M. M. Wolf, F. Verstraete, M. B. Hastings, and J. I. Cirac, Area Laws in Quantum Systems: Mutual Infor- mation and Correlations, Phys. Rev. Lett.100, 070502 (2008)
work page 2008
-
[40]
M. Kliesch, C. Gogolin, M. J. Kastoryano, A. Riera, and J. Eisert, Locality of Temperature, Phys. Rev. X4, 031019 (2014)
work page 2014
-
[41]
T. Kuwahara, ´A. M. Alhambra, and A. Anshu, Improved Thermal Area Law and Quasilinear Time Algorithm for Quantum Gibbs States, Phys. Rev. X11, 011047 (2021)
work page 2021
-
[42]
J. Gray, Quimb: A python package for quantum informa- tion and many-body calculations, Journal of Open Source Software3, 819 (2018)
work page 2018
-
[43]
A Quantum Approximate Optimization Algorithm
E. Farhi, J. Goldstone, and S. Gutmann, A Quan- tum Approximate Optimization Algorithm (2014), arXiv:1411.4028 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[44]
A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, Hardware- efficient variational quantum eigensolver for small molecules and quantum magnets, Nature549, 242 (2017)
work page 2017
-
[45]
M. J. D. Powell, A Direct Search Optimization Method That Models the Objective and Constraint Functions by Linear Interpolation, inAdvances in Optimization and Numerical Analysis, edited by S. Gomez and J.-P. Hen- nart (Springer Netherlands, Dordrecht, 1994) pp. 51–67
work page 1994
- [46]
-
[47]
A. P. Young, Finite-temperature and dynamical proper- ties of the random transverse-field Ising spin chain, Phys. Rev. B56, 11691 (1997)
work page 1997
-
[48]
Watrous,The Theory of Quantum Information(Cam- bridge University Press, Cambridge, 2018)
J. Watrous,The Theory of Quantum Information(Cam- bridge University Press, Cambridge, 2018)
work page 2018
-
[49]
Wolff, Collective Monte Carlo Updating for Spin Sys- tems, Phys
U. Wolff, Collective Monte Carlo Updating for Spin Sys- tems, Phys. Rev. Lett.62, 361 (1989)
work page 1989
-
[50]
S. Hesselmann, Thermal Ising transitions in the vicinity of two-dimensional quantum critical points, Phys. Rev. B93, 10.1103/PhysRevB.93.155157 (2016)
-
[51]
Sachdev,Quantum Phase Transitions, 2nd ed
S. Sachdev,Quantum Phase Transitions, 2nd ed. (Cam- bridge University Press, Cambridge, 2011)
work page 2011
-
[52]
J. J. Wallman, Noise tailoring for scalable quantum com- putation via randomized compiling, Phys. Rev. A94, 10.1103/PhysRevA.94.052325 (2016)
-
[53]
Y. Kim, A. Eddins, S. Anand, K. X. Wei, E. Van Den Berg, S. Rosenblatt, H. Nayfeh, Y. Wu, M. Zaletel, K. Temme, and A. Kandala, Evidence for the utility of quantum computing before fault tolerance, Nature618, 500 (2023)
work page 2023
-
[54]
E. Van Den Berg, Z. K. Minev, A. Kandala, and K. Temme, Probabilistic error cancellation with sparse Pauli–Lindblad models on noisy quantum processors, Nat. Phys.19, 1116 (2023)
work page 2023
-
[55]
E. B¨ aumer, V. Tripathi, D. S. Wang, P. Rall, E. H. Chen, S. Majumder, A. Seif, and Z. K. Minev, Ef- ficient Long-Range Entanglement Using Dynamic Cir- cuits, PRX Quantum5, 030339 (2024)
work page 2024
-
[56]
E. B¨ aumer and S. Woerner, Measurement-based long- range entangling gates in constant depth, Phys. Rev. Re- search7, 10.1103/physrevresearch.7.023120 (2025)
-
[57]
A. Kandala, K. Temme, A. D. C´ orcoles, A. Mezzacapo, J. M. Chow, and J. M. Gambetta, Error mitigation ex- tends the computational reach of a noisy quantum pro- cessor, Nature567, 491 (2019)
work page 2019
-
[58]
J. R. McClean, S. Boixo, V. N. Smelyanskiy, R. Babbush, and H. Neven, Barren plateaus in quantum neural net- work training landscapes, Nat Commun9, 4812 (2018). 16
work page 2018
- [59]
-
[60]
J. M. K¨ ubler, A. Arrasmith, L. Cincio, and P. J. Coles, An Adaptive Optimizer for Measurement-Frugal Varia- tional Algorithms, Quantum4, 263 (2020)
work page 2020
-
[61]
H. R. Grimsley, S. E. Economou, E. Barnes, and N. J. Mayhall, An adaptive variational algorithm for exact molecular simulations on a quantum computer, Nat Com- mun10, 3007 (2019)
work page 2019
- [62]
- [63]
- [64]
- [65]
-
[66]
S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett.69, 2863 (1992)
work page 1992
-
[67]
U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics January 2011 Special Issue,326, 96 (2011)
work page 2011
-
[68]
G. Carleo and M. Troyer, Solving the quantum many- body problem with artificial neural networks, Science 355, 602 (2017)
work page 2017
-
[69]
From architectures to applications: A review of neural quan- tum states,
H. Lange, A. V. de Walle, A. Abedinnia, and A. Bohrdt, From Architectures to Applications: A Review of Neural Quantum States (2024), arXiv:2402.09402 [cond-mat]
-
[70]
J. Martyn and B. Swingle, Product spectrum ansatz and the simplicity of thermal states, Phys. Rev. A100, 032107 (2019)
work page 2019
- [71]
- [72]
-
[73]
T. Baumgratz, D. Gross, M. Cramer, and M. B. Plenio, Scalable Reconstruction of Density Matrices, Phys. Rev. Lett.111, 020401 (2013)
work page 2013
-
[74]
J. Wang, Scalable quantum tomography with fi- delity estimation, Phys. Rev. A101, 10.1103/Phys- RevA.101.032321 (2020)
- [75]
-
[76]
Y. Teng, R. Samajdar, K. Van Kirk, F. Wilde, S. Sachdev, J. Eisert, R. Sweke, and K. Najafi, Learn- ing Topological States from Randomized Measurements Using Variational Tensor-Network Tomography, PRX Quantum6, 040303 (2025)
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.