A&A 416, 375-381 (2004)
DOI: 10.1051/0004-6361:20031703

Dipole allowed and intercombination transitions in K7+ and Ti10+ ions[*]

R. Das1 - N. C. Deb1 - K. Roy 1 - A. Z. Msezane2


1 - Department of Theoretical Physics, Indian Association for the Cultivation of Science, Jadavpur, Calcutta 700 032, India
2 - Department of Physics and Center for Theoretical Studies of Physical Systems, Clark Atlanta University, Atlanta, GA-30314, USA

Received 4 June 2003 / Accepted 17 October 2003

Abstract
We have calculated 87 fine-structure energy levels of K7+ and Ti10+ belonging to the lowest 14 configurations within the n=3 and n=4 complexes using Program CIV3 of Hibbert. Relativistic effects are included through the Breit-Pauli approximation. New improved results of oscillator strengths, transition probabilities and line strengths for all allowed and intercombination transitions among these fine-structure levels are calculated using large-scale configuration interaction wavefunctions. Lifetimes of some relatively longer-lived levels are also computed. Our results are compared with a wide variety of other calculations and available experimental data. It is noted that most of our results compare more favourably with the measurements than the other available theoretical results.

Key words: atomic data

1 Introduction

Stellar coronal spectra are rich in emission lines from high ionization states of the cosmically abundant elements (Brickhouse 2000), such as C, O, Ne, Mg, Si, Fe and Ni, which radiate in the X-ray band. The X-ray Chandra and XMM-Newton observatories have been producing a wealth of spectroscopic information for a wide range of astrophysical X-ray sources, such as accretion driven binary systems and AGNs (Kaastra 2001; Paerels et al. 2000). Accurate and complete atomic data for all stages of ionization for elements at least up to Zn (Kurucz 2002) are required to infer reliably the properties of many cosmic plasmas and, hence for the understanding of the fundamental issues of astrophysics (Savin 2001). However, the field of high-Z atoms poses many open questions and still requires a substantial effort in the atomic physics (Behringer 2000) and some of the available atomic data are riddled with uncertainties or are incomplete. Hence, the present investigation.

Mg-like ions have been observed in rocket-borne measurements (Sofia et al. 1994). More recently, using time resolved X-ray spectroscopy, absorption and emission features from Fe XVI-Fe XIX, Na X, Na XI, F VIII and F IX have been observed in the laboratory (Foord et al. 2001). Reliable dipole allowed and intercombination transition rates and oscillator strengths among the various fine-structure levels of the Mg-like ions are also essential for plasma diagnostics and several diagnostic line ratios for the solar transitions have been determined using transitions in such ions as Fe XV and Si III.

Excitation energies and line strengths for some low-lying excited states of ions in the Mg-like sequence were studied using multiconfiguration Dirac Fock (MCDF) wavefunctions by Cheng & Johnson (1977) and lifetimes were evaluated using multiconfiguration Hartree-Fock wave functions by Fischer & Godefroid (1982). Shorer et al. (1977) calculated, using the relativistic random phase approximation (RRPA) within the frozen core prescription, the excitation energies and absorption oscillator strengths for the $3\rm s^2~ {^1S_0} \rightarrow 3s3p~ ^1P_1$, $ \rm\rightarrow
3s4p~ ^1P_1, \rightarrow3s5p~ ^1P_1$ transitions of the Mg isoelectronic sequence. Fawcett (1983) applied the Cowan-Zelot Hartree-XR computer package (1967, 1968, 1976) to study the weighted oscillator strengths, energy levels and wavelengths of a number of ions in this sequence. Huang & Johnson (1985) used an improved version of the RRPA by introducing a multiconfiguration wavefunction as the reference state to determine excitation energies and oscillator strengths of the resonance transitions in the Mg isoelectronic sequence. Later, Chou et al. (1993) used a large-scale multiconfiguration RRPA (MCRRPA) including excitation channels from core electrons for the calculation of energy levels and oscillator strengths for the intercombination transition $\rm 3s^2~ {^1S_0} \rightarrow
3s3p~ {^3P_1}$ and the resonance transition $\rm 3s^2~ {^1S_0}
\rightarrow3s3p~ {^1P_1}$ in Mg-like ions and thus provided an ab initio treatment of the core-polarization (CP) effects. They showed that the excitation energies for both the intercombination and resonance transitions increased linearly with increasing nuclear charge. Nevertheless, their intercombination excitation energies with CP effects are not in good agreement with experiment and the percentage contributions due to CP effects are far too small compared to the corresponding effects in any CI calculation. They also found that the oscillator strengths for the intercombination transitions increased with increasing nuclear charge, whereas for the resonance transition the oscillator strengths decreased with increasing nuclear charge. Chou et al. (1993) made two important observations regarding oscillator strengths: (i) the differences between the length and the velocity results increased with increasing nuclear charge and (ii) the inclusion of CP effects brought the two forms into better agreement.

Butler et al. (1993) applied a modified R-matrix program to calculate energies and oscillator strengths for some 3l3l'transitions in a number of ions (Z=12-14, 16, 18, 20 and 26). However, this calculation includes no relativistic correction which is expected to contribute to ions of moderately high Z, such as calcium and iron. Relativistic multiconfiguration Dirac-Fock (RMCDF) and configuration interaction Dirac-Fock (CIDF), with and without core-valence correlations, calculations were performed by Stanek et al. (1996) for Mg I through Cl VI. They concluded that intravalence electron correlation is of crucial importance for the ions towards the neutral end of the sequence. However, they (Stanek et al. 1996) have studied only two transitions, namely $\rm 3s^2 (^1S_o)\rightarrow 3s3p (^1P_1)$ and $\rm 3s^2 (^1S_o)\rightarrow 3s3p (^3P_1)$. Valence and core-valence multiconfiguration Dirac-Fock and relativistic configuration interaction calculations including the Breit interaction were performed by Jönsson & Fischer (1997). However, they have studied only the $\rm 3s^2 (^1S)\rightarrow 3s3p (^{1,3}P)$ transitions in some low Z Mg-like ions. They found the Breit interaction to be of great importance both for the $\rm 3s3p$ $\rm {^3P}$ fine-structure splitting and for the intercombination rates. Chen & Cheng (1997) applied B-spline basis techniques to study transition energies in this sequence, and showed that by including the interaction of the valence electrons with this core, along with Breit, QED, and mass-polarization effects, fairly accurate transition energies could be obtained. Correlations with the core are especially important for determining accurate transition rates for the $\rm 3s^2 (^1S_0) \rightarrow 3s3p (^3P_1)$ intercombination transition.

In an attempt to resolve the controvertial energy position of the $\rm 3p^2(^1S_0)$ level, Deb & Msezane (1998) studied the fine- structure levels of Fe XV and confirmed the findings of Sugar & Corliss (1985) with or without the $\rm 3d^2$ configuration included in the calculation. Almaraz et al. (2000) have performed more extensive CIV3 calculations to study the configuration mixing among the levels involved in the $\rm 3s^2 (^1S_0) \rightarrow 3snp
(^{1,3}P_1)$ (n=3-5) transitions in Mg-like ions such as Si III, S V and Fe XV. They also used relativistic quantum defect orbital (RQDO) formalism and MCDF procedures in their investigations. The fine- tuned results for the resonance transition showed a 7% disagreement between length and velocity forms. Interestingly, with a much smaller basis set, level energies for Fe XV from our group (Deb & Msezane 1998; Deb et al. 1999) showed far better agreement with the measurement. However, their (Almaraz et al. 2000) calculated oscillator strengths using MCDF agree closely with their length form of the CIV3 results as do the recent experimental values. Their RQDO methodbased on a simple model potential gave reasonably good oscillator strengths for the ions towards the neutral end of the sequence, but poor agreement with the measured data and both their CIV3 and MCDF for the $\rm 3s^2 (^1S_0) \rightarrow 3s3p (^{1}P_1)$ transition in Fe XV. To obtain the excitation energies, transition probabilities, and lifetimes for the 3l-3l'electric dipole transitions in Mg-like ions (Z=13-100), Safronova et al. (2000) used a relativistic many-body perturbation theory (MBPT). The Breit interaction was also included in their calculation. Their calculated values for allowed transitions agree with experiment within the experimental uncertainties for the 3s3p $\rm (^1P_1)$ level but for many of the other excited levels their agreement with the measured values is rather poor. To extend the data needed for astrophysical applications, Aggarwal et al. (1999, 2000, 2001, 2003) have investigated energy levels, radiative rates, collision strengths and effective collision strengths for Fe XV and compared their results with other available data as well as assessed the reliability of other theoretical calculations.

More recently Biémont et al. (2002) applied the relativistic Hartree-Fock (HFR) and Multiconfiguration Dirac-Fock (MCDF) methods, including both intravalence and corevalence correlations, extensively to calculate the radiative lifetimes and oscillator strengths in Mg-like K VIII. In order to assess the accuracy of theoretical methods for transition rates, it is essential to have reliable experimental values for comparison. A number of comparatively accurate experimental values are available (Curtis 1991) both for the allowed transition $\rm 3s^2(^1S_0)$- $\rm 3s3p(^1P_1)$and for the intercombination transition $\rm 3s^2(^1S_0)$- $\rm 3s3p(^3P_1)$.

The purpose of this work is to study excitation energies, oscillator strengths, transition probabilities and lifetimes of allowed and intercombination transitions of two Mg-like ions, $\rm K^{7+}$ and $\rm Ti^{10+}$, accounting for the electron correlations through an extensive set of configurations and relativistic effects through the Breit-Pauli Hamiltonian.

2 Method of calculation

We use the computer code CIV3 of Hibbert (1975) and of Glass & Hibbert (1978) to optimize various orbitals and to calculate energy levels, oscillator strengths, radiative decay rates and lifetimes for the lowest 87 fine-structure levels of Mg-like K VIII and Ti Xl. Relativistic effects are included through the Breit-Pauli approximation via spin-orbit, spin-other orbit, spin-spin, Darwin and mass correction terms. The basic 14 configurations $\rm 3s^2$, 3s3p, $\rm 3p^2$, 3s3d, 3p3d, $\rm 3d^2$, 3s4s, 3s4p, 3s4d, 3s4f, 3p4s, 3p4p, 3p4d and 3p4f result in the 47 LS states: $\rm 3s^2 (^1S^e),
3s3p (^{1,3}P^o),$ $\rm 3p^2 (^1S^e, ^1D^e,^3P^e)$, $\rm 3s3d
(^{1,3}D^e)$, $\rm 3p3d (^{1,3}P^o$, $\rm ^{1,3}D^o$, $\rm ^{1,3}F^o)$, $\rm 3d^2
(^1S^e, ^1D^e$, $\rm ^1G^e$, $\rm ^3P^e$, $\rm ^3F^e)$, $\rm 3s4s (^{1,3}S^e)$, $\rm 3s4p
(^{1,3}P^o)$, $\rm 3s4d (^{1,3}D^e)$, $\rm 3s4f (^{1,3}F^o)$, $\rm 3p4s
(^{1,3}P^o)$, $\rm 3p4p (^{1,3}S^e$, $\rm ^{1,3}P^e$, $\rm ^{1,3}D^e)$, $\rm 3p4d
(^{1,3}P^o$, $\rm ^{1,3}D^o$, $\rm ^{1,3}F^o)$ and $\rm 3p4f (^{1,3}D^e$, $\rm ^{1,3}F^e$, $\rm ^{1,3}G^e)$. These LS states give rise to 87 fine-structure levels (47 even parity and 40 odd parity), corresponding to various J values.

Table 1: Parameters for the optimized orbitals.

We have taken the core configuration $\rm (1s^2~2s^2~2p^6)$ as frozen but the two valence electrons are distributed among all the possible combinations of the n=3 and n=4 orbitals. To be specific, we have included the $\rm 3s^2, 3p^2, 3d^2, 4s^2, 4p^2, 4d^2, 4f^2,$3s3d, 3s4s, 3s4d, 3p4p, 3p4f, 3d4s, 3d4d, 4s4d, and 4p4f (16 in total) configurations in the generation of the even states and 3s3p, 3s4p, 3s4f, 3p3d, 3p4s, 3p4d, 3d4p, 3d4f, 4s4p, 4s4f, 4p4d, and 4d4f (12 in total) configurations in the generation of the odd states. This set of configurations gives rise to 205 individual configurations, and therefore opening the core configuration and the inclusion of n=5 orbitals has been intentionally avoided to keep the computational requirements to a minimum. This will definitely affect the accuracy of our results (at least for some of the transitions) but leaves room for future improvement.

To construct the configuration-interaction (CI) wave functions we have used the Hartree-Fock radial orbitals for 1s, 2s, 2p and 3s given by Clementi & Roetti (1974) for K VIII and Ti XI. The CIV3 structure code of Hibbert (1975) has been used to optimize the 3p, 3d, 4s, 4p, 4d, and 4f orbitals on the $\rm 3s3p(^1P^o), 3s3d(^3D^e),
3s4s(^3S^e), 3s4p(^1P^o), 3s4d(^3D^e),$ $\rm 3s4f(^3F^o)$ levels, respectively. The explicit form of the wave function and the optimized orbitals are presented elsewhere (e.g. Deb et al. 1999). The parameters of the optimized orbitals for K VIII and Ti XI are given in Table 1.

In our calculation, apart from the inclusion of extensive CI (as described above), relativistic effects have also been accounted for through the Breit-Pauli approximation. For an ion like Mg, for which Z is not very high, this approach leads to satisfactory results. However, in order to improve the accuracy of our results further, we have adjusted the diagonal elements of the Hamiltonian matrix to reproduce the experimental energies as close as possible. This approach, justified further recently (McPeake & Hibbert 2000), leads to a better agreement with the experimental level energies as well as with the splitting between them. Oscillator strengths, transition probabilities and lifetimes are then calculated using these adjusted energies.

3 Result and discussions

New and improved calculations have been performed for the ions $\rm K^{7+}$ and $\rm Ti^{10+}$. Although we have calculated all the 87 fine-structure energy levels (ab initio and fine-tuned) for both the ions, only the levels belonging to the lowest five excited configurations (up to $\rm 3d^2$) are tabulated in Tables 2 and 3. The remaining fine-structure levels for both the ions are available; they can be obtained from the authors on request. For comparison we include the NIST (http://www.nist.gov) recommended data , the measurement of Churilov et al. (1989) and calculations of Fawcett (1983) and Fischer (2002). It is interesting to note that among the lowest 34 excited levels only the $\rm 3s3p(^1P_1)$ level shows the same value in the NIST data and the measurement of Churilov et al. (1989). For all other levels the two sets of results differ by from a few $\rm cm^{-1}$ to a few hundred  $\rm cm^{-1}$.

Table 2: Comparisons of fine-structure energy levels (in cm-1) for K VIII. CSS: measurement of Churilov et al. (1989); NIST: (http://www.nist.gov); CFF: calculation by Fischer (2002) and BF: calculation by Fawcett (1983). (a) and (b) in the present calculation represents fine-tuned and ab initio energies respectively.

During our fine tuning we have adjusted our ab initio energies mostly with the NIST data wherever they are available and with the data from Churilov et al. (1989) wherever the NIST data are not available. Our adjusted energies agree fairly well with both NIST data and the measurement of Churilov et al. (1989). For only a few levels does the calculation of Fawcett (1983) or of Fischer (2002) agree better with the measurement than does ours. For example, the result of Fawcett (1983) for the $\rm 3s3p(^1P_1)$ level merges exactly with both the NIST data and the measured value of Churilov et al. (1989) for both the ions. Similarly Fischer (2002) obtained the closest agreement with the measurement for the $\rm 3p^2 (^1D_2)$ level. The last column in Tables 2 and 3 represents the leading percentage of the configuration mixing. The first number in this column represents the percentage composition of the level on the same row. The next set of numbers of the form M(N) indicates that the next leading percentage is M% of the level number N under the "index" column. We note from the last column of these two tables that many of the levels are either pure or almost pure and some of them couple from weakly to strongly with other levels of the same J values. As an example in $\rm K^{7+}$, 74% of $\rm 3p^2 (^1D_2)$ couples with 24% of $\rm 3s3d(^1D_2)$ and vice versa. Similarly, in $\rm Ti^{10+}$, 62% of $\rm 3p3d(^3D_1)$ strongly couples with the $\rm ^3P_1$ level of the same configuration. It is to be noted here that in Table 2 level numbers 31 and 35 have some coupling with levels 44 and 63, respectively. In the list of 87 levels, level 44 represents $\rm 3s4d(^1D_2)$ and level 63 represents $\rm 3p4p(^1S_0)$. Similarly in Table 3 level 35 is coupled with level 36 which is $\rm 3s4s(^1S_0)$ in the list of 87 levels for this ion. We note that in the present calculation we have used the same set of configurations as used by Deb & Msezane (1999) for Fe XV and Das et al. (2003) for Ni XVII.

Table 3: Comparisons of fine-structure energy levels (in cm-1) for Ti XI. CSS; measurement of Churilov et al. (1989); NIST: (http://www.nist.gov); CFF: calculation by Fischer (2002) and BF: calculation by Fawcett (1983). (a) and (b) in the present calculation represents fine-tuned and ab initio energies respectively.

We present our oscillator strengths $f_{\rm L}$ (length form) and $f_{\rm V}$(velocity form) in Table 4 and also compare them with the corresponding results from other theoretical groups. Unfortunately, most of the other calculations considered only one spin allowed transition $\rm 3s^2 (^1S_o)\rightarrow 3s3p (^1P_1)$. There is generally satisfactory agreement between our results and those of others. Our results in the length and the velocity forms agree within 5% which is an indication of the high accuracy of our wave function generated for the calculation. For the ion K VIII we have compared our oscillator strengths mostly with two available sets of calculations due to Fawcett (1983) and Biémont et al. (2002). While the results of Fawcett (1983) agree very well with those of ours, the results of Biémont et al. (2002) are consistently lower by a few percent. This is interesting because Fawcett (1983) included a limited number of CI only among the n=3 complex, whereas Biémont et al. (2002) included a large number of configurations, promoting one of the valence electrons up to the n=8 complex, while restricting the other electron within the n=3 complex. However, this calculation (Biémont et al. 2002) ignored configurations such as $\rm 4f^2$, $\rm 4s4p$, $\rm 4s4d$, $\rm 4s4f$, $\rm 4p4d$, $\rm 4p4f$, $\rm 4d4f$ etc., which are found to be important in our calculation. In order to minimize the size of the table we have presented oscillator strengths for only a few optically allowed transitions. However, we have calculated oscillator strengths for all possible allowed and intercombination transitions among the lowest 87 fine structure levels of both the ions. Results for these transitions (over eleven hundred) are available from the authors on request.

Table 4: Oscillator strengths for some dipole-allowed transitions, $f_{\rm L}$ and $f_{\rm V}$ are the present length and velocity forms. OT represents other calculations.

In Table 4 we compare our results for the electric dipole transition $\rm 3s^2 (^1S_0) \rightarrow 3s3p (^{1}P_1)$ for both K VIII and Ti XI with five other calculations, where the results of Safronova et al. (2000) have the least agreement with those of ours. While the present calculation and that of Safronova et al. (2000) contain similar relativistic corrections through the Breit-Pauli Hamiltonian, our calculation involves a much larger set of interacting configurations. This could mean that the relatively large differences between the two sets of results are due to the use of a large number of interacting configurations in the present calculation.

Table 5: Lifetimes $(\rm 10^{-9}s)$ for some long-lived levels of K VIII and Ti XI. Superscripts a-e represent the same authors as in Table 4.

In Table 5 we present lifetimes (in nanoseconds) of some relatively longer lived levels of $\rm K^{7+}$ and $\rm Ti^{10+}$. Comparisons are made with other available calculations and the measurements of Curtis (1991) and Biémont et al. (2002). We note that our results give excellent agreement with the measurements. On the other hand the results of Safronova et al. (2000) slightly overestimate and those of Fischer (2002) slightly underestimate the measurement of Curtis (1991), for most of the transitions in Table 5. Finally in Tables 6 and 7 we present oscillator strengths, transition probabilities and line strengths for all allowed and intercombination transitions among the levels of K VIII and Ti XI within the n=3 complex only.

4 Summary and conclusion

The lowest 87 fine-structure levels of the Mg-like ions K VIII and Ti XI have been calculated. We have included large scale interacting configurations and relativistic effects such as spin-orbit, spin-spin, spin-other-orbit, mass correction and Darwin terms. All possible allowed and intercombination transitions have also been incorporated in the calculation of oscillator strengths, transition probabilities and lifetimes of these levels, but only the results of the lowest 35 levels of $\rm K^{7+}$ and $\rm Ti^{10+}$ within the n=3 complex have been presented in this paper. The rest of the 87 levels and associated data for all other allowed and intercombination transitions can be obtained from the authors on request. Present results for the fine-structure levels, oscillator strengths and lifetimes show good agreement with recent measurements. It is concluded that the interactions among the various excited configurations within the n=3 and n=4 complexes are of crucial importance to obtain accurate results.

Acknowledgements
N.C.D. is grateful for financial support from the Department of Science & Technology (DST), Government of India through a research grant SP/S2/L-12/99. A.Z.M. is supported by the U.S. Department of Energy, Division of Chemical Sciences, Office of Basic Energy Sciences, Offices of Energy Research.

References

  
Online Material

Table 6: Oscillator strengths ($f_{\rm L}$ and $f_{\rm V}$), transition probabilities ($A_{\rm L}$) and line strengths (LST) for allowed and intercombination transitions of K7+ among levels of n=3 complex. I and J are initial and final state of the transition as labeled in "Index'' column in Table 2.

Table 7: Oscillator strengths ($f_{\rm L}$ and $f_{\rm V}$), transition probabilities ($A_{\rm L}$) and line strengths (LST) for allowed and intercombination transitions in Ti10+ among its levels of n=3 complex. I and J are initial and final states of the transition as labeled under the "Index'' column in Table 3.



Copyright ESO 2004