  Heavy atom compounds represent a challenge for computational chemistry, due
to the need for simultaneous treatment of relativistic and correlation effects.
Often such systems exhibit also strong correlation which hampers the
application of perturbation theory or single-reference coupled cluster (CC)
methods. As a viable alternative, we have proposed to externally correct 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. {\bf
152}, 174107 (2020)] we have 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 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 dissociation of the
N-U bond, where the system exhibits a 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.
