Abstract:The calculation of non-covalent interaction energies on noisy intermediate-scale quantum (NISQ) computers appears to be challenging with straightforward application of existing quantum algorithms. For example, the use of the standard supermolecular method with the variational quantum eigensolver (VQE) would require extremely precise resolution of the total energies of the fragments to provide for accurate subtraction to the interaction energy. Here we present a symmetry-adapted perturbation theory (SAPT) method that may provide interaction energies with high quantum resource efficiency. Of particular note, we present a quantum extended random-phase approximation (ERPA) treatment of the SAPT second-order induction and dispersion terms, including exchange counterparts. Together with previous work on first-order terms (Chem. Sci., 2022, 13, 3094), this provides a recipe for complete SAPT(VQE) interaction energies up to second order, which is a well established truncation. The SAPT interaction energy terms are computed as first-level observables with no subtraction of monomer energies invoked, and the only quantum observations needed are the VQE one- and two-particle density matrices. We find empirically that SAPT(VQE) can provide accurate interaction energies even with coarsely optimized, low circuit depth wavefunctions from a quantum computer, simulated through ideal statevectors. The errors of the total interaction energy are orders of magnitude lower than the corresponding VQE total energy errors of the monomer wavefunctions. In addition, we present heme-nitrosyl model complexes as a system class for near term quantum computing simulations. They are strongly correlated, biologically relevant and difficult to simulate with classical quantum chemical methods. This is illustrated with density functional theory (DFT) as the predicted interaction energies exhibit a strong sensitivity with respect to the choice of functional. Thus, this work paves the way to obtain accurate interaction energies on a NISQ-era quantum computer with few quantum resources. It is the first step in alleviating one of the major challenges in quantum chemistry, where in-depth knowledge of both the method and system is required a priori to reliably generate accurate interaction energies.

Rigorous Screened Interactions for Realistic Correlated Electron Systems

Rigorous screened interactions for realistic correlated electron systems

Towards a first-principles determination of effective Coulomb interactions in correlated electron materials: Role of intershell interactions

Untangling Strongly and Weakly Interacting Configurations in Many-electron Wave Functions

Effective local potentials for density and density-matrix functional approximations with non-negative screening density

Fast, efficient, and accurate dielectric screening using a local, real-space approach

Fast, efficient, and accurate dielectric screening using a local real-space approach

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Interpolating many-body wave functions for accelerated molecular dynamics on the near-exact electronic surface

Ab Initio Electron-Phonon Interactions in Correlated Electron Systems

Non conventional screening of the Coulomb interaction in low dimensional and finite size system

Dynamical screening in correlated electron systems - from lattice models to realistic materials

Doubly Screened Coulomb Correction Approach for Strongly Correlated Systems

Accurate non-covalent interaction energies on noisy intermediate-scale quantum computers via second-order symmetry-adapted perturbation theory

Variational formulation of dynamical electronic response functions in Bethe-Salpeter, (screened) Hartree-Fock, Hybrid-DFT approaches

Towards Efficient and Unified Treatment of Static and Dynamic Correlations in Generalized Kohn-Sham Density Functional Theory

Toward Efficient and Unified Treatment of Static and Dynamic Correlations in Generalized Kohn-Sham Density Functional Theory

A regularized second-order correlation method from Green's function theory

Electronic structure of low-dimensional inorganic/organic interfaces: Hybrid density functional theory, $G_0W_0$, and electrostatic models

Perturbative Expansion in Reciprocal Space: Bridging Microscopic and Mesoscopic Descriptions of Molecular Interactions