A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians

Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals. For two- and three-center two-electron integrals, the solid-harmonic Hermite scheme is particularly efficient, significantly reducing the cost relative to the Cartesian scheme.     

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.     

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.     

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.     

Evaluation of Hylleraas-CI atomic integrals by integration over the coordinates of one electron. I. Three-electron integrals

A method to evaluate the nonrelativistic electron-repulsion, nuclear attraction and kinetic energy three-electron integrals over Slater orbitals appearing in Hylleraas-CI (Hy-CI) electron structure calculations on atoms is shown. It consists on the direct integration over the interelectronic coordinate r _ij and the sucessive integration over the coordinates of one of the electrons. All the integrals are expressed as linear combinations of basic two-electron integrals. These last are solved in terms of auxiliary two-electron integrals which are easy to compute and have high accuracy. The use of auxiliary three-electron ones is avoided, with great saving of storage memory. Therefore this method can be used for Hy-CI calculations on atoms with number of electrons N ≥ 5. It has been possible to calculate the kinetic energy also in terms of basic two-electron integrals by using the Hamiltonian in Hylleraas coordinates, for this purpose some mathematical aspects like derivatives of the spherical harmonics with respect to the polar angles and recursion relations are treated and some new relations are given.     

The SHARK integral generation and digestion system

In this paper, the SHARK integral generation and digestion engine is described. In essence, SHARK is based on a reformulation of the popular McMurchie/Davidson approach to molecular integrals. This reformulation leads to an efficient algorithm that is driven by BLAS level 3 operations. The algorithm is particularly efficient for high angular momentum basis functions (up to L = 7 is available by default, but the algorithm is programmed for arbitrary angular momenta). SHARK features a significant number of specific programming constructs that are designed to greatly simplify the workflow in quantum chemical program development and avoid undesirable code duplication to the largest possible extent. SHARK can handle segmented, generally and partially generally contracted basis sets. It can be used to generate a host of one- and two-electron integrals over various kernels including, two-, three-, and four-index repulsion integrals, integrals over Gauge Including Atomic Orbitals (GIAOs), relativistic integrals and integrals featuring a finite nucleus model. SHARK provides routines to evaluate Fock like matrices, generate integral transformations and related tasks. SHARK is the essential engine inside the ORCA package that drives essentially all tasks that are related to integrals over basis functions in version ORCA 5.0 and higher. Since the core of SHARK is based on low-level basic linear algebra (BLAS) operations, it is expected to not only perform well on present day but also on future hardware provided that the hardware manufacturer provides a properly optimized BLAS library for matrix and vector operations. Representative timings and comparisons to the Libint library used by ORCA are reported for Intel i9 and Apple M1 max processors.     

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     

应用Monte Carlo平均值法计算STO双中心重迭积分

把Monte Carlo方法引进STO双中心重叠积分的计算中,结果表明,它不仅计算简便、快速、很容易在计算机上实现,而且具有较高的精确度,有望推广应用于更复杂的多中心分子积分中.     

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.     

Fourier transform general formula for systematic potentials

For calculating molecular integrals of systematic potentials, a three‐dimensional (3D) Fourier transform general formula can be derived, by the use of the Abel summation method. The present general formula contains all 3D Fourier transform formulas which are well known as Bethe–Salpeter formulas (Bethe and Salpeter, Handbuch der Physik, Bd. XXXV, 1957) as special cases. It is shown that, in several of the Bethe–Salpeter formulas, the integral does not converge in the meaning of the Riemann integral but converges in the meaning of a hyper function as the Schwartz distribution. For showing an effectiveness of the present general formula, the convergence condition of molecular integrals is derived generally for all of the present potentials. It is found that molecular integrals can be converged in the meaning of the Riemann integral for the present potentials, except for those for extra super singular potentials. It is also found that the convergence condition of molecular integrals over the Slater‐type orbitals is exactly the same as that of the corresponding integrals over the Gaussian‐type orbitals for the present systematic potentials. For showing more effectiveness, the molecular integral over the gauge‐including atomic orbitals is derived for the magnetic dipole‐same‐dipole interaction. © 2012 Wiley Periodicals, Inc.     

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.     

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     

Evaluation of one-electron integrals for arbitrary operators V(r) over Cartesian Gaussians: Application to inverse-square distance and Yukawa operators

A general procedure is presented for generating one-electron integrals over any arbitrary potential operator that is a function of radial distance only. The procedure outlines that for a nucleus centered at point C integrals over Cartesian Gaussians can be written as linear combinations of 1-D integrals. These Cartesian Gaussian functions are expressed in a compact form involving easily computed auxiliary functions. It is well known that integrals over the Coulomb operator can be expressed in terms of F_n(T) integrals, where By means of a substitution for F_n(T) by other simple functions, algorithms that form integrals over an arbitrary function can be generated. Formation of such integrals is accomplished with minor editing of existing code based on the McMurchie–Davidson formalism. Further, the method is applied using the inverse-square distance and Yukawa potential operators V(r) over Cartesian Gaussian functions. Thus, the proposed methodology covers a large class of one-electron integrals necessary for theoretical studies of molecular systems by ab initio calculations. Finally, by virtue of the procedure's recursive nature it provides us with an efficient scheme of computing the proposed class of one-electron integrals. © 1993 John Wiley & Sons, Inc.     

Nonlinear transformation methods for accelerating the convergence of Coulomb integrals over exponential type functions

It is well known that in any ab initio molecular orbital (MO) calculation, the major task involves the computation of molecular integrals, among which the computation of Coulomb integrals are the most frequently encountered. As the molecular system gets larger, computation of these integrals becomes one of the most laborious and time consuming steps in molecular systems calculation. Improvement of the computational methods of molecular integrals would be indispensable to a further development in computational studies of large molecular systems. The atomic orbital basis functions chosen in the present work are Slater type functions. These functions can be expressed as finite linear combinations of B functions which are suitable to apply the Fourier transform method. The difficulties of the numerical evaluation of the analytic expressions of the integrals of interest arise mainly from the presence of highly oscillatory semi-infinite integrals. In this work, we present a generalized algorithm based on the nonlinear transformation of Sidi, for a precise and fast numerical evaluation of molecular integrals over Slater type functions and over B functions. Numerical results obtained for the three-center two-electron Coulomb and hybrid integrals over B functions and over Slater type functions. Comparisons with numerical results obtained using alternatives approaches and an existing code are listed.     

Analytical GIAO and hybrid-basis integral derivatives: application to geometry optimization of molecules in strong magnetic fields

Analytical integral evaluation is a central task of modern quantum chemistry. Here we present a general method for evaluating differentiated integrals over standard Gaussian and mixed Gaussian/plane-wave hybrid orbitals. The main idea is to have a representation of basis sets that is flexible enough to enable differentiated integrals to be reinterpreted as standard integrals over modified basis functions. As an illustration of the method, we report a very simple implementation of Hartree-Fock level geometrical derivatives in finite magnetic fields for gauge-origin independent atomic orbitals, within the London program. As a quantum-chemical application, we optimize the structure of helium clusters and some well-known covalently bound molecules (water, ammonia and benzene) subject to strong magnetic fields.     

On the use of spatial symmetry in atomic-integral calculations: An efficient permutational approach

The minimal number of independent nonzero atomic integrals that occur over arbitrarily oriented basis orbitals of the form ?(r) · Y_lm(Ω) is theoretically derived. The corresponding method can be easily applied to any point group, including the molecular continuous groups C_∞v and D_∞h. On the basis of this (theoretical) lower bound, the efficiency of the permutational approach in generating sets of independent integrals is discussed. It is proved that lobe orbitals are always more efficient than the familiar Cartesian Gaussians, in the sense that GLO s provide the shortest integral lists. Moreover, it appears that the new axial GLO s often lead to a number of integrals, which is the theoretical lower bound previously defined. With AGLO s, the numbers of two-electron integrals to be computed, stored, and processed are divided by factors 2.9 (NH₃), 4.2 (C₅H₅), and 3.6 (C₆H₆) with reference to the corresponding CGTO s calculations. Remembering that in the permutational approach, atomic integrals are directly computed without any four-indice transformation, it appears that its utilization in connection with AGLO s provides one of the most powerful tools for treating symmetrical species.     

Intracule functional models. V. Recurrence relations for two-electron integrals in position and momentum space

The approach used by Ahlrichs [Phys. Chem. Chem. Phys., 2006, 8, 3072] to derive the Obara-Saika recurrence relation (RR) for two-electron integrals over Gaussian basis functions, is used to derive an 18-term RR for six-dimensional integrals in phase space and 8-term RRs for three-dimensional integrals in position or momentum space. The 18-term RR reduces to a 5-term RR in the special cases of Dot and Posmom intracule integrals in Fourier space. We use these RRs to show explicitly how to construct Position, Momentum, Omega, Dot and Posmom intracule integrals recursively.     

New program for molecular calculations with Slater‐type orbitals

A new program for computing all the integrals appearing in molecular calculations with Slater‐type orbitals (STO) is reported. This program follows the same philosophy as the reference pogram previously reported but introduces two main changes: Local symmetry is profited to compute all the two‐electron integrals from a minimal set of seed integrals, and a new algorithm recently developed is used for computing the seed integrals. The new code reduces between one and two orders of magnitude the computational cost in most polyatomic systems. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 148–153, 2001     

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     

Beyond the MNDO model: Methodical considerations and numerical results

It is suggested to improve the MNDO model by the explicit inclusion of valence-shell orthogonalization corrections, penetration integrals, and effective core potentials (ECPs) in the one-center part of the core Hamiltonian matrix. Guided by analytic formulas and numerical ab initio results, the orthogonalization corrections are expressed in terms of the resonance integrals that are represented by a new empirical parametric function. All two-center Coulomb interactions and ECP integrals are evaluated analytically in a Gaussian basis followed by a uniform Klopman–Ohno scaling. One particular implementation of the proposed NDDO SCF approach is described and parameterized for the elements H, C, N, O, and F. In a statistical evaluation of ground-state properties, this implementation shows slight but consistent improvements over MNDO, AM1, and PM3. Significant improvements are found for excited states, transition states, and strong hydrogen bonds. Possible further enhancements of the current implementation are discussed. © 1993 John Wiley & Sons, Inc.