A&A 406, 741-749 (2003)
DOI: 10.1051/0004-6361:20030660

Trends in E2 and M1 transition rates between 3p3/2 and 3p 1/2 levels in 3s2 3pk systems

E. Charro - S. López-Ferrero - I. Martín

Departamento de Química Física, Universidad de Valladolid, 47011 Valladolid, Spain

Received 23 December 2002 / Accepted 2 April 2003

Abstract
The analysis of forbidden lines, such as E2 and M1, in the atomic spectra emitted by certain ions is important for the study of the plasma in astrophysical objects and fusion devices. Transition rates between  $\rm 3p_{3/2}$ and  $\rm 3p_{1/2}$ levels in  $3{\rm s}^2~ 3{\rm p}^k$ ions that belong to the aluminium sequence have been calculated with the Relativistic Quantum Defect Orbital (RQDO) method. The present data are tested by comparison with other theoretical data available in the literature, and the regularity of the transition intensities along the isoelectronic sequence for both (E2 and M1) lines, as well as their competition, are also analysed.

Key words: atomic data - methods: numerical

1 Introduction

Within the last few years, the determination of atomic oscillator strengths and related properties for ionized systems has been an active research area on account of its considerable interest in astrophysics, plasma physics, laser physics and thermonuclear fusion research. In astrophysics, transitions from metastable levels of highly ionized atoms are observed in emission in low-density media. It is well known that line intensities and ratios provide excellent diagnostics for temperature and density and consequently for elemental abundances (Feldman 1981; Herter et al. 1982). In the laboratory, the  $n{\rm p}^k$ configurations of medium-Z elements have attracted considerable attention in recent years as spectroscopic diagnostics in Tokamak plasmas. This is due to the fact that they may give rise to strong lines at much longer wavelengths ( $\lambda = 2000~$Å) than other transitions in the same ion. However, it is still difficult to measure accurately oscillator strengths, particularly for highly ionizad atoms where the available experimental techniques have limitations, and the reproduction in the laboratory of the very extreme conditions to which the constituents of the universe are subjected is extremely difficult.

Experimental transition rates for forbidden lines are scarce and, in most cases, missing when the degree of ionization is higher than 2. It is therefore necessary to rely upon theoretical values. One way to assess the quality of a given theoretical method is a comparison of the so-obtained results with some other measurable quantity which might be available (mostly the observed energies), and by comparison with the results supplied by other computers codes.

On the theoretical side, quantum-mechanical calculations have, thus, become crucial in the accomplishment of these goals. Wiese (1987) long pointed out the usefulness of procedures which account for relativistic effects, the mass production of data for complex atomic systems with increasingly refined semiempirical approaches; and the application of regularities among transition rates, such as the dependence on the nuclear charge, for the scaling of transition rates along isoelectronic sequeces or spectral series.

Many calculations have been done for the ions isoelectronic with Al, most of them corresponding to allowed transitions. In this context, we have reported oscillator strengths corresponding to dipole-allowed transitions for a group of Al-like ions ( Z = 13 - 103) in a previous paper (Lavín et al. l997), in which the Relativistic Quantum Defect Orbital (RQDO) method, was employed. For forbidden transitions, no data have been reported with this method for this sequence before the present calculations.

Certain forbbiden transitions only become possible through the fine-structure splitting of an LS term into a series of J levels. Examples are the M1 and E2 transitions between the $\rm 3s^2 3p$ $ {}^2{\rm P}_{1/2}$ and $\rm 3s^2 3p$ $\rm {}^2P_{3/2}$ levels in the aluminium sequence. In the non-relativistic approximation, these levels are degenerate. Calculations for the 3p $\rm {}^2P_{3/2} \to 3p $ $\rm {}^2P_{1/2}$ transition, corresponding to both E2 and M1 forbidden mechanisms, have presently been performed with the RQDO method, for the aluminium-like elements ranging from Sc IX (Z = 21) to Hg LXVIII (Z = 80). It is well known that fine-structure lines are interesting from a spectroscopic point of view, given their usefulness in spectral analysis in astrophysics and fusion plasma research.

The systematic trend of the A-values for E2 and M1 transitions along the isoelectronic sequence has also been analysed in this work. In previous works (e.g., Lavín et al. l997) we studied the systematic trends of E1 spectral lines. The regular behaviour of oscillator strengths along isoelectronic sequences has proven to be a useful tool for analysing a large body of f-value data. It may also be exploited to obtain additional oscillator strengths by simple interpolation techniques (Martin & Wiese 1996).

Some theoretical E2 and M1 transition rates for aluminium-like elements are available in the literature. The most extensive group of data has been published by Huang (1986), who calculated transition rates for the ions comprised in the $13 \leq Z \leq 106$range with the multiconfiguration Dirac-Fock (MCDF) technique.

In the present work, we deal with heavy ions and in high degree of ionization for which the relativistic contributions to the wavefunctions and energies of the levels involved in the transitions may not be negligible. We have consequently followed a theoretical method that gives explicit account of relativistic effects, the Relativistic Quantum Defect Orbital (RQDO) formalism, as mentioned above.

For over a decade, we have applied the RQDO method (Martín & Karwowski 1991; Karwowski & Martín 1991; Martín et al. 2000) to the calculation of atomic data in several isoelectronic sequences (see, e.g. Martín et al. 1993; Lavín et al. 1997; Charro & Martín 1998; Charro et al. 1996, 1997, 2000, 2001). The RQDO formalism, as opposed to sophisticated and costly self-consistent-field procedures, is a simple but reliable analytical method based on a model Hamiltonian. It has the great advantage that the computational effort is not increased for Z-high systems. The convenience of employing exactly solvable model potentials for calculating atomic transition rates manifests itself not only from a practical point of view but also because of the involved physical implications, when they are capable of achieving a good balance between computational effort and accuracy of results. In the present calculations, a general good accord between the RQDO transition probabilities and the comparative theoretical data is also found.

2 The RQDO method

The Relativistic Quantum Defect Orbital (RQDO) method has been described in detail in previous papers (Karwowski & Martín 1991; Martín et al. 2000). Therefore, we shall only briefly summarise its most fundamental aspects.

The relativistic formulation of the Quantum Defect Orbital (QDO) formalism, as developed by Martín & Karwowski (1991) and Karwowski & Martín (1991) was based on the decoupling of the Dirac second-order equation, and in the interpretation of the resulting solutions that had been carried out by Karwowski & Kobus (1985, 1986). These authors generalised the decoupling procedure of Biedenharn (1962) in such a way that it may be applied to the case of a general atomic spherical potential. They also showed that both the large and small components of the Dirac wavefunction may be recovered without any computational effort.

The quantum defect orbitals are solutions of the Schr $\rm\ddot{o}$dinger equation (Simons 1974; Martín & Simons 1975, 1976)

\begin{displaymath}\left[- \frac{{\rm d}^2}{{\rm d}r^2}~ + \frac{\chi (\chi + 1 ...
...~ \right] {\psi_k}^{\rm QD} = 2 E^{\rm QD} {\psi_k}^{\rm QD},
\end{displaymath} (1)

where $Z_{\rm net}$ is the nuclear charge on the active electron at large r and

\begin{displaymath}\chi = n - \delta + c,
\end{displaymath} (2)

where $\delta$ is the quantum defect and c is an integer chosen to ensure the correct number of nodes and the normalization of the radial wave function. The eigenvalue  $ E^{\rm QD} $ in Eq. (1) depends only upon the nonintegral part of $\chi $ and, hence, is independent of c. The quantum defect is empirically obatained from the following equation:

\begin{displaymath}\ E^{\rm QD} = E^x = - \frac {(Z_{\rm net})^2 } { 2 ( n - {\delta})^2 },
\end{displaymath} (3)

where Ex is the experimental energy.

Analogously to the nonrelativistic case, Eq. (1), the RQDO equation is written as

\begin{displaymath}\left[- \frac{{\rm d}^2}{{\rm d}r^2}~ + \frac{\Lambda(\Lambda...
...~ \right] {\psi_k}^{\rm RD} = 2 e^{\rm RD} {\psi_k}^{\rm RD},
\end{displaymath} (4)

with

\begin{displaymath}\ Z'_{\rm net} = Z_{\rm net} \left( 1 + {\alpha}^2 E^x\right ),
\end{displaymath} (5)

and,

\begin{displaymath}\Lambda = s - 1 - \delta' + c,
\end{displaymath} (6)

when j = l + 1/2, or

\begin{displaymath}\Lambda = - s - \delta' + c,
\end{displaymath} (7)

when j = l - 1/2. $ \delta' $ is the relativistic quantum defect, c is an integer with the same physical significance as that of Eq. (2); $ Z'_{\rm net}$ is the scaled nuclear charge acting on the valence electrons at large radial distances; Ex is also here the experimentally measured energy. Atomic units are used throughout.

The relativistic quantum defect is determined empirically. We have

\begin{displaymath}\ e^{\rm RD} = - \frac {(Z'_{\rm net})^2 } { 2 ( \eta - {\delta'})^2 }
\end{displaymath} (8)

and is has been proven (Martín & Karwowski 1991) that the value of $ \delta' $ may be obtained from:

\begin{displaymath}- \frac {(Z_{\rm net})^2 } { 2 ( \eta - {\delta'})^2 } =
E^...
... {( 1 + {\alpha}^2 E^x / 2 )}{( 1 + {\alpha}^2 E^x )^2 }\cdot
\end{displaymath} (9)

Here $\eta$ is the relativistic principal quantum number, related with n, the principal quantum number, as follows

\begin{displaymath}\eta = n - \mid k\mid + \mid s \mid .
\end{displaymath} (10)

It should be stressed that this formulation is "exact'' in the sense that is equivalent to a four-component formulation based on the standard first-order form of the Dirac equation. All matrix elements, in particular the transition moments, may be expressed in a simple way, using the solutions of the second-order equation. A set of recurrent formulas which are fulfilled by the radial integrals (Martín et al. 2000) allows the formalism to be very simple and compact. Karwowski & Martín (1991) have observed that the relativistic density distribution approximates very well the exact one at large values of r. At small distances, the quality of the density deteriorates, as happens when the nonrelativistic QDO densities are compared with the exact nonrelativistic ones (Simons 1974). Fortunately, the consequences of this drawback are very seldom reflected in the quality of the QDO and RQDO transition rates, because the strongest contribution to radial matrix elements comes, in most cases, from large radial distances.

The most important difference between the RQDO and QDO equations is the explicit dependence of the former on the total angular momentum quantum number k. As a consequence, values of the relativistic quantum defect are determined from the fine structure energies rather than from their centers of gravity. The corresponding relativistic quantum defect orbitals are different for each component of a multiplet and, if c = 0, they retain the nodal structure of the large components of the hydrogenic Dirac wave function.

The Relativistic Quantum Defect Orbital is expressed in terms of Kummer's functions,

$\displaystyle {\psi_k (r)}^{\rm RD} = C~ \exp \left( \frac {-Z'_{\rm net} r}{\e...
...\eta^* - \Lambda - 1 ),
2 \Lambda + 2 , \frac {2 Z'_{\rm net} r}{\eta^*}\right)$     (11)

and C is the normalization constant,

\begin{displaymath}C~ = \left(\frac {2Z'_{\rm net} }{\eta^* }\right)^{1/2} \left...
... - \Lambda )}\right]^{1/2} \frac {1}{ \Gamma (2 \Lambda + 2)}
\end{displaymath} (12)

with $\eta^*$ = $\eta$ - $ \delta' $.

Since the effective Hamiltonian in Eq. (4) includes a screening term, the relativistic quantum defect orbitals are approximately valid in the core region of space. Core polarization effects are implicitly included in the calculations, as they are accounted for in the ${\Lambda}$ parameter of the model Hamiltonian. On the other hand, given the one-electron nature of the RQDO formalism, it can be expected to perform better in highly excited states, where the active electron interacts less with the core and other valence electrons, than in low-lying energy states. The relativistic quantum defect orbitals lead to closed-form analytical expresions for the transition integrals. This allows one to calculate transition probabilities and oscillator strengths by simple algebra and with high computational costs/efficiency ratio.

3 Electric dipole forbidden transitions

The RQDO methodology supplies one-electron radial wavefunctions, characterized by the n, l and j quantum numbers (Martín & Karwowski 1991), that we employ in the transition matrix elements for the initial and final states of the active electron. These correspond to levels of a given L, S and J symmetry in many-electron atoms.

3.1 Transitions via electric quadrupole mechanism

Let us consider electric multipole radiation. The operator responsible for the interaction of this radiation with an atomic system is the electric multipole moment introduced by the general expression

\begin{displaymath}\ Q_m ^{(q)} = - \sum r_i ^q
C_m ^{(q)}(i)
\end{displaymath} (13)

with

\begin{displaymath}\ C_m ^{(q)}(i) = \left( \frac {4\pi } { 2 q + 1 } \right)^{1/2} Y_{qm}(\theta_i,\varphi_i)
\end{displaymath} (14)

where $Y_{qm}(\theta_i,\varphi_i)$ are the spherical harmonics.

In the particular case of q = 2, the transitions take place via the electric quadrupole mechanism, E2. In this context, the electric quadrupole line strength for a transition between two states within the LSJ-coupling, which is the coupling scheme followed throughout in this work, in the notation of Condon & Shortley (1935), is given by the equation

\begin{displaymath}S^{(2)}_{nlj,n'l'j'} ~ = ~ \frac {2}{3} ~
\left\vert\left\la...
...rallel Q^{(2)} \parallel \alpha' J' \right\rangle\right\vert^2
\end{displaymath} (15)

where $\alpha $ and $\alpha'$ are here the sets of quantum numbers which define the lower and upper states, respectively. The matrix-elements have the form
$\displaystyle \langle \alpha J \parallel Q^{(2)}\parallel \alpha' J' \rangle
=\...
...J')
\langle \alpha L \parallel Q^{(2)}\parallel
\alpha' L' \rangle \delta_{SS'}$     (16)

and the Kronecker delta, $\delta_{SS'}$, occurs because the electric multipole operators do not depend on spin. Thus, we define a line factor  $R_{\rm line}$ by
$\displaystyle R_{\rm line}(SLJ,S'L'J') = (2J+1)^{1/2} ~ (2J'+1)^{1/2}~
W(SJL' 2, LJ')$     (17)

where W(SJL' 2, LJ') are the Racah coefficients which can be described in terms of 6j-symbols

\begin{displaymath}W(SJL' 2, LJ') = (-1)^{S+J+L'+2}
\left\{
\begin{array}{ccc}
S & J & L \\
2 & L' & J'
\end{array}\right\}.
\end{displaymath} (18)

The line strength then becomes
                         S(2)nlj,n'l'j' = $\displaystyle \frac {2}{3} ~ (2J+1) ~ (2J'+1)~
\left[W\left(SJL' 2, LJ'\right)\right]^2 ~$  
    $\displaystyle \times \mid \left\langle \alpha L \parallel Q^{(2)}\parallel \alpha '
L' \right \rangle \mid ^2$  
  = $\displaystyle ~ \frac {2}{3} ~ \left[R_{\rm line} ~ \left(SLJ, S' L' J'\right)\right]^2$  
    $\displaystyle \times \mid
\left\langle \alpha L \parallel Q^{(2)}\parallel \alpha 'L' \right\rangle \mid ^2 .$ (19)

As mentioned above, the relevant operator is independent of spin, and hence spin does not change during the transition. The following selection rules apply to E2 transitions: $\Delta S = 0$, $\Delta L$ and $\Delta J = \pm 2 , \pm 1 , 0$, $L + L' \geq 2$, and $J + J' \geq 2$.

Thus, the reduced matrix element $ \langle \alpha L \parallel Q^{(2)}\parallel \alpha '
L' \rangle$ provides the relative strengths of different multiplets. We can write this reduced matrix element as the product of a single-electron reduced matrix element $ \langle n l \parallel Q^{(2)}\parallel n '
l' \rangle$, which depends on only the quantum numbers of the jumping electron, and a factor that we shall call the multiplet factor, $R_{\rm mult}$:

\begin{displaymath}\langle \alpha L \parallel Q^{(2)}\parallel \alpha 'L' \rangl...
... ~ \times \langle n l \parallel Q^{(2)}\parallel n 'l' \rangle
\end{displaymath} (20)

where

\begin{displaymath}\langle n l \parallel Q^{(2)}\parallel n ' l' \rangle ~ = \la...
... \rangle ~ \langle R_{nlj}\vert Q^{(2)}\vert R_{n'l'j'}\rangle
\end{displaymath} (21)

in Eq. (21), $\langle R_{nlj}\vert Q^{(2)}\vert R_{n'l'j'}\rangle$ is the radial transition integral, $ \langle l \parallel C^{(2)} \parallel l' \rangle $ is the pertinent reduced matrix element, which can be evaluated using 3j-symbols by the following expresion

\begin{displaymath}\langle l \parallel C^{(2)} \parallel l' \rangle = [l,l']
\le...
...\begin{array}{ccc}
l & 2 & l' \\
0 & 0 & 0
\end{array}\right)
\end{displaymath} (22)

and $R_{\rm mult}$ may be expressed as follows:
$\displaystyle R_{\rm mult}(\alpha L, \alpha' L' ) ~ = ~ (2L+1)^{1/2} ~ (2L'+1)^{1/2}~
W(L_{\rm c}Ll' 2, l L')$     (23)

where $L_{\rm c}$ refers to the orbital angular momentum of the atomic core, and the last symbol in Eq. (23) is,

\begin{displaymath}W(L_{\rm c}Ll'2,lL') = (-1)^{L_{\rm c}+L+l'+2}
\left\{
\begin...
...y}{ccc}
L_{\rm c} & L & l \\
2 & l' & L'
\end{array}\right\}.
\end{displaymath} (24)

Finally, the line strength takes the form
$\displaystyle S^{(2)}_{nlj,n'l'j'} = \frac {2}{3} ~
\mid R_{\rm line}(SLJ, S'L'...
...allel L' \rangle ~
\langle R_{nlj}\vert Q^{(2)}\vert R_{n'l'j'}\rangle \mid ^2.$     (25)

The total strength of the line strengths of all the lines in a multiplet is therefore:

\begin{displaymath}S({\gamma} L, {\gamma}' L') =
{\sum_{JJ'}} ~ S({\gamma} J, {\gamma}' J').
\end{displaymath} (26)

Thus, the line strength for a multiplet transition is
$\displaystyle S({\gamma} L, {\gamma}' L') =
(2S + 1) ~ \mid \langle \alpha L \p...
... L, \alpha' L') ~ \langle n l \parallel Q^{(2)}\parallel n 'l' \rangle \mid ^2.$     (27)

The relationships between the line strength S(2) (in atomic units, $e^2 a_{\rm o}^4$), the oscillator strength f(2) (dimensionless), and the transition rate A(2) (in s-1) for E2 transitions are given by (Sobelman 1979)

\begin{displaymath}g'A^{(2)} = (8\pi ^2 \hbar \alpha / m \lambda ^2) gf ^{(2)} =
(6.6703 \times 10^{15}/ \lambda ^2) gf^{(2)}
\end{displaymath} (28)


\begin{displaymath}g'A ^{(2)} = (32 \pi ^5 \alpha c a_{\rm o}^4 / 15 {\lambda}^5) S^{(2)}
= (1.11995\times 10^{18}/ \lambda ^5) S^{(2)}
\end{displaymath} (29)

where $\lambda $ is the transition wavelength (in Å), g and g' are the degeneracies of the lower and upper state, respectively, and $\alpha $ is again the fine-structure constant.

3.2 Transitions via magnetic dipole mechanism

Let us consider now the magnetic multipole radiation. The operator responsible for the interacton of this radiation with matter is the magnetic multipole moment introduced by the general expression

\begin{displaymath}\ {\cal M}_m^{(q)} = \sum_i r_i^{q - 1} \sqrt{q (2 q -1 )}
\end{displaymath}


$\displaystyle \left \{ \frac {1}{q + 1}
\left [C_m ^{(q -1)} \times l_i^1 \right ]^{q} +
\left[ C_m ^{(q -1)}(i) \times s_i^1 \right ]_m^{q} \right \}\cdot$     (30)

In the particular case of q = 1, when the transitions are taking place via the magnetic dipole mechanism, M1,

\begin{displaymath}\ {\cal M}_m^{(1)} = \sum_i \frac {1}{2} \left[ l_i^1 + 2 s_i^1 \right]
\end{displaymath} (31)

where l1 and s1 are vector operators of orbital and spin momenta.

In this context, the formulae for the probability of magnetic dipole (M1) transitions do not contain radial integrals, an thus, the line strength, $S^{\rm M1}$, is independent of the frequency of the transition.

The magnetic dipole line strength for a transition between two states within the LSJ-coupling, in the notation of Condon & Shortley (1935), is given by the equation

\begin{displaymath}S^{(\rm M1)}_{nlj,n'l'j'} ~ = ~
\vert\langle \alpha J \parallel {\cal M}^{(1)} \parallel \alpha' J' \rangle\vert^2.
\end{displaymath} (32)

The final expression is
$\displaystyle S^{\rm M1}(LJ,LJ-1) = \frac {1}{ 4J } \left(L+\frac{1}{2}+J+1\rig...
...}-J+1\right)
\left(\frac{1}{2}+J-L
\right) \left(J+L- \frac {1}{2} \right)\cdot$     (33)

The transition rate $A^{\rm M1}$ (in  $\rm s^{-1}$) is related with the line strength  $S^{\rm M1}$, in Bohr magneton units  ${\mu_{\rm B}}^2$, by the expression

\begin{displaymath}A^{M1} = (2.697~ 10^{13}/ \lambda ^3)\frac {S^{\rm M1}}{g'}\cdot
\end{displaymath} (34)

The relationship between the oscillator strength  $f^{\rm M1}$ (dimensionless) and the transition rate  $A^{\rm M1}$ (in  $\rm s^{-1}$) is given by Eq. (28), which is generally valid for E1, E2 and M1 transitions.

4 Regularities along the isoelectronic sequence

Regularities in individual oscillator strengths along an isoelectronic sequence as functions of the nuclear charge have been predicted for E1 transitions from conventional perturbation theory (Cohen & Dalgarno 1966; Dalgarno & Parkinson 1967). We shall now extend this procedure to the analysis of electric quadrupole transitions.

If we denote the generalised radial transition integral

\begin{displaymath}\vert\langle R_{nlj}\vert r^q \vert R_{n'l'j'}\rangle\vert = I_{nlj,n'l'j'}^{(q )},
\end{displaymath} (35)

for a given operator, the variation of a matrix element as a function of Z may be studied by the Rayleigh-Schr$\ddot o$dinger perturbation theory. If we introduce $ \rho = Z.r $ and $\epsilon = E Z ^{-2}$ in a.u., the expansions of $\psi$ and E in powers of the perturbation parameter, Z-1, are:

\begin{displaymath}\psi_{nlj,n'l'j'} = \psi_0 + \psi_1.Z^{-1}+ \psi_2.Z^{-2}+ ...
\end{displaymath} (36)

and

\begin{displaymath}E = Z^2(\epsilon_0 + \epsilon _1 .Z^{-1}+ \epsilon _2 .Z^{-2}+ ...).
\end{displaymath} (37)

The electric dipole (q = 1) transition integral is given by the following Z-expansion:

\begin{displaymath}\vert\langle \psi_{nlj}\vert\rho \vert\psi_{n'l'j'}\rangle\ve...
..._{nlj,n'l'j'}^{(1)} = I_0^{(1)}.Z^{-1} +
I_1^{(1)}.Z^{-2}+ ...
\end{displaymath} (38)

where I0 is the corresponding integral for hydrogen and the superscript (1) refers to E1 or electric dipole transitions.

The dipole line strength, or squared radial integral Eq. (38), may be written as:

Snlj,n'l'j'(1) = S0(1).Z-2 + S1(1).Z-3+ S2(1).Z-4+ ... (39)

and the expression for the E1 oscillator strength will be

fnlj,n'l'j'(1) = f0(1) + f1(1).Z-1+ f2(1).Z-2+ ... (40)

For the line and oscillator strengths, as well as for the transition rate A, it is possible to perform paralell nuclear charge expansion representations in the quadrupole case to study systematic trends of E2 S- or f-values along an isoelectronic sequence. Here the transition operator, $\rho^2$, leads to (Biémont & Godefroid 1978),

Snlj,n'l'j'(2) = S0(2).Z-4 + S1(2).Z-5+ S2(2).Z-6+ ... (41)

with

\begin{displaymath}S_0^{(2)}= \vert< \psi_{nlj}\vert\rho^2 \vert\psi_{n'l'j'}>\vert^{2} .
\end{displaymath} (42)

For the transition rate A, the corresponding expansion is the following,

An'l'j',nlj(2) = A0(2).Z6 + A1(2).Z5 + A2(2).Z4 + A3(2).Z3+ ..., (43)

and the quadrupole oscillator strength may be written, as,

fnlj,n'l'j'(2) = f0(2).Z2 + f1(2).Z + f2(2) + f3(2).Z-1+ ..., (44)

where f0(2), f1(2), f2(2),... are proportional to some power of $\Delta \epsilon _0$, the hydrogenic transition energy corresponding to the transition under study. We may also consider of interest to analyze the behaviour of f(2).Z-2 along the isoelectronic sequence. From Eq. (44) it is inferred:

fnlj,n'l'j'(2) . Z-2 = f0(2) + f1(2).Z-1+ f2(2).Z-2+ ... (45)

When no change of the principal quantum number occurs ( $\triangle n = 0 $) during the transition, if we ignore relativistic effects, $\triangle \epsilon _0 = 0 $. We have thus:

flj,l'j'(2) = f3(2).Z-1+ f4(2).Z-2+ ... (46)

That is, the curve of f(2) versus Z-1 would tend to an asymptotic zero value in the high-Z side of the sequence if no relativistic effects were present. Equation (45) is expected to be a good approximation to the behaviour of the oscillator strength with Z-1 for E2 transitions at least for the first few ions of the sequence of a light element, e.g., Na I, where relativistic effects still are not very important. However, as soon as these effects set in, deviations from Eq. (45), as well as from all the expressions in this section, are expected to occur.

5 Energy data

5.1 Energy levels

In the RQDO context, the quantum defects are often extracted from energy level data. As energy levels, we have employed those provided from the critical compilation by Kelly (1987). These data comprise the elements from Al I to Kr XXIV. Given that the transitions we can study are constrained by the availability of observed energy levels, for the ions for which no data were found, interpolated (e.g. Br XXIII) and extrapolated ( $37 \leq Z \leq 80$) energy values have been used. In order to do this, we first fitted the experimental energy data (belonging to ions between Z = 13 and Z= 36, with the exception of the value for Z = 35) to a polynomic dependence of the nuclear charge, Z. The expression obtained for the 3p $\rm {}^2P_{3/2}$ level is,

$\displaystyle E_{3{\rm p}_{3/2}}=0.136551 ~Z^4 ~- 4.32196 ~Z^3 ~+ 72.6444 ~Z^2
-779.776 Z~+ 3577.09$     (47)

where E is in cm-1 and the correlation factor was found to be 0.9999.

As we advance in the sequence, a Z-value will be reached for which the spin-orbit coupling scheme can no longer be expected be pure LSJ scheme but an intermediate one between that and jj. We have not changed, however, the formulation of either the radial functions or the angular factors employed in the calculation of oscillator strengths at any point in the sequence. We expect that the use of emppirical energy data will, at least to some extent, account for intermediate coupling, as it does, at all values of Z, for configuration interaction.

5.2 Ionization energy

Other input data also needed in the RQDO calculations are the ionization energies of the atomic systems. The corresponding values from Al I to Kr XXIV have been taken from Kelly (1987). Analogously to the energy levels, a fitting of the experimental data (with the exception of Br XXIII) has led to the following expression for the ionization energy in terms of powers of Z

I.E. = -2.07672   - 369447   Z + 17639.8   Z2 - 70.2101   Z3 (48)

where I.E. is given in cm-1, an the correlation factor was equal to 0.9999. With this formula, the gap of values for I.E. can be interpolated (Z = 35) and extrapolated ( $37 \leq Z \leq 80$).

6 Analysis of the results

In Table 1, the wavelengths corresponding to the 3p  $\rm {}^2P_{3/2} \to 3p $ $\rm {}^2P_{1/2}$ spectral line for a large number of Al-like ions are collected. The data employed in the present calculations are compared with those reported by Huang (1986). In Table 2, the presently calculated oscillator strengths are compared with those reported by Huang (1986), which are, to our knowledge, the only data existing in the literature, as regards these transitions and atomic systems. No experimental results have been found.

Table 1: Wavelenght ($\lambda $ in $\AA$) for the $\rm 3p_{1/2}-3p_{3/2}$ spectral line in the aluminium sequence.

The values supplied by Huang (1986), have been calculated with the multiconfiguration Dirac-Fock (MCDF) computer code of Desclaux (1975), which includes contributions from the Breit interaction, as well as those from the Lamb shifts. Huang included all the configurations that arise from within the n = 3 complex in his calculations.

The E2 and M1 f-values collected in Table 2, comprise the ions from Sc IX (Z = 21) up to Hg LXVIII (Z = 80). The magnetic dipole the line strength has been taken as 1.3333  ${\mu_{\rm B}}^2$according to Eq. (33).

Table 2: E2 and M1 oscillator strengths for the $\rm 3p_{1/2}-3p_{3/2}$ spectral line.


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{ALFIG1.eps}
\end{figure} Figure 1: Systematic trends of the transition rates for the 3p $\rm {}^2P_{3/2}$-3p $\rm {}^2P_{1/2}$ spectral line via the electric quadrupole (E2) mechanism along the aluminium sequence.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{ALFIG2.eps}
\end{figure} Figure 2: Systematic trends of the transition rates for the 3p $\rm {}^2P_{3/2}$-3p $\rm {}^2P_{1/2}$ spectral line via the magnetic dipole (M1) mechanism along the aluminium sequence.
Open with DEXTER

The main contribution to the intensity of the spectral line can be analysed by comparing the f-values corresponding to the M1 and E2 transitions for each ion of the sequence. The contribution of the E2 transition turns out to be negligible with respect to that of the M1 transition, the latter being several orders of magnitude larger than the former. An inspection of Table 2 reveals a rather good general agreeement between our results and those reported by Huang (1986).

Additionally, in Figs. 1 and 2, the behaviour of the E2 and M1 transition probabilities, respectively, along the isoelectronic sequence is analysed in graph form. The value of Log A has been plotted against the atomic number, Z, from Z = 21 up to Z = 80. The MCDF (Huang 1986) data have also been included. These figures serve two purposes. One is to show the agreement or deviations between our results and the MCDF calcuations. The second is to reflect the systematic trends complied with by the individual RQDO A-values along the isoelectronic sequence, which have long been considered as a qualitative proof of correctness (Wiese & Kelleher 1998) and can be used for the interpolation or extrapolation of non-calculated data.

In Fig. 3, we show only the present RQDO data in order to show the competition between both mechanism (E2 and M1) for the same transition.

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{ALFIG3.eps}
\end{figure} Figure 3: Competition between the electric quadrupole (E2) and the magnetic dipole (M1) mechanism in the 3p $\rm {}^2P_{3/2}$-3p $\rm {}^2P_{1/2}$ transition along the aluminium sequence.
Open with DEXTER

6.1 Regularities in homologous series

Regularities in individual oscillator strengths along an isoelectronic sequence as functions of the nuclear charge have long been predicted from conventional perturbation theory (Cohen & Dalgarno 1966; Dalgarno & Parkinson 1967), and their usefulness in analysing the f-values has been repeatedly stressed (see, e.g., Martin & Wiese 1996 and references therein).

The oscillator strengths or transition rates for strong analogous ( $nl \to n'l'$) transitions of homologous atoms are expected to be similar on account of the analogous outer electron structure of such elements. Thus, all alkali elements, for example, would be expected to exhibit similar f-values for the leading lines in a Rydberg series, and it is, indeed, well known that the principal resonance lines for the alkalis have f-values close to unity. However, one has to consider that as the elements within a chemical family become heavier, the outer atomic structure becomes modified due to the presence of unfilled electron shells. Generally, the systematic behaviour for homologous atoms is not expected to be as closely adhered to as for the case of a given transition along an isoelectronic sequence, where the electron configuration remains the same and only a scaling of the nuclear charge occurs, but it is observed in some cases. As an illustration, Figs. 4 and 5 display the behaviour of the transition rate for the $n{\rm p}$  ${}^2{\rm P}_{3/2} \to n{\rm p}$ $\rm {}^2P_{1/2}$ spectral line for two sequences, those of aluminium (n=2) and Boron (n=3) against $\xi $. $\xi $ is defined as the charge of the atomic core ( $\xi = Z - N$), where Z is the nuclear charge and N is the number of core electrons. In this way, following the spectroscopic notation, for a neutral atom we have $\xi = 1$. The values for the transition rates via the magnetic dipole (M1) mechanism are shown in Fig. 4, and for the transition via the quadrupole (E2) mechanism in Fig. 5. The RQDO A-values plotted for B-like ions have been taken from Charro et al. (2001).

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{aa3440f4.eps}
\end{figure} Figure 4: Systematic trends for the fine-structure 3p $\rm {}^2P_{3/2}$-3p $\rm {}^2P_{1/2}$ spectral line of the aluminium and boron sequences (n = 3 and 2, respectively). For both sequences, our RQDO transition rates via the magnetic dipole (M1) mechanism are plotted versus $\xi $, the charge on the atomic core or effective nuclear charge.
Open with DEXTER

7 Conclusions

The RQDO procedure has once more proved to be a very useful tool for estimating transition rates and oscillator strengths. In the particular case of the Al-like ions, a general, satisfying, agreement is found between the RQDO results and those obtained with the rather more elaborate MCDF procedure which accounts for a great deal of electron correlation (Huang 1986).

Finally, the RQDO A-values for the E2 and M1 transitions follow the expected systematic trends along the isoelectronic sequence, which have long been considered as a qualitative proof of correctness and can be exploited for the interpolation or extrapolation of non-calculated data. We are confident in that the present RQDO f-values, may be potentially useful in astrophysics and fusion plasma research.


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{ALFIG5.eps}
\end{figure} Figure 5: Systematic trends for the fine-structure 3p $\rm {}^2P_{3/2}$-3p $\rm {}^2P_{1/2}$ spectral line of the aluminium and boron sequences (n = 3 and 2, respectively). For both sequences, our RQDO transition rates via the electric quadrupole (E2) mechanism are plotted versus $\xi $, the charge on the atomic core or effective nuclear charge.
Open with DEXTER

Acknowledgements

This work has been supported by the DGI of the Ministry of Science and Technology within Project No. BQU2001-2935-CO2-O2 and the J.C.L. within Project No. VA119/02.

References



Copyright ESO 2003