key: cord-0977179-8noddf8c
authors: Siddiq, R.; Hameed, M. N.; Zaheer, M. H.; Khan, M. B.; Uddin, Z.
title: Rydberg energies and transition probabilities of Li I for np–ms (m ≤ 5) transitions
date: 2022-03-26
journal: Beni Suef Univ J Basic Appl Sci
DOI: 10.1186/s43088-022-00224-0
sha: d3af69fda5519b28b6a45cce01ec141cde091c5e
doc_id: 977179
cord_uid: 8noddf8c

BACKGROUND: Mathematical modeling provides grounds for understanding scientific systems theoretically. It serves as a guide for experimentalists in determining directions of investigation. Recently, the Covid-19 pandemic has caused disturbances in almost every walk of life. Scientists have played their role and have continued research on the effects of the pandemic. Various mathematical models have been used in different branches of science (Djilali et al. in Phys Scr 96 12 124016, 2021; Math Biosci Eng 18(6):8245–8256, 2021; Zeb et al. in Alex Eng J 61(7):5649–5665). Well-established mathematical models give results close to those obtained by experiments. The Weakest Bound Electron Potential Model is one such model, which explains hydrogen-like atoms and ions. This model has been used extensively for hydrogen-like atoms and ions to calculate energies of Rydberg levels and ionization energies. This model has been used extensively for hydrogen-like atoms and ions to calculate energies of Rydberg levels and ionization energies. RESULTS: This paper presents the energies of the Rydberg series, 2s2ns, and 2s2np of Li I, calculated using WBEPM. The energies are used to calculate transition probabilities from np to 2s, 3s, 4s, and 5s levels. The transition probabilities are compared with corresponding values in published data where available. The agreement with known values is good; most of the transition probabilities calculated in this work are new. A computer program was developed to find the value of the dipole matrix element. The calculations were further verified by calculating the lifetimes of some low-lying levels. CONCLUSIONS: Four series of Li I have been studied, and energies of the Rydberg levels in the series were calculated. The energies then are used to calculate transition probabilities from np to ms transitions, where m = 2, 3, 4, & 5 and n = 1–15. The results are compared where available. An excellent agreement with previously published data shows the reliability of calculations. Most of the transition probabilities are new.

Atomic lithium is less abundant in the universe because of its nuclear instability; it is a part of compounds. Lithium is an ideal atom for the spectroscopic study of astronomical data.; its 6708 Å doublet is present in K-type and cooler stars. The 6708 Å line is also among those with the maximum range of energy distribution. Furthermore, it is present in a part of the spectrum that has little effect from stellar atmospheric absorption or the crowding of lines [4] . More than 300 papers have been published on lithium describing various spectroscopic characteristics; however, only a few papers are on transition probabilities [5] . Heavens used central field approximation to calculate lifetime and transition probabilities between lower excited states of alkali atoms [6] . Fox and Dalgarno used a theoretical method and obtained eigenvalues and eigenfunctions for 1s 2 2p 2 P o , 1s2p 2 2 P, and 1s2p 2 2 D states of lithium sequence. With the help of eigenfunction, they also calculated transition probabilities for these states [7] . Lingard and Nielson performed extensive calculations based on the numerical approximation method to

Beni-Suef University Journal of Basic and Applied Sciences Siddiq et al. Beni-Suef Univ J Basic Appl Sci (2022) 11:42 find dipole transition probabilities, oscillator strengths, lifetimes, and branching ratios of alkali isoelectronic sequences [8] . Kostelecky and Nieto presented a potential function that reproduces spectra in the limits of quantum defect theory. With the help of this potential function, they obtained analytical wavefunctions and calculated transition probabilities of some transitions in lithium and Sodium atoms [9] . Fischer et al. presented multiconfiguration Breit-Pauli energy levels, lifetimes, and transition probabilities for the lithium sequence up to Z = 8 [10] . Zheng et al. employed Weakest bound electron potential model WBEPM theory to calculate transition probabilities for lithium atoms and lithium-like ions [11] . Celik used the weakest bound electron potential model and quantum defect theory to find atomic transition probabilities for the transitions in atomic lithium [12] . 2009 Wiese and Fuhr critically evaluated available literature sources and tabulated atomic transition probabilities for allowed and forbidden transitions in hydrogen, helium, and lithium [13] . Experimentally very few scientists have measured transition probabilities for the lithium atom. Most of the experimental work has measured transition probabilities of 2p and 3p levels of Li I. Anderson et al. measured lifetime and oscillator strength of 2p state of Li I through-beam foil excitation method [14] . Gaupp et al. also found the lifetime and oscillator strength of the 2p state through laser excitation in an atomic beam experiment [15] . In this work, we calculated transition probabilities of 110 transitions using the semi-empirical non-relativistic method.

In atoms, the outermost electrons are weakly bound to the nucleus as compared to the atomic core electrons. Lithium and other atoms have a single electron in their outermost shell; this electron is the weakest bound electron. Zheng et al. developed a theory, known as Weakest Bound Electron Potential Model (WBEPM), for the spectroscopic characteristics of such atoms [16] . The weakest bound electron model theory allows us to deal with many-electron systems as a binary system; the weakest bound electron is considered an electron in a field created by a nucleus dressed up with an atomic core (consists of other electrons surrounding the nucleus). The weakest electron is excited to a higher energy level moves in an orbit that has a larger period; this reduces the coupling between core electrons and the weakest bound electron. Unlike Self consistent field (SCF) method in which each electron moves in average potential arising from other electrons, the WBEPM is based on the consideration of successive dynamic ionization [16] . The potential in which the weakest bound (WB) electron moves is given by The radial part of the Schrodinger equation with this potential is as follows

The last two terms have the same denominator to be combined. Here A = −Z * where Z * is the effective charge of the nucleus, herel * = l − δ n , n * = n − δ n , are effective orbital quantum number and principal quantum number, respectively, δ n is the quantum defect in the quantum numbers and is given as a function of n. The transition probability (A fi ) of a transition for spontaneous emission between levels n f , l f &(n i , l i ) , is given by

In above expression E f > E i and are the energies of upper and lower states, respectively, S represents electric dipole line strength. In lighter atoms, LS coupling dominates; therefore, line strength can be found by

The terms in the bracket contain two 6 J symbols and the matrix element P (1) l i l f , which is given by The integral in the above expression can be found by the formula based on WBEPM theory [11, 12] .

Very few experimental results of transition probabilities and lifetimes are available for the transition studied in this work. In theory, only one study is found in which few transitions probabilities of lithium transitions are presented [8, 17] . Lindgren [8] used numerical coulomb approximation to calculate transition probabilities of various series for which n ≤ 12, l ≤ 4. The radiative transition probabilities of four Rydberg series in lithium are calculated in this work. This work employs a semiempirical non-relativistic weakest bound electron potential model and consists of two parts. In the first part, Rydberg energies are calculated by finding corresponding quantum defects in principle quantum numbers of the Rydberg levels. The experimental values are used to calculate quantum defects and energies of the Rydberg levels and modified quantum numbers n* (Principal) and l* (orbital). The second part uses the energies and quantum defects to calculate transition probabilities for various transitions from np to 2s, 3s, 4s, and 5s levels. The quantum defects for Rydberg series 2s 2 ns and 2s 2 np are found using the following equations,

These quantum defects are then used to calculate the Rydberg energies of both series. WBEPM uses Eq. (10) for the calculations of matrix element, we developed a python program that calculates dipole matrix element

using wavefunction given in Eq. (9) . The python program's reliability is checked by calculating the radii of various ns and np levels, compared with the corresponding values in the literature. An excellent agreement was found between them. A generalized formula for each ns and np levels radii is generated and given in Eqs. (13) (14) (which is a simple quadratic equation in 'n' (principal quantum number)). The energies and the values of matrix elements are used to calculate transition probabilities for the transitions from series np to state 2s, 3s, 4s, and 5s up to n = 30. The results of transition probabilities are given in Tables 1, 2 , 3, and 4. The transition probability results are compared with the corresponding values listed in NIST line data [5] and previously published results [6, [8] [9] [10] [11] [12] . The experimental and theoretical data of transition probabilities is available only for a few transitions, especially for transitions from lower levels; therefore, most of the results are new. The calculation of transition probabilities using WBEPM requires modified values of principal and orbital quantum numbers; the modified values of n * and l * are determined by the quantum defects calculated in the first part. The transition probabilities are calculated using Eq. (7). The radii of ns and np series can be calculated by the simple quadratic equations given below:

The calculated transition probabilities for transitions np to 2s, 3s, 4s, 5s are shown in Tables 1, 2 , 3, and 4. The first and second columns show the configurations of upper and lower levels of the transition, respectively. The third column shows the energies of the 1s 2 np series calculated using quantum defect theory. The transition probabilities calculated using WBEPM for the transitions mentioned above are in the fourth column. In column five, corresponding probabilities published by other scientists are given. Table 1 gives the transition probabilities for the transitions 1s 2 np to 1s 2 2s up to n = 30. Our transition probability results are close to corresponding NIST values; the average difference between is 3.6%, whereas the same difference between Lingard and NIST values is 6.3%. The results of this work are closer to NIST data compared to Zheng's results [11] . Table 2 gives the transition probabilities for transitions 1s 2 np-1s 2 3s. The results of this work are in good agreement with the NIST values, except for one transition probability for 1s 2 4p-1s 2 3s. Due to Cooper's minimum effect expected a large difference between the two values.

Zheng's results [11] for the transitions 1s 2 5p-1s 2 3s and 1s 2 6p-1s 2 3s have large deviations. Table 3 gives the transition probabilities for transitions 1s 2 np-1s 2 4s. The NIST values for transition probability have uncertainties up to 10% in all transitions except 1s 2 4p-1s 2 4s. Only two transition results have large deviations from published data (1s 2 5p-1s 2 4s and 1s 2 6p-1s 2 4s). The transition probabilities for the transitions 1s 2 7p-1s 2 4s and 1s 2 8p to 1s 2 4s are closer to NIST data. Table 4 gives the transition probabilities for transitions 1s 2 np-1s 2 5s. The transition probability for the transition 1s 2 6p-1s 2 5s, calculated in this work, differs from NIST 

The possible lifetimes of Rydberg levels are determined by using the relation τ i = 1 j A ji to check the reliability of the results. The lifetimes of 2p, 3s, 3p, and 4p Rydberg levels are determined and are shown in Table 5 .

The calculated values of lifetimes are compared with the experimental and theoretical results published in [6, 8, 17, [20] [21] [22] [23] [24] [25] . An acceptable agreement is found between the lifetimes calculated in this work and the previously published results.

The dipole matrix elements are zero or have a minimum value for a particular set(s) of principal quantum numbers of upper and lower levels. It occurs due to the overlapping of positive and negative amplitudes of the wavefunctions of these levels. This situation is referred to as Cooper minimum. Due to this effect, the results of transition probabilities are not in good agreement for some of the transitions. Figure 1a -d show transition probabilities for the transition np-ms (m = 2, 3, 4, 5), respectively. A sharp drop can be seen in the transition probabilities in all figures. It shows that the probability of transitioning to neighboring levels is higher than those with larger principal quantum numbers. As principal quantum number increases, the value of transition probability decreases slower than the initial few levels. It results in large lifetimes for higher energy levels and occurs because the transition probability depends on oscillator strength. According to Fano and Cooper [19] , the oscillator strength as a function of effective principal quantum number (n*) drops extremely rapidly towards a minimum, where a reversal of sign in the R integral occurs; this contributes to most of the strength. The rise from the minimum to the second maximum for higher n* is relatively slow. The second maximum has low oscillator strength and is usually not visible in the complete spectrum.

Lithium is a lighter atom, so a non-relativistic approach is suitable to calculate various spectroscopic properties of such atoms. Weakest bound electron potential theory (a semi-empirical method was employed in this work to calculate transition probabilities of various transitions from series 1s 2 np-1s 2 ms ( 2 ≤ n ≤ 30, 2 ≤ m ≤ 5 ). We present 110 transition probabilities, while Heavens [6] , Lingard and Nelson, Kostelecky & Nieto [9] , Fischer [10] , Zheng et al. [11] , and Celik reported 14, 48, 6, 1, 15, 4, and 35 transition probabilities, respectively. Meanwhile, NIST data shows 26 transition probabilities of these series. The results of this work are in good agreement except for four (1s 2 4p-1s 2 3s, 1s 2 5p-1s 2 4s, 1s 2 6p-1s 2 4s, 1s 2 6p-1s 2 5s) of the 110 transitions, and most of the results are new.

Most of the results of transition probabilities in the literature focus on measurement of transition from 2p to 2s (experimental), and very few results of transition probabilities of lithium are reported yet. Lindgard calculated transition probabilities (theoretically) for some series of transitions in lithium for n ≤ 12. we extended the calculation up to n ≤ 30. These values of transition probabilities are helpful in the measurement of oscillator strengths and polarizabilities of Li I. The graphs of transition probabilities show a sharp decrease towards a minimum, followed by a slight rise and eventually a slight decrease. It was observed that the dipole matrix elements have negative values for a few lower transitions that affect oscillator strengths and transition probabilities. It causes transition probabilities to have a minimum known as Cooper minimum. Cooper's minimum effect was first observed in a study of photoionization cross-section of alkali atoms. The calculations are based on the weakest bound electron potential model, which treats lithium as hydrogen-like, but hydrogen does not show zero dipole moment in photoionization. The Cooper minimum reminds us that Weakest bound electron potential model is an approximation; even though lithium is approximately hydrogen-like, it has two additional electrons in the inner shell.

Quantum defect theory was used to calculate quantum defects in principal quantum numbers, which are then used to calculate energies of Rydberg series 1s 2 np, and 1s 2 ns. A generalized formula of quantum defect for each of the series is developed. A computer program in python was developed and used to find the values of the matrix element; this part of the calculation makes our method slightly different from WBEPM theory. A quadratic equation can calculate the radii of the levels in ns and np series in 'n' , which was developed to check the program's reliability. The matrix element, the energies, and quantum defects are used to calculate transition probabilities. For further verification of the results of this work, the lifetimes of some of the Rydberg levels are determined by using calculated transition probabilities. The results are compared with the previously published data, and a good agreement between the two results indicates the reliability of this work. (17) Beam-gas-laser spectroscopy [22] 30.46 Coulomb approximation [26, 27] 30 Central field approximation [6] 30.02 Realistic potentials [14] 30.32 Coulomb approximation [8] 3p 30,925.4 215.578 203 (8) Level-crossing spectroscopy [28] 209 Central field approximation [6] 182 (6) Level crossing spectroscopy [28] 212 Realistic potentials [17] 216.4 Coulomb approximation [8] 4p 36,469.6 416.997 364 Central field approximation [6] 383.9

Realistic potentials [17] 402.6 Coulomb approximation [8] 

Modeling escaping behavior from the herd in different ecological interactions

Global dynamics of alcoholism epidemic model with distributed delays

A robust study of a piecewise fractional order COVID-19 mathematical model

The abundance of lithium and convective mixing in stars of type

Radiative transition probabilities of the lower excited states of the alkali metals

Radiative transition probabilities of the 1s 2p 2 2 P and 1s 2p 2 2 D states of the lithium isoelectronic sequence

Transition probabilities for the alkali isoelectronic sequences Li I, Na I, KI, Rb I, Cs I, Fr I. At Data Nucl Data Tables

Analytical wave functions for atomic quantum-defect theory

Breit-Pauli energies

Transition probability of lithium atom and lithiumlike ions with weakest bound electron wave functions and coupled equations

The calculation of transition probabilities between individual lines for atomic lithium

Accurate atomic transition probabilities for hydrogen, helium, and lithium

Mean life measurements of 2p 2 p 0 levels in the lithium sequence

Accurate lifetime measurements of the lowest 2 P 1/2 states in neutral lithium and sodium

Weakest bound electron potential model theory

Lifetimes of alkali-metal-atom Rydberg states

Radiative lifetime of the 2 P state of lithium

Spectral distribution of atomic oscillator strengths

Accurate time-resolved laser spectroscopy on lithium atoms

Precise atomic radiative lifetime via photoassociative spectroscopy of ultracold lithium

Precision lifetime measurements on alkali atoms and on helium by beam-gas-laser spectroscopy

Lifetime measurements of Lithium D-states

Radiative lifetimes of excited Li I states

Accurate variational calculations of energies of the 2 2 S, 2 2 P, and 3 2 D states and the dipole, quadrupole, and dipolequadrupole polarizabilities and hyperpolarizability of the lithium atom

Lifetimes and polarizabilities of low lying lithium S, P and D states

Level-crossing study of the hyperfine structure of lithium

Hyperfine structure of the 3 2 P and 4 2 P states of lithium and lifetime of the 3 2 P state

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Not applicable.

Not applicable.

Not applicable.

Ethics approval and consent to participate Not applicable.

Not applicable.

The authors declare no competing interests.Received: 26 November 2021 Accepted: 8 March 2022