Overlap integrals of B functions

Two different methods for the evaluation of overlap integrals of B functions with different scaling parameters are analyzed critically. The first method consists of an infinite series expansion in terms of overlap integrals with equal scaling parameters [14]. The second method consists of an integral representation for the overlap integral which has to be evaluated numerically. Bhattacharya and Dhabal [13] recommend the use of Gauss-Legendre quadrature for this purpose. However, we show that Gauss-Jacobi quadrature gives better results, in particular for larger quantum number. We also show that the convergence of the infinite series can be improved if suitable convergence accelerators are applied. Since an internal error analysis can be done quite easily in the case of an infinite series even if it is accelerated, whereas it is very costly in the case of Gauss quadratures, the infinite series is probably more efficient than the integral representation. Overlap integrals of all commonly occurring exponentially declining basis functions such as Slater-type functions, can be expressed by finite sums of overlap integrals of B functions, because these basis functions can be represented by linear combinations of B functions.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Calculation of multicenter electronic attraction,electric field and electric field gradient integrals of Coulomb potential over integer and noninteger <Emphasis Type="Italic">n</Emphasis> Slater orbitals

With the help of complete orthonormal sets of - ETOs, where = 1,0, – 1, – 2, ... a large number of series expansion formulas for the multicenter electronic attraction (EA), electric field (EF) and electric field gradient (EFG) integrals of integer and noninteger n Slater type orbitals (ISTOs and NISTOs) is established through the overlap integrals with the same screening constants and the new central and noncentral interaction potentials depending on the coordinates of the nuclei of a molecule are introduced. The convergence of the series is tested by calculating concrete cases for arbitrary quantum numbers, screening constants and location of ISTOs and NISTOs.     

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     

Evaluation of screened nuclear attraction and electron repulsion molecular integrals over Gaussian basis functions

Analytical solutions to the Yukawa-like screened Coulomb nuclear attraction and electron repulsion molecular basic integrals, as well as to the basic integrals required to compute the virial coefficient, over Gaussian basis functions, are derived and cast into a practical closed form, suitable to interface with modern codes for the calculation of molecular electronic structure. © 1997 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.     

Use of Filter-Steinborn B and Guseinov $${Q_{ns}^q }$$ auxiliary functions in evaluation of two-center overlap integrals over Slater type orbitals

An efficient method for computing overlap integral over Slater type orbitals based on the B Filter-Steinborn and Guseinov auxiliary functions is presented. The final results are expressed through the binomial coefficients with the help of which the overlap integrals can be evaluated efficiently and accurately. The results of calculation are in good agreement with those obtained by other method for arbitrary principal quantum numbers and different screening constants. An erratum to this article can be found at     

Inner and outer radial density functions in many-electron atoms

When any two electrons are considered simultaneously, the radial density function D(r) in many-electron atoms is shown to be rigorously separated into inner D _<(r) and outer D _>(r) radial densities. Accordingly, radial properties such as the electron–nucleus attraction energy V _en and the diamagnetic susceptibility χ _d are the sum of the inner and outer contributions. The electron–electron repulsion energy V _ee has an approximate relation with the minus first moment of the outer density D _>(r). For the 102 atoms He through Lr in their ground states, different characteristics of local maxima in the radial densities D _<(r), D _>(r), and D(r) are reported based on the numerical Hartree-Fock wave functions. Relative contributions of the inner and outer components to V _en and are also discussed for these atoms.     

Fourier analysis,correlation functions and nonadiabatic electron transfer: Wavepackets and exact representations

Summary We investigate the validity of several common approximations in the analysis of nonadiabatic intramolecular electron transfer rate constants. Utilizing the Fourier representation of the golden rule form, we study the evolution of the vibrational correlation function that represents the density-of-states-weighted Franck-Condon factor. In particular, we test the validity of the perturbation theoretic golden rule form and of the Gaussian wavepacket representation for the vibrational wavefunctions against numerically exact quantum mechanical propagations. Although specific cases are found in which both of these break down, for a wide range of conditions (including anharmonic behavior and frequency changes), both the Gaussian wavepacket representation and the golden rule are excellent approximations.     

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.     

Convergence accelerator approach for the evaluation of some three‐electron integrals containing explicit rij factors

A simplified analysis is presented for the evaluation of the three‐electron one‐center integrals of the form ∫rrrrrred r ₁d r ₂d r ₃, for the cases i, j, k, ≥−2, l=−2, m≥−1, n≥−1. These integrals arise in the calculation of lower bounds for energy levels and certain relativistic corrections to the energy when Hylleraas‐type basis sets are employed. Convergence accelerator techniques are employed to obtain a reasonable number of digits of precision, without excessive CPU requirements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 93–99, 1999     

Electron localization functions and local measures of the covariance

The electron localization measure proposed by Becke and Edgecombe is shown to be related to the covariance of the electron pair distribution. Just as with the electron localization function, the local covariance does not seem to be, in and of itself, a useful quantity for elucidating shell structure. A function of the local covariance, however, is useful for this purpose. A different function, based on the hyperbolic tangent, is proposed to elucidate the shell structure encapsulated by the local covariance; this function also seems to work better for the electron localization measure of Becke and Edgecombe. In addition, we propose a different measure for the electron localization that incorporates both the electron localization measure of Becke and Edgecombe and the Laplacian of the electron density; preliminary indications are that this measure is especially good at elucidating the shell structure in valence regions. Methods for evaluating electron localization functions directly from the electron density, without recourse to the Kohn-Sham orbitals, are discussed.     

Calculation of Two-Center Nuclear Attraction Integrals over Integer and Noninteger n-Slater Type Orbitals in Nonlined-Up Coordinate Systems

Two-center nuclear attraction integrals over Slater type orbitals with integer and noninteger principal quantum numbers in nonlined up coordinate systems have been calculated by means of formulas in our previous work (T. Özdoan and M. Orbay, Int. J. Quant. Chem. 87 (2002) 15). The computer results for integer case are in best agreement with the prior literature. On the other hand, the results for noninteger case are not compared with the literature due to the scarcity of the literature, but also compared with the limit of integer case and good agreements are obtained. The proposed algorithm for the calculation of two-center nuclear attraction integrals over Slater type orbitals with noninteger principal quantum numbers in nonlined-up coordinate systems permits to avoid the interpolation procedure used to overcome the difficulty introduced by the presence of noninteger principal quantum numbers. Finally, numerical aspects of the presented formulae are analyzed under wide range of quantum numbers, orbital exponents and internuclear distances.     

Bond functions and many-body effects of the helium trimer

A general scheme for efficient implementation of bond functions in homonuclear triatomics is suggested and applied to the linear and triangular configurations of the helium trimer. It is found that only one set of midbond functions of size 6s3p can provide nearly all of the benefits obtainable from larger sizes as well as 100% of the energy lowering obtained with ten sets of d-functions added at the atom centres. They also enhance the convergence properties of the many body terms at the Hartree-Fock and electron correlation levels. Correct dissociation limits and avoiding spurious minima of potential wells as well as other linear and triangular configurations are taken into account. Received: 5 November 1996 / Accepted: 7 April 1997     

High resolution electron momentum spectroscopy of the valence orbitals of water

The development of a third-generation electron momentum spectrometer with significantly improved energy and momentum resolutions at Tsinghua University (ΔE = 0.45–0.68 eV, Δθ = ±0.53° and Δ? = ±0.84°) has enabled a reinvestigation of the valence orbital electron momentum distributions of H₂O with improved statistical accuracy. The measurements have been conducted at impact energies of 1200 eV and 2400 eV in order to check the validity of the plane wave impulse approximation. The obtained ionization spectra and electron momentum distributions have been compared with the results of computations carried out with Hartree Fock [HF] theory, density functional theory in conjunction with the standard B3LYP functional, one-particle Green’s function [1p-GF] theory along with the third-order algebraic diagrammatic construction scheme [ADC(3)], symmetry adapted cluster configuration interaction [SAC-CI] theory, and a variety of multi-reference [MR-SDCI, MR-RSPT2, MR-RSPT3] theories. The influence of the basis set on the computed momentum distributions has been investigated further, using a variety of basis sets ranging from 6-31G to the almost complete d-aug-cc-pV6Z basis set. A main issue in the present work pertains to a shake-up band of very weak intensity at 27.1 eV, of which the related momentum distribution was analyzed for the first time. The experimental evidences and the most thorough theoretical calculations demonstrate that this band borrows its ionization intensity from the 2a₁ orbital.     

On the correlation between the bifurcation values of the electron localization function and the nuclear independent chemical shift in aromatic and heteroaromatic compounds

In this work, for a representative set of 24 aromatic molecules, that includes hydrocarbons, hetero and charged rings, we explore the correlations between two aromaticity indexes. One of these indexes is based on the bifurcation values of the σ and π electron localization function (ELF), while the other one is the nuclear independent chemical shift (NICS). We observe that the quality of the possible correlations between these two kinds of indexes strongly depends on the kind of rings and on the particular indexes that are considered.