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.     

Simplified Box Orbitals (SBO) for H To Ar atoms: Exact expressions,SBO‐3G approximations,and relations with the ZDO approximation

Simplified Box Orbitals (SBO) are a kind of spatially restricted basis functions. SBOs have a similar use and value to Slater functions but, because they fulfill a version of the zero‐differential overlap approximation, they allow for a drastic reduction in the number of two‐electron integrals to be calculated when dealing with huge systems, and they seem to be specially adapted to study confined systems such as molecules in solution. In a previous study, the mathematical shape of SBOs was discussed and the necessary parameters were obtained by means of the variational method. In the present study, the parameters of each SBO were obtained by applying the condition that it is as similar as possible to the STO that would be used in a basis set without spatial restrictions. We have developed a method to achieve this likeness and deduced simple formulas to describe all the SBOs of any atom. We also present the SBO‐3G expansions of the SBOs obtained, making it possible to use these SBOs with standard quantum chemistry calculation software. Simple formulas were also deduced to directly write the SBOs and SBO‐3G corresponding to the atoms with a Z value of between 1 and 18. Finally, as a first example of the usefulness of this kind of functions, an optimized SBO‐3G basis set is proposed for atoms from H to Cl in molecules. © 2016 Wiley Periodicals, Inc.     

Toward accurate thermochemical models for transition metals: G3Large basis sets for atoms Sc-Zn

An augmented valence triple-zeta basis set, referred to as G3Large, is reported for the first-row transition metal elements Sc through Zn. The basis set is constructed in a manner similar to the G3Large basis set developed previously for other elements (H-Ar, K, Ca, Ga-Kr) and used as a key component in Gaussian-3 theory. It is based on a contraction of a set of 15s13p5d Gaussian primitives to 8s7p3d, and also includes sets of f and g polarization functions, diffuse spd functions, and core df polarization functions. The basis set is evaluated with triples-augmented coupled cluster [CCSD(T)] and Brueckner orbital [BD(T)] methods for a small test set involving energies of atoms, atomic ions, and diatomic hydrides. It performs well for the low-lying s-->d excitation energies of atoms, atomic ionization energies, and the dissociation energies of the diatomic hydrides. The Brueckner orbital-based BD(T) method performs substantially better than Hartree-Fock-based CCSD(T) for molecules such as NiH, where the starting unrestricted Hartree-Fock wavefunction suffers from a high degree of spin contamination. Comparison with available data for geometries of transition metal hydrides also shows good agreement. A smaller basis set without core polarization functions, G3MP2Large, is also defined.     

Cubic-grid Gaussian basis sets for electron scattering calculations. I. Definition and construction

We propose a new type of Gaussian basis sets for use in calculations of electron scattering by molecules. Instead of locating the basis-set functions on the atomic centers of the target molecule, we place primitive s-type Gaussians at the positions of a cubic lattice with a regular grid. The grid and the Gaussian exponent are fixed so as to give the best representation of the plane-wave function. Plane-wave functions and Green functions obtained by means of the cubic-grid basis set are tested graphically against exact functions and functions expressed by means of a conventional Gaussian basis set. © 1995 John Wiley & Sons, Inc.     

Correspondence between GTO and STO molecular basis sets

We present a method for the characterization of the distance between two spaces: one generated by a Gaussian basis set, and another by a Slater basis set. The method is an extension of one previously developed for atoms that has been modified to cover molecular problems. The current version enables us to obtain Slater basis sets capable of reproducing the results (multielectronic wave functions and orbitals) obtained with Gaussian basis sets. The interest of this result arises from the fact that we will be able to profit from the effort invested in the optimization of high‐quality Gaussian basis sets. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1655–1665, 2001     

Linear dependence and energy conservation in Gaussian wavepacket basis sets

We propose a method for dealing with the problem of linear dependence in quantum dynamics simulations employing over-complete Gaussian wavepacket (GWP) basis sets. In particular, by periodically projecting out redundant basis functions using the matching pursuit algorithm whilst simultaneously introducing GWPs which avoid linear dependence with the current basis set, we find that numerical conditioning of the equations-of-motion can be readily controlled. In applications to particle tunnelling in one- and two-dimensional potentials, this method allows us to reproduce the exact quantum-mechanical results with fewer GWP basis functions than similar calculations with non-adaptive basis sets, a result which we trace back to the improved energy conservation of our adaptive approach.     

Calculation of isotropic Compton profiles with Gaussian basis sets

In this paper we present an adaptive algorithm for calculating the isotropic Compton profile (ICP) for any type of Gaussian basis set. The ICP is a measure of the momentum density of electrons and it can be obtained from inelastic X-ray scattering experiments employing synchrotron radiation. We have performed calculations of the ICP for water and helium monomers and dimers using density-functional theory, Hartree-Fock and post-Hartree-Fock methods, with Dunning-type ((d-)aug-)cc-p(C)VXZ basis sets. We have examined the convergence of the Compton profile as a function of the basis set and the level of theory used for the formation of the density matrix. We demonstrate that diffuse basis functions are of utmost importance to the calculation of Compton profiles. Basis sets of at least triple-ζ quality appended by diffuse functions should be used in Compton profile calculations in order to obtain sufficient convergence with regard to the current, experimentally feasible accuracy for systems consisting of light elements.     

Molecule‐adapted basis sets optimized with a quantum Monte Carlo method

The Monte Carlo simulated annealing method is adapted to optimize correlated Gaussian‐type functions in nonrelativistic molecular environments. Starting from an atom‐centered atomic Gaussian basis set, the uncontracted functions are reoptimized in the molecular environments corresponding to the H₂O, CN^?, N₂, CO, BF, NO⁺, CO₂, and CS systems. These new molecular adapted basis sets are used to calculate total energies, harmonic vibrational frequencies, and equilibrium geometries at a correlated level of theory. The present methodology is a simple and effective way to improve molecular correlated wave functions, without the need to enlarge the molecular basis set. Additionally, this methodology can be used to generate hierarchical sequences of molecular basis sets with increasing size, which are relevant to establish complete basis set limits. © 2014 Wiley Periodicals, Inc.     

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.     

Possible criterion for the balance of basis sets in quantum-chemical calculations

We propose to characterize the width of a basis set (BS) by the number of basis functions falling within one electron of the considered atomic or molecular systems. It is established that for atoms, this value (the electron function endowment, or EFE) undergoes drastic changes as the atomic number of a periodic system element increases. It is shown that widening the BS through the addition of valence and polarizational functions increases the imbalance of the basis sets of various atoms in terms of EFE. A scheme of construction is proposed and an example of a BS balanced according to the EFE value is given. The properties of LiH and HF molecules are calculated by the density functional UB3LYP method with the use of standard and molecular-optimized (relaxed) BSes with segmented and general contraction of the Gaussian functions. It is established that there are uniform dependences for the error of calculating the properties of both molecules from the EFE. We conclude that the accuracy of calculating the equilibrium distance, ionization energy, electron affinity, atomization energy, dipole moment, and frequency of normal vibration increases steadily as the EFE value of a molecule rises.     

Evaluation of Gaussian basis sets in ab initio

Ab initio wave functions for NH₃ have been calculated using the STONG program, where N is 2 to 6, and with Gaussian 70, employing the 4–31 and 6–31 basis sets. Electric field gradients at the N atom were computed from the gaussian basis orbitals, and, with STONG, also from the equivalent STO's. Results from the two split valence shell Gaussian 70 calculations are in good agreement with values from extended SCF and SCFCI calculations. Because of adequate agreement with more extended calculations, the use of the 4–31 basis set is chosen to calculate field gradients in molecules consisting of five or more second row atoms. It should be adequate for the prediction of the region for searching for new ¹⁴N, ¹⁷O and ²H NQR spectra. An improved value of the electric quadrupole moment for the ¹⁴N nucleus is proposed.     

Optimization of Gaussian basis sets for Dirac-Hartree-Fock calculations

We investigate the optimization of Gaussian basis sets for relativistic calculations within the framework of the restricted Dirac-Hartree-Fock (DHF) method for atoms. We compare results for Rn of nonrelativistic and relativistic basis set optimizations with a finite nuclear-size. Optimization of separate sets for each spin-orbit component shows that the basis set demands for the lower j component are greater than for the higher j component. In particular, the p _1/2 set requires almost as many functions as the s _1/2 set. This implies that for the development of basis sets for heavy atoms, the symmetry type for which a given number of functions is selected should be based on j, not on l, as has been the case in most molecular calculations performed to date.     

Uniform quality gaussian basis sets for molecular calculations. III. Charge optimized basis sets

Gradient optimized constrained (2s ≠ 2p) and unconstrained (2s ≠ 2p) Gaussian 3G basis sets are reported for the first-row atoms and ions X^O, for Q = ?2 to +4. Analytic equations have been fitted to the logarithm of the exponents as a function of the nuclear charge Z and formal charge Q. Consequently only two parameters Z and Q have to be specified in order to completely define a basis set.     

A study on universal Gaussian basis sets for first-row atoms

Analysis of various optimum and non-optimum Gaussian basis sets for firstrow elements have indicated that with a minimum increase of the basis set size and without loss of accuracy of the calculated total energy, a single universal Gaussian basis set may replace individually optimized Gaussian basis sets for a series of atoms. Such a universal Gaussian basis set may substantially reduce the computational work required for the calculation of molecular integrals in ab initio MO calculations.     

Locally dense basis sets for chemical shift calculations

Calculations of chemical shifts have been carried out using “locally dense” basis sets for the resonant atom of interest, and smaller, attenuated sets on other atoms in the molecule. For carbon, calculations involving a 6-311G(d) triply split valence set with polarization on the resonant atom and 3-21G atomic bases on other heavy atoms result in good agreement with experiment, and are virtually identical to those found employing the larger basis on all atoms. For species such as nitrogen, oxygen, and fluorine where standard balanced basis sets do not agree well with experiment, use of attenuated sets fail as well. The use of locally dense basis sets permits calculations previously impractical, and the successful application to carbon suggests that the chemical shift is most dependent on the local basis set, and less so on whether or not a balanced or unbalanced calculation is being carried out.     

Combined bond-polarization basis sets for accurate determination of dissociation energies. II. Basis set superposition error as a function of the parent basis set

The effect of the parent basis set on the basis set superposition error caused by bond functions is investigated systematically. An important difference between BSSE at the SCF and correlated levels is pointed out. Three new basis sets are defined, denoted 6-311 + G(d,p)B, 6-311 + G(2d,p)B, and 6-311 + G(2df,p)B. BSSE for the first-row hydrides seems to increase uniformly with increasing atomic number of the central atom. Expansion of the valence part of the basis set from 6-31G to 6-311G, as well as adding f functions, has a significant effect on the BSSE. Additional BSSEs incurred by bond functions are less than or equal to 1 kcal/mol for the 6-311 + G(2df,p)B basis set. For the dissociation energies of the first-row hydride species, agreement with experiment within only a few kcal/mol can be obtained even without resorting to isogyric reaction cycles. For high-quality calculations, adding bond functions seems to have definite advantages over expanding the polarization space beyond the [2d1f] level.     

Uniform quality gaussian basis sets for molecular calculations,I. C1 hydrocarbons

The optimality of MO basis sets of Gaussian functions, when constructed from AO basis sets optimized for the neutral atom or for atom ions, is investigated. A formal charge parameter Q is defined and used to adjust the AO basis sets to the molecular environment, by virtue of a simple quadratic expression. Calculations on a series of C₁ hydrocarbons (CH₂, CH₃, CH₃⁺, CH₃^?, CH₄) using 3G basis sets indicate considerable variations in the optimum Q value with the molecular species. The proposed method offers a simple alternative technique to a full molecular basis set optimization.