Unified treatment for the evaluation of arbitrary multielectron multicenter molecular integrals over Slater‐type orbitals with noninteger principal quantum numbers

The expansion formula has been presented for Slater‐type orbitals with noninteger principal quantum numbers (noninteger n‐STOs), which involves conventional STOs (integer n‐STOs) with the same center. By the use of this expansion formula, arbitrary multielectron multicenter molecular integrals over noninteger n‐STOs are expressed in terms of counterpart integrals over integer n‐STOs with a combined infinite series formula. The convergence of the method is tested for two‐center overlap, nuclear attraction, and two‐electron one‐center integrals, due to the scarcity of the literature, and fair uniform convergence and great numerical stability under wide changes in molecular parameters is achieved. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Use of analytical relations in evaluation of exponential integral functions

Analytical formulas through the initial values suitable for numerical computation are developed for the exponential integral functions E_n (x). The relationships obtained are numerically stable for all values of n and for x < 1. Numerical results are also given.PACS No: 31.15.+q, 31.20.Ej AMS subject classification:81-V55, 81V45     

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.     

Analytical evaluation of three-center nuclear-attraction integrals over s-Slater orbitals for asymmetrical conformations of the centers

Analytical formulas for three-center nuclear-attraction integrals over Slater orbitals are given for any location of the three atomic centers. In the mathematical derivations the Neumann expansion has been used and new general auxiliary integrals which depend on the elliptical coordinates of one of the centers are defined. The orbital exponents within the integrals may be different.     

A new algorithm to evaluate bielectronic integrals with 1s Slater‐type orbitals obtained by the integral transforms

This article presents a variation of the integral transform method to evaluate multicenter bielectronic integrals (12|34), with 1s Slater‐type orbitals. It is proved that it is possible to define, out of the expression of (12|34) given by the integral transform method, a function F(q) that has the property of having a unique Q, such that F(Q) = (12|34). Therefore, F(q) may be used to calculate (12|34). It is shown that the evaluation of F(Q) turns out to be simpler than the three‐dimensional integral involved in the calculation of (12|34), and an algorithm is presented to calculate Q. The results show that relative errors on the order of 10^?3 or lower are obtained very efficiently. In addition, it is shown that the proposed algorithm is very stable. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Evaluation of the Boys Function using Analytical Relations

The analytical relations for Boys function F _m (x) are presented. These relations are useful in the fast and more accurate calculations of multicenter molecular integrals over Gaussian type orbitals (GTOs). The formulas obtained are numerically stable for all values of m and x.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals. III. Auxiliary functions \documentclass{article}\pagestyle{empty}\begin{document}$Q_{nn'}^{q}$\end{document} and \documentclass{article}\pagestyle{empty}\begin{document}$G_{-nn'}^{q}$\end{document}

The auxiliary functions $Q_{nn'}^{q}(p,pt)$ and $G_{-nn'}^{q}(p_{a},p,pt)$ which are used in our previous paper [Guseinov, I. I.; Mamedov, B. A. Int J Quantum Chem 2001, 81, 117] for the computation of multicenter electron‐repulsion integrals over Slater‐type orbitals (STOs) are discussed in detail, and the method is given for their numerical computation. The present method is suitable for all values of the parameters p_a, p, and pt. Three‐ and four‐center electron‐repulsion integrals are calculated for extremely large quantum numbers using relations for auxiliary functions obtained in this paper. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Use of Recursion and Analytical Relations in Evaluation of Hypergeometric Functions Arising in Multicenter Integrals with Noninteger <Emphasis Type="Italic">n</Emphasis> Slater Type Orbitals

In this work we present the new recursion and analytical relations for the calculation of hypergeometric functions F(1,b;c;z) occurring in multicenter integrals of noninteger n Slater type orbitals. The formulas obtained are numerically stable for 0 < z < 1 and all integer and noninteger values of parameters b and c The Author cordially congratulates Prof. I.I. Guseinov on his 70th birthday     

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     

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     

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.     

Use of Filter-Steinborn B and Guseinov $${Q_{ns}^q }$$ auxiliary functions in evaluation of two-center overlap integrals over Slater type orbitals

An efficient method for computing overlap integral over Slater type orbitals based on the B Filter-Steinborn and Guseinov auxiliary functions is presented. The final results are expressed through the binomial coefficients with the help of which the overlap integrals can be evaluated efficiently and accurately. The results of calculation are in good agreement with those obtained by other method for arbitrary principal quantum numbers and different screening constants. An erratum to this article can be found at     

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     

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.     

Analytic evaluation of two‐center STO electron repulsion integrals via ellipsoidal expansion

Extant analytic methods for evaluating two‐center electron repulsion integrals in a Slater‐type orbital (STO) basis using ellipsoidal coordinates and the Neumann expansion of 1/r₁₂ have problems of numerical stability that are analyzed in detail using computer‐assisted algebraic techniques. Some of these problems can be eliminated by use of procedures known in this field 40 years ago but seemingly forgotten now. Others can be removed by use of a formulation suitable for small values of the STO screening parameter. A recent attempt at such a formulation is corrected and extended in a way permitting its practical use. The main functions encountered in the integrations over the ellipsoidal coordinate of the range 1 … ∞ are Bessel functions or generalizations thereof, as pointed out here for the first time. This fact is used to motivate the derivation of recurrence relations additional to those previously known. Novel techniques were devised for using these recurrence relations, thereby providing new ways of calculating the quantities that enter the ellipsoidal expansion. The convergence rate of this expansion and the numerical characteristics of several computational strategies are reported in enough detail to identify the ranges where various schemes can be used. This information shows that recent discussions of the “convergence characteristics of [the] ellipsoidal coordinate expansion” are in fact not that, but are instead discussions of an inability to make accurate calculations of the individual terms of the expansion. It is also seen that the parameter range suitable for use of Kotani's well‐known recursive scheme is more limited than seems generally believed. The procedures discussed in this work are capable of yielding accurate two‐center electron repulsion integrals by the ellipsoidal expansion method for all reasonable STO screening parameters, and have been implemented in illustrative public‐domain computer programs. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Evaluation of the auxiliary function $${G_{-ns}^q}$$ using basic overlap integrals over Slater type orbitals

This paper presents a computationally efficient formula in terms of basic overlap integrals over Slater type orbitals (STOs) for the evaluation of auxiliary function which plays a central role in calculations of multicenter molecular integrals. The basic overlap integrals are calculated with the help of recurrence relations. The resulting simple analytical formula for the auxiliary function is completely general for p _a ≤ 1.2 and arbitrary values of parameters p and pt. The efficiency of calculation of auxiliary function is compared with other method.     

Analytical evaluation for two-center nuclear attraction integrals over slater type orbitals by using Fourier transform method

In this study, we shall suggest analytical expressions for two-center nuclear attraction integrals over STO’s with a one-center charge distribution by using Fourier transform method. The derivation is based on partial-fraction decompositions and Taylor expansions of rational functions. Analytical expressions obtained by this method are expressed in terms of Gegenbauer, and binomial coefficients and linear combinations of STO’s. Finally, it is relatively easy to express the Fourier integral representations of two-center nuclear attraction integrals with a one-center charge distribution mentioned above as finite and infinite of series of STO’s and irregular solid harmonics which may be considered to be limiting cases of STO’s.     

Evaluation of Hylleraas-CI atomic integrals by integration over the coordinates of one electron. I. Three-electron integrals

A method to evaluate the nonrelativistic electron-repulsion, nuclear attraction and kinetic energy three-electron integrals over Slater orbitals appearing in Hylleraas-CI (Hy-CI) electron structure calculations on atoms is shown. It consists on the direct integration over the interelectronic coordinate r _ij and the sucessive integration over the coordinates of one of the electrons. All the integrals are expressed as linear combinations of basic two-electron integrals. These last are solved in terms of auxiliary two-electron integrals which are easy to compute and have high accuracy. The use of auxiliary three-electron ones is avoided, with great saving of storage memory. Therefore this method can be used for Hy-CI calculations on atoms with number of electrons N ≥ 5. It has been possible to calculate the kinetic energy also in terms of basic two-electron integrals by using the Hamiltonian in Hylleraas coordinates, for this purpose some mathematical aspects like derivatives of the spherical harmonics with respect to the polar angles and recursion relations are treated and some new relations are given.     

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