Calculation of the Rydberg Energy Levels for Francium Atom

Based on the weakest bound electron potential model theory, the Rydberg energy levels and quantum defects of nsS1/2(n = 8–50), ndD3/2(n = 6–50), and ndD5/2(n = 6–50) spectrum series for francium atom are calculated. The calculated results are in excellent agreement with the 74 known experimentally measured levels (the absolute difference is less than 0.03 cm−1) and 58 energy levels for highly excited states are predicted.


Introduction
Much work [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19] has been carried out to investigate the energy structure in excited states of francium atom, the heaviest alkali metal atom.The experimental investigations of the francium spectrum were pioneered by Liberman et al. in 1978 [1].Later on, the energy levels of 7p 2 P 0 1/2, 3/2 and 8p 2 P 0 1/2, 3/2 were investigated by Doppler-free laser spectroscopy techniques in 1986 and 1987 [2,3].The ns (n = 23, 25-27, 29-31) and nd (n = 22-33) Rydberg levels were observed by two-step laser excitation and detection techniques also in 1987 [4], the ns (n = [10][11][12][13][14][15][16][17][18][19][20][21][22] and nd (n = 8-20) series were determined by stepwise laser excitation in 1989 and 1990 [5,6], and the 9s level was observed in 1996 [7].Meanwhile, the ionization potential was determined in 1988 [8] and refined in 1989 [5].More recently, an interferometric method was used to improve the accuracy of the 7S-7P transition frequency of francium by 1 order of magnitude in 2009 [9].The theoretical investigation of the energy structure of francium could be outlined as follows.In 1995, many-body perturbation theory in screened Coulomb interaction was used to calculate energy levels, E1 transition amplitudes, and the parity-nonconserving E1 amplitude of the 7s-8s transition in francium by Dzuba et al. [10].In 1998, Ritz formalism was utilized to predict the francium Rydberg ns, np, and nd levels up to n = 30 [12].In 1999, removal energies of the n = 7-10 states and hyperfine constants of the n = 7 and 8 states in francium were calculated by Safronova et al. [11], the calculations were based on the relativistic single-double approximation in which single and double excitations of Dirac-Hartree-Fock wave functions are included to all orders in perturbation theory.In 2004, pseudorelativistic Hartree-Fock method including core-polarization effects was used to investigate the radiative parameters for electric dipole transitions in the ions Ra II, Ac III, Th IV, and U VI along the francium isoelectronic sequence by Biémont et al. [13].Also in 2004, relativistic Hartree-Fock method with several correlations was used to perform ab initio calculation of isotope shifts for isotopes of francium (from A = 207 to A = 228) by Dzuba et al. [14].In 2005, radiative corrections to E1 matrix elements for ns-np transitions in the alkalimetal atoms lithium through francium are evaluated by Sapirstein and Cheng [15].In 2006, relativistic third-order and all-order methods were used to calculate the scalar and tensor polarizabilities for the 5 f 5/2 ground state in francium-like ion Th 3+ by Safronova et al. [16].In 2007, some E1 transitions in the francium isoelectronic sequence were computed in the "Dirac-Fock + core-polarization" approximation by Migdalek and Glowacz-Proszkiewicz [17], in which the core-valence electron correlation was treated in a semiclassical picture.Also in 2007, relativistic manybody perturbation theory was applied to study energies of the 7s, 7p, 6d, and 5 f states of francium-like ions with nuclear charges Z = 87-100, and also the transition rates, oscillator strengths, and lifetimes of francium-like ions with Z = 87-92 by Safronova et al. [18].Most of the energy levels for neutral Francium have been compiled by Sansonetti in 2007 [19] In all these theoretical works; however, Rydberg energy levels for highly excited states (n > 31) of francium-like atoms have not been efficiently studied.The reason is that no matter what approximation method (including manybody perturbation theory, relativistic single-double approximation, pseudo-relativistic Hartree-Fock method, relativistic third-order and all-order methods, Dirac-Fock + corepolarization' approximation method, and so on ) is used, it is hard to investigate the highly excited states for such a complicated system that the configuration of the ground state is 1s 2 2s 2 2p 6 3s 2 3p 6 3d 10 4s 2 4p 6 4d 10 4 f 14 5s 2 5p 6 5d 10 6s 2 6p 6 7s.A possible approach to solve this problem is to utilize model potential method.In a previous work [20], we have successfully calculated the Rydberg energy levels and the quantum defects of 1s 2 2s 2 2pns 3 P 0 2,1,0 (n = 3-50) and 1s 2 2s 2 2pnd 3 F 0 3,2 (n = 3-50) spectrum series for Carbon I atom by the weakest bound electron potential model theory (WBEPMT) [21][22][23].In this paper, the WBEPMT is utilized to investigate the Rydberg energy levels and the quantum defects of ns 2 S 1/2 (n = 8-50), nd 2 D 3/2 (n = 6-50) and nd 2 D 5/2 (n = 6-50) spectrum series for francium atom.The calculated results are in excellent agreement with the 74 known experimentally measured levels [19] (the absolute difference is less than 0.03 cm −1 ), and 58 energy levels for highly excited states are predicted, which might be significant for experiment guidance.

Theory of WBEPM
In the weakest bound electron potential model theory (WBEPMT) [21][22][23], the term energy of an atom system is expressed as where T limit is the ionization limit corresponding to a given spectrum series, and E is the energy of the weakest bound electron (WBE), which is determined by the radial equation (in atom unit) 1 2 The potential that suffered by the WBE is modeled as [23] V in which Z stands for the number of effective nuclear charge, and k an adjustable parameter and is not necessarily an integer.The first term of (3) represents the Coulomb potential and the second term represents the dipole potential that originated from polarization.Inserting (3) into (2), the radial equation becomes 1 2 with The solution to ( 4) is known as where n = n + k, A is the normalization constant, and R the Rydberg constant.Let where Z net is the net charge number of the atomic kernel and Z net = 1 for atoms, δ n is the quantum defect, the term energy described by ( 1) is turned into According to Martin's formula about quantum defect [24], δ n can be expressed as where δ 0 is the quantum defect of the lowest energy level in a given series, a i (i = 0, 1, 2, 3) are parameters that can all be obtained through the least-square fitting of (8) with the first few experimental data in a given spectrum series.

Results and Discussion
By using the experimental data of the first four lowest energy levels of a given spectrum series, parameters a i are fitted by least-square method, the results are listed in Table 1.
The Rydberg energy levels of the spectrum series of ns 2 S 1/2 (n = 8-50), nd 2 D 3/2 (n = 6-50), and nd 2 D 5/2 (n = 6-50) for francium atom are calculated by (9).The calculated energy levels, quantum defects, and the differences between calculated and measured energy levels [19] are shown in Table 2 to Table 4, respectively.The unit of energy levels is expressed in cm −1 .
From Tables 2, 3, and 4, it can be seen that the calculated results are in excellent agreement with the known experimental data since the absolute difference is less than 0.03 cm −1 .Therefore, by using the experimental data of the first four lowest energy levels for a given spectrum series, the Rydberg energy levels and the quantum defects of ns 2 S 1/2 (n = 8-50), nd 2 D 3/2 (n = 6-50), and nd 2 D 5/2 (n = 6-50) spectrum series for francium atom have been calculated by virtue of the WBEPMT, the calculated results are in excellent agreement with the 74 known experimentally measured levels, one exception (22d 2 D 5/2 ) is noticed, and 58 energy levels for highly excited states are predicted, which might be significant for experiment guidance.This implies that the WBEPMT is not only reliable to investigate the Rydberg energy levels of CI, but also dependable to probe the Rydberg energy levels of FrI.Compared with many-body perturbation theory, relativistic single-double approximation, pseudo-relativistic Hartree-Fock method, relativistic third-order and all-order methods, and Dirac-Fock + core-polarization' approximation method, the advantage of WBEPMT is that the program code, which can be easily written by the Mathematica language, is simpler and the running time is much shorter.The disadvantage of WBEPMT; however, is that the wave functions of the Rydberg states cannot be obtained, which might be settled by combining the WBEPMT with other theories.

Table 1 :
Spectral coefficients of the four energy series for francium atom by fitting the experimental values.

Table 2 :
Comparison of the calculated and measured values of ns 2 S 1/2 (n = 8-50) energy levels series for francium atoms.

Table 3 :
Comparison of the calculated and measured values of nd 2 D 3/2 (n = 6-50) energy levels series for francium atoms.

Table 4 :
Comparison of the calculated and measured values of nd 2 D 5/2 (n = 6-50) energy levels series for francium atoms.