xilconv: angle-dependent reflection from an ionized disk
This convolution model from Chris Done combines an ionized disk table
model from the XILLVER model of Garcia et al. (2013)
with the
Magdziarz & Zdziarski Compton reflection code. It is a modification of
the rfxconv model described in Kolehmainen, Done &
Diaz Trigo (2011)
which is a modification of the model first
described in Done & Gierlinski (2006).
The algorithm used is as follows.
- Determine the average power-law index of the input spectrum
between 2 and 10 keV. For this index and the other input parameters
interpolate on the table models to generate the reflected spectrum
from the ionized disk.
- Estimate the average power-law index of the reflected spectrum
over the range 12 - 14 keV.
- Iterate over the Compton reflection models changing the
cross-section at 10 keV until a match is found with the index
calculated in the previous step.
- Renormalize the reflection spectrum calculated in step 1 to
match the Compton reflection calculated in step 3 at 14 keV.
- Calculate the final reflection spectrum by using the
renormalized ionized disk spectrum below 14 keV and the Compton
reflection spectrum above 14 keV.
When using this model it is essential to extend the energy range over
which the model is calculated because photons at higher energies are
Compton down-scattered into the target energy range. The energy range
can be extended using the extend command. The upper limit on the
energies should be set above that for which the input spectrum has
significant flux. To speed up the model, calculation of the output
spectrum can be limited to energies below a given value by using xset
to define XILCONV_MAX_E (in units of keV). For instance, suppose that
the original data extends up to 100 keV. To accurately determine the
reflection it may be necessary to extend the energy range up to 500
keV. Now to avoid calculating the output spectrum between 100 and 500
keV use the command xset XILCONV_MAX_E 100.0.
The core of this model is a Greens' function integration with one
numerical integral performed for each model energy. The numerical
integration is done using an adaptive method which continues until
a given estimated fractional precision is reached. The precision can
be changed by setting XILCONV_PRECISION eg xset XILCONV_PRECISION 0.05.
The default precision is 0.01 (ie 1%).
To use different ionized disk table model files than those installed
change the directory searched for these files using xset XILCONV_DIR.
The model parameters are as follows.
par1 = |
the relative reflection normalization. If
is negative then only the reflected component is
returned. |
par2 = |
|
par3 = |
the iron abundance relative to Solar. All other
elements are assumed to have Solar abundance. |
par4 = |
cosine of the inclination angle (degrees). |
par5 = |
the ionization parameter used by the table models. |
par6 = |
the exponential cut-off energy (keV). |
HEASARC Home |
Observatories |
Archive |
Calibration |
Software |
Tools |
Students/Teachers/Public
Last modified: Friday, 23-Aug-2024 13:20:40 EDT
|