Abstract
For general two-electron two-centre integrals over Slater-type orbitals (STOs), the use of the Neumann expansion for the Coulomb interaction potential yields infinite series in terms of few basic functions. In many important cases the number of terms necessary to achieve convergence by a straightforward summation is large and one is forced to calculate the basic integrals of high order. We present a systematic approach to calculation of the higher-order terms in the Neumann series by large-order expansions of the basic integrals. The final expressions are shown to be transparent and straightforward to implement, and all auxiliary quantities can be calculated analytically. Moreover, numerical stability and computational efficiency are also discussed. Results of the present work can be used to speed up calculations of the STOs integral files, but also to study convergence of the Neumann expansion and develop appropriate convergence accelerators.
Similar content being viewed by others
References
J.C. Slater, Phys. Rev. 36, 57 (1930)
J.C. Slater, Phys. Rev. 42, 33 (1932)
S.F. Boys, Proc. Roy. Soc. A200, 542 (1950). (London)
T. Kato, Commun. Pure Appl. Math. 10, 151 (1957)
S. Agmon, Lectures on Exponential Decay of Solutions of Second-Order Elliptic Equations: Bound on Eigenfunctions of N-Body Schrödinger Operators (Princeton University Press, Princeton, 1982)
H. Jones, Int. J. Quantum Chem. 51, 417 (1994)
M.P. Barnett, Int. J. Quantum Chem. 76, 464 (2000)
M.P. Barnett, J. Chem. Phys. 113, 9419 (2000)
L. Berlu, H. Safouhi, P. Hoggan, Int. J. Quantum Chem. 99, 221 (2004)
H. Safouhi, L. Berlu, J. Comput. Phys. 216, 19 (2006)
J. Fernández Rico, J.J. Fernández, I. Ema, R. López, G. Ramírez, Int. J. Quantum Chem. 78, 83 (2000)
J. Fernández Rico, J.J. Fernández, R. López, G. Ramírez, Int. J. Quantum Chem. 78, 137 (2000)
J. Fernández Rico, J.J. Fernández, I. Ema, R. López, G. Ramírez, Int. J. Quantum Chem. 81, 16 (2001)
J. Fernández Rico, R. López, I. Ema, G. Ramírez, J. Comput. Chem. 19, 1284 (1998)
P.E. Hoggan, Int. J. Quantum Chem. 109, 2926 (2009)
P.E. Hoggan, Int. J. Quantum Chem. 110, 98 (2010)
A. Bouferguene, P. E. Hoggan, QCPE, programme no. 667 (1996)
K. Pachucki, Phys. Rev. A 80, 032520 (2009)
K. Pachucki, Phys. Rev. A 85, 042511 (2012)
M. Lesiuk, R. Moszynski, Phys. Rev. A 86, 052513 (2012)
E.N. Maslen, M.G. Trefry, Int. J. Quantum Chem. 37, 51 (1990)
F.E. Harris, Int. J. Quantum Chem. 88, 701 (2002)
M. Belén Ruiz, J. Math. Chem. 43, 701 (2008)
M. Belén Ruiz, J. Math. Chem. 46, 24 (2009)
M. Belén Ruiz, J. Math. Chem. 46, 1322 (2009)
P.E. Hoggan, M. Belén Ruiz, T. Ozdogan, Molecular integrals over slater-type orbitals. From pioneers to recent progress, in Quantum Frontiers of Atoms and Molecules in Physics, Chemistry, and Biology, 1st edn., ed. by Mihai V. Putz (Nova Science Publishers, Inc., New York, 2010), pp. 63–90
M. Belén Ruiz, J. Math. Chem. 49, 2457 (2011)
M. Lesiuk, R. Moszynski, Phys. Rev. E 90, 063318 (2014)
M. Lesiuk, R. Moszynski, Phys. Rev. E 90, 063319 (2014)
M. Lesiuk, M. Przybytek, M. Musial, B. Jeziorski, R. Moszynski, Phys. Rev. A 91, 012510 (2015)
K. Ruedenberg, J. Chem. Phys. 19, 1459 (1951)
K. Ruedenberg, C.C.J. Roothaan, W. Jaunzemis, J. Chem. Phys. 24, 201 (1956)
A.C. Wahl, P.E. Cade, C.C.J. Roothaan, J. Chem. Phys. 41, 2578 (1964)
E.L. Mehler, K. Ruedenberg, J. Chem. Phys. 50, 2578 (1964)
M. Kotani, A. Amemiya, E. Ishiguro, T. Kimura, Table of Molecular Integrals (Maruzen, Tokyo, 1955)
F.E. Harris, J. Chem. Phys. 32, 3 (1960)
J. Yasui, A. Saika, J. Chem. Phys. 76, 468 (1982)
J. Fernández Rico, R. López, G. Ramírez, J. Chem. Phys. 91, 4213 (1989)
J. Fernández Rico, R. López, G. Ramírez, J. Chem. Phys. 97, 7613 (1992)
J. Fernández Rico, R. López, G. Ramírez, C. Tablero, Phys. Rev. A 49, 3381 (1994)
R. López, G. Ramírez, Int. J. Quantum Chem. 49, 11 (1994)
J. Fernández Rico, R. López, G. Ramírez, J. Chem. Phys. 101, 9807 (1994)
M. Abramowitz, I. Stegun (eds.), Handbook of Mathematical Functions (Dover, New York, 1972)
F.W.J. Olver, Asymptotics and Special Functions (A. K. Peters, Wellesley, 1997)
D.S. Jones, Math. Methods Appl. Sci. 24, 369 (2001)
T.M. Dunster, Stud. Appl. Math. 113, 245 (2004)
A. Sidi, P.E. Hoggan, Int. J. Pure Appl. Math. 71, 481 (2011)
P. Hoggan, A. Sidi, D. Pinchon, Adv. Quantum Chem. 68, 43 (2014)
W. Gautschi, SIAM Rev. 9, 24 (1967)
Acknowledgments
This work was supported by the Polish Ministry of Science and Higher Education through the project “Diamentowy Grant”, Number DI2011 012041.
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix 1: Auxiliary integrals
Virtually all final working formulae obtained in the present paper are given in terms of the basic integrals \(E_n(z)\) and \(a_n(z)\). They are defined through the integral representations
where n is an arbitrary integer in the former and a nonnegative integer in the latter. Calculation of \(a_n\) is most easily carried out with help of the Miller algorithm [49] as discussed by Harris [22]. The integral \(E_n\) is usually called the generalised exponential integral. Computation of \(E_n\) differs depending on the sign of n. For a negative integer n the following recursion is completely stable in the upward direction
For positive n and \(z<1\) one uses the series expansion
where \(\varPsi (n)\) is the digamma function at integer argument. The above infinite summations converge to the machine precision in, at most, few tens of terms. Finally, for positive n and \(z>1\) the continued fraction (CF) formula can be applied
To evaluate the CF one can use the Lentz algorithm [43]. The only inconvenience is that consecutive numerators and denominators in the Lentz scheme grow very quickly with the number of terms retained in Eq. (65). Therefore, it is necessary to rescale them from time to time by a small number to avoid numerical overflows. Let us also recall the leading terms of the large-order asymptotic expansions for \(E_n(z)\) and \(a_n(z)\) which read
where \(n>z\).
Appendix 2: Large-order asymptotic formulae for \(i_\mu (\alpha )\) and \(k_\mu (\alpha )\)
According to the work of Sidi et al. [47, 48] the modified spherical Bessel functions posses the following large-order expansions
as \(\mu \rightarrow \infty \) at a fixed z, where
for \(m>0\) and \(b_0(z)=1\). The quantities \(S_{mk}\) in the above expression are the Stirling numbers of the second kind [43] defined recursively as
and we additionally adapt the convention \(S_{mk}=0\) for \(m<k\) or \(m<0\).
Rights and permissions
About this article
Cite this article
Lesiuk, M. Calculation of STOs electron repulsion integrals by ellipsoidal expansion and large-order approximations. J Math Chem 54, 572–591 (2016). https://doi.org/10.1007/s10910-015-0576-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10910-015-0576-5