Abstract Many current and near-future applications of quantum computing utilise parametric families of quantum circuits and variational methods that can suffer from obstacles including non-convergence to the global minimum due to local minima, critical slowing down, or exponential resource scaling. Here we show that quantum imaginary time evolution can overcome these obstacles if the underlying physical system satisfies a set of conditions. This includes many relevant applications such as ground state preparation for local theories in physics or chemistry, combinatorial optimisation problems, or quantum machine learning. In particular, we analyse the quantum imaginary time evolution showing convergence guarantees to the global minimum without critical slowing down and providing a priori estimates on the required evolution time which scale linearly in system size and inverse energy gap. Furthermore, a provided complexity analysis shows that quantum imaginary time evolution can be efficiently compiled into a parametric quantum circuit, finding the optimal parameters included, for a large class of physically relevant problems.
Abstract Many current and near-future applications of quantum computing utilise parametric families of quantum circuits and variational methods that can suffer from obstacles including non-convergence to the global minimum due to local minima, critical slowing down, or exponential resource scaling. Here we show that quantum imaginary time evolution can overcome these obstacles if the underlying physical system satisfies a set of conditions. This includes many relevant applications such as ground state preparation for local theories in physics or chemistry, combinatorial optimisation problems, or quantum machine learning. In particular, we analyse the quantum imaginary time evolution showing convergence guarantees to the global minimum without critical slowing down and providing a priori estimates on the required evolution time which scale linearly in system size and inverse energy gap. Furthermore, a provided complexity analysis shows that quantum imaginary time evolution can be efficiently compiled into a parametric quantum circuit, finding the optimal parameters included, for a large class of physically relevant problems.
The imaginary-time evolution of quantum states is integral to various fields, ranging from natural sciences to classical optimization or machine learning. Since simulating quantum imaginary-time evolution generally requires storing an exponentially large wave function, quantum computers are emerging as a promising platform for this task. However, variational approaches, suitable for near-term quantum computers, struggle with a prohibitive number of measurements and impractical runtimes for relevant system sizes. Here, we suggest a stochastic approach to variational quantum imaginary-time evolution, which allows a significant reduction in runtimes. Our approach allows trading off invested resources and accuracy, which makes it also suitable for ground state preparation, where simulating the exact dynamics is not required. We demonstrate the efficiency of our algorithm in simulations and show a hardware experiment performing the imaginary-time evolution of the transverse field Ising model on 27 qubits.
Abstract Imaginary-time evolution is fundamental for analyzing quantum many-body systems, with applications spanning quantum chemistry, condensed matter physics, and quantum field theory, yet classical simulation requires exponentially growing resources in both system size and evolution time. While quantum approaches reduce the system-size scaling, existing methods rely on heuristic techniques with measurement precision or success probability that deteriorates as evolution time increases. We present a quantum algorithm that prepares normalized imaginary-time evolved states using an adaptive normalization factor to maintain a stable success probability over long imaginary-time intervals. Our algorithm approximates the target state with error polynomially small in the inverse imaginary time using a polynomial number of elementary quantum gates and a single ancilla qubit, with success probability close to one. When the initial state has reasonable overlap with the ground state, this algorithm also achieves polynomial query complexity in the system size. To our knowledge, this is the first quantum algorithm for imaginary-time evolution with provably polynomial resource scaling in evolution time. Numerical experiments validate our theoretical analysis for evolution time up to 50, demonstrating the algorithm’s effectiveness for long-time evolution. Building on this technique, we further develop imaginary-time-evolution-based algorithms for ground-state-related problems and for simulating open quantum systems. These algorithms can reduce circuit depth in certain regimes compared with existing methods, at the expense of higher total query complexity, advancing the practical feasibility of quantum simulation on early fault-tolerant devices.
Preparing ground states of correlated molecules on near-term quantum computers requires algorithms that are both accurate and resource-efficient. We introduce the qubit-adaptive variational quantum imaginary time evolution (QAITE) algorithm, which synthesizes the systematic convergence of imaginary-time evolution with hardware-efficient, adaptive ansatz construction. By integrating the linear-scaling qubit operator pool into the adaptive variational quantum imaginary time evolution framework, QAITE deterministically constructs a compact ansatz without costly classical optimization loops. Numerical simulations for molecular systems (H4, LiH, H2O, BeH2), benchmarked against our own implementations of established adaptive algorithms under identical conditions, show that QAITE achieves chemical accuracy with a 15-45% reduction in CNOT gates compared to the standard adaptive variational quantum imaginary time evolution (AVQITE) method. For hydrogen chains (RHH = 1.5 Å, N = 6-16 qubits), QAITE's circuit complexity grows approximately quadratically and remains the lowest of the four methods at every system size, reaching 31% fewer CNOT gates than AVQITE at 16 qubits─a saving of 1322 gates. This efficiency extends to heterogeneous (LiH)n and (HBe)n chains, where QAITE maintains a CNOT gate ratio below 1 for all sizes studied (e.g., 0.81 for (LiH)2). We further show that the deterministic parameter update is numerically robust: the results are insensitive to the regularization of the quantum geometric tensor over at least 4 orders of magnitude, and chemical accuracy is retained under realistic shot noise applied to both the quantum geometric tensor and the force vector. By outperforming established adaptive algorithms in both absolute gate count and scaling, QAITE establishes a promising, scalable pathway for quantum molecular simulations on noisy intermediate-scale quantum devices.
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.
Abstract Imaginary time evolution is a powerful tool for studying quantum systems. While it is possible to simulate with a classical computer, the time and memory requirements generally scale exponentially with the system size. Conversely, quantum computers can efficiently simulate quantum systems, but not non-unitary imaginary time evolution. We propose a variational algorithm for simulating imaginary time evolution on a hybrid quantum computer. We use this algorithm to find the ground-state energy of many-particle systems; specifically molecular hydrogen and lithium hydride, finding the ground state with high probability. Our method can also be applied to general optimisation problems and quantum machine learning. As our algorithm is hybrid, suitable for error mitigation and can exploit shallow quantum circuits, it can be implemented with current quantum computers.
Variational compression can significantly lower implementation overheads for encoding the time evolution of Hamiltonians into quantum circuits. However, they usually lack global convergence guarantees and well-established scaling behavior. In this work, we provide a recipe for choosing the initial point of such variational optimizations that guarantees convergence to a quantum circuit with near-optimal gate complexity O(Ntpolylog(Nt/ε)) for all local and translationally invariant Hamiltonians. We demonstrate our method by encoding the globally controlled time evolution of a Heisenberg antiferromagnet on a Kagome lattice. For N=48 sites, evolution time t=0.1 and infidelity ε≈1%, the controlled time-evolution circuit requires 960 two-qubit B gates, for which we propose a straightforward implementation scheme for ion-trap setups. Thereby, our recipe extends digital quantum simulators toward system sizes and geometries that are challenging for classical computation.