ISMabs: A high resolution ISM absorption model with variable columns for individual ions.ISMabs is an X-ray photoabsorption model for the interstellar medium that takes into account both neutral and ionized species from H, He, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni and Zn. This model has been developed by Efrain Gatuzz in collaboration with Javier Garcia, Tim Kallman, Claudio Mendoza, and Tom Gorczyca. A complete description of the science behind the model is described in Gatuzz et al. (2014b).
The contents of the tarfile include:
atomic_data/AtomicData.fits -- atomic database binary .fits file. This must reside in the directory atomic_data inside the folder where the model is located.
ismabs.f90 -- source code for ISMabs
lmodel_ismabs.dat -- local model definition file needed by xspec.
compile.sh -- installation script written on bash.
InstallationYou can use the compile.sh file to install the model by doing
sh compile.shOr you can setting up and using this model as described in the xspec manual:
0) You need to have the heasoft package installed on your machine, but it must be built from source. Local models cannot be installed from a binary installation.
1) Untar the tarfile somewhere in your user area.
2) Setup your headas environment (eg. 'setenv HEADAS /path/to/architecture', and 'source \$HEADAS/headas-init.csh').
'initpackage ismabs lmodel.dat path-to-model-directory',where path-to-model-directory is the full path of the model directory.
After the build is complete type
'lmod ismabs path-to-current-directory'In subsequent sessions you don't need to do the initpackage step again, just the lmod.
ParametersThe input parameters include the column densities for the neutral elements as well as the single and double ionized ions. The redshift is also a model parameter.
Note that the He I column density is not included as a free parameter in the model, due to an inherent degeneracy between the relative columns of H, He I, He II. Instead, we impose the following constraint"
He I = 0.1 * H
Atomic DataWith the default set up - that is, if you have run compile.sh, the model will look for the cross-section data file in atomic_data/AtomicData.fits, relative to the directory in which the module is located.
The XSPEC xset command can be used to set the ISMABSROOT variable; if this is set then it is used instead of the path to the module. So after
xset ISMABSROOT /data/ismabs/then the model will use the file /data/ismabs/atomic_data/AtomicData.fits (the ISMABSROOT refers to the directory containing the atomic_data/ directory). Note that ISMABSROOT over-rides any changes made by running compile.sh when building the model.
The location of the file can be found by setting the XSPEC chatter level to 20 or higher - e.g.
before evaluating the model.
We include in our model the following cross sections:
Compared with previous photoabsorption models, which solely rely on neutral species, the inclusion of ions leads to improvements of the spectral fits. Fit parameters comprise the column densities of abundant contributors that allow direct estimates of ionization states. The figure below shows the theoretical photoabsorption cross sections for neutral, single, and double ionized Neon. The vertical lines indicate where the absorption features are observed in the astrophysical observations. The left panel displays the original cross sections, while the right panel shows them after wavelength corrections have been applied (for more details see Gatuzz et al. 2014b).
HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public
Last modified: Thursday, 02-Jul-2015 15:55:05 EDT