A&A 416, 375-381 (2004)
DOI: 10.1051/0004-6361:20031703
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
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
,
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
and the resonance transition
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
and
.
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
transitions in
some low Z Mg-like ions. They found the Breit interaction to be of
great importance both for the
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
intercombination
transition.
In an attempt to resolve the controvertial energy position of the
level, Deb & Msezane (1998) studied the fine-
structure levels of Fe XV and confirmed the findings of Sugar &
Corliss (1985) with or without the
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
(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
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
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
-
and for the intercombination transition
-
.
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,
and
,
accounting for
the electron correlations
through an extensive set of configurations and relativistic
effects through the Breit-Pauli Hamiltonian.
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 ,
3s3p,
,
3s3d, 3p3d,
,
3s4s, 3s4p, 3s4d, 3s4f, 3p4s,
3p4p, 3p4d and 3p4f result in the 47 LS states:
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
and
,
,
.
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
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
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
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.
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
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
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
,
74% of
couples with 24% of
and
vice versa. Similarly, in
,
62% of
strongly
couples with the
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
and level 63 represents
.
Similarly in Table 3
level 35 is coupled with level 36 which is
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
(length form) and
(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
.
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
,
,
,
,
,
,
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,
and
are the present length and velocity
forms. OT represents other calculations.
In Table 4 we compare our results for the electric dipole transition
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
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
and
.
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.
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
and
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.
Table 6:
Oscillator strengths (
and
), transition probabilities (
)
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 (
and
), transition probabilities (
)
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.