Equation-of-motion coupled cluster method for high spin double electron attachment calculations

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.

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.

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.

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.

Non-Valence Anions of Pyridine and the Diazines

The dipole-bound anions of pyridine, pyridazine, and pyrimidine are characterized using equation of motion coupled cluster singles and doubles calculations. These calculations predict that the anions of pyridine, pyrimidine and pyridazine are bound in the Born-Oppenheimer approximation by 0.05, 0.8, and 19.0 meV, respectively. The binding energies of pyrimidine and pyridazine are large enough that the anions will remain bound even when allowing for corrections to the Born-Oppenheimer approximation, while that of pyridine is a borderline case. We were unable to find a stable non-valence correlation-bound anion for pyrazine, which has a zero dipole moment.

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).

Equation-of-motion coupled-cluster theory for double electron attachment with spin-orbit coupling

We report implementation of the equation-of-motion coupled-cluster (EOM-CC) method for double electron-attachment (DEA) with spin-orbit coupling (SOC) at the CC singles and doubles (CCSD) level using a closed-shell reference in this work. The DEA operator employed in this work contains two-particle and three-particle one-hole excitations, and SOC is included in post-Hartree-Fock treatment. Time-reversal symmetry and spatial symmetry are exploited to reduce computational cost. The EOM-DEA-CCSD method with SOC allows us to investigate SOC effects of systems with two-unpaired electrons. According to our results on atoms, double ionization potentials (DIPs), excitation energies (EEs), and SO splittings of low-lying states are calculated reliably using the EOM-DEA-CCSD method with SOC. Its accuracy is usually higher than that of EOM-CCSD for EEs or DIPs if the same target can be reached from single excitations by choosing a proper closed-shell reference. However, performance of the EOM-DEA-CCSD method with SOC on molecules is not as good as that for atoms. Bond lengths for the ground and the several lowest excited states of GaH, InH, and TlH are underestimated pronouncedly, although reasonable EEs are obtained, and splittings of the ³Σ^- state from the π² configuration are calculated to be too small with EOM-DEA-CCSD.

Splittings of d8 configurations of late-transition metals with EOM-DIP-CCSD and FSCCSD methods

Multireference methods are usually required for transition metal systems due to the partially filled d electrons. In this work, the single-reference equation-of-motion coupled-cluster method at the singles and doubles level for double ionization potentials (EOM-DIP-CCSD) is employed to calculate energies of states from the d⁸ configuration of late-transition metal atoms starting from a closed-shell reference. Its results are compared with those from the multireference Fock-space coupled-cluster method at the CCSD level (FSCCSD) for DIP from the same closed-shell reference. Both scalar-relativistic effects and spin-orbit coupling are considered in these calculations. Compared with all-electron FSCCSD results with four-component Dirac-Coulomb Hamiltonian, FSCCSD with relativistic effective core potentials can provide reasonable results, except for atoms with unstable reference. Excitation energies for states in the (n - 1)d⁸ns² configuration are overestimated pronouncedly with these two methods, and this overestimation is more severe than those in the (n - 1)d⁹ns¹ configuration. Error of EOM-CCSD on these excitation energies is generally larger than that of FSCCSD. On the other hand, relative energies of most of the states in the d⁸ configuration with respect to the lowest state in the same configuration are predicted reliably with EOM-DIP-CCSD, except for the ³P₀ state of Hg²⁺ and states in Ir⁺. FSCCSD can provide reasonable relative energies for the several lowest states, while its error tends to be larger for higher states.

Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states

This article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl₂, Br₂, HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Furthermore, the obtained results for the molecular systems agree well with the available experimental values.

Equation-of-motion coupled cluster method for the description of the high spin excited states

The equation-of-motion (EOM) coupled cluster (CC) approach in the version applicable for the excitation energy(EE) calculations has been formulated for high spin components. The EE-EOM-CC scheme based on the restricted Hartree-Fock reference and standard amplitude equations as used in the Davidson diagonalization procedure yields the singlet states. The triplet and higher spin components require separate amplitude equations. In the case of quintets, the relevant equations are much simpler and easier to solve. Out of 26 diagrammatic terms contributing to the R1 and R2 singlet equations in the case of quintets, only R2 operator survives with 5 diagrammatic terms present. In addition all terms engaging three body elements of the similarity transformed Hamiltonian disappear. This indicates a substantial simplification of the theory. The implemented method has been applied to the pilot study of the excited states of the C2 molecule and quintet states of C and Si atoms.