Simple and efficient computational strategies for calculating orbital energies and pair-orbital energies from pCCD-based methods.

We introduce affordable computational strategies for calculating orbital and pair-orbital energies in atomic and molecular systems. Our methods are based on the pair Coupled Cluster Doubles (pCCD) ansatz and its orbital-optimized variant. The computed orbital and pair-orbital energies are then subsequently used to approximate ionization potentials (IPs), electron affinities (EAs), the resulting charge gaps, double ionization potentials (DIPs), and double electron affinities (DEAs). Our methodology builds on the standard Koopmans' theorem and refines it for a pCCD-based wave function. Furthermore, we incorporate pCCD electron correlation effects into the model utilizing canonical Hartree-Fock or natural pCCD-optimized orbitals. The latter represents a diagonal approximation to the (D)IP/D(EA) equation of motion pCCD models. We benchmarked our newly developed models against theoretical and available experimental data for selected atoms in various basis set sizes and a set of 24 organic acceptor molecules. Our numerical results show that the Koopmans' approach based on pCCD natural orbitals provides a balanced treatment of occupied and virtual orbital energies, resulting in reliable predictions of charge gaps at a low computational cost.