Oscillator strengths and effective collision strengths for electron excitation of Mg VI^{⋆}
Department of Physics, Clark Atlanta University, Atlanta, GA 30314, USA
email: stayal@cau.edu
Received: 6 June 2012
Accepted: 2 October 2012
Aims. Electron impact excitation collision strengths and oscillator strengths for the astrophysically important lines in Mg VI are reported. Thermally averaged collision strengths are presented as a function of electron temperature for application to solar and other astrophysical plasmas.
Methods. The collision strengths were calculated in a closecoupling approximation using the Bspline BreitPauli Rmatrix method. The multiconfiguration HartreeFock method with termdependent, nonorthogonal orbitals was employed for an accurate representation of the target wave functions. The closecoupling expansion includes 74 bound levels of Mg VI covering the n = 2 and n = 3 terms. The present calculations led to a total of 2701 forbidden, intercombination, and allowed transitions between finestructure levels. The effective collision strengths were obtained by averaging the electron collision strengths over a Maxwellian distribution of velocities.
Results. The accuracy of present oscillator strengths is evaluated by the agreement between the length and velocity formulations combined with the agreement between the calculated and measured excitation energies. The calculated excitation energies are in excellent agreement with experiments and other extensive configurationinteraction calculations. The oscillator strengths for all E1 transitions are listed. The effective collision strengths are tabulated for all 2701 transitions among the 74 finestructure levels at 11 electron temperatures in the range from 10 000 to 200 000 K. The present results are compared with other closecoupling and distortedwave calculations. Overall good agrement is generally found with the calculations by Ramsbottom & Bell and by Landi & Bhatia for many transitions, but significant differences are also noted for some transitions.
Key words: atomic data / atomic processes / line: formation
Tables 3 and 4 are available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/548/A27
© ESO, 2012
1. Introduction
The prominent Mg VI emission lines have been observed in the ultraviolet (UV), extreme ultraviolet (EUV), and soft Xray spectral ranges in solar and various other astrophysical plasmas. The intensity maps for Mg VI and Fe VII ions are similar, indicating a similar temperature of formation. The maximum ionic abundance of Mg VI in ionization equilibrium of astrophysical plasmas is at around 4 × 10^{5} K. The Mg VI lines due to the –2s2p^{4} ^{2}P_{1/2,3/2}, 2s2p^{4} ^{2}D_{5/2,3/2} transitions at 268.99 and 270.40 have been observed with Hinode/EIS (Young et al. 2007; Brooks et al. 2009) and with SOHO/CDS at transition region temperatures below 10^{6} K. The Mg VI forbidden lines are identified at 1190.09 and 1191.64 Å in the spectra of solar corona by SOHO/SUMER instrument. Several lines of Nlike Mg VI due to the 2s^{2}2p^{3}–2s2p^{4} transitions have been observed with SERTS including lines at 349.11 and 349.18 Å in the EUV spectrum of solar active region (Thomas & Neupert 1994). The Mg VI lines have also been observed in the spectrum of symboitic stars AG Dra and RR Tel by Hubble/STIS instrument (Young et al. 2006).
The Nlike ions can provide some of the best density diagnostics of solar transition region in both quiet and active conditions using intensity ratios from n = 2 forbidden and allowed lines. The forbidden lines among the ground configuration occur in UV spectral region and the allowed lines among the n = 2 levels occur in the EUV spectral region. The Mg VI line intensity ratios between the n = 2 and n = 3 levels are temperature sensitive and can be used as plasma temperature diagnostics over a wide temperature range of the solar corona. There are inconsistencies between the observed and theoretical predictions for line intensities of Nlike Mg VI and Si VIII ions due to the atomic data for these ions (Bhatia & Landi 2003; Young et al. 1998). For example, there are discrepancies between theory and observations from the SERTS for the Mg VI 2s^{2}2p^{3} ^{2}D°–2s2p^{4} ^{2}P, ^{2}D and 2s^{2}2p^{3} ^{4}S°–2s2p^{4} ^{4}P lines which have been attributed to the inaccuracies in the atomic data (Young et al. 1998). In order to make full use of high quality EIS spectroscopic measurements, accurate atomic data of matching quality are needed (Young & Landi 2009; Brooks et al. 2009). Radiative transitions and electronion collisions play an important role in the understanding of the physical processes and conditions in various types of astrophysical plasmas including solar transition region in both quiet and active conditions.
Several theoretical calculations of collision strengths and rates have been performed using distortedwave and closecoupling methods with different approximations for the target wave functions and with different number of target states. To our knowledge, there is no experimental study available for electron impact excitation collision data of Mg VI. The experimental data are important to verify or to calibrate theory. In astrophysical plasma applications, however, the majority of required radiative decay rates and collisional rate coefficients are derived from the calculations. Previous theoretical calculations for electron impact excitation collision strengths include the distortedwave calculations of Bhatia & Mason (1980), Bhatia & Young (1998), Landi & Bhatia (2007), and Zhang & Sampson (1999) and in the Rmatrix calculations of Ramsbottom & Bell (1997). The distortedwave calculations were carried out at a few selected incident electron energies and did not include resonances in the collision strengths. Bhatia & Young (1998) and Ramsbottom & Bell (1997) considered the 12 lowest LS states of the 2s^{2}2p^{3}, 2s2p^{4}, 2p^{5}, and 2s^{2}2p^{2}3s configurations, while Landi & Bhatia (2007) included 74 finestructure levels of the n = 2 and n = 3 states in their calculations.
Previous theoretical calculations have basically been performed in the LScoupling approximation and the collision strengths for finestructure transitions were determined by transforming scattering matrices in LS coupling to LSJ coupling (Saraph 1978). Besides these, the closecoupling expansions did not include all possible n = 3 terms which may be important for the convergence of the closecoupling expansion. Ramsbottom & Bell (1997) included only the four 2s^{2}2p^{2}3s states of the n = 3 complex in their 12state Rmatrix calculation. The purpose of the work reported here is to provide a set of accurate oscillator strengths and effective collision strengths for electron excitation of Mg VI for applications to astrophysical plasmas. Our scattering calculations were performed by using highly accurate target wave functions and by including finestructure effects in the closecoupling expansions directly to make sure that the relativistic effects are adequately accounted for in the scattering calculations. The present calculations were carried out with the Bspline BreitPauli Rmatrix (BSR) approach (Zatsarinny & Tayal 2001; Zatsarinny 2006). One of the distinct features of the present method is the use of termdependent, nonorthogonal orbital sets in the description of the target states. This allows us to optimize the atomic wave functions for different states independently, resulting in a more accurate target description than those used in previous collision calculations.
2. Computational methods
We used nonorthogonal orbitals to represent different Mg VI states. The nonorthogonal orbitals provide much greater flexibility in the choice of wave functions than the orthogonal orbitals and also allow us to include correlation with a reasonable number of configurations and correlated orbitals. Our calculations were performed using the multiconfiguration HartreeFock (MCHF) method (Froese Fischer 2007; Zatsarinny & Froese Fischer 2000). In the MCHF approach the atomic state is represented by an atomic state function (1)where the configuration state functions (CSF) Φ(α_{i}LS) are constructed from oneelectron functions and α_{i} defines the coupling of angular momenta of the electrons.
The Mg VI wave functions exhibit significant correlation corrections and term dependence of the oneelectron orbitals. The nonorthogonal orbitals were optimized for each atomic state separately. We began with the HarteeFock (HF) calculation for the 1s, 2s, and 2p orbitals for the ground 2s^{2}2p^{3} configuration terms and then determined separate sets of excited orbitals for various symmetries of even and odd parities. A set of correlation s, p, d, f, and g orbitals was determined by optimization on the ground 2s^{2}2p^{3} state and another set of s, p, d, f, and g correlation orbitals was optimized on the 2s^{2}2p^{2}3s ^{4}P state. The mean radii of the correlation orbitals are comparable to the spectroscopic orbitals and thus the correlation corrections are represented very well. The spectroscopic and correlation functions are used to construct CI expansions for different atomic states by allowing oneelectron and twoelectron excitations from all the basic configurations 2s^{2}2p^{3}, 2s2p^{4}, 2s^{2}2p^{2}3s, 2s^{2}2p^{2}3p, and 2s^{2}2p^{2}3d to the spectroscopic and correlation orbitals used in our calculations. In the construction of CI expansions for finestructure levels with various J and π, we used configurations generated in this excitation scheme for the atomic LS states and with insignificant configurations with coefficients less than 0.0008 omitted from the expansions. We retained 3314 and 2267 configurations for even and odd parity levels, respectively, in our calculations of oscillator strengths.
We further omitted configurations with weights less than 0.01 in calculations for the description of target wave functions to make the subsequent scattering calculations computationally feasible. We included 226 and 272 configuration state functions in the CI expansions of levels of odd parity and even parity respectively in our final calculation. We used 20, 69, 48, 13, 29, 17, and 30 CSFs respectively for the representation of the , , , , , , and terms of odd parity. For the ^{2}S, ^{2}P, ^{2}D, ^{2}F, ^{2}G, ^{4}P, ^{4}D, and ^{4}F even parity terms we used 31, 74, 51, 33, 10, 44, 15, and 14 CSFs respectively in the CI expansions. The valenceshell correlation effects due to electron excitations to the 2p^{4}3s, 2p^{4}3p, and 2p^{4}3d configurations are essential for the accurate representation of target states. The 2s electron excitation from the 2s^{2}2p^{3} ground configuration to the 3s, 3p, and 3d spectroscopic orbitals is very important. The configurations with oneelectron virtual excitations to the correlation s, p, and d orbitals are also important. Our calculation takes the important correlation and relaxation effects and of the term dependence of the valence orbitals into account.
For the scattering calculations we employed the BSR code (Zatsarinny 2006). Details of this approach can be found in our recent publications on e − K II (Tayal & Zatsarinny 2008) and e − S II (Tayal & Zatsarinny 2010). The distinctive feature of the method is the use of Bsplines as a universal basis to represent the scattering orbitals in the inner region of r ≤ a. Hence, the Rmatrix expansion in this region takes the form (2)\arraycolsep1.75ptHere are the channel functions while the splines B_{j}(r) represent the continuum orbitals. The χ_{i} are (N + 1)electron configurations formed from the oneelectron orbitals used to describe the Nelectrontarget states. These must be included to ensure completeness of the total trial wave function if the continuum orbitals are constructed orthogonally to the bound orbitals. The expansion coefficients a_{ijk} and b_{ik} are determined by diagonalizing the (N + 1)electron BreitPauli Hamiltonian with added Bloch operator. The relativistic effects in the scattering calculations have been incorporated in the BreitPauli Hamiltonian through the use of Darwin, mass correction, and spinorbit operators.
The present scattering approach imposes only limited orthogonality conditions to the continuum orbitals. The use of nonorthogonal orbitals allows us to reduce or even to avoid the introduction of additional (N + 1)electron terms in the Rmatrix expansion. We only require the orthogonality of continuum orbitals to the bound orbitals in the filled 1s shell. No orthogonality constraints to the spectroscopic excited orbitals or the correlated orbitals were imposed. In order to represent the continuum orbitals in internal region, 48 Bsplines of order 8 were used in the present calculations, with the size of the Bspline box a = 7.25 a_{0} (a_{0} = 0.529 × 10^{10} m denoting the Bohr radius). The scattering parameters are found by matching the inner solution at r = a with asymptotic solutions in the outer region. The Seaton’s STGF program as extended by Badnell (1999) has been employed to find the asymptotic solutions and collision parameters. We numerically calculated partialwave contributions up to J = 30. These partial waves are sufficient to achieve convergence for the forbidden transitions under consideration. Higher partial wave contributions are needed for the dipoleallowed transitions and these are estimated in a topup procedure based on the CoulombBethe method (Burgess & Sheorey 1974) or on geometric series approximation. In order to resolve the resonance structure with high accuracy, we used a fine energy mesh of 2.0 × 10^{4} Ry in the closedchannel energy region up to 6.0 Ry and then an energy mesh of 1.0 × 10^{3} Ry up to 10.0 Ry, that allowed us to resolve a vast majority of the narrow resonances in the energy range up through the highest excitation threshold. At higher energies up to 70 Ry where all channels are open and there are no resonances, the cross sections show smooth behavior and we used an energy grid of 0.25 Ry.
In astrophysical applications it is convenient to use excitation rate coefficients or thermally averaged collision strengths as a function of electron temperature. The excitation rates are obtained by averaging collision strengths over a Maxwellian distribution of electron energies. The excitation rate coefficient for a transition from state i to state f at electron temperature T_{e} is given by (3)where g_{i} is the statistical weight of the lower level i, ΔE_{if} = E_{f} − E_{i} is the excitation energy and γ_{if} is a dimensionless quantity called effective collision strength given by (4)where E_{f} is the energy of incident electron with respect to the upper level f. If the collision strength is assumed to be independent of the incident electron energy, we have γ_{if} = Ω_{if}. The effective collision strengths are calculated by integrating collision strengths for finestructure levels over a Maxwellian distribution of electron energies. The integration in Eq. (4) should be carried out using energy dependent collision strengths from threshold to infinity. The collision strengths vary smoothly in high energy region and the energy dependence of collision strengths for high energies have been properly accounted for by using numerical extrapolation technique at E > 70 Ry. In the asymptotic region, the collision strengths follow a high energy limiting behavior for the dipoleallowed transitions (5)where S is the line strength. The collision strengths vary smoothly in the high energy region and exhibit an increasing trend for dipoleallowed transitions. The collision strength increases more rapidly for the stronger dipoleallowed transitions than the weaker transitions.
3. Results and discussion
The target states included in the present scattering calculations are given in Table 1, where we compare our calculated excitation energies with the available experimental values and with recent calculations by Froese Fischer & Tachiev (2004) and Landi & Bhatia (2007). The overall agreement between experiment and present calculation is very satisfactory with an average deviation of about 0.034 Ry. The present excitation energies obtained from wave functions used in the scattering calculation show an average deviations of 0.019, 0.047, 0.053, and 0.017 Ry for the levels of 2s^{2}2p^{3}, 2s2p^{4}, 2p^{5}, and 2s^{2}2p^{2}3l configurations respectively. The accuracy of excitation energies from our extensive MCHF calculation for transition rates is somewhat better with average deviation of about 0.01 Ry. The present structure description represents a substantial improvement over those used in previous Rmatrix calculations by Ramsbottom & Bell (1997) and distortedwave calculations by Landi & Bhatia (2007). We obtained good agreement for all 2s^{2}2p^{3}, 2s2p^{4}, 2p^{5}, and 2s^{2}2p^{2}3l states, indicating that we consistently take correlation corrections in the ground and excited states into account, as our method allows us to use nonorthogonal orbitals, that are specific for given configurations. The accuracy of the present calculations is comparable to the BreitPauli MCHF calculation by Froese Fischer & Tachiev (2004).
Excitation energies (in Ry) of the 74 target levels included in the present transition rates and scattering calculations.
We also present the lifetimes of excited levels in Table 1, where our results are compared with the results of Froese Fischer & Tachiev (2004). The present lifetimes are within 10% of the results of Froese Fischer & Tachiev (2004) for most of the levels. However, larger differences are found for the levels of the 2s^{2}2p^{2}3p , , 2s^{2}2p^{2}3d ^{4}D, ^{2}F, and ^{2}G terms. The present lifetimes were calculated from transition probabilities of all possible E1 transitions between the states considered in our calculation. However, we have not included transition probabilities for forbidden E2 and M1 transitions and, therefore, the lifetimes of the metastable 2s^{2}2p^{3} and levels are not given in Table 1. This may also contribute to some of the discrepancies with the work of Froese Fischer & Tachiev (2004), especially for levels which decay primarily via forbidden transitions. The lifetimes of the 2s2p^{4}, 2p^{5} and 2s^{2}2p^{2}3l configurations varies from a few nanoseconds to fraction of nanosecond with significant lifetimes for the 2s^{2}2p^{2}3p , , , , and 2s^{2}2p^{2}3d ^{2}G levels.
The accuracy of oscillator strengths to some extent relates to the quality of target wave functions and, therefore, to the accuracy and reliability of collision strengths and rate coefficients. For example, collision strengths of the dipoleallowed transitions at high electron energies are directly proportional to the oscillator strengths. We have presented our length values of oscillator strengths and transition probabilities in Table 2 where present results are compared with the calculations of Froese Fischer & Tachiev (2004) and Landi & Bhatia (2007). We used experimental energies in the calculations of oscillator strengths and radiative rates. A rather reasonable agreement between different calculations is noted. Our calculation shows somewhat better agreement with the results of Froese Fischer & Tachiev (2004). We present oscillator strengths and transition probabilities for allowed and intercombination lines between the 74 finestructure levels in Table 3. The relativistic effects were included through the spinorbit, spinotherorbit, spinspin, mass and Darwin BreitPauli operators. We have presented oscillator strengths in both length and velocity formulations and transition probabilities in length form in Table 3. We have also listed wavelengths in this table. The values of oscillator strengths for intercombination lines are usually much lower than those for the dipoleallowed transitions. The agreement between the length and velocity forms of oscillator strengths may to some extent indicate the accuracy of wave functions and convergence of CI expansions. However, it is not a sufficient condition for the accuracy of results. The convergence of results is an important accuracy criterion. We have plotted ratio of the present velocity and length values as a function of length oscillator strengths in Fig. 1. It is clear from Fig. 1 that there is normally a very good agreement between the present length and velocity forms of oscillator strengths, particularly for transitions with significant strengths. The weak transitions, however, are not so accurate. The weak transitions belong mainly to the intercombination transitions. The intercombination transitions occur due to spinorbit mixing of total terms. This mixing strongly depends on the details of calculation and it is difficult to obtain the convergent results in some cases.
Comparison of oscillator strengths and transition probabilities (s^{1}) for some dipole radiative E1 transitions with the calculations of Landi & Bhatia (2007) (LB) and Froese Fischer & Tachiev (2004) (FT).
Fig. 1
Ratio between the present velocity and length oscillator strengths is shown as a function of length oscillator strengths. The dotted lines indicate 20% deviation. 

Open with DEXTER 
Fig. 2
Collision strengths for the forbidden 2s^{2}2p^{3}– (1–2) transition as a function of incident electron energy. Upper panel, resonant collision strengths in the closed channel region; lower panel, nonresonant collision strengths in all open channels region and present results (solid curve) are compared with the calculation of Landi & Bhatia (2007) (solid rectangles). 

Open with DEXTER 
Fig. 3
Collision strengths for the dipoleallowed 2s^{2}2p^{3} –2s2p^{4} (1–6) transition as a function of incident electron energy from threshold to 70 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER 
Fig. 4
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s2p^{4} (3–10) transition as a function of incident electron energy from threshold to 70 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER 
Fig. 5
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4–22) transition as a function of incident electron energy from threshold to 65 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER 
We begin our discussion of collision results for the excitation collision strengths for the forbidden 2s^{2}2p^{3} – (1 − 2) transition in both closedchannel and all openchannel energy regions in Fig. 2. The collision strengths show extensive resonance structures which are typical for lowenergy electron scattering on ions in closedchannel region. These resonances are expected to contribute significantly to the collision rates at low temperatures. Note also that our collision strengths show a richer resonance structure with many narrow resonances. It is partly because our results directly include the finestructure splitting both for the Nelectron target states and for (N + 1)electron resonances. Our calculation predicts very smooth behavior at the abovethreshold region for the forbidden 2s^{2}2p^{3} – (1–2) transition plotted in the lower panel of Fig. 2. Our results for nonresonant collision strengths have been compared with the available distortedwave calculation of Landi & Bhatia (2007) who reported results at five incident electron energies. The present results are shown by solid curve, while the solid rectangles display the results of Landi & Bhatia (2007). There is a good agreement between the two calculations with a maximum deviation of about 10%.
Fig. 6
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3d (1–55) transition as a function of incident electron energy from threshold to 55 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER 
Fig. 7
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3d (3–66) transition as a function of incident electron energy from threshold to 55 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER 
The collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s2p^{4} (1–6), 2s^{2}2p^{3}–2s2p^{4} (3 − 10), 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4–22), 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3d (1–55), and 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3d (3 − 66) transitions are displayed in both resonant and nonresonant regions in Figs. 3 − 7. Comparison of our results with the collision strengths from the distortedwave calculations by Landi & Bhatia (2007) shows an overall good agreement for the background collision strength of 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4–22) and 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3d (1–55) transitions, but some significant differences exist for the other three dipoleallowed 2s^{2}2p^{3}–2s2p^{4} (1–6), 2s^{2}2p^{3}–2s2p^{4} (3–10), and 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3d (3–66) transitions. The collision strengths for the dipoleallowed transitions show expected trend at higher energies. The collision strengths for stronger dipoleallowed transitions increase rapidly at higher energies than for the weaker dipoleallowed transitions. There is an agreement in shape among the two theories but some discrepancies exist in the magnitude of collision strengths. Our results are consistently lower than the calculation of Landi & Bhatia (2007) for the four out of five transitions shown in Figs. 3 − 7. Most likely, the discrepancies between the two calculations arise from the differences in target wave functions because of the strong configuration mixing.
Fig. 8
Effective collision strengths for the forbidden a) 2s^{2}2p^{3}– (1–2), b) 2s^{2}2p^{3}– (1–4), c) 2s^{2}2p^{3}– (4–5), and d) 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3s (1–20) transitions as a function of electron temperature. Present results, solid curve; Ramsbottom & Bell (1997), solid rectangles. 

Open with DEXTER 
Fig. 9
Effective collision strengths for the dipoleallowed a) 2s^{2}2p^{3}–2s2p^{4} (1–6), b) 2s^{2}2p^{3}–2s2p^{4} (3–9), c) 2s^{2}2p^{3}–2s2p^{4} (4–10), and d) 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4 − 22) transitions as a function of electron temperature. Present results, solid curve; Ramsbottom & Bell (1997), solid rectangles. 

Open with DEXTER 
Finally, in Figs. 8, 9 we compare the effective collision strengths as a function of electron temperature for the forbidden and dipoleallowed transitions with the Rmatrix calculations of Ramsbottom & Bell (1997). The effective collision strengths are calculated by averaging the total collision strengths over a Maxwellian distribution of electron energies. The effective collision strengths usually exhibit very smooth dependence on temperature despite strong resonance structures in collision strengths. The effective collision strengths for many transitions, however, are significantly enhanced due to presence of resonances at low electron energies. As seen from panels (a), (b), and (d) of Fig. 8, the present effective collision strengths for the 2s^{2}2p^{3}– (1–2), 2s^{2}2p^{3}– (4–5), and 2s^{2}2p^{3} –2s^{2}2p^{2}(^{3}P)3s (1–20) transitions agree well with the previous Rmatrix calculations at all temperatures except low temperatures. This discrepancy at lower temperatures is caused mainly by the differences in total collision strengths in the nearthreshold region. There are bigger differences among the two calculations for the 2s^{2}2p^{3}– (1–4), transition at all temperatures as shown in Fig. 8b. The differences in the background collision strengths and the resonance structures for the forbidden 2s^{2}2p^{3} – (1–4) transition produce noticeable differences in effective collision strengths. The bigger differences between the calculations are caused by the difference in target wave functions used in two calculations. The effective collision strengths for the dipoleallowed (a) 2s^{2}2p^{3}–2s2p^{4} (1–6); (b) 2s^{2}2p^{3}–2s2p^{4} (3–9); (c) 2s^{2}2p^{3}–2s2p^{4} (4–10); and (d) 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4–22) transitions as a function of electron temperature are displayed in Fig. 9 where our results are compared with the previous Rmatrix calculation of Ramsbottom & Bell (1997). Once again there is a reasonable agreement between the two calculations. It may be noted that Ramsbottom & Bell (1997) gave their results only for two significant digits and a smooth curve could not be obtained in Figs. 8, 9 from their reported results. The differences at lower temperatures are caused by the rich resonance structures with many narrow resonances resolved in our calculation with fine energy grid.
The effective collision strengths between all 74 finestructure levels of the ground 2s^{2}2p^{3} and excited 2s2p^{4}, 2p^{5}, 2s^{2}2p^{2}3l, and 2s2p^{3}3l configurations are presented in the Table 4 for electron temperatures from 10 000 to 200 000 K suitable for astrophysical applications. The indices of the lower and upper levels involved in transitions are given in Table 1. We have calculated all partial waves up to J = 30 in our Bspline BreitPauli Rmatrix approximation and used a topup procedure based on geometric series approximation for the forbidden transitions to estimate the contributions from higher partial waves. The topup procedure based on the CoulombBethe approach was used for the allowed transitions. In order to accurately calculate effective collision strengths at higher temperatures, we also extrapolated our collisions strengths to high energies. The effective collision strengths for the forbidden transitions display a very slow falloff with increasing temperature at higher temperatures. The effective collision strengths for many transitions are significantly enhanced because of the presence of resonances in the total collision strengths at low electron energies.
4. Summary
We have presented radiative parameters and effective collision strengths for all transitions among the 74 LSJ levels of Mg VI. This is most extensive and consistent data set for collisional and radiative parameters for Mg VI to date. The calculations were performed with the BSR code (Zatsarinny 2006) in which a Bspline basis is employed to represent the continuum orbitals. The use of nonorthogonal orbital sets, both constructing the target wave functions and for representing scattering functions, allowed us to optimize different atomic wave functions independently, hence to generate a more accurate description of the target states than those used in previous collision calculations. The radiative parameters for stronger transitions are in good agreement with the recent extensive CI calculations (Froese Fischer & Tachiev 2004).
The effective collision strengths were presented over a wide range of temperatures suitable for use in astrophysical plasmas modeling. Systematic comparisons with other available calculations were made. The overall agreement with previous Rmatrix
calculations of Ramsbottom & Bell (1997) is very reasonable and some discrepancies can be explained on the basis of more accurate target wave functions used in our calculations. Taking the good agreement for oscillator strengths into account, the accuracy for the strong dipoleallowed transitions is estimated to be within 10%. Overall, our results are estimated to be accurate to ~20% or better for most of the transitions.
Acknowledgments
This research work is supported by NASA grant NNX11AB62G from the Solar and Heliophysics program.
References
 Badnell, N. R. 1999, J. Phys. B, 32, 5583 [NASA ADS] [CrossRef] (In the text)
 Bhatia, A. K., & Landi, E. 2003, At. Data Nucl. Data Tables, 85, 317 [NASA ADS] [CrossRef] (In the text)
 Bhatia, A. K., & Mason, H. E. 1980, MNRAS, 190, 925 [NASA ADS] (In the text)
 Bhatia, A. K., & Young, P. R. 1998, At. Data Nucl. Data Tables, 68, 219 [NASA ADS] [CrossRef] (In the text)
 Brooks, D. H., Warren, H. P., Williams, D. R., & Watanabe, T. 2009, ApJ, 705, 1522 [NASA ADS] [CrossRef] (In the text)
 Burgess, A., & Sheorey, V. B. 1974, J. Phys. B, 7, 2403 [NASA ADS] [CrossRef] (In the text)
 Froese Fischer, C. 2007, Comput. Phys. Commun., 176, 559 [NASA ADS] [CrossRef] (In the text)
 Froese Fischer, C., & Tachiev, G. 2004, At. Data Nucl. Data Tables, 87, 1 [NASA ADS] [CrossRef] (In the text)
 Landi, E., & Bhatia, A. K. 2007, At. Data Nucl. Data Tables, 93, 1 [NASA ADS] [CrossRef] (In the text)
 Ramsbottom, C. A., & Bell, K. L. 1997, A&AS, 125, 543 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Saraph, H. E. 1978, Comput. Phys. Commun., 15, 247 [NASA ADS] [CrossRef] (In the text)
 Tayal, S. S., & Zatsarinny, O. 2008, Phys. Rev. A, 78, 012713 [NASA ADS] [CrossRef] (In the text)
 Tayal, S. S., & Zatsarinny, O. 2010, ApJS, 188, 32 [NASA ADS] [CrossRef] (In the text)
 Thomas, R. J., & Neupert, W. M. 1994, ApJS, 91, 461 [NASA ADS] [CrossRef] (In the text)
 Young, P. R., & Landi, E. 2009, ApJ, 707, 173 [NASA ADS] [CrossRef] (In the text)
 Young, P. R., Landi, E., & Thomas, R. J. 1998, A&A, 329, 291 [NASA ADS] (In the text)
 Young, P. R., Dupree, A. K., Espey, B. R., & Kenyon, S. J. 2006, ApJ, 650, 1091 [NASA ADS] [CrossRef] (In the text)
 Young, P. R., Zanna, G. D., Mason, H. E., et al. 2007, Publ. Astron. Soc. Japan, 59, 727 [NASA ADS] (In the text)
 Zatsarinny, O. 2006, Comput. Phys. Commun., 174, 273 [NASA ADS] [CrossRef] (In the text)
 Zatsarinny, O., & Froese Fischer, C. 2000, Comput. Phys. Commun., 124, 247 [NASA ADS] [CrossRef] (In the text)
 Zatsarinny, O., & Tayal, S. S. 2001, J. Phys. B, 34, 1299 [NASA ADS] [CrossRef] (In the text)
 Zhang, H. L., & Sampson, D. H. 1999, At. Data Nucl. Data Tables, 72, 153 [NASA ADS] [CrossRef] (In the text)
All Tables
Excitation energies (in Ry) of the 74 target levels included in the present transition rates and scattering calculations.
Comparison of oscillator strengths and transition probabilities (s^{1}) for some dipole radiative E1 transitions with the calculations of Landi & Bhatia (2007) (LB) and Froese Fischer & Tachiev (2004) (FT).
All Figures
Fig. 1
Ratio between the present velocity and length oscillator strengths is shown as a function of length oscillator strengths. The dotted lines indicate 20% deviation. 

Open with DEXTER  
In the text 
Fig. 2
Collision strengths for the forbidden 2s^{2}2p^{3}– (1–2) transition as a function of incident electron energy. Upper panel, resonant collision strengths in the closed channel region; lower panel, nonresonant collision strengths in all open channels region and present results (solid curve) are compared with the calculation of Landi & Bhatia (2007) (solid rectangles). 

Open with DEXTER  
In the text 
Fig. 3
Collision strengths for the dipoleallowed 2s^{2}2p^{3} –2s2p^{4} (1–6) transition as a function of incident electron energy from threshold to 70 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER  
In the text 
Fig. 4
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s2p^{4} (3–10) transition as a function of incident electron energy from threshold to 70 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER  
In the text 
Fig. 5
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4–22) transition as a function of incident electron energy from threshold to 65 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER  
In the text 
Fig. 6
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3d (1–55) transition as a function of incident electron energy from threshold to 55 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER  
In the text 
Fig. 7
Collision strengths for the dipoleallowed 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3d (3–66) transition as a function of incident electron energy from threshold to 55 Ry. Solid curve, present calculation; solid rectangles, calculation of Landi & Bhatia (2007). 

Open with DEXTER  
In the text 
Fig. 8
Effective collision strengths for the forbidden a) 2s^{2}2p^{3}– (1–2), b) 2s^{2}2p^{3}– (1–4), c) 2s^{2}2p^{3}– (4–5), and d) 2s^{2}2p^{3}–2s^{2}2p^{2}(^{3}P)3s (1–20) transitions as a function of electron temperature. Present results, solid curve; Ramsbottom & Bell (1997), solid rectangles. 

Open with DEXTER  
In the text 
Fig. 9
Effective collision strengths for the dipoleallowed a) 2s^{2}2p^{3}–2s2p^{4} (1–6), b) 2s^{2}2p^{3}–2s2p^{4} (3–9), c) 2s^{2}2p^{3}–2s2p^{4} (4–10), and d) 2s^{2}2p^{3}–2s^{2}2p^{2}(^{1}D)3s (4 − 22) transitions as a function of electron temperature. Present results, solid curve; Ramsbottom & Bell (1997), solid rectangles. 

Open with DEXTER  
In the text 