A new algorithm to evaluate bielectronic integrals with 1s Slater‐type orbitals obtained by the integral transforms

This article presents a variation of the integral transform method to evaluate multicenter bielectronic integrals (12|34), with 1s Slater‐type orbitals. It is proved that it is possible to define, out of the expression of (12|34) given by the integral transform method, a function F(q) that has the property of having a unique Q, such that F(Q) = (12|34). Therefore, F(q) may be used to calculate (12|34). It is shown that the evaluation of F(Q) turns out to be simpler than the three‐dimensional integral involved in the calculation of (12|34), and an algorithm is presented to calculate Q. The results show that relative errors on the order of 10^?3 or lower are obtained very efficiently. In addition, it is shown that the proposed algorithm is very stable. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Molecular integrals over Gaussian-type geminal basis functions

We present formulas for the evaluation of molecular integrals over basis functions with an explicit Gaussian dependence on interelectronic coordinates. These formulas use expansions in Hermite Gaussian functions and represent an extension to the work of McMurchie and Davidson to two-electron basis functions. Integrals that depend on the coordinates of up to four electrons are discussed explicitly. A key feature of this approach is that it allows full exploitation of the shell structure of the orbital part of the basis. Received: 24 February 1997 / Accepted: 4 March 1997     

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.     

Quasi‐random integration in quantum chemistry: Efficiency,stability, and application to the study of small atoms and molecules constrained in spherical boxes

This project consists of two parts. In the first part, a series of test calculations is performed to verify that the integrals involved in the determination of atomic and molecular properties by standard self‐consistent field (SCF) methods can be obtained through Halton, Korobov, or Hammersley quasi‐random integration procedures. Through these calculations, we confirm that all three methods lead to results that meet the levels of precision required for their use in the calculation of properties of small atoms or molecules at least at a Hartree–Fock level. Moreover, we have ensured that the efficiency of quasi‐random integration methods that we have tested is Halton=Korobov>Hammersley?pseudo‐random. We also find that these results are comparable to those yielded by ordinary Monte Carlo (pseudo‐random) integration, with a calculation effort of two orders of smaller magnitude. The second part, which would not have been possible without the integration method previously analyzed, contains a first study of atoms constrained in spherical boxes through SCF calculations with basis functions adapted to the features of the problem: Slater‐type orbitals (STOs) trimmed by multiplying them by a function that yields 1 for 0 < r < (R‐δ), polynomial values for (R‐δ) < r < R and null for r > R, R being the radius of the box and δ a variationally determined interval. As a result, we obtain a equation of state for electrons of small systems, valid just in the limit of low temperatures, but fairly simple. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Evaluation of Multicenter Electronic Attraction, Electric Field and Electric Field Gradient Integrals with Screened and Nonscreened Coulomb Potentials over Integer and Noninteger n Slater Orbitals

By the use of complete orthonormal sets of -exponential-type orbitals, where ( = 1, 0, –1, –2,...) the multicenter electronic attraction (EA), electric field (EF) and electric field gradient (EFG) integrals of nonscreened and Yukawa-like screened Coulomb potentials are expressed through the two-center overlap integrals with the same screening constants and the auxiliary functions introduced in our previous paper (I.I. Guseinov, J. Phys. B, 3 (1970) 1399). The recurrence relations for auxiliary functions are useful for the calculation of multicenter EA, EF and EFG integrals for arbitrary integer and noninteger values of principal quantum numbers, screening constants, and location of slater-type orbitals. The convergence of the series is tested by calculating concrete cases.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals. III. Auxiliary functions \documentclass{article}\pagestyle{empty}\begin{document}$Q_{nn'}^{q}$\end{document} and \documentclass{article}\pagestyle{empty}\begin{document}$G_{-nn'}^{q}$\end{document}

The auxiliary functions $Q_{nn'}^{q}(p,pt)$ and $G_{-nn'}^{q}(p_{a},p,pt)$ which are used in our previous paper [Guseinov, I. I.; Mamedov, B. A. Int J Quantum Chem 2001, 81, 117] for the computation of multicenter electron‐repulsion integrals over Slater‐type orbitals (STOs) are discussed in detail, and the method is given for their numerical computation. The present method is suitable for all values of the parameters p_a, p, and pt. Three‐ and four‐center electron‐repulsion integrals are calculated for extremely large quantum numbers using relations for auxiliary functions obtained in this paper. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Evaluation of two‐center,three‐ and four‐electron integrals over Slater‐type orbitals in elliptical coordinates

An algorithm for evaluation of two‐center, three‐electron integrals with the correlation factors of the type rr and rrr as well as four‐electron integrals with the correlation factors rrr and rrr in the Slater basis is presented. This problem has been solved here in elliptical coordinates, using the generalized and modified form of the Neumann expansion of the interelectronic distance function r for k ≥ ?1. Some numerical results are also included. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

The W algorithm and the D transformation for the numerical evaluation of three‐center nuclear attraction integrals

Three‐center nuclear attraction integrals over exponential‐type functions are required for ab initio molecular structure calculations and density functional theory (DFT). These integrals occur in many millions of terms, even for small molecules, and they require rapid and accurate numerical evaluation. The use of a basis set of B functions to represent atomic orbitals, combined with the Fourier transform method, led to the development of analytic expressions for these molecular integrals. Unfortunately, the numerical evaluation of the analytic expressions obtained turned out to be extremely difficult due to the presence of two‐dimensional integral representations, involving spherical Bessel integral functions. % The present work concerns the development of an extremely accurate and rapid algorithm for the numerical evaluation of these spherical Bessel integrals. This algorithm, which is based on the nonlinear D transformation and the W algorithm of Sidi, can be computed recursively, allowing the control of the degree of accuracy. Numerical analysis tests were performed to further improve the efficiency of our algorithm. The numerical results section demonstrates the efficiency of this new algorithm for the numerical evaluation of three‐center nuclear attraction integrals. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Some exact expressions for the mean and higher moments of functions of sample moments

Examples of exact expressions for the moments (mainly of the mean) of functions of sample moments are given. These provide checks on alternative developments such as asymptotic series for n, and simulation processes. Exact expressions are given for the mean of the square of the sample coefficient of variation, particularly in uniform sampling; Frullani integrals studied by G. H. Hardy arise. It should be kept in mind that exact results for (joint) moment generating functions (mgfs) are of interest as they produce a means of obtaining exact results for (cross) moments—including moments with negative indices. Thus an exact expression for the joint mgf of the 1st two noncentral moments can be used to obtain the mean of the (c.v.)² (but not for the mean of the c..). A general expression is given for the moment generating function of the sample variance. The limitations of Fisher's symbolic formula for the characteristic function of sample moments (or more general statistics) are noted.This research was sponsored by the Applied Mathematical Sciences Research program, Office of Energy Research, U. S. Department of Energy under contract DE-AC0584OR21400 with the Martin Marietta Energy Systems. Inc.     

New analytical expressions,symmetry relations and numerical solutions for the rotational overlap integrals

In this article, extremely simple analytical formulas are obtained for rotational overlap integrals which occur in integrals over two reduced rotation matrix elements. The analytical derivations are based on the properties of the Jacobi polynomials and beta functions. Numerical results and special values for rotational overlap integrals are obtained by using symmetry properties and recurrence relationships for reduced rotation matrix elements. The final results are of surprisingly simple structures and very useful for practical applications. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012