Treatment of spin-orbit coupling with internally contracted multireference coupled cluster theory.

We present the formalism for computing state-interaction matrix elements within the multistate extension of internally contracted multireference coupled-cluster theory. In particular, we focus on the determination of spin-orbit coupling matrix elements. A pilot implementation is presented and tested for the zero-field splitting of atomic 2P and 3P terms as well as for molecular 2Π terms. We also investigate the impact of the non-Hermiticity of the underlying coupled-cluster theory on spin-orbit couplings between non-equivalent states. For this, we consider the coupling between the 3Σ- ground state and the 1Π excited state of NH for several interatomic distances. Only a very small asymmetry is found for this case, but it is also seen that the present theory has to be carefully revisited with respect to the treatment of the response of the reference coefficients in order to avoid artefacts.