© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2211.03343 · Journal of Chemical Theory and Computation

Analytical Formulation of the Second-Order Derivative of Energy for the Orbital-Optimized Variational Quantum Eigensolver: Application to Polarizability

Yuya O. Nakagawa, Jiabao Chen, Shotaro Sudo, Yu‐ya Ohnishi, Wataru Mizukami

Published Mar 28, 2023• 8 citations

Abstract

We develop a quantum-classical hybrid algorithm to calculate the analytical second-order derivative of the energy for the orbital-optimized variational quantum eigensolver (OO-VQE), which is a method to calculate eigenenergies of a given molecular Hamiltonian by utilizing near-term quantum computers and classical computers. We show that all quantities required in the algorithm to calculate the derivative can be evaluated on quantum computers as standard quantum expectation values without using any ancillary qubits. We validate our formula by numerical simulations of quantum circuits for computing the polarizability of the water molecule, which is the second-order derivative of the energy, with respect to the electric field. Moreover, the polarizabilities and refractive indices of thiophene and furan molecules are calculated as a test bed for possible industrial applications. We finally analyze the error scaling of the estimated polarizabilities obtained by the proposed analytical derivative versus the numerical derivative obtained by the finite difference. Numerical calculations suggest that our analytical derivative requires fewer measurements (runs) on quantum computers than the numerical derivative to achieve the same fixed accuracy.

DOI: 10.1021/acs.jctc.2c01176

Related papers

Quantum Equation of Motion with Orbital Optimization for Computing Molecular Properties in Near-Term Quantum Computing

Phillip W. K. Jensen, Erik Rosendahl Kjellgren, Peter Reinholdt +4 more

Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH 2 in the series of STO-3G/6-31G/6-31G* basis sets and of H 4 and H 2 O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H 4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.

arXiv:2312.12386•Journal of Chemical Theory and Computation•May 3, 2024•18 citations

Theory of analytical energy derivatives for the variational quantum eigensolver

Kosuke Mitarai, Yuya O. Nakagawa, Wataru Mizukami

The variational quantum eigensolver (VQE), an quantum-classical algorithm which can be employed to compute properties of molecules and materials, is an appealing candidate for the application of near-term quantum computers. The authors developed an algorithm for computing analytic energy derivatives with respect to some system parameters within the VQE framework. As many physical properties are defined through the energy derivatives, this work widens the applicable range of the VQE.

arXiv:1905.04054•Physical Review Research•Feb 5, 2020•76 citations

An Interface for Variational Quantum Eigensolver based Energy (VQE-E) and Force (VQE-F) Calculator to Atomic Simulation Environment (ASE)

M R Nirmal, Sarkar, Shampa, Manoj Karunakaran Nambiar

The development of quantum algorithms to solve quantum chemistry problems has offered a promising new paradigm of performing computer simulations at the scale of atoms and molecules. Although majority of the research so far has focused on designing quantum algorithms to compute ground and excited state energies and forces, it is useful to run different simulation tasks, such as geometry optimization, with these algorithms as subroutines. Towards this end, we have created an interface for the Variational Quantum Eigensolver based molecular Energy (VQE-E) and molecular Force (VQE-F) code to the Atomic Simulation Environment (ASE). We demonstrate the working of this hybrid quantum-classical interface by optimizing the geometry of water molecule using a native optimizer implemented in ASE. Furthermore, this interface enables one to compare, combine and use quantum algorithms in conjunction with related classical methods quite easily with minimal coding effort.

arXiv:2209.13876•Sep 28, 2022•2 citations

The Cost of Improving the Precision of the Variational Quantum Eigensolver for Quantum Chemistry

Ivana Miháliková, Matej Pivoluska, Martin Plesch +3 more

New approaches into computational quantum chemistry can be developed through the use of quantum computing. While universal, fault-tolerant quantum computers are still not available, and we want to utilize today's noisy quantum processors. One of their flagship applications is the variational quantum eigensolver (VQE)-an algorithm for calculating the minimum energy of a physical Hamiltonian. In this study, we investigate how various types of errors affect the VQE and how to efficiently use the available resources to produce precise computational results. We utilize a simulator of a noisy quantum device, an exact statevector simulator, and physical quantum hardware to study the VQE algorithm for molecular hydrogen. We find that the optimal method of running the hybrid classical-quantum optimization is to: (i) allow some noise in intermediate energy evaluations, using fewer shots per step and fewer optimization iterations, but ensure a high final readout precision; (ii) emphasize efficient problem encoding and ansatz parametrization; and (iii) run all experiments within a short time-frame, avoiding parameter drift with time. Nevertheless, current publicly available quantum resources are still very noisy and scarce/expensive, and even when using them efficiently, it is quite difficult to perform trustworthy calculations of molecular energies.

arXiv:2111.04965•Nanomaterials•Jan 14, 2022•23 citations

Calculating energy derivatives for quantum chemistry on a quantum computer

Thomas E. O’Brien, Bruno Senjean, Ramiro Sagastizabal +5 more

Abstract Modeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within $$0.014$$ 0.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H $${}_{2}$$ 2 molecule, finding agreement at the equilibrium bond length to within $$0.06$$ 0.06 a.u. ( $$2 \%$$ 2% relative error).

arXiv:1905.03742•npj Quantum Information•Dec 11, 2019•101 citations

Quantum computing methods for electronic states of the water molecule

Teng Bian, Daniel Murphy, Rongxin Xia +2 more

We compare recently proposed methods to compute the electronic state energies of the water molecule on a quantum computer. The methods include the phase estimation algorithm based on Trotter decomposition, the phase estimation algorithm based on the direct implementation of the Hamiltonian, direct measurement based on the implementation of the Hamiltonian and a specific variational quantum eigensolver, Pairwise VQE. After deriving the Hamiltonian using STO-3G basis, we first explain how each method works and then compare the simulation results in terms of gate complexity and the number of measurements for the ground state of the water molecule with different O–H bond lengths. Moreover, we present the analytical analyses of the error and the gate-complexity for each method. While the required number of qubits for each method is almost the same, the number of gates and the error vary a lot. In conclusion, among methods based on the phase estimation algorithm, the second-order direct method provides the most efficient circuit implementations in terms of the gate complexity. Moreover, Pairwise VQE serves the most practical method for near-term applications on the current available quantum computers. Finally the possibility of extending the calculation to excited states and resonances is discussed.

arXiv:1804.05453•Molecular Physics•Feb 20, 2019•38 citations

We develop a quantum-classical hybrid algorithm to calculate the analytical second-order derivative of the energy for the orbital-optimized variational quantum eigensolver (OO-VQE), which is a method to calculate eigenenergies of a given molecular Hamiltonian by utilizing near-term quantum computers and classical computers. We show that all quantities required in the algorithm to calculate the derivative can be evaluated on quantum computers as standard quantum expectation values without using any ancillary qubits. We validate our formula by numerical simulations of quantum circuits for computing the polarizability of the water molecule, which is the second-order derivative of the energy, with respect to the electric field. Moreover, the polarizabilities and refractive indices of thiophene and furan molecules are calculated as a test bed for possible industrial applications. We finally analyze the error scaling of the estimated polarizabilities obtained by the proposed analytical derivative versus the numerical derivative obtained by the finite difference. Numerical calculations suggest that our analytical derivative requires fewer measurements (runs) on quantum computers than the numerical derivative to achieve the same fixed accuracy.

Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH 2 in the series of STO-3G/6-31G/6-31G* basis sets and of H 4 and H 2 O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H 4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.

The development of quantum algorithms to solve quantum chemistry problems has offered a promising new paradigm of performing computer simulations at the scale of atoms and molecules. Although majority of the research so far has focused on designing quantum algorithms to compute ground and excited state energies and forces, it is useful to run different simulation tasks, such as geometry optimization, with these algorithms as subroutines. Towards this end, we have created an interface for the Variational Quantum Eigensolver based molecular Energy (VQE-E) and molecular Force (VQE-F) code to the Atomic Simulation Environment (ASE). We demonstrate the working of this hybrid quantum-classical interface by optimizing the geometry of water molecule using a native optimizer implemented in ASE. Furthermore, this interface enables one to compare, combine and use quantum algorithms in conjunction with related classical methods quite easily with minimal coding effort.

New approaches into computational quantum chemistry can be developed through the use of quantum computing. While universal, fault-tolerant quantum computers are still not available, and we want to utilize today's noisy quantum processors. One of their flagship applications is the variational quantum eigensolver (VQE)-an algorithm for calculating the minimum energy of a physical Hamiltonian. In this study, we investigate how various types of errors affect the VQE and how to efficiently use the available resources to produce precise computational results. We utilize a simulator of a noisy quantum device, an exact statevector simulator, and physical quantum hardware to study the VQE algorithm for molecular hydrogen. We find that the optimal method of running the hybrid classical-quantum optimization is to: (i) allow some noise in intermediate energy evaluations, using fewer shots per step and fewer optimization iterations, but ensure a high final readout precision; (ii) emphasize efficient problem encoding and ansatz parametrization; and (iii) run all experiments within a short time-frame, avoiding parameter drift with time. Nevertheless, current publicly available quantum resources are still very noisy and scarce/expensive, and even when using them efficiently, it is quite difficult to perform trustworthy calculations of molecular energies.

Abstract Modeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within

0.014

0.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H

{}_{2}

2 molecule, finding agreement at the equilibrium bond length to within

0.06

0.06 a.u. (

2 \%

2% relative error).

We compare recently proposed methods to compute the electronic state energies of the water molecule on a quantum computer. The methods include the phase estimation algorithm based on Trotter decomposition, the phase estimation algorithm based on the direct implementation of the Hamiltonian, direct measurement based on the implementation of the Hamiltonian and a specific variational quantum eigensolver, Pairwise VQE. After deriving the Hamiltonian using STO-3G basis, we first explain how each method works and then compare the simulation results in terms of gate complexity and the number of measurements for the ground state of the water molecule with different O–H bond lengths. Moreover, we present the analytical analyses of the error and the gate-complexity for each method. While the required number of qubits for each method is almost the same, the number of gates and the error vary a lot. In conclusion, among methods based on the phase estimation algorithm, the second-order direct method provides the most efficient circuit implementations in terms of the gate complexity. Moreover, Pairwise VQE serves the most practical method for near-term applications on the current available quantum computers. Finally the possibility of extending the calculation to excited states and resonances is discussed.