Color-temperature-corrected disc and energetically coupled Comptonization model for AGN. XSPEC model.
AGN spectral energy distributions are complex, but can be phenomenologically fit by a disc, optically thick, low temperature thermal Comptonization (to produce the soft X-ray excess) and an optically thin, high temperature themal Comptonization (to produce the power law emission which dominates above 2 keV). Here we combine these three components together assuming that they are all ultimately powered by gravitational energy released in accretion. We assume that the gravitational energy released in the disc at each radius is emitted as a (color-temperature-corrected) blackbody only down to a given radius, R_corona. Below this radius, we further assume that the energy can no longer completely thermalise, and is distributed between powering the soft excess component and the high energy tail. This imposes an important energetic self consistency on the model. The key aspect of this model is that the optical luminosity constrains the mass accretion rate through the outer disc, M_dot, provided there is an independent estimate of the black hole mass (from, e.g., the H-beta emission line profile). The total luminosity available to power the entire SED is [L_tot=eff M_dot c^2], where the efficiency is set by black hole spin assuming Novikov-Thorne emissivity.
There are two versions of the model: xsoptxagnf is the one recommended for most purposes, and has the color-temperature correction calculated for each temperature from the approximations given in Done et al. (2012); xsoptxagn instead allows the user to define their own color-temperature correction, fcol, which is then applied to annuli with effective temperature > T_scatt. In both models the flux is set by the physical parameters of mass, L/L_Edd, spin and distance, hence the model normalizations MUST be frozen at unity.
This is an additive model component.
|1||mass||Black hole mass in solar masses|
|2||dist||Luminosity distance in Mpc|
|4||astar||Dimensionless black hole spin|
|5||rcor||Coronal radius in R_g=GM/c^2 marking the transition from (color temperature corrected) blackbody emission to a Comptonised spectrum. If this parameter is negative then only the blackbody component is used.|
|6||logrout||Log of the outer radius of the disc in units of R_g; if this is '-ve' the code will use the self gravity radius as calculated from Laor & Netzer (1989)|
|7||kT_E||Electron temperature for the soft Comptonization component (soft excess) in keV|
|8||tau||Optical depth of the soft Comptonization component. If this parameter is negative then only the soft Compton component is used.|
|9||Gamma||Spectral index of the hard Comptonization component ('power law') which has temperature fixed to 100 keV. If this parameter is negative then only the hard Compton component is used.|
|10||fpl||Fraction of the power below rcor which is emitted in the hard Comptonization component|
|14||norm||Normalization, must be frozen|
This information is taken from the XSPEC User's Guide. Version 12.9.1n of the XSPEC models is supplied with CIAO 4.10.
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, 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, xsvnei, 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