Dynamic polarizabilities for open- and closed-shell molecules were obtained by using coupled-cluster (CC) linear response theory with full treatment of singles, doubles, and triples (CCSDT-LR) with large basis sets utilizing the NWChem software suite. By using four approximate CC methods in conjunction with augmented cc-pVNZ basis sets, we are able to evaluate the convergence in both many-electron and one-electron spaces. For systems with primarily dynamic correlation, the results for CC3 and CCSDT are almost indistinguishable. For systems with significant static correlation, the CC3 tends to overestimate the triples contribution, while the PS(T) approximation [J. Chem. Phys. 127, 164105 (2007)] produces mixed results that are heavily dependent on the accuracies provided by noniterative approaches used to correct the equation-of-motion CCSD excitation energies. Our results for open-shell systems show that the choice of reference (restricted open-shell Hartree–Fock versus unrestricted Hartree–Fock) can have a significant impact on the accuracy of polarizabilities. A simple extrapolation based on pentuple-zeta CCSD calculations and triple-zeta CCSDT calculations reproduces experimental results with good precision in most cases.
The PS(T) results reported in this paper differ slightly with respect to the first paper (Ref. 20) for dynamic frequencies due to a bug in the code.