A Comparison of Relativistic Coupled Cluster and Equation of Motion Coupled Cluster Quadratic Response Theory

We present the implementation of relativistic coupled cluster quadratic response theory (QR-CC), following our development of relativistic equation of motion coupled cluster quadratic response theory (QR-EOMCC) [X. Yuan et al., <i>J. Chem. Theory Comput.</i> <b>2023</b>, 19, 9248-9259]. These codes, which can be used in combination with relativistic (2- and 4-component based) as well as nonrelativistic Hamiltonians, are capable of treating both static and dynamic perturbations for electric and magnetic operators. We have employed this new implementation to revisit the calculation of static and frequency-dependent first hyperpolarizabilities of hydrogen halides (HX, X = F-Ts) and the Verdet constant of heavy noble gas atoms (Xe, Rn, Og) and of selected hydrogen halides (HF to HI), in order to investigate the differences and similarities of QR-CC and the more approximate QR-EOMCC. Furthermore, we have determined the relative importance of scalar relativistic effects and spin-orbit coupling to these properties, through a comparison of different Hamiltonians, and extended our calculations to superheavy element species (HTs for hyperpolarizabilities, Og for the Verdet constant). Our results show that as one moves toward the bottom of the periodic table, QR-EOMCC can yield rather different results (hyperpolarizabilities) or perform rather similarly (Verdet constant) to QR-CC. These results underscore the importance of further characterizing the performance of QR-EOMCC for heavy element systems.