Evaluation of one‐electron molecular integrals over complete orthonormal sets of Ψα ‐ETO using auxiliary functions

By the use of expansion and one‐range addition theorems, the one‐electron molecular integrals over complete orthonormal sets of Ψ^α ‐exponential type orbitals arising in Hartree–Fock–Roothaan equations for molecules are evaluated. These integrals are expressed through the auxiliary functions in ellipsoidal coordinates. The comparison is made using Slater‐, Coulomb‐Sturmian‐, and Lambda‐type basis functions. Computation results are in good agreement with those obtained in the literature. The relationships obtained are valid for the arbitrary quantum numbers, screening constants, and location of orbitals. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Auxiliary functions for molecular integrals with Slater‐type orbitals. I. Translation methods

Many types of molecular integrals involving Slater functions can be expressed, with the ζ‐function method in terms of sets of one‐dimensional auxiliary integrals whose integrands contain two‐range functions. After reviewing the properties of these functions (including recurrence relations, derivatives, integral representations, and series expansions), we carry out a detailed study of the auxiliary integrals aimed to facilitate both the formal and computational applications of the ζ‐function method. The usefulness of this study in formal applications is illustrated with an example. The high performance in numerical applications is proved by the development of a very efficient program for the calculation of two‐center integrals with Slater functions corresponding to electrostatic potential, electric field, and electric field gradient. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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     

On the basis‐set dependence of local and integrated electron density properties: Application of a new computer program for quantum‐chemical density analysis

A new computer program for post‐processing analysis of quantum‐chemical electron densities is described. The code can work with Slater‐ and Gaussian‐type basis functions of arbitrary angular momentum. It has been applied to explore the basis‐set dependence of the electron density and its Laplacian in terms of local and integrated topological properties. Our analysis, including Gaussian/Slater basis sets up to sextuple/quadruple‐zeta order, shows that these properties considerably depend on the choice of type and number of primitives utilized in the wavefunction expansion. Basis sets with high angular momentum (l = 5 or l = 6) are necessary to achieve convergence for local properties of the density and the Laplacian. In agreement with previous studies, atomic charges defined within Bader's Quantum Theory of Atoms in Molecules appear to be much more basis‐set dependent than the Hirshfeld's stockholder charges. The former ones converge only at the quadruple‐zeta/higher level with Gaussian/Slater functions. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Auxiliary functions for molecular integrals with Slater‐type orbitals. II. Gauss transform methods

The Gauss transform of Slater‐type orbitals is used to express several types of molecular integrals involving these functions in terms of simple auxiliary functions. After reviewing this transform and the way it can be combined with the shift operator technique, a master formula for overlap integrals is derived and used to obtain multipolar moments associated to fragments of two‐center distributions and overlaps of derivatives of Slater functions. Moreover, it is proved that integrals involving two‐center distributions and irregular harmonics placed at arbitrary points (which determine the electrostatic potential, field and field gradient, as well as higher order derivatives of the potential) can be expressed in terms of auxiliary functions of the same type as those appearing in the overlap. The recurrence relations and series expansions of these functions are thoroughly studied, and algorithms for their calculation are presented. The usefulness and efficiency of this procedure are tested by developing two independent codes: one for the derivatives of the overlap integrals with respect to the centers of the functions, and another for derivatives of the potential (electrostatic field, field gradient, and so forth) at arbitrary points. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Evaluation of integrals over STOs on different centers and the complementary convergence characteristics of ellipsoidal‐coordinate and zeta‐function expansions

One‐electron integrals over three centers and two‐electron integrals over two centers, involving Slater‐type orbitals (STOs), can be evaluated using either an infinite expansion for 1/r₁₂ within an ellipsoidal‐coordinate system or by employing a one‐center expansion in spherical‐harmonic and zeta‐function products. It is shown that the convergence characteristics of both methods are complimentary and that they must both be used if STOs are to be used as basis functions in ab initio calculations. To date, reports dealing with STO integration strategies have dealt exclusively with one method or the other. While the ellipsoidal method is faster, it does not always converge to a satisfactory degree of precision. The zeta‐function method, however, offers reliability at the expense of speed. Both procedures are described and the results of some sample calculation presented. Possible applications for the procedures are also discussed. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 1–13, 1999     

Correspondence between GTO and STO molecular basis sets

We present a method for the characterization of the distance between two spaces: one generated by a Gaussian basis set, and another by a Slater basis set. The method is an extension of one previously developed for atoms that has been modified to cover molecular problems. The current version enables us to obtain Slater basis sets capable of reproducing the results (multielectronic wave functions and orbitals) obtained with Gaussian basis sets. The interest of this result arises from the fact that we will be able to profit from the effort invested in the optimization of high‐quality Gaussian basis sets. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1655–1665, 2001     

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     

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     

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.     

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     

Explicit and implicit multi‐center integrations

We present truncated expansions of multicenter one‐electron nuclear attraction and two‐electron repulsion integrals over localized basis functions in terms of one‐ and two‐center integrals of “Coulomb,” “exchange,” and “hybrid” type. Two variants are discussed: the “Explicit Multi‐center Integrations” and the “Implicit Multi‐Center Integrations” (abbreviated as “EMCI” and “IMCI”, respectively). While EMCI also deals with individual integrals, the IMCI option is the more appealing one: it enables us to evaluate the entire matrix elements of “Restricted Hartree–Fock”‐type in a very effective and chemically meaningful way. Due to the diatomic nature of our expansions, integrations over “Slater‐Type Orbitals” become well‐feasible, too. © 2012 Wiley Periodicals, Inc.     

Radial and angular correlation in heliumlike ions

In the framework of nonrelativistic variational formalism a new type of basis set is proposed, to estimate separately the effect of radial and angular correlations on the ground‐state energy for helium isoelectronic sequence H^? to Ar¹⁶⁺. Effect of radial correlation is incorporated by using multiexponential functions arising from product basis sets suitably formed out of Slater‐type one‐particle orbitals. The angular correlation can be switched on by incorporating an expansion in terms of basis involving interparticle coordinates. With a set of six‐term Slater‐type one‐particle basis and five‐term interparticle expansion, the ground‐state energy of helium is estimated as ?2.9037236 (a.u.) compared with the multiterm variational estimates ?2.9037244 (a.u.) due to Pekeris and Thakkar and Smith and Drake. Matrix elements of different operators in the ground state have been calculated and found to be in good agreement with available accurate results. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Three‐center Coulomb repulsion integrals with Slater functions

The translation method for products of two Slater functions is improved and combined with the short–long range separation in order to develop a robust and efficient algorithm for the evaluation of three‐center Coulomb integrals with Slater functions. Several tests are carried out showing that the algorithm reported here yields integrals with an absolute error below 10^?12 hartree and a computational cost of a few microseconds per integral. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Computation of two‐center Coulomb integrals over Slater‐type orbitals using elliptical coordinates

A general analytic formula is obtained for the two‐center Coulomb integrals over Slater‐type orbitals in elliptical coordinates. Finite series expansions are used in the evaluation of the radial part of the integrals. The analytic formula is expressed in terms of a product of the well‐known auxiliary functions A_k(p) and B_k(p) and incomplete gamma functions. Recursive relations for the computer evaluation of these functions are given as well. The recursive relations are stable and our computer results are in good agreement with the benchmark values given in the literature. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Coupled‐cluster frozen‐density embedding using resolution of the identity methods

Frozen‐density embedding (FDE) is combined with resolution of the identity (RI) Hartree–Fock and a RI‐variant of a second‐order approximate coupled‐cluster singles and doubles (RI‐CC2) to determine solvatochromic shifts for the lowest excitation energy of acetone and pyridazine, respectively, each solvated in different environments with total system sizes of about 2.5 nm diameter. The combination of FDE and RI‐CC2 increases efficiency and enables the calculation of numerous snapshots with 100 to 300 molecules, also allowing for larger basis sets as well as diffuse functions needed for an accurate treatment of properties. The maximum errors in the solvatochromic shifts amount up to 0.2 eV, which are similar to other approximated studies in the literature. © 2014 Wiley Periodicals, Inc.     

Two-center molecular repulsion integrals over slater functions

For the general two-electron two-center integral over Slater functions, use of the Neumann expansion for the electron-electron interaction term yields the standard auxiliary functions. These are expanded and integrated explicitly by two independent methods. The resulting simple analytic formula for the total integral is completely general, requiring only the Slater function quantum numbers and exponents and the internuclear separation. Hence all two-electron hydrid, coulomb, exchange, and one-center integrals are considered. The efficiency of calculation of this expression is compared with those of other methods, indicating an order of magnitude improvement in speed over recursion for the exchange integral.     

Calculation of transition matrix elements by nonsingular orbital transformations

A general strategy is described for the evaluation of transition matrix elements between pairs of full class CI wave functions built up from mutually nonorthogonal molecular orbitals. A new method is proposed for the counter‐transformation of the linear expansion coefficients of a full CI wave function under a nonsingular transformation of the molecular‐orbital basis. The method, which consists in a straightforward application of the Cauchy–Binet formula to the definition of a Slater determinant, is shown to be simple and suitable for efficient implementation on current high‐performance computers. The new method appears mainly beneficial to the calculation of miscellaneous transition matrix elements among individually optimized CASSCF states and to the re‐evaluation of the CASCI expansion coefficients in Slater‐determinant bases formed from arbitrarily rotated (e.g., localized or, conversely, delocalized) active molecular orbitals. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Analytical gradients of the second‐order Møller–Plesset energy using Cholesky decompositions

An algorithm for computing analytical gradients of the second‐order Møller–Plesset (MP2) energy using density fitting (DF) is presented. The algorithm assumes that the underlying canonical Hartree–Fock reference is obtained with the same auxiliary basis set, which we obtain by Cholesky decomposition (CD) of atomic electron repulsion integrals. CD is also used for the negative semidefinite MP2 amplitude matrix. Test calculations on the weakly interacting dimers of the S22 test set (Jure?ka et al., Phys. Chem. Chem. Phys. 2006, 8, 1985) show that the geometry errors due to the auxiliary basis set are negligible. With double‐zeta basis sets, the error due to the DF approximation in intermolecular bond lengths is better than 0.1 pm. The computational time is typically reduced by a factor of 6–7. © 2013 Wiley Periodicals, Inc.