Localized Quantum Chemistry on Quantum Computers

Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today’s quantum devices. Here we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground state energy. Our algorithm, termed “local active space unitary coupled cluster” (LAS-UCC), scales linearly with system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H2)2 into two H2 molecules and by breaking the two double bonds in trans-butadiene and resources estimates are provided for linear chains of up to 20 H2 molecules.

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

We present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) than complete active space second order perturbation theory (CASPT2) at a significantly lower computational cost. In this work we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT we could perform dynamics with an active space that was computationally too expensive for CASPT2.

Stochastic Generalized Active Space Self-Consistent Field: Theory and Application

An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overheads of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS towards fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the triplet over the quintet spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.

Minimal-active-space multistate density functional theory for excitation energy involving local and charge transfer states

Multistate density functional theory (MSDFT) employing a minimum active space (MAS) is presented to determine charge transfer (CT) and local excited states of bimolecular complexes. MSDFT is a hybrid wave function theory (WFT) and density functional theory, in which dynamic correlation is first incorporated in individual determinant configurations using a Kohn–Sham exchange-correlation functional. Then, nonorthogonal configuration-state interaction is performed to treat static correlation. Because molecular orbitals are optimized separately for each determinant by including Kohn–Sham dynamic correlation, a minimal number of configurations in the active space, essential to representing low-lying excited and CT states of interest, is sufficient to yield the adiabatic states. We found that the present MAS-MSDFT method provides a good description of covalent and CT excited states in comparison with experiments and high-level computational results. Because of the simplicity and interpretive capability through diabatic configuration weights, the method may be useful in dynamic simulations of CT and nonadiabatic processes.