Quantum computers have been widely speculated to offer significant advantages in obtaining the ground state of difficult Hamiltonian in chemistry and physics. In this work, we first propose a Lyapunov control-inspired strategy to accelerate the well-established imaginary-time method for ground-state preparation. We also dig for the source of acceleration of the imaginary-time process under Lyapunov control with theoretical understanding and dynamic process visualization. To make the method accessible in the noisy intermediate-scale quantum era, we further propose a variational form of the algorithm that could work with shallow quantum circuits. Through numerical experiments on a broad spectrum of realistic models, including molecular systems, 2D Heisenberg models, and Sherrington-Kirkpatrick models, we show that imaginary-time control may substantially accelerate the imaginary time evolution for all systems and even generate orders of magnitude acceleration (suggesting exponential-like acceleration) for challenging molecular Hamiltonians involving small energy gaps as impressive special cases. Finally, with a proper selection of the control Hamiltonian, the new variational quantum algorithm does not incur additional measurement costs compared to the original variational quantum imaginary-time algorithm.
Quantum computers have been widely speculated to offer significant advantages in obtaining the ground state of difficult Hamiltonian in chemistry and physics. In this work, we first propose a Lyapunov control-inspired strategy to accelerate the well-established imaginary-time method for ground-state preparation. We also dig for the source of acceleration of the imaginary-time process under Lyapunov control with theoretical understanding and dynamic process visualization. To make the method accessible in the noisy intermediate-scale quantum era, we further propose a variational form of the algorithm that could work with shallow quantum circuits. Through numerical experiments on a broad spectrum of realistic models, including molecular systems, 2D Heisenberg models, and Sherrington-Kirkpatrick models, we show that imaginary-time control may substantially accelerate the imaginary time evolution for all systems and even generate orders of magnitude acceleration (suggesting exponential-like acceleration) for challenging molecular Hamiltonians involving small energy gaps as impressive special cases. Finally, with a proper selection of the control Hamiltonian, the new variational quantum algorithm does not incur additional measurement costs compared to the original variational quantum imaginary-time algorithm.
Quantum computers have been widely speculated to offer significant advantages in obtaining the ground state of difficult Hamiltonian in chemistry and physics. In this work, we first propose a Lyapunov control-inspired strategy to accelerate the well-established imaginary-time method for ground-state preparation. We also dig for the source of acceleration of the imaginary-time process under Lyapunov control with theoretical understanding and dynamic process visualization. To make the method accessible in the noisy intermediate-scale quantum era, we further propose a variational form of the algorithm that could work with shallow quantum circuits. Through numerical experiments on a broad spectrum of realistic models, including molecular systems, 2D Heisenberg models, and Sherrington-Kirkpatrick models, we show that imaginary-time control may substantially accelerate the imaginary-time evolution for all systems and even generate orders of magnitude acceleration (suggesting order-of-magnitude acceleration) for challenging molecular Hamiltonians involving small energy gaps as impressive special cases. Finally, with a proper selection of the control Hamiltonian, the new variational quantum algorithm does not incur additional measurement costs compared to the original variational quantum imaginary-time algorithm.
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.
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.
We introduce an adiabatic state preparation protocol which implements quantum imaginary time evolution under the Hamiltonian of the system. Unlike the original quantum imaginary time evolution algorithm, adiabatic quantum imaginary time evolution does not require quantum state tomography during its runtime and, unlike standard adiabatic state preparation, the final Hamiltonian is not the system Hamiltonian. Instead, the algorithm obtains the adiabatic Hamiltonian by integrating a classical differential equation that ensures that one follows the imaginary time evolution state trajectory. We introduce some heuristics that allow this protocol to be implemented on quantum architectures with limited resources. We explore the performance of this algorithm via classical simulations in a one-dimensional spin model and highlight essential features that determine its cost, performance, and implementability for longer times, and compare to the original quantum imaginary time evolution for ground-state preparation. More generally, our algorithm expands the range of states accessible to adiabatic state preparation methods beyond those that are expressed as ground states of simple explicit Hamiltonians. Published by the American Physical Society 2024
Ground-state preparation for a given Hamiltonian is a common quantum-computing task of great importance and has relevant applications in quantum chemistry, computational material modeling, and combinatorial optimization. We consider an approach to simulate dissipative non-Hermitian Hamiltonian quantum dynamics using Hamiltonian simulation techniques to efficiently recover the ground state of a target Hamiltonian. The proposed method facilitates the energy transfer by repeatedly projecting ancilla qubits to the desired state, rendering the effective non-Hermitian Hamiltonian evolution on the system qubits. To make the method more resource friendly in the noisy intermediate-scale quantum (NISQ) and early fault-tolerant era, we combine the non-Hermitian projection algorithm with multiple variational gadgets, including variational module enhancement and variational state recording, to reduce the required circuit depth and avoid the exponentially vanishing success probability for post-selections. We compare our method, the non-Hermitian-variational algorithm, with a pure variational method -- QAOA for solving the 3-SAT problem and preparing the ground state for the transverse-field Ising model. As demonstrated by numerical evidence, the non-Hermitian-variational algorithm outperforms QAOA in convergence speed with improved quantum resource efficiency.
Being able to predict molecular properties and interactions is of utmost interest for academia as well as industry. But the vast complexity of strongly correlated molecular systems limits the performance of classical algorithms. In contrast, quantum computation has the potential to be a game changer in the field of molecular simulations. Despite the hope in quantum computation, the capabilities of current quantum computers are still insufficient for handling molecular systems of interest. In this paper, we propose a variational ansatz for today's noisy quantum computers to calculate the ground state with the help of imaginary time evolution. Although the imaginary time evolution operator is not unitary, it can be implemented on a quantum computer by a linear decomposition and subsequent Taylor series expansion. This has the advantage that only a set of shallow circuits needs to be computed on a quantum computer. The parallel nature of this algorithm can be exploited to speed-up simulations even further, if a privileged access to quantum computers is granted.