Basis set quality versus size II. Approximate GTO wave functions for second row transition metal atoms

Basis sets ranging in size from (16, 10, 7) to (20, 14, 11) have been derived for the atoms Y–Cd. Separate sets represent the energy optimized wave functions for each of the s²dⁿ, s¹dⁿ⁺¹, and s⁰dⁿ⁺² configurations. The energies from the largest sets are within 3 mhartrees of the values obtained in numerical Hartree–Fock calculations. Reasonable Hartree–Fock s²dⁿ– s¹dⁿ⁺¹ and s²dⁿ– s⁰dⁿ⁺² excitation energies may be obtained either using the largest basis sets, or using d-orbitals optimized for the s⁰dⁿ⁺² configurations. The basis sets are slightly unbalanced in favor of the s-functions and in disfavor of the d-functions, but various alternative basis sets may be derived by combining parts of the five parent sets. The convergence of radial expectation values is discussed.     

Gradient optimization of polarization exponents inab initio MO calculations on H2SO → HSOH and CH3SH → CH2SH2

Summary Analytical gradients were used to optimize the polarization function exponents in the 6-31G(d) and 6-31G(d, p) basis sets for the reactants, transition structures and products in the reactions H₂SO HSOH and CH₃SH CH₂SH₂. The optimizedd exponents on the heavy atoms change by ±10% in the course of the reactions and depend on the bonding of the heavy atoms. Thep exponents on the hydrogens change by as much as a factor of 5 and depend on the element to which the hydrogen is bonded and its valency. The effect of exponent optimization on the relative energies is small (±3 kcal/mol). With the 6-31G(d, p) basis set, optimization of the polarization exponents can make some of the bonds significantly more polar, as judged by the Mulliken charges.     

Medium-size Gaussian basis sets for hydrogen through argon

Summary Medium-sized Gaussian basis sets are reoptimized for the ground states of the atoms from hydrogen through argon. The composition of these basis sets is (4s), (5s), and (6s) for H and He, (9s5p) and (12s7p) for the atoms Li to Ne, and (12s8p) and (12s9p) for the atoms Na to Ar. Basis sets for the² P states of Li and Na, and the³ P states of Be and Mg are also constructed since they are useful in molecular calculations. In all cases, our energies are lower than those obtained previously with Gaussian basis sets of the same size.     

Revised relativistic basis sets for the 5d elements Hf–Hg

The relativistic double-zeta (dz) and triple-zeta (tz) basis sets for the 5d elements Hf–Hg have been revised for consistency with the recently optimized 4f basis sets. The new dz basis sets have 24 s functions instead of 22 s functions, and the new tz basis sets have 30 s functions instead of 29 s functions. New contraction patterns have been determined, including the 6p orbital.     

Large atomic natural orbital basis sets for the first transition row atoms

Large atomic natural orbital (ANO) basis sets are tabulated for the Sc to Cu atoms. The primitive sets are taken from the large sets optimized by Partridge, namely (21s13p8d) for Sc and Ti and (20s12p9d) for V to Cu. These primitive sets are supplemented with threep, oned, sixf, and fourg functions. The ANO sets are derived from configuration interaction density matrices constructed as the average of the lowest states derived from the 3d ⁿ 4s ² and 3d ⁿ⁺¹4s ¹ occupations. For Ni, the¹ S(3d ¹⁰) state is included in the averaging. The choice of basis sets for molecular calculations is discussed.     

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.     

Density matrix averaged atomic natural orbital (ANO) basis sets for correlated molecular wave functions

Summary Generally contracted basis sets for first row atoms have been constructed using the Atomic Natural Orbital (ANO) approach, with modifications for allowing symmetry breaking and state averaging. The ANOs are constructed by averaging over several atomic states, positive and negative ions, and atoms in an external electric field. The contracted basis sets give virtually identical results as the corresponding uncontracted sets for the atomic properties, which they have been designed to reproduce. The design objective has been to describe the ionization potential, the electron affinity, and the polarizability as accurately as possible. The result is a set of well-balanced basis sets for molecular calculations. The starting primitive sets are 8s4p3d for hydrogen, 9s4p3d for helium, and 14s9p4d3f for the heavier first row atoms.     

Near-Dirac–Hartree–Fock results for first-row atoms calculated with GTO basis sets

Relativistic basis sets for first-row atoms have been constructed by using the near-Hartree–Fock (nonrelativistic) eigenvectors calculated by Partridge. These bases generate results of near-Dirac–Hartree–Fock quality. Relativistic total and orbital energies, relativistic corrections to the total energy, and magnetic interaction energies for the first-row atoms have been presented. The smallest Gaussian expansions (13s8 p expansions) yield Dirac–Hartree–Fock total energies accurate through six significant digits, while the largest expansions (18s13p expansions) give these energies accurate through seven significant digits. These results are more accurate than some of the results reported earlier, particularly for the open-shell atoms, indicating that the basis employed is reasonably economical for relativistic calculations. © 1995 John Wiley & Sons, Inc.     

Energy-adjusted pseudopotentials for the rare earth elements

Nonrelativistic and quasirelativistic energy-adjusted pseudopotentials and optimized (7s6p5d)/[5s4p3d]-GTO valence basis sets for use in molecular calculations for fixed f-subconfigurations of the rare earth elements, La through Lu, have been generated. Atomic excitation and ionization energies from numerical HF, as well as SCF pseudopotential calculations using the derived basis sets, differ by less than 0.1 eV from numerical HF all-electron results. Corresponding values obtained from CI(SD), CEPA-1, as well as density functional calculations using the quasirelativistic pseudopotentials, are in reasonable agreement with experimental data.     

Theoretical studies of organometallic compounds. III. Structures and bond energies of FeCHn and FeCH (n = 1, 2, 3)

The geometries and dissociation energies for the Fe? C and C? H bonds of FeCH_n and FeCH (n = 1, 2, 3) have been calculated by ab initio quantum mechanical methods using different effective core potential models and Møller–Plesset perturbation theory. The HW3 ECP model, which has a configuration [core] (n?1)s², (n?1)p⁶, (n?1)d¹, (n)s^m for the transition metals, is clearly superior to the larger core LANL1DZ ECP model with the configuration [core] (n?1)d¹, (n)s^m. The Fe? C bond energies calculated at correlated levels using the HW3 ECP are in much better agreement with experiment than the LANL1DZ results. This effect is mainly due to the higher number of correlated electrons rather than the inclusion of the outermost core electrons in the Hartree–Fock calculation. At the PMP4/HW3TZ/6-31G(d)//MP2/HW3TZ/6-31G(d) level, the theoretically predicted Fe? C bond energies for FeCH are in the range of 80% of the experimental values and have nearly the same accuracy as all-electron calculations using large valence basis sets and the MCPF method for the correlation energy. © 1992 by John Wiley & Sons, Inc.     

Small gaussian basis sets for second-row atoms

6s-type and 4p-type gaussian basis sets are obtained for the second row atoms by fitting, using a least squares criterion, to 12s-type and 9p-type gaussian basis sets which are close to the self-consistent field atomic orbital wave functions. The small gaussian expansions are considered to be more suited for molecular calculations using double basis sets. The differences between these sets and the 10s-type, 6p-type and 9s-type, 5p-type are analysed. For molecular calculations using single gaussian basis sets the 10s-type and 6p-type would seem to be the best compromise.

---

Zusammenfassung Ein Basissatz von Gaußfunktionen vom 6s- bzw. 4p-Typ für Atome der zweiten Reihe wird erhalten, indem die Funktionen mit Hilfe des Kriteriums der kleinsten quadratischen Abweichung einem Satz von Gaußfunktionen vom 12s- bzw. 9p-Typ angepaßt werden; dabei ist der letztgenannte Satz der selbstkonsistenten Wellenfunktion aus Atomorbitalen stark angenähert. Die kürzeren Entwicklungen nach Gaußfunktionen werden für geeigneter bei Berechnungen mit zweifachen Basissätzen gehalten. Die Unterschiede zwischen diesen Sätzen und solchen vom 10s- bzw. 6p-Typ sowie vom 9s- und 5p-Typ werden untersucht. Für Molekülrechnungen mit einfachen Basissätzen von Gaußfunktionen scheint der Satz vom 10s- bzw. 6p-Typ den besten Kompromiß darzustellen.

---

Résumé On obtient des bases gaussiennes de type 6s et 4p pour les atomes de la seconde ligne par ajustement selon un critère de moindre carré à des bases gaussiennes de type 12s et 9p proches des orbitales atomiques SCF. Les petits développements en gaussiennes sont plus adaptés à des calculs moléculaires en bases doubles. Analyse des différences entre cas bases et les bases de types 10s et 6p, 9s et 5p. Pour des calculs moléculaires à base simple, 10s et 6p semble le meilleur compromis.

     

Long range interaction energies using Gaussian basis sets and a one center method

A one center method, based on the work of Karplus and Kolker, is discussed and used to calculate the induction energy, through O(R^?8), for the H(ls) – H⁺ interaction employing two types of Gaussian basis sets constructed from functions of the form {r^je^?αr2}. The effective hydrogen atom excitation energies and transition multipole moment matrix elements generated in these calculations are used to calculate the dispersion energy for the H(ls) – H(ls) interaction, through O(R^?10), and the R^?9 triple dipole energy corresponding to the interaction of three H(ls) atoms. The results indicate that Gaussian functions can form good basis sets for obtaining long range forces for a variety of multipole interaction energies.     

Ab initio study of hydrogen bonding in the phenol–water system

Three hydrogen-bonded minima on the phenol-water, C₆H₅OH—H₂O, potential energy surface were located with 3–21G and 6–31G** basis sets at both Hartree–Fock and MP2 levels of theory. MP2 binding energies were computed using large “correlation consistent” basis sets that included extra diffuse functions on all atoms. An estimate of the effect of expanding the basis set to the triple-zeta level (multiple f functions on carbon and oxygen and multiple d functions on hydrogen) was derived from calculations on a related prototype system. The best estimates of the electronic binding energies for the three minima are –7.8, –5.0, and –2.0 kcal/mol. The consequences of uncertainties in the geometries and limitations in the level of correlation recovery are analyzed. It is suggested that our best estimates will likely underestimate the complete basis set, full CI values by 0.1–0.3 kcal/mol. Vibrational normal modes were determined for all three minima, including an MP2/6–31G** analysis for the most strongly bound complex. Computational strategies for larger phenol–water complexes are discussed. © John Wiley & Sons, Inc.     

Energy-adjustedab initio pseudopotentials for the second and third row transition elements

Summary Nonrelativistic and quasirelativisticab initio pseudopotentials substituting the M^(Z–28)+-core orbitals of the second row transition elements and the M^(Z–60)+-core orbitals of the third row transition elements, respectively, and optimized (8s7p6d)/[6s5p3d]-GTO valence basis sets for use in molecular calculations have been generated. Additionally, corresponding spin-orbit operators have also been derived. Atomic excitation and ionization energies from numerical HF as well as from SCF pseudopotential calculations using the derived basis sets differ in most cases by less than 0.1 eV from corresponding numerical all-electron results. Spin-orbit splittings for lowlying states are in reasonable agreement with corresponding all-electron Dirac-Fock (DF) results.     

TiCl, TiH, and TiH+ bond energies: a test of a correlation-consistent Ti basis set

Correlation-consistent basis sets are developed for the Ti atom. The polarization functions are optimized for the average of the ³F and ⁵F states. One series of correlation-consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti ³F–⁵F separation and the dissociation energies of TiCl X⁴Φ, TiH X⁴Φ, and TiH⁺ X³Φ. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas–Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or are computed at the CASSCF level. The Ti³F–⁵F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH⁺ bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation more difficult. Received: 20 January 1999 / Accepted: 26 February 1999 / Published online: 7 June 1999     

A CASPT2 study of the valence and lowest Rydberg electronic states of benzene and phenol

Summary The valence excited states and the 3s, 3p, and 3d (united atom) Rydberg states of benzene and phenol have been obtained by the CASPT2 method, which computes a second-order perturbation correction to complete active space self-consistent field (CASSCF) energies. All non-zero dipole oscillator strengths are also computed, at the CASSCF level. For benzene, 16 singlet and 16 triplet states with excitation energies up to ca. 7.86 eV (63 400 cm^–1) are obtained. Of these, 12 singlet and three triplet energies are experimentally known well enough to allow meaningful comparison. The average error is around 0.1 eV. The highest of these singlet states (2¹ E_2g) is the highest valence * state predicted by elementary -electron theory. Its energy is then considerably lower than has been suggested from laser flash experiments, but in perfect agreement with a reinterpretation of that experiment. For phenol, 27 singlet states are obtained, in the range 4.53–7.84 eV (63 300 cm^–1). Only the lowest has a well-known experimental energy, which agrees with the computed result within 0.03 eV. The ionization energy is in error by 0.05 eV.     

An improved LCAO interpolation scheme for energy band structures. Application to four compounds (ScN,ScP, TiN,ZrN) crystallizing in the sodium chloride structure

An improved version of the LCAO interpolation scheme using metal s-, p-, d-, and non metal s-, and p-basis functions is presented for transition metal compounds with sodium chloride structure. This method enables us to interpolate with reasonable accuracy occupied bands as well as unoccupied energy bands up to 0.9 Rydberg above the Fermi level for the compounds ScN, ScP, TiN and ZrN. Due to the limited basis, problems arise however with bands of predominantly transition metal f or non metal d character lying in this energy range - as is the case for ScP.Optimized parameter sets are presented for the compounds ScN, ScP, TiN and ZrN. They were used for the calculation of the imaginary part of the complex dielectric function, ₂(), as will be shown in two forthcoming papers.     

Vector coupling coefficients for calculations of transition-metal atoms and ions by the SCF coupling operator method

We derived the necessary conditions to which the vector coupling coefficients (VCC ) a and b describing atomic L,S-multiplets of the configurations d^N (1 ≤ N ≤ 9), should satisfy. Special attention is paid to the states of non-Roothaan type for which VCC depend on the choice of degenerate d-orbitals basis set determined within the accuracy up to an orthogonal transformation u. It is shown that for such states the direct sum of matrices ‖a‖ and ‖b‖ must be the non-symmetric matrix. Obtained VCC were used for the ab initio calculations (basis set (14s9p5d)/[8s4p2d] from [15]) on first-row transition atoms (from Sc to Cu) to compare to similar calculations [16], in which the Peterson's VCC have been used, and with calculations [15] carried out by the atomic SCF program [4] as well.