A spin-adapted coupled-cluster based linear response theory for double ionization potentials

Abstract

We developed in this article a spin-adapted formulation of the coupled-cluster based linear response theory (CC-LRT) for computing double-ionization potentials (DIPs), which may be experimentally observed by Auger spectroscopy. CC-LRT is a multireference generalization of the CC theory where the energy differences have no disconnected vacuum (core) diagrams, signifying core-extensivity. For the spin-adaptation of the CC-LRT equations for the singlet and triplet manifolds, we used the Young-Yamanouchi orthogonal spin-eigenfunctions. The orbital version of the CC-LRT equations are then automatically generated by the conjugate projection operators of Young-Yamanouchi spin functions. We illustrated the working of our spin-adaptation procedure by confining our CC-LRT equations to the space of 2h and 1p-3h ionized determinants. As numerical application of our formalism, we computed the Auger kinetic energies of HF and H₂O. We also analyzed the nature of size-extensivity of the DIPs generated by CC-LRT and showed explicitly that when the molecule is composed of two noninteracting fragments the computed DIPs are either DIPs of fragment A or B or a composite DIP depending on both A and B, which are just not sum of ionization potentials (IPs) of A and B. This analysis is done to underscore the fact that DIPs from CC-LRT is only core-extensive and not fully extensive.