A near-linear scaling equation of motion coupled cluster method for ionized states

In this work, a domain-based local pair natural orbital (DLPNO) version of the equation of motion coupled cluster theory with single and double excitations for ionization potentials (IP-EOM-CCSD) equations has been formulated and implemented. The method uses ground state localized occupied and pair natural virtual orbitals and applies the DLPNO machinery to arrive at a linear scaling implementation of the IP-EOM-CCSD method. The accuracy of the method is controllable using ground state truncation parameters. Using default thresholds, the method predicts ionization potential (IP) values with good accuracy (mean absolute error of 0.08 eV). We demonstrate that our code can be used to compute IP values for systems with more than 1000 atoms and 10 000 basis functions.