Molecular integrals over the gauge-including atomic orbitals. II. The Breit-Pauli interaction

Each accompanying coordinate expansion (ACE) formula is derived for each of the orbit-orbit interaction, the spin-orbit coupling, the spin-spin coupling, and the contact interaction integrals over the gauge-including atomic orbitals (GIAOs) by the use of the solid harmonic gradient (SHG) operator. Each ACE formula is the general formula derived at the first time for each of the above molecular integrals over GIAOs. These molecular integrals are arising in the Breit-Pauli two-electron interaction for a relativistic calculation. We may conclude that we can derive a certain ACE formula for any kind of molecular integral over solid harmonic Gaussian-type orbitals by using the SHG operator. The present ACE formulas will be useful, for example, for a calculation of a molecule in a uniform magnetic field, for a relativistic calculation, and so on, with the GIAO as a basis function.     

Evaluation of multicenter one-electron integrals of noninteger u screened Coulomb type potentials and their derivatives over noninteger n Slater orbitals

Multicenter integrals over noninteger n Slater type orbitals with integer and noninteger values of indices u of screened Coulomb type potentials, f(u)(eta,r)=r(u-1)e(-etar), and their first and second derivatives with respect to Cartesian coordinates of the nuclei of a molecule are described. Using complete orthonormal sets of Psi(alpha) exponential type orbitals and rotation transformation of two-center overlap integrals, these integrals are expressed through the noncentral potential functions depending on the molecular auxiliary functions A(k) and B(k). The series expansion formulas derived for molecular integrals of screened Coulomb potentials and their derivatives are especially useful for the computation of multicenter electronic attraction, electric field, and electric field gradient integrals. The convergence of series is tested for arbitrary values of parameters of potentials and orbitals.     

Evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals with integer and noninteger principal quantum numbers

A general formula has been established for the expansion of the product of two normalized associated Legendre functions centered on the nuclei a and b. This formula has been utilized for the evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals (STOs) with integer and noninteger principal quantum numbers. The formulas given in this study for the evaluation of two‐center overlap and nuclear attraction integrals show good rate of convergence and great numerical stability under wide range of quantum numbers, orbital exponents, and internuclear distances. © 2001 Wiley Periodicals, Inc. Int J Quantum Chem, 2001     

Calculation of Multicenter Nuclear Attraction and Electron Repulsion Integrals over Slater Orbitals by Fourier Transform Method Using Gegenbauer Polynomials

In this study, using complete orthonormal sets of exponential type orbitals (ETOs), a single closed analytical relation is derived for a large number of different expansions of overlap integrals over Slater type orbitals (STOs) with the same screening parameters in terms of Gegenbauer coefficients. The general formula obtained for the overlap integrals is utilized for the evaluation of multicenter nuclear attraction and electron repulsion integrals appearing in the Hartree–Fock–Roothaan equations for molecules. The formulas given in this study for the evaluation of these multicenter integrals show good rate of convergence and great numerical stability under wide range of quantum numbers, scaling parameters of STOs and internuclear distances.     

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     

Fast and stable algorithm for the analytical computation of two-center Coulomb and overlap integrals over Slater-type orbitals

Proceeding from analytical expressions for two-center kernel functions that we derived recently, we present new analytical formulas for the two-center Coulomb and overlap integrals over Slater-type orbitals. These formulas are of an exceptionally simple analytical structure and high numerical efficiency. An especially important point is that for the most frequently needed ranges of discrete quantum numbers, the formulas are completely stable in the cases of nearly equal scaling parameters or vanishing interatomic distances, except for one particular case of the Coulomb integral. No special asymptotic formulas are needed any more to compute the two-center integrals over Slater-type orbitals in these case. Furthermore, a largely recursive formulation makes the integral evaluation very economical and fast. In particular, we assess the numerical performance of a new kind of angular momentum recurrences that we have proposed in a previous article [W. Hierse and P.M. Oppeneer, J. Chem. Phys. 99 , 1278 (1993)]. © 1994 John Wiley & Sons, Inc.     

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.     

Evaluation of expansion coefficients for translation of Slater‐type orbitals using complete orthonormal sets of exponential‐type functions

By the use of exponential‐type functions (ETFs) the simpler formulas for the expansion of Slater‐type orbitals (STOs) in terms of STOs at a displaced center are derived. The expansion coefficients for translation of STOs are presented by a linear combination of overlap integrals. The final results are of a simple structure and are, therefore, especially useful for machine computations of arbitrary multielectron multicenter molecular integrals over STOs that arise in the Hartree–Fock–Roothaan approximation and also in the Hylleraas correlated wave function method for the determination of arbitrary multielectron properties of atoms and molecules. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 126–129, 2001     

Unified treatment for the evaluation of arbitrary multielectron multicenter molecular integrals over Slater‐type orbitals with noninteger principal quantum numbers

The expansion formula has been presented for Slater‐type orbitals with noninteger principal quantum numbers (noninteger n‐STOs), which involves conventional STOs (integer n‐STOs) with the same center. By the use of this expansion formula, arbitrary multielectron multicenter molecular integrals over noninteger n‐STOs are expressed in terms of counterpart integrals over integer n‐STOs with a combined infinite series formula. The convergence of the method is tested for two‐center overlap, nuclear attraction, and two‐electron one‐center integrals, due to the scarcity of the literature, and fair uniform convergence and great numerical stability under wide changes in molecular parameters is achieved. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Gaussian expansion method for molecular integrals of molecular properties

The finite Gaussian Expansion method for molecular integrals proposed by Taketa, O-ohata and Huzinaga has been extended to the integrals of molecular properties. The integral formulas of so-called moment, field and field gradient integrals have been derived. It has been numerically shown that in order to evaluate the field and the field gradient integrals based on Slater type orbitals, eight- or ten-term Gaussian expansions are sufficient but this method fails to attain sufficient effective numbers for the moment integrals.     

Symmetry-adapted integral derivatives

A procedure, based on double coset decompositions, is described for reducing formulas for derivatives (with respect to nuclear coordinates) of integrals over symmetry-adapted orbitals to symmetry-distinct integral derivatives over atomic orbitals. The procedure is applicable to any finite point group and to integral derivatives of any order.     

Slater-orbital molecular integrals by numerical Fourier-transform methods. II. Four-center integrals over is orbitals

The four-center nonplanar electron repulsion integrals over 1s Slater-type atomic orbitals are considered by a numerical Fourier-transform method. It is shown that the highly oscillating integrand appearing in the Fourier inversion formula could be successfully treated by using Tchebyscheff quadrature. The resulting formulas are thoroughly discussed with particular emphasis on their numerical features and convergence properties. It follows that the aforementioned integrals may be calculated with a good accuracy with a moderate amount of computing time.     

偶极积分与Boys法分子轨道定域化

定域分子轨道在分子体系的化学图象和物理图象之间充当重要的桥梁作用,它的产生依赖于定域化准则,其中最普遍使用的是Foster-Boys和Edmiston-Ruedenberg(E—R)提出的两种定域化准则。这两种定域化准则是等价的,因而结果也是一致的。但对于E—R定域化来说,由于涉及到大量的多中心积分的计算,计算极为费时,因而远不如Foster-Boys定     

Calculation of molecular electric and magnetic multipole moment integrals of integer and noninteger n Slater orbitals using overlap integrals

Closed formulas are established for the magnetic multipole moment integrals of integer and noninteger n Slater‐type orbitals (ISTOs and NISTOs) in terms of electric multipole moment integrals for which the analytic expressions through the overlap integrals with ISTOs and NISTOs are derived. The overlap integrals are evaluated by the use of auxiliary functions. Using the derived expressions the multipole moment integrals, and therefore the electric and magnetic properties of molecules, can be evaluated most efficiently and accurately. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Spherically symmetrical properties of two-center overlap integrals over arbitrary atomic orbitals and translation coefficients for Slater-type orbitals

The orthogonality relations are derived for the rotation coefficients of two-center overlap integrals over arbitrary atomic orbitals (AAOs) and expansion coefficients for translation of Slater-type orbitals (STOs). Using these formulas, a very interesting theorem regarding the angular dependence is established. If we add the products of all the overlap integrals or all the translation coefficients with the same n and l values, but different m values, the result is independent of orientation. The final results are of a simple structure and are, therefore, especially useful for machine computations of multielectron multicenter molecular integrals by expanding one- and two-center electron charge density over STOs in terms of STOs about a new center.     

Calculation of electric multipole moment integrals with the different screening parameters via the Fourier transform method

Three‐center electric multipole moment integrals over Slater‐type orbitals (STOs) can be evaluated by translating the orbitals on one center to the other and reducing the system to an expansion of two‐center integrals. These are then evaluated using Fourier transforms. The resulting expression depends on the overlap integrals that can be evaluated with the greatest ease. They involve expressions for STO with different screening parameters that are known analytically. This work gives the overall expressions analytically in a compact form, based on Gegenbauer polynomials. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

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     

A new expression for the vibrational overlap integrals between arbitrary one-dimensional harmonic oscillators

A new, compact formula for the vibrational overlap integrals between two harmonic potentials with both arbitrary curvatures and equilibrium positions is derived in a systematic manner with its several derivative formulas. Some mathematical properties of the overlap integrals obtained are discussed.     

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.