Toward a Quantum Computing Formulation of the Electron Nuclear Dynamics Method via Fukutome Unitary Representation

TL;DR

This paper introduces a quantum computing formulation of the electron nuclear dynamics method, enabling simulation of chemical reactions using hybrid quantum-classical schemes with promising accuracy improvements.

Contribution

It develops the first quantum computing version of END using Fukutome unitary representation, facilitating electronic dynamics simulation on quantum hardware.

Findings

Errors decrease linearly with the number of shots

Fukutome matrices factorize into one-qubit rotations

Successful simulation of H2+ electronic dynamics

Abstract

We present the first installment of the quantum computing (QC) formulation of the electron nuclear dynamics (END) method within the variational quantum simulator (VQS) scheme: END/QC/VQS. END is a time-dependent, variational, on-the-flight, and non-adiabatic method to simulate chemical reactions. END represents nuclei with frozen Gaussian wave packets and electrons with a single-determinantal state in the Thouless non-unitary representation. Within the hybrid quantum/classical VQS, END/QC/VQS evaluates the metric matrix M and gradient vector V of the symplectic END/QC equations on a quantum computer, and calculates basis function integrals and time evolution on a classical computer. To adapt END to QC, we substitute the Thouless non-unitary representation with Fukutome unitary representation. We formulate the first END/QC/VQS version for pure electronic dynamics in chemical models consisting of two-electron units. Therein, Fukutome unitary matrices factorize into products of triads of one-qubit rotational matrices, which leads to a QC encoding of one electron per qubit. We design QC circuits to evaluate M and V in one-electron diatomic molecules. In log2-log2 plots, errors and deviations of those evaluations decrease linearly with the number of shots and with slopes = -1/2. We illustrate an END/QC/VQS simulation with the pure electronic dynamics of H2+. We discuss the present results and future END/QC/QVS extensions.