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:2402.11667 · The Journal of Chemical Physics

Fast-forwarding molecular ground state preparation with optimal control on analog quantum simulators

Davide Castaldo, Marta Rosa, Stefano Corni

Published Jul 1, 2024• 0 citations
arXivHTMLPDF

Abstract

We show that optimal control of the electron dynamics is able to prepare molecular ground states, within chemical accuracy, with evolution times approaching the bounds imposed by quantum mechanics. We propose a specific parameterization of the molecular evolution only in terms of interaction already present in the molecular Hamiltonian. Thus, the proposed method solely utilizes quantum simulation routines, retaining their favorable scalings. Due to the intimate relationships between variational quantum algorithms and optimal control, we compare, when possible, our results with state-of-the-art methods in the literature. We found that the number of parameters needed to reach chemical accuracy and algorithmic scaling is in line with compact adaptive strategies to build variational Ansätze. The algorithm, which is also suitable for quantum simulators, is implemented by emulating a digital quantum processor (up to 16 qubits) and tested on different molecules and geometries spanning different degrees of electron correlation.

DOI
10.1063/5.0204618

Related papers

Qubit-Adaptive Variational Quantum Imaginary Time Evolution: A Resource-Efficient Algorithm for Ground-State Preparation of Correlated Molecules on Near-Term Quantum Computers

Thitiwuth Chaiyatho, Zhenyu Li, Jin‐Shi Xu +2 more

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.

Journal of Chemical Theory and Computation•Sep 26, 2026•0 citations

Nonvariational ADAPT algorithm for quantum simulations

Ho Lun Tang, Yanzhu Chen, Prakriti Biswas +3 more

We explore a nonvariational quantum state preparation approach combined with the ADAPT operator selection strategy in the application of preparing the ground state of a desired target Hamiltonian. In this algorithm, energy gradient measurements determine both the operators and the gate parameters in the quantum circuit construction. We compare this nonvariational algorithm with ADAPT-VQE and with feedback-based quantum algorithms in terms of the rate of energy reduction, the circuit depth, and the measurement cost in molecular simulation. We find that, despite using deeper circuits, this new algorithm reaches chemical accuracy at a similar measurement cost to ADAPT-VQE. Since it does not rely on a classical optimization subroutine, it may provide robustness against circuit parameter errors due to imperfect control or gate synthesis.

arXiv:2411.09736•Physical Review Research•May 12, 2025•7 citations

Pulse based Variational Quantum Optimal Control for hybrid quantum computing

Robert de Keijzer, Oliver T.C. Tse, S. J. J. M. F. Kokkelmans

This work studies pulse based variational quantum algorithms (VQAs), which are designed to determine the ground state of a quantum mechanical system by combining classical and quantum hardware. In contrast to more standard gate based methods, pulse based methods aim to directly optimize the laser pulses interacting with the qubits, instead of using some parametrized gate based circuit. Using the mathematical formalism of optimal control, these laser pulses are optimized. This method has been used in quantum computing to optimize pulses for quantum gate implementations, but has only recently been proposed for full optimization in VQAs. Pulse based methods have several advantages over gate based methods such as faster state preparation, simpler implementation and more freedom in moving through the state space. Based on these ideas, we present the development of a novel adjoint based variational method. This method can be tailored towards and applied in neutral atom quantum computers. This method of pulse based variational quantum optimal control is able to approximate molecular ground states of simple molecules up to chemical accuracy and is able to compete with the gate based variational quantum eigensolver in terms of total number of quantum evaluations. The total evolution time T and the form of the control Hamiltonian Hc are important factors in the convergence behavior to the ground state energy, both having influence on the quantum speed limit and the controllability of the system.

arXiv:2202.08908•Quantum•Jan 26, 2023•31 citations

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Niladri Gomes, Anirban Mukherjee, Feng Zhang +5 more

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

arXiv:2102.01544•Advanced Quantum Technologies•Oct 29, 2021•84 citations

Shortcuts for adiabatic and variational algorithms in molecular simulation

Julián Ferreiro-Vélez, Iñaki Iriarte-Zendoia, Yue Ban +1 more

Quantum algorithms offer a promising route toward computational advantage, but current implementations remain constrained by limited coherence and gate errors, making circuit complexity and noise sensitivity critical considerations. Here, we develop a shortcuts-to-adiabaticity strategy for molecular ground-state preparation and apply it to both digitized adiabatic quantum computing and the variational quantum eigensolver. In the adiabatic framework, we construct an approximate counterdiabatic (CD) Hamiltonian using a nested-commutator expansion. For the molecular systems considered, the CD correction accelerates state preparation and improves energy convergence in the fast-evolution regime, enabling a target accuracy to be reached with a lower implemented-operator count. In the variational setting, we introduce the adiabatic gauge ansatz ($\mbox{AGA}$), obtained by promoting the Pauli-string components generated by the approximate adiabatic gauge potential (AGP) to independent variational generators. We further propose a reduced variant, $\mbox{AGAR}$, that retains only Pauli generators with weight not larger than two, substantially reducing the ansatz size while maintaining competitive accuracy over the regimes considered. Benchmarks on $\mathrm{LiH}$ and $\mathrm{BeH}_2$ show that the proposed ansätze achieve accuracies competitive with established approaches such as UCCSD. Experiments on the \texttt{IBM\_basquecountry} superconducting quantum processor, together with simulations under depolarizing noise, show that the reduced CD-inspired ansatz yields smaller energy discrepancies than UCCSD under the noise conditions considered. Overall, our results demonstrate how AGP-inspired operator structures can be used to balance accuracy and implementation cost in molecular simulations on current noise-limited quantum hardware.

arXiv:2407.20957•Physical Review A•Oct 7, 2026•2 citations

An adaptive variational algorithm for exact molecular simulations on a quantum computer

Harper R. Grimsley, Sophia E. Economou, Edwin Barnes +1 more

Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.

arXiv:1812.11173•Nature Communications•Jul 8, 2019•1,029 citations
We show that optimal control of the electron dynamics is able to prepare molecular ground states, within chemical accuracy, with evolution times approaching the bounds imposed by quantum mechanics. We propose a specific parameterization of the molecular evolution only in terms of interaction already present in the molecular Hamiltonian. Thus, the proposed method solely utilizes quantum simulation routines, retaining their favorable scalings. Due to the intimate relationships between variational quantum algorithms and optimal control, we compare, when possible, our results with state-of-the-art methods in the literature. We found that the number of parameters needed to reach chemical accuracy and algorithmic scaling is in line with compact adaptive strategies to build variational Ansätze. The algorithm, which is also suitable for quantum simulators, is implemented by emulating a digital quantum processor (up to 16 qubits) and tested on different molecules and geometries spanning different degrees of electron correlation.
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.
We explore a nonvariational quantum state preparation approach combined with the ADAPT operator selection strategy in the application of preparing the ground state of a desired target Hamiltonian. In this algorithm, energy gradient measurements determine both the operators and the gate parameters in the quantum circuit construction. We compare this nonvariational algorithm with ADAPT-VQE and with feedback-based quantum algorithms in terms of the rate of energy reduction, the circuit depth, and the measurement cost in molecular simulation. We find that, despite using deeper circuits, this new algorithm reaches chemical accuracy at a similar measurement cost to ADAPT-VQE. Since it does not rely on a classical optimization subroutine, it may provide robustness against circuit parameter errors due to imperfect control or gate synthesis.
This work studies pulse based variational quantum algorithms (VQAs), which are designed to determine the ground state of a quantum mechanical system by combining classical and quantum hardware. In contrast to more standard gate based methods, pulse based methods aim to directly optimize the laser pulses interacting with the qubits, instead of using some parametrized gate based circuit. Using the mathematical formalism of optimal control, these laser pulses are optimized. This method has been used in quantum computing to optimize pulses for quantum gate implementations, but has only recently been proposed for full optimization in VQAs. Pulse based methods have several advantages over gate based methods such as faster state preparation, simpler implementation and more freedom in moving through the state space. Based on these ideas, we present the development of a novel adjoint based variational method. This method can be tailored towards and applied in neutral atom quantum computers. This method of pulse based variational quantum optimal control is able to approximate molecular ground states of simple molecules up to chemical accuracy and is able to compete with the gate based variational quantum eigensolver in terms of total number of quantum evaluations. The total evolution time T and the form of the control Hamiltonian Hc are important factors in the convergence behavior to the ground state energy, both having influence on the quantum speed limit and the controllability of the system.
Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.
Quantum algorithms offer a promising route toward computational advantage, but current implementations remain constrained by limited coherence and gate errors, making circuit complexity and noise sensitivity critical considerations. Here, we develop a shortcuts-to-adiabaticity strategy for molecular ground-state preparation and apply it to both digitized adiabatic quantum computing and the variational quantum eigensolver. In the adiabatic framework, we construct an approximate counterdiabatic (CD) Hamiltonian using a nested-commutator expansion. For the molecular systems considered, the CD correction accelerates state preparation and improves energy convergence in the fast-evolution regime, enabling a target accuracy to be reached with a lower implemented-operator count. In the variational setting, we introduce the adiabatic gauge ansatz ($\mbox{AGA}$), obtained by promoting the Pauli-string components generated by the approximate adiabatic gauge potential (AGP) to independent variational generators. We further propose a reduced variant, $\mbox{AGAR}$, that retains only Pauli generators with weight not larger than two, substantially reducing the ansatz size while maintaining competitive accuracy over the regimes considered. Benchmarks on LiH\mathrm{LiH}LiH and BeH2\mathrm{BeH}_2BeH2​ show that the proposed ansätze achieve accuracies competitive with established approaches such as UCCSD. Experiments on the \texttt{IBM\_basquecountry} superconducting quantum processor, together with simulations under depolarizing noise, show that the reduced CD-inspired ansatz yields smaller energy discrepancies than UCCSD under the noise conditions considered. Overall, our results demonstrate how AGP-inspired operator structures can be used to balance accuracy and implementation cost in molecular simulations on current noise-limited quantum hardware.
Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.