Physically Motivated Improvements of Variational Quantum Eigensolvers.

The adaptive derivative-assembled pseudo-Trotter variational quantum eigensolver (ADAPT-VQE) has emerged as a pivotal promising approach for electronic structure challenges in quantum chemistry with noisy quantum devices. Nevertheless, to surmount existing technological constraints, this study endeavors to enhance ADAPT-VQE's efficacy. Leveraging insights from the electronic structure theory, we concentrate on optimizing state preparation without added computational burden and guiding ansatz expansion to yield more concise wave functions with expedited convergence toward exact solutions. These advancements culminate in shallower circuits and, as demonstrated, reduced measurement requirements. This research delineates these enhancements and assesses their performance across mono, di, and tridimensional arrangements of H4 models, as well as in the water molecule. Ultimately, this work attests to the viability of physically motivated strategies in fortifying ADAPT-VQE's efficiency, marking a significant stride in quantum chemistry simulations.

Quantum Simulation of Molecular Response Properties in the NISQ Era.

Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this article, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offer a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate the state-specific response properties. Further, through noiseless quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster-based linear response models due to the improved quality of the ground-state wave function obtained using the ADAPT-VQE algorithm.

Generalization of the Tensor Product Selected CI Method for Molecular Excited States.

In a recent paper [JCTC, 2020, 16, 6098], we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact, local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excited-state character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons with respect to the extrapolated result (i.e., near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 doubly excited (biexcitonic) states.

Quantum self-consistent equation-of-motion method for computing molecular excitation energies, ionization potentials, and electron affinities on a quantum computer.

Near-term quantum computers are expected to facilitate material and chemical research through accurate molecular simulations. Several developments have already shown that accurate ground-state energies for small molecules can be evaluated on present-day quantum devices. Although electronically excited states play a vital role in chemical processes and applications, the search for a reliable and practical approach for routine excited-state calculations on near-term quantum devices is ongoing. Inspired by excited-state methods developed for the unitary coupled-cluster theory in quantum chemistry, we present an equation-of-motion-based method to compute excitation energies following the variational quantum eigensolver algorithm for ground-state calculations on a quantum computer. We perform numerical simulations on H2, H4, H2O, and LiH molecules to test our quantum self-consistent equation-of-motion (q-sc-EOM) method and compare it to other current state-of-the-art methods. q-sc-EOM makes use of self-consistent operators to satisfy the vacuum annihilation condition, a critical property for accurate calculations. It provides real and size-intensive energy differences corresponding to vertical excitation energies, ionization potentials and electron affinities. We also find that q-sc-EOM is more suitable for implementation on NISQ devices as it is expected to be more resilient to noise compared with the currently available methods.

New Local Explorations of the Unitary Coupled Cluster Energy Landscape.

The recent quantum information boom has effected a resurgence of interest in the unitary coupled cluster (UCC) theory. Our group's interest in local energy landscapes of unitary ansätze prompted us to investigate the classical approach of truncating the Taylor series expansion (instead of a perturbative expansion) of the UCC with singles and doubles (UCCSD) energy at the second order. This amounts to an approach where the electron correlation energy is estimated by taking a single Newton-Raphson step from Hartree-Fock toward UCCSD. Such an approach has been explored previously, but the accuracy was not extensively studied. In this paper, we investigate the performance and observe similar pathologies to the linearized coupled cluster with singles and doubles. We introduce the use of derivatives of order three or greater to help partially recover the variational lower bound of true UCCSD, restricting these derivatives to those of the "unmixed" category in order to simplify the model. By testing the approach on several potential energy surfaces and reaction energies, we find this "diagonal" approximation to higher order terms to be effective at reducing sensitivity near singularities for strongly correlated regimes, while not significantly diminishing the accuracy of weakly correlated systems.

Symmetry Breaking Slows Convergence of the ADAPT Variational Quantum Eigensolver.

Because quantum simulation of molecular systems is expected to provide the strongest advantage over classical computing methods for systems exhibiting strong electron correlation, it is critical that the performance of VQEs be assessed for strongly correlated systems. For classical simulation, strong correlation often results in symmetry breaking of the Hartree-Fock reference, leading to Löwdin's well-known "symmetry dilemma", whereby accuracy in the energy can be increased by breaking spin or spatial symmetries. Here, we explore the impact of symmetry breaking on the performance of ADAPT-VQE using two strongly correlated systems: (i) the "fermionized" anisotropic Heisenberg model, where the anisotropy parameter controls the correlation in the system, and (ii) symmetrically stretched linear H4, where correlation increases with increasing H-H separation. In both of these cases, increasing the level of correlation of the system leads to spontaneous symmetry breaking (parity and S^2, respectively) of the mean-field solutions. We analyze the role that symmetry breaking in the reference states and orbital mappings of the fermionic Hamiltonians have in the compactness and performance of ADAPT-VQE. We observe that improving the energy of the reference states by breaking symmetry has a deleterious effect on ADAPT-VQE by increasing the length of the ansatz necessary for energy convergence and exacerbating the problem of "gradient troughs".

ONIOM Method with Charge Transfer Corrections (ONIOM-CT): Analytic Gradients and Benchmarking.

Hybrid methods such as ONIOM that treat different regions of a large molecule using different methods are widely used to investigate chemical reactions in a variety of materials and biological systems. However, there are inherent sources of significant errors due to the standard treatment of the boundary between the regions using hydrogen link atoms. In particular, an unbalanced charge distribution in the chemically important model region is a potential source of such problems. We have previously suggested ONIOM-CT (ONIOM with charge transfer corrections) which addresses this issue by applying a potential in the form of point charges to obtain a desired charge redistribution. The metric for charge redistribution relies on the type of population analysis used to obtain the charges. ONIOM-CT has been implemented using Mulliken and Löwdin population analyses and has been shown to improve computed reaction energies for illustrative chemical reactions. In this work, we derive and implement the analytic gradients for ONIOM-CT that requires solving two sets of coupled-perturbed self-consistent equations, one each for the model system and the full system. However, both are needed only at the low level of theory, allowing for an efficient formulation and implementation for both Mulliken and Löwdin population analyses. Benchmarking and illustrative geometry optimizations have been carried out for a previously studied set of reactions involving a single link atom between regions. Additionally, we have generalized our method for the treatment of model systems involving multiple link atoms to enable applications for a broader set of problems. The generalized methods are illustrated for both charge models. Furthermore, we have studied a set of three proton transfer reactions and demonstrate that significant improvement is achieved by ONIOM-CT over ONIOM using both Mulliken and Löwdin population analyses.

Coupled Electron Pair-Type Approximations for Tensor Product State Wave Functions.

Size extensivity, defined as the correct scaling of energy with system size, is a desirable property for any many-body method. Traditional configuration interaction (CI) methods are not size extensive, hence the error increases as the system gets larger. Coupled electron pair approximation (CEPA) methods can be constructed as simple extensions of a truncated CI that ensures size extensivity. One of the major issues with the CEPA and its variants is that singularities arise in the amplitude equations when the system starts to be strongly correlated. In this work, we extend the traditional Slater determinant based coupled electron pair approaches like CEPA-0, averaged coupled-pair functional, and average quadratic coupled-cluster to a new formulation based on tensor product states (TPS). We show that a TPS basis can often be chosen such that it removes the singularities that commonly destroy the accuracy of CEPA based methods. A suitable TPS representation can be formed by partitioning the system into separate disjoint clusters and forming the final wave function as the tensor product of the many body states of these clusters. We demonstrate the application of these methods on simple bond breaking systems such as CH4 and F2 where determinant based CEPA methods fail. We further apply the TPS-CEPA approach to stillbene isomerization and few planar π-conjugated systems. Overall, the results show that the TPS-CEPA method can remove the singularities and provide improved numerical results compared to common electronic structure methods.

Revealing the Contest between Triplet-Triplet Exchange and Triplet-Triplet Energy Transfer Coupling in Correlated Triplet Pair States in Singlet Fission.

Understanding the separation of the correlated triplet pair state 1(TT) intermediate is critical for leveraging singlet fission to improve solar cell efficiency. This separation mechanism is dominated by two key interactions: (i) the exchange interaction (K) between the triplets which leads to the spin splitting of the biexciton state into 1(TT),3(TT) and 5(TT) states, and (ii) the triplet-triplet energy transfer integral (t) which enables the formation of the spatially separated (but still spin entangled) state 1(T···T). We develop a simple ab initio technique to compute both the biexciton exchange (K) and biexciton transfer coupling. Our key findings reveal new conditions for successful correlated triplet pair state dissociation. The biexciton exchange interaction needs to be ferromagnetic or negligible to the triplet energy transfer for favorable dissociation. We also explore the effect of chromophore packing to reveal geometries where these conditions are achieved for tetracene.

Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference.

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 1066 determinants for the exact FCI solution.

Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules.

We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.

Selected Configuration Interaction in a Basis of Cluster State Tensor Products.

Selected configuration interaction (SCI) methods are currently enjoying a resurgence due to several recent developments which improve either the overall computational efficiency or the compactness of the resulting SCI vector. These recent advances have made it possible to get full CI (FCI) quality results for much larger orbital active spaces compared to conventional approaches. However, due to the starting assumption that the FCI vector has only a small number of significant Slater determinants, SCI becomes intractable for systems with strong correlation. This paper introduces a method for developing SCI algorithms in a way which exploits local molecular structure to significantly reduce the number of SCI variables. The proposed method is defined by first grouping the orbitals into clusters over which we can define many-particle cluster states. We then directly perform the SCI algorithm in a basis of tensor products of cluster states instead of Slater determinants. While the approach is general for arbitrarily defined cluster states, we find significantly improved performance by defining cluster states through a Tucker decomposition of the global (and sparse) SCI vector. To demonstrate the potential of this method, called tensor product selected configuration interaction (TPSCI), we present numerical results for a diverse set of examples: (1) modified Hubbard model with different inter- and intracluster hopping terms, (2) less obviously clusterable cases of bond breaking in N2 and F2, and (3) ground state energies of large planar π-conjugated systems with active spaces of up to 42 electrons in 42 orbitals. These numerical results show that TPSCI can be used to significantly reduce the number of SCI variables in the variational space, thus paving a path for extending these deterministic and variational SCI approaches to a wider range of physical systems.

Spin-Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach.

We present a practical approach for computing the Breit-Pauli spin-orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin-orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner-Eckart theorem. A mean field spin-orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin-orbit coupling constants, despite the fact that the efficient wave function approximations we explore neglect the bulk of dynamical correlation.

Is the Trotterized UCCSD Ansatz Chemically Well-Defined?

The variational quantum eigensolver (VQE) has emerged as one of the most promising near-term quantum algorithms that can be used to simulate many-body systems such as molecular electronic structures. Serving as an attractive ansatz in the VQE algorithm, unitary coupled cluster (UCC) theory has seen a renewed interest in recent literature. However, unlike the original classical UCC theory, implementation on a quantum computer requires a finite-order Suzuki-Trotter decomposition to separate the exponentials of the large sum of Pauli operators. While previous literature has recognized the nonuniqueness of different orderings of the operators in the Trotterized form of UCC methods, the question of whether or not different orderings matter at the chemical scale has not been addressed. In this Letter, we explore the effect of operator ordering on the Trotterized UCCSD ansatz, as well as the much more compact k-UpCCGSD ansatz recently proposed by Lee et al. [ J. Chem. Theory Comput. , 2019 , 15 , 311 . arXiv, 2019 , quant-ph:1909.09114. https://arxiv.org/abs/1909.09114 ]. We observe a significant, system-dependent variation in the energies of Trotterizations with different operator orderings. The energy variations occur on a chemical scale, sometimes on the order of hundreds of kcal/mol. This Letter establishes the need to define not only the operators present in the ansatz but also the order in which they appear. This is necessary for adhering to the quantum chemical notion of a "model chemistry", in addition to the general importance of scientific reproducibility. As a final note, we suggest a useful strategy to select out of the combinatorial number of possibilities, a single well-defined and effective ordering of the operators.

Simple and Efficient Truncation of Virtual Spaces in Embedded Wave Functions via Concentric Localization.

We present a strategy to generate "concentrically local orbitals" for the purpose of decreasing the computational cost of wave function-in-density functional theory (WF-in-DFT) embedding. The concentric localization is performed for the virtual orbitals by first projecting the virtual space onto atomic orbitals centered on the embedded atoms. Using a one-particle operator, these projected orbitals are then taken as a starting point to iteratively span the virtual space, recursively creating virtual orbital "shells" with consecutively decreasing correlation energy recovery at each iteration. This process can be repeated to convergence, allowing for tunable accuracy. Assessment of the proposed scheme is performed by application to the potential energy diagram of the Menshutkin reaction of chloromethane and ammonia inside a segment of a carbon nanotube and the torsional potential of a simplified version of the retinal chromophore.

An adaptive variational algorithm for exact molecular simulations on a quantum computer.

Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.

A Combined Spin-Flip and IP/EA Approach for Handling Spin and Spatial Degeneracies: Application to Double Exchange Systems.

Many multiconfigurational systems, such as single-molecule magnets, are difficult to study using traditional computational methods due to the simultaneous existence of both spin and spatial degeneracies. In this work, a new approach termed n-spin-flip ionization potential/electron affinity ( nSF-IP or nSF-EA) is introduced which combines the spin-flip method of Anna Krylov with particle-number changing IP/EA methods. We demonstrate the efficacy of the approach by applying it to the strongly correlated N2+, as well as several double exchange systems. We also demonstrate that when these systems are well-described by a double exchange model Hamiltonian, only 1SF-IP/EA is required to extract the double exchange parameters and accurately predict energies for the low-spin states. This significantly reduces the computational effort for studying such systems. The effects of including additional excitations (using a RAS- nSF-IP/EA scheme) are also examined, with particular emphasis on hole and particle excitations.

Automatic Partition of Orbital Spaces Based on Singular Value Decomposition in the Context of Embedding Theories.

We present a simple approach for orbital space partitioning to be employed in the projection-based embedding theory developed by Goodpaster and co-workers [ Manby et al. J. Chem. Theory Comput. 2012 , 8 , 2564 ]. Once the atoms are assigned to the desired subspaces, the molecular orbitals are projected onto the atomic orbitals centered on active atoms and then singular value decomposed. The right singular vectors are used to rotate the initial molecular orbitals, taking the largest gap in the singular values spectrum to define the most suitable partition of the occupied orbital space. This scheme is free from numerical parameters, contrary to the Mulliken charge threshold or the completeness criterion previously used. The performance of this new prescription is assessed in a test set of several distinct reactions, the deprotonation of decanoic acid, the Diels-Alder reaction of 1,3-butadiene and octadecanonaene, the torsional potential of a retinal derivative, and the critical points along the reaction coordinate of an example of the Menshutkin SN2 reaction inside a carbon nanotube.

Simple Rule To Predict Boundedness of Multiexciton States in Covalently Linked Singlet-Fission Dimers.

Because of the potential for increasing solar cell efficiencies, significant effort has been spent understanding the mechanism of singlet fission. We provide a simple connectivity rule to predict whether the through-bond coupling will be stabilizing or destabilizing for the 1(TT) state in covalently linked singlet-fission chromophores. By drawing an analogy between the chemical system and a simple spin-lattice, one is able to determine the ordering of the multiexciton spin state via a generalized usage of Ovchinnikov's rule. This allows one to predict (without any computation) whether the 1(TT) multiexciton state will be bound or unbound with respect to the separated triplets in covalently linked singlet-fission dimers. To test our hypothesis, we have performed ab initio calculations on a systematic series of covalently linked singlet-fission dimers. Numerical examples are given, and the limitations of the proposed theory are explored.