Recurrence relations for the evaluation of electron repulsion integrals over spherical Gaussian functions

Recurrence relations are derived for the evaluation of two-electron repulsion integrals (ERIs) over Hermite and spherical Gaussian functions. Through such relations, a generic ERI or ERI derivative may be reduced to “basic” integrals, i.e., true and auxiliary integrals involving only zero angular momentum functions. Extensive use is made of differential operators, in particular, of the spherical tensor gradient ??(?). Spherical Gaussians, being nonseparable in the x, y, and z coordinates, were not included in previous formulations. The advantages of using spherical Gaussians instead of Cartesian or Hermite Gaussians are briefly discussed. © 1993 John Wiley & Sons, Inc.     

New recurrence relations for the rapid evaluation of electron repulsion integrals based on the accompanying coordinate expansion formula

We present an algorithm for the rapid computation of electron repulsion integrals (ERIs) over Gaussian basis functions based on the accompanying coordinate expansion (ACE) formula. The present algorithm uses equations termed angular momentum reduced expressions and introduces two types of recurrence relations to ACE formulas. Numerical efficiencies are assessed for (p pmid R:p p) and (sp spmid R:sp sp) ERIs by using the floating-point operation count. The algorithm is suitable for calculating ERIs for the same exponents but different angular momentum functions, such as L shells and derivatives of ERIs. The present algorithm is also capable of calculating ERIs with highly contracted Gaussian basis functions.     

A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians

Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals. For two- and three-center two-electron integrals, the solid-harmonic Hermite scheme is particularly efficient, significantly reducing the cost relative to the Cartesian scheme.     

Analytical calculations of molecular integrals for multielectron R-matrix methods

Closed-form analytical expressions for one- and two-electron integrals between Cartesian Gaussians over a finite spherical region of space are developed for use in ab initio molecular scattering calculations. In contrast with some previous approaches, the necessary integrals are formulated solely in terms of finite summations involving standard functions. The molecular integrals evaluated over the finite region of space are computed by subtracting the contributions outside the region from the integrals over all space. The latter integrals can be efficiently and accurately obtained from existing bound-state algorithms. Our approach incorporates molecular scattering calculations into current quantum chemistry programs and facilitates the unification of bound- and continuum-state calculations for both diatomic and polyatomic molecules. Multidimensional Monte Carlo numerical integrations validate the high accuracy of our closed form results for the two-electron integrals.     

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.     

Use of STOs in Hartree–Fock calculations: Error analysis and variance-minimized pseudospectral method

In this study it is demonstrated that STO (Slater-type orbital) basis sets are particularly well suited to pseudospectral Hartree–Fock calculations. The reduction of two-electron integrals, to ones that are (at worst) equivalent to a one-electron integral over three centers, eliminates the need for slowly convergent one-center expansions. This allows all integrals to be calculated quickly and accurately in either spherical or ellipsoidal coordinates. A new variance-minimized variant of the pseudospectral method is derived and applied to a number of small closed-shell molecules. The performance of the algorithm is assessed relative to purely spectral calculations employing STO and GTO (Gaussian-type orbital) basis sets. The pseudospectral operator is used to assess the errors contained in solutions found by the purely spectral method. The suitability of a number of different de-aliasing set types is also examined. Orthogonal sets of hydrogen-like eigenfunctions were found to be optimal. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1537–1548, 1999     

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.     

Molecular integrals for Gaussian and exponential‐type functions: Shift operators

Basis functions with arbitrary quantum numbers can be attained from those with the lowest numbers by applying shift operators. We derive the general expressions and the recurrence relations of these operators for Cartesian basis sets with Gaussian and exponential radial factors. In correspondence, the expressions of molecular integrals involving functions with arbitrary quantum numbers can be obtained by applying these operators on the integrals with the lowest quantum numbers. Since the original form of the shift operators is not appropriate to deal with integrals, we give their representation in terms of derivatives with respect to the parameters on which these integrals explicitly depend. Moreover, we translate the recurrence relations to the new representation and, finally, we analyze the general expressions ot the molecular integrals. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 137–145, 2000     

Recurrence formulas for molecular integrals over Laguerre Gaussian-type functions

Recurrence formulas for overlap, nuclear attraction, and electron-repulsion integrals over Laguerre Gaussian-type functions are presented. They have been derived using compact recurrence relations for homogeneous solid spherical harmonic operators but are rather lengthy as compared to those over Cartesian Gaussian-type functions. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 273–279, 1998     

Calculation of overlap integrals over Slater-type orbitals using translational and rotational transformations for spherical harmonics

Using translation and rotation formulas for spherical harmonics the finite sums through the basic overlap integrals and spherical harmonics are derived for the arbitrary overlap integrals over Slater-type orbitals (STOs). The recurrence relations for the evaluation of basic overlap integrals have been established recently [Guseinov II, Mamedov BA (1999) J Mol Struct (THEOCHEM) 465:1]. By the use of the derived expressions the overlap integrals can be calculated most efficiently and accurately, especially for large quantum numbers of STOs. Received: 2 May 2000 / Accepted: 31 May 2000 / Published online: 11 September 2000     

On the use of spatial symmetry in atomic-integral calculations: An efficient permutational approach

The minimal number of independent nonzero atomic integrals that occur over arbitrarily oriented basis orbitals of the form ?(r) · Y_lm(Ω) is theoretically derived. The corresponding method can be easily applied to any point group, including the molecular continuous groups C_∞v and D_∞h. On the basis of this (theoretical) lower bound, the efficiency of the permutational approach in generating sets of independent integrals is discussed. It is proved that lobe orbitals are always more efficient than the familiar Cartesian Gaussians, in the sense that GLO s provide the shortest integral lists. Moreover, it appears that the new axial GLO s often lead to a number of integrals, which is the theoretical lower bound previously defined. With AGLO s, the numbers of two-electron integrals to be computed, stored, and processed are divided by factors 2.9 (NH₃), 4.2 (C₅H₅), and 3.6 (C₆H₆) with reference to the corresponding CGTO s calculations. Remembering that in the permutational approach, atomic integrals are directly computed without any four-indice transformation, it appears that its utilization in connection with AGLO s provides one of the most powerful tools for treating symmetrical species.     

Application of second-order Møller-Plesset perturbation theory with resolution-of-identity approximation to periodic systems

Efficient periodic boundary condition (PBC) calculations by the second-order M?ller-Plesset perturbation (MP2) method based on crystal orbital formalism are developed by introducing the resolution-of-identity (RI) approximation of four-center two-electron repulsion integrals (ERIs). The formulation and implementation of the PBC RI-MP2 method are presented. In this method, the mixed auxiliary basis functions of the combination of Poisson and Gaussian type functions are used to circumvent the slow convergence of the lattice sum of the long-range ERIs. Test calculations of one-dimensional periodic trans-polyacetylene show that the PBC RI-MP2 method greatly reduces the computational times as well as memory and disk sizes, without the loss of accuracy, compared to the conventional PBC MP2 method.     

Derivative of electron repulsion integral using accompanying coordinate expansion and transferred recurrence relation method for long contraction and high angular momentum

In this study, an early‐working algorithm is designed to evaluate derivatives of electron repulsion integrals (DERIs) for heavy‐element systems. The algorithm is constructed to extend the accompanying coordinate expansion and transferred recurrence relation (ACE‐TRR) method, which was developed for rapid evaluation of electron repulsion integrals (ERIs) in our previous article (M. Hayami, J. Seino, and H. Nakai, J. Chem. Phys. 2015, 142, 204110). The algorithm was formulated using the Gaussian derivative rule to decompose a DERI of two ERIs with the same sets of exponents, different sets of contraction coefficients, and different angular momenta. The algorithms designed for segmented and general contraction basis sets are presented as well. Numerical assessments of the central processing unit time of gradients for molecules were conducted to demonstrate the high efficiency of the ACE‐TRR method for systems containing heavy elements. These heavy elements may include a metal complex and metal clusters, whose basis sets contain functions with long contractions and high angular momenta.     

Comment on “Analytical evaluation for two-center nuclear attraction integrals over Slater type orbitals by using Fourier transform method” by S. ?zcan and E. ?ztekin (J. Math. Chem. DOI 10.1007/s10910-008-9398-z)

S. ?zcan and E. ?ztekin, (J. Math. Chem. doi:) published formulas for evaluating the two-center nuclear attraction integrals over Slater type orbitals. It is shown that the analytical relations for these integrals through the expansion coefficients of the electron charge density for the one-center case and the overlap integrals presented in Sect. 3 of this work can easily be derived by means of a simple algebra from the formulas published in our papers (I.I. Guseinov, J Mol Struct (Theochem) 417:117, 1997; J Math Chem 42:415, 2007 and B.A. Mamedov, Chin J Chem 22:545, 2004). It should be noted that the formulas of overlap integrals presented by E. ?ztekin et al., in previous paper (E. ?ztekin, M. Yavuz, Ş. Atalay, J Mol Struct (Theochem) 544:69, 2001) for the calculation of two-center nuclear attraction integrals also are obtained from our papers (see Comment: I.I. Guseinov, J Mol Struct (Theochem) 638:235, 2003).     

Cubic-grid Gaussian basis sets for electron scattering calculations. I. Definition and construction

We propose a new type of Gaussian basis sets for use in calculations of electron scattering by molecules. Instead of locating the basis-set functions on the atomic centers of the target molecule, we place primitive s-type Gaussians at the positions of a cubic lattice with a regular grid. The grid and the Gaussian exponent are fixed so as to give the best representation of the plane-wave function. Plane-wave functions and Green functions obtained by means of the cubic-grid basis set are tested graphically against exact functions and functions expressed by means of a conventional Gaussian basis set. © 1995 John Wiley & Sons, Inc.     

Molecular symmetry and ab initio calculations. II. Symmetry-matrix and symmetry-supermatrix in the Dirac-Fock method

The concepts of symmetry-matrix and symmetry-supermatrix introduced in article I [J. Comput. Chem., 10, 957 (1989)] can be generalized to the Dirac-Fock method. By using the semidirect product decomposition of O_h and the linear vector space theory, the irreducible representation basis of O_h for any molecular system (O_h or its subgroups) can be deduced analytically in the nonorthonormal Cartesian Gaussian basis. This method is extended to discuss the double-valued representations of O_h* in the complex Cartesian Gaussian spinor basis. In the double-valued irreducible representation basis of D₂*, the matrix of kinetic operator c(OVERLINE)σ(/OVERLINE)·(OVERLINE)p(/OVELINE) in the Dirac-Fock equation can be reduced into a real symmetric and can be grouped into classes under the operations in D_3d. Therefore, the symmetry-matrix and symmetry-supermatrix can also be used in the Dirac-Fock method to reduce the storage of two electron integrals and calculations of Fock matrix during iterations by a factor of ca. g² (g is the order of the molecular symmetry group). In addition, a method to deal with the nonorthonormal space is presented. © 1996 by John Wiley & Sons, Inc.     

Integral evaluation algorithms and their implementation

Three different algorithms for the calculation of many center electron-repulsion integrals are discussed, all of which are considered to be economic in terms of the number of arithmetic operations. The common features of the algorithms are as follows: Cartesian Gaussian functions are used, integrals are calculated by blocks (a block being defined as the set of integrals obtainable from four given exponents on four given centers), and functions may be adopted to R(3). Adaption to molecular point group symmetry is not considered. Tables are given showing the minimum number of operations for a selection of block types allowing one to identify the theoretically most economic, and the corresponding salient features. Comments concerning the computer implementations are also given both on sealar and vector processors. In particular, the Cyber 205 is considered, a vector processor on which we have implemented what we believe to be the most efficient algorithm.     

Auxiliary functions for molecular integrals with Slater‐type orbitals. II. Gauss transform methods

The Gauss transform of Slater‐type orbitals is used to express several types of molecular integrals involving these functions in terms of simple auxiliary functions. After reviewing this transform and the way it can be combined with the shift operator technique, a master formula for overlap integrals is derived and used to obtain multipolar moments associated to fragments of two‐center distributions and overlaps of derivatives of Slater functions. Moreover, it is proved that integrals involving two‐center distributions and irregular harmonics placed at arbitrary points (which determine the electrostatic potential, field and field gradient, as well as higher order derivatives of the potential) can be expressed in terms of auxiliary functions of the same type as those appearing in the overlap. The recurrence relations and series expansions of these functions are thoroughly studied, and algorithms for their calculation are presented. The usefulness and efficiency of this procedure are tested by developing two independent codes: one for the derivatives of the overlap integrals with respect to the centers of the functions, and another for derivatives of the potential (electrostatic field, field gradient, and so forth) at arbitrary points. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Two-electron integral evaluation for uncontracted geometrical-type Gaussian functions

A new algorithm for efficient evaluation of two-electron repulsion integrals (ERIs) using uncontracted geometrical-type Gaussian basis functions is presented. Integrals are evaluated by the Habitz and Clementi method. The use of uncontracted geometrical basis sets allows grouping of basis functions into shells (s, sp, spd, or spdf) and processing of integrals in blocks (shell quartets). By utilizing information common to a block of integrals, this method achieves high efficiency. This technique has been incorporated into the KGNMOL molecular interaction program. Representative timings for a number of molecules with different basis sets are presented. The new code is found to be significantly faster than the previous program. For ERIs involving only s and p functions, the new algorithm is a factor of two faster than previously. The new program is also found to be competitive when compared with other standard molecular packages, such as HONDO-8 and Gaussian 86.     

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