A vector processing algorithm of auxiliary integral evaluation for two-electron Gaussian integrals

A six-term auxiliary integral expression for the two-electron Gaussian integral is derived on the basis of the Chebyshev polynomial approximation instead of the seven-term Taylor expansion. This expression and the related recurrence formula enable us to perform a high-speed calculation on a vector processing computer.     

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     

Semi-asymptotic evaluation of certain three-center two-electron integrals

Two-electron repulsion integrals between a two-center charge distribution and a charge distribution about a third center, which do not appreciably interpenetrate, are shown to be given to useful accuracy by numerical differentiation of certain three-center one-electron integrals. This method also may be used to evaluate integrals of this type for which the Mulliken or Sklar approximations are inapplicable.Supported by Contract SD-102 with the Advanced Research Projects Agency.     

Analytic derivatives for the Cholesky representation of the two-electron integrals

We propose a formalism for calculating analytic derivatives of the electronic energy with respect to nuclear coordinates using Cholesky decomposition of the two-electron integrals. The formalism is derived by exploiting the equivalence of Cholesky decomposition and density fitting when a suitable auxiliary basis set is used for expanding atomic orbital product densities in the latter. An implementation of gradients at the nonhybrid density functional theory level is presented, and sample calculations demonstrate that the errors in equilibrium geometries due to the Cholesky representation of the integrals can be controlled by adjusting the decomposition threshold.     

Symmetry properties of one- and two-electron molecular integrals

The maximum numbers of distinct one- and two-electron integrals that arise in calculating the electronic energy of a molecule are discussed. It is shown that these may be calculated easily using the character table of the symmetry group of the set of basis functions used to express the wave function. Complications arising from complex group representations and from a conflict of symmetry between the basis set and the nuclear configuration are considered and illustrated by examples.     

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.     

Pre-processing two-electron integrals for efficient utilization in many-electron self-consistent field calculations

A procedure for pre-processing a random list of two-electron integrals is described. Construction of Fock-matrices with the resulting PK-file is shown to be three times faster than use of a file containing the random two-electron integrals, indices, and a branching code.     

An economic storage and processing method for two-electron integrals in LCAO-MO calculations

An economic algorithm for storing and processing two-electron integrals arising in LCAO-MO calculations is presented. The integrals are sorted prior to the SCF iterative scheme, classified according to equivalences in the orbital indices and finally stored on separate files that contain only integrals of one type. The novel approach of physically separating the integrals according to category is shown to be more efficient than random storage. Actual computing times for the new technique are tabulated for a representative number of molecular systems and compared with times obtained using previously reported methods.     

Multipole-based integral estimates for the rigorous description of distance dependence in two-electron integrals

We derive multipole-based integral estimates (MBIE) as rigorous and tight upper bounds to four-center two-electron integrals in order to account for the 1/R distance decay between the charge distributions, which is missing in the Schwarz screening commonly used in ab initio methods. Our screening criteria are valid for all angular momenta and can be formulated for any order of multipoles. We have found the expansion limited to dipoles to be sufficiently tight for estimating the integrals in Hartree-Fock and density-functional theories, while the screening effort is negligible. For, e.g., a DNA fragment with 1052 atoms and 10,674 basis functions (6-31G*) the exchange part is faster by a factor of 2.1 as compared to the Schwarz screening both within our linear exchange scheme, whereas a smaller factor of 1.3 is gained for the Coulomb part within the continuous fast multipole method. Most importantly, our new MBIE screening is perfectly suited to exploit the strong distance decay of electron-correlation effects of at least 1/R4 in atomic-orbital-based formulations of correlation methods.     

Careful neglect of small two-electron integrals applied to borazane and diazirine

A formulation is given to systematically neglect two-electron integrals in a gaussian lobe basis set smaller than a prescribed threshold. With a threshold of 10^?4 au the energies are unchanged to seven significant figures and the dipole moments to four for borazane and diazirine. This technique offers a 15-20% saving in computation time for these integrals with a negligible loss of accuracy.     

Some aspects of parallelization of two-electron integrals in molecular orbital programs

Evaluation of two-electron integrals forms a substantial part of the CPU time for any ab initio molecular orbital program. This part of the package, “MICROMOL”, is parallelized. However, this parallelization leads to only sublinear speedups (typically 3 on a 4-node machine). In view of these results, the task of development of an efficient program for two-electron integrals suitable for the parallel environment has been taken up. The program is written in FORTRAN considering specific symmetry features and application of rigorous bounds. This program is further parallelized with a good load balancing strategy. The molecules used as the test cases are: trans-butadiene, benzene, nitrobenzene, naphtalene and cytosine, with 3G and 4–31G basis sets. The results indicate that the parallel version of this program gives a typical speedup of 3.6 for a 3G basis set and approximately 3.4 for a 4–31G basis set for all the molecules tested. The sequential version of this program is ～1.2 times faster than the sequential version of MICROMOL, whereas the parallel version is ～1.4 times faster than the parallelized MICROMOL.     

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     

Self consistent molecular orbital calculations on polyatomic molecules: Gaussian approximation for two-electron integrals

A MO –LCAO –SCF treatment is performed on water, ammonia and methane using a recently proposed approximate method. The procedure is found capable of predicting total and orbital energies in close agreement with the results of accurate computations using double ζ basis sets of Slater type orbitals. A comparison is made with the results of similar approximate ab initio procedures.     

Many-Electron,multicenter integrals in superposition of correlated configurations method. I. one- and two-electron integrals

The calculations by means of the superposition of correlated configurations method (Hylleraas-CI ), that is, the combination of configuration interaction with the Hylleraas-type correlation factors, needs the effective evaluation of some nontrivial integrals. This series of papers gives the formulas for all types of integrals needed for molecular calculations when Gaussian lobe functions are used as a basis set. The formulas for two-electron integrals are given in the present paper. The preliminary results for two-electron systems are presented.