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

Analytical Nuclear Gradients and Hessians on Quantum Hardware via Orbital-Optimized VQE with Error Mitigation

Renato Olarte Hernandez, Karl Michael Ziems, Erik Kjellgren, Jacob Kongsted, Sonia Coriani, Stephan P. A. Sauer

physics.chem-phquant-phSubmitted Aug 9, 2026• Announced Aug 11, 2026• 0 citations
arXivHTMLPDF

Abstract

Nuclear gradients and Hessians are fundamental quantities in computational chemistry, essential for a wide range of applications including geometry optimization, vibrational spectroscopy, and molecular property calculations. In this work, we present their analytical implementation on quantum hardware. The methodology is formulated within an active-space framework combining orbital optimization and linear-response theory. On the quantum-computing side, the approach employs the tiled unitary product state (tUPS) ansatz to directly evaluate the tensor elements required for solving the response equations. Moreover, the expectation values are corrected using an adapted confusion-matrix error-mitigation scheme in combination with post-selection criteria. The resulting workflow is assessed on molecular hydrogen and on water through the calculation of potential energy surfaces, nuclear gradients, Hessians, and vibrational frequencies, enabling the evaluation of both its capabilities and current limitations. The results demonstrate good performance for the hydrogen molecule, whereas the water molecule provides a more demanding test of quantum-hardware resources and highlights the trade-offs associated with error-mitigation strategies. The quantified analysis of the results identify the main sources of errors, suggesting improvement directions for more accurate quantum computer applications.

Comments
37 pages, 13 figures

Related papers

Analytical Ground- and Excited-State Gradients for Molecular Electronic Structure Theory from Hybrid Quantum/Classical Methods

Robert M. Parrish, Gian-Luca Anselmetti, Christian Gogolin

We develop analytical gradients of ground- and excited-state energies with respect to system parameters including the nuclear coordinates for the hybrid quantum/classical multistate contracted variational quantum eigensolver (MC-VQE) applied to fermionic systems. We show how the resulting response contributions to the gradient can be evaluated with a quantum effort similar to that of obtaining the VQE energy and independent of the total number of derivative parameters (e.g. number of nuclear coordinates) by adopting a Lagrangian formalism for the evaluation of the total derivative. We also demonstrate that large-step-size finite-difference treatment of directional derivatives in concert with the parameter shift rule can significantly mitigate the complexity of dealing with the quantum parameter Hessian when solving the quantum response equations. This enables the computation of analytical derivative properties of systems with hundreds of atoms, while solving an active space of their most strongly correlated orbitals on a quantum computer. We numerically demonstrate the exactness the analytical gradients and discuss the magnitude of the quantum response contributions.

arXiv:2110.05040•Oct 11, 2021•5 citations

The impacts of optimization algorithm and basis size on the accuracy and efficiency of variational quantum eigensolver

Xian‐Hu Zha, Zhang, Chao, Dengdong Fan +4 more

Variational quantum eigensolver (VQE) is demonstrated to be the promising methodology for quantum chemistry based on near-term quantum devices. However, many problems are yet to be investigated for this methodology, such as the influences of optimization algorithm and basis size on the accuracy and efficiency for quantum computing. To address these issues, five molecules (H2, LiH, HF, N2 and F2) are studied in this work based on the VQE method using unitary coupled cluster (UCC) ansatz. The performance of the gradient optimization L-BFGS-B is compared with that of the direct search method COBYLA. The former converges more quickly, but the accuracy of energy surface is a little lower. The basis set shows a vital influence on the accuracy and efficiency. A large basis set generally provides an accurate energy surface, but induces a significant increase in computing time. The 631g basis is generally required from the energy surface of the simplest H2 molecule. For practical applications of VQE, complete active space (CAS) is suggested based on limited quantum resources. With the same number of qubits, more occupied orbitals included in CAS gives a better accuracy for the energy surface and a smaller evaluation number in the VQE optimization. Additionally, the electronic structure, such as filling fraction of orbitals, the bond strength of a molecule and the maximum nuclear charge also influences the performance of optimization, where half occupation of orbitals generally requires a large computation cost.

arXiv:2006.15852•Jun 29, 2020•3 citations

The impacts of optimization algorithm and basis size on the accuracy and efficiency of variational quantum eigensolver based on unitary coupled cluster ansatz

Xian‐Hu Zha, C. Zhang, Dengdong Fan +3 more

Variational quantum eigensolver (VQE) is demonstrated to be the promising methodology for quantum chemistry based on near-term quantum devices. However, many problems are yet to be investigated for this methodology, such as the influences of optimization algorithm and basis size on the accuracy and efficiency for quantum computing. To address these issues, five molecules (H2, LiH, HF, N2 and F2) are studied in this work based on the VQE method using unitary coupled cluster (UCC) ansatz. The performance of the gradient optimization L-BFGS-B is compared with that of the direct search method COBYLA. The former converges more quickly, but the accuracy of energy surface is a little lower. The basis set shows a vital influence on the accuracy and efficiency. A large basis set generally provides an accurate energy surface, but induces a significant increase in computing time. The 631g basis is generally required from the energy surface of the simplest H2 molecule. For practical applications of VQE, complete active space (CAS) is suggested based on limited quantum resources. With the same number of qubits, more occupied orbitals included in CAS gives a better accuracy for the energy surface and a smaller evaluation number in the VQE optimization. Additionally, the electronic structure, such as filling fraction of orbitals, the bond strength of a molecule and the maximum nuclear charge also influences the performance of optimization, where half occupation of orbitals generally requires a large computation cost.

•Jun 29, 2020•0 citations

Implementing QESEM's High-Accuracy Error Mitigation on a Quantum Computer: a Water Potential Energy Surface Study

Renato Olarte Hernandez, Emanuele Rossi, Sonia Coriani +5 more

Quantum error mitigation (QEM) is essential for extracting chemically accurate results from near-term quantum hardware. Many widely used QEM methods rely on uncontrolled heuristics whose bias depends on the specific circuit and noise realization. In this work, we employ QESEM---a characterization-based, unbiased quasi-probabilistic mitigation method---on IBM's Aachen quantum processor to compute the ground-state potential energy surface (PES) of the symmetrically-stretched water molecule. We consider a classically-optimized, single-layer perfect-pairing tiled unitary product state ansatz. We map this ansatz to an 8-qubit register corresponding to a (4,4) active space and the STO-3G basis set. Compared to the statevector reference, we find that raw QPU results typically overestimate the ground-state energy by around 500~mHa across the scanned geometries. On the other hand, QESEM-mitigated results fall approximately within 100~mHa, 30~mHa, or ``chemical accuracy'' ($\sim$1.5~mHa), depending on the target precision. We benchmark QESEM at both loose (0.1~Ha) and tight (0.01~Ha) precision targets, evaluating both individual and merged batches of runs. As the precision target is tightened, the accuracy improves systematically, matching or exceeding results reported in the literature. We further quantify the sampling cost of these results, reporting the number of shots required at each precision level. Our results show how, with the current levels of hardware error, reaching the highest accuracies demands substantial QPU time. The results demonstrate that the characterization-based, unbiased error mitigation provided by QESEM allows to measure quantitatively meaningful potential energy surfaces on current quantum hardware. Concurrently, they highlight the sampling overhead associated with QEM, which remains a central bottleneck en route to larger chemical problems and higher precision.

quant-pharXiv:2609.07284•Sep 8, 2026•0 citations

Accurate Vibrational Frequency Calculations for Quantum Computing via an Analytic Second-Order Energy Derivative Framework

Juntao Lai, Qiang Fu, Zhenyu Li +1 more

Quantum computing has emerged as a promising paradigm for tackling electronic structure problems, with most efforts to date focused on molecular energies and, more recently, first-order derivatives. However, the extension to second-order energy derivatives with respect to nuclear coordinates─essential for predicting vibrational spectra and identifying transition states─has remained relatively limited. Here, we present an analytic implementation for computing nuclear Hessians within the variational quantum eigensolver framework. Our approach produces harmonic vibrational frequencies and normal modes in excellent agreement with full configuration interaction (FCI) benchmarks, even for systems with challenging cases of orbital degeneracy or involving weak intermolecular interactions. We further assess the quantum measurement cost of these second-order derivative calculations and compare it with that of variational quantum eigensolver (VQE) energy calculations. Additionally, we demonstrate how point group symmetry can be incorporated to reduce measurement cost without loss of accuracy. This work extends quantum computing capabilities toward advanced quantum chemical simulations, such as the characterization of noncovalent interactions via low-frequency vibrational signatures and the identification of transition states.

Journal of Chemical Theory and Computation•Apr 17, 2026•0 citations

Limitations of quantum hardware for molecular energy estimation using VQE

Abel Carreras, Román Orús, David Casanova

Variational quantum eigensolvers (VQEs) are among the most promising quantum algorithms for solving electronic structure problems in quantum chemistry, particularly during the noisy intermediate-scale quantum (NISQ) era. In this study, we investigate the capabilities and limitations of VQE algorithms implemented on current quantum hardware for determining molecular ground-state energies, focusing on the adaptive derivative-assembled pseudo-Trotter ansatz VQE (ADAPT-VQE). To address the significant computational challenges posed by molecular Hamiltonians, we explore various well known strategies to simplify the Hamiltonian, optimize the ansatz, and improve classical parameter optimization through modifications of the COBYLA optimizer. These enhancements are integrated into a tailored quantum computing implementation designed to minimize the circuit depth and computational cost. Using benzene as a benchmark system, we demonstrate the application of these optimizations on an IBM quantum computer. Despite these improvements, our results highlight the limitations imposed by current quantum hardware, particularly the impact of quantum noise on state preparation and energy measurement. The noise levels in today's devices prevent meaningful evaluations of molecular Hamiltonians with sufficient accuracy to produce reliable quantum chemical insights. Finally, we extrapolate the requirements for future quantum hardware to enable practical and scalable quantum chemistry calculations using VQE algorithms. This work provides an assessment of current quantum algorithms for molecular modeling on real quantum hardware, highlighting the impact of noise and hardware limitations on the achievable accuracy.

arXiv:2506.03995•Physical Chemistry Chemical Physics•Jan 1, 2026•5 citations
Nuclear gradients and Hessians are fundamental quantities in computational chemistry, essential for a wide range of applications including geometry optimization, vibrational spectroscopy, and molecular property calculations. In this work, we present their analytical implementation on quantum hardware. The methodology is formulated within an active-space framework combining orbital optimization and linear-response theory.
On the quantum-computing side, the approach employs the tiled unitary product state (tUPS) ansatz to directly evaluate the tensor elements required for solving the response equations. Moreover, the expectation values are corrected using an adapted confusion-matrix error-mitigation scheme in combination with post-selection criteria. The resulting workflow is assessed on molecular hydrogen and on water through the calculation of potential energy surfaces, nuclear gradients, Hessians, and vibrational frequencies, enabling the evaluation of both its capabilities and current limitations.
The results demonstrate good performance for the hydrogen molecule, whereas the water molecule provides a more demanding test of quantum-hardware resources and highlights the trade-offs associated with error-mitigation strategies. The quantified analysis of the results identify the main sources of errors, suggesting improvement directions for more accurate quantum computer applications.
We develop analytical gradients of ground- and excited-state energies with respect to system parameters including the nuclear coordinates for the hybrid quantum/classical multistate contracted variational quantum eigensolver (MC-VQE) applied to fermionic systems. We show how the resulting response contributions to the gradient can be evaluated with a quantum effort similar to that of obtaining the VQE energy and independent of the total number of derivative parameters (e.g. number of nuclear coordinates) by adopting a Lagrangian formalism for the evaluation of the total derivative. We also demonstrate that large-step-size finite-difference treatment of directional derivatives in concert with the parameter shift rule can significantly mitigate the complexity of dealing with the quantum parameter Hessian when solving the quantum response equations. This enables the computation of analytical derivative properties of systems with hundreds of atoms, while solving an active space of their most strongly correlated orbitals on a quantum computer. We numerically demonstrate the exactness the analytical gradients and discuss the magnitude of the quantum response contributions.
Variational quantum eigensolver (VQE) is demonstrated to be the promising methodology for quantum chemistry based on near-term quantum devices. However, many problems are yet to be investigated for this methodology, such as the influences of optimization algorithm and basis size on the accuracy and efficiency for quantum computing. To address these issues, five molecules (H2, LiH, HF, N2 and F2) are studied in this work based on the VQE method using unitary coupled cluster (UCC) ansatz. The performance of the gradient optimization L-BFGS-B is compared with that of the direct search method COBYLA. The former converges more quickly, but the accuracy of energy surface is a little lower. The basis set shows a vital influence on the accuracy and efficiency. A large basis set generally provides an accurate energy surface, but induces a significant increase in computing time. The 631g basis is generally required from the energy surface of the simplest H2 molecule. For practical applications of VQE, complete active space (CAS) is suggested based on limited quantum resources. With the same number of qubits, more occupied orbitals included in CAS gives a better accuracy for the energy surface and a smaller evaluation number in the VQE optimization. Additionally, the electronic structure, such as filling fraction of orbitals, the bond strength of a molecule and the maximum nuclear charge also influences the performance of optimization, where half occupation of orbitals generally requires a large computation cost.
Variational quantum eigensolver (VQE) is demonstrated to be the promising methodology for quantum chemistry based on near-term quantum devices. However, many problems are yet to be investigated for this methodology, such as the influences of optimization algorithm and basis size on the accuracy and efficiency for quantum computing. To address these issues, five molecules (H2, LiH, HF, N2 and F2) are studied in this work based on the VQE method using unitary coupled cluster (UCC) ansatz. The performance of the gradient optimization L-BFGS-B is compared with that of the direct search method COBYLA. The former converges more quickly, but the accuracy of energy surface is a little lower. The basis set shows a vital influence on the accuracy and efficiency. A large basis set generally provides an accurate energy surface, but induces a significant increase in computing time. The 631g basis is generally required from the energy surface of the simplest H2 molecule. For practical applications of VQE, complete active space (CAS) is suggested based on limited quantum resources. With the same number of qubits, more occupied orbitals included in CAS gives a better accuracy for the energy surface and a smaller evaluation number in the VQE optimization. Additionally, the electronic structure, such as filling fraction of orbitals, the bond strength of a molecule and the maximum nuclear charge also influences the performance of optimization, where half occupation of orbitals generally requires a large computation cost.
Quantum error mitigation (QEM) is essential for extracting chemically accurate results from near-term quantum hardware. Many widely used QEM methods rely on uncontrolled heuristics whose bias depends on the specific circuit and noise realization. In this work, we employ QESEM---a characterization-based, unbiased quasi-probabilistic mitigation method---on IBM's Aachen quantum processor to compute the ground-state potential energy surface (PES) of the symmetrically-stretched water molecule. We consider a classically-optimized, single-layer perfect-pairing tiled unitary product state ansatz. We map this ansatz to an 8-qubit register corresponding to a (4,4) active space and the STO-3G basis set. Compared to the statevector reference, we find that raw QPU results typically overestimate the ground-state energy by around 500~mHa across the scanned geometries. On the other hand, QESEM-mitigated results fall approximately within 100~mHa, 30~mHa, or ``chemical accuracy'' ($\sim$1.5~mHa), depending on the target precision. We benchmark QESEM at both loose (0.1~Ha) and tight (0.01~Ha) precision targets, evaluating both individual and merged batches of runs. As the precision target is tightened, the accuracy improves systematically, matching or exceeding results reported in the literature. We further quantify the sampling cost of these results, reporting the number of shots required at each precision level. Our results show how, with the current levels of hardware error, reaching the highest accuracies demands substantial QPU time. The results demonstrate that the characterization-based, unbiased error mitigation provided by QESEM allows to measure quantitatively meaningful potential energy surfaces on current quantum hardware. Concurrently, they highlight the sampling overhead associated with QEM, which remains a central bottleneck en route to larger chemical problems and higher precision.
Quantum computing has emerged as a promising paradigm for tackling electronic structure problems, with most efforts to date focused on molecular energies and, more recently, first-order derivatives. However, the extension to second-order energy derivatives with respect to nuclear coordinates─essential for predicting vibrational spectra and identifying transition states─has remained relatively limited. Here, we present an analytic implementation for computing nuclear Hessians within the variational quantum eigensolver framework. Our approach produces harmonic vibrational frequencies and normal modes in excellent agreement with full configuration interaction (FCI) benchmarks, even for systems with challenging cases of orbital degeneracy or involving weak intermolecular interactions. We further assess the quantum measurement cost of these second-order derivative calculations and compare it with that of variational quantum eigensolver (VQE) energy calculations. Additionally, we demonstrate how point group symmetry can be incorporated to reduce measurement cost without loss of accuracy. This work extends quantum computing capabilities toward advanced quantum chemical simulations, such as the characterization of noncovalent interactions via low-frequency vibrational signatures and the identification of transition states.
Variational quantum eigensolvers (VQEs) are among the most promising quantum algorithms for solving electronic structure problems in quantum chemistry, particularly during the noisy intermediate-scale quantum (NISQ) era. In this study, we investigate the capabilities and limitations of VQE algorithms implemented on current quantum hardware for determining molecular ground-state energies, focusing on the adaptive derivative-assembled pseudo-Trotter ansatz VQE (ADAPT-VQE). To address the significant computational challenges posed by molecular Hamiltonians, we explore various well known strategies to simplify the Hamiltonian, optimize the ansatz, and improve classical parameter optimization through modifications of the COBYLA optimizer. These enhancements are integrated into a tailored quantum computing implementation designed to minimize the circuit depth and computational cost. Using benzene as a benchmark system, we demonstrate the application of these optimizations on an IBM quantum computer. Despite these improvements, our results highlight the limitations imposed by current quantum hardware, particularly the impact of quantum noise on state preparation and energy measurement. The noise levels in today's devices prevent meaningful evaluations of molecular Hamiltonians with sufficient accuracy to produce reliable quantum chemical insights. Finally, we extrapolate the requirements for future quantum hardware to enable practical and scalable quantum chemistry calculations using VQE algorithms. This work provides an assessment of current quantum algorithms for molecular modeling on real quantum hardware, highlighting the impact of noise and hardware limitations on the achievable accuracy.