Efficiency of numerical basis sets for predicting the binding energies of hydrogen bonded complexes: evidence of small basis set superposition error compared to Gaussian basis sets

Binding energies of selected hydrogen bonded complexes have been calculated within the framework of density functional theory (DFT) method to discuss the efficiency of numerical basis sets implemented in the DFT code DMol3 in comparison with Gaussian basis sets. The corrections of basis set superposition error (BSSE) are evaluated by means of counterpoise method. Two kinds of different numerical basis sets in size are examined; the size of the one is comparable to Gaussian double zeta plus polarization function basis set (DNP), and that of the other is comparable to triple zeta plus double polarization functions basis set (TNDP). We have confirmed that the magnitudes of BSSE in these numerical basis sets are comparative to or smaller than those in Gaussian basis sets whose sizes are much larger than the corresponding numerical basis sets; the BSSE corrections in DNP are less than those in the Gaussian 6-311+G(3df,2pd) basis set, and those in TNDP are comparable to those in the substantially large scale Gaussian basis set aug-cc-pVTZ. The differences in counterpoise corrected binding energies between calculated using DNP and calculated using aug-cc-pVTZ are less than 9 kJ/mol for all of the complexes studied in the present work. The present results have shown that the cost effectiveness in the numerical basis sets in DMol3 is superior to that in Gaussian basis sets in terms of accuracy per computational cost.     

Newly developed basis sets for density functional calculations

Optimized contracted Gaussian basis sets of double-zeta valence polarized (DZVP) quality for first-row transition metals are presented. The DZVP functions were optimized using the PWP86 generalized gradient approximation (GGA) functional and the B3LYP hybrid functional. For a careful analysis of the basis sets performance the transition metal atoms and cations excitation energies were calculated and compared with the experimental ones. The calculated values were also compared with those obtained using the previously available DZVP basis sets developed at the local-density functional level. Because the new basis sets work better than the previous ones, possible reasons of this behavior are analyzed. The newly developed basis sets also provide a good estimation of other atomic properties such as ionization energies.     

Bond functions,covalent potential curves,and the basis set superposition error

In the current practice of quantum chemistry, it is not clear whether corrections for basis set superposition errors should be applied to the calculation of potential energy curves, in order to improve agreement with experimental data. To examine this question, spectroscopic parameters derived from theoretical potential curves are reported for the homonuclear diatomics C₂, N₂, O₂, and F₂, using a configuration interaction method. Three different basis sets were used, including double zeta plus polarization, triple zeta plus double polarization, and double zeta polarization augmented by bond functions. The bond function basis sets, which were optimized in the preceding paper to obtain accurate dissociation energies, also gave the most accurate parameters. The potential curves were then corrected for basis set superposition error using the counterpoise correction, and the spectroscopic parameters were computed again. The BSSE-corrected curves showed worse agreement with experiment for all properties than the original (uncorrected) curves. The reasons for this finding are discussed. In addition to the numerical results, some problems in the application of the BSSE correction to basis sets containing bond functions are shown. In particular, there is an overcounting of the lowering due to the bond functions, regardless of which type of correction is applied. Also, genuine BSSE affects cannot be separated from energy-lowering effects due to basis set incompleteness, and we postulate that it is the latter which is strongly dominant in the calculation of covalent potential curves. Based on these arguments, two conclusions follow: (1) application of BSSE corrections to potential curves should not be routinely applied in situations where the bonding is strong, and (2) appropriate use of bond functions can lead to systematic improvement in the quality of potential curves.     

Quality of contracted Gaussian-type function basis sets

The valence quality of contracted (C) Gaussian-type function (GTF) basis sets in molecular calculations is discussed for the first- through fourth-row atoms. The split-valence basis sets derived from minimal-type CGTF sets are compared with those derived from primitive (P) GTF sets. Using F, Cl, Br, and I atoms and their homonuclear diatomics as test species, we find that the split-valence CGTF sets have almost the same quality as PGTF sets with larger s and p expansion terms: for example, the (53/5), (533/53), (5333/533/5), and (53 333/5333/53) CGTF sets correspond approximately to the [9/5], [15/9], [19/15/5], and [22/18/7] PGTF sets for the first- to fourth-row atoms, respectively, where the slash separates the s, p, and d symmetries. For the main group atoms of the four rows, we recommend using the above-mentioned CGTFs or larger.     

Property-optimized gaussian basis sets for molecular response calculations

With recent advances in electronic structure methods, first-principles calculations of electronic response properties, such as linear and nonlinear polarizabilities, have become possible for molecules with more than 100 atoms. Basis set incompleteness is typically the main source of error in such calculations since traditional diffuse augmented basis sets are too costly to use or suffer from near linear dependence. To address this problem, we construct the first comprehensive set of property-optimized augmented basis sets for elements H-Rn except lanthanides. The new basis sets build on the Karlsruhe segmented contracted basis sets of split-valence to quadruple-zeta valence quality and add a small number of moderately diffuse basis functions. The exponents are determined variationally by maximization of atomic Hartree-Fock polarizabilities using analytical derivative methods. The performance of the resulting basis sets is assessed using a set of 313 molecular static Hartree-Fock polarizabilities. The mean absolute basis set errors are 3.6%, 1.1%, and 0.3% for property-optimized basis sets of split-valence, triple-zeta, and quadruple-zeta valence quality, respectively. Density functional and second-order M?ller-Plesset polarizabilities show similar basis set convergence. We demonstrate the efficiency of our basis sets by computing static polarizabilities of icosahedral fullerenes up to C(720) using hybrid density functional theory.     

Minimal basis sets in calculations of intermolecular interaction energies

Several minimal (7, 3/3) Gaussian basis sets have been used to calculate the energies and some other properties of CH₄ and H₂O. Improved basis sets developed for these molecules have been extended to NH₃ and HF and employed to H₂CO and CH₃OH. Interaction energies between XH_n molecules have been calculated using the old and the new minimal basis sets. The results obtained with the new basis sets are comparable in accuracy to those calculated with significantly more extended basis sets involving polarization functions. Binding energies calculated using the counterpoise method are not much different for the new and the old minimal basis sets, and are likely to be more accurate than the results of much more extended calculations.     

Performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene

This work compares the performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene. Silylene, its cation and anion have been studied in ¹ A ₁, ² A ₁ and ² B ₁ states, respectively, in the gas phase and C_2v symmetry. The methods considered are second-order Møller-Plesset perturbation theory (MP2), the density functional theory (DFT), Gaussian-2 (G2) and complete basis set methods (CBS-4M and CBS-Q). The basis sets used are 6-31G(d,p), 6-311G(d,p), 6-31++G(d,p) and 6-311++G(d,p). The functional used for the DFT method is B3LYP. Silylene and its cation and anion have been optimised using the MP2 and DFT methods and the named basis sets. Single-point energy calculations (G2, CBS-4M and CBS-Q) were performed using MP2/6-311++G(d,p) structures and these energies have been used to calculate atomisation energy, ionisation energy and adiabatic electron affinity. Frequency calculations were also done and the raw vibrational frequencies were assigned. It is interesting to note the close similarity between the predicted parameters and some of the available literature values. The results obtained are consistent and converge with different basis sets with improved size and quality. However, the parameters obtained are very much method dependent.     

Intramolecular basis set superposition errors

The effects of intramolecular basis set superposition errors are less well documented than the corresponding intermolecular effects. The intramolecular basis set superposition errors are examined, using the approach of Jensen, for several basis sets developed by Pople and his co‐workers, which are widely used in studies of larger molecules. Prototype calculations are reported for the ground state of the water molecule using both the matrix Hartree–Fock method and the many‐body perturbation expansion for the correlation energy taken through second order. A similar investigation is carried out for some of the correlation consistent basis sets published by Dunning and his collaborators. Specifically, the following aspects are investigated: (i) the magnitude of the intramolecular basis set superposition error, (ii) the nonadditivity of intramolecular counterpoise corrections when applied in a pairwise fashion, and (iii) the use of multiple “ghost” centers. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 282–292, 2001     

Systematically convergent basis sets for explicitly correlated wavefunctions: the atoms H, He, B-Ne, and Al-Ar

Correlation consistent basis sets have been optimized for use with explicitly correlated F12 methods. The new sets, denoted cc-pVnZ-F12 (n=D,T,Q), are similar in size and construction to the standard aug-cc-pVnZ and aug-cc-pV(n+d)Z basis sets, but the new sets are shown in the present work to yield much improved convergence toward the complete basis set limit in MP2-F12/3C calculations on several small molecules involving elements of both the first and second row. For molecules containing only first row atoms, the smallest cc-pVDZ-F12 basis set consistently recovers nearly 99% of the MP2 valence correlation energy when combined with the MP2-F12/3C method. The convergence with basis set for molecules containing second row atoms is slower, but the new DZ basis set still recovers 97%-99% of the frozen core MP2 correlation energy. The accuracy of the new basis sets for relative energetics is demonstrated in benchmark calculations on a set of 15 chemical reactions.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

Approximations to complete basis set-extrapolated, highly correlated non-covalent interaction energies

The first-principles calculation of non-covalent (particularly dispersion) interactions between molecules is a considerable challenge. In this work we studied the binding energies for ten small non-covalently bonded dimers with several combinations of correlation methods (MP2, coupled-cluster single double, coupled-cluster single double (triple) (CCSD(T))), correlation-consistent basis sets (aug-cc-pVXZ, X = D, T, Q), two-point complete basis set energy extrapolations, and counterpoise corrections. For this work, complete basis set results were estimated from averaged counterpoise and non-counterpoise-corrected CCSD(T) binding energies obtained from extrapolations with aug-cc-pVQZ and aug-cc-pVTZ basis sets. It is demonstrated that, in almost all cases, binding energies converge more rapidly to the basis set limit by averaging the counterpoise and non-counterpoise corrected values than by using either counterpoise or non-counterpoise methods alone. Examination of the effect of basis set size and electron correlation shows that the triples contribution to the CCSD(T) binding energies is fairly constant with the basis set size, with a slight underestimation with CCSD(T)∕aug-cc-pVDZ compared to the value at the (estimated) complete basis set limit, and that contributions to the binding energies obtained by MP2 generally overestimate the analogous CCSD(T) contributions. Taking these factors together, we conclude that the binding energies for non-covalently bonded systems can be accurately determined using a composite method that combines CCSD(T)∕aug-cc-pVDZ with energy corrections obtained using basis set extrapolated MP2 (utilizing aug-cc-pVQZ and aug-cc-pVTZ basis sets), if all of the components are obtained by averaging the counterpoise and non-counterpoise energies. With such an approach, binding energies for the set of ten dimers are predicted with a mean absolute deviation of 0.02 kcal/mol, a maximum absolute deviation of 0.05 kcal/mol, and a mean percent absolute deviation of only 1.7%, relative to the (estimated) complete basis set CCSD(T) results. Use of this composite approach to an additional set of eight dimers gave binding energies to within 1% of previously published high-level data. It is also shown that binding within parallel and parallel-crossed conformations of naphthalene dimer is predicted by the composite approach to be 9% greater than that previously reported in the literature. The ability of some recently developed dispersion-corrected density-functional theory methods to predict the binding energies of the set of ten small dimers was also examined.     

Sensitivity of the properties of ruthenium "blue dimer" to method, basis set, and continuum model

The ruthenium "blue dimer" [(bpy)(2)Ru(III)OH(2)](2)O(4+) is best known as the first well-defined molecular catalyst for water oxidation. It has been subject to numerous computational studies primarily employing density functional theory. However, those studies have been limited in the functionals, basis sets, and continuum models employed. The controversy in the calculated electronic structure and the reaction energetics of this catalyst highlights the necessity of benchmark calculations that explore the role of density functionals, basis sets, and continuum models upon the essential features of blue-dimer reactivity. In this paper, we report Kohn-Sham complete basis set (KS-CBS) limit extrapolations of the electronic structure of "blue dimer" using GGA (BPW91 and BP86), hybrid-GGA (B3LYP), and meta-GGA (M06-L) density functionals. The dependence of solvation free energy corrections on the different cavity types (UFF, UA0, UAHF, UAKS, Bondi, and Pauling) within polarizable and conductor-like polarizable continuum model has also been investigated. The most common basis sets of double-zeta quality are shown to yield results close to the KS-CBS limit; however, large variations are observed in the reaction energetics as a function of density functional and continuum cavity model employed.     

Auxiliary basis sets for density-fitted correlated wavefunction calculations: weighted core-valence and ECP basis sets for post-d elements

We report optimised auxiliary basis sets for the resolution-of-the-identity (or density-fitting) approximation of two-electron integrals in second-order M?ller-Plesset perturbation theory (MP2) and similar electronic structure calculations with correlation-consistent basis sets for the post-d elements Ga-Kr, In-Xe, and Tl-Rn. The auxiliary basis sets are optimised such that the density-fitting error is negligible compared to the one-electron basis set error. To check to which extent this criterion is fulfilled we estimated for a test set of 80 molecules the basis set limit of the correlation energy at the MP2 level and evaluated the remaining density-fitting and the one-electron basis set errors. The resulting auxiliary basis sets are only 2-6 times larger than the corresponding one-electron basis sets and lead in MP2 calculations to speed-ups of the integral evaluation by one to three orders of magnitude. The density-fitting errors in the correlation energy are at least hundred times smaller than the one-electron basis set error, i.e. in the order of only 1-100 μH per atom.     

Gaussian basis sets for the fifth row elements,Mo-Cd,and the sixth row elements W-RN

Gaussian basis sets, consisting of 17 s-type, 12 p-type, and 8 d-type functions, for the fifth row elements, Mo? Cd, and 19 s-type, 14 p-type (16 p-type), 10 d-type and 5 f-type functions for the sixth row elements, W? Rn, are presented. The basis sets are of double zeta quality, and are optimized to .002 a.u. in the energy. The energies are compared with D.Z. STO basis sets.     

Optimized Slater-type basis sets for the elements 1-118

Seven different types of Slater type basis sets for the elements H (Z = 1) up to E118 (Z = 118), ranging from a double zeta valence quality up to a quadruple zeta valence quality, are tested in their performance in neutral atomic and diatomic oxide calculations. The exponents of the Slater type functions are optimized for the use in (scalar relativistic) zeroth-order regular approximated (ZORA) equations. Atomic tests reveal that, on average, the absolute basis set error of 0.03 kcal/mol in the density functional calculation of the valence spinor energies of the neutral atoms with the largest all electron basis set of quadruple zeta quality is lower than the average absolute difference of 0.16 kcal/mol in these valence spinor energies if one compares the results of ZORA equation with those of the fully relativistic Dirac equation. This average absolute basis set error increases to about 1 kcal/mol for the all electron basis sets of triple zeta valence quality, and to approximately 4 kcal/mol for the all electron basis sets of double zeta quality. The molecular tests reveal that, on average, the calculated atomization energies of 118 neutral diatomic oxides MO, where the nuclear charge Z of M ranges from Z = 1-118, with the all electron basis sets of triple zeta quality with two polarization functions added are within 1-2 kcal/mol of the benchmark results with the much larger all electron basis sets, which are of quadruple zeta valence quality with four polarization functions added. The accuracy is reduced to about 4-5 kcal/mol if only one polarization function is used in the triple zeta basis sets, and further reduced to approximately 20 kcal/mol if the all electron basis sets of double zeta quality are used. The inclusion of g-type STOs to the large benchmark basis sets had an effect of less than 1 kcal/mol in the calculation of the atomization energies of the group 2 and group 14 diatomic oxides. The basis sets that are optimized for calculations using the frozen core approximation (frozen core basis sets) have a restricted basis set in the core region compared to the all electron basis sets. On average, the use of these frozen core basis sets give atomic basis set errors that are approximately twice as large as the corresponding all electron basis set errors and molecular atomization energies that are close to the corresponding all electron results. Only if spin-orbit coupling is included in the frozen core calculations larger errors are found, especially for the heavier elements, due to the additional approximation that is made that the basis functions are orthogonalized on scalar relativistic core orbitals.     

Accurate extrapolation of electron correlation energies from small basis sets

A new two-point scheme is proposed for the extrapolation of electron correlation energies obtained with small basis sets. Using the series of correlation-consistent polarized valence basis sets, cc-pVXZ, the basis set truncation error is expressed as deltaE(X) proportional, variant(X + xi(i))(-gamma). The angular momentum offset xi(i) captures differences in effective rates of convergence previously observed for first-row molecules. It is based on simple electron counts and tends to values close to 0 for hydrogen-rich compounds and values closer to 1 for pure first-row compounds containing several electronegative atoms. The formula is motivated theoretically by the structure of correlation-consistent basis sets which include basis functions up to angular momentum L = X-1 for hydrogen and helium and up to L = X for first-row atoms. It contains three parameters which are calibrated against a large set of 105 reference molecules (H, C, N, O, F) for extrapolations of MP2 and CCSD valence-shell correlation energies from double- and triple-zeta (DT) and triple- and quadruple-zeta (TQ) basis sets. The new model is shown to be three to five times more accurate than previous two-point schemes using a single parameter, and (TQ) extrapolations are found to reproduce a small set of available R12 reference data better than even (56) extrapolations using the conventional asymptotic limit formula deltaE(X) proportional, variantX(-3). Applications to a small selection of boron compounds and to neon show very satisfactory results as well. Limitations of the model are discussed.