Towards the analytic theory of Potential Energy Curves for diatomic molecules. Studying He${}_2^+$ and LiH diatomics as illustration
Following the first principles the elements of the analytic theory of potential curves for diatomic molecules (diatomics) are presented. It is based on matching the perturbation theory at small internuclear distances $R$ and multipole expansion at large distances, modified by instanton-induced trans-series for homonuclear case, with addition of the phenomenologically described equilibrium configuration, if exists. It leads to a new class of (generalized) meromorphic potentials (modified by exponential terms) with difference in degrees of polynomials in numerator and denominator equal to 4 (6) for positively charged (neutral) diatomics. As examples the He$_2^+$ and LiH diatomics in Born-Oppenheimer (adiabatic) approximation are considered. For ${}^{4}$He$_2^+$ (${}^{3}$He$_2^+$) diatomics it is found the approximate analytic expression for the potential energy curves (analytic PEC) $V(R)$ for the ground state $X^2 Σ_u^+$ and the first excited state $A^2 Σ_g^+$. It provides 3-4 s.d. correctly for distances $R \in [1, 10]$\,a.u. with some irregularities for $A^2 Σ_g^+$ PEC at small distances (much smaller than equilibrium distances) probably related to level crossings which may occur there. The analytic PEC for the ground state $X^2 Σ_u^+$ supports 829 (626) rovibrational states with 3-4 s.d. of accuracy in energy, which is only by 1 state less (more) than 830 (625) reported in the literature. In turn, the analytic PEC for the excited state $A^2 Σ_g^+$ supports all 9 reported weakly-bound rovibrational states. Entire rovibrational spectra is found in a single calculation using the code based on the Lagrange mesh method.