Equation-of-motion coupled-cluster method for ionized states with spin-orbit coupling

We report implementation of the equation of motion coupled-cluster approach for ionized states (EOMIP-CC) with spin-orbit coupling (SOC) using closed-shell state as reference in this work. Ionization potentials (IPs) are calculated in the ionized 1h and 2h1p space with EOM at the CC singles (CCS) as well as the CC singles and doubles levels (CCSD). In this EOMIP-CC approach, SOC is included either in both the CC and EOM steps or only in the EOM step. It should be noted that IPs provided by the EOMIP-CC approach with SOC included only in the EOM step are not size-intensive. Time-reversal symmetry and spatial symmetry are exploited for D(2h) and its subgroups to reduce computational effort. All these approaches have been shown to be able to afford acceptable estimates for SOC splittings. The EOMIP-CCSD with SOC included only in the EOM step can provide reasonable IPs for systems containing up to 5th row elements. On the other hand, the EOMIP-CCS approach with SOC included in both CC and EOM steps could not predict a bounded (2)∑(g) (+) state for I(2) (+) and should be used with care.