Platinum, gold, and silver standards of intermolecular interaction energy calculations

High-accuracy noncovalent interaction energies are indispensable as data points for potential energy surfaces and as benchmark values for improving and testing more approximate approaches. The preferred algorithm (the gold standard) for computing these energies has been the coupled-cluster method with singles, doubles, and perturbative triples [CCSD(T)] converged to the complete basis set (CBS) limit. However, gold-standard calculations are expensive as correlated interaction energies converge slowly with the basis set size, and establishing the CBS limit to better than 0.05 kcal/mol typically requires a CCSD(T) calculation in a basis set of at least triple-zeta quality. If an even higher accuracy is required (for example, for the assignment of complicated high-resolution spectra), establishing a superior platinum standard requires both a precisely converged CCSD(T)/CBS limit and the corrections for the core correlation, relativistic effects, and higher-order coupled-cluster terms at least through the perturbative quadruple excitations. On the other hand, if a triple-zeta CCSD(T) calculation is not feasible but a double-zeta one is, it is worthwhile to look for a silver standard that provides the most accurate and consistent approximation to the gold standard at a reduced computational cost. We review the recent developments aimed at (i) increasing the breadth and diversity of the available collection of gold-standard benchmark interaction energies, (ii) evaluating the best computational strategies for platinum-standard calculations and producing beyond-CCSD(T) potential energy surfaces for spectroscopic and scattering applications of the highest precision, and (iii) improving the accuracy of the silver-standard, double-zeta-level CCSD(T)/CBS estimates through the use of explicit correlation and midbond basis functions. We also outline the remaining challenges in the accurate ab initio calculations of noncovalent interaction energies.High-accuracy noncovalent interaction energies are indispensable as data points for potential energy surfaces and as benchmark values for improving and testing more approximate approaches. The preferred algorithm (the gold standard) for computing these energies has been the coupled-cluster method with singles, doubles, and perturbative triples [CCSD(T)] converged to the complete basis set (CBS) limit. However, gold-standard calculations are expensive as correlated interaction energies converge slowly with the basis set size, and establishing the CBS limit to better than 0.05 kcal/mol typically requires a CCSD(T) calculation in a basis set of at least triple-zeta quality. If an even higher accuracy is required (for example, for the assignment of complicated high-resolution spectra), establishing a superior platinum standard requires both a precisely converged CCSD(T)/CBS limit and the corrections for the core correlation, relativistic effects, and higher-order coupled-cluster terms at least through the per...