Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods

We report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states' character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. The numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.