Journal of Chemistry

Journal of Chemistry / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 1598951 | https://doi.org/10.1155/2017/1598951

Ebru Çopuroğlu, "Evaluation of Self-Friction Three-Center Nuclear Attraction Integrals with Integer and Noninteger Principal Quantum Numbers over Slater Type Orbitals", Journal of Chemistry, vol. 2017, Article ID 1598951, 6 pages, 2017. https://doi.org/10.1155/2017/1598951

Evaluation of Self-Friction Three-Center Nuclear Attraction Integrals with Integer and Noninteger Principal Quantum Numbers over Slater Type Orbitals

Academic Editor: Demeter Tzeli
Received02 Nov 2016
Revised07 Feb 2017
Accepted21 Feb 2017
Published06 Mar 2017

Abstract

We have proposed a new approach to evaluate self-friction (SF) three-center nuclear attraction integrals over integer and noninteger Slater type orbitals (STOs) by using Guseinov one-range addition theorem in standard convention. A complete orthonormal set of Guseinov exponential type orbitals (-ETOs, ) has been used to obtain the analytical expressions. The overlap integrals with noninteger quantum numbers occurring in SF three-center nuclear attraction integrals have been evaluated using auxiliary functions. The accuracy of obtained formulas is satisfactory for arbitrary integer and noninteger principal quantum numbers.

1. Introduction

It is well known that Roothaan open-shell Hartree-Fock theory (HFR) and its extensions are not, in general, applicable to any state of a single configuration, which has any symmetry of open shells [110]. In a recent paper [11], Guseinov has eliminated the insufficiencies arising in HFR theory and has suggested a combined open-shell Hartree-Fock-Roothaan theory (CHFR) which includes the arbitrary number and symmetry of closed- and open-shell cases. It is clear that the evaluation of multicenter molecular integrals arising in CHFR theory has prime importance in many branches of physics particularly quantum physics, atomic and molecular physics, and nuclear physics. Using different approximations, considerable efforts have been made for calculating these integrals over STOs with noninteger principal quantum numbers (NISTOs). Notice that the NISTOs describe the physical systems more accurately than do ISTOs with integer principal quantum numbers [1217]. But the main problem of using NISTOs basis in molecular calculations arises in the evaluation of the multicenter integrals. Various methods have been proposed in the literature for improving evaluation for the use of NISTOs in molecular calculations. It should be noted that most of the studies in the literature on the evaluation of multicenter molecular integrals with NISTOs were not entirely successful [1820].

According to classical electrodynamics, the Lagrangian for a system of charges does not include the Lorentz self-friction (SF) potentials [2123]. By considering this deficiency, Guseinov had suggested a new centrally symmetric potential of the SF field and on the basis of this idea Guseinov proposed the new complete orthonormal sets of -ETOs. The indices arising from the use of total potential and occurring in the Guseinov -ETOs are the frictional quantum number. For the evaluation of multicenter molecular integrals of integer Slater type orbitals (ISTOs) and noninteger Slater type orbitals (NISTOs) appearing in the CHFR approximation, Guseinov derived one-range addition theorems by the use of complete orthonormal sets of Guseinov -ETOs [24].

In this case the analytical formulas for the evaluation of multicenter molecular integrals are directly depending on the SF quantum number. Thus, considering SF interactions, the calculation results of physical and chemical properties of atoms and molecules obtained from CHFR equations are closer to their real values. Therefore, the new type of multicenter molecular integrals arising in CHFR equations should be called SF multicenter molecular integrals. SF three-center nuclear attraction integral is one of the most important molecular integrals arising in CHFR equations.

The aim of this paper is to provide the general analytical expressions of SF three-center nuclear attraction integrals over ISTOs and NISTOs. With the help of Guseinov one-range addition theorems and -ETOs, we have established new algorithm for calculating SF three-center nuclear attraction integrals. We noticed that the SF three-center nuclear attraction integrals are expressed in terms of two-center overlap integrals and two-center nuclear attraction integrals over NISTOs. The convergence, accuracy, and CPU time have been tested by our previous studies.

2. Definitions and Basic Formulas

The ISTOs can be written as follows [25, 26]:where is an integer principal quantum number, is the screening constant, and is a complex or real spherical harmonic (see [2729] for exact definitions).

The NISTOs with noninteger values of the principal quantum number are defined byHere, is the gamma function.

To obtain the expansion of NISTO’s about a new center by the use of complete orthonormal sets of -ETOs we can use the method for the expansion of ISTOs. Finally the formulas for the expansion of NISTOs in terms of ISTOs at a new origin are defined as [24]where , , and the coefficients are determined as follows:In (3) and (4), and the quantities are the overlap integrals between the NISTOs and ISTOs:

For the evaluation of overlap integrals we use the formula defined in terms of auxiliary functions [30]. As can be seen from (3), the coefficients for expansion of NISTOs for are expressed through the noninteger overlap integrals with the same screening constants defined by (5). It is clear that the obtained expansion formulas from this paper can also be used in the case of integer values of . Equation (3) becomes the series expansion formulas for translation of ISTOs (see, (5) of [31]).

3. Evaluation of SF Three-Center Nuclear Attraction Integrals

The SF three-center nuclear attraction integrals over NISTOs are defined by where and .

Taking into account (3) in (6) we obtain the following relations:Here is the two-center nuclear attraction integral over NISTOs determined by [32]

In order to calculate the two-center nuclear attraction integrals we use the Laplace expansion of Coulomb potential. Then we find the following analytical expression for the two-center nuclear attraction integrals: here , , and . See Appendix for exact formulas of , and quantities and functions occurring in (9). The quantities and in (9) are well-known complete and incomplete gamma functions defined by, respectively, These functions have already been investigated by numerous authors with different algorithms (see, e.g., [33, 34] and references therein).

As can be seen from (7) and (9), by the use of Guseinov one-range addition theorems, the SF three-center nuclear attraction integrals are expressed through the two-center overlap integrals, gamma functions, auxiliary functions, and coefficients. The obtained formula is basic and valid for arbitrary integer and noninteger principal quantum numbers.

4. Numerical Results and Discussion

We have presented a new analytical method for the calculation of self-friction three-center nuclear attraction integrals over STOs and NISTOs with the help of Guseinov one-range addition theorems. The computer programs of established formulas have been performed by using Mathematica 7.0 programming language in Intel (R) Core (TM) 2.5 GHz computer. On the basis of (7), the calculation results of SF three-center nuclear attraction integrals over STOs and NISTOs have been given in Tables 1, 2, and 3. In Table 1, for and , the comparison results of SF three-center nuclear attraction integrals over STOs and NISTOs have been presented. As can be seen from Table 1 our calculation results are satisfactory. Also, for different values of self-friction quantum number , the convergence of the series expansion relations of SF three-center nuclear attraction integral over NISTOs as a function of summation limits of had been shown in Tables 2 and 3.


Eq. (7)



Eq. (7) for

50.336896384834373
60.339487682453034
70.340270532993467
80.340116042929251
90.339661315498115
100.339243752352591
110.338980006668039
120.338868607202027
130.338861299222622
140.338908803548408
150.338971794137852
160.339028513191853
170.339068958094987
180.339092542090420
190.339102386725724
200.339103160558990


Eq. (7) for

50.34020570374536
60.34023045324266
70.33972836080133
80.33926825347282
90.33899764955237
100.33890266956856
110.33891035305428
120.33896166770645
130.33901588628170
140.33905691933072
150.33908024198451
160.33909004101075
170.33909111299376
180.33908862095281
190.33908529189431
200.33908278090016

The coefficients , , and are repeatedly needed in a computation of the three-center nuclear attraction integrals. We believe that a common storage scheme for the and coefficients and overlap integrals with the same selection rule, as proposed in this study, will give important contributions in reducing requirements of computer time for computation of multicenter integrals which arise in the Hartree-Fock-Roothaan and Hylleraas approaches. For quick calculations, the and coefficients and integrals are stored in the memory of the computer. In order to put these coefficients and integrals into or to get them back from the memory, the positions of certain coefficients , , and integrals are determined by the following relations.For binomial coefficientsFor coefficientsFor overlap integrals

To demonstrate the effectiveness of storage method for given coefficients and integrals we calculated three-center nuclear attraction integral with noninteger quantum numbers. The computer time required for the calculation of three-center nuclear attraction integrals is not given in the tables due to the fact that the comparison cannot be made with the different computers used in the literature. For instance, for with sets , , , , , , , and CPU time takes about 5.897 ms and 2.57 ms by using usual and storage methods, respectively.

In Figure 1 we plot the convergence of the series (see (7)) of for with different values of parameter as a function of the indices . The series accuracy is shown in Figure 1 where the quantities are the values of integral for . As can be seen from Figure 1, for given , the convergence rates are satisfactory for arbitrary values of parameter . Greater accuracy is attainable by the use of more terms of expansions (7).

5. Conclusion

In conclusion, by considering self-friction effects, we have proposed an analytical efficient method for the calculation of three-center nuclear attraction integrals over STOs and NISTOs using Guseinov one-range addition theorems. This allows us to take into account the self-friction field effect in atomic and molecular electronic structure evaluations.

Appendix

The generalized Gaunt coefficients and occurring in (9) are determined by the following relationships [32]:where See [32] for the exact definition of quantities and .

In this work is the complex or real spherical harmonics. The definition of phases for the complex spherical harmonics differs from the Condon-Shortley phases [27] by the sign factor and can be defined aswhereHere and .

Competing Interests

The author declares that they have no competing interests.

References

  1. C. C. Roothaan, “Self-consistent field theory for open shells of electronic systems,” Reviews of Modern Physics, vol. 32, pp. 179–185, 1960. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  2. S. Huzinaga, “Coupling operator method in the SCF theory,” The Journal of Chemical Physics, vol. 51, no. 9, pp. 3971–3975, 1969. View at: Publisher Site | Google Scholar
  3. K. Hirao and H. Nakatsuji, “General SCF operator satisfying correct variational condition,” The Journal of Chemical Physics, vol. 59, no. 3, pp. 1457–1462, 1973. View at: Publisher Site | Google Scholar
  4. T. Kagawa, “Relativistic Hartree-Fock-Roothaan theory for open-shell atoms,” Physical Review A, vol. 12, no. 6, pp. 2245–2256, 1975. View at: Publisher Site | Google Scholar
  5. H. J. Silverstone, “Hartree-Fock equations for open shell states,” The Journal of Chemical Physics, vol. 67, no. 9, pp. 4172–4180, 1977. View at: Google Scholar
  6. R. Carbo and J. M. Riera, Lecture Notes in Chemistry: A General SCF Theory, vol. 5, Springer, Berlin, Germany, 1978.
  7. G. Fieck and J. Wirsich, “Electrostatic interaction for arbitrary point groups, configurations, and coupling schemes,” International Journal of Quantum Chemistry, vol. 18, no. 3, pp. 753–773, 1980. View at: Publisher Site | Google Scholar
  8. M. Haser, “Molecular point-group symmetry in electronic structure calculations,” The Journal of Chemical Physics, vol. 95, no. 11, pp. 8259–8265, 1991. View at: Publisher Site | Google Scholar | MathSciNet
  9. K. Koc and Y. Ishikawa, “Single-Fock-operator method for matrix Dirac-Fock self-consistent-field calculations on open-shell atoms,” Physical Review A, vol. 49, no. 2, pp. 794–798, 1994. View at: Publisher Site | Google Scholar
  10. P. M. Kozlowski and E. R. Davidson, “Construction of open shell perturbation theory invariant with respect to orbital degeneracy,” Chemical Physics Letters, vol. 226, no. 5-6, pp. 440–446, 1994. View at: Publisher Site | Google Scholar
  11. I. I. Guseinov, “Combined open shell Hartree-Fock theory of atomic-molecular and nuclear systems,” Journal of Mathematical Chemistry, vol. 42, no. 2, pp. 177–189, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  12. I. I. Guseinov, B. A. Mamedov, M. Ertürk, H. Aksu, and E. Şahin, “Application of combined Hartree-Fock-Roothaan theory to atoms with an arbitrary number of closed and open shells of any symmetry,” Few-Body Systems, vol. 41, no. 3-4, pp. 211–220, 2007. View at: Publisher Site | Google Scholar
  13. I. I. Guseinov, E. Sahin, and M. Erturk, “An improvement on-exponential type orbitals for atoms in standard convention,” Molecular Physics, vol. 112, no. 1, pp. 35–40, 2014. View at: Publisher Site | Google Scholar
  14. M. Erturk, “An improvement on noninteger n-generalized exponential type orbitals with hyperbolic cosine in atomic calculations,” Bulletin of the Chemical Society of Japan, vol. 85, no. 7, pp. 753–757, 2012. View at: Publisher Site | Google Scholar
  15. T. Koga and K. Kanayama, “Noninteger principal quantum numbers increase the efficiency of Slater-type basis sets: singly charged cations and anions,” Journal of Physics B: Atomic, Molecular and Optical Physics, vol. 30, no. 7, p. 1623, 1997. View at: Publisher Site | Google Scholar
  16. T. Koga and K. Kanayama, “Noninteger principal quantum numbers increase the efficiency of Slater-type basis sets: heavy atoms,” Chemical Physics Letters, vol. 266, no. 1-2, pp. 123–129, 1997. View at: Publisher Site | Google Scholar
  17. T. Koga, T. Shimazaki, and T. Satoh, “Noninteger principal quantum numbers increase the efficiency of Slater- type basis sets: double-zeta approximation,” Journal of Molecular Structure: THEOCHEM, vol. 496, pp. 95–100, 2000. View at: Publisher Site | Google Scholar
  18. K. G. Kay and H. J. Silverstone, “Analytical evaluation of multicenter integrals of r 12− 1 with Slater‐type atomic orbitals. V. four‐center integrals by Fourier‐transform method,” The Journal of Chemical Physics, vol. 51, pp. 4287–4304, 1969. View at: Publisher Site | Google Scholar | MathSciNet
  19. E. O. Steinborn and E. Filter, “Translations of fields represented by spherical-harmonic expansions for molecular calculations,” Theoretica Chimica Acta, vol. 38, no. 4, pp. 247–260, 1975. View at: Publisher Site | Google Scholar | MathSciNet
  20. E. Filter and E. O. Steinborn, “A matrix representation of the translation operator with respect to a basis set of exponentially declining functions,” Journal of Mathematical Physics, vol. 21, no. 12, pp. 2725–2736, 1980. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  21. H. A. Lorentz, The Theory of Electrons, Dover, New York, NY, USA, 1953.
  22. W. Heitler, The Quantum Theory of Radiation, Oxford University Press, Oxford, UK, 1950.
  23. L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, Pergamon, New York, NY, USA, 1987.
  24. I. I. Guseinov, “Unified treatment of one-range addition theorems for complete orthonormal sets of generalized exponential-type orbitals and noninteger N slater functions,” Bulletin of the Chemical Society of Japan, vol. 87, no. 10, pp. 1101–1103, 2014. View at: Publisher Site | Google Scholar
  25. P.-O. Löwdin, “Expansion theorems for the total wave function and extended Hartree-Fock schemes,” Reviews of Modern Physics, vol. 32, pp. 328–334, 1960. View at: Publisher Site | Google Scholar | MathSciNet
  26. J. C. Slater, Quantum Theory of Atomic Structure, vol. 2, McGraw-Hill, New York, NY, USA, 1960.
  27. E. U. Condon and G. H. Shortley, Theory of Atomic Spectra, Cambridge University Press, London, UK, 1963.
  28. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, New York, NY, USA, 1980. View at: MathSciNet
  29. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 1972.
  30. I. I. Guseinov and B. A. Mamedov, “Evaluation of overlap integrals with integer and noninteger n Slater-type orbitals using auxiliary functions,” Journal of Molecular Modeling, vol. 8, no. 9, pp. 272–276, 2002. View at: Publisher Site | Google Scholar
  31. I. I. Guseinov, “One-range addition theorems for complete sets of modified exponential type orbitals and noninteger n slater functions in standard convention,” Journal of the Chinese Chemical Society, vol. 61, no. 4, pp. 477–480, 2014. View at: Publisher Site | Google Scholar
  32. I. I. Guseinov, “Analytical evaluation of two-centre Coulomb, hybrid and one-electron integrals for Slater-type orbitals,” Journal of Physics B: Atomic and Molecular Physics, vol. 3, no. 11, pp. 1399–1412, 1970. View at: Publisher Site | Google Scholar
  33. E. Çopuroğlu and T. Mehmetoğlu, “Full analytical evaluation of the einstein relation for disordered semiconductors,” IEEE Transactions on Electron Devices, vol. 62, no. 5, pp. 1580–1583, 2015. View at: Publisher Site | Google Scholar
  34. I. I. Guseinov and B. A. Mamedov, “Evaluation of incomplete gamma functions using downward recursion and analytical relations,” Journal of Mathematical Chemistry, vol. 36, no. 4, pp. 341–346, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

Copyright © 2017 Ebru Çopuroğlu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views585
Downloads275
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.