On the difference between the transition properties calculated with linear response- and equation of motion-CCSD approaches

In this work, we quantitatively investigate the difference between the linear response (LR) and the equation of motion (EOM) coupled cluster (CC) approaches in the calculation of transition properties, namely, dipole and oscillator strengths, for the most widely used truncated CC wave function, which includes single and double excitation operators. We compare systems of increasing size, where the size-extensivity may be important. Our results suggest that, for small molecules, the difference is small even with large basis sets. The difference increases with the size of the system, but it is numerically small until hundreds of electron pairs are correlated. Although these calculations may be possible in a few years, at present the EOM approach is more advantageous, albeit more approximate, because it is computationally less demanding.