next previous
Up: On the interpolation

2. Fitting isoelectronic data

The reduced charge number tex2html_wrap_inline1147 is defined by
equation223
where Z is the atomic number (i.e. nuclear charge in atomic units) for which the range is tex2html_wrap_inline1151. tex2html_wrap_inline1153 is usually equal to the number of bound electrons N although on occasions it is more convenient if tex2html_wrap_inline1157 as in the second example (Fig. 2 (click here)) where tex2html_wrap_inline1159. Here we omit the neutral atom case and consequently obtain a much better fit to the positive ion data points. The parameter C is adjusted in order to optimise the spline fit. The plots we use are either of type 2 (tex2html_wrap_inline1163 for large Z) or type 3 (tex2html_wrap_inline1167 for large Z). The different types of plot are defined and discussed by Burgess & Tully (1992).

  figure236
Figure 2: Fluorine sequence: tex2html_wrap_inline1171, tex2html_wrap_inline1173, C = 1.8

2.1. Type 2 plots

Figure 1 (click here) concerns temperatures of maximum abundance tex2html_wrap_inline1179 for fluorine-like ions in conditions of coronal ionization equilibrium. From Arnaud & Rothenflug's (1985) tabulation we obtain estimates of log(tex2html_wrap_inline1181) for 10 ions in the sequence. Our spline fit can be used to complete their tabulation for the intermediate ions which they did not consider. Here we input tex2html_wrap_inline1183 and treat it as a type 2 case. The optimised fit has C = 4.4 and rms error 0.38%. Since we only wish to interpolate the data but not extrapolate them, we are not concerned with the high Z limit point. In fact it does not exist since there is a weak (logarithmic) divergence in this case.

Figure 2 (click here) deals with the ionization energy I of fluorine-like ions. Since I increases like tex2html_wrap_inline1193 as Z tends to tex2html_wrap_inline1197, we treat this as type 2 by inputting tex2html_wrap_inline1199, with I in tex2html_wrap_inline1203. The limit point at tex2html_wrap_inline1205 corresponding to tex2html_wrap_inline1207 is the hydrogenic value tex2html_wrap_inline1209, while the optimised fit (rms error 0.16%) is obtained with C = 1.8.

  table254
Table 1: Spline fitting parameters for the curves shown in the figures

  figure263
Figure 3: Aluminium sequence: tex2html_wrap_inline1235, tex2html_wrap_inline1237, C=7.8

2.2. Type 3 plots

Figure 3 (click here) shows our way of interpolating the ground term magnetic dipole transition energies for the aluminium sequence. It is instructive to compare this way of plotting the data with that used by Edlén (1942) in his Fig. 1 (click here). We invert the observed spin-orbit splitting tex2html_wrap_inline1243 of the tex2html_wrap_inline1245 term, which Edlén (1942) gives in tex2html_wrap_inline1247 in Table 2 (click here), and take the square root. This yields a quantity tex2html_wrap_inline1249 at high Z. The value of the limit point for this type 3 plot is

  table281
Table 2: (a) Log(tex2html_wrap_inline1253 from the spline fit shown in Fig. 1. tex2html_wrap_inline1255 is the temperature of maximum coronal abundance for F-like ions. (b) Ionization energy in tex2html_wrap_inline1257 for F-like ions deduced from the spline fit shown in Fig. 2. The results given here are calculated using data from Table 1

  table312
Table 3: Transition energy in tex2html_wrap_inline1299 for the fine structure splitting in the ground configuration of aluminium and Al-like ions. (a) Using the spline interpolation data given in Table 1. (b) Using Edlén's (1964) extrapolation formula

  table336
Table 4: Fine structure collision strengths tex2html_wrap_inline1331 at zero energy for the ground tex2html_wrap_inline1333 term in C-like ions. (a) tex2html_wrap_inline1335; (b) tex2html_wrap_inline1337; (c) tex2html_wrap_inline1339. Results obtained using the interpolating spline fit given in Table 1


equation367
and is deduced from the well-known expression for the spin-orbit splitting given by Edlén (1964, p. 167), viz.
equation375


4 where the screening parameter s' > 0 is a finite quantity which varies slowly with Z. The optimised spline fit (rms error of 0.05%) is obtained when C = 7.8. Edlén's (1964) extrapolation formula for Z-s' as a function of Z is given by the pair of equations
displaymath1509

displaymath1511
which reduces to finding the roots of the cubic equation (Maple 1995)
equation386
where
equation388
and
equation390
Only one of the roots of (5) is a physical solution, and its range of validity as a function of Z is limited, since for tex2html_wrap_inline1515 the screening parameter s' becomes negative.

Figure 4 (click here) shows the compacting and interpolation of Blaha's (1969) distorted wave results for three fine structure collision strengths at threshold energy in the carbon sequence. The transitions are between the lowest three levels. Blaha gives results for neutral carbon and 7 ions in the sequence. Here we take tex2html_wrap_inline1519 and input tex2html_wrap_inline1521 which decreases like tex2html_wrap_inline1523 as tex2html_wrap_inline1525. The high Z limits for tex2html_wrap_inline1529 are from Saraph et al. (1969).

2.3. Spline fitting parameters

The five knot values tex2html_wrap_inline1533 and C parameter for each of the spline curves shown in Figs. 1 to 4 are given in Table 1 (click here). By means of the program FUNCTION SPLINE tex2html_wrap_inline1537, (see Burgess & Tully's 1992 Appendix), it is possible to interpolate the reduced data points shown in the figures. Notice that tex2html_wrap_inline1539 and by using the definition of tex2html_wrap_inline1541 in terms of tex2html_wrap_inline1543 and C, and knowing the type of plot one can obtain A(Z).

The interpolated results given in Tables 2 (click here), 3 (click here) and 4 (click here) are calculated using data in Table 1 (click here) and the spline program from Burgess & Tully (1992).

  figure411
Figure 4: Carbon sequence: tex2html_wrap_inline1549, tex2html_wrap_inline1551, tex2html_wrap_inline1553, tex2html_wrap_inline1555, tex2html_wrap_inline1557


next previous
Up: On the interpolation

Copyright EDP Sciences
web@edpsciences.com