DOI: 10.18710/klqvok
title: Replication Data for: Static polarizabilities at the basis set limit: A benchmark of 124 species

description: &lt;h3&gt;Introduction&lt;/h3&gt;
This Dataverse entry contains replication data for our journal article “Static polarizabilities at the basis set limit:
A benchmark of 124 species” published in Journal of Chemical Theory and Computation. It contains highly precise static polarizabilities computed in multiwavelet basis in combination with density functional theory (DFT, PBE functional). In addition, the d/Preliminaryata set contains analysis tools (Jupyter Notebooks with Python3 code) for generating the figures in the journal article.

&lt;h3&gt;How to use&lt;/h3&gt;
Because our multiwavelet data is guaranteed to be at the complete basis set limit (to within the specified limit), it is suitable as a benchmark reference in studies of static polarizabilities where the basis set convergence is important. With multiwavelets we don't have to assume that the computed property is at the basis set limit, as is the case with Gaussian type orbital (GTO) basis sets, and it is therefore possible to confirm whether the property of interest computed basis is sufficiently converged with respect to the complete basis set limit. Our benchmark reference can also be used in the development of new methodology that requires accurate training data.

&lt;h3&gt;Running the Jupyter Notebooks&lt;/h3&gt;
The Anaconda Python distribution is usually recommended for obtaining Jupyter Notebook. It can be downloaded from here: https://www.anaconda.com/distribution/
&lt;br&gt;
&lt;br&gt;
The simplest way to run the notebooks is to download all files in this DataverseNO dataset. That will preserve the directory structure, which is absolutely necessary to avoid errors. Then start your Jupyter Notebook session, navigate to the data set directory, and open the desired notebook.

&lt;h3&gt;Journal article&lt;/h3&gt;
&lt;a href="https://doi.org/10.1021/acs.jctc.0c00128"&gt;Brakestad et al. "Static polarizabilities at the basis set limit: A benchmark of 124 species". J. Chem. Theory Comput. (2020)&lt;/a&gt;&lt;br&gt;&lt;br&gt;

&lt;h3&gt;Abstract from journal article&lt;/h3&gt;
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 artefacts 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.
description: ORCA, 4.1.2
description: MRChem, 0
description: Python, 3.6