Frequency-dependent nonlinear optical properties with explicitly correlated coupled-cluster response theory using the CCSD(R12) model

Response theory up to infinite order is combined with the explicitly correlated coupled-cluster singles and doubles model including linear-r(12) corrections, CCSD(R12). The additional terms introduced by the linear-r(12) contributions, not present in the conventional CCSD calculation, are derived and discussed with respect to the extra costs required for their evaluation. An implementation is presented up to the cubic response function for one-electron perturbations, i.e., up to frequency-dependent second hyperpolarizabilities. As first applications the authors computed the electronic polarizabilities and second hyperpolarizabilities of BH, N(2), and formaldehyde and show that the improvement in the one-electron basis set convergence known from the R12 method for ground state energies is retained for higher-order optical properties. Frequency-dependent results are presented for the second hyperpolarizability of N(2).