# Calculation of the Rydberg Energy Levels for Francium Atom

Authors: Huang Shizhong; Sun Qiufeng
Journal: Physics Research International
 (2010)

---

## Abstract

Based on the weakest bound electron potential model theory, the Rydberg energy levels and quantum defects ofnsS21/2(n=8–50), ndD23/2(n=6–50), and ndD25/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.


---

## Body

## 1. Introduction

Much work [1–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 7p2P1/2,3/20 and 8p2P1/2,3/20 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–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 initiocalculation 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 alkali-metal 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 5f5/2 ground state in francium-like ion Th3+ 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 many-body perturbation theory was applied to study energies of the 7s, 7p, 6d, and 5f 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 many-body perturbation theory, relativistic single-double approximation, pseudo-relativistic Hartree-Fock method, relativistic third-order and all-order methods, Dirac-Fock + core-polarization’ 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 1s22s22p63s23p63d104s24p64d104f145s25p65d106s26p67s. Apossible 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 1s22s22pnsP32,1,00 (n=3–50) and 1s22s22pndF33,20 (n=3–50) spectrum series for Carbon I atom by the weakest bound electron potential model theory (WBEPMT) [21–23]. In this paper, the WBEPMT is utilized to investigate the Rydberg energy levels and the quantum defects of nsS1/22 (n=8–50), ndD3/22  (n=6–50) and ndD5/22 (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.


## 2. Theory of WBEPM

In the weakest bound electron potential model theory (WBEPMT) [21–23], the term energy of an atom system is expressed as(1)T=Tlimit+E,
where Tlimit 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)(2)12d2Rdr2+1rdRdr+[E-V(r)-l(l+1)2r2]R=0.
The potential that suffered by the WBE is modeled as [23](3)V(r)=-Z′r+k(k+1)+2kl2r2,
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(4)12d2Rdr2+1rdRdr+[E+Z′r-l′(l′+1)2r2]R=0,
with(5)l′=l+k.The solution to (4) is known as(6)Rnl(r)=Aexp⁡(-Z′rn′)rl'Ln-l-12l'+1(2Z′rn′),E=-R(Z′n′)2,
where n′=n+k, A is the normalization constant, and R the Rydberg constant. Let(7)Z′n′=Znetn-δn,
where Znet is the net charge number of the atomic kernel and Znet=1 for atoms, δn is the quantum defect, the term energy described by (1) is turned into(8)T=Tlimit+E=Tlimit-R(Znetn-δn)2.According to Martin’s formula about quantum defect [24], δn can be expressed as(9)δn=∑i=03ai(n-δ0)2i,
where δ0 is the quantum defect of the lowest energy level in a given series, ai (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.


## 3. Results and Discussion

By using the experimental data of the first four lowest energy levels of a given spectrum series, parametersai are fitted by least-square method, the results are listed in Table 1.Table 1
Spectral coefficients of the four energy series for francium atom by fitting the experimental values.
a0a1a2a3ns2S1/25.070850.32392110.454696.28960nd2D3/23.423660.101560−0.18468−1.68691nd2D5/23.402620.133411−0.117487−0.89660The Rydberg energy levels of the spectrum series ofnsS1/22 (n=8–50), ndD3/22 (n=6–50), and ndD5/22 (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.Table 2
Comparison of the calculated and measured values ofnsS1/22 (n=8–50) energy levels series for francium atoms.
nTcalTexpDiff.δcalnTcalTexpDiff.δcal819739.9819739.980.005.106703032672.2932672.290.005.07131925671.0025671.000.005.089983132685.645.071281028310.61728310.6170.005.082633232697.545.071251129718.90929718.9090.0005.078833332708.195.071221230559.49030559.5040.0145.076623432717.745.071201331101.51831101.5390.0225.075233532726.365.071171431471.44331471.4650.0245.074293632734.165.071151531735.15831735.1820.0215.073633732741.235.071141631929.76831929.7890.0145.073153832747.675.071121732077.47832077.4920.0165.072793932753.555.071101832192.23532192.2510.0185.072504032758.935.071091932283.16232283.1800.0165.072284132763.865.071082032356.42832356.4440.0115.072104232768.405.071072132416.32932416.3400.0085.071954332772.595.071062232465.92932465.9370.0105.071834432776.465.071042332507.4632507.470.005.071734532780.045.071042432542.5932542.590.005.071644632783.365.071032532572.5532572.560.015.071564732786.455.071022632598.3332598.340.015.071504832789.335.071012732620.6632620.670.015.071444932792.015.071002832640.1432640.140.005.071395032794.515.071002932657.2232657.220.005.07135From Tables2, 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.Table 3
Comparison of the calculated and measured values ofndD3/22  (n=6–50) energy levels series for francium atoms.
nTcalTexpDiff.δcalnTcalTexpDiff.δcal616229.8716229.870.003.430352932681.1132681.110.003.42381724244.0324244.030.003.428873032693.5032693.500.003.42379827600.65727600.6570.0003.427313132704.573.42378929316.49729316.4970.0003.426303232714.493.423781030309.96430309.9620.0023.425643332723.423.423771130936.33030936.3250.0053.425193432731.493.423761231356.51231356.5060.0063.424883532738.813.423751331652.00731652.0000.0073.424653632745.463.423751431867.68731867.6820.0053.424483732751.533.423741532029.91032029.9090.0013.424353832757.083.423741632154.98932154.9790.0103.424253932762.173.423731732253.45532253.4490.0063.424174032766.853.423731832332.35632332.3540.0023.424104132771.153.423731932396.55332396.5520.0013.424054232775.133.423722032449.48432449.4830.0013.424004332778.813.423722132493.6432493.640.003.423974432782.223.423722232530.8632530.860.003.423944532785.393.423712332562.5232562.520.003.423914632788.343.423712432589.6832589.680.003.423884732791.083.423712532613.1532613.150.003.423864832793.653.423712632633.5732633.570.003.423854932796.043.423702732651.4432651.440.003.423835032798.293.423702832667.1832667.180.003.42382Table 4
Comparison of the calculated and measured values ofndD5/22 (n=6–50) energy levels series for francium atoms.
nTcalTexpDiff.δcalnTcalTexpDiff.δcal616429.6416429.640.003.414762932681.3932681.390.003.4028724332.9324332.930.003.410283032693.7532693.750.003.40280827645.37327645.3730.0003.407713132704.793.40279929341.81729341.8170.0003.406213232714.693.402781030325.61130325.6050.0063.405283332723.603.402771130946.65130946.6430.0083.404673432731.663.402761231363.66731363.6550.0123.404253532738.963.402751331657.16831657.1550.0133.403943632745.603.402741431871.52931871.5140.0153.403713732751.653.402731532032.84732032.8210.0263.403543832757.193.402731632157.28332157.2740.0093.403403932762.273.402721732255.28132255.2750.0063.403304032766.943.402721832333.83432333.8270.0073.403214132771.243.402711932397.76532397.7610.0043.403144232775.213.402712032450.49032450.4880.0023.403084332778.883.402702132494.4832494.480.003.403034432782.293.402702232531.5732531.570.003.402994532785.453.402692332563.1332563.130.003.402954632788.403.402692432590.2032590.200.003.402924732791.143.402692532613.6032613.600.003.402894832793.703.402682632633.9732633.970.003.402874932796.093.402682732651.8032651.790.013.402855032798.333.402682832667.4932667.490.003.40283Therefore, 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 ofnsS1/22 (n=8–50), ndD3/22 (n=6–50), and ndD5/22 (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 (22dD5/22) 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.



---

*Source: 203497-2010-12-16.xml*