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:2605.03017

Preparing High-Fidelity Thermofield Double States

Brian J. J. Khor, Nadie LiTenn, Martin Sasieta, Brian Swingle

quant-phcond-mat.stat-mechhep-thSubmitted May 4, 2026• Announced May 6, 2026• 0 citations
arXivHTMLPDF

Abstract

A major promise of quantum computers is the controlled preparation of many-body quantum states beyond the reach of efficient classical computation. Among the most important targets are thermal mixed states and their thermofield double (TFD) purifications, which play central roles in quantum many-body physics and quantum gravity. For target systems with a bounded energy spectrum that obey the eigenstate thermalization hypothesis (ETH), we present a parent Hamiltonian built from two copies of the target Hamiltonian and ultra-local couplings between the copies, which we argue is gapped with a ground state that approximates a TFD state of the target Hamiltonian. By adiabatically evolving down from strong coupling, we can thus prepare a high-fidelity TFD state. We study two variants of the parent Hamiltonian using numerical methods in two classes of models: mixed field Ising models in one and two dimensions and non-local "spin Sachdev-Ye-Kitaev'' models. In the simpler variant, the parent Hamiltonian ground state has high overlap with a TFD for system sizes accessible to near-term quantum devices. However, the global overlap decays exponentially with the number of qubits, with a small error per degree of freedom. The second variant introduces an additional penalty term which can be tuned to reduce or remove the decay of the overlap with system size. Together with a general ETH-based analysis, these results suggest a broadly applicable method for TFD preparation that is not limited to particular temperatures or geometric locality.

Comments
13 pages + appendices

Related papers

Product spectrum ansatz and the simplicity of thermal states

John M. Martyn, Brian G. Swingle

Calculating the physical properties of quantum thermal states is a difficult problem for classical computers, rendering it intractable for most quantum many-body systems. A quantum computer, by contrast, would make many of these calculations feasible in principle, but it is still nontrivial to prepare a given thermal state or sample from it. It is also not known how to prepare special simple purifications of thermal states known as thermofield doubles, which play an important role in quantum many-body physics and quantum gravity. To address this problem, we propose a variational scheme to prepare approximate thermal states on a quantum computer by applying a series of two-qubit gates to a product mixed state. We apply our method to a nonintegrable region of the mixed field Ising chain and the Sachdev-Ye-Kitaev model. We also demonstrate how our method can be easily extended to large systems governed by local Hamiltonians and the preparation of thermofield double states. By comparing our results with exact solutions, we find that our construction enables the efficient preparation of approximate thermal states on quantum devices. Our results can be interpreted as implying that the details of the many-body energy spectrum are not needed to capture simple thermal observables.

arXiv:1812.01015•Physical Review A•Sep 6, 2019•64 citations

Preparing a Thermofield Double State with Feedback Quantum Algorithms

Guilherme E. L. Pexe, Lucas A. M. Rattighieri, Felipe F. Fanchini +2 more

The efficient preparation of correlated thermal states, such as the Thermofield Double (TFD) state, is a fundamental prerequisite for simulating quantum gravity models and many-body thermodynamics on quantum processors. In this work, we investigate the ground state preparation of the Two Coupled Sachdev-Ye-Kitaev model, known as the Maldacena-Qi model, which is dual to a traversable wormhole in $AdS_2$, utilizing feedback-based quantum algorithms. We demonstrate that the standard feedback-based quantum algorithm (FALQON) and its time-rescaled variant (TR-FALQON) face severe kinetic limitations in this system, failing to converge to the highly entangled ground state when initialized in trivial product states. To overcome these barriers, we propose the hybrid ITE-TR-FALQON protocol, which integrates the imaginary-time evolution present in imaginary-time-enhanced FALQON (ITE-FALQON) with the time-rescaling mechanism. Our numerical results indicate that the introduction of non-unitary dynamics is strictly necessary to break symmetry traps and filter out excited states, while time-rescaling drastically accelerates algorithm convergence. The proposed method achieves fidelities close to unity and reproduces the von Neumann and Rényi entropy spectra of the exact TFD state with high precision.

hep-thCross-listedarXiv:2607.01653•Physical review. D/Physical review. D.•Jul 3, 2026•0 citations

Variational preparation of thermofield double states for SYK models via multi-angle QAOA: sequential angle pruning for circuit reduction

Haji Muhammad Husnain Ashfaq, Moongul Byun, Keun-Young Kim

Variational preparation of thermofield double (TFD) states can require deep quantum circuits, particularly for interacting many-body systems. Reducing these circuits while retaining high fidelity is therefore crucial for TFD-state preparation on noisy quantum processors. We study this problem by applying the multi-angle quantum approximate optimization algorithm (ma-QAOA) to TFD-state preparation and introducing two top-down sequential angle-pruning algorithms. Starting from the optimized initial ma-QAOA circuit, both algorithms sequentially remove Pauli-string evolutions with small optimized angles and reoptimize the remaining parameters after each removal. We apply these algorithms to Gaussian and binary Sachdev--Ye--Kitaev (SYK) models in both dense and sparse cases. We find that ma-QAOA prepares the target TFD states with high fidelity and that sequential small-angle pruning retains high fidelity while reducing the circuit depth, particularly at low temperature. Moreover, using the post-reoptimization cost in sequential small-angle pruning further improves the fidelity. For the binary sparse $N=10$ SYK model at $β=10$, $88.8\%$--$92.1\%$ of the nonlocal Pauli-string evolutions are removed while retaining an average fidelity of approximately $95\%$. Finally, we propose extensions of the sequential pruning algorithms toward quantum--classical hybrid implementation.

quant-pharXiv:2609.02793•Sep 3, 2026•0 citations

Cooling the Sachdev-Ye-Kitaev model using thermofield double states

Thomas Schuster, Bryce Kobrin, Vincent P. Su +2 more

We analyze a simple and efficient experimental protocol to cool the Sachdev-Ye-Kitaev (SYK) model to low temperatures. The protocol utilizes local couplings between two copies of an SYK model to create a gapped adiabatic path, between a high temperature product state and a low temperature thermofield double state. By smoothly varying the coupling strength between these two limits, one efficiently cools the SYK model. We support these predictions-and demonstrate fast cooling to the low-temperature gravitational regime-via exact numerical solutions to the large-N equations of motion that govern the ground state and dynamical properties of the coupled system. Finally, we present a theoretical framework based upon eigenstate thermalization that provides a microscopic explanation for the efficacy of the cooling protocol; intriguingly, this suggests that the protocol may be applicable for cooling strongly-interacting quantum Hamiltonians more broadly.

quant-pharXiv:2511.09620•Nov 14, 2025•0 citations

Double-bracket quantum algorithms for thermal state preparation

Andrew Wright, Reyhaneh Aghaei Saem, Supanut Thanasilp +2 more

We propose quantum algorithms for preparing thermal states via the simulation of the thermofield double states. The key idea is to leverage double-bracket quantum algorithms to implement imaginary-time evolution on thermofield double states, whose reduced state realizes the Gibbs state. Our method, termed double-bracket thermofield double (DB-TFD), introduces two variants. The first, the vanilla DB-TFD algorithm, directly implements imaginary-time evolution using double-bracket quantum imaginary-time evolution. The second, poly DB-TFD, employs double-bracket quantum signal processing to approximate the imaginary-time evolution operator via a polynomial transformation. We demonstrate that the complexity of the poly DB-TFD algorithm scales exponentially with the inverse temperature in a broad practical regime. This scaling is consistent with existing methods, and numerical simulations support the corresponding theoretical bound. We further demonstrate the utility of DB-TFD in quantum Boltzmann machines for generative modeling, achieving improved performance compared with variational imaginary-time evolution approaches. These results establish DB-TFD as a promising route for thermal state preparation in the near-term and early-fault-tolerant regimes.

quant-pharXiv:2606.05947•Jun 5, 2026•0 citations

Variational Thermal Quantum Simulation via Thermofield Double States

Jingxiang Wu, Timothy H. Hsieh

We present a variational approach for quantum simulators to realize finite temperature Gibbs states by preparing thermofield double (TFD) states. Our protocol is motivated by the quantum approximate optimization algorithm and involves alternating time evolution between the Hamiltonian of interest and interactions which entangle the system and its auxiliary counterpart. As a simple example, we demonstrate that thermal states of the 1D classical Ising model at any temperature can be prepared with perfect fidelity using L/2 iterations, where L is system size. We also show that a free fermion TFD can be prepared with nearly optimal efficiency. Given the simplicity and efficiency of the protocol, our approach enables near-term quantum platforms to access finite temperature phenomena via preparation of thermofield double states.

arXiv:1811.11756•Physical Review Letters•Nov 26, 2019•135 citations
A major promise of quantum computers is the controlled preparation of many-body quantum states beyond the reach of efficient classical computation. Among the most important targets are thermal mixed states and their thermofield double (TFD) purifications, which play central roles in quantum many-body physics and quantum gravity. For target systems with a bounded energy spectrum that obey the eigenstate thermalization hypothesis (ETH), we present a parent Hamiltonian built from two copies of the target Hamiltonian and ultra-local couplings between the copies, which we argue is gapped with a ground state that approximates a TFD state of the target Hamiltonian. By adiabatically evolving down from strong coupling, we can thus prepare a high-fidelity TFD state. We study two variants of the parent Hamiltonian using numerical methods in two classes of models: mixed field Ising models in one and two dimensions and non-local "spin Sachdev-Ye-Kitaev'' models. In the simpler variant, the parent Hamiltonian ground state has high overlap with a TFD for system sizes accessible to near-term quantum devices. However, the global overlap decays exponentially with the number of qubits, with a small error per degree of freedom. The second variant introduces an additional penalty term which can be tuned to reduce or remove the decay of the overlap with system size. Together with a general ETH-based analysis, these results suggest a broadly applicable method for TFD preparation that is not limited to particular temperatures or geometric locality.
Calculating the physical properties of quantum thermal states is a difficult problem for classical computers, rendering it intractable for most quantum many-body systems. A quantum computer, by contrast, would make many of these calculations feasible in principle, but it is still nontrivial to prepare a given thermal state or sample from it. It is also not known how to prepare special simple purifications of thermal states known as thermofield doubles, which play an important role in quantum many-body physics and quantum gravity. To address this problem, we propose a variational scheme to prepare approximate thermal states on a quantum computer by applying a series of two-qubit gates to a product mixed state. We apply our method to a nonintegrable region of the mixed field Ising chain and the Sachdev-Ye-Kitaev model. We also demonstrate how our method can be easily extended to large systems governed by local Hamiltonians and the preparation of thermofield double states. By comparing our results with exact solutions, we find that our construction enables the efficient preparation of approximate thermal states on quantum devices. Our results can be interpreted as implying that the details of the many-body energy spectrum are not needed to capture simple thermal observables.
The efficient preparation of correlated thermal states, such as the Thermofield Double (TFD) state, is a fundamental prerequisite for simulating quantum gravity models and many-body thermodynamics on quantum processors. In this work, we investigate the ground state preparation of the Two Coupled Sachdev-Ye-Kitaev model, known as the Maldacena-Qi model, which is dual to a traversable wormhole in AdS2AdS_2AdS2​, utilizing feedback-based quantum algorithms. We demonstrate that the standard feedback-based quantum algorithm (FALQON) and its time-rescaled variant (TR-FALQON) face severe kinetic limitations in this system, failing to converge to the highly entangled ground state when initialized in trivial product states. To overcome these barriers, we propose the hybrid ITE-TR-FALQON protocol, which integrates the imaginary-time evolution present in imaginary-time-enhanced FALQON (ITE-FALQON) with the time-rescaling mechanism. Our numerical results indicate that the introduction of non-unitary dynamics is strictly necessary to break symmetry traps and filter out excited states, while time-rescaling drastically accelerates algorithm convergence. The proposed method achieves fidelities close to unity and reproduces the von Neumann and Rényi entropy spectra of the exact TFD state with high precision.
Variational preparation of thermofield double (TFD) states can require deep quantum circuits, particularly for interacting many-body systems. Reducing these circuits while retaining high fidelity is therefore crucial for TFD-state preparation on noisy quantum processors. We study this problem by applying the multi-angle quantum approximate optimization algorithm (ma-QAOA) to TFD-state preparation and introducing two top-down sequential angle-pruning algorithms. Starting from the optimized initial ma-QAOA circuit, both algorithms sequentially remove Pauli-string evolutions with small optimized angles and reoptimize the remaining parameters after each removal. We apply these algorithms to Gaussian and binary Sachdev--Ye--Kitaev (SYK) models in both dense and sparse cases. We find that ma-QAOA prepares the target TFD states with high fidelity and that sequential small-angle pruning retains high fidelity while reducing the circuit depth, particularly at low temperature. Moreover, using the post-reoptimization cost in sequential small-angle pruning further improves the fidelity. For the binary sparse N=10N=10N=10 SYK model at β=10β=10β=10, 88.8%88.8\%88.8%--92.1%92.1\%92.1% of the nonlocal Pauli-string evolutions are removed while retaining an average fidelity of approximately 95%95\%95%. Finally, we propose extensions of the sequential pruning algorithms toward quantum--classical hybrid implementation.
We analyze a simple and efficient experimental protocol to cool the Sachdev-Ye-Kitaev (SYK) model to low temperatures. The protocol utilizes local couplings between two copies of an SYK model to create a gapped adiabatic path, between a high temperature product state and a low temperature thermofield double state. By smoothly varying the coupling strength between these two limits, one efficiently cools the SYK model. We support these predictions-and demonstrate fast cooling to the low-temperature gravitational regime-via exact numerical solutions to the large-N equations of motion that govern the ground state and dynamical properties of the coupled system. Finally, we present a theoretical framework based upon eigenstate thermalization that provides a microscopic explanation for the efficacy of the cooling protocol; intriguingly, this suggests that the protocol may be applicable for cooling strongly-interacting quantum Hamiltonians more broadly.
We propose quantum algorithms for preparing thermal states via the simulation of the thermofield double states. The key idea is to leverage double-bracket quantum algorithms to implement imaginary-time evolution on thermofield double states, whose reduced state realizes the Gibbs state. Our method, termed double-bracket thermofield double (DB-TFD), introduces two variants. The first, the vanilla DB-TFD algorithm, directly implements imaginary-time evolution using double-bracket quantum imaginary-time evolution. The second, poly DB-TFD, employs double-bracket quantum signal processing to approximate the imaginary-time evolution operator via a polynomial transformation. We demonstrate that the complexity of the poly DB-TFD algorithm scales exponentially with the inverse temperature in a broad practical regime. This scaling is consistent with existing methods, and numerical simulations support the corresponding theoretical bound. We further demonstrate the utility of DB-TFD in quantum Boltzmann machines for generative modeling, achieving improved performance compared with variational imaginary-time evolution approaches. These results establish DB-TFD as a promising route for thermal state preparation in the near-term and early-fault-tolerant regimes.
We present a variational approach for quantum simulators to realize finite temperature Gibbs states by preparing thermofield double (TFD) states. Our protocol is motivated by the quantum approximate optimization algorithm and involves alternating time evolution between the Hamiltonian of interest and interactions which entangle the system and its auxiliary counterpart. As a simple example, we demonstrate that thermal states of the 1D classical Ising model at any temperature can be prepared with perfect fidelity using L/2 iterations, where L is system size. We also show that a free fermion TFD can be prepared with nearly optimal efficiency. Given the simplicity and efficiency of the protocol, our approach enables near-term quantum platforms to access finite temperature phenomena via preparation of thermofield double states.