Non-equilibrium ionization collisional plasma model with variable abundances. XSPEC model.
Non-equilibrium ionization collisional plasma model. This assumes a constant temperature and single ionization parameter. It provides a characterisation of the spectrum but is not a physical model. The references for this model can be found in the help file for the xsequil model.
See also the xsnei and xsvvnei models.
Several versions are available. To switch between them use the set_xsxset command: set_xsxset('NEIVERS', '1.0') gives the version from XSPEC v11.1; set_xsxset('NEIVERS', '1.1') uses updated calculations of ionization fractions using dielectronic recombination rates from Mazzotta et al (1998); set_xsxset('NEIVERS', '2.0') uses the same ionization fractions as 1.1 but uses APED to calculate the resulting spectrum; and set_xsxset('NEIVERS', '3.0') uses the beta release from AtomDB 3.0. Note that versions 1.x have no emission from Ar. The default is version 3.0.
This is an additive model component.
|1||kT||plasma temperature (keV)|
|2||H||hydrogen density in cm^-3|
|3-14||(element)||Abundances for He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni with respect to Solar, as defined by the set_xsabund command.|
|15||Tau||ionization timescale (s/cm^3)|
|17||norm||10^-14 / (4 pi (D_A*(1+z))^2) Int n_e n_H dV, where D_A is the angular diameter distance to the source (cm), n_e is the electron density (cm^-3), and n_H is the Hydrogen density (cm^-3)|
This information is taken from the XSPEC User's Guide. Version 12.9.1n of the XSPEC models is supplied with CIAO 4.10.
Using NEIVERS 2.0
The data files needed when version 2.0 of NEIVERS are not provided with Sherpa (or X-Spec) and must be downloaded separately: XSPEC data for APEC NEI Models (19Mb). Once unpacked, it will create the directory heasoft-6.22.1/spectral/modelData/ and the two files can either be moved to $ASCDS_INSTALL/ots/spectral/modelData/ or the NEIAPECROOT variable set to
sherpa> set_xsxset('NEIAPECROOT', '/path/to/heasoft-6.22.1/spectral/modelData/APEC_nei_v11')
where "/path/to" should be replaced by the location where the tar file was unpacked.
Checking it worked
To check that this works, try:
sherpa> set_xschatter(25) sherpa> set_xsxset('NEIVERS', '2.0') sherpa> create_model_component('xsnei', 'mdl') sherpa> mdl([0.1, 0.2, 0.3, 0.4, 0.5])
which will evaluate the xsnei model over the bins 0.1-0.2, 0.2-0.3, 0.3-0.4, and 0.4-0.5 keV, and - as the X-Spec chatter level has been increased to 25 - will also print messages to the screen showing what files were used or if there was an error loading them.
For a list of known bugs and issues with the XSPEC models, please visit the XSPEC bugs page.
To check the X-Spec version used by Sherpa, use the get_xsversion routine from the xspec module:
sherpa> from sherpa.astro.xspec import get_xsversion sherpa> get_xsversion() '12.9.1n'
- absorptionedge, absorptiongaussian, absorptionlorentz, absorptionvoigt, accretiondisk, atten, bbody, bbodyfreq, beta1d, beta2d, blackbody, box1d, box2d, bpl1d, bremsstrahlung, brokenpowerlaw, ccm, const1d, const2d, cos, delta1d, delta2d, dered, devaucouleurs2d, disk2d, edge, emissiongaussian, emissionlorentz, emissionvoigt, erf, erfc, exp, exp10, fm, gauss1d, gauss2d, hubblereynolds, jdpileup, linebroad, list_model_components, list_models, lmc, load_xscflux, load_xsgsmooth, load_xsireflect, load_xskdblur, load_xskdblur2, load_xskerrconv, load_xslsmooth, load_xspartcov, load_xsrdblur, load_xsreflect, load_xssimpl, load_xszashift, load_xszmshift, log, log10, logabsorption, logemission, logparabola, lorentz1d, lorentz2d, models, normbeta1d, normgauss1d, normgauss2d, opticalgaussian, poisson, polynom1d, polynom2d, polynomial, powerlaw, powlaw1d, recombination, scale1d, scale2d, schechter, seaton, sersic2d, shell2d, sigmagauss2d, sin, sm, smc, sqrt, stephi1d, steplo1d, tablemodel, tan, xgal, xs, xsabsori, xsacisabs, xsagauss, xsapec, xsbapec, xsbbody, xsbbodyrad, xsbexrav, xsbexriv, xsbkn2pow, xsbknpower, xsbmc, xsbremss, xsbvapec, xsbvvapec, xsc6mekl, xsc6pmekl, xsc6pvmkl, xsc6vmekl, xscabs, xscemekl, xscevmkl, xscflow, xscompbb, xscompls, xscompmag, xscompps, xscompst, xscomptb, xscompth, xscomptt, xsconstant, xsconvolve, xscplinear, xscutoffpl, xscyclabs, xsdisk, xsdiskbb, xsdiskir, xsdiskline, xsdiskm, xsdisko, xsdiskpbb, xsdiskpn, xsdust, xsedge, xseplogpar, xseqpair, xseqtherm, xsequil, xsexpabs, xsexpdec, xsexpfac, xsezdiskbb, xsgabs, xsgadem, xsgaussian, xsgnei, xsgrad, xsgrbm, xsheilin, xshighecut, xshrefl, xskerrbb, xskerrd, xskerrdisk, xslaor, xslaor2, xslogpar, xslorentz, xslyman, xsmeka, xsmekal, xsmkcflow, xsnei, xsnotch, xsnpshock, xsnsa, xsnsagrav, xsnsatmos, xsnsmax, xsnsmaxg, xsnsx, xsnteea, xsnthcomp, xsoptxagn, xsoptxagnf, xspcfabs, xspegpwrlw, xspexmon, xspexrav, xspexriv, xsphabs, xsplabs, xsplcabs, xsposm, xspowerlaw, xspshock, xspwab, xsraymond, xsredden, xsredge, xsrefsch, xsrnei, xssedov, xssirf, xssmedge, xsspexpcut, xsspline, xssrcut, xssresc, xssss_ice, xsstep, xsswind1, xstbabs, xstbgrain, xstbvarabs, xsuvred, xsvapec, xsvarabs, xsvbremss, xsvequil, xsvgadem, xsvgnei, xsvmcflow, xsvmeka, xsvmekal, xsvnpshock, xsvphabs, xsvpshock, xsvraymond, xsvrnei, xsvsedov, xsvvapec, xsvvgnei, xsvvnei, xsvvnpshock, xsvvpshock, xsvvrnei, xsvvsedov, xswabs, xswndabs, xsxion, xszagauss, xszbabs, xszbbody, xszbremss, xszdust, xszedge, xszgauss, xszhighect, xszigm, xszpcfabs, xszphabs, xszpowerlw, xszredden, xszsmdust, xsztbabs, xszvarabs, xszvfeabs, xszvphabs, xszwabs, xszwndabs, xszxipcf