*Chin. Phys B*23 (2014) 093203]

*Chin. Phys B*23 (2014) 093203]

† Corresponding author. E-mail:

Recently, S. Aggarwal [*Chin. Phys. B* ^{2}2p^{5}, 2s2p^{6}, and 2s^{2}2p^{4}3*ℓ* configurations of F-like tungsten. There is no discrepancy for his calculated energies for the levels and the radiative rates for the limited number of E1 transitions, but the reported results for the lifetimes are highly inaccurate. According to our calculations, errors in his reported lifetimes are up to 6 orders of magnitude for several levels. Here we report the correct lifetimes for future comparisons and applications, and also explain the reason for the discrepancies.

Atomic data, including energy levels and oscillator strengths or radiative decay rates, are required for many ions in order to model plasmas and/or for diagnostics, such as to determine elemental density, temperature and chemical composition. Generally, lighter ions (with *Z* ≤ 28) are important for astrophysical applications whereas heavier ones are more important for lasing and fusion plasmas. An important element of particular interest for studies of fusion plasmas is tungsten (W), because apart from being an important constituent of the reactor walls, it radiates at several of its ionization stages. Considering the requirement of the developing ITER project, atomic data for its ions are in high demand. However, it is not any data, but an accurate and complete set of data which is required for modelling and diagnostics.^{[1]} Recently, S. Aggarwal,^{[2]} henceforth to be referred to as SA, has reported results for energy levels, oscillator strengths (*f*-values), radiative rates (*A*-values), line strengths (*S*-values), and lifetimes (*τ*) for F-like W LXVI. For his calculations, the modified version of the GRASP (general-purpose relativistic atomic structure package) code was adopted. This code is known as GRASP0 and is freely available at the website *j j* coupling scheme, and includes further relativistic corrections arising from the Breit interaction and QED (quantum electrodynamics) effects.

SA reported energies and lifetimes for the lowest 60 levels of the 2s^{2}2p^{5}, 2s2p^{6}, and 2s^{2}2p^{4}3*ℓ* configurations of W LXVI, although he included *configuration interaction* (CI) among 359 levels. The additional 299 levels arise from the 2s^{2}2p^{4}4*ℓ*, 2s^{2}2p^{4}5*ℓ* (*ℓ* ≤ 3), 2s2p^{5}3*ℓ*, 2s2p^{5}4*ℓ*, and 2p^{6}3*ℓ* configurations, i.e. 23 in total. However, he listed *A*-values (and other related parameters) for only electric dipole (E1) transitions from the lowest two levels of the ground configuration *τ A*-values are also required for the electric quadrupole (E2), magnetic dipole (M1), and magnetic quadrupole (M2) transitions. Because of the limitations of his data, we notice that while his reported *A*-values are correct, the results for *τ* are not. Therefore, in this communication we report the correct values of the lifetimes.

For consistency, we have employed the same version of the GRASP code as adopted by SA.^{[2]} Similarly, for compatibility we use the same option of *extended average level*, as adopted by him. To make comparisons, we have performed three sets of calculations, namely GRASP1: with 113 levels of 11 configurations, i.e., 2s^{2}2p^{5}, 2s2p^{6}, 2s^{2}2p^{4}3*ℓ*, 2s2p^{5}3*ℓ*, and 2p^{6}3*ℓ*; GRASP2: which includes the same 23 configurations and 359 levels as by SA (see the list in Section 1); and finally GRASP3: which includes a further 142 levels of the 2s^{2}2p^{4}5g, 2s2p^{5}5*ℓ*, 2p^{6}4*ℓ*, and 2p^{6}5*ℓ* configurations. Therefore, GRASP3 calculations include 501 levels from 38 configurations in total. This exercise is performed mainly to (i) demonstrate that the inclusion of CI additional to the basic 11 configurations of GRASP1 makes a nominal impact on the calculations of atomic parameters, including *τ*, and (ii) levels of higher configurations intermix with the lowest 113. There are no (major) differences for energy levels, as may be noted from Table ^{2}2p^{4}4s configurations intermix (see Table 8 of Aggarwal and Keenan^{[3]}). Similarly, there are no differences in *A*-values between our calculations and those of SA. However, we note here that energies obtained by SA with the *Flexible Atomic Code* (FAC) of Gu^{[4]} and listed in his Table 2 are *incorrect* for some of the levels, particularly the last 5, because these differ with his GRASP calculations by nearly 1 Ryd. We have performed similar calculations with FAC and do not observe any such discrepancies. Furthermore, it is clear from Table *τ* for which there are large discrepancies.

The lifetime *τ* of a level *j* is determined as *τ*_{j} = 1/∑_{i} *A*_{ji}, i.e., the summation is over *all* types of transitions from lower levels *i* to higher *j*. In Table *τ* from GRASP1, GRASP2, and GRASP3 calculations, along with those of SA (GRASP4) for a ready comparison. However, results are listed for only the lowest 92 levels, because beyond those the intermixing of the levels of the 2s^{2}2p^{4}4s configuration affect the calculations of *τ*. However, our calculated *τ* for the last two levels of SA, namely (2s^{2}2p^{4}3d) ^{4}F_{3/2} and ^{2}F_{5/2}, not included in our Table ^{−16} and 2.195×10^{−12} s, respectively, and correspond to levels 101 and 106 listed in Table 8 of Aggarwal and Keenan.^{[3]} For these two levels the respective *τ* by SA are 8.65×10^{−10} and 1.00×10^{−11} s, i.e., differing by up to 6 orders of magnitude.

It is clear from Table *τ* within a satisfactory agreement, and the discrepancy, if any, is below 10% — see level 33 for example. Therefore, as stated above, the effect of additional CI included in GRASP2 and GRASP3 calculations is nominal on the determination of atomic parameters. On the other hand, the *τ* results of SA are invariably *higher*, by up to 6 orders of magnitude for several levels – see for example 14, 41, 52, 56, and 59. Since SA included the same CI as in our GRASP2 calculations, his results for *τ* should have been comparable with ours. Unfortunately, this is not the case, and the discrepancy cannot be understood because of the limited listing of *A*-values by SA. In this short communication it is not possible for us to list *A*-values for all transitions, but in Table *A*-values which contribute to over 20% to ∑_{i} *A*_{ji}. This will help us to understand the differences between ours and the *τ* of SA. We also note that for level 60 (2s^{2}2p^{4}3d ^{2}D_{5/2}), several transitions make substantial contributions but are individually below 20% of the sum total. Therefore, for this level those transitions are listed which contribute to over 10%. Finally, all *A*-values for a wide range of transitions are provided in a separate paper.^{[3]}

SA mentioned the calculations for M2 transitions (only in the abstract but not in the text of the paper), and nothing about E2 and M1. Additionally, he presented no results for any other types except E1. However, it is clear from Table *τ* for level 2, the 1–2 M1 transition dominates. Since there is no discrepancy between the two sets of results for this level (see Table *must* have included the contributions of M1 transitions. However, he did not calculate *A*-values for E2 transitions, because for 7–9 levels, their *A*-values clearly dominate, hence the discrepancies between the GRASP2 and GRASP4 results for *τ*. On the other hand, for level 20 the 1–20 M2 and 10–20 E1 transitions contribute almost equally. As there is a discrepancy of over an order of magnitude in *τ* for this level, it is clear that he made no calculations for M2 transitions. However, the 10–20 E1 *A*-value (9.891×10^{10} s^{−1}) should give *τ* = 1.01×10^{−11} s and not 9.36×10^{−10} s, as listed by him. This is because, his calculation is based on *A*-value (1.068×10^{9} s^{−1}) for the 8–20 E1 transition alone. Among the E1 transitions, the one which dominates is 12–20 with *A* = 2.747×10^{10} s^{−1}, not listed in Table

For level 14, the 6–14 E1 *A*-value (1.476×10^{12} s^{−1}) contributes to over 76%, but the *τ* of SA is perhaps based on 3–14 (*A* = 3.80×10^{9} s^{−1}) alone, and hence the discrepancy of four orders of magnitude. Similarly, for level 56 the *A*-value for the 38–56 E1 transition contributes to about half, but the calculation of SA is probably based on the *A*-value for the 8–56 E1 (*A* = 2.446×10^{6} s^{−1}) transition alone. Therefore, it is apparent that his calculations are totally random and hence highly unreliable. It may be worth noting here that W LXVI is not the only ion for which we have found discrepancies, errors and anomalies in his (and the group members) reported work. Several of their calculations have been demonstrated as unreliable — see for example W XL,^{[5]}, W LVIII,^{[6]} and W LXII.^{[7]}

In this work, we have calculated energy levels, radiative rates and lifetimes for the lowest 501 levels/transitions of F-like tungsten, i.e., W LXVI, by adopting the widely used GRASP code. However, results are presented here only for the lowest 92 levels to make comparisons with the recent work of SA,^{[2]} who has reported similar results. His results for energy levels and *A*-values, only for limited levels/transitions, are found to be correct, but the corresponding values of *τ* have large errors, of up to 6 orders of magnitude, for several levels. This is because he has randomly selected transitions for the calculations, which are not always the most contributing ones. Therefore, for the benefit of users (mainly experimentalists) as well as for future workers, we have listed the correct values of *τ*, and have also explained the reasons for the discrepancies.

Since SA^{[2]} reported *A*-values only for E1 transitions from the ground levels, these are insufficient for applications, because a *complete* set of data for *all* transitions is required for any modelling application. Additionally, he did not consider all levels of the calculations as several are missing (see Table ^{[3]}

**Reference**

1 | |

2 | Chin. Phys. 23 093203 |

3 | |

4 | Can. J. Phys. 86 675 |

5 | Can. J. Phys. 92 545 |

6 | Can. J. Phys. 92 1166 |

7 | Chin. Phys. 24 123201 |