Dirac Rmatrix calculation for electronimpact excitation of S xiii
^{1}
Key Laboratory of Optical Astronomy, National Astronomical Observatories,
Chinese Academy of Sciences, 100012
Beijing, PR China
email: [gyliang;gzhao]@bao.ac.cn
^{2}
Pakistan Atomic Energy Commission, PO Box 1114, 44000
Islamabad,
Pakistan
Received:
2
May
2013
Accepted:
28
May
2013
Context. Sulfur emission lines in the soft Xray and extremeultraviolet regions are observed in a variety of laboratory and astrophysical spectra. But accurate electron impact excitation data for S xiii for stateoftheart NLTE spectral models are scarce.
Aims. We calculated electronimpact excitation collision strengths and effective collision strengths of S xiii for transitions among the lowestlying 98 finestructure states 1s^{2}2lnl′ corresponding to principal quantum numbers n = 2,3,4. The effective collision strengths for these transitions were computed over a wide temperature range (log _{10}T_{e} (K) = 4.53–7.53) for various astrophysical plasma conditions.
Methods. We used the fullyrelativistic parallel Dirac Rmatrix code to calculate collision strengths. To generate target wavefunctions and energy levels for scattering calculations, we employed the GRASP0 multiconfiguration DiracFock code for states up to n = 5.
Results. The wavefunctions are generated from 27 configurations – 1s^{2}2lnl′(n = 2,3,4,5) – giving rise to 166 jj energy levels. The collision and effective collision strengths among the lowest 98 finestructure levels are compared with the previous theoretical calculations. The collision strengths for most transitions agree well at higher incident electron energies.
Conclusions. The resonant contributions to effective collision strengths are most dominant at lower temperatures.
Key words: atomic data / atomic processes / plasmas / line: identification
© ESO, 2013
1. Introduction
Sulfur is an abundant element that widely exists in various astrophysical environments (Feldman 1992). In solar and stellar coronae, rich emission lines emanating from sulfur ions have been detected over a wide wavelength range from Xray to extremeultraviolet (EUV) bands (Acton et al. 1985; Raassen et al. 2002; Li et al. 2013). Belike ions are of great interest for modelling hightemperature astrophysical and laboratory plasmas, because this isoelectronic sequence provides line ratios suitable for temperature diagnostics (Landi et al. 2001; Kavanagh et al. 2010). The emission lines emanating from S xiii have been used as a test bed to infer electron densities and plasma temperatures (Keenan & Berrington 1985). The emission lines from S xiii are mainly formed for temperatures in the range of 1.3–16 MK, with a peak formation temperature of 2.5 MK (Bryans et al. 2009). Keenan et al. (1988) theoretically investigated the emission line strength of S xiii and compared them with XUV solar flares observations. Moreover, spectral emission lines of S xiii have also been studied extensively in laboratory plasmas. Goldsmith et al. (1974) studied the spectra of S xiii and S xiv and identified twenty new transitions in the 25–40 Å wavelength range. Khardi et al. (1994) studied the spectra of highly ionised Beisoelectronic ions (Al, Si, S, and Ar) with a beamfoil experiment in the 50–400 Å wavelength range. Lepson et al. (2005), with the aid of the Lawrence Livermore electronbeam ion traps EBITI and EBITII, observed a large number of sulfur emission lines in the soft Xray and EUV regions. For astrophysical spectra, the authors detected ten strong emission features from S xiii in the 20–75 Å wavelength range.
Previous calculations about atomic structure data for S xiii have been reported by several authors. Muehlethaler & Nussbaumer (1976) calculated the transition probabilities for several ions of the Beisoelectronic sequence within an n = 2 complex intermediatecoupling scheme. Tully et al. (1990) reported energy levels, oscillator strengths, and photoionisation crosssections for opacity calculations of fifteen members of the Beisoelectronic sequence in the Russell Saunders (LS) coupling scheme. Recently, Mania et al. (2009) calculated the excitation energies, lifetimes, and oscillator strengths of levels belonging to the configurations of 1s^{2}2lnl′ (n = 2,3,4,5) of S xiii using the multiconfiguration HartreeFock method, and carried out an extensive comparison with the experimental results. On the experimental side, oscillator strengths, transition rates, and lifetimes within the n = 2 complex of S xiii have been studied in a number of experimental works (Pegg et al. 1977; Träbert & Heckmann 1980; Yang et al. 2005; Du et al. 2005; Bhattacharya et al. 1998; Schippers et al. 2012).
Accurate electronimpact excitation data for S xiii are scarse in the literature. The very few previous electronimpact excitation calculations on S xiii have been carried out using nonperturbative CoulombBorn and distorted wave (DW) methods. An early calculation of collision strengths for 17 Belike ions (10 ≤ Z ≤ 74) was performed by Sampson et al. (1984) using the CoulombBornexchange method between levels of n = 2,3 configurations. Their collision strengths were calculated at nine electron impact energies (1 ≤ ϵ ≤ 15 in threshold units) in the intermediatecoupling scheme. (Zhang & Sampson 1992, hereafter ZS92) provided the relativistic DW collision strengths and oscillator strengths for 85 Belike ions (8 ≤ Z ≤ 92) for all of the possible 45 △ n = 0 transitions among 10 levels corresponding to 1s^{2}2l^{2} configurations. This method does not include contributions from resonances. A partial relativistic CoulombBethe (PRCBe) approximation was used to estimate the contributions from large partial waves, which may not be accurate for high incident electron energies (Zhang & Fontes 2013). Until recently, the most extensive DW collision strength calculation for S xiii was performed by (Landi & Bhatia 2008, hereafter LB08) using the University College London SUPERSTRUCTURE (SS) code (Eissner et al. 1974). The collision strengths were calculated for lowest eighteen configurations (2l^{2}, 2lnl′, n = 3,4,5 with l = s,p and l′ = s,p,d) in the CoulombBethe approximation. The authors reported the collision strengths at seven (10, 20, 45, 90, 135, 180, and 225 Ry) and five (45, 90, 135, 180, and 225 Ry) incident electron energies for the transitions within the three lowest n = 2 configurations and for the transitions between the lowest five levels as well as the levels belonging to the n = 3,4,5 configurations. This data set is employed in astrophysical spectral modelling packages, such as the current version of the Chianti database (Landi et al. 2012). Resonant contributions to the collision strength in this collisional data set were omitted, which are pertinent for lowtemperature plasmas but their omission is not justified in any accurate electron collisional data. However, higher levels arising from the configurations 2s4f and 2p4f are clearly not included in their calculations, while it is necessary to include them for accurate practical spectral modelling techniques. In addition, the above collision strength calculations focused on a limited number of transitions between levels of n = 2 configurations and the lowest five levels, corresponding to n = 3,4,5 configurations.
There are no previous Rmatrix calculations for electronimpact excitation data on S xiii apart from the work of (Keenan 1988, hereafter K88), to the best of our knowledge. In their Rmatrix calculation, electronimpact excitation rates for S xiii were calculated by interpolating effective collision strength of three Belike ions (Ne, Si, Ca) for all transitions amongst 1s^{2}2l^{2}, and results were presented in the form of polynomial fits to the effective collision strengths. Albeit less accurate, it is the only collisional data set on S xiii pertaining to resonance contributions available in the literature to date. In the present work, we primarily intend to perform fully relativistic collision strength calculations for S xiii up to n = 4 by employing the parallel Dirac Rmatrix code (DARC) by taking into account the effect of resonances, since they can lead to a considerable enhancement of the collision rates. All the 98 jjlevels corresponding to the configurations 1s^{2}2l^{2} and 1s^{2}2lnl (n = 3,4) are incorporated to delineate the target structure properly. Then Maxwellianaveraged effective collision strengths are evaluated for electron temperatures log _{10}T_{e} (K) = 4.53–7.53.
2. Computational methods
To accomplish collision strength calculations with DARC, we initially need to calculate atomic structure parameters (electronic orbitals, energy levels, and energy states) to delineate the S xiii target ion. In our calculation, we employed the GRASP0 (generalpurpose relativistic atomic structure package) code (Grant et al. 1980), which is a fully relativistic jjcoupling code based on the multiconfiguration DiracFock method. The relativistic corrections arising from the Breit interaction and quantum electrodynamics (QED) effects (vacuum polarisation and Lamb shift) have also been included in GRASP0. For the atomic structure calculations, we optimised the energy levels by using the extended averagelevel (EAL) mode, in which a weighted (proportional to 2j + 1) trace of the Hamiltonian matrix is minimised. The Lagrange multipliers between orbitals of the same generalised occupation number are set to zero to achieve convergence in a selfconsistent field (SCF) process. The wave functions generated by GRASP0 were used as input data for further scattering calculations. The DARC is based on the jjcoupling scheme and uses the DiracCoulomb Hamiltonian in the Rmatrix approach. Relativistic effects were considered in a systematic way in both the target description and the scattering model (colliding electron+target ions). Inclusion of finestructure in the channel coupling increases the size of the corresponding Dirac Hamiltonian matrix substantially. Therefore, it is common practice to perform largescale calculations in the atomic structure part and then reduce the size of the target state wavefunctions for the scattering part to solve the ionpluselectron (or (N + 1)electron) colliding system.
2.1. Energy levels
The wave functions were generated for the 27 lowestlying configurations 1s^{2}2lnl′ (n = 2,3,4,5), which give rise to 166 finestructure levels. The lowest 98 levels in Rydbergs from states up to n = 4 displayed in Table 1 were employed to construct the target in the scattering calculations. In this table, we also list the experimental results obtained from the Chianti v7^{1} atomic database and the NIST v5^{2} atomic spectra database to estimate the accuracy of wavefunctions used in the representation of target states. LB08 calculated electronimpact collision strengths, energy levels, oscillator strengths, and transition rates of the 92 finestructure levels belonging to the configurations (1s^{2}2l^{2},1s^{2}2lnl′ (n = 3,4,5,l′ = s,p,d)) by using the SS code. Their calculations lack 2s4f and 2p4f configurations, compared with our atomic structure model. Our calculations use the same atomic model as Mania et al. (2009) who employed the multiconfigurational HartreeFock (MCHF) method.
Level energies (Ryd) of S xiii from different calculations, along with the compilations of experimental values and other theoretical values.
A comparison of our calculated level energies with those of NIST and other theoretical results is shown in Fig. 1. To provide an immediate overview of the differences between theory and experiment, we plot the ratio of the level energies between the available observed values together with other theoretical results and our calculations. One can see from Fig. 1 and Table 1 that the GRASP0 ten lowestlying levels (2l^{2}) show a considerable discrepancy with NIST values, which is still below 2.5% however. We differ by 2.2% at maximum for the 2s2p^{1}P_{1} level from NIST. However, our theoretical model yields more accurate energies of higher excited states, and our results agree within less than 0.2% for higher configurations with those of NIST. This good agreement between our calculations and the observed values along with the MCDHF/MCHF calculations from the NIST database is consistent for all higherlying states (see Fig. 1). It is obvious from Fig. 1 that the level energies of our model and the DW energies of LB08 agree within 0.2%. Furthermore, we employed more states originating from two additional configurations (2s4f, 2p4f), than LB08. The addition of these correlation functions improve the accuracy of the target model for the scattering part. Overall, the energy differences seen between our calculations and those of LB08 with the observed energies are very similar. These two sets of calculations agree excellently with each other within 0.2%, as depicted in Fig. 1. The energies of the MCHF calculations of Mania et al. (2009) show a better agreement with the NIST observed values than the LB08 calculations because of energy corrections via a leastsquares procedure. Such adjustments by Mania et al. (2009) were made to optimise the electronic parameters of the oscillator strengths and lifetimes of S xiii. Their resulting level energies agree excellently within 0.2% with the observed energies found from the NIST database, as can be seen in Table 1. The agreement between our predicted energies and those of the MCHF calculations by Mania et al. (2009) is also lower than 2.1% for the lowestlying ten states, while the energy differences for the higherlying levels are lower than 0.2%. It should be noted that our atomic model is equal in size to that calculated by Mania et al. (2009). Because of the similarity of these two models, we have listed only a few energy levels in Table 1 from Mania et al. (2009). Table 1 clearly shows that our calculated energies and those of other theoretical calculations agree very well. Based on this analysis, we conclude that there is no (major) discrepancy (better than 2.5%) between theory and the experimental values for the energy levels of S xiii. This confirms the reliability of our target used for scattering calculations.
Fig. 1 Ratio between our level energies and previous calculations: black triangle – DW calculations of Landi & Bhatia (2008), blue circle – experimental values, red cross – both the NIST’s MCDHF compilations and MCHF results from Mania et al. (2009). 

Open with DEXTER 
The energyorder of the levels of all three independent calculations (SS/MCHF/GRASP0) is same with a few exceptions. Among the available observed energies, our GRASP0 energyorder differs only in a few cases, namely for levels (45–46, 65–67). Similarly, the level labelling in the SS calculations of LB08 varies for some levels in our GRASP0 and the MCHF calculations of Mania et al. (2009). In Table 1, we present our energy levels by closely following the energy ordering of LB08 to facilitate comparing our results. We expect that energy ordering is accurate due to the small energy differences specifically for strong transitions. Therefore, the overall magnitude of collision strengths is probably not be affected because the level ordering is correct with a few exceptions. It is also helpful to compare the oscillator strengths and Avalues to better appreciate the quality of wavefunctions used in the target description.
2.2. Oscillator strengths
For a given transition i → j, the absorption oscillator strength (f_{ij}) and emission transition probabilities A_{ji} (in s^{1}) are related by the following expression: (1)where m and e are the electron mass and charge, respectively, c is the velocity of light, λ_{ji} is the transition energy/wavelength in Å and ω_{i} and ω_{j} are the statistical weight factors of the initial (i) and final (j) states, respectively.
We computed electricdipole (E1), magneticdipole (M1), electricquadrupole (E2), and magneticquadrupole (M2) transitions amongst the 98 finestructure levels using the GRASP0 code. In Table 2, a sample of the weighted oscillator strengths (gfvalues) and Avalues of the transitions among the lowest ten levels are compared with the results from the LB08 and MCDHF compilations from the NIST, as well as available experimental gfvalues. Our Avalues agree excellently with the results of LB08 and MCDHF compilations.
In Table 2, all transitions obtained from the three sets of calculations agree within 25%, while the agreement for the E1 transitions is better than 15%. Our gfvalues agree with the experimental values within 10% expect for a single transition, namely that of 4–8 (2s2p–2p^{2} [entity ! #x20 !] ^{3}P_{2}). Our calculated gfvalue for this transition shows a maximum discrepancy of 27% with the corresponding experimental value, but is consistent with other theoretical calculations within 4%. In Fig. 2, we also compare the gfvalue with the L08 and MCDHF collections from NIST. Only the transitions of the levels up to the n = 3 configurations are provided in the MCDHF database. Figure 2 illustrates that the oscillator strength of about 89% of the transitions agree with the L08 calculations within 20%. Our calculation also agrees with the values from the MCDHF calculations within 20% for most of the transitions. To conclude, we may state that for almost all strong transitions (gf ≥ 0.01), Our radiative rates are accurate to better than 20%. However, for some weak transitions (gf ≤ 0.01), our calculation differs significantly (by more than a factor of 2) from the results of LB08. Overall, we expect that these atomic structure parameters are reliable because of the small differences with observations. The uncertainties in the collision strength data due to some inaccuracies in the target structure data are correspondingly small.
2.3. Lifetimes
The lifetime τ of an excited level j can be defined as follows: (2)Lifetimes of excited states can be measured experimentally, and several efforts in the literature have been devoted to measure the lifetimes of excited levels corresponding to the n = 2 complex of S xiii. To access the accuracy of our calculations, we compared our calculated lifetimes with available experimental values in Table 3. We also derived lifetimes of some levels of S xiii based on the Avalues of LB08. The lifetimes calculated by Mania adjusted with experimental values are closer to the experimental values than our lifetime predictions with GRASP0. Our lifetimes agree well with the recent experimental measurements (Du et al. 2005) within the measurement uncertainty except for two excited levels (2s2p , 2p^{2} [entity ! #x20 !] ^{3}P_{0}), which are lower by 6.25% and 7.4%, respectively. Our calculated lifetime of the 2p^{2} [entity ! #x20 !] ^{3}P_{0} level is 15% lower than the measured value of Bhattacharya et al. (1998). As for the lifetime for the level 2s2p [entity ! #x20 !] ^{1}P_{1}, our calculation is 23% lower than the measured value of Pegg et al. (1977). Overall, the discrepancy between our calculated lifetimes and all measured ones is within 10% except for lifetimes of the above two levels. The measured lifetimes agree with each other within the experimental uncertainties. Moreover, theoretical lifetimes are closer to the experimental measurements of Du et al. (2005), which has the smallest experimental uncertainties among the available measured lifetimes. Based on this comparison, we may state that the accuracy of our calculated lifetimes is better than 10%.
Comparison of weighted oscillator strengths and transition rates (Avalues) with other theoretical and/or experimental results.
Fig. 2 Ratio between our gfvalues and other theoretical calculations: open circle – DW calculations of Landi & Bhatia (2008), solid triangle – NIST’s MCDHF collections. 

Open with DEXTER 
2.4. Scattering calculations
The calculations of electronimpact excitation collision strength (Ω) were performed using a set of parallel DARC Rmatrix code^{3} (Ballance & Griffin 2006), originally developed by Norrington (Norrington & Grant 1987). The Rmatrix boundary radius is about 7.36 au, which is determined automatically by the code to enclose most of the target wavefunctions. We employed 30 continuum orbitals for each angular momentum in the expansion of wavefunctions. The allowed energy range for Ω can be up to 160 Ryd, which is about three times the ionisation potential of S xiii. Contributions from partial waves up to J = 41.5 are included. This is sufficient for most transitions in the energy range ≲160 Ryd. However, a large number of the Jvalue would be preferred for some allowed transitions above the resonance region. To account for the contribution of higher partial waves, we included a “topup”, based on the CoulombBethe method (Burgess & Sheorey 1974) in our calculation. In the resonance region, we used a fine energy mesh of 8.317 × 10^{6} z^{2} Ryd and 10^{4} energy points (up to 47.9 Ryd), where z is the effective or residual charge (12) of the ion. This is sufficient to resolve the resonance structures in the closed channel region. Above 47.9 Ryd and up to three times the ionisation potential of this ion, a coarser energy mesh (6.655 × 10^{4} z^{2} Ryd) and 2000 energy points were employed.
Comparison of the lifetimes with available experimental results.
In astrophysical and laboratory plasma applications, often the excitation rate coefficients (in units of cm^{3} s^{1}) are used. In this context, the excitation rate coefficient for a transition i → j can be expressed as follows: (3)where T_{e} is the electron temperature in K, ΔE_{ij} is the transition energy, and k is the Boltzmann constant, the Υ_{ij} is the the dimensionless thermally averaged effective collision strength, which can be defined as follows: (4)where E_{j} is the incident energy with respect to the level j. The effective collision strengths were calculated in a wide temperature range (3.38 × 10^{4}–3.38 × 10^{7} K) to encompass all temperatures of importance in astrophysical and plasma applications. The highenergy Born/Bethe limits were calculated using the AUTOSTRUCTURE (AS) code (Badnell 1986). These limits were then used to extrapolate the collision strength at the energy above 160 Ryd to obtain the effective collision strength at high temperatures.
3. Results and discussions
Collision strengths and the Maxwellianaveraged effective collision strengths among 98 finestructure levels of S xiii for a total of 4753 transitions were calculated adopting a suite of parallel Dirac Rmatrix programs. Therefore, it is not feasible to tabulate the collision strengths and the effective collision strengths because of the large number of transitions and the incident electron energies involved. A complete set of collision strengths and effective collision strengths for all transitions involved can be obtained online. In this section, we solely focus on the collision strength and the corresponding effective collision strength of some typical transitions. Figures 3–7 display the collision strengths and/or the associated effective collision strengths for the selected transitions. Comparisons of these electronimpact collisional parameters are mainly made from the earlier Rmatrix work of K88 and the two independent DW calculations by ZS92 and LB08. The selection of transitions shown in Figs. 3–7 are based on the kinds of agreement and the extent of differences we observe with previous theoretical calculations. Here, we also present Table 4 to compare a sample of the current effective collision strengths with those from the Rmatrix calculations of K88 and the DW calculation of LB08. Only the values at three temperatures (log _{10}T_{e} (K) = 6.23,6.53,6.93) are displayed in Table 4.
Fig. 3 Collision strength for the 1–5 (2s^{2} [entity ! #x20 !] ^{1}S_{0}–2s2p [entity ! #x20 !] ^{1}P_{1}) transition: solid black curve – our calculation without “topup”, dashed red curve – our calculation with “topup”, solid circle – DW calculations of Landi & Bhatia (2008). 

Open with DEXTER 
Fig. 4 Collision strength 1–9 and 1–10: solid black curve – our DARC results, solid blue circle – DW calculation of Zhang & Sampson (1992), open red circle – DW calculation of Landi & Bhatia (2008). The insets in the graphs compare the corresponding effective collision strength (Υ) as a function of the electron temperature: open black circle – our DARC results, open star – Rmatrix interpolated work of Keenan (1988), open upward triangle – DW calculations from Landi & Bhatia (2008). 

Open with DEXTER 
Fig. 5 Collision strengths for the optically allowed transitions 1–13 and 2–7, optically forbidden transition 1–4, and the intercombination transition 1–11: solid curve – our DARC calculation, open red square – DW calculations from Landi & Bhatia (2008), solid blue circle – calculations from Zhang & Sampson (1992). The insets show the details of the collision strengths near the background. 

Open with DEXTER 
Fig. 6 Comparison of our DARC effective collision strengths: open circle – our calculation, open red upward triangle – DW calculation of Landi & Bhatia (2008), open blue star – interpolated Rmatrix values of Keenan (1988). 

Open with DEXTER 
In Fig. 3, we compare the collision strength for the optically allowed transition 1–5 (2s^{2}^{1}S_{0}–2s2p ^{1}P_{1}) with the available DW collision strengths of LB08, who reported the collision strengths of S xiii at seven and five incident electron energies, as described above. In Fig. 3, our results with “topup” procedure are shown with the dashed red curve, while the solid black curve displays our results without the “topup” method. The background values of the collision strength calculated with the “topup” procedure agree well with LB08, but the collision strength without the “topup” is lower than LB08 by 27%, especially above the threshold excitation. From Fig. 3, we can see that if we neglect this “topup” procedure to our collision strength calculations, it will lead to inaccuracies in the calculated results for higher partial waves, which ultimately will affect the calculations of the effective collision strengths. Furthermore, we see from the inset in Fig. 3 that the “topup” contributions near the resonance threshold to the collision strengths of the dipoleallowed transition are negligible and the DW calculations agree well with the background level in this lowenergy resonance region. We also compared all available transitions reported in LB08 and ZS92. Our results for the background of the collision strength agree well with LB08 and ZS92 for the vast majority of transitions with a few exceptions. For example, LB08 noticed that their collision strengths for the transition 1–9 (2s^{2}^{1}S_{0}–2p^{2}^{1}D_{2}) and 1–10 (2s^{2}^{1}S_{0}–2p^{2}^{1}S_{0}) are higher by a factor of 1.5–2.0 than the results of ZS92. For the transition 1–9, the background level of our collision strength is slightly lower than LB08, but higher than ZS92. This difference is caused by the configuration interaction effects, which have been explored in LB08.
Fig. 7 Comparison of our effective collision strengths with those of Landi & Bhatia (2008): dashed lines indicate agreement within 20%, dotted lines indicate agreement within a factor of 2. 

Open with DEXTER 
It is obvious from Fig. 4 that the background level of our collision strength for the transition 1–10 is higher by a factor of 1.5–2.0 than those of LB08 in the threshold region, but converges in the highenergy region. The effect of differing collisions between our results and the LB08 calculations on the resulting effective collision strengths are shown in the insets in Fig. 4. These two insets clearly show the impact of resonances on the effective collision strength for both transitions. The available effective collision strength data for the two transitions from the interpolated work of K88 and the DW calculations of LB08 over a range of temperature exhibit strong deviations from our results. There is even little or no agreement for all three independent calculations for the 1–9 transition with the exception of a single temperature value of LB08 in the highertemperature domain, where the background level of our collision strength and that of LB08 agree well. For the lowtemperature region, the height of the lowlying resonances enhances the effective collision strength significantly. For the 1−10 transition, the effective collision strengths produced by K88 are approximately 3%–78% lower, while the LB08 predictions 61%–78% are lower than our DARC values over the temperature range. Based on this analysis, both of these earlier calculations are not sophisticated enough because they exclude resonances and produce an underestimated effective collision strength for the two transitions. The same trend is observed for the two calculations with our Rmatrix results for the 1–10 transition. The K88 results are employed by the Chianti database because they include the resonance effects. However, as can be seen in the inset panel of Fig. 4, the K88’s effective collision strength for the transitions 1–9 is lower than our results by about 30% and lower than the LB08 results. For most transitions, the effective collision strengths of K88 are lower than our calculation, but higher than that of LB08, which clearly indicates that the resonance contributions are underestimated.
Effective collision strength for the first 10 levels in S xiii at three temperatures (log _{10}T_{e} (K) = 6.23, 6.53, and 6.93).
In Fig. 5, we consider the collision strengths for the optically forbidden transition 1–4 (2s^{2}^{1}S_{0}–2s2p ^{3}P_{2}), the spinexchange transition 1–11 (2s^{2}^{1}S_{0}–2s3s ^{3}S_{1}), and the two optically allowed transitions 1–13 (2s^{2}^{1}S_{0}–2s3p ^{1}P_{1}) and 2–7 (2s2p ^{3}P_{0}–2p^{2}^{3}P_{1}). Here, we again compare our collision strength values for the four transitions with the two previous DW calculations of ZS92 and LB08. It is clear from Fig. 5 that our background collision strengths for all four transitions agree well with them. The insets in Fig. 5 show that the resonance effects dominate the collision strength in the lowenergy region and are larger than the background level of the collision strength in the highenergy region. Our results agree excellently with the previous DW results for all transitions with a few exceptions. The resonant effects can be up to 2–3 orders of magnitude except for transition 2–7, which is only a few times higher than the background.
The corresponding effective collision strengths for the four transitions (see Fig. 5) are plotted in Fig. 6. For the optically forbidden transition 1–4, our results are systematically higher than those of LB08 and K88 across the entire temperature range because of the clear resonance effects and a slightly higher background in the highenergy range. The resonance effects mainly enhance the effective collision strength for the optically allowed transition 1–13 over the temperature range 10^{4} K–10^{5} K. However, our effective collision strength for the optically allowed transition 2–7 agrees excellently with that of LB08. This is expected because of the weak resonance of the collision strength. For the spinexchange transition 1–11, the resonance contribution is extremely high at lower temperatures, and the effective collision strengths converge to the values of LB08 at high temperatures.
The differences between our Rmatrix effective collision strength and the results from LB08 for all transitions are also visually analysed in Fig. 7. We plot the ratios of effective collision strengths between our calculations and the results from LB08 for all 333 transitions available at three different temperatures (log _{10}T_{e} (K) = 6.23,6.53, and 6.93). One can see that our effective collision strengths are systematically higher than those of LB08, which can be explained by the resonance effects. There is generally a very good agreement between our collision strengths and the previous DW results. This agreement of Υ is better at higher temperatures (blue solid triangles) than at the lower temperatures (black solid squares), which is also expected from the behaviour of the collision strength due to the lack of resonances at higher incident electron energies. Table 4 contains a set of effective collision strengths for a sample of three electron temperatures for all transitions among the first ten finestructure levels arising from the n = 2 complex of S xiii. In Table 4, We compare our DARC effective strengths with the results of K88 and LB08. As shown in Table 4, our effective collision strengths show an average difference of 29.6%, 23.6%, and 14.3% at log _{10}T_{e} (K) 6.23, 6.53, and 6.93, respectively, with the DW calculations of LB08. Overall, our effective collision strengths are consistent and closer at higher temperatures with the DW calculations due to the absence of significant contributions of resonances. The background of our collision strength agrees excellently with the results of LB08. By further considering the accuracy of gfvalues involved in scattering calculations, we estimate our collision strength to be accurate to better than 20%. For effective collision strengths, our results are estimated to be accurate to better than 15%.
4. Conclusions
Electronimpact excitation data for Belike sulfur were calculated with the DARC code among the 98 finestructure levels (54 LS terms) belonging to the configurations 2s^{x}2p^{y} (x + y = 2) and 2lnl′ (l = s,p;n = 3,4;l′ = s,p,d and f). The target wavefunctions were generated by the GRASP0 code. Our energy levels, oscillator strengths, and transitions rates agree excellently with previous theoretical calculations and experimental measurements. The background values of our collision strength agree very well with previous DW calculations. The effective collision strengths were calculated over a temperature range 3.38 × 10^{4} K–3.38 × 10^{7} K important to astrophysically and plasma applications, and cover the whole equilibrium temperatures for Belike sulfur. Regarding the effective collision strengths, we found strong differences with values available in the literature because these exclude resonances. The resonance contributions to the effective collision strength can be as high as several times (2–10) for forbidden transitions, in particular at low temperatures. However, the resonance contribution is not significant for most of the optically allowed transitions, especially at high temperatures. In conclusion, we provided a robust and complete set of Rmatrix data of collision strengths and effective collision strengths for S xiii, which are useful for astrophysical and plasma applications.
Acknowledgments
The authors would like to thank C. P. Ballance for many fruitful discussions and valuable suggestions. This work was supported by the National Natural Science Foundation of China under grant Nos. 11273032, 11233004, 11243004, and 11103040. Dr. G. Y. Liang acknowledges the support from the OneHundredTalents programme of the Chinese Academy of Sciences (CAS). Dr. M. A. Bari acknowledges the support by the Chinese Academy of Sciences Fellowship for Young International Scientists.
References
 Acton, L. W., Bruner, M. E., Brown, W. A., et al. 1985, ApJ, 291, 865 [NASA ADS] [CrossRef] [Google Scholar]
 Badnell, N. R. 1986, J. Phys. B Atom. Mol. Phys., 19, 3827 [NASA ADS] [CrossRef] [Google Scholar]
 Ballance, C. P., & Griffin, D. C. 2006, J. Phys. B Atom. Mol. Phys., 39, 3617 [NASA ADS] [CrossRef] [Google Scholar]
 Bhattacharya, N., Bapat, B., Rangwala, S. A., Kumar, S. V. K., & Krishnakumar, E. 1998, Eur. Phys. J. D, 2, 125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bryans, P., Landi, E., & Savin, D. W. 2009, ApJ, 691, 1540 [NASA ADS] [CrossRef] [Google Scholar]
 Burgess, A., & Sheorey, V. B. 1974, J. Phys. B Atom. Mol. Phys., 7, 2403 [NASA ADS] [CrossRef] [Google Scholar]
 Du, S.B., Yang, Z.H., Chang, H.W., & Su, H. 2005, Chin. Phys. Lett., 22, 1638 [NASA ADS] [CrossRef] [Google Scholar]
 Eissner, W., Jones, M., & Nussbaumer, H. 1974, Comput. Phys. Commun., 8, 270 [NASA ADS] [CrossRef] [Google Scholar]
 Feldman, U. 1992, Phys. Scr., 46, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Goldsmith, S., Oren, L., & Cohen, L. 1974, ApJ, 188, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Grant, I. P., McKenzie, B. J., Norrington, P. H., Mayers, D. F., & Pyper, N. C. 1980, Comput. Phys. Comm., 21, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Kavanagh, A. P., Watanabe, H., Li, Y. M., et al. 2010, Phys. Rev. A, 81, 022712 [NASA ADS] [CrossRef] [Google Scholar]
 Keenan, F. P. 1988, Phys. Scr., 37, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Keenan, F. P., & Berrington, K. A. 1985, Sol. Phys., 99, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Keenan, F. P., McCann, S. M., & Widding, K. G. 1988, Sol. Phys., 117, 69 [NASA ADS] [CrossRef] [Google Scholar]
 Khardi, S., BuchetPoulizac, M. C., Buchet, J. P., et al. 1994, Phys. Scr., 49, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Landi, E., & Bhatia, A. K. 2008, Atomic Data and Nuclear Data Tables, 94, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Landi, E., Doron, R., Feldman, U., & Doschek, G. A. 2001, ApJ, 556, 912 [NASA ADS] [CrossRef] [Google Scholar]
 Landi, E., Del Zanna, G., Young, P. R., Dere, K. P., & Mason, H. E. 2012, ApJ, 744, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Lepson, J. K., Beiersdorfer, P., Behar, E., & Kahn, S. M. 2005, ApJ, 625, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Li, F., Liang, G. Y., & Zhao, G. 2013, ApJ, 762, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Mania, A. J., Luna, F. R. T., & Hernandes, J. A. 2009, J. Quant. Spec. Radiat. Transf., 110, 82 [NASA ADS] [CrossRef] [Google Scholar]
 Muehlethaler, H. P., & Nussbaumer, H. 1976, A&A, 48, 109 [NASA ADS] [Google Scholar]
 Norrington, P. H., & Grant, I. P. 1987, J. Phys. B Atom. Mol. Phys., 20, 4869 [NASA ADS] [CrossRef] [Google Scholar]
 Pegg, D. J., Forester, J. P., Vane, C. R., et al. 1977, Phys. Rev. A, 15, 1958 [NASA ADS] [CrossRef] [Google Scholar]
 Raassen, A. J. J., Mewe, R., Audard, M., et al. 2002, A&A, 389, 228 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sampson, D. H., Goett, S. J., & Clark, R. E. H. 1984, Atomic Data and Nuclear Data Tables, 30, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Schippers, S., Bernhardt, D., Müller, A., et al. 2012, Phys. Rev. A, 85, 012513 [NASA ADS] [CrossRef] [Google Scholar]
 Träbert, E., & Heckmann, P. H. 1980, Phys. Scr., 22, 489 [NASA ADS] [CrossRef] [Google Scholar]
 Tully, J. A., Seaton, M. J., & Berrington, K. A. 1990, J. Phys. B Atom. Mol. Phys., 23, 3811 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, Z.H., Du, S.B., Chang, H.W., et al. 2005, Chin. Phys. Lett., 22, 1099 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, H. L., & Fontes, C. J. 2013, Atomic Data and Nuclear Data Tables, 99, 416 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, H. L., & Sampson, D. H. 1992, Atomic Data and Nuclear Data Tables, 52, 143 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Level energies (Ryd) of S xiii from different calculations, along with the compilations of experimental values and other theoretical values.
Comparison of weighted oscillator strengths and transition rates (Avalues) with other theoretical and/or experimental results.
Effective collision strength for the first 10 levels in S xiii at three temperatures (log _{10}T_{e} (K) = 6.23, 6.53, and 6.93).
All Figures
Fig. 1 Ratio between our level energies and previous calculations: black triangle – DW calculations of Landi & Bhatia (2008), blue circle – experimental values, red cross – both the NIST’s MCDHF compilations and MCHF results from Mania et al. (2009). 

Open with DEXTER  
In the text 
Fig. 2 Ratio between our gfvalues and other theoretical calculations: open circle – DW calculations of Landi & Bhatia (2008), solid triangle – NIST’s MCDHF collections. 

Open with DEXTER  
In the text 
Fig. 3 Collision strength for the 1–5 (2s^{2} [entity ! #x20 !] ^{1}S_{0}–2s2p [entity ! #x20 !] ^{1}P_{1}) transition: solid black curve – our calculation without “topup”, dashed red curve – our calculation with “topup”, solid circle – DW calculations of Landi & Bhatia (2008). 

Open with DEXTER  
In the text 
Fig. 4 Collision strength 1–9 and 1–10: solid black curve – our DARC results, solid blue circle – DW calculation of Zhang & Sampson (1992), open red circle – DW calculation of Landi & Bhatia (2008). The insets in the graphs compare the corresponding effective collision strength (Υ) as a function of the electron temperature: open black circle – our DARC results, open star – Rmatrix interpolated work of Keenan (1988), open upward triangle – DW calculations from Landi & Bhatia (2008). 

Open with DEXTER  
In the text 
Fig. 5 Collision strengths for the optically allowed transitions 1–13 and 2–7, optically forbidden transition 1–4, and the intercombination transition 1–11: solid curve – our DARC calculation, open red square – DW calculations from Landi & Bhatia (2008), solid blue circle – calculations from Zhang & Sampson (1992). The insets show the details of the collision strengths near the background. 

Open with DEXTER  
In the text 
Fig. 6 Comparison of our DARC effective collision strengths: open circle – our calculation, open red upward triangle – DW calculation of Landi & Bhatia (2008), open blue star – interpolated Rmatrix values of Keenan (1988). 

Open with DEXTER  
In the text 
Fig. 7 Comparison of our effective collision strengths with those of Landi & Bhatia (2008): dashed lines indicate agreement within 20%, dotted lines indicate agreement within a factor of 2. 

Open with DEXTER  
In the text 