Minimally augmented Karlsruhe basis sets

We propose an extension of the basis sets proposed by Ahlrichs and coworkers at Karlsruhe (these basis sets are designated as the second-generation default or “def2” basis sets in the Turbomole program). The Karlsruhe basis sets are very appealing because they constitute balanced and economical basis sets of graded quality from partially polarized double zeta to heavily polarized quadruple zeta for all elements up to radon (Z = 86). The extension consists of adding a minimal set of diffuse functions to a subset of the elements. This yields basis sets labeled minimally augmented or with “ma” as a prefix. We find that diffuse functions are not quite as important for the def2 basis sets as they are for Pople basis sets, but they are still necessary for good results on barrier heights and electron affinities. We provide assessments and validations of this extension for a variety of data sets and representative cases. We recommend the new ma-TZVP basis set for general-purpose applications of density functional theory.     

Comments on relativistic basis sets

 Adding the tight and diffuse Gaussian-type functions (GTFs), Faegri's variationally determined double-zeta-quality basis sets for molecular relativistic calculations are examined. An example atom is Cm. When the tight s-type GTF is added the total energy increases, whereas when diffuse GTFs are added the total energy decreases. The reasons for these findings are clarified. It is also pointed out that not only the Faegri's sets but also other variationally determined basis sets would show similar behavior so far as the expansion terms are not sufficient. Received: 22 July 2002 / Accepted: 21 October 2002 / Published online: 31 January 2003 Correspondence to: H. Tatewaki e-mail: htatewaki@nsc.nagoya.cu.ac.jp     

Gaussian basis sets for calculation of spin densities in first-row atoms

Summary The suitability of Gaussian basis sets for ab initio calculation of Fermi contact spin densities is established by application to the prototype first-row atoms B-F having open shell p electrons. Small multiconfiguration self-consistent-field wave functions are used to describe relevant spin and orbital polarization effects. Basis sets are evaluated by comparing the results to highly precise numerical grid calculations previously carried out with the same wave function models. It is found that modest contracted Gaussian basis sets developed primarily for Hartree-Fock calculations can give semiquantitative results if augmented by diffuse functions and if further uncontracted in the outer core-inner valence region.     

Reactivity of haloethanes with hydroxyl radicals: Effects of basis set and correlation energy on reaction energetics

Ab initio calculations on fluoroethane reactions with the hydroxyl radical have been carried out at different levels of theory. The convergence of reaction barriers and reaction enthalpies has been systematically investigated with respect to the size and quality of the basis set and the treatment of correlation energy. The G2 and MP2 barrier heights and reaction enthalpies show the best agreement with the experimental data. The split valence basis sets of triple-zeta quality supplemented by diffuse and polarization functions are necessary to reproduce experimental values for barrier heights and reaction enthalpies at the MP2 level of theory. The full counterpoise correction was used to calculate the basis set superposition error for several standard basis sets, including polarization and diffuse functions. The smallest counterpoise corrections are associated with basis sets that contain polarization and diffuse functions, the diffuse functions being the most effective in reducing BSSE. However, in our case, the uncorrected barrier heights are in better agreement with experimental results than the counterpoise-corrected data. Thus, at the MP2 level of theory, which seems to be dictated for larger electronic systems of chemical interest, the optimal approach is to increase the basis set to the maximum size affordable and to use results without counterpoise corrections for the calculation of reaction barriers. A viable alternative is the use of G2 theory because its results for the barrier heights and reaction enthalpies are in excellent agreement with the experimental data. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1190–1199     

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.     

Uniform quality gaussian basis sets for molecular calculations. IV. Gradient and charged optimized basis sets for CH4

Comparison of the molecular Q-optimized and molecular gradient optimized carbon basis sets for CH 4 showed that molecular Q optimization is an excellent substitute to the more expensive molecular gradient optimization. The parameter Q of the Q optimization is related to the population (i.e., net charge) on the atom.     

Basis sets for geometry optimizations of second-row transition metal inorganic and organometallic complexes

Modest-sized basis sets for the second-row transition metal atoms are developed for use in geometry optimization calculations. Our method is patterned after previous work on basis sets for first-row transition metal atoms. The basis sets are constructed from the minimal basis sets of Huzinaga and are augmented with a set of diffuse p and d functions. The exponents of these diffuse functions are chosen to minimize both the difference between the calculated and experimental equilibrium geometries and the total molecular energies for several second-row transition metal inorganic and organon etallic complexes. Slightly smaller basis sets, based on the same Huzinaga minimal sets but augmented with a set of diffuse s and p functions rather than diffuse p and d functions, are also presented. The performance of these basis sets is tested on a wide variety of second-row transition metal inorganic and organometallic complexes and is compared to pseudopotential basis sets incorporating effective core potentials.     

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.     

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.     

Cubic-grid Gaussian basis sets for electron scattering calculations. II. Matrix elements

In this article, we report an efficient computational procedure for electron scattering matrix elements in the previously developed cubic-grid Gaussian basis sets. The Green function matrix elements derived for the cubic-grid basis set are simpler and easier to calculate than are those available in the literature for conventional Gaussian basis sets. Special features of the cubic-grid basis sets may also be exploited for a very efficient computation of Coulomb and exchange integrals. Inelastic scattering amplitudes for vibrational excitations may be efficiently calcualted in the harmonic approximation by numerical differention of the T-matrix elements. © 1995 John Wiley & Sons, Inc.     

Choosing GTO basis sets for periodic HF calculations

We investigate numerical linear dependencies of Gaussian-type orbital basis sets employed in the framework of the Hartree-Fock self-consistent field method for periodic structures, which so far have hampered the use of extended basis sets for non-ionic crystals. These linear dependencies occur when diffuse basis functions are included in a basis set in an uncontrolled manner. We use the condition number of the overlap matrix to lead us in the construction of extended basis sets for periodic structures which avoid numerical linear dependencies. Extended basis sets of high quality are optimized for a number of periodic structures (fcc He, α-Be, α-BN, and B1 NaF) with respect to the energy of the constituent atoms or ions. The results obtained with our basis sets, which do not require reoptimization in the crystal environment, compare favorably with those obtained with other extended basis sets reported in the literature. Received: 20 July 1998/Accepted: 21 August 1998 / Published online: 19 October 1998     

Development and testing of a compact basis set for use in effective core potential calculations on rhodium complexes

We present a set of effective core potential (ECP) basis sets for rhodium atoms which are of reasonable size for use in electronic structure calculations. In these ECP basis sets, the Los Alamos ECP is used to simulate the effect of the core electrons while an optimized set of Gaussian functions, which includes polarization and diffuse functions, is used to describe the valence electrons. These basis sets were optimized to reproduce the ionization energy and electron affinity of atomic rhodium. They were also tested by computing the electronic ground state geometry and harmonic frequencies of [Rh(CO)₂μ‐Cl]₂, Rh(CO)₂ClPy, and RhCO (neutral and its positive, and negative ions) as well as the enthalpy of the reaction of [Rh(CO)₂μ‐Cl]₂ with pyridine (Py) to give Rh(CO)₂ClPy, at different levels of theory. Good agreement with experimental values was obtained. Although the number of basis functions used in our ECP basis sets is smaller than those of other ECP basis sets of comparable quality, we show that the newly developed ECP basis sets provide the flexibility and precision required to reproduce a wide range of chemical and physical properties of rhodium compounds. Therefore, we recommend the use of these compact yet accurate ECP basis sets for electronic structure calculations on molecules involving rhodium atoms. © 2012 Wiley Periodicals, Inc.