ISIS Spectra from Badenes' Models

Badenes test input: layer 110 of DDTcb model

The test input files from Badenes are given here - they describe a single layer of a model at several different times.

• DDTcb.02.2.00_layer110_200.dat
• DDTcb.02.2.00_layer110_430.dat
• DDTcb.02.2.00_layer110_1000.dat
• DDTcb.02.2.00_layer110_2000.dat

Create ISIS parameters from the layer information

This is done by reading the Badenes' data file into an IDL structure and then outputting a file for each layer with variables for input to ISIS modeling set:

• IDL> bad_name = 'DDTcb.02.2.00_layer110_2000'
• IDL> bd = read_badenes(bad_name)
• IDL> ilay = 1 ; starts at 0
• IDL> badenes_to_ions, bd, ilay, DKPC=10.0, BAD_NAME=bad_name, /OUTPUT

The output parameter files are:

DDTcb...200_L110.i

DDTcb...430'

DDTcb...1000'

DDTcb...2000'

The parameters in these files are:

• mttemp - temp in degrees K (of the electrons)
• mtels - the Z values of elements present in the layer
• mtabs - the abundance of the elements with respect to AG89 "solar" values.
• mtnorm - the "norm" of the emission in the usual "n_e n_H dV / D^2" sense with appropriate factors and constants. The distance assumed is list in the file header, 10 kpc was used here.
• NOTE: for both mtnorm and mtabs the value of "n_H" used is fictional, that is, it is defined as: n_H = n_e/1.20. See comments in the badenes_to_ions routine for specifics.
• allels, allion, allifs --- These arrays give the "ion fraction" (allifs) for the element-ion-s (allels, allion) that are present (up to 11 per element included.) These fractions are used to scale the contribution of each ion to the spectrum based on the ratio of this fraction to the fraction that is assumed for the CEI case spectrum.

Create ISIS spectrum for each layer

This is done in ISIS making use of the previous parameter files:

• [ isis> .load start.i % one time "startup" file... includes much more than is needed for this. ]
• isis>
• isis> variable bad_layer = 110;
• isis> variable bad_name = "DDTcb.02.2.00_layer110_200";
• isis> ()=evalfile("badenes_ions_calc.i");

The output spectrum files are:

DDTcb...200_L110spect.dat

DDTcb...430'

DDTcb...1000'

DDTcb...2000'

They contain two columns:

• wavelngth - in A ; the range and spacing can be changed in 'ions_calc.
• flux - in photons/s/cm^2/bin

Sum all layers together

In general if calculating the flux from all layers of a model then the various 'LNNNspect.dat files would be summed and put into a 'spectlam.dat file. For this single-layer case just copy the 'L110spect.dat files into 'spectlam.dat ... (In general the simple IDL routine badenes_assemble.pro is used... Like I said, it's not optimized!)

Convert to spectrum vs Energy; include instrument effects

Back in IDL read in the wavelength-flux file and convert it to Energy units, multiply by approximate EPIC MOS effective area and smooth to represent CCD resolution. Output plots and energy spectra files:

• IDL> bad_name = 'DDTcb.02.2.00_layer110_2000'
• IDL> badenes_espect, bad_name

The output energy-spectrum files are:

DDTcb...200_spectvse.dat

DDTcb...430'

DDTcb...1000'

DDTcb...2000'

The energy-spectrum files contain two columns:

• energy - in keV
• flux - in counts/s/keV

Plots of the wavelength spectrum and the smoothed energy spectrum are available:

DDTcb...200_spectplots.pdf

DDTcb...430'

DDTcb...1000'

DDTcb...2000'

Note that the high-resolution wavelength spectrum (top panel) has been converted to energy units and plotted at high resolution for reference with the smoothed energy version (bottom panel.)

Make a comparison plot...

The four spectra above are plotted together in: 'ddSpectra.pdf ,
(created by badenes_compare.pro). This can be compared with Badenes's plot: 'Spectra.pdf.
Hmmm... There are differences... but the Si complex around 2 keV looks very similar for the 2000 year case...