Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy
30 Aug 2005-Physical Chemistry Chemical Physics (The Royal Society of Chemistry)-Vol. 7, Iss: 18, pp 3297-3305
TL;DR: A large set of more than 300 molecules representing all elements-except lanthanides-in their common oxidation states was used to assess the quality of the bases all across the periodic table, and recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
Abstract: Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 molecules representing (nearly) all elements-except lanthanides-in their common oxidation states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, density functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
Citations
More filters
TL;DR: The revised DFT-D method is proposed as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.
Abstract: The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.
32,589 citations
TL;DR: A new density functional of the generalized gradient approximation (GGA) type for general chemistry applications termed B97‐D is proposed, based on Becke's power‐series ansatz from 1997, and is explicitly parameterized by including damped atom‐pairwise dispersion corrections of the form C6 · R−6.
Abstract: A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.
23,058 citations
TL;DR: It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT‐D methods has only a minor impact on the quality of the results and BJ‐damping seems to provide a physically correct short‐range behavior of correlation/dispersion even with unmodified standard functionals.
Abstract: It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.
14,151 citations
TL;DR: A series of auxiliary basis sets to fit Coulomb potentials for the elements H to Rn (except lanthanides) is presented and computation times for the Coulomb part are reduced by a factor of ca.15 kJ mol(-1) per atom.
Abstract: A series of auxiliary basis sets to fit Coulomb potentials for the elements H to Rn (except lanthanides) is presented. For each element only one auxiliary basis set is needed to approximate Coulomb energies in conjunction with orbital basis sets of split valence, triple zeta valence and quadruple zeta valence quality with errors of typically below ca. 0.15 kJ mol−1 per atom; this was demonstrated in conjunction with the recently developed orbital basis sets of types def2-SV(P), def2-TZVP and def2-QZVPP for a large set of small molecules representing (nearly) each element in all of its common oxidation states. These auxiliary bases are slightly more than three times larger than orbital bases of split valence quality. Compared to non-approximated treatments, computation times for the Coulomb part are reduced by a factor of ca. 8 for def2-SV(P) orbital bases, ca. 25 for def2-TZVP and ca. 100 for def2-QZVPP orbital bases.
4,876 citations
TL;DR: The construction of transferable, hierarchical basis sets are demonstrated, allowing the calculation to range from qualitative tight-binding like accuracy to meV-level total energy convergence with the basis set, since all basis functions are strictly localized.
2,178 citations
References
More filters
TL;DR: In this article, a semi-empirical exchange correlation functional with local spin density, gradient, and exact exchange terms was proposed. But this functional performed significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
Abstract: Despite the remarkable thermochemical accuracy of Kohn–Sham density‐functional theories with gradient corrections for exchange‐correlation [see, for example, A. D. Becke, J. Chem. Phys. 96, 2155 (1992)], we believe that further improvements are unlikely unless exact‐exchange information is considered. Arguments to support this view are presented, and a semiempirical exchange‐correlation functional containing local‐spin‐density, gradient, and exact‐exchange terms is tested on 56 atomization energies, 42 ionization potentials, 8 proton affinities, and 10 total atomic energies of first‐ and second‐row systems. This functional performs significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
87,732 citations
TL;DR: In this paper, a detailed study of correlation effects in the oxygen atom was conducted, and it was shown that primitive basis sets of primitive Gaussian functions effectively and efficiently describe correlation effects.
Abstract: In the past, basis sets for use in correlated molecular calculations have largely been taken from single configuration calculations. Recently, Almlof, Taylor, and co‐workers have found that basis sets of natural orbitals derived from correlated atomic calculations (ANOs) provide an excellent description of molecular correlation effects. We report here a careful study of correlation effects in the oxygen atom, establishing that compact sets of primitive Gaussian functions effectively and efficiently describe correlation effects i f the exponents of the functions are optimized in atomic correlated calculations, although the primitive (s p) functions for describing correlation effects can be taken from atomic Hartree–Fock calculations i f the appropriate primitive set is used. Test calculations on oxygen‐containing molecules indicate that these primitive basis sets describe molecular correlation effects as well as the ANO sets of Almlof and Taylor. Guided by the calculations on oxygen, basis sets for use in correlated atomic and molecular calculations were developed for all of the first row atoms from boron through neon and for hydrogen. As in the oxygen atom calculations, it was found that the incremental energy lowerings due to the addition of correlating functions fall into distinct groups. This leads to the concept of c o r r e l a t i o n c o n s i s t e n t b a s i s s e t s, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation consistent sets are given for all of the atoms considered. The most accurate sets determined in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding ANO sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estimated that this set yields 94%–97% of the total (HF+1+2) correlation energy for the atoms neon through boron.
26,705 citations
TL;DR: Numerical results for atoms, positive ions, and surfaces are close to the exact correlation energies, with major improvements over the original LM approximation for the ions and surfaces.
Abstract: Langreth and Mehl (LM) and co-workers have developed a useful spin-density functional for the correlation energy of an electronic system. Here the LM functional is improved in two ways: (1) The natural separation between exchange and correlation is made, so that the density-gradient expansion of each is recovered in the slowly varying limit. (2) Uniform-gas and inhomogeneity effects beyond the randomphase approximation are built in. Numerical results for atoms, positive ions, and surfaces are close to the exact correlation energies, with major improvements over the original LM approximation for the ions and surfaces.
16,378 citations
TL;DR: In this article, various contracted Gaussian basis sets for atoms up to Kr are presented which have been determined by optimizing atomic self-consistent field ground state energies with respect to all basis set parameters, i.e., orbital exponents and contraction coefficients.
Abstract: Various contracted Gaussian basis sets for atoms up to Kr are presented which have been determined by optimizing atomic self‐consistent field ground state energies with respect to all basis set parameters, i.e., orbital exponents and contraction coefficients.
8,258 citations
TL;DR: In this article, a triple zeta valence (TZV) basis set is presented for Li to Kr. The TZV bases are characterized by typically including a single contraction to describe inner shells and three basis functions for valence shells.
Abstract: Contracted Gaussian basis sets of triple zeta valence (TZV) quality are presented for Li to Kr. The TZV bases are characterized by typically including a single contraction to describe inner shells and three basis functions for valence shells. All parameters—orbital exponents and contraction coefficients—have been determined by minimization of atomic self‐consistent field ground state energies. Advantages and necessary modifications of TZV basis sets are discussed for simple test calculations of molecular energies and nuclear magnetic resonance (NMR) chemical shieldings in treatments with and without inclusion of electron correlation.
7,859 citations