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     

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     

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.     

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     

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     

Computation of overlap integrals over Slater-type orbitals using auxiliary functions

Analytical expressions through the binomial coefficients and recursive relations are derived for the expansion coefficients of overlap integrals in terms of a product of well-known auxiliary functions A_k and B_k. These formulas are especially useful for the calculation of overlap integrals for large quantum numbers. Accuracy of the computer results is satisfactory for the values of quantum numbers up to 50 and for the arbitrary values of screening constants of atomic orbitals and internuclear distances. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 199–204, 1998     

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     

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.     

Auxiliary functions for Slater molecular integrals

Summary Some types of recurrence relations are modified to overcome the cases in which their conventional application is unstable in both the forward and backward directions. The original recurrence relations — connecting adjacent elements — are replaced by more general ones, where the non-adjacent elements are connected by coefficients obtained by new sets of relations derived from the original ones. This modification can be helpful for the calculation of the complicated molecular integrals with Slater Type Orbitals (STOs).As a simple test we prove that some auxiliary functions — previously evaluated by expensive expansions — appearing in two-center two-electron integrals can be thus calculated with very low computer cost and high accuracy.     

Libcint: An efficient general integral library for Gaussian basis functions

An efficient integral library Libcint was designed to automatically implement general integrals for Gaussian‐type scalar and spinor basis functions. The library is able to evaluate arbitrary integral expressions on top of p, r and σ operators with one‐electron overlap and nuclear attraction, two‐electron Coulomb and Gaunt operators for segmented contracted and/or generated contracted basis in Cartesian, spherical or spinor form. Using a symbolic algebra tool, new integrals are derived and translated to C code programmatically. The generated integrals can be used in various types of molecular properties. To demonstrate the capability of the integral library, we computed the analytical gradients and NMR shielding constants at both nonrelativistic and 4‐component relativistic Hartree–Fock level in this work. Due to the use of kinetically balanced basis and gauge including atomic orbitals, the relativistic analytical gradients and shielding constants requires the integral library to handle the fifth‐order electron repulsion integral derivatives. The generality of the integral library is achieved without losing efficiency. On the modern multi‐CPU platform, Libcint can easily reach the overall throughput being many times of the I/O bandwidth. On a 20‐core node, we are able to achieve an average output 8.3 GB/s for C₆₀ molecule with cc‐pVTZ basis. © 2015 Wiley Periodicals, Inc.     

Improved partition–expansion of two‐center distributions involving slater functions

The calculation of the electronic structure of large systems is facilitated by the substitution of the two‐center distributions by their projections on auxiliary basis sets of one‐center functions. An alternative is the partition–expansion method in which one first decides what part of the distribution is assigned to each center, and next expands each part in spherical harmonics times radial factors. The method is exact, requires neither auxiliary basis sets nor projections, and can be applied to Gaussian and Slater basis sets. Two improvements in the partition–expansion method for Slater functions are reported: general expressions valid for arbitrary quantum numbers are derived and the efficiency of the procedure is increased giving analytical solutions to integrals previously computed by numerical quadrature. The efficiency of the new version is assessed in several molecules and the advantages over the projection methods are pointed out. © 2013 Wiley Periodicals, Inc.     

Generalized spherical functions in the theory of many-electron atoms

A new method for finding non-relativistic and relativistic wave-functions of an electron moving in the field of a nuclear charge in the jj coupling scheme is proposed. It is based on the usage of generalized spherical functions. The mathematical apparatus necessary to find the expressions for matrix elements of the non-relativistic and relativistic energy or electron transition operators is developed. The formulas obtained for these matrix elements are more convenient than those usually used in jj coupling scheme; only their radial integrals and some phase multipliers depend on orbital quantum numbers.     

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     

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.     

A general formulation for the efficient evaluation of n-electron integrals over products of Gaussian charge distributions with Gaussian geminal functions

In this work, we present a general formulation for the evaluation of many-electron integrals which arise when traditional one particle expansions are augmented with explicitly correlated Gaussian geminal functions. The integrand is expressed as a product of charge distributions, one for each electron, multiplied by one or more Gaussian geminal factors. Our formulation begins by focusing on the quadratic form that arises in the general n-electron integral. Using the Rys polynomial method for the evaluation of potential energy integrals, we derive a general formula for the evaluation of any n-electron integral. This general expression contains four parameters ω, θ, v, and h, which can be evaluated by an examination of the general quadratic form. Our analysis contains general expressions for any n-electron integral over s-type functions as well as the recursion needed to build up arbitrary angular momentum. The general recursion relation requires at most n + 1 terms for any n-electron integral. To illustrate the general method, we develop explicit expressions for the evaluation of two, three, and four particle electron repulsion integrals as well as two and three particle overlap and nuclear attraction integrals. We conclude our exposition with a discussion of a preliminary computational implementation as well as general computational requirements. Implementation on parallel computers is briefly discussed.     

Application of graphical methods of spin algebras to limited CI approaches. I. Closed shell case

The time independent diagrammatic technique based on the mathematical methods of quantum electrodynamics (second quantization, Wick's theorem, Feynman-like diagrams) is combined with graphical techniques of spin algebras to derive general expressions for the matrix elements of spin independent one- and two-particle operators between spin symmetry adapted ground, mono- and bi-excited configurations of a closed shell system. Two coupling schemes are considered for bi-excited states and their relative merits are discussed. Finally, the results are used to derive compact expressions for the coupling coefficients of the direct configuration interaction from molecular integrals (CIMI ) method.     

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.     

A simple O(4,2) approximation of hydrogenic coulomb integrals

Repulsion integrals for electrons with hydrogenic wavefunctions are expressed in terms of operations of the Lie algebra of the group O(4,2) × O(4,2). The results is used to obtain approximate expressions for Coulomb integrals as functions of hydrogenic quantum numbers     

Gaussian matrix elements of the operator X: Reduction to confluent hypergeometric functions

A new method has been proposed for calculating the Gaussian matrix elements of the interaction operator represented by an arbitrary degree of interelectron distance X. The method is based on the expansion of two-electron integrals as the sum of one-electron integrals which in turn admit compact operator representation in terms of confluent hypergeometric functions. The generating differential operator has been shown to be related to the modified Hermitian polynomials. The standard structure of the special functions encountered in this approach is useful in studying the analytical behavior of the integrals and makes it possible to obtain for these integrals recurrence relations, direct algebraic expressions in the forms of finite sum of confluent hypergeometric functions, integral representations, and asymptotic properties. Unlike the usual methods based on integral transformation of the interaction operator, the proposed approach has a wider field of application, and in addition, leads to compact and convenient analytical expressions. The idea of using differential properties of integrals to simplify the integrand structure gives the proposed approach a certain resemblance to that suggested by Boys but not developed in detail in his pioneer work.