Fast Emulation of Fermionic Circuits with Matrix Product States.

We describe a matrix product state (MPS) extension for the Fermionic Quantum Emulator (FQE) software library. We discuss the theory behind symmetry-adapted MPSs for approximating many-body wave functions of spin-1/2 Fermions, and we present an open-source, MPS-enabled implementation of the FQE interface (MPS-FQE). The software uses the open-source pyblock3 and block2 libraries for most elementary tensor operations, and it can largely be used as a drop-in replacement for FQE that allows for more efficient but approximate emulation of larger Fermionic circuits. Finally, we show several applications relevant to both near-term and fault-tolerant quantum algorithms where approximate emulation of larger systems is expected to be useful: characterization of state preparation strategies for quantum phase estimation, the testing of different variational quantum eigensolver ansätze, the numerical evaluation of Trotter errors, and the simulation of general quantum dynamics problems. In all these examples, approximate emulation with MPS-FQE allows us to treat systems that are significantly larger than those accessible with a full statevector emulator.

Quantum harmonic free energies for biomolecules and nanomaterials.

Obtaining the free energy of large molecules from quantum mechanical energy functions is a long-standing challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of large size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocrystals, and show that the error per atom decreases with system size in the limit of large systems. We further show that we can obtain the vibrational contributions to the binding free energies of prototypical protein-ligand complexes where exact computation is too expensive to be practical. Our work raises the possibility of routine quantum mechanical estimates of thermodynamic quantities in complex systems.

Fast Exchange with Gaussian Basis Set Using Robust Pseudospectral Method.

In this article, we present an algorithm to efficiently evaluate the exchange matrix in periodic systems when a Gaussian basis set with pseudopotentials is used. The usual algorithm for evaluating exchange matrix scales cubically with the system size because one has to perform O(N2) fast Fourier transform (FFT). Here, we introduce an algorithm that retains the cubic scaling but reduces the prefactor significantly by eliminating the need to do FFTs during each exchange build. This is accomplished by representing the products of Gaussian basis function using a linear combination of an auxiliary basis the number of which scales linearly with the size of the system. We store the potential due to these auxiliary functions in memory, which allows us to obtain the exchange matrix without the need to do FFT, albeit at the cost of additional memory requirement. Although the basic idea of using auxiliary functions is not new, our algorithm is cheaper due to a combination of three ingredients: (a) we use a robust pseudospectral method that allows us to use a relatively small number of auxiliary basis to obtain high accuracy; (b) we use occ-RI exchange, which eliminates the need to construct the full exchange matrix; and (c) we use the (interpolative separable density fitting) ISDF algorithm to construct these auxiliary basis sets that are used in the robust pseudospectral method. The resulting algorithm is accurate, and we note that the error in the final energy decreases exponentially rapidly with the number of auxiliary functions.

Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers.

An accurate assessment of how quantum computers can be used for chemical simulation, especially their potential computational advantages, provides important context on how to deploy these future devices. To perform this assessment reliably, quantum resource estimates must be coupled with classical computations attempting to answer relevant chemical questions and to define the classical algorithms simulation frontier. Herein, we explore the quantum computation and classical computation resources required to assess the electronic structure of cytochrome P450 enzymes (CYPs) and thus define a classical-quantum advantage boundary. This is accomplished by analyzing the convergence of density matrix renormalization group plus n-electron valence state perturbation theory (DMRG+NEVPT2) and coupled-cluster singles doubles with noniterative triples [CCSD(T)] calculations for spin gaps in models of the CYP catalytic cycle that indicate multireference character. The quantum resources required to perform phase estimation using qubitized quantum walks are calculated for the same systems. Compilation into the surface code provides runtime estimates to compare directly to DMRG runtimes and to evaluate potential quantum advantage. Both classical and quantum resource estimates suggest that simulation of CYP models at scales large enough to balance dynamic and multiconfigurational electron correlation has the potential to be a quantum advantage problem and emphasizes the important interplay between classical computations and quantum algorithms development for chemical simulation.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Conservation laws in coupled cluster dynamics at finite temperature.

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137-6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest's theorem) for all one-particle properties while remaining energy conserving for time-independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of non-equilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

A coupled cluster framework for electrons and phonons.

We describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard-Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron-phonon coupling. We also report the implementation of electron-phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

Recent developments in the PySCF program package.

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This paper explains the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems.

We discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas (UEG) at warm, dense conditions, and some simple materials. The performance of model systems at high temperatures is encouraging: for the one-dimensional Hubbard model, FT-CCSD provides a qualitatively accurate description of finite-temperature correlation effects even at U = 8, and it allows for the computation of systematically improvable exchange-correlation energies of the warm, dense UEG over a wide range of conditions. We highlight the obstacles that remain in using the method for realistic ab initio calculations on materials.

Generalized single excitation configuration interaction: an investigation into the impact of the inclusion of non-orthogonality on the calculation of core-excited states.

In this paper, we investigate different non-orthogonal generalizations of the configuration interaction with single substitutions (CIS) method for the calculation of core-excited states. Fully non-orthogonal CIS (NOCIS) has been described previously for species with singlet and doublet ground states, and this paper reports the extension to molecules in their triplet ground state. In addition to NOCIS, we present a novel method, one-center NOCIS (1C-NOCIS), for open-shell molecules which is intermediate between NOCIS and the computationally less demanding static exchange approximation (STEX). We explore this hierarchy of spin-pure methods for core excitations of molecules with singlet, doublet, and triplet ground states. We conclude that, while NOCIS provides the best results and preserves the spatial symmetry of the wavefunction, 1C-NOCIS retains much of the accuracy of NOCIS at a dramatically reduced cost. For molecules with closed-shell ground states, STEX and 1C-NOCIS are identical.

Time-Dependent Coupled Cluster Theory on the Keldysh Contour for Nonequilibrium Systems.

We leverage the Keldysh formalism to extend our implementation of finite temperature coupled cluster theory [J. Chem. Theory Comput. 2018, 14, 5690-5700] to thermal systems that have been driven out of equilibrium. The resulting Keldysh coupled cluster theory is discussed in detail. We describe the implementation of the equations necessary to perform Keldysh coupled cluster singles and doubles calculations of finite-temperature dynamics, and we apply the method to some simple systems including a Hubbard model with a Peierls phase and an ab initio model of warm-dense silicon subject to an ultrafast XUV pulse.

Non-Orthogonal Configuration Interaction with Single Substitutions for Core-Excited States: An Extension to Doublet Radicals.

In this paper, we present an open-shell extension of the non-orthogonal configuration interaction singles (NOCIS) method for the calculation of core-excited states, intended for peak assignment in XAS spectra of doublet radicals. This extension requires the consideration of additional configurations due to the singly occupied open-shell orbital, and the addition of essential orbital relaxation effects is found to provide a significant improvement on standard CIS, while maintaining the desirable properties of spin purity, variationality, and size consistency. We apply this method to the calculation of core excitations for several open-shell molecules and demonstrate that it performs competitively with other available methods, despite a lack of dynamic correlation. In particular, relative to CVS-ADC(2)-x, RMS error is reduced by a factor of 6 over usual orthogonal CIS and is comparable to time-dependent density functional theory with the best short-range corrected functionals.

A Time-Dependent Formulation of Coupled-Cluster Theory for Many-Fermion Systems at Finite Temperature.

We present a time-dependent formulation of coupled cluster theory. This theory allows for direct computation of the free energy of quantum systems at finite temperature by imaginary time integration and is closely related to the thermal cluster cumulant theory of Mukherjee and co-workers [ Chem. Phys. Lett. 1992 , 192 , 55 -61 ; Phys. Rev. E 1993 , 48 , 3373 -3389 ; Chem. Phys. Lett. 2001 , 335 , 281 -288 ; Chem. Phys. Lett. 2002 , 352 , 63 -69 ; Int. J. Mod. Phys. B 2003 , 17 , 5367 -5377 ]. Our derivation of the finite-temperature theory highlights connections to perturbation theory and to zero-temperature coupled cluster theory. We show explicitly how the finite-temperature coupled cluster singles and doubles amplitude equations can be derived in analogy with the zero-temperature theory and how response properties can be efficiently computed using a variational Lagrangian. We discuss the implementation for realistic systems and showcase the potential utility of the method with calculations of the exchange correlation energy of the uniform electron gas under warm dense matter conditions.

Non-orthogonal configuration interaction with single substitutions for the calculation of core-excited states.

In this paper, we present the non-orthogonal configuration interaction singles (NOCIS) method for calculating core-excited states of closed-shell molecules. NOCIS is a black-box variant of NOCI, which uses A different core-ionized determinants for a molecule with A atoms of a given element to form single substitutions. NOCIS is a variational, spin-pure, size-consistent ab initio method that dramatically improves on standard CIS by capturing essential orbital relaxation effects, in addition to essential configuration interaction. We apply it to the calculation of core-excitations for several smaller molecules and demonstrate that it performs competitively with other Hartree-Fock and DFT-based methods. We also benchmark it in several basis sets.

Second order Møller-Plesset and coupled cluster singles and doubles methods with complex basis functions for resonances in electron-molecule scattering.

The method of complex basis functions is applied to molecular resonances at correlated levels of theory. Møller-Plesset perturbation theory at second order and equation-of-motion electron attachment coupled-cluster singles and doubles (EOM-EA-CCSD) methods based on a non-Hermitian self-consistent-field reference are used to compute accurate Siegert energies for shape resonances in small molecules including N2-, CO-, CO2-, and CH2O-. Analytic continuation of complex 𝜃-trajectories is used to compute Siegert energies, and the 𝜃-trajectories of energy differences are found to yield more consistent results than those of total energies. The ability of such methods to accurately compute complex potential energy surfaces is investigated, and the possibility of using EOM-EA-CCSD for Feshbach resonances is explored in the context of e-helium scattering.

Stabilizing potentials in bound state analytic continuation methods for electronic resonances in polyatomic molecules.

The computation of Siegert energies by analytic continuation of bound state energies has recently been applied to shape resonances in polyatomic molecules by several authors. We critically evaluate a recently proposed analytic continuation method based on low order (type III) Padé approximants as well as an analytic continuation method based on high order (type II) Padé approximants. We compare three classes of stabilizing potentials: Coulomb potentials, Gaussian potentials, and attenuated Coulomb potentials. These methods are applied to a model potential where the correct answer is known exactly and to the Πg2 shape resonance of N2- which has been studied extensively by other methods. Both the choice of stabilizing potential and method of analytic continuation prove to be important to the accuracy of the results. We conclude that an attenuated Coulomb potential is the most effective of the three for bound state analytic continuation methods. With the proper potential, such methods show promise for algorithmic determination of the positions and widths of molecular shape resonances.