Substituent effects in second-row molecules: Basis set performance in calculations of normal valency phosphorus and sulfur compounds

Ab inito molecular orbital calculations of the phosphorus- and sulfur-containing series PH₂X, PH₃X⁺, SHX, and SH₂X⁺ (X = H, CH₃, NH₂, OH, F) have been carried out over a range of Gaussian basis sets and the results (optimized geometrical structures, relative energies, and electron distributions) critically compared. As in first-row molecules there are large discrepancies between substituent interaction energies at different basis set levels, particularly in electron-rich molecules; use of basis sets lower than the supplemented 6-31G basis incurs the risk of obtaining substituent stabilizations with large errors, including the wrong sign. Only a small part of the discrepancies is accounted for by structural differences between the optimized geometries. Supplementation of low level basis sets by d functions frequently leads to exaggerated stabilization energies for π-donor substituents. Poor performance also results from the use of split valence basis sets in which the valence shell electron density is too heavily concentrated in diffuse component of the valence shell functions, again likely to occur in electron-rich molecules. Isodesmic reaction energies are much less sensitive to basis set variation, but d function supplementation is necessary to achieve reliable results, suggesting a marginal valence role for d functions, not merely polarization of the bonding density. Optimized molecular geometries are relatively insensitive to basis set and electron population analysis data, for better-than-minimal bases, are uniform to an unexpected degree.     

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.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for divalent and tetravalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modeling divalent (5f ⁿ⁺¹ occupation with n = 5–13 for Pu–No) respectively tetravalent (5f ^n-1 occupation with n = 1–9 for Th–Cf) actinides together with corresponding core-polarization potentials have been generated. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets as well as 2f1g correlation functions have been derived and contracted to polarized double, triple, and quadruple zeta quality. Corresponding smaller (4s4p) and (5s5p) respectively (4s4p3d) and (5s5p4d) basis sets suitable for calculations on actinide(II) respectively actinide(IV) ions in crystalline solids form subsets of these basis sets designed for calculations on molecules. Results of Hartree–Fock test calculations for actinide di- and tetrafluorides show a satisfactory agreement with calculations using 5f-in-valence pseudopotentials. Electronic Supplementary Material The online version of this article doi: contains supplementary material, which is available to authorized users.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for trivalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modelling trivalent actinides, corresponding to a near-integral 5f ⁿ occupation (n = 0–14 for Ac–Lr), have been generated. Energy-optimized (6s5p4d), (7s6p5d), and (8s7p6d) primitive valence basis sets contracted to polarized double to quadruple zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller basis sets (4s4p3d), (5s5p4d), and (6s6p5d) suitable for calculations on actinide(III) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Results of Hartree–Fock test calculations for actinide(III) monohydrates and actinide trifluorides show a satisfactory agreement with corresponding calculations using 5f-in-valence pseudopotentials. Even in the beginning of the actinide series, where the 5f shell is relatively diffuse, only quite acceptable small deviations occur as long as the 5f-shell does not participate significantly in covalent bonding. Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Quasirelativistic energy-consistent 4f-in-core pseudopotentials for tetravalent lanthanide elements

Quasirelativistic energy-consistent 4f-in-core pseudopotentials modeling tetravalent lanthanides (4f ^n?1 occupation with n = 1, 2, 3, 8, 9 for Ce, Pr, Nd, Tb, Dy) have been adjusted. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets contracted to polarized double- to quadruple-zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller (4s4p3d) and (5s5p4d) basis sets suitable for calculations on lanthanide(IV) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Calculations for lanthanide tetrafluorides using the 4f-in-core pseudopotentials at the Hartree–Fock level show satisfactory agreement with calculations using 4f-in-valence pseudopotentials. For cerium tetrafluoride the experimental bond length is well reproduced using the 4f-in-core pseudopotential at the coupled-cluster level with single and double excitation operators and a perturbative estimate of triple excitations. For cerium dioxide 4f-in-core and 4f-in-valence pseudopotential calculations agree quite well, if a proper f basis set instead of f polarization functions is applied.     

Pseudopotential study of the rare earth monohydrides,monoxides and monofluorides

Nonrelativistic and quasirelativistic energy-adjusted pseudopotentials for fixed 4f subconfigurations of the rare earth elements La through Lu together with corresponding optimized valence basis sets have been used in SCF and CI(SD) calculations to determine the spectroscopic constants for the energetically low lying superconfigurations of the lanthanide monohydrides, monoxides and monofluorides. The experimentally observed trends in dissociation energies, bond lengths and vibrational frequencies for the ground states of the calculated superconfigurations of the monoxides and monofluorides are well reproduced. The results for the monohydrides are mainly predictions.     

The expansion of hydrogen states in Gaussian orbitals

The convergence properties of Gaussian orbitals are studied by considering a very simple system, the hydrogen atom. We have variationally optimized even-tempered basis sets containing up to 60 s functions for the ground state and the first excited S state of the hydrogen atom, to an accuracy of 10^–15E_h. In addition, we have freely optimized the exponents in basis sets containing up to 12 Gaussians. We have studied the convergence of the total energy, the kinetic energy, the extent of the atom as measured by r², and the Fermi-contact interaction at the nucleus in these basis sets as well as in basis sets augmented with additional diffuse or steep functions.     

Contracted Gaussian-type basis functions revisited. IV. Atoms Rb to Xe

We report minimal-type contracted Gaussian-type function (GTF) sets, #n=(n3333/n33/n3) with n=5 and 6, #7= (74333/743/74), and #8= (84333/843/75), for the fourth-row atoms from Rb to Xe. Test calculations are performed on the Ag₂ molecule. Spectroscopic constants given by split valence sets derived from #5 and #6 are a little contaminated by basis set superposition error. However, we find that the fully valence split #8 set, (8433111/84111/711111), yields essentially the same results as a large GTF set, (22s15p12d), with a general contraction, when p-, d-, and f-type polarization functions are augmented. The present #7 and #8 CGTF sets are recommended for ab initio molecular calculations including fourth-row atoms. Received: 15 January 2002 / Accepted: 16 April 2002 / Published online: 24 June 2002     

A note on the importance of d orbitals in the calculation of effective interactions and the excitation spectra of atoms and molecules

The focal point of our discussion is the examination of truncated basis sets used in obtaining an accurate first principles clculation of the effective valence shell Hamiltonian by the canonical transformation-cluster expansion approasch. Subsequent diagonalization of this effecitve valence shell hamiltonian yields the valence shell transition energies. A detailed analysis of numerical results obtained using a number of different basis sets of hydrogen-like orbitals together with rigorous symmetry arguments celarly demonstrates the special role played by d orbitals in computing the ³P → ¹D transition energy in carbon. The failure of early attempts to calculate the effective Hamiltonian for ethylene from first principles is examined in the light of recent ab initio calculations on ethylene involving d orbitals and the computations reported in this paper. We conclude that accurate calculations of the effective valence shell Hamiltonian for molecules must consider d orbitals in the excited orbital basis set.     

HF and MP2 calculations on CN−, N2, AlF,SiO, PN,SC, ClB,and P2 using correlated molecular wave functions

Contracted basis sets of double zeta valence quality plus polarization functions (DZP) and augmented DZP basis sets, which were recently constructed for the first‐ and second‐row atoms, are applied to study the electronic ground states of the diatomic molecules CN^?, N₂, AlF, SiO, PN, SC, ClB, and P₂. At the Hartree–Fock (HF) and/or Møller–Plesset second‐order (MP2) levels, total and molecular orbital energies, dissociation energies, bond lengths, harmonic vibrational frequencies, and dipole moments are calculated and compared with available experimental data and with the results obtained from correlation consistent polarized valence basis sets of Dunning's group. For N₂, calculations of polarizabilities at the HF and MP2 levels with the sets presented above are also done and compared with results reported in the literature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Core correlating basis functions for elements 31–118

Gaussian functions for correlation of all core shells of elements from Z = 31 to Z = 118 have been optimized in relativistic singles and doubles CI calculations, performed on the shell of highest angular momentum for each principal quantum number. The SCF functions were derived from the double-zeta, triple-zeta, and quadruple-zeta basis sets previously optimized by the author. Only those Gaussian functions that are not represented in the SCF basis sets were optimized. The functions are available from the Dirac program web site, .     

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.     

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.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for pentavalent and hexavalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modeling pentavalent (5f ^n?2 occupation with n = 2–6 for Pa–Am) and hexavalent (5f ^n?3 occupation with n = 3–6 for U–Am) actinides have been adjusted. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets contracted to polarized double- to quadruple-zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller basis sets (4s4p3d) and (5s5p4d) suitable for calculations on actinide(V) and actinide(VI) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Calculations using the Hartree–Fock and the coupled-cluster method with single and double excitation operators and a perturbative estimate of triple excitations for actinide pentafluorides show satisfactory agreement with calculations using 5f-in-valence pseudopotentials and experimental data, respectively. However, in the hexavalent case the 5f-in-core approximation seems to reach its limitations except for hexavalent uranium (5f⁰), where results for both uranium hexafluoride and the uranyl ion deviate only slightly from the 5f-in-valence reference data.