Energy Levels and Lifetimes of 1 s 2 and 1 snl ( n = 2-5 ) for Neutral Helium

In this paper, we present calculations for some of the lowest energy levels and lifetimes for neutral helium. The FAC (Flexible-Atomic-Code) is a reliable code for calculating 49 energy levels and their lifetimes. The calculation is performed up to n=5 including a series of configuration of 1s2 and 1snl. Comparison has been made with similar data published in the NIST database. A good agreement of less than1% was found for most levels expect the 1s2p P2 level. This proves the reliability of our results. New values for lifetimes are presented for the first time. article History Received: 15 December 2017 Accepted: 28 December 2017

introduction Atomic data are very much needed for the modeling of plasmas 1 .They are useful for many applications in astrophysics and nuclear fusion Tokamak.Moreover, the ITER project needs accurate atomic data for a wide area of ions.The simplest multi-electron system, He-like can play an important role in providing the needed accurate atomic data.In previous work, we provided calculations of He-like neon 2 .
We just focus on this study in neutral Helium.For this Nobel gas, the atom contains two protons and two electrons.The energy structure of the 1snl configuration is mainly dominated by the electron-nucleus and electron-electron Coulomb interactions.The separations of levels having belonging to the same n multiplicity and having l = s, p, or d are mainly determined by Spin-orbit interactions between the electrons.The final state is given by the 2S+1 L notation where n,l,S,L are the standard quantum numbers.The quantity 2S + 1 is the multiplicity of the term.The S and L momentums are coupled to obtain, J = S + L, for a given level.The level is denoted as 2S+1 L J .J being the total angular momentum.
Many experimental and theoretical studies were performed for this element.Energy levels have been published by the National Institute of Standards and Technology (NIST) and are available at their website 3 .Sow et al. 4 performed calculations for the 2s 2 1 S, 2p 2 1 D, 3s 2 1 S, 3p 2 1 D, 3d 2 1 G, 4p 2 1 D, 4d 2 1 D, 4f 2 1 I doubly excited states of 2≤ Z ≤15.The energy eigen-values of doubly excited 2pnp( 1 P e ) (n=3-8) and 2pnp ( 3 P e ) (n=2-7) bound states of a neutral helium atom were calculated under the weakly coupled plasma screening by K. Saha et al. 5 .Doubly excited states in helium are calculated by E. Lindroth 6 with a finite discrete spectrum for states with electrons in the n=2 and n=3 states.An approach to calculating the energies and widths of resonances for neutral Helium was developed on the basis of the stabilization method, the energies of 28 resonances of n S symmetry with the spin multiplicities n=1, 2, 3, and 4 were calculated by I. A. Misurkin et al. 7 Therefore, in this paper, we just interest on calculations of the singly excited energy levels for neutral Helium, namely He I. We employed the fully relativistic code Flexible Atomic Code (FAC) of Gu 8 .FAC code provides many different atomic parameters such as energy levels, transitions rates and lifetimes.In the rest of paper, we shall give, for He I, energy levels and lifetimes of the lowest 49 levels belonging to 1s 2 and 1snl (with n ≤ 5; 0 ≤ l ≤ n-1).

Theoretical Method
We employed for our calculations the widely used FAC code of Gu 8 .We simply give here a short description of the theoretical method used by FAC.
By diagonalizing the relativistic Hamiltonian H, we get the energy levels of an N electrons atom 9 where H D (i) is the single-electron Dirac Hamiltonian.The basis states ν φ , which are usually referred to as configuration state functions (CSF), are antisymmetric sums of products of N one-electron Dirac spinors m nκ ϕ .... (2)   where X km is the usual spin-angular function.n is the principal quantum number, k is the relativistic angular quantum number and m is the z-component of the total angular momentum j.
The approximate atomic state functions are obtained by mixing the basis states using the same symmetries where b v are the mixing coefficients obtained from diagonalizing the total Hamiltonian.

Choice of Local Central Potential
To build the Hamiltonian matrix, the one-electron radial orbital must be known.According to the standard Dirac-Fock-Slater method, the large and small components, P nk and Q nk , satisfy the coupled Dirac equation for a local central field V (r) where a is the fine structure constant and ε nk are the energy eigenvalues of the radial orbitals.

Solution of Dirac Equations
The radial orbitals sought have a direct influence on the potential, so the Eq. ( 4) requires a self-consistent iteration.In each iteration, the orbitals from the previous step are used to derive the potential.Consequently, solving the eigenvalue problem using known potential is sufficient.As is standard, we convert Eq. ( 4) into a Schrodinger-like equation in two steps: eliminating the small component and performing the transformation 10 Under this transformation, we have .... (6)   where U(r) is an effective potential defined as .... (7)   We use the standard Numerov method to solve Eq. ( 6).However, it is customary to perform another transformation before seeking the solution .... (8)   where t(r) as a function of radial distance is suitably chosen so that a uniform grid can be used in the new variable t, and the corresponding transformation on the wave function is to bring the differential equation for G a (t) to a Schrodinger-like form, i.e., without the first derivative term .... (9)   However, for free orbitals with sufficiently high energy, solving Eq. ( 9) in a conventional way becomes impractical.We shall use a different approach for continuum states, namely, the phase amplitude method.
The minimum and maximum radial distances, r min and r max , in setting up the radial grid are chosen as .... (10)   where z eff is the residual charge of the atomic ion that the electrons experience at large r.The low-n and high-n states are treated differently.The dividing n 0 is determined by the choice of r max , specifically, For n ≤ n 0 , the orbitals were found by outward and inward integration of Eq. ( 9) with zero amplitudes at both ends, and matching at the outer classical turning point.Node counting is used to pick out the appropriate solution corresponding to the quantum numbers n and l.The wave functions are then normalized by numerical integration.For n > n 0 , Eq. ( 9) is integrated outward until r = r core , where the potential has reached its asymptotic Coulomb value.For r > r core , the wave function is the exponentially decaying Whit taker function where , and λ = l in the nonrelativistic limit 8 .In the relativistic case, the asymptotic behavior of the effective potential is modified according to Eq. ( 7), and corresponds to ..

..(12)
To normalize the wave function, we note that the correct normalization is given by 11  where µ is the quantum defect.The quantum defect of a Rydberg atom refers to a correction applied to the equations governing Rydberg atom behavior.For a non-hydrogen atom -alkali for example-the binding energy of the Rydberg states is ε = R y /( nδ) 2 where δ is the quantum defect and R y is the Rydberg constant.For high n states we are concerned with, (ν n ) = 1 is a very good approximation.

Results and Discussions
The calculation is performed up to n=5 which generate up to 49 levels.The 1s 2 , 1snl (n=2-5) configurations are given in Table 1.We compare in this table our energy levels from FAC code with the level energies published by NIST 3 .The differences presented by percentage in the Table are very small.Our calculated energy levels agree within 1%, the only case where the difference is larger being for the 1s2p 3 P 2 level.One can see that results from these two calculations match well for most of the levels and proves that our results are consistent.We can state with confidence that the results are in good agreement to the other published values for the energy levels of He I. New values of lifetimes are presented in the same Table 1.

Conclusion
We have presented in this paper results for energy levels among the lowest 49 levels, for He I. Based on the experimental published results in NIST database, our energy levels are accurate to better than 1%.Moreover, we presented the results for the lifetimes for almost levels of study.Good agreement between our calculated energy levels for He I and the available NIST data reflects the quality of calculation of the wave functions.
As we don't have other results of lifetimes, we expect that the present set of results will be highly useful for comparison with other future experimental work.

acknowledgments
This work has been realised with the financial support of the Tunisian Ministry of Higher Education and Scientific Research.