Molecular symmetry and ab initio calculations: IV. Symmetry-matrix and symmetry-supermatrix in calculations of two-electron repulsion integrals

The product of two Gaussians having different centers is itself a one-center Gaussian, thus multicenter integrals with a Cartesian Gaussian basis can be reduced to one-center integrals. Recurrence relations for overlap integrals and electron repulsion integrals (ERIs) are derived at these centers. The calculations of overlap integrals and ERIs are carried out step by step from the highest symmetry case (one center) to required cases (different centers) by using the translation of Cartesian Gaussians. Full exploitation of symmetry in calculation processes can result in optimal use of these recurrence relations. Compared with the recently published algorithms, based on the recurrence relations derived by Obara and Saika [J. Chem. Phys., 84 , 3963 (1986)], the floating point operations (FLOPs) for ERI calculations (having four different centers) can be reduced by a factor of ca. 2. A significant extra saving in calculations and storage can be obtained if atoms, linear, or planar molecules are discussed. © 1997 John Wiley & Sons, Inc.     

Coupled-cluster theory in a projected atomic orbital basis

We present a biorthogonal formulation of coupled-cluster (CC) theory using a redundant projected atomic orbital (PAO) basis. The biorthogonal formulation provides simple equations, where the projectors involved in the definition of the PAO basis are absorbed in the integrals. Explicit expressions for the coupled-cluster singles and doubles equations are derived in the PAO basis. The PAO CC equations can be written in a form identical to the standard molecular orbital CC equations, only with integrals that are related to the atomic orbital integrals through different transformation matrices. The dependence of cluster amplitudes, integrals, and correlation energy contributions on the distance between the participating atomic centers and on the number of involved atomic centers is illustrated in numerical case studies. It is also discussed how the present reformulation of the CC equations opens new possibilities for reducing the number of involved parameters and thereby the computational cost.     

Calculation of 4-Center Coulomb Repulsion Matrix Elements in a Basis of Exponential Type Spherical AO Using 9-Dimensional Polyspherical Harmonics

A method for calculating 4-center Coulomb repulsion integrals in a basis of exponential type AO with regular sectorial harmonics as angular terms is proposed. The initial integrals are represented as a partial differentiation operator with respect to the Cartesian coordinates of the centers of AO, acting on the scalar function which is a 4-center integral of s functions. Differentiation is performed by calculating the Fourier transform of this scalar function in 9-dimensional Euclidean space with the help of the sectorial harmonic argument summing theorem. Thus compact representation of quantum-chemical multicenter integrals is obtained in a basis of exponential type functions with arbitrary angular parts.     

Unified analytical treatment of one-electron two-center integrals with noninteger n Slater-type orbitals

A unified treatment of one-electron two-center integrals over noninteger n Slater-type orbitals is described. Using an appropriate prolate spheroidal coordinate system with the two atomic centers as foci, all the molecular integrals are expressed by a single analytical formula which can be readily and compactly programmed. The analysis of the numerical performance of the computational algorithm is also presented. Received: 1 April 1999 / Accepted: 2 July 1999 / Published online: 2 November 1999     

Calculation of 3-Center One-Electron Integrals \left\langle {AB|r_c^{ - 1} } \right\rangle and Coulomb Repulsion Matrix Elements \left\langle {AB|CC'} \right\rangle in an Exponential Type Spherical AO Basis

A new method for calculating 3-center one-electron integrals and matrix elements of electron interaction of $\left\langle {AB|CC'} \right\rangle $ type in MO LCAO theory, where A, B, and C are the centers of exponential type AO including Slater and hydrogen-like basis functions, has been developed. Integrals of this type are reduced to double integrals over the square with edge 1. This provides the grounds for effective calculation of quantum-chemical 3-center integrals in a basis of exponential type spherical functions.     

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     

Integral evaluation algorithms and their implementation

Three different algorithms for the calculation of many center electron-repulsion integrals are discussed, all of which are considered to be economic in terms of the number of arithmetic operations. The common features of the algorithms are as follows: Cartesian Gaussian functions are used, integrals are calculated by blocks (a block being defined as the set of integrals obtainable from four given exponents on four given centers), and functions may be adopted to R(3). Adaption to molecular point group symmetry is not considered. Tables are given showing the minimum number of operations for a selection of block types allowing one to identify the theoretically most economic, and the corresponding salient features. Comments concerning the computer implementations are also given both on sealar and vector processors. In particular, the Cyber 205 is considered, a vector processor on which we have implemented what we believe to be the most efficient algorithm.     

Simplified expansion of Slater orbitals about displaced centers

A new and very simple one-range expansion of the 0s function is derived and employed as the starting point of three recurrence relations which allow the expansion of arbitrary Slater functions over displaced centers. Convergence of the expansions, both pointwise and in norm, are analyzed, and three-center nuclear attraction integrals are chosen for a further test of the formal developments and the numerical behavior of these expansions.     

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     

Analytic first derivatives for explicitly correlated,multicenter, Gaussian geminals

Variational calculations utilizing the analytic gradient of explicitly correlated Gaussian molecular integrals are presented for the ground state of the hydrogen molecule. Preliminary results serve to motivate the need for general formulas for analytic first derivatives of molecular integrals involving multicenter, explicitly correlated Gaussian geminals with respect to Gaussian exponents and coordinates of the orbital centers. Explicit formulas for analytic first derivatives of Gaussian functions containing correlation factors of the form exp(-βr_ij²) are derived and discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 991–999, 1997     

Multipole expansions for numerical orbital products

The product of two numerically defined atomic angular momentum orbitals at different centers is considered. The product can be expanded about a third center on the line segment joining the two centers. A numerical procedure for evaluating the expansion functions is developed. The application of the expansion to the evaluation of four‐center electron–electron repulsion integrals is discussed. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Energy levels of the weakly interacting radicals

Energy levels of the weakly interacting radicals are estimated on the basis of planar methyl radical interaction as a model. These levels depend on the given separation r between the centers of mass. The appropriate Schrödinger equation is solved by using the Brillouin-Wigner series method. Analytical forms for the integrals used to estimate matrix elements are derived. The principle of total momentum conservation is strictly obeyed. Some energy levels cannot be estimated as a result of divergence.     

Theoretical investigation of the magnetic interactions of Ni9 complexes

On the basis of density-functional theory (DFT) calculations, a theoretical analysis of the exchange interactions in Ni9L2(O2CMe)8{(2-py)2CO2}4, was performed, where L is a bridging ligand, OH- (1) or N3- (2). Each magnetic interaction between the Ni spin centers is analyzed for 1 and 2 in terms of exchange integrals (J values), orbital overlap integrals (T values) and natural orbitals. It was found that a J3 interaction, which is a magnetic interaction via the bridging ligand orbitals, mainly controls the whole magnetic properties, and the dominant interaction is a sigma-type orbital interaction between Ni dz2 orbitals. Further investigations on the magnetostructural correlations are performed on the J3 interactions using simplest Ni-L-Ni models. These models reproduced the magnetic interactions qualitatively well not only for the Ni9 complexes but also for other inorganic complexes. Strong correlations have been found between the magnetic orbital overlaps (T values) and the Ni-L-Ni angle. These results revealed that the difference of the magnetic properties between OH- and N3- is caused by the orbital overlap integral (T values) of the sigma-type J3 interaction pathway. The magnetic interactions are also discussed from a Hubbard model by evaluating the transfer integral (t) and on-site Coulomb integrals (U), in relation to the Heisenberg picture.     

The limits of local correlation theory: electronic delocalization and chemically smooth potential energy surfaces

Local coupled-cluster theory provides an algorithm for measuring electronic correlation quickly, using only the spatial locality of localized electronic orbitals. Previously, we showed [J. Subotnik et al., J. Chem. Phys. 125, 074116 (2006)] that one may construct a local coupled-cluster singles-doubles theory which (i) yields smooth potential energy surfaces and (ii) achieves near linear scaling. That theory selected which orbitals to correlate based only on the distances between the centers of different, localized orbitals, and the approximate potential energy surfaces were characterized as smooth using only visual identification. This paper now extends our previous algorithm in three important ways. First, locality is now based on both the distances between the centers of orbitals as well as the spatial extent of the orbitals. We find that, by accounting for the spatial extent of a delocalized orbital, one can account for electronic correlation in systems with some electronic delocalization using fast correlation methods designed around orbital locality. Second, we now enforce locality on not just the amplitudes (which measure the exact electron-electron correlation), but also on the two-electron integrals themselves (which measure the bare electron-electron interaction). Our conclusion is that we can bump integrals as well as amplitudes, thereby gaining a tremendous increase in speed and paradoxically increasing the accuracy of our LCCSD approach. Third and finally, we now make a rigorous definition of chemical smoothness as requiring that potential energy surfaces not support artificial maxima, minima, or inflection points. By looking at first and second derivatives from finite difference techniques, we demonstrate complete chemical smoothness of our potential energy surfaces (bumping both amplitudes and integrals). These results are significant both from a theoretical and from a computationally practical point of view.     

Use of STOs in Hartree–Fock calculations: Error analysis and variance-minimized pseudospectral method

In this study it is demonstrated that STO (Slater-type orbital) basis sets are particularly well suited to pseudospectral Hartree–Fock calculations. The reduction of two-electron integrals, to ones that are (at worst) equivalent to a one-electron integral over three centers, eliminates the need for slowly convergent one-center expansions. This allows all integrals to be calculated quickly and accurately in either spherical or ellipsoidal coordinates. A new variance-minimized variant of the pseudospectral method is derived and applied to a number of small closed-shell molecules. The performance of the algorithm is assessed relative to purely spectral calculations employing STO and GTO (Gaussian-type orbital) basis sets. The pseudospectral operator is used to assess the errors contained in solutions found by the purely spectral method. The suitability of a number of different de-aliasing set types is also examined. Orthogonal sets of hydrogen-like eigenfunctions were found to be optimal. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1537–1548, 1999     

Neglect of four- and approximation of one-, two-, and three-center two-electron integrals in a symmetrically orthogonalized basis

A new integral approximation for use in molecular electronic structure calculations is proposed as an alternative to the traditional neglect of diatomic differential overlap models. The similarity between the symmetrically orthogonalized and the original basis functions (assumed orthonormal within each atomic set but nonorthogonal between different centers) is used to construct a robust approximation for the two-electron integrals, with the error being quadratic in the deviation between the products of the functions. Invariance properties of this procedure are rigorously proved. Numerical studies on a representative set of molecules at valence-only minimal basis Hartree-Fock level show that the approximation introduces relatively small errors, encouraging its future application in the semiempirical field.     

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.     

A New Method of Calculating Multicenter Matrix Elements in MO LCAO Theory Using an Exponential Type Basis

Using integral representation of the product of reduced Bessel functions (RBF) specified on different centers and a new generalized integral identity for RBF one can prove that the 4-center integral of Coulomb repulsion in an exponential type AO basis may be expressed as a three-dimensional integral over the volume of a cube with an edge 1. A new method of calculating the multicenter matrix elements of quantum chemistry in an exponential AO basis is suggested based on this representation. Numerical calculations of a number of multicenter integrals using this algorithm illustrate the efficiency of the method.     

Evaluation of Hylleraas-CI atomic integrals by integration over the coordinates of one electron. II. Four-electron integrals

An alternative procedure to the classical method for evaluating the four-electron Hylleraas-CI integrals is given. The method consists of direct integration over the r ₁₂ coordinate and integration over the coordinates of one of the electrons, reducing the integrals to lower order. The method based on the earlier work of Calais and L?wdin and of Perkins is extended to the general angular case. In this way it is possible to solve all of the four-electron integrals appearing in the Hylleraas-CI method. The four-electron integrals are expanded in three-electron ones which are in turn expanded in two-electron integrals. Finally the two-electron integrals are expanded into two-electron auxiliary integrals which usually have one negative power. The use of three- and four-electron electron auxiliary integrals is avoided. Some special cases lead to one- and two-electron auxiliary integrals with negative powers which do not converge individually but do converge in combination with others. These relations and their solutions are presented, together with results of various kinds of integrals.