Master formulas for two‐ and three‐center one‐electron integrals involving Cartesian GTO,STO, and BTO

As a first application of the shift operators method we derive master formulas for the two‐ and three‐center one‐electron integrals involving Gaussians, Slater, and Bessel basis functions. All these formulas have a common structure consisting in linear combinations of polynomials of differences of nuclear coordinates. Whereas the polynomials are independent of the type (GTO, BTO, or STO) of basis functions, the coefficients depend on both the class of integral (overlap, kinetic energy, nuclear attraction) and the type of basis functions. We present the general expression of polynomials and coefficients as well as the recurrence relations for both the polynomials and the whole integrals. Finally, we remark on the formal and computational advantages of this approach. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 83–93, 2000     

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     

One-center expansion for pseudopotential matrix elements

Semilocal pseudopotential operators can be expressed as a linear combination of nonlocal (projection) operators. Pseudopotential operator integrals over a molecular basis set are therefore reduced to linear combinations of overlap integrals products. Molecular calculations indicate that sufficient precision can be achieved with a limited number of nonlocal operators. Analytic derivatives of pseudopotential integrals are easily deduced and implemented in a standard quantum chemistry program.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals I. Single‐center expansion method

The multicenter charge‐density expansion coefficients [I. I. Guseinov, J Mol Struct (Theochem) 417 , 117 (1997)] appearing in the molecular integrals with an arbitrary multielectron operator were calculated for extremely large quantum numbers of Slater‐type orbitals (STOs). As an example, using computer programs written for these coefficients, with the help of single‐center expansion method, some of two‐electron two‐center Coulomb and four‐center exchange electron repulsion integrals of Hartree–Fock–Roothaan (HFR) equations for molecules were also calculated. Accuracy of the results is quite high for the quantum numbers, screening constants, and location of STOs. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 146–152, 2000     

An accurate numerical multicenter integration for molecular orbital theory

A powerful and accurate numerical three‐dimensional integration scheme was developed especially for molecular orbital calculations. A multicenter integral is decomposed into the sum of single‐center integrals using nuclear weight functions and calculated using Gaussian quadrature rules. The decomposed single‐center integrands show strong anisotropy. With a careful selection of the Gaussian quadrature rule according to the anisotropy, it is possible to obtain an accuracy of 13 digits with a small number of integration points for the overlap integrals, normalization integrals, and molecular integrals for the hydrogen molecule. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 509–523, 1999     

Attractive electron–electron interactions within robust local fitting approximations

An analysis of Dunlap's robust fitting approach reveals that the resulting two‐electron integral matrix is not manifestly positive semidefinite when local fitting domains or non‐Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four‐center two‐electron integrals based on the resolution‐of‐the‐identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair‐atomic resolution‐of‐the‐identity (PARI) approach, atomic‐orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product. Numerical tests indicate that in 1% or less of all Hartree–Fock and Kohn–Sham calculations, the indefinite integral matrix causes nonconvergence in the self‐consistent‐field iterations. In these cases, the two‐electron contribution to the total energy becomes negative, meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the standard Coulomb‐metric RI method. The speed‐up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple‐zeta quality basis set. A positive semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky‐decomposition techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calculations. © 2013 Wiley Periodicals, Inc.     

A full‐pivoting algorithm for the Cholesky decomposition of two‐electron repulsion and spin‐orbit coupling integrals

A significant reduction in the computational effort for the evaluation of the electronic repulsion integrals (ERI) in ab initio quantum chemistry calculations is obtained by using Cholesky decomposition (CD), a numerical procedure that can remove the zero or small eigenvalues of the ERI positive (semi)definite matrix, while avoiding the calculation of the entire matrix. Conversely, due to its antisymmetric character, CD cannot be directly applied to the matrix representation of the spatial part of the two‐electron spin‐orbit coupling (2e‐SOC) integrals. Here, we present a computational strategy to achieve a Cholesky representation of the spatial part of the 2e‐SOC integrals, and propose a new efficient CD algorithm for both ERI and 2e‐SOC integrals. The proposed algorithm differs from previous CD implementations by the extensive use of a full‐pivoting design, which allows a univocal definition of the Cholesky basis, once the CD δ threshold is made explicit. We show that is the upper limit for the errors affecting the reconstructed 2e‐SOC integrals. The proposed strategy was implemented in the ab initio program Computational Emulator of Rare Earth Systems (CERES), and tested for computational performance on both the ERI and 2e‐SOC integrals evaluation. © 2017 Wiley Periodicals, Inc.     

On numerical solution of the integral Hartree-Fock equations for molecules by the method of MO LCAO in the momentum space

To develop a numerical solution of mentioned equations the method of factorized projection of integral operator kernel is applied. All matrix elements of the method are calculated analytically, being expressed in terms of two types of standard integrals: the overlap integrals and one-electron Coulomb integrals. To calculate the integrals we used the O(4)-symmetry of hydrogen-like atomic orbitals as well as operational technique of differentiation with respect to scalar and vector parameters.     

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     

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.     

Energy and energy gradient matrix elements with N-particle explicitly correlated complex Gaussian basis functions with L=1

In this work we consider explicitly correlated complex Gaussian basis functions for expanding the wave function of an N-particle system with the L=1 total orbital angular momentum. We derive analytical expressions for various matrix elements with these basis functions including the overlap, kinetic energy, and potential energy (Coulomb interaction) matrix elements, as well as matrix elements of other quantities. The derivatives of the overlap, kinetic, and potential energy integrals with respect to the Gaussian exponential parameters are also derived and used to calculate the energy gradient. All the derivations are performed using the formalism of the matrix differential calculus that facilitates a way of expressing the integrals in an elegant matrix form, which is convenient for the theoretical analysis and the computer implementation. The new method is tested in calculations of two systems: the lowest P state of the beryllium atom and the bound P state of the positronium molecule (with the negative parity). Both calculations yielded new, lowest-to-date, variational upper bounds, while the number of basis functions used was significantly smaller than in previous studies. It was possible to accomplish this due to the use of the analytic energy gradient in the minimization of the variational energy.     

Symmetry reduction of the matrix elements of a two-particle operator

Local coordinate systems are chosen for each quadruple of atoms relative to a four-center integral, in order to avoid linear combinations of orbitals when symmetry operations perform on an orbital. This choice can utilize the complete molecular symmetry to attain the optimal number of symmetry-unique integrals and to construct two-particle matrix elements by multiplying symmetry-unique integrals, called the “standard four-center integrals,” by the corresponding coefficients, called the “C coefficients.” A simple algorithm to use the complete molecular symmetry to reduce calculations of molecular matrix elements is outlined for general highly symmetric molecules. A tetrahedral molecule is analyzed. © 1996 John Wiley & Sons, Inc.     

Calculation of molecular integrals for systems confined by hard spherical walls: Use of the single‐center expansion of floating spherical gaussians

Assuming a gaussian basis set representation of atomic and molecular wave functions, the single‐center expansion of off‐centered spherical gaussian orbitals is exploited to calculate the one and two‐electron integrals for multielectronic atoms and molecules confined within hard spherical walls. As a validating test, the ground‐state energy of a helium atom positioned off‐center in a spherical box is calculated by applying the simplest form of the floating spherical gaussian orbital (FSGO) scheme, i.e., the use of a primitive basis set consisting of a single FSGO per electron pair. Comparison with corresponding recent accurate calculations gives supporting evidence of the adequacy of the method for its application to more elaborate gaussian‐type basis set representations for confined atoms and molecules. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 271–278, 2001     

Spin-orbit (core) and core potential integrals

Molecular integral formulas and corresponding computational algorithms are developed for the relativistic spin-orbit and core potential operators that are obtained from atomic relativistic calculations by means of the effective core potential procedure. Much use is made of earlier work on core potential integrals by McMurchie and Davidson. The resulting computer code has been made part of the ARGOS (Argonne, Ohio State) program from the C?OLUMBUS suite of programs, which computes the needed integrals over symmetry-adapted combinations of generally contacted Gaussian atomic orbitals.     

Near‐exact supersymmetric partner potentials: Construction and exploitation

Precise supersymmetric (SUSY) partner potentials can be generated only for exactly solvable problems of the stationary Schrödinger equation. This is a severe restriction, as most problems are not amenable to exact solutions. We employ here a linear variational strategy to explicitly construct approximate SUSY partners of a few common, not exactly solvable potentials and subsequently examine their properties to explore the advantages in practical implementation. The efficacy of our proposed scheme is commendable. We demonstrate that, for symmetric potentials, the constructed partners may be so good that the overall recipe has the nicety of generating the whole eigenspectrum by employing only half of the full Hilbert space functions. A similar strategy is shown to work for the odd states too, with proper boundary conditions. Pilot calculations involve a number of low‐lying states of some mixed oscillator and double‐well potentials. Analysis of the results reveals a few interesting features of the problem of construction of approximate SUSY partners and their practical use. Particularly, we identify places where the operator‐level approximations are involved and how far they affect the bounding properties of energies that are obtained as eigenvalues of a matrix diagonalization problem associated with linear variations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Hartree-Fock soft Coulomb hole to estimate correlation energies in atoms, molecules and polymers

We have shown that the empirical correction introduced into the Hartree-Fock method to calculate correlation energies for atoms and therefore to remove the error caused by the so-called Coulomb hole can be extended from atoms to molecules and polymers. A reformulation was required of the necessary parameter representation. The reparametrization has been performed staying as close as possible to the original expressions for atoms reported by Chakravorty and Clementi (S.J. Chakravorty and E. Clementi, Phys. Rev. A, 39 (1989) 2290). In addition to their work, where the correlation energy has been calculated with the self-consistent Hartree-Fock wavefunction and the correction integrals, we have performed investigations, including the perturbation operator in the Fock operator, so that the total energy also contains the correlation energy. The applications of this approach to atoms and molecules show that the total electron correlation energies and ionization potentials calculated as differences of total energies can be obtained very satisfactorily. On the basis of the reported calculations it turns out that one obtains better agreement with reference values of more sophisticated calculations when the correction integrals are used to build up the Fock matrix. Furthermore we have found that the magnitude of the correlation energy depends only weakly on the size of the basis sets, which makes this empirical method very attractive for its application to large molecular and polymeric systems.     

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.     

A quartic force field coordinate substitution scheme using hyperbolic sine coordinates

Quartic force fields (QFF) are currently the most cost‐effective method for the approximation of potential energy surfaces for the calculation of anharmonic vibrational energies. It is known, although, that its performance can be less than satisfactory due to limitations related to slow convergence of the series. In this article, we present a coordinate substitution scheme using a combination of Morse and sinh coordinates, well adapted for its use with cartesian normal coordinates. We derive expressions for analytical integrals for use in VSCF and VCI calculations and show that the simultaneous substitution of symmetric and antisymmetric normal coordinates by Morse and sinh coordinates, respectively, significantly improves the vibrational transition frequencies for these modes in a well‐balanced fashion. The accuracy of this substitution scheme is demonstrated by comparing one and two‐dimensional sections of substituted and unsubstituted QFF with ab initio potential energy grids, as well as with vibrational energy calculations using as test cases two well‐studied benchmark molecules: water and formaldehyde. We conclude that the coordinate substitution scheme presented constitutes a very attractive alternative to simple QFFs in the context of cartesian normal coordinates.