**Journal of Applied Mathematics and Physics**

Vol.08 No.01(2020), Article ID:97634,15 pages

10.4236/jamp.2020.81007

Modified Atomic Orbital Calculations of Energy of the (2s^{2} ^{1}S) Ground-State, the (2p^{2} ^{1}D); (3d^{2} ^{1}G) and (4f^{2} ^{1}I) Doubly Excited States of Helium Isoelectronic Sequence from H^{−} to Ca^{18+}

Malick Sow^{1*}, Ibrahima Sakho^{2}, Boubacar Sow^{1}, Abdou Diouf^{1}, Youssou Gning^{1}, Babou Diop^{1}, Matabara Dieng^{3}, Abdourahmane Diallo^{4}, Mamadou Diouldé Ba^{4}, Jean Kouhissoré Badiane^{4}, Mamadi Biaye^{5}, Ahmadou Wagué^{1}^{ }

^{1}Department of Physics, Atoms Lasers Laboratory, Faculty of Sciences and Technologies, University Cheikh Anta Diop, Dakar, Senegal

^{2}Department of Physics, UFR Sciences et Technologies, University of Thiès, Thiès, Senegal

^{3}UFR SATIC, University Alioune DIOP de Bambey, Bambey, Senegal

^{4}Department of Physics, UFR of Sciences and Technologies, University Assane Seck of Ziguinchor, Ziguinchor, Senegal

^{5}Department of Physics and Chemistry, Faculty of Sciences and Technologies of Training and Education, University Cheikh Anta Diop, Dakar, Senegal

Copyright © 2020 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 3, 2019; Accepted: January 3, 2020; Published: January 6, 2020

ABSTRACT

We report in this paper the ground-state energy 2s^{2} ^{1}S and total energies of doubly excited states 2p^{2} ^{1}D, 3d^{2} ^{1}G, 4f^{2} ^{1}I of the Helium isoelectronic sequence from H^{−} to Ca^{18+}. Calculations are performed using the Modified Atomic Orbital Theory (MAOT) in the framework of a variational procedure. The purpose of this study required a mathematical development of the Hamiltonian applied to Slater-type wave function [1] combining with Hylleraas-type wave function [2]. The study leads to analytical expressions which are carried out under special MAXIMA computational program. This first proposed MAOT variational procedure, leads to accurate results in good agreement as well as with available other theoretical results than experimental data. In the present work, a new correlated wave function is presented to express analytically the total energies for the 2s^{2} ^{1}S ground state and each doubly 2p^{2} ^{1}D, 3d^{2} ^{1}G, 4f^{2} ^{1}I excited states in the He-like systems.The present accurate data may be a useful guideline for future experimental and theoretical studies in the (nl^{2}) systems.

**Keywords:**

Modified Atomic Orbital Theory, Variational Calculations, Correlated Wavefunction, Energies, Ground-State, Doubly Excited States, Helium Isoelectronic Sequence, Atoms and Ions

1. Introduction

The resolution of the Schrödinger equation gives for the ground-state of the Helium atom the value E = 108.8 eV. The experimental result being equal to E_{0} = −79.0 eV, one conceives that one takes into account the electronic correlation term, overestimating, numerous studies provided evidence of the importance of the electronic correlation in the ground-state and in the doubly excited state of He-like series [3] [4] [5].

For the ground states, energy calculations of Helium isoelectronic sequence are performed by using several correlated wavefunction and an analytical technique calculation. Thus, Moore [6], Radzig and Smirnov [7] and Arnaud [8] measured the ground state energies E (1s^{2} ^{1}S) of two-electron systems with atomic number Z (2 ≤ Z ≤ 10). Combining the perturbation theory [9] to the Ritz variational method [10], Sakho et al., [11] set in motion a technic of analytical calculation of the ground state energy E (^{1}S_{0}), the first ionization energy J (^{1}S_{0}) and the radial correlation expectation value
$\langle {r}_{12}^{-1}\rangle $ (^{1}S_{0}) for the Helium-like ions from Hydrogen ion H^{−} (Z = 1) to silicon Si^{12+} (Z = 14). Utpal and Talukdar [12] used an analytical approach to also calculate the ground-state energies of helium isoelectronic sequence from Hydrogen ion H^{−} (Z = 1) up to the silicon ion Si^{12+} (Z = 14). For some methods using variational wavefunction like Hylleraas’s one, a good approximation of the eigenvalues is obtained when the minima of the

function ( $\frac{\text{d}E}{\text{d}\alpha}$, with α a variational parameter, E is the energy) converges with

the increasing values of the dimension D and when the function $E=f\left(\alpha ,D\right)$ exhibits a plateau [2]. For other methods using a non variational wavefunction [13], the ground-state energy of He-like ions is determined by the use of proper core boundary conditions correct behavior for ${r}_{12}\to 0$ and ${r}_{12}\to \infty $ and by taking recourse to a perturbative method. For some methods using an analytical technique calculation, some authors were interested in the setting in work of calculation techniques permitting one to succeed to an analytical expression of the ground-state energy of the He-like ions. Thus, an analytical calculation for the ground-state energy and radial expectation values of Helium isoelectronic sequence is managed by using a wavefunction of the type of Bhattacharyya [12] [14]. Besides, developing the orbital atomic theory, Slater in [15] introduced the notions of screen constant σ and effective quantum number n* for the calculation of the energy of an electronic configuration given containing N electrons.

On the basis of his theory, Slater expresses analytically the total energy of an atomic system of N electrons according to σ and n^{*} determined from rules that he established. The analytical formula of Slater permits the simple calculation of the ground-state energy of He-like ions for which N = 2, σ = 0.30 and n^{*} = 1.

For the doubly excited states in He-like ions, since the early experiment [16] and theoretical explanation [17], Doubly Excited States (DES) of Helium isoelectronic sequence have been the target of a number of theoretical approaches. Greatest attention have been concentrated on studying symmetric DES (nl^{2}) with excited electrons having equal values of principal quantum number n (intrashell states) where the electronic correlation effect may be predominant [18]. The investigations of the intrashell states of two-electrons systems are advanced and due to the group theoretical method [19] which allowed intrashell states to be approximatively classified and some of these properties studied [20]. Theoretical investigations of (nl^{2}) doubly excited states are performed by using various method. The projection operator method and group theoretical methods [20] have been used for energies calculations of the 2s^{2}, 2p^{2} states in Helium-like ions. Time independent variational perturbation [21] was applied for total energies calculations of the 2s^{2}, 2p^{2} and 3d^{2} states in He, Li^{+}, Be^{2+} and B^{3+}. The correlation part of the energies for the 2s^{2}, 2p^{2}, 3s^{2}, 3p^{2}, 3d^{2} states in He isoelectronic series have been investigated by using perturbation theory [4]. The Screening Constant by Unit Nuclear charge (SCUNC) method [22] used a semi-empirical procedure to calculate (ns^{2}) ^{1}S^{e}, (np^{2}) ^{1}D^{e} and (Nsnp) ^{1}P° excited state of He-like ions. Recently the Modified Atomic Orbital Theory (MAOT) has been applied successfully in the studies of high lying ^{1,3}P° of He-like ions [1]. In this paper, we apply the first MAOT variational procedure, to calculate the ground-state energy 2s^{2} ^{1}S and the total energies of the singlet DES 2p^{2} ^{1}D, 3d^{2} ^{1}G and 4f^{2} ^{1}I of He isoelectronic sequence from H^{−} to Ca^{18+}. In addition, for the first time in our knowledge, we have also calculated theoretical the screening constant
${\sigma}_{theo}$ which is compared with experimental Slater screening constant (
${\sigma}_{exp}$ ) determined from his rule. Our present procedure leads to analytical expression which are carried out under MAXIMA computational program. Our energies positions are compared to other available theoretical and experimental data.

2. Theory

2.1. Brief Description of the MAOT Formalism

In the framework of Modified Atomic Orbital Theory (MAOT), total energy of (νℓ)-given orbital is expressed in the form [1] [23] :

$E\left(\upsilon \mathcal{l}\right)=-\frac{{\left[Z-\sigma \left(\mathcal{l}\right)\right]}^{2}}{{\upsilon}^{2}}$ (1)

For an atomic system of several electrons M, the total energy is given by (in Rydberg):

$E=-{\displaystyle \underset{i=1}{\overset{M}{\sum}}\frac{{\left[Z-{\sigma}_{i}\left(\mathcal{l}\right)\right]}^{2}}{{\upsilon}_{i}^{2}}}$

With respect to the usual spectroscopic notation $\left(N\mathcal{l},N{\mathcal{l}}^{\prime}\right){}^{2S+1}L{}^{\pi}$, this equation becomes

$E=-{\displaystyle \underset{i=1}{\overset{M}{\sum}}\frac{{\left[Z-{\sigma}_{i}\left({}^{2S+1}L{}^{\pi}\right)\right]}^{2}}{{\upsilon}_{i}^{2}}}$ (2)

In the photoionization study, energy resonances are generally measured relatively to the E_{∞} converging limit of a given (^{2S+1}L_{J})nl-Rydberg series. For these states, the general expression of the energy resonances is given by the formula [24] presented previously (in Rydberg units):

$\begin{array}{c}{E}_{n}={E}_{\infty}-\frac{1}{{n}^{2}}\{Z-{\sigma}_{1}\left({}^{2S+1}{L}_{J}\right)-{\sigma}_{2}\left({}^{2S+1}{L}_{J}\right)\times \frac{1}{n}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{{\sigma}_{2}^{\alpha}\left({}^{2S+1}{L}_{J}\right)\times \left(n-m\right)\times \left(n-q\right){\displaystyle \underset{k}{\sum}\frac{1}{{f}_{k}\left(n,m,q,s\right)}}\}}^{2}\end{array}$ (3)

In this equation m and q (m < q) denote the principal quantum numbers of the (^{2S+1}L_{J})nl-Rydberg series of the considered atomic system used in the empirical determination of the σ_{i}(^{2S+1}L_{J})-screening constants, s represents the spin of the nl-electron (s = 1/2), E_{∞} is the energy value of the series limit generally determined from the NIST atomic database, E_{n} denotes the corresponding energy resonance, and Z represents the nuclear charge of the considered element. The only problem that one may face by using the MAOT formalism is linked to the

determination of the $\underset{k}{\sum}\frac{1}{{f}_{k}\left(n,m,q,s\right)}$ term. The correct expression of this term

is determined iteratively by imposing general Equation (3) to give accurate data with a constant quantum defect values along all the considered series. The value of α in the ${\sigma}_{2}$ of the last term is fixed to 1 and 2 during the iteration. The quantum defect $\delta $ is calculated from the standard formula

${E}_{n}={E}_{\infty}-\frac{R{Z}_{core}^{2}}{{\left(n-\delta \right)}^{2}}\Rightarrow \delta =n-{Z}_{core}\sqrt{\frac{R}{\left({E}_{\infty}-{E}_{n}\right)}}$ (3 bis)

In this Equation (3 bis), R is the Rydberg constant, Z_{core} represents the electric charge of the core ion.

Z_{core} is directly obtained by the photoionization process from an atomic X system
$X+h\nu \to {X}^{p+}+p{e}^{-}$

${f}_{k}={f}_{k}\left({}^{2S+1}L{}_{J},n,s,m,q\right)$ are screening constants to be evaluated empirically with k taking the values from 1 to q.

L_{J}: denote the considered quantum state (S, P, D, F, ...).

2.2. Variational Procedure of Calculations

For the 2s^{2} ^{1}S ground state and each doubly excited states 2p^{2} ^{1}D, 3d^{2} ^{1}G, 4f^{2} ^{1}I, we constructed the basis wave functions below by combining Slater-type wave function [1] and Hylleraas-type wave functions [2] :

${\varphi}_{jkmn}\left({r}_{1},{r}_{2}\right)={\left({r}_{1}{r}_{2}\right)}^{\left(\upsilon -1\right)}\times \mathrm{exp}\left(-\left(\frac{Z-\sigma}{\upsilon \times {a}_{0}}\right)\left({r}_{1}+{r}_{2}\right)\right){\left({r}_{1}+{r}_{2}\right)}^{j}{\left({r}_{1}-{r}_{2}\right)}^{k}{\left|{r}_{1}-{r}_{2}\right|}^{m}$ (4)

We considered parameter

$\xi =\frac{Z-\sigma}{\upsilon \times {a}_{0}}$ (5)

and $Z,\sigma ,\upsilon $ and ${a}_{0}$ are respectively the nucleus charge number, the screening constant, the principal quantum number and Bohr’s radius

${\varphi}_{jkmn}\left({r}_{1},{r}_{2}\right)={\left({r}_{1}{r}_{2}\right)}^{\left(\upsilon -1\right)}\times \mathrm{exp}\left(-\xi \left({r}_{1}+{r}_{2}\right)\right){\left({r}_{1}+{r}_{2}\right)}^{J}{\left({r}_{1}-{r}_{2}\right)}^{K}{\left|{r}_{1}-{r}_{2}\right|}^{M}$ (6)

With $\left|{r}_{1}-{r}_{2}\right|=\sqrt{{r}_{1}^{2}+{r}_{2}^{2}-2{r}_{1}{r}_{2}\mathrm{cos}{\theta}_{12}}$

where ${r}_{1}$ and ${r}_{2}$ denote the positions of the two electrons;

${r}_{1}$ and ${r}_{2}$ are respectively used for $\left|{r}_{1}\right|$ and $\left|{r}_{2}\right|$,$J,K,M$ are Hylleraas parameters with $\left(J,K,M\ge 0\right)$.

J takes into account the distance of the two electrons from the nucleus.

K takes into account the approximation of the two electrons from the nucleus.

M takes into account the distance between the two electrons.

From Equation (5) the screening constant can be expressed

$\sigma =Z-{a}_{0}\upsilon \xi $ (7)

The final form of the wave function of the singlet doubly excited state can be written as follow:

${\Psi}_{n}\left({r}_{1},{r}_{2}\right)={\displaystyle \underset{jkm}{\sum}{\beta}_{jkm}{\varphi}_{jkmn}\left({r}_{1},{r}_{2}\right)}$ (8)

where the coefficients ${\beta}_{jkm}$ are determined by solving the Schrödinger equation:

$H{\Psi}_{n}\left({r}_{1},{r}_{2}\right)=E{\Psi}_{n}\left({r}_{1},{r}_{2}\right)$ (9)

where the Hamiltonian operator H has the form:

$H=T+C+W$ with (10)

$T=\frac{-{\hslash}^{2}}{2m}\left({\nabla}_{1}+{\nabla}_{2}\right)$ ; $C=-\left(\frac{Z{e}^{2}}{{r}_{1}}+\frac{Z{e}^{2}}{{r}_{2}}\right)$ ; $W=\frac{{e}^{2}}{\left|{r}_{1}-{r}_{2}\right|}$ (11)

where T is the kinetic energy, C is the Coulomb potential between the atomic nucleus and the two electrons, W is the Coulomb interaction between electrons.

In the Hamilton operator we neglected all magnetic and relativistic effects together with the motion of the atomic nucleus.

In this Equation (11), Z is the nuclear charge ∇_{1} is the Laplacian with reference to the coordinates of the vector radius
${r}_{1}$ which detect the position of the electron 1. ∇_{2} Laplacian defines the coordinates of the vector radius
${r}_{2}$ which detect the position of the electron 2 and
$\left|{r}_{1}-{r}_{2}\right|$ inter-electronic distance.

The representation of the Schrödinger equation on the non-orthogonal basis leads to the general eigenvalue equation [25] ;

$\underset{\left(j,k,m\right),\left({j}^{\prime}{k}^{\prime}{m}^{\prime}\right)}{\sum}\left({H}_{JKMnl}-E{N}_{JKMnl}\right)}=0$ (12)

With $J=j+{j}^{\prime}$ ; $K=k+{k}^{\prime}$ ; $M=m+{m}^{\prime}$

${N}_{JKMnl}=\langle {\varphi}_{jkmnl}|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}nl}\rangle $ : is the normalization factor (13)

${H}_{JKMnl}=\langle {\varphi}_{jkmnl}\left|H\right|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}nl}\rangle $ : is the matrix elements of Hamilton operator (14)

E: is the eigenvalue of the energy

${H}_{JKMnl}={T}_{jkm,{j}^{\prime}{k}^{\prime}{m}^{\prime}nl}+{C}_{JKMnl}+{W}_{JKMnl}$ (15)

${T}_{jkm,{j}^{\prime}{k}^{\prime}{m}^{\prime}nl}$ : is the matrix elements of the kinetic Energy operator of the two electrons

${C}_{JKMnl}$ : is the matrix elements of the Coulombian interaction Energy operator between the nucleus and the two electrons

${W}_{JKMnl}$ : is the matrix elements of the Coulombian interaction Energy operator between the two electrons

Thus following the form of the basis wave function above (4), we have constructed for each state a special wave function and then calculated the matrix elements ${N}_{JKMnl}$,${C}_{JKMnl}$,${W}_{JKMnl}$,${T}_{jkm,{j}^{\prime}{k}^{\prime}{m}^{\prime}nl}$

For example on 2p^{2} states:
$n=\upsilon =2;l=1$

The wave function is written as follow:

${\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)=\left({r}_{1}{r}_{2}\right)\times \mathrm{exp}\left(-\xi \left({r}_{1}+{r}_{2}\right)\right){\left({r}_{1}+{r}_{2}\right)}^{J}{\left({r}_{1}-{r}_{2}\right)}^{K}{\left|{r}_{1}-{r}_{2}\right|}^{M}$ (16)

The matrix elements of the normalization factor is written as follow:

${N}_{JKM21}=\langle {\varphi}_{jkm21}|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\rangle $ (17)

${N}_{JKM21}={\displaystyle \iiint {\text{d}}^{3}{r}_{1}}{\text{d}}^{3}{r}_{2}{\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)\times {\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)$ (18)

${N}_{JKM21}={{\displaystyle \iiint {\text{d}}^{3}{r}_{1}{\text{d}}^{3}{r}_{2}\left({r}_{1}{r}_{2}\right)}}^{2}{\left({r}_{1}+{r}_{2}\right)}^{J}{\left({r}_{1}-{r}_{2}\right)}^{K}{\left|{r}_{1}-{r}_{2}\right|}^{M}\times \mathrm{exp}\left(-2\xi \left({r}_{1}+{r}_{2}\right)\right)$ (19)

With $\iiint {\text{d}}^{3}{r}_{i}}={\displaystyle \int {r}_{i}^{2}\text{d}{r}_{i}{\displaystyle \underset{0}{\overset{\pi}{\int}}\mathrm{sin}{\theta}_{i}\text{d}{\theta}_{i}{\displaystyle \underset{0}{\overset{2\pi}{\int}}\text{d}{\phi}_{i}}}}=4\pi {\displaystyle \int {r}_{i}^{2}\text{d}{r}_{i}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(i=1,2\right)$ (20)

The matrix elements of the Coulombian interaction Energy operator between the nucleus and the two electrons is written as follow:

${C}_{JKMnl}=\langle {\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)\left|C\right|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)\rangle $ (21)

${C}_{JKM21}=\langle {\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)|-Z{e}^{2}\left(\frac{1}{{r}_{1}}+\frac{1}{{r}_{2}}\right)|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)\rangle $ (22)

${C}_{JKM21}=-Z{e}^{2}{\displaystyle \iiint {\text{d}}^{3}{r}_{1}}{\text{d}}^{3}{r}_{2}{\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)\left(\frac{1}{{r}_{1}}+\frac{1}{{r}_{2}}\right){\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)$ (23)

The matrix elements of the Coulombian interaction Energy operator between the two electrons is expressed as follow:

${W}_{JKM21}=\langle {\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)|\frac{1}{\left|{r}_{1}-{r}_{2}\right|}|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)\rangle $ (24)

${W}_{JKM21}={\displaystyle \iiint {\text{d}}^{3}{r}_{1}}{\text{d}}^{3}{r}_{2}{\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)\left(\frac{1}{\left|{r}_{1}-{r}_{2}\right|}\right){\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)$ (25)

The matrix elements of the kinetic Energy operator of the two electrons is expressed as follow:

${T}_{jkm,{j}^{\prime}{k}^{\prime}{m}^{\prime}21}=\langle {\varphi}_{jkm21}\left({r}_{1},{r}_{2}\right)\left|\frac{-{\hslash}^{2}}{2m}\left({\nabla}_{1}+{\nabla}_{2}\right)\right|{\varphi}_{{j}^{\prime}{k}^{\prime}{m}^{\prime}21}\left({r}_{1},{r}_{2}\right)\rangle $ (26)

Calculations for the other matrix elements of other states were obtained in the same procedure.

The mathematical development of the Hamiltonian applied to each wave function of each states leads to simplified analytical expressions which are carried out under MAXIMA computational program.

Concerning the screening constant σ and the variational parameter ξ, the procedure is as follow:

From the Slater condition $0.3\le \sigma \le 1$, and taken into account Equation (5), the parameter ξ can be expressed as follow:

$\frac{Z-1}{n}\le \xi \le \frac{Z}{n}$ (27)

In order to obtain the theoretical screening constant ${\sigma}_{theo}$, the exponential parameter ξ and the minimum eigenvalue in which we are interested, the analytical expressions of each matrix elements of each state are carried out by our self MAXIMA computational program.

MAXIMA is a computer algebra system for the manipulation of symbolic and numerical expressions, including differentiation, integration, ordinary differential equations, and matrix elements. MAXIMA yields high precision numeric results by using exact fractions, arbitrary precision integers, and variable precision floating point numbers. Our MAXIMA source code is compiled on windows systems.

At the beginning, the variational parameter ξ is determined as shown (27). For each value of Z and n, we know a limited area of the variational parameter ξ. Thus for each value of Z, n and ξ, the program calculated directly the screening constant σ as shown the Equation (7) and then the eigenvalue E. To obtain the minimum eigenvalue and the theoretical screening constant in which we are interested and quoted in tables 1 - 4, the variational parameter ξ and the Hylleraas parameters (J, K, M) are slightly varied that exhibit a plateau for the energy.

3. Results and Discussions

The main results of our calculations for the theoretical screening constant, the variational parameter
$\xi $ and energies concerning the 1s^{2} ^{1}S ground-state and the Doubly Excited State for 2p^{2} ^{1}D, 3d^{2} ^{1}G and 4f^{2} ^{1}I of helium isoelectronic sequence are quoted respectively in Tables 1-4. Our present results are compared with other theoretical calculations and experimental data. Then our results are converted into Rydberg for direct comparison by using the infinite Rydberg 1Ry = 0.5 a.u = 13.605698 eV. Table 1 shows a comparison of the present calculation for the 1s^{2} ^{1}S ground-state energy with the experimental data [6] [7] [8], and the theoretical results [9] [10] [11] [12] [15]. In addition the theoretical screening

Table 1. Variational parameter ξ and ground energies for 1s^{2} ^{1}S states of Helium isoelectronic sequence. The values are reported in Rydberg (Ry). 1Ry = 13.605698 eV. The screening constant sigma (σ_{theo} = 0.3101).

E^{MAOT}: Energy E of the Modified Atomic Orbital Theory (MAOT), present work. E^{sakho} (Sakho, 2006) [11] ; E^{e} (Pekeris, 1962); E^{c} (Drake, 1988) [9] ; E^{Slate}^{r} (Minkine, 1982) [15] ; E^{a}^{,b} Experimental data (Radzig, 1985 [7] ; Arnaud, 1993) [8] ; E^{ex}^{p} (Moore, 1971) [6] ; E^{d} (Utpal and Talukdar, 1999) [12] ; E^{f} (Roothaan et al., 1960) [10].

constant σ_{the}_{o} = 0.3101 is also presented and agree very good with the experimental screening constant of Slater σ_{Slater} = 0.30 determined from his rule (σ_{Slate}_{r} = 0.30 for 1s^{2} state, n = 1). For H^{−} (Z = 1), our results at −1.0098 Ry are compared with those of [10] [11] [12], respectively at −1.14069 Ry; −1.01648 Ry and −0.97459 Ry and it can see the good agreement. Comparison indicates that MAOT results agree well with each other up to Z = 10. For 10 < Z ≤ 14, our results are compared with those of [10] [11] [15]. Comparison indicates that for Z = 14, our total energy at −374.9353 Ry agree with [11] at −374.89072 Ry, with [9] at −374.88558 Ry, with [15] at −375.38022 Ry and with [10] at −374.72241 Ry. Generally, this good agreement enables to expect our results with MAOT calculation for ground-state energy up to Z = 20 to be accurate.

In Table 2, we presented the theoretical screening constant σ_{theo} = 0.3512 and reported the value of the variational parameter ξ and total energies for 2p^{2} ^{1}D

Table 2. Parameter ξ and total energies E for doubly excited 2p^{2} ^{1}D states. The values are reported in Rydberg (Ry). 1Ry = 13.605698 eV. The screening constant sigma (σ_{theo} = 0.3512).

E^{MAOT}: Energy E of the Modified Atomic Orbital Theory (MAOT), present work. SCUNC (Sakho, 2008) [22] ; E^{BIA} (Biaye et al., 2005) [2] ; E^{Ray} Experimental data (Ray et al., 1991) [21] ; E^{c} (Roy et al., 1997) [26] ; E^{g} (Ho and Bhatia, 1991) [27] ; E^{h} (Ivanov and Safronova, 1993) [4] ; E^{e} Experimental data (Hicks and Comer, 1975) [28] ; E^{f} Experimental data (Diehl et al., 1999) [29].

doubly excited states of He-like ions up to Z = 20. For this state, comparison shows that MAOT results agree well with the theoretical results of [2] [4] [21] [22] [26] [27] and also the experimental data of [28] [29]. For Z ≤ 10, comparison indicates that results agree well with each other. In addition, we calculated the energy of the ion Hydrogen H^{−} where there is no available data for comparison. For Z = 10, our results at −46.5156 Ry are compared with those of [2] [4] [22] respectively at −46.49816 Ry; −46.53271 Ry and −46.46435 Ry. Thus, this good agreement allowed us to expect our results up to Z = 20 to be also accurate.

Table 3 indicates the present MAOT calculation for the 3d^{2} ^{1}G DES that are compared with those of [2] [4] [21] [30] [31]. Here, except the ion Hydrogen H^{−} (Z = 1) where there is no available results, the agreements between the calculation are generally good up to Z = 5. Thus for 5 < Z ≤ 10, our results are compared with

Table 3. Parameter ξ and total energies E for doubly excited (3d^{2} ^{1}G) states. The values are reported in Rydberg (Ry). 1Ry = 13.605698 eV. The screening constant sigma (σ_{theo} = 0.3512).

E^{MAOT}: Energy E of the Modified Atomic Orbital Theory (MAOT), present work. E^{sakho} (Sakho,2010) [30] ; E^{a} (Bachau et al., 1991) [31] ; E^{b} (Biaye,2005) [2] ; E^{c} (Ivanov and Safronova, 1993) [4] ; E^{d} (Ho, 1989) as quoted in (Biaye, 2005) [2] ; E^{Ray} (Ray et al., 1991) [21].

those of [2] [4] [30] [31] and it can be seen that the present MAOT results agrees well with each other. Comparison indicates that (Z = 10) our results at −20.5969 Ry agrees very well with them respectively at −20.56099 Ry, −20.52000 Ry, −20.63354 Ry and −20.52658 Ry for [4]. In addition, the theoretical screening constant is also calculated (σ_{theo} = 0.3512). This good agreement allows us to expect our results with MAOT calculation for 3d^{2} ^{1}G doubly excited state up to Z = 20 to be accurate.

In Table 4, we present our calculation for the 4f^{2} ^{1}I doubly excited state and compare our results with theoretical results of [2] [30] and [32]. In our knowledge, experimental data are not available yet and also for the ion Hydrogen H^{−} (Z = 1). Here the agreements between the calculations are seen to be very satisfactory. As far as comparisons with the SCUNC results of [2] [30] and [32] are

Table 4. Parameter ξ and total energies E for doubly excited 4f^{2} ^{1}I states. The values are reported in Rydberg (Ry). 1Ry = 13.605698 eV. The screening constant sigma (σ_{theo} = 0.3508).

E^{MAOT}: Energy E of the Modified Atomic Orbital Theory (MAOT), present work. E^{Sakho} (Sakho, 2010) [30] ; E^{Biaye} (Biaye et al., 2005) [2] ; E^{Ho} (Ho, 1989) as quoted in (Biaye et al., 2005) [2].

concerned, our results for Z = 2, at −0.32906 Ry, agree well with them respectively at −0.32380 Ry; −0.35913 Ry and −0.34080 Ry. For 2< Z ≤ 10, it can be seen the present MAOT results agrees with each other, and the agreement can be seen up to Z = 10 between our results at −11.56322 Ry and that from theoretical calculations [30] and [2] respectively at −11.41048 Ry, and −11.61005. As a result, our listed data up to Z = 20 are expected to be also accurate. In addition the screening constant is also calculated σ_{theo} = 0.3508. The good agreement between the results quoted in Tables 1-4 indicate the validity of this MAOT variational procedure and his merit to calculate the (nl^{2}) ground-state and doubly excite state energies.

4. Summary and Conclusion

In this work, the variational procedure of the Modified Atomic Orbital Theory (MAOT) has been applied for the first time to the calculations of the ground-state energy 2s^{2} ^{1}S and excitation energies of the doubly 2p^{2} ^{1}D, 3d^{2} ^{1}G, 4f^{2} ^{1}I excited states of the Helium isoelectronic sequence from Hydrogen ion H^{−} to Calcium ion Ca^{18+}. It has demonstrated the possibilities to construct a new correlated wave function adapted to the correct description of the electron-electron correlations phenomena in the ground and in the doubly excited nl^{2} states of the He-like systems. These very important results obtained in this work indicate the possibility to apply the MAOT variational procedure to the treatment of atomic spectra in two electron systems and probably in more complex atomic systems. The good results give also the possibility to analyze resonance energies via a very MAOT flexible procedure, in contrast to the complex procedures of experimental and theoretical methods based on the determination of the photoionization cross-section.

Acknowledgements

The authors are grateful to the Abdus Salam International Centre for Theoretical Physics, Trieste, Italy, the Orsay Institute of Molecular Sciences (OIMS), Paris, France and the Swedish International Development Agency (SIDA) for support.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Sow, M., Sakho, I., Sow, B., Diouf, A., Gning, Y., Diop, B., Dieng, M., Diallo, A., Ba, M.D., Badiane, J.K., Biaye, M. and Wagué, A. (2020) Modified Atomic Orbital Calculations of Energy of the (2s^{2} ^{1}S) Ground-State, the (2p^{2} ^{1}D); (3d^{2} ^{1}G) and (4f^{2} ^{1}I) Doubly Excited States of Helium Isoelectronic Sequence from H^{−} to Ca^{18+}. Journal of Applied Mathematics and Physics, 8, 85-99. https://doi.org/10.4236/jamp.2020.81007

References

- 1. Sakho, I. (2012) Modified Atomic Orbital Theory Applied to the Calculation of High-Lying 2(K,T)n
^{±1,3}P^{0}Rydberg Series of He-Like Ions. Natural Science, 4, 73-78. https://doi.org/10.4236/ns.2012.41011 - 2. Biaye, M., Konté, A., Faye, N.A.B. and Wague, A. (2005) The Singlet Doubly Excited (nl)
^{2}States of Helium-Like Systems. Physica Scripta, 71, 39. https://doi.org/10.1088/0031-8949/71/1/006 - 3. Wague, A. (1990) Electron Correlation Processes in Atomic Doubly Excited States Ionization. LAMP Seminar ICTP, Trieste, July 1990.
- 4. Ivanov, A.I. and Safronova, I.U. (1993) Calculation of the Correlation Part of the Energy of Two-Electron System. Optics and Spectroscopy, 75, 506.
- 5. Ivanov, A.I. (1995) Analytical Properties of the Exact Energy of the Ground-State of a Two Electrons Atoms as Function of 1/Z. Physical Review A, 52, 3. https://doi.org/10.1103/PhysRevA.52.1942
- 6. Moore, C.E. (1971) Atomic Energy Levels (National Standard Reference Data Series, No. 35 U. S. GPO, Washington, D.C). Vol. 1, 4.
- 7. Radzig, A.A. and Smirnov, M.B. (1985) Reference Data on Atoms, Molecules and Ions. Springer, Berlin. https://doi.org/10.1007/978-3-642-82048-9
- 8. Arnaud, P. (1993) Cours de chimie physique. 3rd Edition, Dunod, Paris, Chapter 11, 88.
- 9. Drake, G.W. (1988) Theoretical Energies for the n = 1 and 2 States of the Helium Isoelectronic Sequence up to Z = 100. Canadian Journal of Physics, 66, 586-611. https://doi.org/10.1139/p88-100
- 10. Roothaan, C.C.J., Sachs, L.M. and Weiss, A.W. (1960) Analytical Self-Consistent Field Functions for the Atomic Configurations 1s
^{2}, 1s^{2}2s, and 1s^{2}2s^{2}. Reviews of Modern Physics, 32, 186. - 11. Sakho, I., et al. (2006) Calculation of the Ground-State Energy, the First Ionization Energy and the Radial Correlation Expectation Value for He-Like Atoms. Physica Scripta, 74, 180-186. https://doi.org/10.1088/0031-8949/74/2/005
- 12. Utpal, R. and Talukdar, B. (1999) Electron Correlation for Helium-Like Atoms. Physica Scripta, 59, 133. https://doi.org/10.1238/Physica.Regular.059a00133
- 13. Kleinekathöfer, U., Patil, S.H., Tang, K.T. and Toenies, J.P. (1996) Boundary-Condition-Determined Wave Function for the Ground State of Helium and Isoelectronic Ions. Physical Review A, 54, 2840. https://doi.org/10.1103/PhysRevA.54.2840
- 14. Bhattacharyya, S., Bhattacharyya, A., Talukdar, B. and Deb, N.C. (1994) Analytical Approach to the Helium-Atom Ground State Using Correlated Wavefunctions. Journal of Physics B: Atomic, Molecular and Optical Physics, 29, L147. https://doi.org/10.1088/0953-4075/29/5/003
- 15. Minkine, V., Simkine, B. and Miniaev, R. (1982) Théorie de la Structure Moléculaire. Mir, Paris, Chapter 3, 65.
- 16. Codling, K., Madden, R.P. and Ederer, D.L. (1967) Resonances in the Photo-Ionization Continuum of Ne I (20-150 eV). Physical Review, 155, 26. https://doi.org/10.1103/PhysRev.155.26
- 17. Cooper, J.W., Fano, U. and Prats, F. (1963) Classification of Two-Electron Excitation Levels of Helium. Physical Review Letters, 10, 518. https://doi.org/10.1103/PhysRevLett.10.518
- 18. Fano, U. (1983) Correlations of Two Excited Electrons. Reports on Progress in Physics, 46, 97. https://doi.org/10.1088/0034-4885/46/2/001
- 19. Rau, P.R.A. (1990) Group Theoretical Treatments of Strongly Correlated Atomic Dynamics. Reports on Progress in Physics, 53, 181. https://doi.org/10.1088/0034-4885/53/2/002
- 20. Herrick, R.D., Kellman, E.M. and Poliak, D.R. (1980) Supermultiplet Classification of Higher Intrashell Doubly Excited States of H
^{-}and He. Physical Review A, 22, 1517. https://doi.org/10.1103/PhysRevA.22.1517 - 21. Ray, D. and Mukherjee, P.K. (1991) Doubly Excited
^{1}S^{e},^{1}D^{e}and^{1}G^{e}States of He, Li^{+}, Be^{2+}and B^{3+}. Journal of Physics B: Atomic, Molecular and Optical Physics, 19, 1241-1249. https://doi.org/10.1088/0953-4075/24/6/013 - 22. Sakho, I. (2008) Screening Constant by Unit Nuclear Charge Calculations for (n s
^{2})^{1}S^{e}, (np^{2})^{1}D^{e}and (N s n p)^{1}P○ Excited States of He-Like Systems. The European Physical Journal D, 47, 37-44. https://doi.org/10.1140/epjd/e2008-00018-2 - 23. Sakho, I. (2013) Photoabsorption of H
^{-}and He via the Modified Atomic Orbital Theory: Application to the^{1}P-Rydberg State. Chinese Journal of Physics, 51, 209. https://doi.org/10.6122/CJP.51.209 - 24. Sakho, I. (2014) General Formalism of the Modified Atomic Orbital theory for the Rydberg Series of Atoms and Ions: Application to the Photoionization of Ne+. Journal of Atomic and Molecular Sciences, 5, 206. https://doi.org/10.4208/jams.110813.021414a
- 25. Schmid, E.W., Spitz, G. and Losch, W. (1987) Theoretical Physics on the Personal Computer. Springer-Verlag, Berlin, Heidelberg, New York, London, Paris, Tokyo.
- 26. Roy, K.A., Singh, R. and Deb, M.B. (1997) Density-Functional Calculations for Doubly Excited States of He, Li
^{+}, Be^{2+}and B^{3+}(^{1,3}S^{e},^{3}P○,^{1,3}D^{e},^{1,3}P○,^{1}G^{e}). Journal of Physics B, 30, 4763. https://doi.org/10.1088/0953-4075/30/21/014 - 27. Ho, Y.K. and Bathia, A.K. (1991) Complex-Coordinate Calculation of
^{1,3}D Resonances in Two-Electron Systems. Physical Review A, 44, 2895. https://doi.org/10.1103/PhysRevA.44.2895 - 28. Hicks, P.J. and Comer, J. (1975) Ejected Electron Spectroscopy of Autoionizing States Excited by Low Energy Electron Impact. Journal of Physics B, 8, 1866. https://doi.org/10.1088/0022-3700/8/11/022
- 29. Diehl, S., Cubaynes, D., Bizau, J.P., Wuilleumier, F.J., Kenedy, E.T., Mosnier, J.P. and Morgan, T.J. (1999) New High-Resolution Measurements of Doubly Excited States of Li
^{+}. Physica B, 32, 4193. https://doi.org/10.1088/0953-4075/32/17/305 - 30. Sakho, I. (2010) Calculations of (nl)
^{2}and (3lnl’) Autoionizing States in Two-Electron Systems. Physica Scripta, 82, Article ID: 035301. https://doi.org/10.1088/0031-8949/82/03/035301 - 31. Bachau, H., Martin, F., Riera, A. and Yanez, M. (1991) Resonance Parameters and Properties of Helium-Like Doubly Excited States: 2 ≤ Z ≤ 10. Atomic Data and Nuclear Data Tables, 48, 167-212. https://doi.org/10.1016/0092-640X(91)90006-P
- 32. Ho, Y.K. (1989) Resonances in Helium Atoms Associated with the N=4 and N=5 He
^{+}Thresholds. Zeitschrift für Physik D, 11, 277-281. https://doi.org/10.1007/BF01438499

List of the Symbols

σ_{theo}: (theo = theory): screening constant σ determined theoretically

E: Energy

H: Hamilton operator

N: Normalization factor

T: Kinetic Energy

W: Coulombian interaction potential operator between the two electrons

C: Coulombian interaction potential between the nucleus and the two electrons

∇_{i}: denotes the Laplacian operator ∇ with reference to the coordinates of the vector radius r_{i} which detect the position of the electron i.

∇_{1}: is the Laplacian with reference to the coordinates of the vector radius r_{1} which detect the position of the electron 1.

∇_{2}: Laplacian defines the coordinates of the vector radius r_{2} which detect the position of the electron 2.

r_{i} = denotes the position r of the electron i relative to the nucleus

r_{1} and r_{2} = denote respectively the position r of the electron 1 and the electron 2 relative to the nucleus.

$\left|{r}_{1}-{r}_{2}\right|$ : is an inter-electronic distance.

${\theta}_{i};{\phi}_{i}$ = angular variables of the electron i

S: total spin

$\hslash $ = reduced Planck constant such as $\hslash =h/2\pi $ with h the Planck constant

N = principal quantum number

l = orbitalar quantum number

Θ_{12}: angle between electron 1 and electron 2.

$\xi $ : variational parameter zeta, a Greek alphabet

H^{−}: Hydrogen ion

He: Helium atom

Li^{+}: Lithium ion

Be^{2+}: Beryllium ion

B^{3+}: Boron ion

C^{4+}: Carbon ion

N^{5+}: Nitrogen ion

O^{6+:} Oxygen ion

F^{7+}: Fluorine ion

Ne^{8+}: Neon ion

Na^{9+}: Sodium ion

Mg^{10+}: Magnesium ion

Al^{11+}: Aluminium ion

Si^{12+:} Silicon ion

P^{13+}: Phosphorus ion

S^{14+}: Sulfur ion

Cl^{15+}: Chlorine ion

Ar^{16+}: Argon ion

K^{17+}: Potassium ion

Ca^{18+}: Calcium ion