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.     

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

Generally contracted basis sets for the first row transition metal atoms Sc-Zn 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 the three electronic configurationsd ⁿ ,d ^n–1 s, andd ^n–2 s ² for the neutral atom as well as the ground state for the cation and the ground state atom in an external electric field. The primitive sets are 21s15p10d6f4g. Contraction to 6s5p4d3f2g yields results that are virtually identical to those obtained with the corresponding uncontracted basis sets for the atomic properties, which they have been designed to reproduce. Slightly larger deviations are obtained with the 5s4p3d2f1g for the polarizability, while energetic properties still have only small errors. The design objective has been to describe the ionization potential, the polarizability and the valence spectrum as accurately as possible. The result is a set of well-balanced basis sets for molecular calculations, which can be used together with basis sets of the same quality for the first and second row atoms.     

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.     

The excited states of pyrazine: A basis set study

Summary Multiconfigurational second order perturbation theory (CASSCF/CASPT2) has been used to investigate the dependence of computed valence excitation energies and transition moments on the basis sets. Pyrazine has been selected as the test molecule. Atomic normal orbital (ANO) type basis sets are used throughout. Contractions of the structure (4s3p1d/2s) are found to be an optimal compromise between the quality and the size of the calculations and are capable of yielding results virtually identical to more extended basis sets.     

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.     

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.     

Reliability of atomic natural orbital basis sets in calculations involving pseudopotentials

The performance of Atomic Natural Orbital (ANO) basis sets for calculations involving nonempirical core pseudopotentials has been studied by comparing the results for atomic and molecular nitrogen obtained using contracted ANO basis sets with those obtained using both the primitive set and a segmented one. The primitive set has been optimized at the SCF level for atomic N treated as a five-electron pseudo-atom, and consists of 7s and 7p primitive GTOs supplemented by 2d and 1f GTOs optimized at the CI level. From this primitive set three contracted [3s 3p 2d 1f] sets have been obtained. The first one has been derived from the ANOs of the neutral atom, the second has been obtained from an averaged density matrix and the third one is a segmented set. For the atom, the segmented set gives a zero contraction error at the SCF level as it must be in valence-only calculations. The ANO basis sets show some small contraction error at the SCF level but perform better in CI calculations. However, for the diatomic N₂ molecule the ANO basis sets exhibit a rather large contraction error in the calculated SCF energy. A detailed analysis of the origin of this error is reported, which shows that the conventional strategy used to derive ANO basis sets does not work very well when pseudopotentials are involved.     

Configuration interaction calculations on the 2P ground state of boron atom and C+ using Slater orbitals

Configuration Interaction (CI) calculations on the ground ²P state of boron atom are presented using a wave function expansion constructed with L‐S eigenfunction configurations of s‐, p‐, and d‐Slater orbitals. Two procedures of optimization of the orbital exponents have been investigated. First, CI(SD) calculations including few types of configurations and full optimization of the orbital exponents led to the energy ?24.63704575 a.u. Second, full‐CI (FCI) calculations including a large number of configuration types using a fixed set of orbital exponents for all configurations gave ?24.63405222 a.u. using the basis [4s3p2d] and 2157 configurations, and to an improved result of ?24.64013999 a.u. for 3957 configurations and a [5s4p3d] basis. This last result is better than earlier calculations of Schaefer and Harris (Phys Rev 1968, 167, 67), and compares well with the recent ones from Froese Fischer and Bunge (personal communication). In addition, using the same wave functions, CI calculations of the boron isoelectronic ion C⁺ have been performed obtaining an energy of ?37.41027598 a.u. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Model potentials suitable for calculations with slater-type basis for Sc through Zn

Model potentials appropriate for molecular calculations with Slater-type (ST) basis sets have been generated for the first-row transition-metal atoms. Two sets of model potentials are presented. The first one has been optimized using standard 2ζ ST basis sets. The second is consistent with a reduced ST basis set. The reduced bases have been obtained by means of a new algorithm, whose results are compared to those found with the method of Y. Sakai and S. Huzinaga (J. Chem. Phys.76, 2537 (1982)). The comparison shows that the new approach leads to significant improvements in the overall results. Two different valence shells have been investigated, one formed by the 3s, 3p, 3d, and 4s AOs (SPDS), and the other one formed by the 3d and 4s AOs (DS). The model potentials presented here describe these valence shells in good agreement with the all-electron calculations taken as reference, the SPDS calculations being uniformly more accurate. Special attention has been paid to the transferability of the model potentials to electronic states with different orbital occupation. It is shown that the reported potentials have a wide transferability of this sort, being thus suitable for calculating the electronic structure of transition-metal compounds with STO 2ζ quality.     

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.     

Coupling coefficients for systems with two open electronic shells: Transition metal ions with pM dN configuration

We derived the necessary conditions to which the vector coupling coefficients (VCC) a_mn and b_mn, describing atomic L, S-multiplets of p^Md^N and d^Ns¹ configurations (1 ? N ? 9, 1 ? M ? 5), should satisfy. It is shown that for two-open-shell systems under consideration the unknown VCC should satisfy not only usual restrictions resulting from the spheric symmetry, but also some additional equation introduced in the present paper in the form of a postulate. VCC obtained were used for the ab initio calculations (by the general SCF coupling operator method) of several transitionmetal atoms and ions with electronic configurations 3d¹4p¹, 3p⁴3d³, and 3p⁵3d³. To check the presented theory, we carried out a detailed comparison between these results and analogous data, obtained by the atomic Roothaan-Hartree-Fock method [2]. © 1993 John Wiley & Sons, Inc.     

Compact contracted basis sets for third-row atoms: Ga–Kr

The (14s11p5d) primitive basis set of Dunning for the third-row main group atoms Ga-Kr has been contracted [6s4p1d]. The core functions have been relatively highly contracted while those which represent the valence region have been left uncontracted to maintain flexibility. Calculations with the [6s4p1d] contraction are reported for a variety of molecules involving third-row atoms. This basis set is found to satisfactorily reproduce experimental properties such as geometric configurations, dipole moments, and vibrational frequencies for a range of molecules. Comparisons are made with the performance of the uncontracted basis set. Polarization functions for the contracted basis set are reported and performance of the basis set with and without polarization functions is examined. A relaxation of the [6s4p1d] contraction to [9s6p2d] for higher level evergy calculations is also presented.     

Coupling between the convergence behavior of basis set and electron correlation: a quantitative study

The coupling between improvement of the basis set and the valence electron correlation method has been studied quantitatively for the total atomization energies (TAEs) of a number of small molecules, using basis sets of up to [7s6p5d4f3g2h/5s4p3d2f1g] quality. Very significant coupling is found to exist. Using a scaled basis set extrapolation beyond [6s5p4d3f2g/ 4s3p2d1f] at the MP2 or CCSD level, mean absolute errors of 0.18 and 0.15 kcal/mol, respectively, can be obtained for the TAEs of a number of small polyatomic molecules, compared to 0.12 kcal/mol using CCSD(T) throughout. Received: 7 February 1997 / Accepted: 6 May 1997