Two-component relativistic equation-of-motion coupled cluster for electron ionization

Abstract

We present an implementation of the relativistic ionization-potential (IP) equation-of-motion coupled-cluster (EOMCC) with up to 3-hole–2-particle (3h2p) excitations that makes use of the molecular mean-field exact two-component framework and the full Dirac–Coulomb–Breit Hamiltonian. The closed-shell nature of the reference state in an X2C-IP-EOMCC calculation allows for accurate predictions of spin–orbit splittings in open-shell molecules without breaking degeneracies, as would occur in an excitation-energy EOMCC calculation carried out directly on an unrestricted open-shell reference. We apply X2C-IP-EOMCC to the ground and first excited states of the HCCX+ (X = Cl, Br, I) cations, where it is demonstrated that a large basis set (i.e., quadruple-zeta quality) and 3h2p correlation effects are necessary for accurate absolute energetics. The maximum error in calculated adiabatic IPs is on the order of 0.1 eV, whereas spin–orbit splittings themselves are accurate to ≈0.01 eV, as compared to experimentally obtained values.