Articles

CHIANTI—AN ATOMIC DATABASE FOR EMISSION LINES. XIII. SOFT X-RAY IMPROVEMENTS AND OTHER CHANGES

, , , , and

Published 2013 January 11 © 2013. The American Astronomical Society. All rights reserved.
, , Citation E. Landi et al 2013 ApJ 763 86 DOI 10.1088/0004-637X/763/2/86

0004-637X/763/2/86

ABSTRACT

The CHIANTI spectral code consists of two parts: an atomic database and a suite of computer programs in Python and IDL. Together, they allow the calculation of the optically thin spectrum of astrophysical objects and provide spectroscopic plasma diagnostics for the analysis of astrophysical spectra. The database includes atomic energy levels, wavelengths, radiative transition probabilities, collision excitation rate coefficients, ionization, and recombination rate coefficients, as well as data to calculate free–free, free–bound, and two-photon continuum emission. Version 7.1 has been released, which includes improved data for several ions, recombination rates, and element abundances. In particular, it provides a large expansion of the CHIANTI models for key Fe ions from Fe viii to Fe xiv to improve the predicted emission in the 50–170 Å wavelength range. All data and programs are freely available at http://www.chiantidatabase.org and in SolarSoft, while the Python interface to CHIANTI can be found at http://chiantipy.sourceforge.net.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The CHIANTI database includes assessed atomic parameters and transition rates necessary for the calculation of the emissivity of optically thin, collisionally dominated plasmas. IDL-based software was initially developed to allow for synthetic spectra calculation and plasma diagnostics. In the first release (Dere et al. (1997), CHIANTI was designed to provide synthetic spectra at wavelengths longer than 50 Å. In subsequent releases, many additions have been made to the database to extend its wavelength coverage, improve the models of the ions included in CHIANTI, as well as add new physical processes and ions. Also, CHIANTI has been interfaced and included in the standard software of virtually all major solar missions, such as the Solar and Heliospheric Observatory (SOHO; Domingo et al. 1995), Transition Region and Coronal Explorer (TRACE; Handy et al. 1999), Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI; Lin et al. 2002), Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008), Hinode (Kosugi et al. 2007), and the Solar Dynamics Observatory (SDO; Pesnell et al. 2012). As part of the CHIANTI project, systematic comparisons between CHIANTI emissivities of individual ions or overextended wavelength ranges with observed solar spectra have been carried out to ensure the quality of the CHIANTI data and identify areas where improvement was needed. CHIANTI is one of the most widespread, complete, and accurate spectral codes available for optically thin emission in the 1–2000 Å wavelength range. At longer wavelengths the coverage is not as complete, but many important ions that emit at these wavelengths are also included.

The present paper describes the latest version of the CHIANTI database (version 7.1), which has been developed to include improved data for several ions, as well as greatly expanded models for Fe ions from Fe viii to Fe xiv. Updated recombination rates have also been included, and a new set of ionization equilibria have been calculated. Several new element abundance data sets have also been added. This version of CHIANTI has been designed to improve the database and synthetic spectrum in the 50–170 Å range, in order to meet the analysis requirements of the SDO instruments Extreme ultraviolet Variability Experiment and Atmospheric Imaging Assembly (AIA) 94 Å and 131 Å soft X-ray bands.

The CHIANTI database includes several files for each ion: the information on the atomic model and energy level is stored in the "ELVLC" file (e.g., si_9.elvlc for Si ix), wavelengths and radiative data are stored in the "WGFA" file (e.g., si_9.wgfa for Si ix), while electron–ion collision rates are stored in the "SPLUPS" file (e.g., si_9.splups for Si ix). We will use the ELVLC, WGFA, and SPLUPS conventions throughout the text.

2. DATABASE IMPROVEMENTS

2.1. Be-like Ions

2.1.1. O v

The previous model for O v included data for the entire n = 2 and n = 3 complexes for a total of 46 fine-structure levels. CHIANTI version 7.1 extends the model to include all the n = 4 and n = 5 configurations, for a total of 166 fine-structure levels.

Observed energy levels have been taken from version 5 of the National Institute of Standards and Technology (NIST; Kramida et al. 2012) for levels belonging to the n = 2, 3 complexes, and from the compilation of Moore (1993) for the levels in the n = 4, 5 configurations.

Theoretical energy levels, radiative transition rates, and collision strengths have been taken from Bhatia & Landi (2012), who used the Flexible Atomic Code (FAC; Gu 2003) to calculate the full set of data included in CHIANTI. Bhatia & Landi (2012) included all configurations of complexes up to n = 9 in the target model but the 2snl (n = 6–9) configurations have been omitted here. Levels from these configurations had lower energies than those in the 2p5l configurations, but their omission has no consequence either for the total radiative decay from the 2p5l levels, or on the overall O v population calculations under solar coronal conditions. Collision strengths were calculated at six values of the incident electron energy using the distorted wave (DW) approximation ranging from threshold to 41 Ry. However, resonant excitation has important effects on the O v Maxwellian-averaged collision strengths, so the R-matrix calculations included in the previous CHIANTI version (described in Landi et al. 2006), along with the energy separation in the SPLUPS file, have been retained for the n = 2, 3 levels.

The extended O v model allows the prediction of lines in the soft X-ray range between 110 Å and 170 Å which can be observed by high-resolution spectrometers and also fall within the passband of the SDO/AIA 131 Å channel.

2.2. B-like Ions

2.2.1. C ii

The previous model for C ii was based on the R-matrix collision strengths of Blum & Pradhan (1992) that consisted of 18 fine-structure levels. The current model is based on the calculations of Tayal (2008) that provide oscillator strengths, A-values, and electron effective collision strengths for an atom with 37 fine-structure levels. Observed energies and wavelengths are taken from the compilation of NIST (version 5; Kramida et al. 2012). In addition, five observed wavelengths between the 2s22p2P and the 2s2p2 4P levels are taken from the analysis of Young et al. (2011).

2.2.2. O iv

The electron collision data for O iv have been largely unchanged since CHIANTI 1 (Dere et al. 1997), for which R-matrix data of Zhang & Sampson (1994) were used for transitions among the n = 2 levels and unpublished data of D.H. Sampson for transitions involving n = 3 levels. These data are completely replaced here with the data from Aggarwal & Keenan (2008), who provided effective collision strengths for all transitions between the 75 levels of the configurations 2s22p, 2s2p2, 2p3, 2s23l (l = s, p, d), 2s2p3l (l = s, p, d), and 2s24l (l = s, p, d, f). The previous CHIANTI model contained 57 levels from the 2p23l (l = s, p, d) configurations, which are dropped from the new model. These levels have high energies (many are above the O iv ionization limit) and do not produce significant emission. Coupled with the fact that the data for these levels were never published, it has been judged better to remove them.

The Aggarwal & Keenan (2008) effective collision strengths were computed with the DARC code and are tabulated for 12 temperatures between 3 × 103 and 1 × 106 K. As the tabulated temperatures were not evenly distributed when projected onto a logarithmic scale, they were interpolated to provide additional values at 1.7 × 104 and 1.7 × 105 K. Nine-point spline fits were then performed to the effective collision strengths tabulated at 14 temperature points. For allowed transitions, high-temperature limit points (Burgess & Tully 1992) were derived using theoretical energies and oscillator strengths provided by Aggarwal & Keenan (2008). For a number of allowed transitions, the scaled effective collision strengths did not tend toward the high-temperature limit and in some cases it was necessary to remove one or two of the highest temperature points in order to get a good fit. For other transitions it was necessary to remove between one and three of the lowest temperature points to get a good fit. Care was taken to ensure that low-temperature points were not removed for transitions that are important for low temperature, photoionized plasmas, however. The fits reproduce the original data to better than 1.85% in the 104–5 × 105 K temperature range.

The only O iv levels with significant population for astrophysical applications are the energetically lowest five levels, and so only electron collision data for transitions involving these levels are included.

Three sources for radiative decay rates are used. For all transitions among levels of the 2s22p, 2s2p2, 2p3, 2s23l (l = s, p, d) configurations (levels 1–20 in the CHIANTI indexing), except the ground transition, data from Corrégé & Hibbert (2004) are used. The ground transition decay rate is from Tachiev & Froese Fischer (2000), although we note that this value is in very good agreement with the earlier calculations of Flower & Nussbaumer (1975) and Galavís et al. (1998). Rates for decays from the 2s2p3s4PJ levels (numbers 21–23 in the CHIANTI indexing) are also from Tachiev & Froese Fischer (2000). Decay rates for all other transitions are from Aggarwal & Keenan (2008); however, these values were computed using the authors' theoretical energies and so they have been scaled using experimental energies to yield more accurate values.

Experimental energies are available for all of the 75 levels of the CHIANTI model. The 2s22p2P3/2 level (number 2 in the CHIANTI model) energy is from Feuchtgruber et al. (1997), and the 2s2p2 4PJ (numbers 3–5 in the CHIANTI model) energies are from Young et al. (2011). All other energies are from version 4.1 of the NIST database (Ralchenko et al. 2011).

Giunta et al. (2012) identified a problem with the previous CHIANTI model whereby transitions from the 2s23s2S1/2 level were underpredicted. They demonstrated that using the Aggarwal & Keenan (2008) atomic data fixes this problem, and so the new CHIANTI model resolves this issue. This has consequences for studies of the λλ279.63, 279.93 emission lines in Hinode/EIS spectra (Culhane et al. 2007; Muglach et al. 2010) which are emitted from this level.

2.3. C-like Ions

2.3.1. N ii

The N ii lines observed in the visible provide a useful density diagnostic at densities above 104 cm−3. The previous model for N ii included 23 fine-structure levels and was based on the effective collision strength calculations of Stafford et al. (1994). The new model includes 58 fine-structure levels. Observed energies and wavelengths are taken from the compilation of NIST (version 5; Kramida et al. 2012). In addition, two observed wavelengths for transitions within the 2s22p2 configuration and one transition between the 2s22p2 and the 2s2p3 configurations are taken from the analysis of Young et al. (2011). The A-values are from Tayal (2011) with the exception of the forbidden transitions within the ground configuration, which are from Tachiev & Froese Fischer (2001). In addition, wavelengths for three transitions within the ground configuration are taken from Young et al. (2011). The effective collision strengths were calculated by Tayal (2011) over a temperature range of 5.0 × 102 K to 105 K.

2.4. N-like Ions

2.4.1. O ii

The O ii lines at 3729 Å and 3726 Å provide an important density diagnostic ratio that is often used in the analysis of astrophysical plasmas such a planetary nebulae (Wang et al. 2004). The previous model for O ii consisted of 15 fine-structure levels. The current model has been expanded to 35 fine-structure levels with observed energy levels and wavelengths taken from the NIST compilation (version 4; Ralchenko et al. 2011). The A-values for the ground configuration were calculated by Zeippen (1987) and those for allowed transitions were calculated by Tayal (2007). Effective collision strengths were calculated by Tayal (2007) over a temperature range from 2 × 103 K to 105 K. Following the suggestion of Kisielius et al. (2009), the collision strength indices for levels 4 and 5 were reversed.

2.4.2. Ca xiv

The radiative rates included in the previous version of CHIANTI for Ca xiv were taken from Bhatia & Landi (2005). These have been retained, but an error was corrected in the numbering of the levels in the WGFA file that caused the rates of levels 48 and 52 to be exchanged, and those of level 49 to be omitted from the file itself. All the rest of the data are unchanged.

2.5. Mg-like Ions

2.5.1. Ca ix

Ca ix data were left unchanged since the first release of CHIANTI in 1996, and consisted of a 16-level model where all data were taken from the DW calculation of Christensen et al. (1986). The entire data set has been superseded in CHIANTI version 7.1 by the FAC calculations of Landi & Bhatia (2012), who provided the energy levels, A-values, and collision strengths between the 283 fine-structure levels of the n = 3, n = 4, and n = 5 configurations.

Energy levels were taken from a number of sources (Parker & Phillips 1940; Fawcett 1970, 1976; Ekberg 1971; Edlen & Boden 1976; Litzen & Redfors 1987; Redfors 1988; Churilov et al. 1989; Curdt et al. 2004). Collision strengths were calculated at six values of the incident electron energy ranging from threshold to up to 57.9 Ry. Only the collision data for transitions originating in the lowest four levels in the atomic model (corresponding to the ground level and the 3s3p3P triplet) have been retained, as they were the only ones relevant for the level population calculation for astrophysical applications.

2.6. Iron Ions

The data available in the older versions of CHIANTI for Fe viii to Fe xiv came from a variety of sources, whose main focus was to calculate accurate atomic data and transition rates for transitions within the n = 3 complex, which gives rise to the very bright Fe lines dominating the solar and stellar spectra at wavelengths in the ≈170–450 Å wavelength range. Only a few data sets were available in the literature for transitions from the n = 4 and n = 5 complexes in these ions emitted at wavelengths between 50 Å and 170 Å. As a consequence, the CHIANTI spectrum for this range was clearly incomplete, as shown by e.g., Testa et al. (2012). In this version of CHIANTI, we have expanded the atomic models of all iron ions from Fe viii to Fe xiv to include a large number of levels from high-energy configurations in the n = 3–5 complexes, to account for the missing lines in the 50–170 Å range.

The CHIANTI models for Fe viii to Fe xiv have been deeply expanded, although almost all of the data from the previous CHIANTI versions was retained for the configurations emitting spectral lines in the 170–450 Å range. The strategy of this version of CHIANTI has been to add data for high-energy configurations to the already existing atomic models of these ions. The limited effect of resonant excitation for these configurations indicates that the use of the DW approximation is reasonable. The new DW calculations for Fe viii and Fe ix we include in version 7.1 have recently been presented in O'Dwyer et al. (2012). They carried out calculations using a new DW development of the autostructure code, described in Badnell (2011). Atomic data and transition rates for the high-energy configurations of Fe x to Fe xiv have been taken from Landi & Dere (2013), and are based on DW calculations carried out with the FAC code.

The wavelengths in the CHIANTI WGFA files are calculated using the experimental energies in the ELVLC files when available, otherwise theoretical energies are used. The latter provide wavelengths which can sometimes be several Å apart from their real value and so they are of limited use for spectroscopic analysis. To improve such wavelengths, theoretical energies are sometimes corrected using factors determined by inspecting the difference between the available experimental energies and their theoretical counterparts. These corrected energies are reported in the third pair of columns in the ELVLC files. Although this is an ad hoc procedures and correction factors may vary from configuration to configuration within the same ion, the accuracy of the wavelengths calculated using corrected energy levels is much improved. For the Fe ions we expanded in this version of CHIANTI, only a few experimental energy values were available for the new configurations added in the CHIANTI model. It was then decided to calculate correction factors for each ion, using as a guide the available experimental values. These factors have been either averaged (in case the differences were constant) or fitted with a straight line; the correction has been then put in the third energy columns and used to calculate improved energies. Table 1 reports the method used for each of the ions, and the levels from which it has been derived and to which it has been applied. Since no experimental energy values were available for the additional levels of Fe xiv, no correction was applied to the theoretical energies of this ion.

Table 1. Method of Energy Correction for Each of the Fe Ions Modified in CHIANTI Version 7.1

Ion Total No. of Levels Levels Method
Fe viii 104 3–55 Linear fit
    56–77 Linear fit
    78–104 Average
Fe ix 379 1–140 Average (version 7)
    141–379 Average
Fe x 825 4–48 Average
    49–825 Average
Fe xi 999 10–48 Linear fit
    49–999 Average
Fe xii 898 6–41 Linear fit
    42–898 Average
Fe xiii 950 29–114 Linear
    115–950 Average
Fe xiv 739   None

Download table as:  ASCIITypeset image

This new data set for Fe viii to Fe xiv vastly improve CHIANTI performances in the soft X-ray range between ≈80 Å and ≈120 Å. Figures 1 and 2 show a comparison between spectral observations and CHIANTI synthetic spectra calculated using version 7.0 and the present version. Observations come from the quiet-Sun spectrum from Manson (1972; see Del Zanna 2012a for details) and the Chandra spectrum of Procyon analyzed by Testa et al. (2012). In both cases, the new Fe calculations account for a large number of transitions giving rise to observed spectral lines in both spectra.

Figure 1.

Figure 1. Comparison between the quiet-Sun spectrum observed by Manson (1972, top panel) with synthetic spectra from versions 7.0 and 7.1 of CHIANTI (bottom panel). For details of the synthetic spectrum calculation, see Del Zanna (2012a).

Standard image High-resolution image
Figure 2.

Figure 2. Comparison between the spectrum of the solar-like star Procyon (top panel) with synthetic spectra from versions 7.0 and 7.1 of CHIANTI. The synthetic spectra calculation has been done using the DEM and parameters described in Testa et al. (2012).

Standard image High-resolution image

A future version of CHIANTI will include extensive R-matrix calculations (see, e.g., Del Zanna et al. 2012a), being completed under the auspices of the Atomic Processes for Astrophysical Plasmas network (http://www.apap-network.org/). These calculations will replace the DW data of all the levels added in this version of CHIANTI with R-matrix results, as well as improve on the already available R-matrix results for the lower energy configurations emitting in the 170–450 Å wavelength range, and add even more high-energy levels.

2.6.1. Fe viii

The version 7.0 data for levels in the 3s23p6nl (n = 5, 6, ln − 1) and 3s23p67l' (l' = s, p, d, f) configurations have been replaced with the calculations from O'Dwyer et al. (2012).

Experimental energy levels come from a combination of measurements from NIST (version 2; Martin et al. 1999), Ekberg & Feldman (2003), Young & Landi (2009), Del Zanna (2009), and Landi & Young (2010). These have been used to determine a correction to predicted energies which has been used to calculate improved wavelengths for transitions where experimental energy levels were not available.

The R-matrix calculation by Griffin et al. (2000) was retained for n = 3, 4 levels. As shown in Del Zanna (2009), significant limitations in the target adopted by Griffin et al. (2000) were present. To improve them, the collision strengths of the dipole-allowed transitions were scaled by Del Zanna (2009) according to the ratio of the gf values. These scaled values were also adopted in the previous version 7.0 of CHIANTI. Radiative data for the n = 3, 4 levels are taken from Del Zanna (2009).

Collisional and radiative data were calculated by O'Dwyer et al. (2012) for levels up to n = 7 using autostructure. These replace old calculations from Czyzak & Krueger (1966).

2.6.2. Fe ix

The CHIANTI model for Fe ix was expanded to include many more configurations and levels. The original model included the 3s23p6, 3s23p53d, 3s3p63d, 3s23p43d2, and 3s23p54s and 4p configurations for a total of 140 fine-structure levels, whose data were taken from the calculations of Storey et al. (2002). The present model adds the 3s3p53d2, 3p63d2, 3s23p5nl (n = 4–6, ln − 1), and 3s3p64l' (l' = s, p, d, f) configurations, reaching 379 fine-structure levels.

The data from Storey et al. (2002), already used in the previous version of CHIANTI, were retained for all levels, with the only exception of radiative data for the 3s23p54s and 4p configurations which are replaced by the Landi & Young (2009) calculation.

The Fe ix energy level and radiative data files were also updated with revised energies and wavelengths for transitions involving levels 3s23p4(1D)3d2 3D1, 2, 3 (with indices 110–112 in the CHIANTI model). The energies are from Young & Landi (2009) who gave new experimental energies for levels 110 and 111, and an updated theoretical energy for level 112. All other experimental energies for the new levels come from version 3 of NIST (Ralchenko et al. 2005). In the benchmark study of Del Zanna (2012a) some discrepancies between predictions and observations (from Manson 1972) in several Fe ix lines were noted, and alternative identifications have been proposed. New laboratory measurements of these lines will be needed to assess which identifications are correct. Observed energies have been used to calculate corrections to the energies of levels where no experimental entry was available; these corrected energies have been used to calculate improved values of unobserved wavelengths.

Radiative data and DW collision strengths have been taken from the autostructure calculations of O'Dwyer et al. (2012). Recently, also Foster & Testa (2011) carried out Fe ix DW calculations using the FAC code. The actual collision strengths obtained by O'Dwyer et al. (2012) with autostructure and those from FAC were very similar (within 10% in most cases), although some differences in the line intensities were present, due to different modeling.

2.6.3. Fe x

The new CHIANTI model for Fe x includes 825 fine-structure levels from 27 configurations. Energy levels, A-values, and collision data for the lower 54 levels come from version 7.0 of CHIANTI, and consist of the same R-matrix calculation from Del Zanna et al. (2004) described in Landi et al. (2006).

Data for the additional levels 55–825 were calculated by Landi & Dere (2013) and are from the 3s23p33d2, 3s23p23d3, 3s3p53d, 3s3p43d2, 3s3p33d3, 3p63d, 3p53d2, 3s23p44l, 3s23p33d4l, 3s3p54l, and 3s23p45l' (with l = s, p, d, f and l' = s, p, d, f, g). A much larger set of configurations was adopted by Landi & Dere (2013) for the Fe x target, which included 31,598 levels from the complete n = 3 complex, as well as from the n = 4, 5 and n = 6 complexes.

Observed energy levels are taken from Fawcett et al. (1972), with additional levels from Del Zanna et al. (2012a) and Del Zanna (2012a). Differences between observed and calculated levels have been used to determine a correction to the predicted energy of levels for which an experimental energy was not available; the resulting corrected energies were used to calculate improved wavelengths. A-values were calculated for all transitions within the 825 levels. Collision strengths were calculated for six values of the incident electron energy ranging from threshold to up to 202 Ry; we retain in CHIANTI the collisional data for transitions originating in the ground 3s23p5 2P doublet as well as from the metastable levels in the 3s23p43d configuration, whose population is large enough to influence the overall ion population in astrophysical plasmas.

2.6.4. Fe xi

The Fe xi model has been extended from 365 to 999 fine-structure levels with the addition of 3s23p3d3, 3s3p23d3, 3p43d2, 3s23p34l (l = p, d, f), 3s33p23d4s, and 3s3p44l' (l' = s, p, d, f). The lowest levels of some of these configurations had lower energies than the highest levels in the version 7.0 model, so that some renumbering was necessary. However, all the version 7.0 data available for the 365 levels in the version 7.0 model have been retained.

New data for the additional levels have been taken from the FAC calculations of Landi & Dere (2013). The Fe xi atomic target model adopted by Landi & Dere (2013) included the complete set of configurations in the n = 3 and n = 4 complexes, as well as the 3s23p35l configuration, for a total of 38,850 levels. The observed energies for these configurations were taken from Fawcett et al. (1972); additional energy levels have been taken from Del Zanna (2012a). Differences between observed and calculated level energies have been used to determine a correction to the predicted energy of levels for which an experimental value was not available; the resulting corrected energies were used to calculate improved wavelengths. A-values have been computed for all transitions within the additional levels, and between the latter and those in the configurations already included in version 7.0 of the database.

Collision strengths have been calculated at six values of the incident electron energy, ranging from threshold to up to 90 Ry using the DW approximation. Effective collision strengths have been calculated and retained only for transitions originating from the ground configuration and from the some metastable levels with significant population at log Ne = 12.0: 3s23p33d5D4, 3G4, 5, 1G4, (2D)3F4, and (2P)3F4.

2.6.5. Fe xii

The version 7.0 model of Fe xii included 143 levels from the lowest six configurations in the n = 3 complex. Version 7.1 extends it to 898 fine-structure levels by including additional levels from 28 additional configurations: 3s23d3, 3s3p23d2, 3s3p3d3, 3p43d,3p33d2, 3s23p24l, 3s23p3d4l, 3s3p34l (l = s, p, d, f), and 3s23p25l' (l' = s, p, d, f, g). The new Fe xii model retains all the data from version 7.0 for the lowest 143 levels, and uses the Landi & Dere (2013) calculations for the additional levels.

Landi & Dere (2013) used FAC on an Fe xii model which included the complete n = 3–5 complexes, for a total of 33,680 fine-structure levels. Experimental energies for the lowest 41 levels are the same as in version 7.0, coming from Del Zanna & Mason (2005); those for the higher-energy configurations are taken from Fawcett et al. (1972), with additional values from Del Zanna (2012a) and Del Zanna et al. (2012b). The differences between experimental and calculated energies were used to determine a correction to the energy of all other levels; corrected energies were used to calculate transition wavelengths. A-values were provided by Landi & Dere (2013) for all transitions included in the present version of CHIANTI.

Collision strengths were calculated at six values of the incident electron energy ranging from threshold to 141 Ry using the DW approximation. Effective collision strengths were calculated and retained for all levels within the ground configurations.

2.6.6. Fe xiii

The 7-configuration, 114-level model available in version 7.0 has been expanded to include data from 34 additional configurations: 3s3p3d2, 3s3d3, 3p23d2, 3p3d3, 3s23p4l, 3s23d4l, 3s3p24l, 3s3p3d4l, 3p34l, 3s23p5l', and 3s3p25l' with l = s, p, d, f and l' = s, p, d, f, g. The new CHIANTI model includes 950 fine-structure levels. The data in version 7.0 for the lowest 114 levels have been retained.

Experimental values for the new energy levels were taken from Fawcett et al. (1972); additional values were taken from Del Zanna (2012a) and Del Zanna & Storey (2012). A correction to the calculated energies was determined using the difference between the available experimental values and their theoretical counterpart; the corrected energies were used to provide more accurate wavelengths for the transitions involving the unobserved levels. The target ion model adopted by Landi & Dere (2013) included the entire n = 3–5 complexes, as well as the additional configurations 3s23pnl, 3s23dnl, and 3s3p2nl, with n = 6, 7 and ln − 1. A-values were provided for all transitions involving the 950 levels in the ELVLC file.

Collision strengths were calculated by Landi & Dere (2013) using the FAC code and the DW approximation. Six values of the electron energy ranging from transition threshold to 119 Ry. Effective collision strengths were calculated and retained for transitions originating in the five levels of the ground configuration, as well as from metastable level 18 (3s23p3d3F4).

2.6.7. Fe xiv

The version 7.0 model of Fe xiv, including data from 16 configurations and 197 fine-structure levels from Liang et al. (2010), has been expanded to include additional 542 levels from 39 additional configurations: 3s3p4f, 3s3d4l, 3p24l, 3s3d4l, 3s3p5l', 3s3p6s, and 3s2nl'', where l = s, p, d, f, l' = s, p, d, f, g, n = 4–7, and l'' ⩽ n − 1. The total number of levels in CHIANTI is now 739. All the data in the previous version of CHIANTI have been retained.

Experimental energy levels for the additional configurations were not available, while we added the value for levels 122, 125, and 136, where the new identifications proposed by Del Zanna (2012a) have been included. The lack of observed energies in the additional configurations has prevented the calculation of corrected wavelengths.

All data for the additional configurations come from Landi & Dere (2013), who used FAC with an Fe xiv atomic model that included the entire n = 3 to n = 7 complexes, for a total of 5079 levels. A-values were calculated for all transitions involving the new configurations. Collision strengths for the new levels were calculated by Landi & Dere (2013) using the DW approximation at six incident electron energies ranging from transition threshold to 256 Ry. Effective collision strengths were calculated and retained for transitions originating in the ground configuration.

3. IMPROVED RECOMBINATION RATES AND IONIZATION EQUILIBRIA

Recently, two new sets of calculations of radiative and dielectronic recombination rates have become available. Abdel-Naby et al. (2012) have provided rates for recombination of the aluminum-like ions to silicon-like ions and Nikolic et al. (2010) have provided rates for recombination of the argon-like ions to potassium-like ions. In both cases, these new rates generally supplant the rates of Shull & Van Steenberg (1982) that we have used in the past. Exceptions are the ions where the rates for Fe ix and Fe xiv were taken from Badnell (2006).

A revised version of the ionization equilibrium for all elements between H and Zn were then calculated at temperatures between 104 and 109 K. These are based on the ionization rates of Dere (2007), the new recombination rates and the rates listed by Dere et al. (2009). Figure 3 illustrates the differences for ions Ni ix through Ni xii between the current values (solid line) and those from version 6 (dashed line) resulting from the updated recombination rates for argon-like ion Ni xi.

Figure 3.

Figure 3. Ion fractions of Ni ix to Ni xii calculated using the ionization and recombination rates in CHIANTI. Solid line: present CHIANTI version (version 7.1). Dashed line: CHIANTI version 6 (Dere et al. 2009).

Standard image High-resolution image

4. IMPROVED AIA TEMPERATURE RESPONSE FUNCTIONS

For narrowband imaging instruments, such as SOHO/EIT, TRACE, STEREO/EUVI, Hinode/XRT, and SDO/AIA, it is crucial for scientists to have an accurate model of the spectral composition of the radiation that enters into the instrument. In particular, the spectral composition determines the temperature response functions (TRFs) that are needed for estimating physical quantities such as emission measure, temperature, and density. The TRF of an imager indicates, for every temperature T0, the amount of data numbers DN that the instrument returns when it observes an isothermal plasma at T0 per unit emission measure.

All the imaging instruments devoted to solar physics make use of CHIANTI to determine their TRFs, and the completeness of the database and the accuracy of its data are of critical importance. The bandpasses for SOHO/EIT, TRACE, STEREO/EUVI, and Hinode/XRT have been thoroughly tested by the CHIANTI team in several studies (Young et al. 1998; Landi et al. 2002; Landi & Phillips 2006; Chifor et al. 2007; Del Zanna 2012b), which have shown that CHIANTI synthetic spectra are both complete and accurate. However, Del Zanna et al. (2011) and Testa et al. (2012) showed that the bandpasses of the SDO/AIA 94 Å and 131 Å channels were incomplete, resulting in an underestimation of the TRF for those two channels. In an independent study, Aschwanden & Boerner (2011) also found that the SDO/AIA 94 Å TRF was underestimated and, with some assumptions on the thermal structure of the plasma, determined an empirical correction factor of the order 6.7 ± 1.7 to be multiplied to the standard TRF for temperatures lower than log T = 6.3.

Since one of the main goals of this version of CHIANTI is to address the shortcoming of the CHIANTI-based TRF for the AIA 94 Å and 131 Å channels, we have tested the effects of the new Fe calculations on the TRF of both bandpasses. The official AIA TRF is calculated by the AIA team; here we only show the new curves to illustrate the effects of the new CHIANTI data on those curves. The AIA effective areas were taken from SolarSoft. CHIANTI isothermal spectra were calculated per unit emission measure at a grid of temperatures in the log T = 4.0 = 8.0 range; both line and continuum emission have been included at a constant density of log Ne = 9.0. To verify that our software correctly calculated the AIA TRF curves, we first calculated them using version 6 of CHIANTI, the coronal element abundances of Feldman et al. (1992; file sun\coronal\ext.abund) and version 6.0 chianti.ioneq ionization equilibria. The TRF curves were in excellent agreement with those in the standard AIA software in SolarSoft; for details on the AIA temperature responses also see Del Zanna et al. (2011). We then used the present version of CHIANTI (including the revised ionization equilibria) and calculated new TRF curves. Results are shown in Figures 4 and 5 for the 94 Å and 131 Å channels, respectively; in each of these figures, the top panel shows the absolute values of the TRF as calculated with the AIA software (red line), and CHIANTI version 7.1 (blue line) with the same settings as in the AIA software. Figures 4 and 5 show that both channels are affected by the new data in the present version of CHIANTI. The difference is the 131 Å bandpass is due to changes in the equilibrium ion abundance for of Fe ions, especially for Fe viii, as well as new Fe xi lines from the new n = 4 levels added in this version of CHIANTI. The difference in the 94 Å bandpass is due to the presence of many more spectral lines in the present version of CHIANTI. Figures 6 and 7 show the comparison between the synthetic spectra of version 7.0 with those in version 7.1 under quiet Sun (with log Ne = 8.3) and active region (with log Ne = 8.7) conditions, respectively. New transitions from Fe ix, Fe xii, and Fe xiv, and the enhanced ion abundance of Fe viii (with some decrease in the ion abundance of Fe x) account for the large differences in the AIA TRF for channel 94 Å, as also discussed by Del Zanna (2012a).

Figure 4.

Figure 4. Comparison of SDO/AIA temperature response function (TRF) of the 94 Å channel calculated using version 7.1 of CHIANTI, with those available in the AIA software. Top: absolute values of the TRF. Bottom: ratio of TRF calculated using version 7.1 and the AIA software.

Standard image High-resolution image
Figure 5.

Figure 5. Same as Figure 4, for the 131 Å channel.

Standard image High-resolution image
Figure 6.

Figure 6. Synthetic spectra calculated from CHIANTI versions 7.0 (blue curve) and 7.1 (black curve) in the wavelength range covered by the AIA 94 Å channel for quiet-Sun conditions. Black labels indicate new lines included in version 7.1, blue labels indicate lines present in both CHIANTI versions. Additional lines from Fe xii and Fe xiv are predicted by version 7.1 only, at ≃ 93.6 Å, which contribute significantly to the total intensity of the 93.6 Å feature.

Standard image High-resolution image
Figure 7.

Figure 7. Same as Figure 6, for active region conditions.

Standard image High-resolution image

5. ELEMENT ABUNDANCES

Recently, new sets of element abundances have been proposed by different authors, and we expanded the available selection in CHIANTI to include them. A new set of coronal abundances has been taken from Schmelz et al. (2012), and two new sets of photospheric abundances come from Lodders et al. (2009) and Caffau et al. (2011). Caffau et al. (2011) provided new abundances for 12 elements, 8 of which are included in CHIANTI; they also recommended that the values of Lodders et al. (2009) be used for the remaining elements, so these two sets are combined in the CHIANTI file with the Caffau et al. (2011) file.

6. SUMMARY

In this paper we have summarized the new features of version 7.1 of the CHIANTI database. They consist of updated data for several ions, as well as greatly expanded data sets for several key Fe ions. CHIANTI version 7.1 is a significant improvement over the previous version in the 50–170 Å spectral range as it allows a more complete calculation of the synthetic spectrum by adding a large number of previously unavailable Fe lines. However, some of the problems of CHIANTI in this wavelength range outlined by Testa et al. (2012) due to a few ions of other elements than Fe still need being addressed. This will be done in the next version of CHIANTI. The changes we made in version 7.1 have significant effects in the TRF of the 94 Å and 131 Å channels of the SDO/AIA narrowband imager.

The entire database as well as the standard suite of IDL computer programs are distributed through SolarSoft; they are also freely available at the CHIANTI Web site http://www.chiantidatabase.org. ChiantiPy can be found at http://chiantipy.sourceforge.net.

The work of E.L. is supported by the NNX11AC20G and NNX10AQ58G NASA grants, and by NSF grant AGS-1154443. G.D.Z. acknowledges support from STFC via the Advanced Fellowship program and the UK APAP Network. The work of P.R.Y. was performed under contract with the Naval Research Laboratory and was funded by NASA. K.P.D. acknowledges funding from the NASA grant NNX11AC20G. H.E.M. acknowledges support from STFC. This research has made use of NASA's Astrophysics Data System. We warmly thank all the authors that provided us with the data necessary to make this version of the database possible, and all CHIANTI users that helped us improve the database and the software with their comments.

Please wait… references are loading.
10.1088/0004-637X/763/2/86