Theoretical energy level spectra and transition data for 4p{}^{6}4d, 4p{}^{6}4f and 4p{}^{5}4d{}^{2} configurations of W{}^{37+} ion

Theoretical energy level spectra and transition data for 4p4d, 4p4f and 4p4d configurations of W ion


The ab initio quasirelativistic Hartree-Fock method developed specifically for the calculation of spectral parameters of heavy atoms and highly charged ions was applied to determine atomic data for tungsten ions. The correlation effects were included by adopting configuration interaction method. The Breit-Pauli approximation for quasirelativistic Hartree-Fock radial orbitals was employed to take into account relativistic effects. The energy level spectra, radiative lifetimes, Lande factors were calculated for the , and configurations of W ion. The atomic data, namely, the transition wavelengths, spontaneous emission rates and oscillator strengths for the electric dipole, electric quadrupole and magnetic dipole transitions among and within the levels of these configurations are tabulated.




[cor1]Corresponding author.


1 Introduction

Due to its unique features, the metallic tungsten is widely employed in devices which operate at extremely high temperatures. These properties facilitate adoption of tungsten in fusion reactors including internationally developed ITER tokamak [[1], [2]]. High temperatures within fusion plasma cause tungsten to evaporate. Furthermore, the tungsten atoms in plasma are ionized to very high degrees. This process has a negative consequences on plasma since it makes plasma temperatures to drop significantly.

To perform plasma processes modeling and diagnostics one needs huge amount of spectroscopic data describing various ionization stages of tungsten and other atoms [[3], [4]]. Atomic spectroscopy can help to determine important properties of plasma [[5], [6], [7]]. Nevertheless, the spectra of the tungsten ions present in fusion plasma devices are still not investigated properly. This conclusion is evident from review [[8]] where all available experimental wavelengths and energy level spectra for ions from W III to W LXXIV are provided.

It is evident from the analysis of [[8]] and references within, that only a few data for the spectra of tungsten ions with open 4d- shells are available. In present work we investigate a relatively simple system, namely W. This ion has only one valence 4d- electron in its ground state. It follows from [[8]] that only 11 excited levels from the possible 48 ones for the configurations under consideration have been determined experimentally. On the other hand, a theoretical study of energy level spectra was presented in [[9]]. However, this work suffers from the fact that the data are provided just for some fifty percent of existing levels representing the configurations under consideration. Furthermore, the transition data were determined for the highly excited levels only.

In present work we employ the quasirelativistic method [[10], [11], [12]] developed specifically for the ab initio calculation of spectral parameters of heavy atoms and highly charged ions. A preliminary analysis of the suitability of this method for the tungsten ions with open 4d- shell was presented in [[13]]. It was demonstrated that our quasirelativistic approach is appropriate for the inclusion of the relativistic effects, and the further improvement of accuracy of the calculated data can be achieved with inclusion of the correlation effects.

In Section 2 we provide a short description of our calculation method. The more detailed description of the applied approach and the analysis of the convergence within the configuration interaction (CI) approximation applied for the inclusion of correlation effects is presented in [[14]]. Obtained results are discussed in Section 3.

2 Calculation method

2.1 Quasirelativistic approximation

We use a quasirelativistic approximation for ab initio calculations of ion energy spectra. This approach significantly differs from widely used method described in [[15]]. In our method, the radial orbitals (RO) are obtained by solving the quasirelativistic equations of the following form:


The first line of Eq.(1) represents the traditional Hartree-Fock equations, where denotes the exchange part of the potential and represents the direct part of the potential including the interaction of an electron with nucleus and with other electrons. We take into account the finite size of a nucleus within the nucleus potential [[16]]. This allow us to represent the radial orbitals in powers of a radial variable in the nucleus region. Next two terms with the multiplier describe the relativistic correction of the mass-velocity dependence. The last term of Eq.(1) represents the potential of the electron contact interaction with the nucleus. In our approach we include the contact interaction with the nucleus not only for the s- electrons but also some part of it for the p- electrons:


A detailed discussion of the particular features of Eq.(1) is given in [[10], [11]], whereas their solution techniques are described in [[17], [18]].

Concluding the description of applied approximation, we want to articulate the unique features of our quasirelativistic method which significantly differs from widely used approach described in [[15]]. The main differences arise from our adopted set of Hartree-Fock quasirelativistic equations (QRHF) featuring several distinctive properties, namely:

  1. No statistical potentials are used. There are only conventional self-consistent field direct and exchange potentials in QRHF.

  2. The finite size of the nucleus is taken into account to determine the potential while solving the quasirelativistic equations [[16], [18]].

  3. The mass-velocity term is divided into two parts - the direct potential and the exchange one.

  4. The contact interaction term contains only the nucleus potential derivative in the numerator. There are no two-electron potentials in the numerator.

  5. Only the direct part of the potential is included into the denominator of the contact interaction with a nucleus term.

  6. The contact interaction with a nucleus is taken into account not only for s-electrons, but also for p-electrons with some additional corrections made [[11]].

One must look in [[10], [11]] for the most complete description of the way to determine QRHF equations employed in this work. The methods to calculate the energy level spectra were discussed extensively in [[12]].

For the energy level spectra calculation, we include all two-electron interactions in the same way as it is done in conventional Breit-Pauli approximation. This similarity makes it possible to apply widely used code mchf breit [[19]] for the angular integration of the Breit-Pauli Hamiltonian matrix elements. We employ the computer program mchf mltpol [[20]] to determine the matrix elements of transition operators along with the code mchf lsjtr [[21]] which has been adopted for use with the quasirelativistic radial orbitals.

2.2 Correlation effects

We use configuration interaction (CI) approximation to include electron correlation effects. In present work, depending on the type of one-electron state, two sorts of radial orbitals are exploited in the CI expansion of wavefunctions. For the configurations under consideration (adjusted configurations), we employ one-electron radial orbitals determined from solutions of Eq.(1). This type of RO was applied for all electrons with and . We adopt the transformed radial orbitals (TRO) to describe virtually excited electrons included in CI expansion.

The method of transformed radial orbitals to account for the configuration interaction effects was developed in [[22]] and successfully applied for the non-relativistic radial orbitals (see, e.g. [[23], [24], [25], [26]]). The comparison of TRO properties with those of the solutions of the multiconfigurational Hartree-Fock-Jucys equations [[27]] demonstrates their close similarity [[28], [29]]. Such a way to include correlation effects within our quasirelativistic approach was successfully applied for the studies of the tungsten spectral parameters in [[30]].

In current work the TRO were employed to describe the admixed configurations, having virtually excited electrons with the principal quantum number and with all allowed values of orbital quantum number . These orbitals are determined from the RO of the configurations under consideration by ensuring the effective inclusion of the correlation effects. For the present calculations, we employ the form of TRO with two free variational parameters and :


Here the multiplier ensures the normalization of determined TRO, the first term in parenthesis performs the necessary transformation, and the second term provides the orthogonality of all employed radial orbitals. The parameters and are varied in order to ensure the maximum of averaged energy correction to the energy of the adjusted configuration in the second order of perturbation theory:


We apply the same set of radial orbitals to describe both the even and odd configurations. This approach enables us to avoid any problems caused by the non-orthogonality of RO in calculation of the electron transition parameters.

As a consequence of adopted TRO, the basis of radial orbitals and the number of possible admixed configurations increases rapidly. Therefore one needs to perform a selection of the admixed configurations by including only those which have the greatest influence on the configurations being adjusted. As a selection criteria we apply the mean weight of the admixed configuration within the CI expansion of the adjusted configuration wavefunction determined in the second order of perturbation theory:


We include only those configurations which have larger than the specified small parameter in our calculation of energy levels and transition parameters. The methods and computer codes to calculate the mean weight from Eq. (5) have been described in [[31], [32]]. Such a way to select the interacting configurations in the non-relativistic approximation was sucessfully exploited (see, e.g., earlier mentioned [[23], [24], [25], [26]] and many others).

The configuration state functions (CSFs) of the selected configurations are used to form the energy operator matrices. By diagonalizing these Hamiltonian matrices, we can determine the energy eigenvalues and eigenfunctions for the investigated configurations and exploit them for calculation of transition parameters. More details of employed approach one can find in [[12]] and references therein, whereas some examples of application of our method are presented in [[30], [33]].

In present work we have selected only those admixed configurations which have the mean weights in the wave functions of the configurations under consideration larger than . This procedure has allowed us to reduce the number of admixed configurations in the CI expansion of the wavefunctions to 102 for even configurations and to 217 for odd configurations. So here we use an approximation denoted as F in [[14]]. It was demonstrated in [[14]] that any further expansion of the admixed configurations basis does not lead to significant improvement in the calculated energy level spectra accuracy.

We have determined energy eigenvalues and eigenfunctions by diagonalizing Hamiltonian matrices obtained in the -coupling scheme in Breit-Pauli approximation. Using these results, we have calculated the transition parameters. We have computed data for the electric dipole (E1), electric quadrupole (E2) and magnetic dipole (M1) operators for the transitions both between the configurations under consideration and within these configurations. The emission transition rates of all above mentioned types were used to determine the radiative lifetimes of excited states.

3 Results and discussion

We present the energy level spectra with the percentage contribution of eigenfunctions together with Lande factors and the radiative lifetimes for the investigated configurations in Table 2. It was demonstrated in [[14]] that deviations of our theoretical level energy values from the experimental spectra data do not exceed in most cases. Consequently, we can conclude that these discrepancies in transition energies do not affect significantly the accuracy of calculated values of the transition rates or oscillator strengths. In Table 1 we compare our calculated energy levels with data obtained from the completely relativistic calculations RELAC [[9]] and with available experimental data [[7]] presented in compilation [[8]]. This evaluation favorably demonstrates a better agreement of our data with the experimental level energies, especially for highly-excited states, comparing to the data from relativistic calculations RELAC. This is caused mainly by considerably improved inclusion of correlation effects.

Such a good agreement with experimental data leads to conclusion that our quasirelativistic approach adequately accounts for the relativistic effects even for highly-ionized heavy atoms. It is worth to mention a very good agreement of theoretical fine-structure level splitting for the ground term 4p4d D pointing to close spin-orbit interaction values obtained in our work and in fully relativistic calculations [[9]].

As it was mentioned before, all our calculations were performed in -coupling scheme. It is well established fact that the spin-orbit interaction induces a very strong mixing of -terms for highly-charged ions. Therefore in Table 1 we present only the level numbers (indices) taken from Table 2 where one can find a corresponding CI wavefunction expansion. It is easy to notice from the CSFs contributions presented in Table 2 that the strong mixing occurs both for the -terms of the same configuration and for terms of different configurations.

A well-established case of a strong interaction between the and configurations (see [[34]]) can be confirmed from the analysis of the wave function percentage contributions. For example, a contribution from the configuration for the level in Table 2 is merely . Moreover, the largest contribution for the level comes from the configuration although the level itself clearly must be attributed to the configuration .

In Table 3 we present the transition wavelengths , the emission transition rates and the weighted oscillator strength values . We use only the level indices and the total angular momentum from Table 2 to describe both the initial and final levels along with the type of transition operator (E1, E2 or M1) applied to calculate the transition parameters. We present the transitions from the excited level to lower-lying levels in decreasing order of their rates. We limit the number of presented transitions by removing those having the transition rate values by 3 orders of magnitude lower comparing to the strongest transition from the particular energy level. This condition would remove completely the E2 transitions. Nonetheless, we have placed some most important E2 transitions in Table 3.

The magnetic dipole M1 transitions within configuration were computed using a completely relativistic approximation in [[35]]. As one can see from Table 3, their value of transition rate is very close to the value of from our calculation.

From analysis of data presented in Table 3, we can notice that the electric dipole transition rates are the most strong ones for all excited configuration levels when the total angular momentum . Furthermore, the strong E1 transitions define the values of radiative lifetimes for these levels presented in Table 2. The only exemption from this rule demonstrates the level with total angular momentum . Even if the E1 transitions from this level have the largest rates , the values of two M1 transition rates are of the same order of magnitude as for the E1 transitions. This causes the significant decrease of the radiative lifetime value for this energy level.

The electric dipole transitions from the excited levels with are forbidden because the largest value of total angular momentum of the ground state . Consequently, the radiative lifetimes for these levels are determined by the magnetic dipole transitions since the E2 transition rates are significantly smaller.

This work, partially supported by the European Communities under the contract of Association between EURATOM/LEI FU07-CT-2007-00063, was carried out within the framework of the European Fusion Development Agreement.
Term [[8]] Level Name [[9]] EXP [[8]] QRHF RELAC [[9]]
1 4p4d D (4p)(4d) 1.5 0 0 0
2 4p4d D (4p)(4d) 2.5 1546 1538 1537
3 4p4d F (3/2,2) 0.5 12276 12311
4 4p4d F (3/2,2) 1.5 12276 12312
5 (4p)(4p)(4d) 2.5 12591 12569
6 (4p)(4p)(4d) 3.5 12752 12720
10 (4p)(4p)(4d)(4d) 4.5 13922 13861
13 (4p)(4p)(4d)(4d) 5.5 14255 14201
16 (4p)(4p)(4d)(4d) 3.5 14631 14646
17 (4p)(4p)(4d)(4d) 2.5 14680 14691
18 (4p)(4p)(4d)(4d) 1.5 14756 14794
19 (4p)(4p)(4d)(4d) 4.5 14853 14861
20 4p4d G (3/2,4) (4p)(4p)(4d)(4d) 2.5 15085 15220 15220
21 (4p)(4p)(4d) 5.5 15500 15443
22 4p4d F (3/2,3) (4p)(4p)(4d)(4d) 1.5 15426 15618 15680
24 (4p)(4p)(4d) 3.5 15671 15649
25 (4p)(4p)(4d) 4.5 15854 15851
28 (4p)(4p)(4d) 3.5 16561 16567
30 4p4d F (3/2,4) (4p)(4p)(4d) 2.5 17315 17535 17630
32 4p4f F (4p)(4f) 2.5 17581 17752 17840
33 4p4f F (4p)(4f) 3.5 17696 17889 17951
34 4p4d F (1/2,2) (4p)(4p)(4d) 2.5 20145 20319 20396
44 4p4d D (1/2,2) 1.5 21843 22069
46 4p4d G (1/2,4) (4p)(4p)(4d)(4d) 3.5 23255 23274 23488
47 (4p)(4p)(4d)(4d) 1.5 23309 23525
48 (4p)(4p)(4d)(4d) 2.5 23369 23347
Table 1: Comparison of calculated energy levels (in 100 cm) of W with available experimental data (EXP) from [[8]] and relativistic calculations RELAC from [[9]]. Level identifications from both [[8]] and [[9]] are used where available, and level numbers are from Table 2.

Table 1. Energy levels and radiative lifetimes of and configurations of W.

Throughout the table we present all the eigenfunction components which have their percentage contributions higher than

The energy level number (index)
The level energy in cm
The Lande factor . The notation a(b) means
The radiative lifetime in  s
The total angular momentum
Contribution The percentage contribution of CSF in level eigenfunction

Table 2. Electron transitions among and within the , and configurations of W.

The transitions from the upper level to the lower ones are presented in the descending order of the transition rate values . The number of presented transitions is limited to those having the transition rate values less than 3 orders of magnitude smaller comparing to the strongest transition from the particular energy level.

The total angular momentum of the initial level
The initial energy level number
The total angular momentum of the final level
The final energy level number
Type The electron transition type
(Å) The transition wavelength in Å
(s The emission transition rate in s. The notation aEb means
The weighted oscillator strength . The notation aEb means
1 0 0.80 1.5 98 4d D
2 1538 1.20 2.55(4) 2.5 98 4d D
3 12311 1.00 8.44(1) 1.5 42 4p4d(F) D +23 4p4d(F) F +10 4p4d(D) D
4 12312 0.12 2.44(1) 0.5 84 4p4d(F) D +11 4p4d(D) P
5 12591 1.03 3.83(1) 2.5 36 4p4d(F) F +18 4p4d(D) D +12 4p4d(F) G +10 4p4d(F) D
6 12752 1.03 1.75(1) 3.5 32 4p4d(F) G +31 4p4d(F) G +23 4p4d(D) F +7 4p4d(F) F
7 13284 1.36 5.37(0) 1.5 28 4p4d(F) P +21 4p4d(S) P +15 4p4d(P) D
8 13636 1.30 1.95(1) 2.5 32 4p4d(F) D +22 4p4d(P) P +15 4p4d(P) D
9 13915 1.29 2.88(1) 3.5 35 4p4d(F) F +27 4p4d(P) D +16 4p4d(F) D +10 4p4d(D) F
10 13922 1.18 3.04(4) 4.5 70 4p4d(F) G +15 4p4d(F) G +13 4p4d(F) F
11 14052 1.09 1.29(1) 1.5 40 4p4d(F) D +14 4p4d(F) D +14 4p4d(D) D +11 4p4d(D) P
12 14104 1.23 4.29(1) 2.5 40 4p4d(P) P +17 4p4d(F) F +15 4p4d(G) F
13 14255 1.17 2.68(6) 5.5 57 4p4d(G) H +42 4p4d(F) G
14 14290 1.06 1.15(4) 3.5 29 4p4d(F) G +15 4p4d(F) F +15 4p4d(F) F +15 4p4d(G) G +14 4p4d(G) F
15 14408 1.66 4.37(2) 0.5 33 4p4d(P) P +29 4p4d(P) S +19 4p4d(D) P
16 14631 1.17 3.20(1) 3.5 20 4p4d(P) D +13 4p4d(D) F +13 4p4d(G) G +12 4p4d(G) F +11 4p4d(F) D
17 14680 1.22 1.60(1) 2.5 36 4p4d(P) D +25 4p4d(P) D +15 4p4d(D) D
18 14756 1.39 2.77(0) 1.5 30 4p4d(P) S +25 4p4d(P) D +17 4p4d(P) P
19 14853 1.13 1.47(5) 4.5 38 4p4d(G) G +24 4p4d(F) F +18 4p4d(F) G +17 4p4d(G) H
20 15220 0.88 6.50(3) 2.5 41 4p4d(G) F +28 4f F +13 4p4d(F) F
21 15500 1.20 1.70(4) 5.5 57 4p4d(F) G +42 4p4d(G) H
22 15618 0.85 2.16(3) 1.5 32 4p4d(F) F +22 4p4d(P) D +12 4p4d(F) D
23 15619 1.21 1.60(0) 2.5 24 4p4d(D) D +18 4p4d(P) P +15 4p4d(F) D +13 4p4d(P) D
24 15671 1.22 5.32(2) 3.5 34 4p4d(D) F +29 4p4d(P) D +13 4fF
25 15854 1.17 1.20(4) 4.5 38 4p4d(F) F +30 4p4d(G) G +18 4p4d(F) G +13 4p4d(G) H
26 15972 1.08 6.58(3) 0.5 65 4p4d(D) P +14 4p4d(P) P +12 4p4d(P) S
27 16350 0.76 1.49(3) 0.5 52 4p4d(P) P +24 4p4d(P) D +18 4p4d(P) S
28 16551 1.18 4.48(2) 3.5 17 4f F +16 4p4d(F) F +15 4p4d(D) F +14 4p4d(F) D +12 4p4d(P) D
29 16897 1.27 5.82(1) 1.5 50 4p4d(S) P +14 4p4d(P) D
30 17535 1.09 1.02(3) 2.5 29 4p4d(F) D +15 4p4d(F) F +13 4f F
31 17553 1.36 1.58(3) 1.5 20 4p4d(S) P +17 4p4d(P) S +16 4p4d(P) P +13 4p4d(D) D +12 4p4d(P) D
32 17752 0.93 1.21(3) 2.5 49 4f F +18 4p4d(G) F
33 17889 1.12 1.72(3) 3.5 63 4f F +18 4p4d(G) G
34 20319 0.77 8.80(4) 2.5 39 4p4d(F) G +28 4p4d(F) F +15 4p4d(D) F
35 21314 1.02 7.13(1) 3.5 46 4p4d(F) G +19 4p4d(F) G +14 4p4d(F) F
36 21387 1.07 1.70(2) 2.5 18 4p4d(F) F +16 4p4d(D) F +14 4p4d(F) G +14 4p4d(P) D +12 4p4d(F) D
37 21478 1.07 9.10(3) 1.5 37 4p4d(P) D +16 4p4d(F) D
38 21547 1.96 1.39(3) 0.5 41 4p4d(P) P +33 4p4d(P) S +13 4p4d(P) P
39 21555 0.43 9.08(4) 0.5 56 4p4d(P) D +21 4p4d(P) P +13 4p4d(S) P
40 21754 1.04 1.39(2) 4.5 45 4p4d(G) H +19 4p4d(G) G +16 4p4d(F) G
41 22069 0.88 3.03(4) 1.5 43 4p4d(F) D +17 4p4d(G) D +12 4p4d(F) F
42 22794 1.29 2.73(3) 3.5 45 4p4d(F) D +32 4p4d(F) F +18 4p4d(F) F
43 23116 1.15 3.88(4) 2.5 34 4p4d(F) D +12 4p4d(P) D +10 4p4d(F) D
44 23159 1.20 4.30(3) 1.5 33 4p4d(D) P +26 4p4d(D) D +12 4p4d(P) D +11 4p4d(P) P
45 23221 1.11 1.30(2) 4.5 31 4p4d(F) G +23 4p4d(G) H +17 4p4d(F) G +15 4p4d(F) F +11 4p4d(G) G
46 23274 1.06 5.17(4) 3.5 53 4p4d(G) F +29 4p4d(G) G
47 23309 1.48 4.69(4) 1.5 52 4p4d(P) P +20 4p4d(P) S +11 4p4d(P) P
48 23369 1.11 4.68(1) 2.5 35 4p4d(D) F +24 4p4d(P) D +17 4p4d(P) D
49 24058 0.67 3.13(2) 0.5 75 4p4d(S) P +13 4p4d(P) D
Table 2: (continued)
Type ( Å)  (s)
2.5 2 1.5 1 M1 650.2 3.93E04 1.49E05
1.5 1 E2 650.2 1.89E01 7.19E09
1.5 3 1.5 1 E1 81.23 1.15E09 4.56E03
2.5 2 E1 92.83 3.38E07 1.75E04
0.5 4 1.5 1 E1 81.22 4.09E09 8.09E03
2.5 5 1.5 1 E1 79.42 2.37E09 1.34E02
2.5 2 E1 90.47 2.47E08 1.82E03
3.5 6 2.5 2 E1 89.18 5.72E07 5.45E04
1.5 7 2.5 2 E1 85.14 1.78E08 7.74E04
1.5 1 E1 75.28 8.25E06 2.80E05
2.5 8 1.5 1 E1 73.34 5.13E07 2.48E04
2.5 2 E1 82.66 3.11E04 1.91E07
3.5 9 2.5 2 E1 80.80 3.47E09 2.72E02
4.5 10 3.5 6 M1 854.6 3.29E04 3.60E05
3.5 6 E2 854.6 2.13E00 2.33E09
1.5 11 1.5 1 E1 71.17 5.76E09 1.75E02
2.5 2 E1 79.91 1.97E09 7.56E03
2.5 12 2.5 2 E1 79.58 2.31E09 1.31E02
1.5 1 E1 70.90 2.44E07 1.11E04
5.5 13 4.5 10 M1 2999 3.68E02 5.95E06
3.5 6 E2 665.1 5.94E00 4.73E09
3.5 14 2.5 2 E1 78.42 4.94E04 3.65E07
2.5 5 M1 588.5 2.32E04 9.64E06
3.5 6 M1 650.0 1.42E04 7.21E06