Calculate a simple table model (grid18)
The xstar webpage contains links to pre-calculated table models that can be downloaded and directly used in xspec. While these tables are good for quicklook and preliminary analyses we generally recommend users recalculate tables for their specific data analysis need to make sure they have full control and understanding of the parameter space and guarantee the use of the most recent version of xstar. In this example we recalculate the table model ‘grid18’ using xstar2xspec interactively.
After initializting HEAsoft, in an empty directory, we create a copy
of all required parameter files and reset $PFILES environment
varriable. This step is optional but highly recommended.
$ cp $HEADAS/syspfiles/xstar.par .
$ cp $HEADAS/syspfiles/xstinitable.par .
$ cp $HEADAS/syspfiles/xstar2table.par .
$ export PFILES=`pwd`
Now run xstar2xspec from the command line and with the following prompted parameters.
$ xstar2xspec
rm: xstinitable.lis: No such file or directory
rm: xstinitable.fits: No such file or directory
rm: xstar2xspec.log: No such file or directory
XSTAR2XSPEC: Initialize the Multiple XSTAR run...
spectrum type? [pow]
radiation temperature or alpha soft maximum? [-1.2] -1
covering fraction soft maximum (0.:1.) [1.]
density soft maximum (cm**-3) (1.e+4:1.E21) [1.E+12]
luminosity soft maximum (/10**38 erg/s) (1.e-20:1.E10) [1.e+6]
column density soft maximum (cm**-2) (1.e+10:1.E25) [1.e+24] 1e23
column density interpolation type (0:1) [1]
column density soft minimum (1.:1.E25) [1.E+21] 1e19
column density number of steps (1:20) [7] 5
log(ionization parameter) soft maximum (erg cm/s) (-10.:+10.) [5.] 4
log(ionization parameter) interpolation type (0:1) [0]
log(ionization parameter) soft minimum (-10.0:+10.0) [1.] -4
log(ionization parameter) number of steps (1:40) [9] 17
model name [template] XSTAR_fine_grid_1-3
abundance table [xdef]
XSTAR2XSPEC: Entering XSTAR2XSPEC Main Processing Loop...
XSTAR2XSPEC: Executing Pass #1
xstar returns: xstar returns: xstar version 2.59cj
pass number= 1 -1
log(r) delr/r log(N) log(xi) x_e log(n) log(t) h-c(%) h-c(%) log(tau)
fwd rev
18.00 -36.00 -10.00 -4.00 0.37 12.00 3.93 -0.24 0.00 -10.00 -10.00 78
18.00 -36.00 -10.00 -4.00 0.37 12.00 3.93 0.22 0.88 -10.00 -10.00 32
18.00 -12.92 17.08 -4.00 0.31 12.00 3.62 8.38 1.00 -0.30 -10.00 71
[...]
XSTAR2XSPEC: Executing Pass #85
xstar returns: xstar returns: xstar version 2.59cj
pass number= 1 -1
log(r) delr/r log(N) log(xi) x_e log(n) log(t) h-c(%) h-c(%) log(tau)
fwd rev
14.00 -36.00 -10.00 4.00 1.21 12.00 7.59 -0.00 0.00 -10.00 -10.00 26
14.00 -36.00 -10.00 4.00 1.21 12.00 7.59 -0.00 -1.21 -10.00 -10.00 1
14.00 -3.00 23.00 4.00 1.21 12.00 7.59 -0.00 -1.21 -5.28 -10.00 1
final print: 0
xstar: Prepping to write spectral data
xstar: Done writing spectral data
total time 42.034025847911835
xstar2table returns: xstar2table v1.0
Compiled: Mar 15 2024 13:12:37
1) column : Value is interpolated.
Logarithmically sampled at 5 points:
1e+19 1e+20 1e+21 1e+22 1e+23
2) rlogxi : Value is interpolated.
Linearly sampled at 17 points:
-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4
Model Control Parameters:
Spectrum Type: pow
Spectrum File: spect.dat
Spectrum Units: Energy
Abundance Table:
Is redshift a parameter?: Yes
Number of Steps: 3
Number of Iterations: 99
Write Switch: Yes
Print Switch: Yes
Step Size Choice Switch: 0
Number of Passes: 1
Constant Pressure Switch: No
Courant multiplier (emult) 0.50
taumax for Courant step: 5.00
Minimum electron fraction: 0.10
Critical ion abundance: 0.00
Turbulent velocity (km/s): 300.00
radius exponent: 0.00
number of energy bins: 9999
Model Name: XSTAR_fine_grid_1-3
Energy Range: 100.00 - 20000.00 eV
SliceEnergySpectra: Selected 3414 bins in range ( 99.891- 19939) eV with indices (4452-7866).
xstar2table: Install a check of the ENERGIES extension here?
xstar2table: Does it match the previous ENERGIES extension?
xstar2table: Successful Completion
xstar2table returns: xstar2table v1.0
Compiled: Mar 15 2024 13:12:37
1) column : Value is interpolated.
Logarithmically sampled at 5 points:
1e+19 1e+20 1e+21 1e+22 1e+23
2) rlogxi : Value is interpolated.
Linearly sampled at 17 points:
-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4
Model Control Parameters:
Spectrum Type: pow
Spectrum File: spect.dat
Spectrum Units: Energy
Abundance Table:
Is redshift a parameter?: Yes
Number of Steps: 3
Number of Iterations: 99
Write Switch: Yes
Print Switch: Yes
Step Size Choice Switch: 0
Number of Passes: 1
Constant Pressure Switch: No
Courant multiplier (emult) 0.50
taumax for Courant step: 5.00
Minimum electron fraction: 0.10
Critical ion abundance: 0.00
Turbulent velocity (km/s): 300.00
radius exponent: 0.00
number of energy bins: 9999
Model Name: XSTAR_fine_grid_1-3
Energy Range: 100.00 - 20000.00 eV
SliceEnergySpectra: Selected 3414 bins in range ( 99.891- 19939) eV with indices (4452-7866).
xstar2table: Install a check of the ENERGIES extension here?
xstar2table: Does it match the previous ENERGIES extension?
xstar2table: Successful Completion
XSTAR2XSPEC: 0 runs remaining to process.
rm: xout_detail.fits: No such file or directory
rm: xout_step.lis: No such file or directory
XSTAR2XSPEC: Execution time: 63464.00 seconds
XSTAR2XSPEC: Processing Complete!