Numerous papers for the transition dealt with here have appeared over the past 10 years or so. These have all been based on essentially the same R-matrix suite of programs as the one we use, but with different target wavefunctions. Obviously it is important to compare some of this earlier work with ours, especially in those cases where there exist significant differences between the results. Note that the comparisons we make concern calculations in which the target wavefunctions have more than two terms.
The three ions we have selected for comparison, namely those corresponding to Z = 14, 26 and 28, have been studied by Mohan and collaborators: Si+5 (Mohan & Le Dourneuf 1990), Fe+17 (Mohan et al. 1987), Ni+19 (Mohan et al. 1990). Each of these raises questions which we now comment on.
Mohan & Le Dourneuf (1990) tabulate the effective collision strength for Si VI as a function of temperature. They calculated using the collision strength they had obtained earlier by means of a 16 state close-coupling expansion in LS coupling. As can be seen in Fig. 11, where we plot as a function of log, their tabulated results rise dramatically with temperature when this exceeds about . Burgess et al. (1991) were able to make a reliable spline fit to Mohan & Le Dourneuf's data using the program OMEUPS (see Burgess & Tully 1992), the r.m.s. error of the fit being 2.2%. In order to get a fit to this accuracy Burgess et al. (1991) omitted the lowest two temperature points given in the table. Mohan (private communication) had previously informed one of us (JAT) that the published value at the lowest temperature, i.e. , was unreliable. However he gave no explanation why this was so.
Saraph & Tully (1994) state that at higher temperatures (not energies, as printed incorrectly in their paper) the results of Mohan & Le Dourneuf (1990) indicate a steep rise in the effective collision strength that is due to resonances converging on higher target states (i.e. on terms above ). The present investigation shows, however, that the reason for the increase manifested by the starred data points in Fig. 11 is principally due to the fact that Mohan & Le Dourneuf (1990) used the trapezoidal rule to carry out the integration over energy. Since there is a big gap in their tabulation of between about 10 Ry and 48 Ry, this numerical method should not have been used; see Burgess et al (1997b) for a detailed explanation. The gap is apparent in Fig. 12a where we plot as a function of energy: the collision strength is taken from the atomic databank at The Queen's University of Belfast (QUB). The gap is spanned in the figure by a dashed straight line. The resonances above 48 Ry in Fig. 12a, i.e. the highest target term included by Mohan & Le Dourneuf, are clearly responsible for the peak in near seen in Fig. 11 (curves a and b). The origin and physical significance, if any, of these resonances have not been explained by Mohan and Le Dourneuf. Fig. 12b shows that, as is to be expected, from the IRON Project has no structure in the region beyond about 12 Ry. Consequently the resulting curve c in Fig. 11 does not have a peak in the region of .
Figure 11: : Using
Mohan & Le Dourneuf's (1990)
from the QUB atomic databank and a) the trapezoidal rule,
b) the linear interpolation
method. : Mohan & Le Dourneuf (1990);
c) IRON Project
Figure 12: : a) by
Mohan & Le Dourneuf (1990) from
the QUB atomic databank; b) IRON Project
Saraph & Tully (1994) state that the effective collision strength given by Mohan et al. (1987) rises steeply at higher energies. This is wrong on two counts. Firstly, since the effective collision strength is a function of T they should have said "at higher temperatures'' and not "at higher energies''. Secondly, the effective collision strength that Mohan et al. (1987) tabulate does not vary enormously with temperature. As can be seen in Fig. 13, where we show versus log from Mohan et al. (1987) and the IRON Project, rises only slightly at temperatures above one million degrees. The effective collision strength from Mohan et al. (1987) lies about 10% lower than the present one. This may be due to the fact that they include less target terms in their close-coupling expansion and use a much coarser energy mesh than we do.
Figure 13: : a) Mohan et al. (1987);
b) IRON Project
In Figs. 14 and 15 we compare the collision strength by Mohan et al. (1987) from the QUB atomic databank with the present calculation.
Figure 14: : a) by
Mohan et al. (1987) from the QUB
atomic databank; b) IRON Project
Figure 15: : a)
by Mohan et al. (1987) from the QUB
atomic databank; b) IRON Project
Mohan et al. (1990) tabulate for this ion as a function of T. The range covered is from T = 104 to T = 107. Now according to Arnaud & Rothenflug (1985) Ni XX has its maximum coronal abundance at about , so the tabulation by Mohan et al. does not cover the temperature range of interest fully. Figure 16 shows as a function of log. The starred points () correspond to data taken from Table 1 in Mohan et al. (1990), while the full line curve represents the present results. Both sets of results have a maximum at a temperature below log. Although the maxima occur at approximately the same value of log, they differ in size somewhat. At temperatures above 106, the present exhibits a second peak which, although not present in the tabulation of Mohan et al. (1990), does show when we thermally average their collision strength from the QUB atomic databank. We conclude that Mohan et al. (1990) may have truncated the integration over energy too soon when calculating .
Figure 16: : a)
Mohan et al. (1990); b) IRON Project