About Chandra Archive Proposer Instruments & Calibration Newsletters Data Analysis HelpDesk Calibration Database NASA Archives & Centers Chandra Science Links

Search In:   

Home

Features

Download

Sources

Physics

Issues
 

Comparing APEC to Other Plasma Codes


ATOMDB v1.2.0 vs. Raymond-Smith | ATOMDB v1.2.0 vs. Mekal/SPEX | ATOMDB v1.2.0 vs. ATOMDB v1.1.0

Differences between ATOMDB Version 1.2.0 and Raymond-Smith

The APEC outputs should be viewed as the replacement for the Raymond-Smith code. At high spectral resolution, the Raymond-Smith is not practicable because the emission lines are not calculated in fine structure (i.e. they are calculated via line bundling). Radiated power as a function of temperature (i.e. the cooling curves) are not significantly different. However, even at moderate (CCD) resolution, there are significant differences between APEC V1.2.0 and the Raymond-Smith model implemented in public fitting packages, due to different models for ionization balance, fine structure and atomic data. The most important differences are for iron, shown for binned spectra at different temperatures in these figures (ps): We compare element-by-element line spectra from APEC V.1.2.0 and SPEX V1.1.0. (Although newer versions of SPEX may be available, the V1.1.0 models are similar to the MEKAL model currently available in XSPEC and Sherpa.) These comparison figures (ps) provide valuable information for X-ray spectral analysis near 1 keV. Examples from these figures (below) illustrate the major differences.  Examples of noticeable differences between ATOMDB and MEKAL/SPEX at moderate (CCD) resolution

a. ATOMDB does not yet include fluorescence lines, whereas SPEX does. The most significant implication is that APEC models do not include the Fe K 6.4 keV (1.9 Angstrom) line (see other CAVEATS).

b. For Oxygen the most significant differences occur in the O VIII Ly alpha/ O VIII Ly beta ratio, as shown by Smith et al. (2001) ( ADS).

Comparison of Oxygen Line Spectra at 40 Million K

Oxygen Line Comparison
Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 

c. Other H-like ions show a similar difference in Ly alpha/Ly beta ratio. Silicon also shows large differences in He-like emissivities at low temperatures, presumably caused by differences in the ionization balance. SPEX assumes Arnaud and Rothenflug (1985, A&AS, 60, 425), while APEC assumes Mazzotta et al. (1998, A&AS, 133, 403).

Comparison of Silicon Line Spectra at 2.5 Million K

Silicon Line Comparison
Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 

d. Fe L-shell data show significant differences. Some effects can be attributed to the differences in ionization balance (SPEX assumes Arnaud and Raymond (1992, ApJ, 398, 394), while APEC assumes Mazzotta et al.). Other effects can be attributed to scaling laws used in SPEX, which are not used in APEC (see, e.g. Gu et al. 1999, ApJ, 518, 1002). For example, the large differences at low temperatures are from Fe XVII only, and must be related to the scaling laws. Differences at higher temperature appear to be a combination of scaling laws and ionization balance effects for Fe XXI to XXII. (The APEC calculations use the collision strengths produced by HULLAC and have been tested relative to the original calculations.)

Comparisons of Iron Line Spectra

Iron Line Comparison
Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Iron Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 

Examples of noticeable differences between ATOMDB and MEKAL/SPEX at high resolution

a. The DR satellite lines appearing near the resonance lines show significant spectral structure not captured by assumptions introduced in SPEX (e.g. near O VII 21.7 Angstroms at 3 million K and Fe XVII 15.0 Angstroms at 2 million K). Proper treatment of the DR lines is important for distinguishing between line and continuum emission (see Brickhouse et al. 1995, ApJS, 97, 551).

Comparison of DR Satellite Line Spectra

Oxygen Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Iron Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

b. Fine structure lines are split up in the ATOMDB, but not always in the SPEX list. For example, the O VIII H-alpha lines near 102.5 Angstroms should be treated carefully when trying to measure the density diagnostic Fe XXI 102.22 Angstrom line.

O VIII H alpha Fine Structure

Oxygen Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

c. Wavelength differences are also apparent. The ATOMDB has included the measurements of Brown et al. (1998, ApJ, 502, 1015) and Brown et al. (1999, Internal Report LLNL.JC-136647B). H- and He-like wavelengths in the ATOMDB are from relativistic quantum mechanics calculations and include fine structure.

d. SPEX includes inner shell excitation lines which are not currently included in the ATOMDB (see other CAVEATS), for example, the Fe line at 17.2 Angstroms; however, this line strength in SPEX is much stronger than in Raymond-Smith models and needs investigation. Note also wavelength differences noted in (c) above.

Fe XVII Comparison: Wavelengths and Inner Shell Excitation

Iron Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

e. The ATOMDB includes more lines at high principal quantum number than SPEX (see Brickhouse et al. 2000, ApJ, 530, 387). This is apparent near 10 Angstroms at 6 x 10^6 K). However, APEC models are still not necessarily complete enough in detail to achieve good spectral fits at high resolution, depending on the spectrum to fit. We do not believe APEC is complete enough yet at soft X-ray wavelengths to perform global fits to LETG spectra.

Iron Comparison: High n Lines

Iron Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Differences between ATOMDB Version 1.2.0 and ATOMDB Version 1.1.0

a. Major problems with Ni ion wavelengths in ATOMDB V 1.1.0 have been corrected.

Ni Wavelength Comparisons at 4 Temperatures

Nickel comparison a
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Nickel comparison b
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Nickel comparison c
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Nickel comparison d
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

b. The ATOMDB now replaces the CHIANTI calculations for Li-like ions with calculations by D. Liedahl using the HULLAC. The new calculations include levels up to principal quantum number n=7. The mutual lines are in good agreement. The new lines are apparent at the short wavelength end of the sample spectra.

Comparison of Lines from Li-like Ions

Neon Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Silicon
Line Comparison Plot
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 


Last modified:16 July 2002


The Chandra X-Ray Center (CXC) is operated for NASA by the Smithsonian Astrophysical Observatory.
60 Garden Street, Cambridge, MA 02138 USA.    Email: cxcweb@head.cfa.harvard.edu
Smithsonian Institution, Copyright © 1998-2004. All rights reserved.