Abstract
Heavy atom compounds represent a challenge for computational chemistry due to the need for simultaneous treatment of relativistic and correlation effects. Often such systems also exhibit strong correlation, which hampers the application of perturbation theory or single-reference coupled cluster (CC) methods. As a viable alternative, we have proposed externally correcting the CC method using the density matrix renormalization group (DMRG) wave functions, yielding the DMRG-tailored CC method. In a previous paper [J. Chem. Phys. 2020, 152, 174107], we reported a first implementation of this method in the relativistic context, which was restricted to molecules with real double group symmetry. In this work, we present a fully general implementation of the method, covering complex and quaternion double groups as well. The 4c-TCC method thus becomes applicable to polyatomic molecules, including heavy atoms. For the assessment of the method, we performed calculations of the chiral uranium compound NUHFI, which was previously studied in the context of the enhancement of parity violation effects. In particular, we performed calculations of a cut of the potential energy surface of this molecule along the stretching of the N-U bond, where the system exhibits strong multireference character. Since there are no experimental data for NUHFI, we have performed also an analogous study of the (more symmetric) NUF(3) molecule, where the vibrational frequency of the N-U bond can be compared with spectroscopic data.