Benchmark variational coupled cluster doubles results

Exploring Ground and Excited States Via Single Reference Coupled-Cluster Theory and Algebraic Geometry

The exploration of the root structure of coupled cluster (CC) equations holds both foundational and practical significance for computational quantum chemistry. This study provides insight into the intricate root structures of these nonlinear equations at both the CCD and CCSD level of theory. We utilize computational techniques from algebraic geometry, specifically the monodromy and parametric homotopy continuation methods, to calculate the full solution set. We compare the computed CC roots against various established theoretical upper bounds, shedding light on the accuracy and efficiency of these bounds. We hereby focus on the dissociation processes of four-electron systems such as (H2)2 in both D2h and D∞h configurations, H4 symmetrically distorted on a circle, and lithium hydride. We moreover investigate the ability of single-reference CC solutions to approximate excited state energies. We find that multiple CC roots describe energies of excited states with high accuracy. Our investigations reveal that for systems like lithium hydride, CC not only provides high-accuracy approximations to several excited state energies but also to the states themselves.

A new "gold standard": Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models

The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.

S-Diagnostic─An a Posteriori Error Assessment for Single-Reference Coupled-Cluster Methods

We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the S-diagnostic. We provide a derivation of the S-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the S-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of systems with varying numerical difficulty, and (3) the square-planar copper complexes [CuCl4]2-, [Cu(NH3)4]2+, and [Cu(H2O)4]2+. Throughout the numerical investigations, the S-diagnostic is compared to other SRCC diagnostic procedures, that is, the T1, D1, max T2, and D2 diagnostics as well as different indices of multideterminantal and multireference character in coupled-cluster theory. Our numerical investigations show that the S-diagnostic outperforms the T1, D1, max T2 and D2 diagnostics and is comparable to the indices of multideterminantal and multireference character in coupled-cluster theory in their individual fields of applicability. The experiments investigating the performance of the S-diagnostic for geometry optimizations using SRCC reveal that the S-diagnostic correlates well with different error measures at a high level of statistical relevance. The experiments investigating the performance of the S-diagnostic for electronic correlation simulations show that the S-diagnostic correctly predicts strong multireference regimes. The S-diagnostic, moreover, correctly detects the successful SRCC computations for [CuCl4]2-, [Cu(NH3)4]2+, and [Cu(H2O)4]2+, which have been known to be misdiagnosed by T1 and D1 diagnostics in the past. This shows that the S-diagnostic is a promising candidate for an a posteriori diagnostic for SRCC calculations.

Two Algorithms for Excited-State Quantum Solvers: Theory and Application to EOM-UCCSD

Near-term quantum devices promise to revolutionize quantum chemistry, but simulations using the current noisy intermediate-scale quantum (NISQ) devices are not practical due to their high susceptibility to errors. This motivated the design of NISQ algorithms leveraging classical and quantum resources. While several developments have shown promising results for ground-state simulations, extending the algorithms to excited states remains challenging. This paper presents two cost-efficient excited-state algorithms inspired by the classical Davidson algorithm. We implemented the Davidson method into the quantum self-consistent equation-of-motion unitary coupled-cluster (q-sc-EOM-UCC) excited-state method adapted for quantum hardware. The circuit strategies for generating desired excited states are discussed, implemented, and tested. We demonstrate the performance and accuracy of the proposed algorithms (q-sc-EOM-UCC/Davidson and its variational variant) by simulations of H2, H4, LiH, and H2O molecules. Similar to the classical Davidson scheme, q-sc-EOM-UCC/Davidson algorithms are capable of targeting a small number of excited states of the desired character.

Benchmarking ionization potentials using the simple pCCD model

The electron-detachment energy is measured by the ionization potential (IP). As a result, it is a fundamental, observable and important molecular electronic signature in photoelectron spectroscopy. A precise theoretical prediction of electron-detachment energies or ionization potentials is essential for organic optoelectronic systems like transistors, solar cells, or light-emitting diodes. In this work, we benchmark the performance of the recently presented IP variant of the equation-of-motion pair coupled cluster doubles (IP-EOM-pCCD) model to determine IPs. Specifically, the predicted ionization energies are compared to experimental results and higher-order coupled cluster theories based on statistically assessing 201 electron-detached states of 41 organic molecules for three different molecular orbital basis sets and two sets of particle-hole operators. While IP-EOM-pCCD features a reasonable spread and skewness of ionization energies, its mean error and standard deviation differ by up to 1.5 eV from reference data. Our study, thus, highlights the importance of dynamical correlation to reliably predict IPs from a pCCD reference function in small organic molecules.

Hybrid Auxiliary Field Quantum Monte Carlo for Molecular Systems

We propose a quantum Monte Carlo approach to solve the many-body Schrödinger equation for the electronic ground state. The method combines optimization from variational Monte Carlo and propagation from auxiliary field quantum Monte Carlo in a way that significantly alleviates the sign problem. In application to molecular systems, we obtain highly accurate results for configurations dominated by either dynamic or static electronic correlation.

Quantum Simulation of Molecules in Solution

Quantum chemical calculations on quantum computers have been focused mostly on simulating molecules in the gas phase. Molecules in liquid solution are, however, most relevant for chemistry. Continuum solvation models represent a good compromise between computational affordability and accuracy in describing solvation effects within a quantum chemical description of solute molecules. In this work, we extend the variational quantum eigensolver to simulate solvated systems using the polarizable continuum model. To account for the state dependent solute-solvent interaction we generalize the variational quantum eigensolver algorithm to treat non-linear molecular Hamiltonians. We show that including solvation effects does not impact the algorithmic efficiency. Numerical results of noiseless simulations for molecular systems with up to 12 spin-orbitals (qubits) are presented. Furthermore, calculations performed on a simulated noisy quantum hardware (IBM Q, Mumbai) yield computed solvation free energies in fair agreement with the classical calculations.

Benchmark Density Functional Theory Approach for the Calculation of Bond Dissociation Energies of the M-O2 Bond: A Key Step in Water Splitting Reactions

A very fascinating aspect in quantum chemical research is to determine the accurate and cost-effective methods for the calculation of electronic and structural properties through a benchmark study. The current study focuses on the performance evaluation of density functional theory methods for the accurate measurement of bond dissociation energies (BDEs) of chemically important M-O₂ bonds in water splitting reactions. The BDE measurement has got noteworthy attention due to its importance in all areas of chemistry. For BDE measurements of M-O₂ bonds in five metal complexes with oxygen molecules, 14 density functionals (DFs) are chosen from seven classes of DFs with two series of mixed basis sets. A combination of pseudopotential and Pople basis sets [LANL2DZ & 6-31G(d) and SDD & 6-31+G(d)] are used as a series of mixed basis sets. The B3LYP-GD3BJ functional with LANL2DZ & 6-31G(d) gives outstanding results due to low deviations, error, and the best Pearson's correlation (R) between the experimental and theoretical data. Our study suggested an efficient, low-cost, precise, and accurate B3LYP-GD3BJ/LANL2DZ & 6-31G(d) level of theory for BDE of the M-O₂ bond, which may be useful for chemists working in the field of energy generation and utilization.

Regularized second-order correlation methods for extended systems

Second-order Møller-Plesset perturbation theory (MP2) constitutes the simplest form of many-body wavefunction theory and often provides a good compromise between efficiency and accuracy. There are, however, well-known limitations to this approach. In particular, MP2 is known to fail or diverge for some prototypical condensed matter systems like the homogeneous electron gas (HEG) and to overestimate dispersion-driven interactions in strongly polarizable systems. In this paper, we explore how the issues of MP2 for metallic, polarizable, and strongly correlated periodic systems can be ameliorated through regularization. To this end, two regularized second-order methods (including a new, size-extensive Brillouin-Wigner approach) are applied to the HEG, the one-dimensional Hubbard model, and the graphene-water interaction. We find that regularization consistently leads to improvements over the MP2 baseline and that different regularizers are appropriate for the various systems.

Electron-molecule collisions with explicit rovibrational resolution at MRCI level and using even tempered basis sets

A method for calculating the generalized oscillator strengths (GOSs) and differential cross section (DCS) with vibration and rotation resolution is presented. The importance of accounting for the rotational contribution is to be emphasized since it has not previously been considered in GOS calculations. Although largely neglected due to its small effect on various properties, the rotational resolution proved to be fundamental in the study of certain phenomena, such as the interference between rotational states in a molecule. As the general goal of this work is to obtain theoretical values comparable to high resolution experiments, special care was taken on the calculation of the electronic part of the scattering amplitude, particularly in what concerns the choice of the atomic basis set. Accordingly, even-tempered basis sets have proved to lead to good results. The helium atom was taken as a model system for this aspect of the problem. Then, GOS and DCS, for explicit vibrational and rotational transitions, were calculated for hydrogen and nitrogen molecules. For higher accuracy, a non-Franck-Condon approach was used to obtain transitions involving vibrational states. The resultant values have shown good agreement with the available experimental data.

Variational coupled cluster for ground and excited states

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes that are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations that may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC and we compare it with its (projected) traditional version in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear H₄ molecule and the continuous deformation of the square H₄ molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of variational pCCD (VpCCD) and traditional pCCD is gauged against their configuration interaction (CI) equivalent, known as doubly occupied CI, for reference Slater determinants made of ground- or excited-state Hartree-Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated.

Excited States from State-Specific Orbital-Optimized Pair Coupled Cluster

The pair coupled cluster doubles (pCCD) method (where the excitation manifold is restricted to electron pairs) has a series of interesting features. Among others, it provides ground-state energies very close to what is obtained with doubly occupied configuration interaction (DOCI), but with a polynomial cost (compared with the exponential cost of the latter). Here, we address whether this similarity holds for excited states by exploring the symmetric dissociation of the linear H4 molecule. When ground-state Hartree-Fock (HF) orbitals are employed, pCCD and DOCI excited-state energies do not match, a feature that is assigned to the poor HF reference. In contrast, by optimizing the orbitals at the pCCD level (oo-pCCD) specifically for each excited state, the discrepancies between pCCD and DOCI decrease by 1 or 2 orders of magnitude. Therefore, the pCCD and DOCI methodologies still provide comparable energies for excited states, but only if suitable, state-specific orbitals are adopted. We also assessed whether a pCCD approach could be used to directly target doubly excited states, without having to resort to the equation-of-motion (EOM) formalism. In our Δoo-pCCD model, excitation energies are extracted from the energy difference between separate oo-pCCD calculations for the ground state and the targeted excited state. For a set comprising the doubly excited states of CH+, BH, nitroxyl, nitrosomethane, and formaldehyde, we found that Δoo-pCCD provides quite accurate excitation energies, with root-mean-square deviations (with respect to full configuration interaction results) lower than those of CC3 and comparable to those of EOM-CCSDT, two methods with a much higher computational cost.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Quantum orbital-optimized unitary coupled cluster methods in the strongly correlated regime: Can quantum algorithms outperform their classical equivalents?

The Coupled Cluster (CC) method is used to compute the electronic correlation energy in atoms and molecules and often leads to highly accurate results. However, due to its single-reference nature, standard CC in its projected form fails to describe quantum states characterized by strong electronic correlations and multi-reference projective methods become necessary. On the other hand, quantum algorithms for the solution of many-electron problems have also emerged recently. The quantum unitary variant of CC (UCC) with singles and doubles (q-UCCSD) is a popular wavefunction Ansatz for the variational quantum eigensolver algorithm. The variational nature of this approach can lead to significant advantages compared to its classical equivalent in the projected form, in particular, for the description of strong electronic correlation. However, due to the large number of gate operations required in q-UCCSD, approximations need to be introduced in order to make this approach implementable in a state-of-the-art quantum computer. In this work, we evaluate several variants of the standard q-UCCSD Ansatz in which only a subset of excitations is included. In particular, we investigate the singlet and pair q-UCCD approaches combined with orbital optimization. We show that these approaches can capture the dissociation/distortion profiles of challenging systems, such as H₄, H₂O, and N₂ molecules, as well as the one-dimensional periodic Fermi-Hubbard chain. These results promote the future use of q-UCC methods for the solution of challenging electronic structure problems in quantum chemistry.

Hermitian second-order methods for excited electronic states: Unitary coupled cluster in comparison with algebraic-diagrammatic construction schemes

Employing an intermediate state representation (ISR) approach, Hermitian second-order methods for the calculation of electronic excitation energies are presented and compared in detail. These comprise the algebraic-diagrammatic construction scheme for the polarization propagator, a hybrid second-order ISR scheme based on traditional coupled-cluster theory as well as two similar approaches based on a unitary coupled-cluster (UCC) ansatz. Although in a strict perturbation-theoretical framework all prove to be identical, differences emerge when the corresponding converged cluster amplitudes are used and depending on how the similarity-transformed UCC Hamiltonian is evaluated. The resulting excitation energies, however, do not significantly differ for systems well described by means of perturbation theory.

What Levels of Coupled Cluster Theory Are Appropriate for Transition Metal Systems? A Study Using Near-Exact Quantum Chemical Values for 3d Transition Metal Binary Compounds

Transition metal compounds are traditionally considered to be challenging for standard quantum chemistry approximations like coupled cluster (CC) theory, which are usually employed to validate lower level methods like density functional theory (DFT). To explore this issue, we present a database of bond dissociation energies (BDEs) for 74 spin states of 69 diatomic species containing a 3d transition metal atom and a main group element, in the moderately sized def2-SVP basis. The presented BDEs appear to have an (estimated) 3σ error less than 1 kJ/mol relative to the exact solutions to the nonrelativistic Born-Oppenheimer Hamiltonian. These benchmark values were used to assess the performance of a wide range of standard single reference CC models, as the results should be beneficial for understanding the limitations of these models for transition metal systems. We find that interactions between metals and monovalent ligands like hydride and fluoride are well described by CCSDT. Similarly, CCSDTQ appears to be adequate for bonds between metals and nominally divalent ligands like oxide and sulfide. However, interactions with polyvalent ligands like nitride and carbide are more challenging, with even CCSDTQ(P)_Λ yielding errors on the scale of a few kJ/mol. We also find that many perturbative and iterative approximations to higher order terms either yield disappointing results or actually worsen the performance relative to the baseline low level CC method, indicating that complexity does not always guarantee accuracy.

Quantum Chemistry in the Age of Quantum Computing

Practical challenges in simulating quantum systems on classical computers have been widely recognized in the quantum physics and quantum chemistry communities over the past century. Although many approximation methods have been introduced, the complexity of quantum mechanics remains hard to appease. The advent of quantum computation brings new pathways to navigate this challenging and complex landscape. By manipulating quantum states of matter and taking advantage of their unique features such as superposition and entanglement, quantum computers promise to efficiently deliver accurate results for many important problems in quantum chemistry, such as the electronic structure of molecules. In the past two decades, significant advances have been made in developing algorithms and physical hardware for quantum computing, heralding a revolution in simulation of quantum systems. This Review provides an overview of the algorithms and results that are relevant for quantum chemistry. The intended audience is both quantum chemists who seek to learn more about quantum computing and quantum computing researchers who would like to explore applications in quantum chemistry.

Unrestricted treatment for the direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions

This work extends to the unrestricted orbital approach the procedure described in our previous report [Alcoba et al., J. Chem. Phys. 148, 024105 (2018)] for determining variationally the two-electron reduced density matrix arising from doubly occupied-configuration-interaction wave functions by imposing two- and three-index N-representability conditions. An analysis of the numerical results obtained in selected systems, from both restricted and unrestricted treatments, allows one to assess the performance of these methodologies as well as to show the influence of the P, Q, G, T1, and T2 positivity conditions. We highlight the satisfactory results obtained within the unrestricted scheme.

Orbital-invariant spin-extended approximate coupled-cluster for multi-reference systems

We present an approximate treatment of spin-extended coupled-cluster (ECC) based on the spin-projection of the broken-symmetry coupled-cluster (CC) ansatz. ECC completely eliminates the spin-contamination of unrestricted CC and is therefore expected to provide better descriptions of dynamical and static correlation effects, but introduces two distinct problems. The first issue is the emergence of non-terminating amplitude equations, which are caused by the de-excitation effects inherent in symmetry projection operators. In this study, we take a minimalist approach and truncate the Taylor series of the exponential ansatz at a certain order such that the approximation safely recovers the traditional CC without spin-projection. The second issue is that the nonlinear equations of ECC become underdetermined, although consistent, yielding an infinitude of solutions. This problem arises because of the redundancies in the excitation manifold, as is common in other multi-reference approaches. We remove the linear dependencies in ECC by employing an orthogonal projection manifold. We also propose an efficient solver for our method, in which the components are usually sparse but not diagonal-dominant. It is shown that our approach is rigorously orbital-invariant and provides more accurate results than its configuration interaction and linearized CC analogues for chemical systems.

Spin polynomial similarity transformation for repulsive Hamiltonians: interpolating between coupled cluster and spin-projected unrestricted Hartree–Fock

Our overarching goal is to be able to describe both weak and strong correlation with a single, computationally affordable method without sacrificing important qualities of the wavefunction, e.g. symmetries of the Hamiltonian.

H4: A challenging system for natural orbital functional approximations

The correct description of nondynamic correlation by electronic structure methods not belonging to the multireference family is a challenging issue. The transition of D(2h) to D(4h) symmetry in H4 molecule is among the most simple archetypal examples to illustrate the consequences of missing nondynamic correlation effects. The resurgence of interest in density matrix functional methods has brought several new methods including the family of Piris Natural Orbital Functionals (PNOF). In this work, we compare PNOF5 and PNOF6, which include nondynamic electron correlation effects to some extent, with other standard ab initio methods in the H4 D(4h)/D(2h) potential energy surface (PES). Thus far, the wrongful behavior of single-reference methods at the D(2h)-D(4h) transition of H4 has been attributed to wrong account of nondynamic correlation effects, whereas in geminal-based approaches, it has been assigned to a wrong coupling of spins and the localized nature of the orbitals. We will show that actually interpair nondynamic correlation is the key to a cusp-free qualitatively correct description of H4 PES. By introducing interpair nondynamic correlation, PNOF6 is shown to avoid cusps and provide the correct smooth PES features at distances close to the equilibrium, total and local spin properties along with the correct electron delocalization, as reflected by natural orbitals and multicenter delocalization indices.

Density cumulant functional theory from a unitary transformation: N-representability, three-particle correlation effects, and application to O4(+)

A new approach to density cumulant functional theory is developed that derives density cumulant N-representability conditions from an approximate Fock space unitary transformation. We present explicit equations for the third- and fourth-order two-particle cumulant N-representability, as well as the second-order contributions that depend on the connected three-particle density cumulant. These conditions are used to formulate the ODC-13 method and the non-iterative (λ3) correction that employ an incomplete description of the fourth-order two-particle cumulant N-representability and the second-order three-particle correlation effects, respectively. We perform an analysis of the ODC-13 N-representability description for the dissociation of H2 and apply the ODC-13 method and the (λ3) correction to diatomic molecules with multiple bond character and the symmetry-breaking tetraoxygen cation (O4(+)). For the O4(+) molecule, the vibrational frequencies of the ODC-13(λ3) method do not exhibit spatial symmetry breaking and are in a good agreement with the recent infrared photodissociation experiment. We report the O4(+) equilibrium structure, harmonic frequencies, and dissociation energy computed using ODC-13(λ3) with a diffuse, core-correlated aug-cc-pCVTZ basis set.