Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Imaginary Time Evolution and Ground State Preparation Using Unitary Multi-Copy Protocols

Tal Schwartzman1,*, Torsten V. Zache2,3,†, Hannes Pichler2,3, and H. R. Sadeghpour1

  • *Contact author: tal.schwartzman@cfa.harvard.edu
  • †Contact author: torsten.zache@uibk.ac.at

PRX Quantum 7, 033033 – Published 18 August, 2026

DOI: https://doi.org/10.1103/1tf6-bc55

Abstract

Efficient low-energy state preparation is a key objective in quantum computation and quantum simulation. Quantum imaginary-time evolution replaces real-time dynamics with imaginary-time dynamics, exponentially suppressing higher-energy eigenstates. We introduce deterministic unitary protocols that approximate imaginary-time evolution for ground state preparation. The protocols require multiple copies of the system, real-time evolution under the system Hamiltonian, and controlled-swap operations (or more general swap-generated unitaries). Our analysis focuses on two concrete circuit families: a tree architecture with provable polynomial-in-depth convergence but rapidly growing width, and a compact “hedge” architecture that achieves comparable accuracy with only polynomial width in a heuristic construction supported by numerics. Numerical evidence indicates that mid-circuit post-selection can accelerate convergence with practical success probabilities. Separately, we demonstrate that circuit volume can be traded for the shot complexity of post-circuit observable estimation in the ground state preparation setting. Finally, we outline concrete platform-specific implementations in which multi-copy registers and swap-mediated couplings are natural, illustrating how these hybrid analog-digital circuits can complement existing state-preparation methods in the near term.

View figure in article

Physics Subject Headings (PhySH)

Popular Summary

Article Text

References (83)

  1. S. Sachdev, Quantum phase transitions, 2nd ed. (Cambridge University Press, Cambridge, England, 2011).
  2. S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry, Rev. Mod. Phys. 92, 015003 (2020).
  3. E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, Quantum computation by adiabatic evolution, arXiv:Quant-Ph/0001106.
  4. D. Aharonov, W. van Dam, J. Kempe, Z. Landau, S. Lloyd, and O. Regev, Adiabatic quantum computation is equivalent to standard quantum computation, SIAM Rev. 50, 755 (2008).
  5. J. Kempe, A. Kitaev, and O. Regev, The complexity of the local Hamiltonian problem, in International Conference on Foundations of Software Technology and Theoretical Computer Science (Springer, New York, 2004), pp. 372–383.
  6. S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett. 69, 2863 (1992).
  7. U. Schollwöck, The density-matrix renormalization group in the age of matrix product states, Ann. Phys. (Amsterdam) 326, 96 (2011).
  8. R. Orús, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Ann. Phys. (Amsterdam) 349, 117 (2014).
  9. Y. Ge, J. Tura, and J. I. Cirac, Faster ground state preparation and high-precision ground energy estimation with fewer qubits, J. Math. Phys. 60, 022202 (2019).
  10. A. Gilyén, Y. Su, G. H. Low, and N. Wiebe, Quantum singular value transformation and beyond: Exponential improvements for quantum matrix arithmetics, in Proceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing (Association for Computing Machinery, New York, 2019), pp. 193–204.
  11. T. Albash and D. A. Lidar, Adiabatic quantum computation, Rev. Mod. Phys. 90, 015002 (2018).
  12. L. Lin and Y. Tong, Near-optimal ground state preparation, Quantum 4, 372 (2020).
  13. Y. Dong, L. Lin, and Y. Tong, Ground-state preparation and energy estimation on early fault-tolerant quantum computers via quantum eigenvalue transformation of unitary matrices, PRX Quantum 3, 040305 (2022).
  14. F. Verstraete, M. M. Wolf, and J. Ignacio Cirac, Quantum computation and quantum-state engineering driven by dissipation, Nat. Phys. 5, 633 (2009).
  15. J. Cotler, S. Choi, A. Lukin, H. Gharibyan, T. Grover, M. E. Tai, M. Rispoli, R. Schittko, P. M. Preiss, A. M. Kaufman, M. Greiner, H. Pichler, and P. Hayden, Quantum virtual cooling, Phys. Rev. X 9, 031013 (2019).
  16. S. Roy, J. T. Chalker, I. V. Gornyi, and Y. Gefen, Measurement-induced steering of quantum systems, Phys. Rev. Res. 2, 033347 (2020).
  17. L. J. Schulman, T. Mor, and Y. Weinstein, Physical limits of heat-bath algorithmic cooling, Phys. Rev. Lett. 94, 120501 (2005).
  18. S. Raeisi and M. Mosca, Asymptotic bound for heat-bath algorithmic cooling, Phys. Rev. Lett. 114, 100404 (2015).
  19. H. Fan, F.-L. Wu, L. Wang, S.-Q. Liu, and S.-Y. Liu, Algorithmic cooling based on quantum measurement, Phys. Rev. A 112, 013104 (2025).
  20. T. S. Cubitt, Dissipative ground state preparation and the dissipative quantum eigensolver, arXiv:2303.11962.
  21. X. Mi, A. Michailidis, S. Shabani, K. Miao, P. Klimov, J. Lloyd, E. Rosenberg, R. Acharya, I. Aleiner, T. Andersen et al., Stable quantum-correlated many-body states through engineered dissipation, Science 383, 1332 (2024).
  22. Z. Ding, C.-F. Chen, and L. Lin, Single-ancilla ground state preparation via Lindbladians, Phys. Rev. Res. 6, 033147 (2024).
  23. 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).
  24. Z. Ding, Y. Zhan, J. Preskill, and L. Lin, End-to-end efficient quantum thermal and ground state preparation made simple,arXiv:2508.05703.
  25. Y. Zhan, Z. Ding, J. Huhn, J. Gray, J. Preskill, G. K.-L. Chan, and L. Lin, Rapid quantum ground state preparation via dissipative dynamics, Phys. Rev. X 16, 011004 (2026).
  26. Y. Suzuki, M. Gluza, J. Son, B. H. Tiang, N. H. Ng, and Z. Holmes, Grover’s algorithm is an approximation of imaginary-time evolution, arXiv:2507.15065.
  27. P. Ginsparg, Applied conformal field theory, arXiv:Hep-th/9108028.
  28. P. Calabrese and J. Cardy, Quantum quenches in extended systems, J. Stat. Mech.:Theory Exp. 2007, P06008 (2007).
  29. J. Maldacena, The large-N limit of superconformal field theories and supergravity, Int. J. Theor. Phys. 38, 1113 (1999).
  30. J. Maldacena, Eternal black holes in anti-de Sitter, J. High Energy Phys. 03 (2003) 021.
  31. J. Maldacena and L. Maoz, Wormholes in ads, J. High Energy Phys. 04 (2004) 053.
  32. V. Balasubramanian, P. Hayden, A. Maloney, D. Marolf, and S. F. Ross, Multiboundary wormholes and holographic entanglement, Classical Quantum Gravity 31, 185015 (2014).
  33. A. R. Brown, H. Gharibyan, S. Leichenauer, H. W. Lin, S. Nezami, G. Salton, L. Susskind, B. Swingle, and M. Walter, Quantum gravity in the lab. I. Teleportation by size and traversable wormholes, PRX Quantum 4, 010320 (2023).
  34. S. Nezami, H. W. Lin, A. R. Brown, H. Gharibyan, S. Leichenauer, G. Salton, L. Susskind, B. Swingle, and M. Walter, Quantum gravity in the lab. II. Teleportation by size and traversable wormholes, PRX Quantum 4, 010321 (2021).
  35. T. Schuster, B. Kobrin, P. Gao, I. Cong, E. T. Khabiboulline, N. M. Linke, M. D. Lukin, C. Monroe, B. Yoshida, and N. Y. Yao, Many-body quantum teleportation via operator spreading in the traversable wormhole protocol, Phys. Rev. X 12, 031013 (2022).
  36. 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).
  37. 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).
  38. Y. Mao, M. Chaudhary, M. Kondappan, J. Shi, E. O. Ilo-Okeke, V. Ivannikov, and T. Byrnes, Measurement-based deterministic imaginary time evolution, Phys. Rev. Lett. 131, 110602 (2023).
  39. Y. Huang, Y. Shao, W. Ren, J. Sun, and D. Lv, Efficient quantum imaginary time evolution by drifting real-time evolution: An approach with low gate and measurement complexity, J. Chem. Theory Comput. 19, 3868 (2023).
  40. K. Hejazi, M. Motta, and G. K.-L. Chan, Adiabatic quantum imaginary time evolution, Phys. Rev. Res. 6, 033084 (2024).
  41. M. Gluza, J. Son, B. H. Tiang, R. Zander, R. Seidel, Y. Suzuki, Z. Holmes, and N. H. Ng, Double-bracket quantum algorithms for quantum imaginary-time evolution, Phys. Rev. Lett. 136, 020601 (2024).
  42. R. Zander, R. Seidel, L. Xiaoyue, and M. Gluza, Role of Riemannian geometry in double-bracket quantum imaginary-time evolution, in International Conference on Geometric Science of Information (Springer, New York, 2025), pp. 105–114.
  43. M. Alghadeer, K. U. Giang, S. Cao, S. D. Fasciati, M. Piscitelli, N. Ng, P. J. Leek, M. Gluza, and M. Bakr, Double-bracket algorithmic cooling, arXiv:2510.00302.
  44. Y. Suzuki, B. H. Tiang, J. Son, N. H. Ng, Z. Holmes, and M. Gluza, Double-bracket algorithm for quantum signal processing without post-selection, Quantum 9, 1954 (2025).
  45. M. Robbiati, E. Pedicillo, A. Pasquale, X. Li, O. Kiss, A. Wright, R. Farias, K. U. Giang, J. Son, J. Knörzer et al., Double-bracket quantum algorithms for high-fidelity ground state preparation, arXiv:2408.03987.
  46. M. Gluza, Double-bracket quantum algorithms for diagonalization, Quantum 8, 1316 (2024).
  47. The Hamiltonian here is the Hermitian operator −i[ϕβ,H].

  48. S. Lloyd, M. Mohseni, and P. Rebentrost, Quantum principal component analysis, Nat. Phys. 10, 631 (2014).
  49. S. Kimmel, C. Y.-Y. Lin, G. H. Low, M. Ozols, and T. J. Yoder, Hamiltonian simulation with optimal sample complexity, npj Quantum Inf. 3, 13 (2017).
  50. M. Kjaergaard, M. Schwartz, A. Greene, G. Samach, A. Bengtsson, M. O’Keeffe, C. McNally, J. Braumüller, D. Kim, P. Krantz et al., Demonstration of density matrix exponentiation using a superconducting quantum processor, Phys. Rev. X 12, 011005 (2022).
  51. F. Wei, Z. Liu, G. Liu, Z. Han, X. Ma, D.-L. Deng, and Z. Liu, Realizing non-physical actions through Hermitian-preserving map exponentiation, arXiv:2308.07956.
  52. J. Son, M. Gluza, R. Takagi, and N. H. Y. Ng, Quantum dynamic programming, Phys. Rev. Lett. 134, 180602 (2025).
  53. Acting on copies with different entries of c→ is still a well-defined unitary operation. Such unbalanced gates may be useful as heuristic or variational additions, for example near the end of the circuit, but they are not part of the constructions analyzed here.

  54. We note that ϕβ, defined in Eq. (2), is pure, and therefore the expression for the fidelity simplifies to the trace of the product of the two states.

  55. With our sign convention, the tree circuit described here takes the first copy to the bookkeeping value −n. This choice is arbitrary: one could instead follow the opposite output branch at each layer, which would take, for example, the last copy to +n. One can also follow both branches simultaneously, as in the four-copy example of Eq. (16), so that the final bookkeeping values of the first and last copies are (−n,n). Although these two output copies are generally correlated, each individual reduced state follows the recurrence of Eq. (19), with the sign of ϵ in U12 determined by the branch being followed. No measurement or reinitialization is involved in this deterministic construction; these operations are only introduced later in the post-selected variant.

  56. N. A. McMahon, M. Pervez, and C. Arenz, Equating quantum imaginary time evolution, Riemannian gradient flows, and stochastic implementations, Phys. Rev. Res. 8, 023024 (2025).
  57. M. Pervez, A. Haqq, N. A. McMahon, and C. Arenz, Riemannian gradient descent-based quantum algorithms for ground state preparation with guarantees, arXiv:2512.13401.
  58. P. Hauke, H. G. Katzgraber, W. Lechner, H. Nishimori, and W. D. Oliver, Perspectives of quantum annealing: Methods and implementations, Rep. Prog. Phys. 83, 054401 (2020).
  59. M. Cerezo, A. Arrasmith, R. Babbush, S. C. Benjamin, S. Endo, K. Fujii, J. R. McClean, K. Mitarai, X. Yuan, L. Cincio et al., Variational quantum algorithms, Nat. Rev. Phys. 3, 625 (2021).
  60. J. Tilly, H. Chen, S. Cao, D. Picozzi, K. Setia, Y. Li, E. Grant, L. Wossnig, I. Rungger, G. H. Booth et al., The variational quantum Eigensolver: A review of methods and best practices, Phys. Rep. 986, 1 (2022).
  61. L. Zhou, S.-T. Wang, S. Choi, H. Pichler, and M. D. Lukin, Quantum approximate optimization algorithm: Performance, mechanism, and implementation on near-term devices, Phys. Rev. X 10, 021067 (2020).
  62. S. Ebadi, A. Keesling, M. Cain, T. T. Wang, H. Levine, D. Bluvstein, G. Semeghini, A. Omran, J.-G. Liu, R. Samajdar et al., Quantum optimization of maximum independent set using Rydberg atom arrays, Science 376, 1209 (2022).
  63. L. Bombieri, Z. Zeng, R. Tricarico, R. Lin, S. Notarnicola, M. Cain, M. D. Lukin, and H. Pichler, Quantum adiabatic optimization with Rydberg arrays: Localization phenomena and encoding strategies, PRX Quantum 6, 020306 (2025).
  64. The hedge circuit requires one small adjustment: within each layer, any copy not acted on by a V gate is evolved for an additional ϵ real time. This padding of the circuit equates the real time of the copies on which the V gates act, such that also the acquired phases agree and the copies are symmetric with respect to swap (at leading order).

  65. W. J. Huggins, S. McArdle, T. E. O’Brien, J. Lee, N. C. Rubin, S. Boixo, K. B. Whaley, R. Babbush, and J. R. McClean, Virtual distillation for quantum error mitigation, Phys. Rev. X 11, 041036 (2021).
  66. H. Pichler, L. Bonnes, A. J. Daley, A. M. Läuchli, and P. Zoller, Thermal versus entanglement entropy: A measurement protocol for fermionic atoms with a quantum gas microscope, New J. Phys. 15, 063003 (2013).
  67. This is the optimal gate count regardless of the circuit structure. To see this, let us represent the bookkeeping vector after the mth use of Uij as c→(m). Consider the following function S(c→(m))=(c1(m))2+(c2(m))2+⋯, where ci(m) are the elements of c→(m). As Uij only acts on two elements if they are equal, and it adds 1 to one of them and −1 to the other, we understand that S(c→(m+1))−S(c→(m))=2. Therefore, S(c→(m))=2m, and the minimal number M to reach a state is given by S(c→(m))/2. For the state c→f=(−n,…,−1,1,…,n), we obtain the claimed result.

  68. C. Gross and I. Bloch, Quantum simulations with ultracold atoms in optical lattices, Science 357, 995 (2017).
  69. R. Islam, R. Ma, P. M. Preiss, M. Eric Tai, A. Lukin, M. Rispoli, and M. Greiner, Measuring entanglement entropy in a quantum many-body system, Nature (London) 528, 77 (2015).
  70. A. M. Kaufman, M. E. Tai, A. Lukin, M. Rispoli, R. Schittko, P. M. Preiss, and M. Greiner, Quantum thermalization through entanglement in an isolated many-body system, Science 353, 794 (2016).
  71. I. Bloch, J. Dalibard, and S. Nascimbène, Quantum simulations with ultracold quantum gases, Nat. Phys. 8, 267 (2012).
  72. C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Feshbach resonances in ultracold gases, Rev. Mod. Phys. 82, 1225 (2010).
  73. H. Pichler, G. Zhu, A. Seif, P. Zoller, and M. Hafezi, Measurement protocol for the entanglement spectrum of cold atoms, Phys. Rev. X 6, 041033 (2016).
  74. D. González-Cuadra, D. Bluvstein, M. Kalinowski, R. Kaubruegger, N. Maskara, P. Naldesi, T. V. Zache, A. M. Kaufman, M. D. Lukin, H. Pichler et al., Fermionic quantum processing with programmable neutral atom arrays, Proc. Natl. Acad. Sci. U. S. A. 120, e2304294120 (2023).
  75. F. Gkritsis, D. Dux, J. Zhang, N. Jain, C. Gogolin, and P. M. Preiss, Simulating chemistry with fermionic optical superlattices, PRX Quantum 6, 010318 (2025).
  76. A. Browaeys and T. Lahaye, Many-body physics with individually controlled Rydberg atoms, Nat. Phys. 16, 132 (2020).
  77. D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner et al., A quantum processor based on coherent transport of entangled atom arrays, Nature (London) 604, 451 (2022).
  78. R. Acharya et al., Quantum error correction below the surface code threshold, Nature (London) 638, 920 (2025).
  79. R. Blatt and C. F. Roos, Quantum simulations with trapped ions, Nat. Phys. 8, 277 (2012).
  80. A. Ransford, M. Allman, J. Arkinstall, J. Campora III, S. F. Cooper, R. D. Delaney, J. M. Dreiling, B. Estey, C. Figgatt, A. Hall et al., Helios: A 98-qubit trapped-ion quantum computer, arXiv:2511.05465.
  81. T. I. Andersen, N. Astrakhantsev, A. H. Karamlou, J. Berndtsson, J. Motruk, A. Szasz, J. A. Gross, A. Schuckert, T. Westerhout, Y. Zhang et al., Thermalization and criticality on an analogue–digital quantum simulator, Nature (London) 638, 79 (2025).
  82. M. Foss-Feig, G. Pagano, A. C. Potter, and N. Y. Yao, Progress in trapped-ion quantum simulation, Annu. Rev. Condens. Matter Phys. 16, 145 (2025).
  83. We note that the sign of ϵ here is opposite to the definition of U in Eq. (11). This is because we (arbitrarily) chose to evolve the first copy from ϕ0 to ϕβ, and not to ϕ−β in the error bound derivation.

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation