JMPJournal of Modern Physics2153-1196Scientific Research Publishing10.4236/jmp.2023.1412094JMP-129265ArticlesPhysics&Mathematics Variational Calculations for (ns<sup>2</sup>) <sup>1</sup>Se, (np<sup>2</sup>) <sup>1</sup>D<sup>e</sup> and (nd<sup>2</sup>) <sup>1</sup>G<sup>e</sup> Resonance States for He Isoelectronic Sequences below the n = 2, 3 and 4 Hydrogenic Thresholds MatabaraDieng1*Department of Physics, UFR-SATIC, University Alioune Diop of Bambey, Bambey, Senegal0911202314121617163312, September 202321, November 2023 24, November 2023© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

This work presents results of the different parameters which characterize the nonrelativistic Hamilton operator for the helium atoms allowing us to solve the Schr &ouml;dinger equation. The total energy is decomposed into three terms allowing to separate the kinetic energy, the electrons-nucleus interaction energy and the electron-electron interaction energy of the (2s 2, 3s 2 and 4s 2) 1S e, (2p 2, 3p 2 and 4p 2) 1D e and (3d 2 and 4d 2) 1G e resonance singlet states of the helium isoelectronic sequences. The states have been defined by using special forms of the Hylleraas type wave functions. The calculations have been carried out in the framework of the variational method using configuration interaction basis states with a real Hamiltonian. The agreement of the energy value of other states between the present theoretical values available in the literature is excellent. But as for the comparison of the kinetic energies, the electrons-nucleus energies interaction and the electron-electron interaction energies, we note a slight difference with the theoretical values common in literature.

Hylleraas Method Helium-Like Ions Systems Singlet Doubly-Excited States Matrix Elements Autoionization1. Introduction

Double photoionization occupies an important place because of the crucial role that electron-electron interaction plays in the ejection of two electrons by the photoabsorption of a single photon. Electron-electron correlation is also an essential feature in the process of simultaneous excitation and ionization when an electron is ejected and the residual ion is in excited state. Many years ago, in helium atoms, all of the spectra transitions emitted by the atoms were visible and only one electron was in excited state. No spectra transition corresponding to the de-excitation of atoms with two electrons having been excited was observed. At the time, we did not have any light sources possessing photons with sufficient energies for double excitation of the atoms. It was in 1963, at the National Bureau of Standards (NBS) in Gaithersburg, that Madden and Codling [1] used the first light source with a continuous UV spectrum that was able to describe doubly excited helium atoms. The experimental observations of Madden and Codling allowing the identification of the painfully excited states of helium [1] [2] and the qualitative study of the correlation between two doubly excited electrons provided by the work of Cooper et al. [3] have been sources of deep inspiration for many researchers working in the study of doubly excited states (DES) of helium and helium-likes ions [4] - [9] .

Over the years, several theoretical calculations have been developed for the descriptions of the doubly excited states, especially for the doubly excited states of helium and helium-like ions for which many experimental data have confirmed. A recent photoionization experiments have been performed by Morgan and Ederer [10] to observe resonance with excellent precision and Bizau et al. [11] have studied also the photoionization descriptions with excitation above the n = 2 threshold. On the theoretical part, different approaches have been used for that purpose by the theorists. Among these theoretical methods we can note: Berrington et al. [12] have used a six-state pseudostate close-coupling method, the photoionization process has been studied by a Feshbach projection method by Bhatia and Temkin [13] , Salomonson et al. [14] have employed the many-body perturbation theory to study photoionization with excitation, Ojha [15] has applied multichannel quantum-defect theory to study these processes near the n=2 threshold, the variational method [16] [17] , the screening constant by unit nuclear charge method [18] etc.

In general, combinations of functions are necessary in order to give a good description of doubly excited states. The wave functions are used in the incomplete hydrogenic product combined with the truncated diagonalization method, the Slater wave functions in the complex rotation method, the Hylleraas-type wave functions in the complex rotation method, the variational method and the diagonalization method. In the case of the complex rotation method, many authors use the Hylleraas-type wave functions for the intrashell doubly excited energies and width calculations. For the intershell doubly excited energies and width calculations of two electron systems, they use the Slater wave functions [5] .

The present work is an extension of the earlier calculations of two-photon excitation and ionization energies of the Rydberg helium [7] [8] . Special forms of the Hylleraas-type wave functions constructed in the previous paper have been used for the intrashell doubly excited energies calculations of the helium-like ions. In this communication we reported, without the Slater wave functions, the total energy and the various parameters as the kinetic energy, the electrons-nucleus interaction energy and the electron-electron interaction energy of the (2s2, 3s2 and 4s2) 1Se, (2p2, 3p2 and 4p2) 1De and (3d2 and 4d2) 1Ge resonance states of the helium isoelectronic sequences with Z = 2 - 10, by using the same special forms of Hylleraas type wave functions. The calculations have been carried out in the framework of the variational method using configuration interaction basis states with a real Hamiltonian. In the calculations we have neglected the Feshbach shifts. The aim of this work is to show the possibility of these special forms of Hylleraas-type wave functions to describe some intrashell singlet of doubly excited states of two electron atoms. This work allowed the use of small basis sets (13 terms) in the calculations.

The rest of this paper is arranged as follows: Section 2 presents the theoretical procedures used in this work. Section 3 is about the comparisons of our results with the other theoretical calculations data, through curves. Here also, we have no experimental data and there is not much theoretical data available, for comparisons. Finally, conclusion is drawn.

2. Method of Numerical Calculations2.1. Hamiltonian and Wave Functions

The Schrödinger equation for the relative motion of the helium-like ion, which interacts with each other by a spherically symmetric potential, can be written as:

H Φ ( r 1 , r 2 ) = E Φ ( r 1 , r 2 ) (1)

where H is the Hamilton operator for the helium atom we shall use the approximation is given by:

H = − ℏ 2 m ( Δ 1 + Δ 2 ) − Z e 2 r 1 − Z e 2 r 2 + e 2 | r 1 − r 2 | , (2)

The vector r 1 and r 2 denote the positions of the two electrons, m the mass of an electron, e the elementary charge and Z the nuclear charge number, Δ 1 and Δ 2 the Laplace operator of the two electrons.

The Hamilton operator can be consists of three parts:

H = T + C + W , (3)

where T, C and W are respectively the kinetic energy operator of the two electrons, the Coulomb interaction operator between the atomic nucleus and the two electrons and the Coulomb interaction operator between the two electrons:

T = − ℏ 2 m ( Δ 1 + Δ 2 ) , (4)

C = − Z e 2 r 1 − Z e 2 r 2 , (5)

W = e 2 | r 1 − r 2 | , (6)

In this case of the Hamilton operator, all magnetic and relativistic effects together with the motion of the atomic nucleus are neglected.

In this article, Φ ( r 1 , r 2 ) are the trials non-orthogonal of two-electron wave functions that we have considered for the description of the intrashell singlet doubly excited states of the helium-like ions. There are special constructions of the incomplete hydrogenic wave functions and Hylleraas type wave functions as follows:

Φ j , k , m ( r 1 , r 2 ) = χ j , k , m ( r 1 , r 2 ) φ j , k , m ( r 1 , r 2 ) , (7)

With

φ j , k , m ( r 1 , r 2 ) = ( r 1 + r 2 ) j ( r 1 − r 2 ) k | r 1 − r 2 | m exp [ − λ ( r 1 + r 2 ) ] , (8)

which are Hylleraas ground state wave functions of the helium-like ions [19] , including electron correlation effects; j, k, m are Hylleraas parameters with (j, k, m ≥ 0), λ is a coefficient defined by

λ = Z α n a 0 , (9)

where Z, α, a0 and n are respectively the nucleus charge number, variational parameters, Bohr’s radius and the principal quantum number.

These wave functions φ j , k , m ( r 1 , r 2 ) are not orthogonal.

The set of parameters (j, k, m) define the basis states (i.e. the configurations).

The even values of k define the symmetric wave functions describing the singlet states.

The wave functions Φ j , k , m ( r 1 , r 2 ) , are incomplete hydrogenic wave functions and can be expressed as follows:

Φ n , l ( r 1 , r 2 ) = ( 2 r 1 2 r 2 ) l { ∑ υ = 0 υ = n − l − 1 [ ( n 2 a 0 2 λ 2 ) ( 2 r 1 2 r 2 ) ] υ } ( r 1 + r 2 ) j ( r 1 − r 2 ) k       × | r 1 − r 2 | m exp [ − λ ( r 1 + r 2 ) ] , (10)

With

n = ( j , k , m ) : principal quantum number.

l : orbital angular moment.

The interesting feature in the wave functions Φ j , k , m ( r 1 , r 2 ) , is that they contain an electron correlation term: | r 1 − r 2 | , which represents the angular part of the wave functions instead of the spherical harmonic in the other Hylleraas type wave functions.

This electron correlation term plays an important role in our trial wave functions for the description of the intrashell singlet doubly excited states. The wave functions Φ j , k , m ( r 1 , r 2 ) have also the advantage that, in the eigenvalue calculations E, the exhibition of a plateau and the convergence of the minima of the functions ( d E / d α = 0 ) arise quickly for small basis set (13 terms).

2.2. Calculation Procedures

The final form of the wave functions of the intrashell singlet doubly excited state including the correlation effects due to the mixing of configurations can be expressed as follows:

ψ n , l ( r 1 , r 2 ) = ∑ n a n Φ n , l , (11)

where a n are the eigenvectors which can be determined by solving the Schrödinger equation.

H ψ n , l ( r 1 , r 2 ) = E ψ n , l ( r 1 , r 2 ) , (12)

The representation of the Schrodinger equation on the non-orthogonal basis leads to the general eigenvalue equation:

∑ n ' ( H n , n ' − E N n , n ' ) a n ' = 0 , (13)

with:

N n , n ' = N J , K , M = 〈 ψ n , l | ψ n ' , l 〉 , (14)

H n , n ' = H J , K , M = 〈 ψ n , l | H | ψ n ' , l 〉 , (15)

H n , n ' = H J , K , M = 〈 ψ n , l | T + C + W | ψ n ' , l 〉 , (16)

H n , n ' = 〈 ψ n , l | T | ψ n ' , l 〉 + 〈 ψ n , l | C | ψ n ' , l 〉 + 〈 ψ n , l | W | ψ n ' , l 〉 , (17)

T n , n ' = T J , K , M = 〈 ψ n , l | T | ψ n ' , l 〉 = 〈 T 〉 , (18)

C n , n ' = C J , K , M = 〈 ψ n , l | C | ψ n ' , l 〉 = 〈 C 〉 , (19)

W n , n ' = W J , K , M = 〈 ψ n , l | W | ψ n ' , l 〉 = 〈 W 〉 , (20)

H J , K , M = T J , K , M + C J , K , M + W J , K , M = 〈 E 〉 = 〈 T 〉 + 〈 C 〉 + 〈 W 〉 , (21)

Wherein N J , K , M are the matrix elements of normalisation factor, H J , K , M the matrix elements of Hamilton operator, T J , K , M the matrix elements of kinetic energy operator of the two electrons, C J , K , M the matrix elements of electrons-nucleus interaction energy operator and W J , K , M the matrix elements of electron-electron interaction energy operator.

For example, we present the result of the different parameters of 2s2 state:

- Matrix elements of normalisation factor:

N J , K , M = 2 π 2 M + 2 ( 1 K + 1 − 1 K + 3 − 1 K + M + 3 + 1 K + M + 5 ) ( J + K + M + 5 ) ! ( 1 2 λ ) J + K + M + 6                             + 4 π 2 M + 2 ( n a 0 λ ) 2 ( 1 K + 1 − 2 K + 3 + 1 K + 5 − 1 K + M + 3 + 2 K + M + 5                             − 1 K + M + 7 ) ( J + K + M + 7 ) ! ( 1 2 λ ) J + K + M + 8                             + 2 π 2 M + 2 ( n a 0 λ ) 4 ( 1 K + 1 − 3 K + 3 + 3 K + 5 − 1 K + 7 − 1 K + M + 3 + 3 K + M + 5                             − 3 K + M + 7 + 1 K + M + 9 ) ( J + K + M + 9 ) ! ( 1 2 λ ) J + K + M + 10       for     K = 0 , 2 , 4 , 6 , ⋯ N J , K , M = 0     for     K = 1 , 3 , 5 , 7 , ⋯ (22)

- Matrix elements of electrons-nucleus interaction energy:

C J , K , M = 2 π 2 M + 2 ( 1 K + 1 − 1 K + M + 3 ) ( J + K + M + 4 ) ! ( 1 2 λ ) J + K + M + 5                           + 16 π 2 M + 2 ( n a 0 λ ) 2 ( 1 K + 1 − 1 K + 3 − 1 K + M + 3 + 1 K + M + 5 ) ( J + K + M + 6 ) ! ( 1 2 λ ) J + K + M + 7                           + 8 π 2 M + 2 ( n a 0 λ ) 4 ( 1 K + 1 − 2 K + 3 + 1 K + 5 − 1 K + M + 3 + 2 K + M + 5                           − 1 K + M + 7 ) ( J + K + M + 8 ) ! ( 1 2 λ ) J + K + M + 9       for     K = 0 , 2 , 4 , 6 , ⋯ C J , K , M = 0       for     K = 1 , 3 , 5 , 7 , ⋯ (23)

- Matrix elements of electron-electron interaction energy:

W J , K , M = e 2 N J , K , M − 1 (24)

- Matrix elements of kinetic energy:

T J , K , M = 2 ( λ 2 N J , K , M − J λ N J − 1 , K , M + j j ' N J − 2 , K , M + k k ' N J , K − 2 , M + m m ' N J , K , M − 2 )                         + 1 2 [ − M λ ( C J , K , M − C J , K + 2 , M − 2 ) + ( m j '   + j m ' ) ( C J − 1 , K , M − C J − 1 , K + 2 , M − 2 )                         + ( m k '   + k m ' ) ( C J + 1 , K , M − 2 − C J − 1 , K , M ) ] (25)

All the other states are calculated with this same way.

The intrashell singlet doubly excited wave functions were found in the basis containing the configurations with the following condition for the Hylleraas parameters j + k + m ≤ 3, corresponding to the basis dimension D = 13.

In order to obtain the minimum eigenvalue in which we are interested in the calculations are carried out for various values of the parameter α.

The eigenvalues E obtained in the present calculations follow the Hylleraas-Undheim theorem [20] and do not include the Feshbach shifts because of the use of the incomplete basis sets of the wave functions.

According to the Hylleraas-Undheim theorem [20] , a good approximation for the eigenvalues is obtained when the minima of the functions ( d E / d α = 0 ) converge with increasing values of the dimension D and when the functions exhibit a plateau.

3. Results and Discussions

In Figures 1-8, we show the variation, of the kinetic energies 〈 T 〉 , the electrons-nucleus interaction energies 〈 C 〉 , the electron-electron interaction energies 〈 W 〉 and the total energies fof the present work for (2s2) 1Se, (3s2) 1Se, (4s2) 1Se, (2p2) 1De, (3p2) 1De, (4p2) 1De, (3d2) 1Ge and (4d2) 1Ge states of helium-like ions with Z = 2 - 10.

We compare in the Figure 9 and Figure 10, our results of the kinetic energies 〈 T 〉 for (2s2) 1Se and (3s2) 1Se states of helium-like ions (Z = 2 - 10) with the

theoretical values of Gning et al. [21] , Konté et al. [22] which have use respectively the complex rotation method implemented in Scilab and the method of stabilization using correlated Hylleraas-type wave functions. We note that the kinetic energy increases when Z is large and it varies like a parabolic function. We note also a little disagreement between our results and the other theoretical calculations of Konté et al. [22] when Z increases and a big difference with the results of Gning et al. [21] . The increase of the nucleus charge number Z leads more kinetic energy.

Figure 11 and Figure 12 compare our results of the electrons-nucleus energies 〈 C 〉 for (2s2) 1Se and (3s2) 1Se states of helium-like ions (Z = 2 - 10) with the theoretical values of Gning et al. [21] . and Konté et al. [22] . We find that the variational of 〈 C 〉 looks like a hyperbolic function and it is substantially opposite to twice 〈 T 〉 . It is seen that; a little discrepancy appears when the charge number increases between our results and the theoretical calculations of Konté et al. [22] and a big gap with Gning et al. [21] results. The absolute value of electron-nucleus interaction energy increase as the nucleus charge number Z.

Figure 13 and Figure 14 contain comparison electron-electron interaction energies 〈 W 〉 of the present work with results of Gning et al. [21] and Konté et al. [22] for (2s2) 1Se and (3s2) 1Se states of helium-like ions (Z = 2 - 10). We see 〈 W 〉 is linear function of Z. Comparison done with results of Gning et al. [21] and Konté et al. [22] shows a discrepancy when Z increases. This is the only

comparison where we have a gread disagreement with the results of Konté et al. [22] . Electron-electron correlation is essential for the process of excitation and ionization and is gread when the nucleus charge number Z increase.

A comparison is made in Figure 15 and Figure 16 with the total energies of

Gning et al. [21] and Konté et al. [22] for (2s2) 1Se and (3s2) 1Se states of helium-like ions up to Z = 10. We note generally a good agreement between our results with those of Gning et al. [21] and Konté et al. [22] .

4. Conclusion

We have estimated separately, using special forms of Hylleraas-type wave functions, the kinetic energies, the electrons-nucleus interaction energies, the electron-electron interaction energies and the total energies for (ns2) 1Se, (np2) 1De and (nd2) 1Ge resonance states for He isoelectronic sequences below the n = 2, 3 and 4 hydrogenic thresholds up to Z = 10. The calculations have been done in the framework of the variation method using configuration interaction basis states with a real Hamiltonian. Ours results for total energies are in good agreement with cited theoretical literatures values with other methods. For the kinetic energies and the electrons-nucleus interaction energies, we note a slight disagreement between our results and those of the other calculations. For the electron-electron interaction energies, this is the only comparison where we have a great disagreement with the other results. In a general way, we have presented in this paper satisfactory results which will be beneficial for theoretical and experimental research. Examples for other states as n l n l ' ( l ≠ l ' ) will be considered in the future.

Data Availability Statement

This manuscript has no associated data or the data will not be deposited. [Authors’ comment: Data are available upon request.]

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Dieng, M. (2023) Variational Calculations for (ns2) 1Se, (np2) 1De and (nd2) 1Ge Resonance States for He Isoelectronic Sequences below the n = 2, 3 and 4 Hydrogenic Thresholds. Journal of Modern Physics, 14, 1617-1633. https://doi.org/10.4236/jmp.2023.1412094

ReferencesMadden, R.P. and Codling, K. (1963) Physics Letters, 10, 516-518. https://doi.org/10.1103/PhysRevLett.10.516Madden, R.P. and Codling, K. (1965) The Astrophysical Journal, 141, 364-375. https://doi.org/10.1086/148132Cooper, J.W., Fano, U. and Prats, F. (1963) Physical Review Letters, 10, 518-521. https://doi.org/10.1103/PhysRevLett.10.518Herrick, D.R. and Sinanoglu, O. (1975) Physical Review A, 11, 97. https://doi.org/10.1103/PhysRevA.11.97Ho, Y.K. (1996) Zeitschrift für Physik D, 38, 191-196. https://doi.org/10.2307/1270431Bachau, H., Martin, F., Riera, A. and Yanez, M. (1991) Atomic Data and Nuclear Data Tables, 48, 167-212. https://doi.org/10.1016/0092-640X(91)90006-PBiaye, M., Dieng, M., Sakho, I., Ndao, A.S. and Wagué, A. (2009) Chinese Journal of Physics, 47, 166.Dieng, M., Sakho, I., Biaye, M. and Wagué, A. (2010) Chinese Journal of Physics, 48, 38.Dieng, M., Sakho, I., Biaye, M. and Wagué, A. (2014) Radiation Physics and Chemistry, 96, 19-26. https://doi.org/10.1016/j.radphyschem.2013.08.009Morgan, H.D. and Ederer, D.L. (1984) Physical Review A, 29, 1901.Bizau, J.H., Wuilleumier, F., Ederer, D.L., Chang, T.N., Krummacher, S. and Schmidt, V. (1982) Physical Review Letters, 48, 588. https://doi.org/10.1103/PhysRevLett.48.588Berrington, K.A., Burke, P.G., Fon, W.C. and Taylor, K.T. (1982) Journal of Physics B: Atomic, Molecular and Optical Physics, 15, L603. https://doi.org/10.1088/0022-3700/15/17/006Bathia, A.K. and Temkin, A. (1984) Physical Review A, 29, 1895.Salomonson, S., Carter, S.L. and Kelly, H.P. (1985) Journal of Physics B, 18, L149. https://doi.org/10.1088/0022-3700/18/6/001Ojha, P.C. (1984) Journal of Physics B, 17, 1807. https://doi.org/10.1088/0022-3700/17/9/019Ray, D. and Mukherjee, P.K. (1991) Journal of Physics B, 24, 1241. https://doi.org/10.1088/0953-4075/24/6/013Ray, D., Kundu, B., Mukherjee, P.K., Ohtsuki, K. and Ohno, K. (1989) Physics Letters A, 136, 423-427. https://doi.org/10.1016/0375-9601(89)90427-1Sakho, I., Konté, A., Ndao, A.S., Biaye, M. and Wagué, A. (2010) Physica Scripta, 82, Article ID: 035301. https://doi.org/10.1088/0031-8949/82/03/035301Schimid, E.W., Spitz, G. and Losch, W. (1987) Theoretical Physics on Personal Computer. Springer-Verlag, Berlin.Hylleraas, E. and Undheim, B. (1930) Zeitschrift für Physik, 65, 759-772. https://doi.org/10.1007/BF01397263Gning, Y., Sow, M., Traoré, A., Dieng, M., Diakhaté, B., Biaye, M. and Wagué, A. (2015) Radiation Physics and Chemistry, 106, 1-6.Konté, A., Ndao, A.S., Biaye, M. and Wagué, A. (2006) Physica Scripta, 74, 605. https://doi.org/10.1088/0031-8949/74/6/002