Expansion Formulae for Two-center Integer and Noninteger <Emphasis Type="Italic">n</Emphasis> STO Charge Densities and their Use in Evaluation of Multi-center Integrals

Using complete orthonormal sets of Ψ^α-exponential type orbitals (Ψ^α-ETOs, α =1, 0, −1, −2, ...) introduced by the author, the series expansion formulae are derived for the two-center integer and noninteger n STO (ISTO and NISTO) charge densities in terms of integer n STOs at a third center. The expansion coefficients occurring in these relations are presented through the two-center overlap integrals between STOs with integer and noninteger principal quantum numbers. The general formulae obtained for the STO charge densities are utilized for the evaluation of two-center Coulomb and hybrid integrals of NISTOs appearing in the Hartee–Fock–Roothaan approximation. The final results are expressed in terms of both the overlap integrals and the one-center basic integrals over integer n STOs. It should be noted that the result for the multi-center multielectron integrals with two-center noninteger n STO charge densities presented in this paper were not appeared in our past publications.     

Gaussian vs. Slater representations of d orbitals: An information theoretic appraisal based on both position and momentum space properties

A detailed appraisal of Gaussian-type orbital (GTO) and Slater-type orbital (STO) expansions of 3d orbitals is carried out for the ²S state of copper—a case that should be maximally unfavorable for STOs. The appraisal is based on a wide variety of both position and momentum space properties and utilizes an information theoretic quality assessment technique. It is found that GTO expansions are not as useful as STO expansions for the prediction of 〈p⁸〉, 〈p⁷〉, and 〈r^?6〉 because these properties probe the functional deficiencies of GTOs at small r and large p. On the other hand, GTO expansions can predict accurate values of large r properties like 〈r⁸〉 despite the fact that their position space asymptotic decay is too fast. Unlike the case of s orbitals in helium, there does not seem to be any consistent ordering between accuracy in position space and accuracy in momentum space. The quality measures are found to be very useful for pinpointing the deficiencies of various expansions. This information enables us to construct easily a new GTO and a new STO expansion that are more accurate than any of the others in the literature. It is suggested that one STO is worth no more than two GTOs in the case of d orbitals.     

Calculation of electric multipole moment integrals with the different screening parameters via the Fourier transform method

Three‐center electric multipole moment integrals over Slater‐type orbitals (STOs) can be evaluated by translating the orbitals on one center to the other and reducing the system to an expansion of two‐center integrals. These are then evaluated using Fourier transforms. The resulting expression depends on the overlap integrals that can be evaluated with the greatest ease. They involve expressions for STO with different screening parameters that are known analytically. This work gives the overall expressions analytically in a compact form, based on Gegenbauer polynomials. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

On the use of Gaussian shell type basis orbitals for single-center expansions. I. Evaluation of integrals

Formulas are derived for all Hamiltonian integrals required for molecular computations using a novel basis for single-center expansions. The basis orbitals depend exponentially upon α(r ? ρ)² where r and ρ are the distance from center to electron and to a variationally scaled spherical shell, respectively. Comparisons are made between these so-called Gaussian shell orbitals (GSO ) and the conventional GTO and STO bases for single-center calculations. A preliminary comparison on H using a single GSO , a non-integer STO , and a GTO gives the optimized energies: ?0.51089 a.u., ?0.50504 a.u., and ?0.50422 a.u., respectively.     

Calculation of two-center overlap integral in molecular coordinate system over Slater type orbital using Löwdin α-radial and Guseinov rotation–angular functions

By use of Löwdin and Guseinov relations for the radial and angular part of two-center overlap integrals, respectively, the computer calculations of overlap integrals over Slater type orbitals (STOs) in molecular coordinate system are performed. The results of calculations are valid for arbitrary principal quantum numbers, screening constants and location of STOs. Excellent agreement with benchmark results and stability of the technique are demonstrated.     

Evaluation of integrals over STOs on different centers and the complementary convergence characteristics of ellipsoidal‐coordinate and zeta‐function expansions

One‐electron integrals over three centers and two‐electron integrals over two centers, involving Slater‐type orbitals (STOs), can be evaluated using either an infinite expansion for 1/r₁₂ within an ellipsoidal‐coordinate system or by employing a one‐center expansion in spherical‐harmonic and zeta‐function products. It is shown that the convergence characteristics of both methods are complimentary and that they must both be used if STOs are to be used as basis functions in ab initio calculations. To date, reports dealing with STO integration strategies have dealt exclusively with one method or the other. While the ellipsoidal method is faster, it does not always converge to a satisfactory degree of precision. The zeta‐function method, however, offers reliability at the expense of speed. Both procedures are described and the results of some sample calculation presented. Possible applications for the procedures are also discussed. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 1–13, 1999     

Four‐center integrals for Gaussian and exponential functions

The shift operator technique is used for deriving, in a unified manner, the master formulas for the four‐center repulsion integrals involving Gaussian (GTO), Slater (STO), and Bessel (BTO) basis functions. Moreover, for the two classes of exponential‐type functions (ETO), i.e., STO and BTO, we give the expressions corresponding to both the Gauss and Fourier transforms. From the comparison of the master formulas of GTO and ETO, we conclude that ETO can perform more efficiently than GTO, and we remark the points where the effort must be focused to carry out this possibility. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 16–28, 2001     

Spherically symmetrical properties of two-center overlap integrals over arbitrary atomic orbitals and translation coefficients for Slater-type orbitals

The orthogonality relations are derived for the rotation coefficients of two-center overlap integrals over arbitrary atomic orbitals (AAOs) and expansion coefficients for translation of Slater-type orbitals (STOs). Using these formulas, a very interesting theorem regarding the angular dependence is established. If we add the products of all the overlap integrals or all the translation coefficients with the same n and l values, but different m values, the result is independent of orientation. The final results are of a simple structure and are, therefore, especially useful for machine computations of multielectron multicenter molecular integrals by expanding one- and two-center electron charge density over STOs in terms of STOs about a new center.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals. II. Two‐center expansion method

Using expansion formulas for the charge‐density over Slater‐type orbitals (STOs) obtained by the one of authors [I. I. Guseinov, J Mol Struct (Theochem) 1997, 417, 117] the multicenter molecular integrals with an arbitrary multielectron operator are expressed in terms of the overlap integrals with the same screening parameters of STOs and the basic multielectron two‐center Coulomb or hybrid integrals with the same operator. In the special case of two‐electron electron‐repulsion operator appearing in the Hartree–Fock–Roothaan (HFR) equations for molecules the new auxiliary functions are introduced by means of which basic two‐center Coulomb and hybrid integrals are expressed. Using recurrence relations for auxiliary functions the multicenter electron‐repulsion integrals are calculated for extremely large quantum numbers. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 117–125, 2001     

Overlap Integrals Between Irregular Solid Harmonics and STOs Via the Fourier Transform Methods

In this paper, a unified analytical and numerical treatment of overlap integrals between Slater type orbitals (STOs) and irregular Solid Harmonics (ISH) with different screening parameters is presented via the Fourier transform method. Fourier transform of STOs is probably the simplest to express of overlap integrals. Consequently, it is relatively easy to express the Fourier integral representations of the overlap integrals as finite sums and infinite series of STOs, ISHs, Gegenbauer, and Gaunt coefficients. The another mathematical tools except for Fourier transform have used partial-fraction decomposition and Taylor expansions of rational functions. Our approach leads to considerable simplification of the derivation of the previously known analytical representations for the overlap integrals between STOs and ISHs with different screening parameters. These overlap integrals have also been calculated for extremely large quantum numbers using Gegenbauer, Clebsch-Gordan and Binomial coefficients. The accuracy of the numerical results is quite high for the quantum numbers of Slater functions, irregular solid harmonic functions and for arbitrary values of internuclear distances and screening parameters of atomic orbitals.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals I. Single‐center expansion method

The multicenter charge‐density expansion coefficients [I. I. Guseinov, J Mol Struct (Theochem) 417 , 117 (1997)] appearing in the molecular integrals with an arbitrary multielectron operator were calculated for extremely large quantum numbers of Slater‐type orbitals (STOs). As an example, using computer programs written for these coefficients, with the help of single‐center expansion method, some of two‐electron two‐center Coulomb and four‐center exchange electron repulsion integrals of Hartree–Fock–Roothaan (HFR) equations for molecules were also calculated. Accuracy of the results is quite high for the quantum numbers, screening constants, and location of STOs. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 146–152, 2000     

Evaluation of expansion coefficients for translation of Slater‐type orbitals using complete orthonormal sets of exponential‐type functions

By the use of exponential‐type functions (ETFs) the simpler formulas for the expansion of Slater‐type orbitals (STOs) in terms of STOs at a displaced center are derived. The expansion coefficients for translation of STOs are presented by a linear combination of overlap integrals. The final results are of a simple structure and are, therefore, especially useful for machine computations of arbitrary multielectron multicenter molecular integrals over STOs that arise in the Hartree–Fock–Roothaan approximation and also in the Hylleraas correlated wave function method for the determination of arbitrary multielectron properties of atoms and molecules. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 126–129, 2001     

Use of Cartesian coordinates in evaluation of multicenter multielectron integrals over Slater type orbitals and their derivatives

Using addition theorems for interaction potentials and Slater type orbitals (STOs) obtained by the author, and the Cartesian expressions through the binomial coefficients for complex and real regular solid spherical harmonics (RSSH) and their derivatives presented in this study, the series expansion formulas for multicenter multielectron integrals of arbitrary Coulomb and Yukawa like central and noncentral interaction potentials and their first and second derivatives in Cartesian coordinates were established. These relations are useful for the study of electronic structure and electron-nuclei interaction properties of atoms, molecules, and solids by Hartree–Fock–Roothaan and correlated theories. The formulas obtained are valid for arbitrary principal quantum numbers, screening constants and locations of STOs.     

Calculation of Multicenter Nuclear Attraction and Electron Repulsion Integrals over Slater Orbitals by Fourier Transform Method Using Gegenbauer Polynomials

In this study, using complete orthonormal sets of exponential type orbitals (ETOs), a single closed analytical relation is derived for a large number of different expansions of overlap integrals over Slater type orbitals (STOs) with the same screening parameters in terms of Gegenbauer coefficients. The general formula obtained for the overlap integrals is utilized for the evaluation of multicenter nuclear attraction and electron repulsion integrals appearing in the Hartree–Fock–Roothaan equations for molecules. The formulas given in this study for the evaluation of these multicenter integrals show good rate of convergence and great numerical stability under wide range of quantum numbers, scaling parameters of STOs and internuclear distances.     

Importance of the basis set for the spin-state energetics of iron complexes

We have performed a systematic investigation of the influence of the basis set on relative spin-state energies for a number of iron compounds. In principle, with an infinitely large basis set, both Slater-type orbital (STO) and Gaussian-type orbital (GTO) series should converge to the same final answer, which is indeed what we observe for both vertical and relaxed spin-state splittings. However, we see throughout the paper that the STO basis sets give consistent and rapidly converging results, while the convergence with respect to the basis set size is much slower for the GTO basis sets. For example, the large GTO basis sets that give good results for the vertical spin-state splittings of compounds 1-3 (6-311+G**, Ahlrichs VTZ2D2P) fail for the relaxed spin-state splittings of compound 4 (where 1 is Fe-(PyPepS)2 (PyPepSH 2 = N-(2-mercaptophenyl)-2-pyridinecarboxamide), 2 is Fe(tsalen)Cl (tsalen = N, N'-ethylenebis-(thiosalicylideneiminato)), 3 is Fe(N(CH2-o-C6H4S) 3)(1-Me-imidazole), and 4 is FeFHOH). Very demanding GTO basis sets like Dunning's correlation-consistent (cc-pVTZ, cc-pVQZ) basis sets are needed to achieve good results for these relaxed spin states. The use of popular (Pople-type) GTO, effective core potentials basis set (ECPB), or mixed ECPB(Fe):GTO(rest) basis sets is shown to lead to substantial deviations (2-10 kcal/mol, 14-24 kcal/mol for 3-21G), in particular for the high spin states that are typically placed at too low energy. Moreover, the use of an effective core potential in the ECPB basis sets results in spin-state splittings that are systematically different from the STO-GTO results.     

New program for molecular calculations with Slater‐type orbitals

A new program for computing all the integrals appearing in molecular calculations with Slater‐type orbitals (STO) is reported. This program follows the same philosophy as the reference pogram previously reported but introduces two main changes: Local symmetry is profited to compute all the two‐electron integrals from a minimal set of seed integrals, and a new algorithm recently developed is used for computing the seed integrals. The new code reduces between one and two orders of magnitude the computational cost in most polyatomic systems. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 148–153, 2001     

Calculation of Two-center Nuclear Attraction Integrals over Slater Type Orbitals in Molecular Coordinate System

A closed analytical relation is derived for the two-center nuclear attraction integrals over Slater type orbitals (STOs) in terms of binomial coefficients. This formula can be used in highly accurate calculations of the nuclear attraction integrals. The relationships obtained are valid for arbitrary values of quantum numbers and screening constants of STOs and location of nuclei.     

Calculation of three-center electric and magnetic multipole moment integrals using translation formulas for Slater-type orbitals

 By the use of translation formulas for the expansion of Slater-type orbitals (STOs) in terms of STOs at a new origin, three-center electric and magnetic multipole moment integrals are expressed in terms of two-center multipole moment integrals for the evaluation of which closed analytical formulas are used. The convergence of the series is tested by calculating concrete cases. Computer results with an accuracy of 10⁻⁷ are obtained for 2^ν– pole electric and magnetic multipole moment integrals for 1≤ν≤5 and for arbitrary values of screening constants of atomic orbitals and internuclear distances. Received: 28 October 1999 / Accepted: 15 February 2000 / Published online: 5 June 2000     

Unified treatment of multicenter integrals of integer and noninteger u Yukawa-type screened Coulomb type potentials and their derivatives over Slater orbitals

Multicenter integrals with integer and noninteger values of indices u of Yukawa-type screened Coulomb-type potentials (SCTPs) f(u)(eta,r)=r(u-1)e(-etar) and their derivatives over Slater-type orbitals (STOs) are evaluated using series expansion formulas obtained from the expansions in terms of complete orthonormal sets of Psi(alpha)-exponential-type orbitals (Psi(alpha)-ETOs, alpha=1, 0, -1, -2,...). The final results obtained are valid for the arbitrary values of parameters of SCTPs and STOs.