A semiempirical model for the two-center repulsion integrals in the NDDO approximation

A self-consistent formalism is proposed for the two-center electron repulsion integrals in the NDDO approximation, based on their expansion in terms of multipole-multipole interactions and free from adjustable parameters.     

Some remarks on the semiempirical one-center electron repulsion integrals

The identification of the stage of ionization for various kinds of one-center electron repulsion integrals, occurring when nonbonding or lone-pair electrons are considered explicitly as well as π-electrons, is discussed for conjugated organic molecules containing heteroatoms N. It is concluded that the value for the negative ions should be used for (π_Cπ_C | π_Cπ_C) in all the states but for (π_Nπ_N | π_Nπ_N) only in the π-π states. In the n-π states, the appropriate value of (π_Nπ_N | π_Nπ_N) is that of the neutral atom if the molecule contains only one N atom. If more than one N atom is involved in the molecule, some weighted mean of the values for the negative ion and for the neutral atom should be used. The value for the neutral atom is most adequate for one-center repulsion integrals other than the (ππ | ππ) type in both the π-π and the n-π states. The method of determining these integrals is also discussed. It is concluded that they are to be determined from the consideration of appropriate electron-transfer reactions except for exchange integrals. The exchange integrals are shown to have to be determined from the Slater–Condon parameters derived from the analysis of the experimental atomic energy levels. Illustrative calculations are given for the lower singlet levels of the formaldehyde, pyrazine, pyridine, and the p-benzoquinone molecule. It is found that the calculated energies of the n-π transitions become much too low unless the (ππ | ππ) values of the heteroatoms in the molecule are chosen differently in the n-π states and in the π-π states as pointed out theoretically in this article.     

Methods for rapid evaluation of electron repulsion integrals in large-scale LCGO calculations

A method is proposed for efficient use of molecular symmetry in the evaluation of two-electron integrals. This provides a means of avoiding the recalculation of symmetry-redundant integrals, and of symmetry-blocking matrices and supermatrices without the usual time-consuming transformation procedures. Various methods for speeding up the calculation of integrals are also discussed. Integral calculation times are given for some representative molecules.     

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.     

Horizontal vectorization of electron repulsion integrals

We present an efficient implementation of the Obara–Saika algorithm for the computation of electron repulsion integrals that utilizes vector intrinsics to calculate several primitive integrals concurrently in a SIMD vector. Initial benchmarks display a 2–4 times speedup with AVX instructions over comparable scalar code, depending on the basis set. Speedup over scalar code is found to be sensitive to the level of contraction of the basis set, and is best for quartets when l_D = 0 or , which makes such a vectorization scheme particularly suitable for density fitting. The basic Obara–Saika algorithm, how it is vectorized, and the performance bottlenecks are analyzed and discussed. © 2016 Wiley Periodicals, Inc.     

Parallelization of three-center electron repulsion integrals

The parallelization of the three-center electron repulsion integrals arising from the variational fitting of the Coulomb potential is presented. A scheme for dynamical load balancing of the corresponding loop structure is discussed. The implementation in the density functional theory program deMon using the message passing interface is described. The efficiency of the parallelization is analyzed by selected benchmark calculations     

Evaluation of screened nuclear attraction and electron repulsion molecular integrals over Gaussian basis functions

Analytical solutions to the Yukawa-like screened Coulomb nuclear attraction and electron repulsion molecular basic integrals, as well as to the basic integrals required to compute the virial coefficient, over Gaussian basis functions, are derived and cast into a practical closed form, suitable to interface with modern codes for the calculation of molecular electronic structure. © 1997 John Wiley & Sons, Inc.     

One-Center coulomb,two-center hybrid,and two-center coulomb integrals over STP functions

Diatomic integral formulas involving Slater-transform-Preuss functions are obtained for the Coulomb and hybrid types. One-center Coulomb integrals are also considered. All results have finite sums containing a one-dimensional numerical integration.     

Efficient electronic integrals and their generalized derivatives for object oriented implementations of electronic structure calculations

For the new parallel implementation of electronic structure methods in ACES III (Lotrich et al., in preparation) the present state-of-the-art algorithms for the evaluation of electronic integrals and their generalized derivatives were implemented in new object oriented codes with attention paid to efficient execution on modern processors with a deep hierarchy of data storage including multiple caches and memory banks. Particular attention has been paid to define proper integral blocks as basic building objects. These objects are stand-alone units and are no longer tied to any specific software. They can hence be used by any quantum chemistry code without modification. The integral blocks can be called at any time and in any sequence during the execution of an electronic structure program. Evaluation efficiency of these integral objects has been carefully tested and it compares well with other fast integral programs in the community. Correctness of the objects has been demonstrated by several application runs on real systems using the ACES III program.     

Recursion formula for electron repulsion integrals over Hermite polynomials

A method for computing electron repulsion integrals over contracted Gaussian functions is described in which intermediate integrals over Hermite polynomials are generated by a “pre‐Hermite” recursion (PHR) step before the conversion to regular integrals. This greatly reduces the floating‐point operation counts inside the contraction loops, where only simple “scaling”‐type operations are required, making the method efficient for contracted Gaussians, particularly of high angular momentum. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Poisson equation for molecular exchange,hybrid and coulomb electron repulsion integrals

The various multicenter exchange, hybrid and Coulomb electron repulsion integrals that occur in molecular quantum mechanics are shown to satisfy a Poisson equation in which an overlap integral plays the role of a source distribution function. Two-, three-and four-center exchange integrals arise from four-center source functions; two- and three-center hybrid integrals arise from three-center distributions; and one- and two-center Coulomb integrals have two-center sources.     

Computation of electron repulsion integrals using the rys quadrature method

Following an earlier proposal to evaluate electron repulsion integrals over Gaussian basis functions by a numerical quadrature based on a set of orthogonal polynomials (Rys polynomials), a computational procedure is outlined for efficient evaluation of the two-dimensional integrals I_x, I_y, and I_z. Compact recurrence formulas for the integrals make the method particularly fitted to handle high-angular-momentum basis functions. The technique has been implemented in the HONDO molecular orbital program.     

Three-center expansion of electron repulsion integrals with linear combination of atomic electron distributions

A new systematic way of constructing auxiliary basis functions for approximating the evaluation of electron repulsion integrals is proposed and applied to SCF and MCSCF wavefunction calculations. In the approximation, the one-electron density is expanded in terms of a linear combination of atomic electron distributions (LCAD), and the four-center two-electron repulsion integrals are reduced to the three- and two-center quantities. This results in a high-accuracy approximation as well as a large reduction in disk storage and input/output requirement, proportional to N³ rather than N⁴, N being the number of basis functions. Numerical results indicate that the error from the present approximation decreases as the size of molecular basis functions increases and that the LCAD version of MCSCF calculations requires only a fractional amount of the CPU time required in the conventional procedure without loss of accuracy.     

Wavelets for electronic structure calculations

Molecular electronic structure calculations have a multi‐scale character through the presence of a set of singularities corresponding to atomic nuclei, and thus there exists a potential to improve the efficiency of these calculations using fast wavelet transform techniques. We report on the development of a one dimensional prototype benchmark problem of sufficient complexity to capture the features of 3‐D problems that are being solved today in quantum electronics calculations. Theoretical estimates of decay across scales and spatial distribution of wavelet coefficients for the solutions of the 1‐D and 3‐D problems are derived and verified experimentally. Equivalence in a multi‐resolution context of the solutions of the 1‐D prototype and the 3‐D problem is established. This revised version was published online in July 2006 with corrections to the Cover Date.     

On approximating electron repulsion integrals with linear combination of atomic-electron distributions

We recently investigated some features of atomic-electron distributions which are products of atomic orbitals without linear redundancies. Then, the atomic-electron distributions were employed as expansion functions of molecular charge distributions in the three-center expansion method for electron repulsion integrals. It turned out that the method is very accurate even after substantial reductions of required integrals. We shall illustrate the method and discuss its applicability to many electron theories. © 1996 John Wiley & Sons, Inc.