Accurate singlet-triplet gaps in biradicals via the spin averaged anti-Hermitian contracted Schrödinger equation

Exact Ansatz of Fermion-Boson Systems for a Quantum Device

We present an exact Ansatz for the eigenstate problem of mixed fermion-boson systems that can be implemented on quantum devices. Based on a generalization of the electronic contracted Schrödinger equation (CSE), our approach guides a trial wave function to the ground state of any arbitrary mixed Hamiltonian by directly measuring residuals of the mixed CSE on a quantum device. Unlike density functional and coupled cluster theories applied to electron-phonon or electron-photon systems, the accuracy of our approach is not limited by the unknown exchange-correlation functional or the uncontrolled form of the exponential Ansatz. To test the performance of the method, we study the Tavis-Cummings model, commonly used in polaritonic quantum chemistry. Our results demonstrate that the CSE is a powerful tool in the development of quantum algorithms for solving general fermion-boson many-body problems.

Forte: A suite of advanced multireference quantum chemistry methods

Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements.

Many-Body Excited States with a Contracted Quantum Eigensolver

Calculating ground and excited states is an exciting prospect for near-term quantum computing applications, and accurate and efficient algorithms are needed to assess viable directions. We develop an excited-state approach based on the contracted quantum eigensolver (ES-CQE), which iteratively attempts to find a solution to a contraction of the Schrödinger equation projected onto a subspace and does not require a priori information on the system. We focus on the anti-Hermitian portion of the equation, leading to a two-body unitary ansatz. We investigate the role of symmetries, initial states, constraints, and overall performance within the context of the model strongly correlated rectangular H4 system. We show that the ES-CQE achieves near-exact accuracy across the majority of states, covering regions of strong and weak electron correlation, while also elucidating challenging instances for two-body unitary ansatz.

Quantum simulation of conical intersections

We explore the simulation of conical intersections (CIs) on quantum devices, setting the groundwork for potential applications in nonadiabatic quantum dynamics within molecular systems. The intersecting potential energy surfaces of H3+ are computed from a variance-based contracted quantum eigensolver. We show how the CIs can be correctly described on quantum devices using wavefunctions generated by the anti-Hermitian contracted Schrödinger equation ansatz, which is a unitary transformation of wavefunctions that preserves the topography of CIs. A hybrid quantum-classical procedure is used to locate the seam of CIs. Additionally, we discuss the quantum implementation of the adiabatic to diabatic transformation and its relation to the geometric phase effect. Results on noisy intermediate-scale quantum devices showcase the potential of quantum computers in dealing with problems in nonadiabatic chemistry.

Universal Generalization of Density Functional Theory for Static Correlation

A major challenge for density functional theory (DFT) is its failure to treat static correlation, yielding errors in predicted charges, band gaps, van der Waals forces, and reaction barriers. Here we combine one- and two-electron reduced density matrix (1- and 2-RDM) theories with DFT to obtain a universal O(N^{3}) generalization of DFT for static correlation. Using the lowest unitary invariant of the cumulant 2-RDM, we generate a 1-RDM functional theory that corrects the convexity of any DFT functional to capture static correlation in its fractional orbital occupations. Importantly, the unitary invariant yields a predictive theory by revealing the dependence of the correction's strength upon the trace of the two-electron repulsion matrix. We apply the theory to the barrier to rotation in ethylene, the relative energies of the benzynes, as well as an 11-molecule, dissociation benchmark. By inheriting the computational efficiency of DFT without sacrificing the treatment of static correlation, the theory opens new possibilities for the prediction and interpretation of significant quantum molecular effects and phenomena.

Molecules in Environments: Toward Systematic Quantum Embedding of Electrons and Drude Oscillators

We develop a quantum embedding method that enables accurate and efficient treatment of interactions between molecules and an environment, while explicitly including many-body correlations. The molecule is composed of classical nuclei and quantum electrons, whereas the environment is modeled via charged quantum harmonic oscillators. We construct a general Hamiltonian and introduce a variational Ansatz for the correlated ground state of the fully interacting molecule-environment system. This wave function is optimized via the variational Monte Carlo method and the ground state energy is subsequently estimated through the diffusion Monte Carlo method. The proposed scheme allows an explicit many-body treatment of electrostatic, polarization, and dispersion interactions between the molecule and the environment. We study solvation energies and excitation energies of benzene derivatives, obtaining excellent agreement with explicit ab initio calculations and experiments.

Threshold Photoelectron Spectrum of m-Benzyne

Due to their unusual electronic structure, the biradical m-benzyne, C₆H₄, and its cation are of considerable interest in chemistry. Here, the photoion mass-selected threshold photoelectron spectrum of the m-benzyne biradical is presented. An adiabatic ionization energy of 8.65 ± 0.015 eV is derived, while a vibrational progression of 0.10 eV is assigned to the ν₉⁺ ring breathing mode, in excellent agreement with computations. The experimental spectrum was reproduced well by Franck-Condon spectral modeling of the ²A₁ ← X ¹A₁ transition, in which the cation retains a monocyclic C₆ framework. The energetically close-lying bicyclic ²A₂ cation state exhibits low Franck-Condon factors, due to the large change in geometry, and thus cannot be observed.

Reducing the Quantum Many-Electron Problem to Two Electrons with Machine Learning

An outstanding challenge in chemical computation is the many-electron problem where computational methodologies scale prohibitively with system size. The energy of any molecule can be expressed as a weighted sum of the energies of two-electron wave functions that are computable from only a two-electron calculation. Despite the physical elegance of this extended "aufbau" principle, the determination of the distribution of weights─geminal occupations─for general molecular systems has remained elusive. Here we introduce a new paradigm for electronic structure where approximate geminal-occupation distributions are "learned" via a convolutional neural network. We show that the neural network learns the N-representability conditions, constraints on the distribution for it to represent an N-electron system. By training on hydrocarbon isomers with only 2-7 carbon atoms, we are able to predict the energies for isomers of octane as well as hydrocarbons with 8-15 carbons. The present work demonstrates that machine learning can be used to reduce the many-electron problem to an effective two-electron problem, opening new opportunities for accurately predicting electronic structure.

Variational determination of the two-electron reduced density matrix within the doubly occupied configuration interaction framework: Treatments of triplet N-electron systems

This work performs variational determinations of two-electron reduced density matrices corresponding to open-shell N-electron systems within the framework of the doubly occupied configuration interaction treatment, traditionally limited to studies of closed-shell systems. The procedure has allowed us to describe satisfactorily molecular systems in triplet states following two methods. One of them adds hydrogen atoms at an infinite distance of the triplet system studied, constituting a singlet supersystem. Energies and reduced density matrices of the triplet system are obtained by removing the contributions of the added atoms from the singlet supersystem results. The second procedure determines variationally the two-electron reduced density matrices corresponding to the triplet systems by means of adequate couplings of basis-set functions. Both models have been managed by imposing N-representability conditions on the reduced density matrix calculations. Results obtained from these methods for molecular systems in triplet ground states are reported and compared with those provided by benchmark methods.

Interplay of Electronic and Geometric Structure Tunes Organic Biradical Character in Bimetallic Tetrathiafulvalene Tetrathiolate Complexes

The synthesis and design of organic biradicals with tunable singlet-triplet gaps have become the subject of significant research interest, owing to their possible photochemical applications and use in the development of molecular switches and conductors. Recently, tetrathiafulvalene tetrathiolate (TTFtt) has been demonstrated to exhibit such organic biradical character in doubly ionized bimetallic complexes. In this article we use high-level ab initio calculations to interrogate the electronic structure of a series of TTFtt-bridged metal complexes, resolving the factors governing their biradical character and singlet-triplet gaps. We show that the degree of biradical character correlates with a readily measured experimental predictor, the central TTFtt C-C bond length, and that it may be described by a one-parameter model, providing valuable insight for the future rational design of TTFtt based biradical compounds and materials.

Elucidating the molecular orbital dependence of the total electronic energy in multireference problems

The accurate resolution of the chemical properties of strongly correlated systems, such as biradicals, requires the use of electronic structure theories that account for both multi-reference and dynamic correlation effects. A variety of methods exist that aim to resolve the dynamic correlation in multi-reference problems, commonly relying on an exponentially scaling complete-active-space self-consistent-field (CASSCF) calculation to generate reference molecular orbitals (MOs). However, while CASSCF orbitals provide the optimal solution for a selected set of correlated (active) orbitals, their suitability in the quest for the resolution of the total correlation energy has not been thoroughly investigated. Recent research has shown the ability of Kohn-Shan density functional theory to provide improved orbitals for coupled cluster (CC) and Møller-Plesset perturbation theory (MP) calculations. Here, we extend the search for optimal and more cost effective MOs to post-configuration-interaction [post-(CI)] methods, surveying the ability of the MOs obtained with various density functional theory (DFT) functionals, as well as Hartree-Fock and CC and MP calculations to accurately capture the total electronic correlation energy. Applying the anti-Hermitian contracted Schrödinger equation to the dissociation of N₂, the calculation of biradical singlet-triplet gaps, and the transition states of bicylobutane isomerization, we demonstrate that DFT provides a cost-effective alternative to CASSCF in providing reference orbitals for post-CI dynamic correlation calculations.

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg [J. Chem. Phys. 107, 432 (1997)]. This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple-ζ basis set predicts the spin splittings of [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ to be -35.7 and -17.1 kcal mol^-1, respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.