AGN Model

class lightning.agn.AGNModel

Bases: BaseEmissionModel

An implementation of the Stalevski (2016) SKIRTOR models.

The broken power law accretion disk model is the SKIRTOR default. Future versions may include an option to swap in the Schartman+(2005) power law model as in X-Cigale. Optionally includes the X-Cigale recipe for polar dust extinction.

Parameters:

filter_labelslist, str

List of filter labels.

redshiftfloat

Redshift of the model.

wave_gridnp.ndarray, (Nwave,), float, optional

If set, the spectra are interpreted to this wavelength grid.

polar_dustbool

If True, AGN polar dust extinction and re-emission are implemented following the recipe from X-Cigale Note that even if this keyword is set to False, the polar dust optical depth remains a parameter of the model - it just doesn’t do anything, and should held at 0 in any fitting. (Default: True)

Attributes:

Lnu_restnumpy.ndarray, (29, 7, 132), float

High-res spectral model grid. First dimension covers U, the second q_PAH, and the third covers wavelength.

Lnu_obsnumpy.ndarray, (29, 7, 132), float

(1 + redshift) * Lnu_rest

mean_Lnunumpy.ndarray, (29, 7, Nfilters), float

The Lnu_obs grid integrated against the filters.

Lbolnumpy.ndarray, (29, 7), float

Total luminosity of each model in the grid.

wave_grid_restnumpy.ndarray, (132,), float

Rest-frame wavelength grid for the models.

wave_grid_obs

nu_grid_rest

nu_grid_obs

Methods

get_model_lnu(params[, exptau])	Produce the AGN model as observed in the given filters.
get_model_lnu_hires(params[, exptau])	Produce the high-resolution AGN spectral model.
print_params([verbose])	If verbose, print a nicely formatted table of the models, their parameters, and the description of the parameters.

References

• https://ui.adsabs.harvard.edu/abs/2016MNRAS.458.2288S/abstract

• https://sites.google.com/site/skirtorus

get_model_lnu(params, exptau=None)

Produce the AGN model as observed in the given filters.

Given a set of parameters, produce the observed AGN SED, optionally attenuating it with the ISM dust attenuation model and a polar dust model, a là CIGALE.

Parameters:

paramsnp.ndarray, (Nmodels, 3) or (3,) float

The AGN model parameters.

exptaunp.ndarray, (Nmodels, Nwave) or (Nwave,) float

The ISM attenuation curve.

Returns:

lnu_hiresnp.ndarray, (Nmodels, Nfilters) or (Nfilters,) float

The high resolution AGN models.

get_model_lnu_hires(params, exptau=None)

Produce the high-resolution AGN spectral model.

Given a set of parameters, produce the AGN spectrum, optionally attenuating it with the ISM dust attenuation model and a polar dust model, a là CIGALE.

Parameters:

paramsnp.ndarray, (Nmodels, 3) or (3,) float

The AGN model parameters.

exptaunp.ndarray, (Nmodels, Nwave) or (Nwave,) float

The ISM attenuation curve.

Returns:

lnu_hiresnp.ndarray, (Nmodels, Nwave) or (Nwave,) float

The high resolution AGN models.