Lah Number Mediated Analytical Two‐center Two‐electron Integrals of Coulomb Green Function of H‐like Orbitals for ¹Σ States

Abstract

Energetics of two‐center two‐electron (2c–2e) systems carry challenges in theoretical understanding of Schrödinger equation (SE) for well‐known divergence of Coulomb interactions and nuclear separation (R) in modified H‐like AOs, Slater type orbitals (STOs), Gaussian type orbitals (GTOs), B‐spline, Sturmian function and etc. employed to VBT and MOT. Certain elegant computational and analytical techniques were developed for STO, GTO and other square integrable basis set within Born‐Oppenheimer (BO) approximation. STOs and GTOs have an essential limitation of absence of radial nodes. Thus, analytical treatment has become an urge for H‐like AOs. We have considered the diatomic molecules only for the sake of simplicity. Employing Sheffer identity in associated Laguerre polynomial/Whittaker‐M function forms of H‐like AOs and transforming integrals into elliptic coordinates with two nuclei on two foci furnishes exact, analytical and simple Coulomb integrals (Js) in terms of R. Lah number originated from for nuclear coordinates only due to Sheffer identity shows that energetics of diatomic molecules can be anticipated as extremum function of R. Therefore, the optimization of potential energy surface (PES) of electrons as gradient of R may lead to σ‐bond formation. In this paper, we have developed diagonal Js for bound states of H₂ molecule.