Multiwavelets are emerging as an attractive alternative to traditional basis sets such as Gaussian-type orbitals and plane waves. One of their distinctive properties is the ability to reach the basis set limit (often a chimera for traditional approaches) reliably and consistently by fixing the desired precision ε. We present our multiwavelet implementation of the linear response formalism, applied to static magnetic properties, at the self-consistent field level of theory (both for Hartree–Fock and density functional theories). We demonstrate that the multiwavelets consistently improve the accuracy of the results when increasing the desired precision, yielding results that have four to five digits precision, thus providing a very useful benchmark which could otherwise only be estimated by extrapolation methods. Our results show that magnetizabilities obtained with the augmented quadruple-ζ basis (aug-cc-pCVQZ) are practically at the basis set limit, whereas absolute nuclear magnetic resonance shielding tensors are more challenging: even by making use of a standard extrapolation method, the accuracy is not substantially improved. In contrast, our results provide a benchmark that: (1) confirms the validity of the extrapolation ansatz; (2) can be used as a reference to achieve a property-specific extrapolation scheme, thus providing a means to obtain much better extrapolated results; (3) allows us to separate functional-specific errors from basis-set ones and thus to assess the level of cancellation between basis set and functional errors often exploited in density functional theory.

Benchmarking molecular properties with Gaussian-type orbital (GTO) basis sets can be challenging, because one has to assume that the computed property is at the complete basis set (CBS) limit, without a robust measure of the error. Multiwavelet (MW) bases can be systematically improved with a controllable error, which eliminates the need for such assumptions. In this work, we have used MWs within Kohn–Sham density functional theory to compute static polarizabilities for a set of 92 closed-shell and 32 open-shell species. The results are compared to recent benchmark calculations employing the GTO-type aug-pc4 basis set. We observe discrepancies between GTO and MW results for several species, with open-shell systems showing the largest deviations. Based on linear response calculations, we show that these discrepancies originate from artifacts caused by the field strength and that several polarizabilies from a previous study were contaminated by higher order responses (hyperpolarizabilities). Based on our MW benchmark results, we can affirm that aug-pc4 is able to provide results close to the CBS limit, as long as finite difference effects can be controlled. However, we suggest that a better approach is to use MWs, which are able to yield precise finite difference polarizabilities even with small field strengths.

Transition metal-catalyzed reactions invariably include steps where ligands associate or dissociate. In order to obtain reliable energies for such reactions, sufficiently large basis sets need to be employed. In this paper, we have used high-precision multiwavelet calculations to compute the metal–ligand association energies for 27 transition metal complexes with common ligands, such as H₂, CO, olefins, and solvent molecules. By comparing our multiwavelet results to a variety of frequently used Gaussian-type basis sets, we show that counterpoise corrections, which are widely employed to correct for basis set superposition errors, often lead to underbinding. Additionally, counterpoise corrections are difficult to employ when the association step also involves a chemical transformation. Multiwavelets, which can be conveniently applied to all types of reactions, provide a promising alternative for computing electronic interaction energies free from any basis set errors.