Comment on “Symbolic Calculation of Two‐Center Overlap Integrals Over Slater‐type Orbitals”

Published by Gümü?. and Özdo?an formulas for the evaluation of two‐center overlap integrals (Gümü?, S.; Özdo?an, T. J. Chin. Chem. Soc. 2004, 51, 243) are critically analyzed. It is demonstrated that the formulas presented in this work are not original and they can easily be derived from the relationships contained in our papers (Guseinov, I. I. J. Phys. B 1970 , 3, 1399; Phys. Rev. A 1985 , 32, 1864; J. Mol Struct. (Theochem) 1995 , 336, 17) by changing the summation indices and application of a simple algebra. It should be noted that the symbolic results of overlap integrals between different combinations of quantum numbers given in Table 1 and 2 can also be obtained from the use of established in above mentioned our papers general formulas or presented in the literature relations for overlap integrals in terms of the products of molecular auxiliary functions A_n(p) and B_n(pt) (see, e.g., Lofthus, A. Mol. Phys. 1962 , 5, 105).     

Overlap integrals of B functions

Two different methods for the evaluation of overlap integrals of B functions with different scaling parameters are analyzed critically. The first method consists of an infinite series expansion in terms of overlap integrals with equal scaling parameters [14]. The second method consists of an integral representation for the overlap integral which has to be evaluated numerically. Bhattacharya and Dhabal [13] recommend the use of Gauss-Legendre quadrature for this purpose. However, we show that Gauss-Jacobi quadrature gives better results, in particular for larger quantum number. We also show that the convergence of the infinite series can be improved if suitable convergence accelerators are applied. Since an internal error analysis can be done quite easily in the case of an infinite series even if it is accelerated, whereas it is very costly in the case of Gauss quadratures, the infinite series is probably more efficient than the integral representation. Overlap integrals of all commonly occurring exponentially declining basis functions such as Slater-type functions, can be expressed by finite sums of overlap integrals of B functions, because these basis functions can be represented by linear combinations of B functions.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Multicenter integrals over polarization potential operators

Analytic expressions for multicenter integrals over the general one‐particle operator x^n′y^l′z^m′| r |^−k(1−exp(−αr²))ⁿ(n′, m′, l′, n≥0, k>2, α>0), employing Cartesian Gaussians, are presented. While until now only P. Schwerdtfeger and H. Silberbach (Phys Rev A 1988, 37, 2834) have succeeded in finding such expressions, using a Laplace transform, we shall show that one can also get them according to the method of L. E. McMurchie and E. R. Davidson (J Comp Phys 1978, 26, 218; J Comp Phys 1981, 44, 289). ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 403–416, 1999     

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     

Correlation energy extrapolation by intrinsic scaling. II. The water and the nitrogen molecule

The extrapolation method for determining benchmark quality full configuration-interaction energies described in preceding paper [L. Bytautas and K. Ruedenberg, J. Chem. Phys. 121, 10905 (2004)] is applied to the molecules H(2)O and N(2). As in the neon atom case, discussed in preceding paper [L. Bytautas and K. Ruedenberg, J. Chem. Phys. 121, 10905 (2004)] remarkably accurate scaling relations are found to exist between the correlation energy contributions from various excitation levels of the configuration-interaction approach, considered as functions of the size of the correlating orbital space. The method for extrapolating a sequence of smaller configuration-interaction calculations to the full configuration-interaction energy and for constructing compact accurate configuration-interaction wave functions is also found to be effective for these molecules. The results are compared with accurate ab initio methods, such as many-body perturbation theory, coupled-cluster theory, as well as with variational calculations wherever possible.     

Molecular symmetry and ab initio calculations: IV. Symmetry-matrix and symmetry-supermatrix in calculations of two-electron repulsion integrals

The product of two Gaussians having different centers is itself a one-center Gaussian, thus multicenter integrals with a Cartesian Gaussian basis can be reduced to one-center integrals. Recurrence relations for overlap integrals and electron repulsion integrals (ERIs) are derived at these centers. The calculations of overlap integrals and ERIs are carried out step by step from the highest symmetry case (one center) to required cases (different centers) by using the translation of Cartesian Gaussians. Full exploitation of symmetry in calculation processes can result in optimal use of these recurrence relations. Compared with the recently published algorithms, based on the recurrence relations derived by Obara and Saika [J. Chem. Phys., 84 , 3963 (1986)], the floating point operations (FLOPs) for ERI calculations (having four different centers) can be reduced by a factor of ca. 2. A significant extra saving in calculations and storage can be obtained if atoms, linear, or planar molecules are discussed. © 1997 John Wiley & Sons, Inc.     

Nonlinear transformation methods for accelerating the convergence of Coulomb integrals over exponential type functions

It is well known that in any ab initio molecular orbital (MO) calculation, the major task involves the computation of molecular integrals, among which the computation of Coulomb integrals are the most frequently encountered. As the molecular system gets larger, computation of these integrals becomes one of the most laborious and time consuming steps in molecular systems calculation. Improvement of the computational methods of molecular integrals would be indispensable to a further development in computational studies of large molecular systems. The atomic orbital basis functions chosen in the present work are Slater type functions. These functions can be expressed as finite linear combinations of B functions which are suitable to apply the Fourier transform method. The difficulties of the numerical evaluation of the analytic expressions of the integrals of interest arise mainly from the presence of highly oscillatory semi-infinite integrals. In this work, we present a generalized algorithm based on the nonlinear transformation of Sidi, for a precise and fast numerical evaluation of molecular integrals over Slater type functions and over B functions. Numerical results obtained for the three-center two-electron Coulomb and hybrid integrals over B functions and over Slater type functions. Comparisons with numerical results obtained using alternatives approaches and an existing code are listed.     

Evaluation of Multicenter Electronic Attraction, Electric Field and Electric Field Gradient Integrals with Screened and Nonscreened Coulomb Potentials over Integer and Noninteger n Slater Orbitals

By the use of complete orthonormal sets of -exponential-type orbitals, where ( = 1, 0, –1, –2,...) the multicenter electronic attraction (EA), electric field (EF) and electric field gradient (EFG) integrals of nonscreened and Yukawa-like screened Coulomb potentials are expressed through the two-center overlap integrals with the same screening constants and the auxiliary functions introduced in our previous paper (I.I. Guseinov, J. Phys. B, 3 (1970) 1399). The recurrence relations for auxiliary functions are useful for the calculation of multicenter EA, EF and EFG integrals for arbitrary integer and noninteger values of principal quantum numbers, screening constants, and location of slater-type orbitals. The convergence of the series is tested by calculating concrete cases.     

Computation of two‐center Coulomb integrals over Slater‐type orbitals using elliptical coordinates

A general analytic formula is obtained for the two‐center Coulomb integrals over Slater‐type orbitals in elliptical coordinates. Finite series expansions are used in the evaluation of the radial part of the integrals. The analytic formula is expressed in terms of a product of the well‐known auxiliary functions A_k(p) and B_k(p) and incomplete gamma functions. Recursive relations for the computer evaluation of these functions are given as well. The recursive relations are stable and our computer results are in good agreement with the benchmark values given in the literature. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Analytical Evaluation of Multicenter Multielectron Integrals of Central and Noncentral Interaction Potentials over Slater Orbitals Using Overlap Integrals and Auxiliary Functions

Using expansion formulas for central and noncentral interaction potentials (CIPs and NCIPs, respectively) in terms of Slater type orbitals (STOs) obtained by the author (I.I. Guseinov, J. Mol. Model., in press), the multicenter multielectron integrals of arbitrary interaction potentials (AIPs) are expressed through the products of overlap integrals with the same screening parameters and new auxiliary functions. For auxiliary functions, the analytic and recurrence relations are derived. The relationships obtained for multicenter multielectron integrals of AIDs are valid for the arbitrary quantum numbers, screening parameters and location of orbitals.     

Analytic first derivatives for explicitly correlated,multicenter, Gaussian geminals

Variational calculations utilizing the analytic gradient of explicitly correlated Gaussian molecular integrals are presented for the ground state of the hydrogen molecule. Preliminary results serve to motivate the need for general formulas for analytic first derivatives of molecular integrals involving multicenter, explicitly correlated Gaussian geminals with respect to Gaussian exponents and coordinates of the orbital centers. Explicit formulas for analytic first derivatives of Gaussian functions containing correlation factors of the form exp(-βr_ij²) are derived and discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 991–999, 1997     

Comment on “Calculation of Two-Center Nuclear Attraction Integrals over Integer and Noninteger n-Slater-Type Orbitals in Nonlined-Up Coordinate Systems”

Recently published formulas for the calculation of two-center nuclear attraction integrals (T. Özdoan, S. Gümü and M. Kara, J. Math. Chem. 33 (2003) 181) are critically analyzed. It is shown that the analytical relations presented in this work are not original and they can easily be derived by means of a simple algebra from the formulas for overlap integrals, their rotation coefficients and expansion of the product of two normalized associated Legendre functions in elliptical coordinates published in our papers (I.I. Guseinov, J. Phys. B 3 (1970) 1399; Phys. Rev. A 32 (1985) 1864; J. Mol. Struct.: Theochem 336 (1995) 17).     

Simplified ab-initio calculations for molecular systems

A method is described for reducing a large part of the arithmetic of exact ab-initio SCF molecularorbital calculations based on Slater-type-orbitals without noticeable loss of numerical accuracy. The procedure involves the transformation to Löwdin orthogonalized orbitals and then invoking the NDDO approximation. The three- and four-centre two-electron integrals required are estimated by a truncated Ruedenberg expansion. All one-electron integrals are evaluated exactly. No empirical parameters are employed. Numerical tests on CO, OF₂, O₃ and ONF show that the NDDO approximation is very accurate for Löwdin functions and that the Ruedenberg expansion is arithmetically satisfactory for the SCF MO calculations.     

The W algorithm and the D transformation for the numerical evaluation of three‐center nuclear attraction integrals

Three‐center nuclear attraction integrals over exponential‐type functions are required for ab initio molecular structure calculations and density functional theory (DFT). These integrals occur in many millions of terms, even for small molecules, and they require rapid and accurate numerical evaluation. The use of a basis set of B functions to represent atomic orbitals, combined with the Fourier transform method, led to the development of analytic expressions for these molecular integrals. Unfortunately, the numerical evaluation of the analytic expressions obtained turned out to be extremely difficult due to the presence of two‐dimensional integral representations, involving spherical Bessel integral functions. % The present work concerns the development of an extremely accurate and rapid algorithm for the numerical evaluation of these spherical Bessel integrals. This algorithm, which is based on the nonlinear D transformation and the W algorithm of Sidi, can be computed recursively, allowing the control of the degree of accuracy. Numerical analysis tests were performed to further improve the efficiency of our algorithm. The numerical results section demonstrates the efficiency of this new algorithm for the numerical evaluation of three‐center nuclear attraction integrals. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A method of computing multicenter integrals

A transformation analogous to the transform to the center-of-mass system is used to derive simple formulas for these integrals. A method of computation is used which involves expansion of the radial wave functions with respect to functions of an isotropic harmonic oscillator. It is shown that all the multicenter integrals in the case of an isotropic-oscillator function are expressed as finite sums, and that it is sufficient to use the known values for the unitary unimolecular groups SU₃ and SU₂ in order to compute these integrals.     

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     

Improved algorithm for the calculation of one-electron two-center integrals with STOs

In a previous article (J. Fernández Rico, R. López and G. Ramírez, J. Comp. Chem., 9 , 790 (1988)) we have proposed the calculation of molecular integrals involving STOs by means of some recurrence relations which use two sets ( h and H ) of overlap integrals (basic matrices). In the present paper, we derive explicit expressions of these integrals employing the two-range expansion of the 0s-function. This approach yields equations for the elements of the two basic matrices in terms of two further matrices, k (x,y) and i (x,y), and some auxiliary functions. Relations between the elements of these matrices and the functions are thoroughly explored and numerical tests are included for illustrating the behavior of the method.     

Evaluation of one-electron integrals for arbitrary operators V(r) over Cartesian Gaussians: Application to inverse-square distance and Yukawa operators

A general procedure is presented for generating one-electron integrals over any arbitrary potential operator that is a function of radial distance only. The procedure outlines that for a nucleus centered at point C integrals over Cartesian Gaussians can be written as linear combinations of 1-D integrals. These Cartesian Gaussian functions are expressed in a compact form involving easily computed auxiliary functions. It is well known that integrals over the Coulomb operator can be expressed in terms of F_n(T) integrals, where By means of a substitution for F_n(T) by other simple functions, algorithms that form integrals over an arbitrary function can be generated. Formation of such integrals is accomplished with minor editing of existing code based on the McMurchie–Davidson formalism. Further, the method is applied using the inverse-square distance and Yukawa potential operators V(r) over Cartesian Gaussian functions. Thus, the proposed methodology covers a large class of one-electron integrals necessary for theoretical studies of molecular systems by ab initio calculations. Finally, by virtue of the procedure's recursive nature it provides us with an efficient scheme of computing the proposed class of one-electron integrals. © 1993 John Wiley & Sons, Inc.     

Matrix elements for Ĵ2 and Ĵz operators over explicitly correlated Cartesian Gaussian functions

General formalism for evaluation of multiparticle integrals involving J?² and J?_z operators over explicitly correlated Cartesian Gaussian functions is presented. The integrals are expressed in terms of the general overlap integrals. An explicitly correlated Cartesian Gaussian function is a product of spherical orbital Gaussian functions, powers of the Cartesian coordinates of the particle, and exponential Gaussian factors, which depend on interparticular distances. This development is relevant to both adiabatic and nonadiabatic calculations of energy and properties of multiparticle systems. © 1995 John Wiley & Sons, Inc.     

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