Efficient evaluation of Coulomb integrals in a mixed Gaussian and plane‐wave basis

An efficient way of calculation is presented for matrix elements between two plane waves interacting with a molecular Coulombic field. In concurrence with the absolute value of the momentum transfer vector, K = k₁ ? k ₂ , the most effective method of calculation is selected. The case of K = 0 requires special treatment. For 0 < |K| ≤ 0.3, it is profitable to evaluate the integrals by means of the multipole expansion, and for |K| > 0.3 the density fitting can be applied. For the large |K| the electronic part of the integral is much smaller than the nuclear part and the integral may be approximated by the nuclear contribution only. Some examples for testing the accuracy and time saved are presented. The primary purpose of this paper is to accelerate electron scattering calculations, but it also may be profitable for the electronic structure theory in attempts to use mixed Gaussian and plane‐wave basis sets. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density fitting for derivatives of Coulomb integrals in ab initio calculations using mixed Gaussian and plane‐wave basis

This article presents an efficient way for evaluation of Coulomb integrals of the type (k₁ (1)k₂ (1)|g₁ (2)g₂ (2)) and their derivatives with respect to nuclear coordinates by means of density fitting. Symbols k₁ andk₂ stand for plane‐wave functions and g₁ and g₂ for gaussians. The study was undertaken with the objective to accelerate electronic structure and electron scattering calculations in which a mixed Gaussian and plane‐wave basis set is used. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Efficient evaluation of one-center three-electron Gaussian integrals

 An algorithm is presented for the efficient evaluation of two types of one-center three-electron Gaussian integrals. These integrals are required to avoid the resolution-of-identity (RI) approximation in explicitly correlated linear R12 methods. Without the RI approximation, it is possible to enforce rigorously the strong orthogonality of the second-order M?ller–Plesset R12 ansatz. A test calculation is performed using atomic Gaussian-type orbitals of the neon atom. Received: 21 November 2000 / Accepted: 6 April 2001 / Published online: 9 August 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     

Half-numerical evaluation of pseudopotential integrals

A half-numeric algorithm for the evaluation of effective core potential integrals over Cartesian Gaussian functions is described. Local and semilocal integrals are separated into two-dimensional angular and one-dimensional radial integrals. The angular integrals are evaluated analytically using a general approach that has no limitation for the l-quantum number. The radial integrals are calculated by an adaptive one-dimensional numerical quadrature. For the semilocal radial part a pretabulation scheme is used. This pretabulation simplifies the handling of radial integrals, makes their calculation much faster, and allows their easy reuse for different integrals within a given shell combination. The implementation of this new algorithm is described and its performance is analyzed.     

Rapid evaluation of molecular integrals with ETOs

In this article, we discuss a way in which the theory of hyperspherical harmonics may be used for rapid evaluation of difficult molecular integrals when exponential‐type orbitals (ETOs) are used as a basis. One of us (J.E.A.) has implemented the method, and programs are available for general use. As a byproduct of this work, we are also able to evaluate generalized scattering factors for ETOs which allow first‐order density matrices to be measured experimentally using high‐quality X‐ray diffraction data. © 2015 Wiley Periodicals, Inc.     

Computation of some new two-electron Gaussian integrals

Summary The evaluation of a new form of two-electron integrals is required if the interelectronic distancer ₁₂ is used as a variable in then-electron functions of electron correlation methods. The McMurchie-Davidson algorithm for the generation of molecular integrals over Gaussian-type functions is ideally suited to this. The new Gaussian integrals are formed from Hermite integrals overr ₁₂ (rather than 1/r ₁₂) by standard techniques. The Hermite integrals overr ₁₂ itself are generated by a simple procedure with negligible computational effort. The key results are discussed in the context of general recursion formulas. On leave from: Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, W-4630 Bochum, Germany     

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.     

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     

Method for evaluation of density functional integrals in molecular calculations

Numerical methods for computing variationally optimized molecular orbitals within the Hartree–Fock approximation are augmented to include correlation functionals of the density in the energy and the numerical methods for carrying this out are described. The approach is applied explicitly to the Colle–Salvetti correlation energy functional. It is found that the gradient terms in the Colle–Salvetti functional present numerical problems associated with the low-density behavior, but also that they make a relatively small contribution to the correlation energy. In the three cases considered, HF, H₂O and N₂, it is found that the Colle–Salvetti correction considerably underestimates the correlation energies obtained in coupled-cluster theory.     

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.     

Gaussian expansions of orbital products for the evaluation of two-electron integrals

A method is described for evaluating multicenter integrals over contracted gaussian-type orbitals by the use of gaussian expansion of orbital products. The expansions are determined by the method of non-linear least squares with constraints. There is no restriction upon the symmetry of the orbital product and the method is applicable to all products arising from s, p and d-type orbitals. Results are given to indicate the accuracy of the method.     

Convergence accelerator approach for the evaluation of some three‐electron integrals containing explicit rij factors

A simplified analysis is presented for the evaluation of the three‐electron one‐center integrals of the form ∫rrrrrred r ₁d r ₂d r ₃, for the cases i, j, k, ≥−2, l=−2, m≥−1, n≥−1. These integrals arise in the calculation of lower bounds for energy levels and certain relativistic corrections to the energy when Hylleraas‐type basis sets are employed. Convergence accelerator techniques are employed to obtain a reasonable number of digits of precision, without excessive CPU requirements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 93–99, 1999     

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     

The reduced multiplication scheme of the Rys-Gauss quadrature for 1st order integral derivatives

Summary An implementation of the reduced multiplication scheme of the Rys-Gauss quadrature to compute the gradients of electron repulsion integrals is discussed. The study demonstrates that the Rys-Gauss quadrature is very suitable for efficient utilization of simplifications as offered by the direct computation of symmetry adapted gradients and the use of the translational invariance of the integrals. The introduction of the so-called intermediate products is also demonstrated to further reduce the floating point operation count. Two prescreening techniques based on the 2nd order density matrix in the basis of the uncontracted Gaussian functions is proposed and investigated in the paper. This investigation gives on hand that it is not necessary to employ the Cauchy-Schwarz inequality to achieve efficient prescreening. All the features mentioned above were demonstrated by their implementation into the gradient programalaska. The paper offers a theoretical and practical assessment of the modified Rys-Gauss quadrature in comparison with other methods and implementations and a detailed analysis of the behavior of the method as suggested above as a function of changes with respect to symmetry, basis set quality, molecular size, and prescreening threshold.     

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     

Gaussian and finite-element Coulomb method for the fast evaluation of Coulomb integrals

The authors propose a new linear-scaling method for the fast evaluation of Coulomb integrals with Gaussian basis functions called the Gaussian and finite-element Coulomb (GFC) method. In this method, the Coulomb potential is expanded in a basis of mixed Gaussian and finite-element auxiliary functions that express the core and smooth Coulomb potentials, respectively. Coulomb integrals can be evaluated by three-center one-electron overlap integrals among two Gaussian basis functions and one mixed auxiliary function. Thus, the computational cost and scaling for large molecules are drastically reduced. Several applications to molecular systems show that the GFC method is more efficient than the analytical integration approach that requires four-center two-electron repulsion integrals. The GFC method realizes a near linear scaling for both one-dimensional alanine alpha-helix chains and three-dimensional diamond pieces.     

Strategies for an efficient implementation of the Gauss-Bessel quadrature for the evaluation of multicenter integral over STFs

In a previous work, a new Gauss quadrature was introduced with a view to evaluate multicenter integrals over Slater-type functions efficiently. The complexity analysis of the new approach, carried out using the three-center nuclear integral as a case study, has shown that for low-order polynomials its efficiency is comparable to the SD. The latter was developed in connection with multi-center integrals evaluated by means of the Fourier transform of B functions. In this work we investigate the numerical properties of the Gauss-Bessel quadrature and devise strategies for an efficient implementation of the numerical algorithms for the evaluation of multi-center integrals in the framework of the Gaussian transform/Gauss-Bessel approach. The success of these strategies are essential to elaborate a fast and reliable algorithm for the evaluation of multi-center integrals over STFs.