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     

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     

Efficiency of the algorithms for the calculation of Slater molecular integrals in polyatomic molecules

The performances of the algorithms employed in a previously reported program for the calculation of integrals with Slater-type orbitals are examined. The integrals are classified in types and the efficiency (in terms of the ratio accuracy/cost) of the algorithm selected for each type is analyzed. These algorithms yield all the one- and two-center integrals (both one- and two-electron) with an accuracy of at least 12 decimal places and an average computational time of very few microseconds per integral. The algorithms for three- and four-center electron repulsion integrals, based on the discrete Gauss transform, have a computational cost that depends on the local symmetry of the molecule and the accuracy of the integrals, standard efficiency being in the range of eight decimal places in hundreds of microseconds.     

Direct quantum chemical integral evaluation

Performing quantum chemical integral evaluation directly, without recursion and without direct coupling of angular momenta according to the rotation group is analyzed. The rotation group limits the structure of these closed‐form expressions. The result of all cross differentiation is a rotational invariant. Closed‐form expressions are obtained for the general three‐ and four‐center Gaussian integral. The solid harmonic addition formula can be used to express these integrals as sums of products of an exponent‐independent (angular) factor and a molecular‐orientation‐independent (exponential) factor in a variety of ways. The results are products of two such factors summed over the set of distinct, relevant polynomials of the exponents. The coefficients of these polynomials, angular factors, are complicated but common to all n‐center matrix elements and independent of any type of contraction. Derivatives must be obtained using the product rule. An implementation in the Solid Spherical Harmonic Gaussian (SSHG) computer code is outlined and preliminary comparison is made. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 373–383, 2001     

Reference program for molecular calculations with Slater-type orbitals

A program for computing all the integrals appearing in molecular calculation with Slater-type orbitals is reported. The program is mainly intended as a reference for testing and comparing other algorithms and techniques. An analysis of the performance of the program is presented, paying special attention to the computational cost and the accuracy of the results. Results are also compared with others obtained with Gaussian basis sets of similar quality. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1284–1293, 1998     

Evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals with integer and noninteger principal quantum numbers

A general formula has been established for the expansion of the product of two normalized associated Legendre functions centered on the nuclei a and b. This formula has been utilized for the evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals (STOs) with integer and noninteger principal quantum numbers. The formulas given in this study for the evaluation of two‐center overlap and nuclear attraction integrals show good rate of convergence and great numerical stability under wide range of quantum numbers, orbital exponents, and internuclear distances. © 2001 Wiley Periodicals, Inc. Int J Quantum Chem, 2001     

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     

Auxiliary functions for molecular integrals with Slater‐type orbitals. I. Translation methods

Many types of molecular integrals involving Slater functions can be expressed, with the ζ‐function method in terms of sets of one‐dimensional auxiliary integrals whose integrands contain two‐range functions. After reviewing the properties of these functions (including recurrence relations, derivatives, integral representations, and series expansions), we carry out a detailed study of the auxiliary integrals aimed to facilitate both the formal and computational applications of the ζ‐function method. The usefulness of this study in formal applications is illustrated with an example. The high performance in numerical applications is proved by the development of a very efficient program for the calculation of two‐center integrals with Slater functions corresponding to electrostatic potential, electric field, and electric field gradient. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A full‐pivoting algorithm for the Cholesky decomposition of two‐electron repulsion and spin‐orbit coupling integrals

A significant reduction in the computational effort for the evaluation of the electronic repulsion integrals (ERI) in ab initio quantum chemistry calculations is obtained by using Cholesky decomposition (CD), a numerical procedure that can remove the zero or small eigenvalues of the ERI positive (semi)definite matrix, while avoiding the calculation of the entire matrix. Conversely, due to its antisymmetric character, CD cannot be directly applied to the matrix representation of the spatial part of the two‐electron spin‐orbit coupling (2e‐SOC) integrals. Here, we present a computational strategy to achieve a Cholesky representation of the spatial part of the 2e‐SOC integrals, and propose a new efficient CD algorithm for both ERI and 2e‐SOC integrals. The proposed algorithm differs from previous CD implementations by the extensive use of a full‐pivoting design, which allows a univocal definition of the Cholesky basis, once the CD δ threshold is made explicit. We show that is the upper limit for the errors affecting the reconstructed 2e‐SOC integrals. The proposed strategy was implemented in the ab initio program Computational Emulator of Rare Earth Systems (CERES), and tested for computational performance on both the ERI and 2e‐SOC integrals evaluation. © 2017 Wiley Periodicals, Inc.     

Symbolic Calculation of Two‐Center Overlap Integrals Over Slater‐Type Orbitals

Two‐center overlap integrals over Slater type orbitals (STOs) have been expressed in terms of the well‐known Mulliken's integrals B_n(pt) using Rodrigues's formula for normalized associated Legendre functions. A computer program is written in Mathematica 4.0 for the evaluation of two‐center overlap integrals over STOs. Using this computer program, symbolic tables are presented for two‐center overlap integrals up to quantum numbers 1 ≤ n,n′ ≤ 3, 0 ≤ l,l′ ≤ 2, ?2 ≤ m,m′ ≤ 2. Numerical results of this work, for some quantum sets, have also been compared with prior literature and best agreement achieved with recent works of Barnett while some discrepancies were obtained with works of Öztekin et al. and Guseinov et al.     

Three‐center Coulomb repulsion integrals with Slater functions

The translation method for products of two Slater functions is improved and combined with the short–long range separation in order to develop a robust and efficient algorithm for the evaluation of three‐center Coulomb integrals with Slater functions. Several tests are carried out showing that the algorithm reported here yields integrals with an absolute error below 10^?12 hartree and a computational cost of a few microseconds per integral. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Atomic and molecular cubature grids

This article presents cubature grids of the Gaussian type that are adapted for the purpose of wave function calculation on atoms and molecules. The problems of the singularity at the nucleus, the derivatives in the kinetic energy, and the presence of two‐electron integrals are shown to be resolved. Each grid has a definite degree of accuracy so that it reproduces the exact values of all the integrals in a defined class. Seventh‐degree accuracy can be obtained from a grid of 143 nodes. The grids are applied, as simple illustrations, to well‐known self‐consistent field (SCF) calculations on helium. Grids for homonuclear diatomics are also discussed and an illustrative application given to a homonuclear diatomic molecule. A comparison between a molecular grid and the union of two unmodified atomic grids shows that overlaps and distortions in weights can occur to the extent that this is not practical. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

New translation method for STOs and its application to calculation of two-center two-electron integrals

The new translation method for Slater-type orbitals (STOs) previously tested in the case of the overlap integral is extended to the calculation of two-center two-electron molecular integrals. The method is based on the exact translation of the regular solid harmonic part of the orbital followed by the series expansion of the residual spherical part in powers of the radial variable. Fair uniform convergence and stability under wide changes in molecular parameters are obtained for all studied two-center hybrid, Coulomb, and exchange repulsion integrals. Ten-digit accuracy in the final numerical results is achieved through multiple precision arithmetic calculation of common angular coefficients and Gaussian numerical integration of some of the analytical formulas resulting for the radial integrals. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 91–100, 2000     

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     

Master formulas for two‐ and three‐center one‐electron integrals involving Cartesian GTO,STO, and BTO

As a first application of the shift operators method we derive master formulas for the two‐ and three‐center one‐electron integrals involving Gaussians, Slater, and Bessel basis functions. All these formulas have a common structure consisting in linear combinations of polynomials of differences of nuclear coordinates. Whereas the polynomials are independent of the type (GTO, BTO, or STO) of basis functions, the coefficients depend on both the class of integral (overlap, kinetic energy, nuclear attraction) and the type of basis functions. We present the general expression of polynomials and coefficients as well as the recurrence relations for both the polynomials and the whole integrals. Finally, we remark on the formal and computational advantages of this approach. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 83–93, 2000     

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     

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     

Nuclear attraction and electron interaction integrals of exponentially decaying functions and the Poisson equation

Summary The technique proposed by O-Ohata and Ruedenberg (J Math Phys 7:547 (1966)) and by Silver and Ruedenberg (J Chem Phys 49:4306 (1968)) of computing nuclear attraction and electron interaction integrals by solving an inhomogeneous Laplace equation can also be applied ifB functions (Filter E, Steinborn EO (1978) Phys Rev A 18:1) are used as basis functions in atomic and molecular calculations. It is shown that because of the remarkable mathematical properties ofB functions the derivation of compact explicit expressions for the multicenter integrals mentioned above is particularly simple. These results are also of interest in the context of other exponentially decaying functions, since all the other commonly occurring exponentially decaying functions as, for instance, Slater functions or bound state hydrogen eigenfunctions can be expressed as simple linear combinations ofB functions. Consequently, their multicenter integrals can also be expressed in terms of multicenter integrals ofB functions.Dedicated to Prof. Klaus Ruedenberg on the occasion of his 70th birthday     

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     

Correspondence between GTO and STO molecular basis sets

We present a method for the characterization of the distance between two spaces: one generated by a Gaussian basis set, and another by a Slater basis set. The method is an extension of one previously developed for atoms that has been modified to cover molecular problems. The current version enables us to obtain Slater basis sets capable of reproducing the results (multielectronic wave functions and orbitals) obtained with Gaussian basis sets. The interest of this result arises from the fact that we will be able to profit from the effort invested in the optimization of high‐quality Gaussian basis sets. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1655–1665, 2001