VAMPyR-A high-level Python library for mathematical operations in a multiwavelet representation.

Wavelets and multiwavelets have lately been adopted in quantum chemistry to overcome challenges presented by the two main families of basis sets: Gaussian atomic orbitals and plane waves. In addition to their numerical advantages (high precision, locality, fast algorithms for operator application, linear scaling with respect to system size, to mention a few), they provide a framework that narrows the gap between the theoretical formalism of the fundamental equations and the practical implementation in a working code. This realization led us to the development of the Python library called VAMPyR (Very Accurate Multiresolution Python Routines). VAMPyR encodes the binding to a C++ library for multiwavelet calculations (algebra and integral and differential operator application) and exposes the required functionality to write a simple Python code to solve, among others, the Hartree-Fock equations, the generalized Poisson equation, the Dirac equation, and the time-dependent Schrödinger equation up to any predefined precision. In this study, we will outline the main features of multiresolution analysis using multiwavelets and we will describe the design of the code. A few illustrative examples will show the code capabilities and its interoperability with other software platforms.

Full Breit Hamiltonian in the Multiwavelets Framework.

New techniques in core-electron spectroscopy are necessary to resolve the structures of oxides of f-elements and other strongly correlated materials that are present only as powders and not as single crystals. Thus, accurate quantum chemical methods must be developed to calculate core spectroscopic properties in such materials. In this contribution, we present an important development in this direction, extending our fully adaptive real-space multiwavelet basis framework to tackle the four-component Dirac-Coulomb-Breit Hamiltonian. We show that multiwavelets can reproduce one-dimensional grid-based approaches. They are however a fully three-dimensional approach which can later be extended to molecules and materials. Our multiwavelet implementation attained precise results irrespective of the chosen nuclear model, provided that the error threshold is tight enough and that the chosen polynomial basis is sufficiently large. Furthermore, our results confirmed that in two-electron species, the magnetic and Gauge contributions from s-orbitals are identical in magnitude and can account for the experimental evidence from K and L edges.

Scalar Relativistic Effects with Multiwavelets: Implementation and Benchmark.

The importance of relativistic effects in quantum chemistry is widely recognized, not only for heavier elements but throughout the periodic table. At the same time, relativistic effects are strongest in the nuclear region, where the description of electrons through a linear combination of atomic orbitals becomes more challenging. Furthermore, the choice of basis sets for heavier elements is limited compared with lighter elements where precise basis sets are available. Thanks to the framework of multiresolution analysis, multiwavelets provide an appealing alternative to overcoming this challenge: they lead to robust error control and adaptive algorithms that automatically refine the basis set description until the desired precision is reached. This allows one to achieve a proper description of the nuclear region. In this work, we extended the multiwavelet-based code MRChem to the scalar zero-order regular approximation framework. We validated our implementation by comparing the total energies for a small set of elements and molecules. To confirm the validity of our implementation, we compared both against a radial numerical code for atoms and the plane-wave-based code EXCITING.

Quantifying Intramolecular Basis Set Superposition Errors.

We show that medium-sized Gaussian basis sets lead to significant intramolecular basis set superposition errors at Hartree-Fock and density functional levels of theory, with artificial stabilization of compact over extended conformations for a 186 atom deca-peptide. Errors of ∼80 and ∼10 kJ/mol are observed, with polarized double zeta and polarized triple zeta quality basis sets, respectively. Two different procedures for taking the basis set superposition error into account are tested. While both reduce the error, it appears that polarized quadruple zeta basis sets are required to reduce the error below a few kJ/mol. Alternatively, the basis set superposition error can be eliminated using multiresolution methods based on Multiwavelets.

Cavity-Free Continuum Solvation: Implementation and Parametrization in a Multiwavelet Framework.

We present a multiwavelet-based implementation of a quantum/classical polarizable continuum model. The solvent model uses a diffuse solute-solvent boundary and a position-dependent permittivity, lifting the sharp-boundary assumption underlying many existing continuum solvation models. We are able to include both surface and volume polarization effects in the quantum/classical coupling, with guaranteed precision, due to the adaptive refinement strategies of our multiwavelet implementation. The model can account for complex solvent environments and does not need a posteriori corrections for volume polarization effects. We validate our results against a sharp-boundary continuum model and find a very good correlation of the polarization energies computed for the Minnesota solvation database.

MRChem Multiresolution Analysis Code for Molecular Electronic Structure Calculations: Performance and Scaling Properties.

MRChem is a code for molecular electronic structure calculations, based on a multiwavelet adaptive basis representation. We provide a description of our implementation strategy and several benchmark calculations. Systems comprising more than a thousand orbitals are investigated at the Hartree-Fock level of theory, with an emphasis on scaling properties. With our design, terms that formally scale quadratically with the system size in effect have a better scaling because of the implicit screening introduced by the inherent adaptivity of the method: all operations are performed to the requested precision, which serves the dual purpose of minimizing the computational cost and controlling the final error precisely. Comparisons with traditional Gaussian-type orbitals-based software show that MRChem can be competitive with respect to performance.

Multiwavelets applied to metal-ligand interactions: Energies free from basis set errors.

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 H2, 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.

Static Polarizabilities at the Basis Set Limit: A Benchmark of 124 Species.

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.

The Elephant in the Room of Density Functional Theory Calculations.

Using multiwavelets, we have obtained total energies and corresponding atomization energies for the GGA-PBE and hybrid-PBE0 density functionals for a test set of 211 molecules with an unprecedented and guaranteed µHartree accuracy. These quasi-exact references allow us to quantify the accuracy of standard all-electron basis sets that are believed to be highly accurate for molecules, such as Gaussian-type orbitals (GTOs), all-electron numeric atom-centered orbitals (NAOs), and full-potential augmented plane wave (APW) methods. We show that NAOs are able to achieve the so-called chemical accuracy (1 kcal/mol) for the typical basis set sizes used in applications, for both total and atomization energies. For GTOs, a triple-ζ quality basis has mean errors of ∼10 kcal/mol in total energies, while chemical accuracy is almost reached for a quintuple-ζ basis. Due to systematic error cancellations, atomization energy errors are reduced by almost an order of magnitude, placing chemical accuracy within reach also for medium to large GTO bases, albeit with significant outliers. In order to check the accuracy of the computed densities, we have also investigated the dipole moments, where in general only the largest NAO and GTO bases are able to yield errors below 0.01 D. The observed errors are similar across the different functionals considered here.

Magnetic properties with multiwavelets and DFT: the complete basis set limit achieved.

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.