Excited states with pair coupled cluster doubles tailored coupled cluster theory

Complete Active Space Iterative Coupled Cluster Theory

In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through a similarity transformation of the Hamiltonian with those CC amplitudes that are not encompassed by the active space. We denote this approach as the complete active space iterative coupled cluster (CASiCC) ansatz. We investigate its efficiency and accuracy in the singles and doubles approximation by studying the prototypical molecules H4, H8, H2O, and N2. Our results demonstrate that CASiCC systematically improves on the single-reference CCSD and the externally corrected CCSD methods across entire potential energy curves while retaining modest computational costs. However, the tailored coupled cluster method shows superior performance in the strong correlation regime, suggesting that its accuracy is based on error compensation. We find that the iterative versions of externally corrected and tailored coupled cluster methods converge to the same results.

ModelHamiltonian: A Python-scriptable library for generating 0-, 1-, and 2-electron integrals

ModelHamiltonian is a free, open source, and cross-platform Python library designed to express model Hamiltonians, including spin-based Hamiltonians (Heisenberg and Ising models) and occupation-based Hamiltonians (Pariser-Parr-Pople, Hubbard, and Hückel models) in terms of 1- and 2-electron integrals, so that these systems can be easily treated by traditional quantum chemistry software programs. ModelHamiltonian was originally intended to facilitate the testing of new electronic structure methods using HORTON but emerged as a stand-alone research tool that we recognize has wide utility, even in an educational context. ModelHamiltonian is written in Python and adheres to modern principles of software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. While we anticipate that most users will use ModelHamiltonian as a Python library, we include a graphical user interface so that models can be built without programming, based on connectivity/parameters inferred from, for example, a SMILES string. We also include an interface to ChatGPT so that users can specify a Hamiltonian in plain language (without learning ModelHamiltonian's vocabulary and syntax). This article marks the official release of the ModelHamiltonian library, showcasing its functionality and scope.

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively "cheap" and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

Delving into the catalytic mechanism of molybdenum cofactors: a novel coupled cluster study

In this work, we use modern electronic structure methods to model the catalytic mechanism of different variants of the molybdenum cofactor (Moco). We investigate the dependence of various Moco model systems on structural relaxation and the importance of environmental effects for five critical points along the reaction coordinate with the DMSO and NO3- substrates. Furthermore, we scrutinize the performance of various coupled-cluster approaches for modeling the relative energies along the investigated reaction paths, focusing on several pair coupled cluster doubles (pCCD) flavors and conventional coupled cluster approximations. Moreover, we elucidate the Mo-O bond formation using orbital-based quantum information measures, which highlight the flow of σM-O bond formation and σN/S-O bond breaking. Our study shows that pCCD-based models are a viable alternative to conventional methods and offer us unique insights into the bonding situation along a reaction coordinate. Finally, this work highlights the importance of environmental effects or changes in the core and, consequently, in the model itself to elucidate the change in activity of different Moco variants.

Tailored and Externally Corrected Coupled Cluster with Quantum Inputs

We propose to use wave function overlaps obtained from a quantum computer as inputs for the classical split-amplitude techniques, tailored and externally corrected coupled cluster, to achieve balanced treatment of static and dynamic correlation effects in molecular electronic structure simulations. By combining insights from statistical properties of matchgate shadows, which are used to measure quantum trial state overlaps, with classical correlation diagnostics, we can provide quantum resource estimates well into the classically no longer exactly solvable regime. We find that rather imperfect wave functions and remarkably low shot counts are sufficient to cure qualitative failures of plain coupled cluster singles doubles and to obtain chemically precise dynamic correlation energy corrections. We provide insights into which wave function preparation schemes have a chance of yielding quantum advantage, and we test our proposed method using overlaps measured on Google's Sycamore device.

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces

This paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC 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.

Two determinant distinguishable cluster

A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.

Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry

Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled ab initio quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations. The best that one can do in a basis is the 'full CI,' the exact solution of the non-relativistic Schrödinger equation or, if need be, the relativistic Dirac equation. With due regard to converging the basis set and adequate consideration of higher clusters and relativity in a calculation, virtually predictive results can be obtained. But in addition to its numerical performance, coupled-cluster theory also offers a conceptually new, many-body foundation for the theory that should be appreciated by all practitioners. The latter is emphasized in this perspective, leading to the 'evolution toward simplicity' in the title. The ultimate theory will benefit from the several features that are uniquely exact in coupled-cluster theory and its equation-of-motion (EOM-CC) extensions.

Variational quantum eigensolver for closed-shell molecules with non-bosonic corrections

The realization of quantum advantage with noisy-intermediate-scale quantum (NISQ) machines has become one of the major challenges in computational sciences. Maintaining coherence of a physical system with more than ten qubits is a critical challenge that motivates research on compact system representations to reduce algorithm complexity. Toward this end, the variational quantum eigensolver (VQE) used to perform quantum simulations is considered to be one of the most promising algorithms for quantum chemistry in the NISQ era. We investigate reduced mapping of one spatial orbital to a single qubit to analyze the ground state energy in a way that the Pauli operators of qubits are mapped to the creation/annihilation of singlet pairs of electrons. To include the effect of non-bosonic (or non-paired) excitations, we introduce a simple correction scheme in the electron correlation model approximated by the geometrical mean of the bosonic (or paired) terms. Employing it in a VQE algorithm, we assess ground state energies of H2O, N2, and Li2O in good agreement with full configuration interaction (FCI) models respectively, using only 6, 8, and 12 qubits with quantum gate depths proportional to the squares of the qubit counts. With the adopted seniority-zero approximation that uses only one half of the qubit counts of a conventional VQE algorithm, we find that our non-bosonic correction method reaches reliable quantum chemistry simulations at least for the tested systems.