Seniority-based coupled cluster theory

Toward Reliable Dipole Moments without Single Excitations: The Role of Orbital Rotations and Dynamical Correlation

The dipole moment is a crucial molecular property linked to a molecular system's bond polarity and overall electronic structure. To that end, the electronic dipole moment, which results from the electron density of a system, is often used to assess the accuracy and reliability of new electronic structure methods. This work analyses electronic dipole moments computed with the pair coupled cluster doubles (pCCD) ansätze and its linearized coupled cluster (pCCD-LCC) corrections using the canonical Hartree-Fock and pCCD-optimized (localized) orbital bases. The accuracy of pCCD-based dipole moments is assessed against experimental and CCSD(T) reference values using relaxed and unrelaxed density matrices and different basis set sizes. Our test set comprises molecules of various bonding patterns and electronic structures, exposing pCCD-based methods to a wide range of electron correlation effects. Additionally, we investigate the performance of pCCD-in-DFT dipole moments of some model complexes. Finally, our work indicates the importance of orbital relaxation in the pCCD model and shows the limitations of the linearized couple cluster corrections in predicting electronic dipole moments of multiple-bonded systems. Most importantly, pCCD with a linearized CCD correction can reproduce the dipole moment surfaces in singly bonded molecules, which are comparable to the multireference ones.

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.

Coupled cluster-inspired geminal wavefunctions

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

Correlated pair ansatz with a binary tree structure

We develop an efficient algorithm to implement the recently introduced binary tree state (BTS) ansatz on a classical computer. BTS allows a simple approximation to permanents arising from the computationally intractable antisymmetric product of interacting geminals and respects size-consistency. We show how to compute BTS overlap and reduced density matrices efficiently. We also explore two routes for developing correlated BTS approaches: Jastrow coupled cluster on BTS and linear combinations of BT states. The resulting methods show great promise in benchmark applications to the reduced Bardeen-Cooper-Schrieffer Hamiltonian and the one-dimensional XXZ Heisenberg Hamiltonian.

Single Reference Treatment of Strongly Correlated H4 and H10 Isomers with Richardson-Gaudin States

Richardson-Gaudin (RG) states are employed as a variational wave function ansatz for strongly correlated isomers of H4 and H10. In each case, a single RG state describes the seniority-zero sector quite well. Simple natural orbital functionals offer a cheap and reasonable approximation of the outstanding weak correlation in the seniority-zero sector, while systematic improvement is achieved by performing a configuration interaction in terms of RG states.

Seniority and Hierarchy Configuration Interaction for Radicals and Excited States

Hierarchy configuration interaction (hCI) has recently been introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number and has been shown to efficiently recover both dynamic and static correlations for closed-shell molecular systems [ J. Phys. Chem. Lett. 2022, 13, 4342]. Here we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and excited states in a state-specific way. We gauge this route against excitation-based CI (eCI) and seniority-based CI (sCI) by evaluating how different ground-state properties of radicals converge to the full CI limit. We find that hCI outperforms or matches eCI, whereas sCI is far less accurate, in line with previous observations for closed-shell molecules. Employing second-order Epstein-Nesbet (EN2) perturbation theory as a correction significantly accelerates the convergence of hCI and eCI. We further explore various hCI and sCI models to calculate the excitation energies of closed- and open-shell systems. Our results underline that the choice of both the reference determinant and the set of orbitals drives the fine balance between correlation of ground and excited states. State-specific hCI2 and higher-order models perform similarly to their eCI counterparts, whereas lower orders of hCI deliver poor results unless supplemented by the EN2 correction, which substantially improves their accuracy. In turn, sCI1 produces decent excitation energies for radicals, encouraging the development of related seniority-based coupled-cluster methods.

Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials

We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.

Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure

A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for ground state energy with high accuracy. While classical quantum chemistry is a relatively mature field, the accurate and scalable prediction of strongly correlated states found, e.g., in bond breaking and polynuclear transition metal compounds remains an open problem. Within the context of a variational quantum eigensolver, we propose a new family of ansatzes which provides a more physically appropriate description of strongly correlated electrons than a unitary coupled cluster with single and double excitations (qUCCSD), with vastly reduced quantum resource requirements. Specifically, we present a set of local approximations to the unitary cluster Jastrow wavefunction motivated by Hubbard physics. As in the case of qUCCSD, exactly computing the energy scales factorially with system size on classical computers but polynomially on quantum devices. The local unitary cluster Jastrow ansatz removes the need for SWAP gates, can be tailored to arbitrary qubit topologies (e.g., square, hex, and heavy-hex), and is well-suited to take advantage of continuous sets of quantum gates recently realized on superconducting devices with tunable couplers. The proposed family of ansatzes demonstrates that hardware efficiency and physical transparency are not mutually exclusive; indeed, chemical and physical intuition regarding electron correlation can illuminate a useful path towards hardware-friendly quantum circuits.

Excited states with pair coupled cluster doubles tailored coupled cluster theory

It is known that some non-dynamic effects of electron correlation can be included in coupled cluster theory using a tailoring technique that separates the effects of non-dynamic and dynamic correlations. Recently, the simple pCCD (pair coupled cluster doubles) wavefunction was shown to provide good results for some non-dynamic correlation problems, such as bond-breaking, in a spin-adapted way with no active space selection. In this paper, we report a study of excited states using "tailored coupled cluster singles and doubles," to attempt to use pCCD as a kernel for more complete coupled-cluster singles and doubles (CCSD) results for excited states. Several excited states are explored from those primarily due to single excitations to those dominated by doubly excited states and from singlet-triplet splittings for some diradical states. For the first two situations, tailored pCCD-TCCSD offers no improvement over equation of motion-CCSD. However, when we explore the singlet-triplet gap of diradical molecules that are manifestly multi-reference, a pCCD kernel provides improved results, particularly with generalized valence bond orbitals.

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.

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations

In this perspective, the various measures of electron correlation used in wave function theory, density functional theory and quantum information theory are briefly reviewed. We then focus on a more traditional metric based on dominant weights in the full configuration solution and discuss its behavior with respect to the choice of the N-electron and the one-electron basis. The impact of symmetry is discussed, and we emphasize that the distinction among determinants, configuration state functions and configurations as reference functions is useful because the latter incorporate spin-coupling into the reference and should thus reduce the complexity of the wave function expansion. The corresponding notions of single determinant, single spin-coupling and single configuration wave functions are discussed and the effect of orbital rotations on the multireference character is reviewed by analyzing a simple model system. In molecular systems, the extent of correlation effects should be limited by finite system size and in most cases the appropriate choices of one-electron and N-electron bases should be able to incorporate these into a low-complexity reference function, often a single configurational one.

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen-Cooper-Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. The method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree-Fock based electronic structure methods.

A configuration interaction correction on top of pair coupled cluster doubles

Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N₂ and F₂ dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.

Fanpy: A python library for prototyping multideterminant methods in ab initio quantum chemistry

Fanpy is a free and open-source Python library for developing and testing multideterminant wavefunctions and related ab initio methods in electronic structure theory. The main use of Fanpy is to quickly prototype new methods by making it easier to convert the mathematical formulation of a new wavefunction ansätze to a working implementation. Fanpy is designed based on our recently introduced Flexible Ansatz for N-electron Configuration Interaction (FANCI) framework, where multideterminant wavefunctions are represented by their overlaps with Slater determinants of orthonormal spin-orbitals. In the simplest case, a new wavefunction ansatz can be implemented by simply writing a function for evaluating its overlap with an arbitrary Slater determinant. Fanpy is modular in both implementation and theory: the wavefunction model, the system's Hamiltonian, and the choice of objective function are all independent modules. This modular structure makes it easy for users to mix and match different methods and for developers to quickly explore new ideas. Fanpy is written purely in Python with standard dependencies, making it accessible for various operating systems. In addition, it adheres to principles of modern software development, including comprehensive documentation, extensive testing, quality assurance, and continuous integration and delivery protocols. This article is considered to be the official release notes for the Fanpy library.

Reduced density matrices/static correlation functions of Richardson–Gaudin states without rapidities

Seniority-zero geminal wavefunctions are known to capture bond-breaking correlation. Among this class of wavefunctions, Richardson–Gaudin states stand out as they are eigenvectors of a model Hamiltonian. This provides a clear physical picture, clean expressions for reduced density matrix (RDM) elements, and systematic improvement (with a complete set of eigenvectors). Known expressions for the RDM elements require the computation of rapidities, which are obtained by first solving for the so-called eigenvalue based variables (EBV) and then root-finding a Lagrange interpolation polynomial. In this paper, we obtain expressions for the RDM elements directly in terms of the EBV. The final expressions can be computed at the same cost as the rapidity expressions. Therefore, except, in particular, circumstances, it is entirely unnecessary to compute rapidities at all. The RDM elements require numerically inverting a matrix, and while this is usually undesirable, we demonstrate that it is stable, except when there is degeneracy in the single-particle energies. In such cases, a different construction would be required.

Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry

In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.

Resolving the π-assisted U-N σf-bond formation using quantum information theory

We model the potential energy profiles of the UO₂ (NCO)Cl₂^- → NUOCl₂^- + CO₂ reaction pathway [Y. Gong, V. Vallet, M. del Carmen Michelini, D. Rios and J. K. Gibson, J. Phys. Chem. A, 2014, 118, 325-330] using different pair coupled-cluster doubles (pCCD) methods. Specifically, we focus on pCCD and pCCD-tailored coupled cluster models in predicting relative energies for the various intermediates and transition states along the reaction coordinate. Furthermore, we augment our study on energetics with an orbital-pair correlation analysis of the complete reaction pathway that features two distinct paths. Our analysis of orbital correlations sheds new light on the formation and breaking of respective bonds between the uranium, oxygen, and nitrogen atoms along the reaction coordinates where the "yl" bond is broken and a nitrido compound formed. Specifically, the strengthening of the U-N σ_f-bond is assisted by a π-type interaction that is delocalized over the C-N-U backbone of the UO₂ (NCO)Cl₂^- complex.

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.

Density matrices of seniority-zero geminal wavefunctions

Scalar products and density matrix elements of closed-shell pair geminal wavefunctions are evaluated directly in terms of the pair amplitudes, resulting in an analog of Wick’s theorem for fermions or bosons. This expression is, in general, intractable, but it is shown how it becomes feasible in three distinct ways for Richardson–Gaudin (RG) states, the antisymmetrized geminal power, and the antisymmetrized product of strongly orthogonal geminals. Dissociation curves for hydrogen chains are computed with off-shell RG states and the antisymmetrized product of interacting geminals. Both are near exact, suggesting that the incorrect results observed with ground state RG states (a local maximum rather than smooth dissociation) may be fixable using a different RG state.

Near-exact treatment of seniority-zero ground and excited states with a Richardson-Gaudin mean-field

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian, Richardson-Gaudin (RG) states, are used as a variational wavefunction ansatz for strongly correlated electronic systems. These states are geminal products whose coefficients are solutions of non-linear equations. Previous results showed an un-physical apparent avoided crossing in ground state dissociation curves for hydrogen chains. In this paper, it is shown that each seniority-zero state of the molecular Coulomb Hamiltonian corresponds directly to an RG state. However, the seniority-zero ground state does not correspond to the ground state of a reduced BCS Hamiltonian. The difficulty is in choosing the correct RG state. The systems studied showed a clear choice, and we expect that it should always be possible to reason physically which state to choose.

Hierarchy Configuration Interaction: Combining Seniority Number and Excitation Degree

We propose a novel partitioning of the Hilbert space, hierarchy configuration interaction (hCI), where the excitation degree (with respect to a given reference determinant) and the seniority number (i.e., the number of unpaired electrons) are combined in a single hierarchy parameter. The key appealing feature of hCI is that each hierarchy level accounts for all classes of determinants whose number shares the same scaling with system size. By surveying the dissociation of multiple molecular systems, we found that the overall performance of hCI usually exceeds or, at least, parallels that of excitation-based CI. For higher orders of hCI and excitation-based CI, the additional computational burden related to orbital optimization usually does not compensate the marginal improvements compared with results obtained with Hartree-Fock orbitals. The exception is orbital-optimized CI with single excitations, a minimally correlated model displaying the qualitatively correct description of single bond breaking at a very modest computational cost.

Assessing the Accuracy of Tailored Coupled Cluster Methods Corrected by Electronic Wave Functions of Polynomial Cost

Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N₂, F₂, C₂, CN⁺, CO, BN, BO⁺, and Cr₂) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.

Open-shell extensions to closed-shell pCCD

The pair coupled cluster doubles (pCCD) ansatz represents an inexpensive but accurate single-reference method to describe multi-reference problems. By construction, pCCD remains, however, applicable to closed-shell systems. For the first time, we present extensions to pCCD that allow us to target open-shell molecules with up to 4 unpaired electrons. Although requiring only modest computational cost, our methods approach chemical accuracy for some challenging cases, while their performance is comparable to more expensive models like DMRG or CCSD(T).

Accurate full configuration interaction correlation energy estimates for five- and six-membered rings

Following our recent work on the benzene molecule [P.-F. Loos, Y. Damour, and A. Scemama, J. Chem. Phys. 153, 176101 (2020)], motivated by the blind challenge of Eriksen et al. [J. Phys. Chem. Lett. 11, 8922 (2020)] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for 12 five- and six-membered ring molecules (cyclopentadiene, furan, imidazole, pyrrole, thiophene, benzene, pyrazine, pyridazine, pyridine, pyrimidine, s-tetrazine, and s-triazine) in the standard correlation-consistent double-ζ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with an estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction calculations performed with the configuration interaction using a perturbative selection made iteratively algorithm. Having at our disposal these accurate reference energies, the respective performance and convergence properties of several popular and widely used families of single-reference quantum chemistry methods are investigated. In particular, we study the convergence properties of (i) the Møller-Plesset perturbation series up to fifth-order (MP2, MP3, MP4, and MP5), (ii) the iterative approximate coupled-cluster series CC2, CC3, and CC4, and (iii) the coupled-cluster series CCSD, CCSDT, and CCSDTQ. The performance of the ground-state gold standard CCSD(T) as well as the completely renormalized CC model, CR-CC(2,3), is also investigated. We show that MP4 provides an interesting accuracy/cost ratio, while MP5 systematically worsens the correlation energy estimates. In addition, CC3 outperforms CCSD(T) and CR-CC(2,3), as well as its more expensive parent CCSDT. A similar trend is observed for the methods including quadruple excitations, where the CC4 model is shown to be slightly more accurate than CCSDTQ, both methods providing correlation energies within 2 millihartree of the FCI limit.

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.

GQCP: The Ghent Quantum Chemistry Package

The Ghent Quantum Chemistry Package (GQCP) is an open-source electronic structure software package that aims to provide an intuitive and expressive software framework for electronic structure software development. Its high-level interfaces (accessible through C++ and Python) have been specifically designed to correspond to theoretical concepts, while retaining access to lower-level intermediates and allowing structural run-time modifications of quantum chemical solvers. GQCP focuses on providing quantum chemical method developers with the computational "building blocks" that allow them to flexibly develop proof of principle implementations for new methods and applications up to the level of two-component spinor bases.

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 H₄ 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.

Transition density matrices of Richardson-Gaudin states

Recently, ground state eigenvectors of the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian, Richardson-Gaudin (RG) states, have been employed as a wavefunction ansatz for strong correlation. This wavefunction physically represents a mean-field of pairs of electrons (geminals) with a constant pairing strength. To move beyond the mean-field, one must develop the wavefunction on the basis of all the RG states. This requires both practical expressions for transition density matrices and an idea of which states are most important in the expansion. In this contribution, we present expressions for the transition density matrix elements and calculate them numerically for half-filled picket-fence models (reduced BCS models with constant energy spacing). There are no Slater-Condon rules for RG states, though an analog of the aufbau principle proves to be useful in choosing which states are important.

Doubly Occupied Pair Coupled Cluster F12 Approach

Inspired by reports of the good performance of the doubly occupied pair coupled cluster (pCCD) theory in describing static electron correlation, we have introduced and implemented a variant thereof that includes single excitations and explicitly treats the dynamic electron correlation using the F12 methodology (pCCSD-F12). This drastically reduces the computation scaling with respect to the standard method using the full double-excitation operator (CCSD-F12). Slater-type geminals as a correlation factor, together with fixed cusp conditions, were used, which is known as the SP-ansatz. For sample model systems, we have investigated the performance of reference states constructed from either canonical or localized molecular orbitals. Finaly, the employment of Brueckner orbitals has been tested, which causes the single excitations to naturally vanish from the wave function expansion (B-pCCD-F12). Our test systems include different-sized rings of hydrogen atoms and dissociation curves for small molecules such as HF, N₂, and CO₂; and comparison with CCSD-F12 is presented for a series of reaction enthalpies.

Reduced density matrices of Richardson-Gaudin states in the Gaudin algebra basis

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer Hamiltonian have recently been employed as a variational wavefunction ansatz in quantum chemistry. This wavefunction is a mean-field of pairs of electrons (geminals). In this contribution, we report optimal expressions for their reduced density matrices in both the original physical basis and the basis of the Richardson-Gaudin pairs. Physical basis expressions were originally reported by Gorohovsky and Bettelheim [Phys. Rev. B 84, 224503 (2011)]. In each case, the expressions scale like O(N⁴), with the most expensive step being the solution of linear equations. Analytic gradients are also reported in the physical basis. These expressions are an important step toward practical mean-field methods to treat strongly correlated electrons.

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.

Downfolding of many-body Hamiltonians using active-space models: Extension of the sub-system embedding sub-algebras approach to unitary coupled cluster formalisms

In this paper, we discuss the extension of the recently introduced subsystem embedding subalgebra coupled cluster (SES-CC) formalism to unitary CC formalisms. In analogy to the standard single-reference SES-CC formalism, its unitary CC extension allows one to include the dynamical (outside the active space) correlation effects in an SES induced complete active space (CAS) effective Hamiltonian. In contrast to the standard single-reference SES-CC theory, the unitary CC approach results in a Hermitian form of the effective Hamiltonian. Additionally, for the double unitary CC (DUCC) formalism, the corresponding CAS eigenvalue problem provides a rigorous separation of external cluster amplitudes that describe dynamical correlation effects-used to define the effective Hamiltonian-from those corresponding to the internal (inside the active space) excitations that define the components of eigenvectors associated with the energy of the entire system. The proposed formalism can be viewed as an efficient way of downfolding many-electron Hamiltonian to the low-energy model represented by a particular choice of CAS. In principle, this technique can be extended to any type of CAS representing an arbitrary energy window of a quantum system. The Hermitian character of low-dimensional effective Hamiltonians makes them an ideal target for several types of full configuration interaction type eigensolvers. As an example, we also discuss the algebraic form of the perturbative expansions of the effective DUCC Hamiltonians corresponding to composite unitary CC theories and discuss possible algorithms for hybrid classical and quantum computing. Given growing interest in quantum computing, we provide energies for H₂ and Be systems obtained with the quantum phase estimator algorithm available in the Quantum Development Kit for the approximate DUCC Hamiltonians.

Polynomial-product states: A symmetry-projection-based factorization of the full coupled cluster wavefunction in terms of polynomials of double excitations

Our goal is to remedy the failure of symmetry-adapted coupled-cluster theory in the presence of strong correlation. Previous work along these lines has taken us from a diagram-level analysis of the coupled-cluster equations to an understanding of the collective modes which can occur in various channels of the coupled-cluster equations to the exploration of non-exponential wavefunctions in efforts to combine coupled-cluster theory with symmetry projection. In this manuscript, we extend these efforts by introducing a new, polynomial product wavefunction ansatz that incorporates information from symmetry projection into standard coupled-cluster theory in a way that attempts to mitigate the effects of the lack of size extensivity and size consistency characteristic of symmetry-projected methods. We describe the new approach in detail within the context of our previous efforts, explore some illustrative calculations, and consider one route for reducing the computational cost of the new method.

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer

In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO–LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.

Assessing the accuracy of simplified coupled cluster methods for electronic excited states in f0 actinide compounds

We scrutinize the performance of different variants of equation of motion coupled cluster (EOM-CC) methods to predict electronic excitation energies and excited state potential energy surfaces in closed-shell actinide species.