Disentangling accretion disk and dust emissions in the infrared spectrum of type 1 AGN

We use a semi-empirical model to reproduce the 0.1-10um spectral energy distribution (SED) of a sample of 85 luminous quasars. In the model, the continuum emission from the accretion disk as well as the nebular lines are represented by a single empirical template (disk), where differences in the optical spectral index are reproduced by varying the amount of extinction. The near- and mid-infrared emission of the AGN-heated dust is modelled as the combination of two black-bodies (dust). The model fitting shows that the disk and dust components are remarkably uniform among individual quasars, with differences in the observed SED largely accounted for by varying levels of obscuration in the disk as well as differences in the relative luminosity of the disk and dust components. By combining the disk-subtracted SEDs of the 85 quasars, we generate a template for the 1-10um emission of the AGN-heated dust. Additionally, we use a sample of local Seyfert 1 galaxies with full spectroscopic coverage in the 0.37um to 39um range to demonstrate a method for stitching together spectral segments obtained with different PSF and extraction apertures. We show that the disk and dust templates obtained from luminous quasars also reproduce the optical-to-mid-infrared spectra of local Seyfert 1s when the contribution from the host galaxy is properly subtracted.


OBTAINING ACCRETION DISK AND DUST TEMPLATES FROM LUMINOUS QUASARS
We use a sample of 85 luminous quasars (νL ν [3µm]>10 45.5 erg s −1 ) selected for their spectroscopic coverage (AKARI and/or Spitzer/IRS) in the rest-frame 2.5-10 µm range. In addition to the AKARI and Spitzer spectroscopy we obtain optical photometry from the Sloan Digital Sky Survey (SDSS) Data Release 12, near-infrared (NIR) photometry from the Two Micron All Sky Survey (2MASS), the VISTA Hemisphere Survey (VHS), and the UKIRT Infrared Deep Sky Survey (UKIDSS), and mid-infrared photometry in 4 bands from WISE [see 1, for details].
In such luminous quasars the optical emission of the AGN easily outshines that of the host galaxy. The emission from dust is expected to be negligible at λ 0.85µm because the maximum temperature of dust grains is limited by sublimation to ∼1500 K [2]. As a consequence, only the accretion disk and the emission lines from the broad line region (BLR) and narrow line region (NLR) contribute significantly to the rest-frame 0.1-0.85µm spectrum.
In Hernán-Caballero et al. [1] we showed that a single empirical quasar template like that of Shen [3] suffices to model the rest-frame UV-optical (0.1-0.85µm) SED of luminous quasars if we allow for an adjustable extinction to reproduce the diversity in optical spectral indices among individual quasars. To extend the template into the NIR and MIR ranges, we assumed that the disk emission follows a powerlaw with the theoretical slope α=1/3 predicted for a locally heated optically thick disk [e.g. 4,5] and confirmed through polarized light observations [6]. We also added to the new template the NIR nebular lines extracted from the quasar template of Glikman et al. [7].
We fit the rest-frame 0.1-10µm SED of the 85 quasars with a two component disk+dust model. The disk component (plus nebular lines) is represented by the template described above, modified by an adjustable level of extinction with a wavelength dependency following the Small Magellanic Cloud Bar extinction law [8]. This law is often used to de-redden quasars [e.g. 9, 10] since it lacks the 2175Å absorption feature. The AGN-heated dust component is represented by the linear combination of two black-bodies (hot and warm) at adjustable temperatures within the intervals 850-2000 K and 150-900 K, respectively. Figure 1 shows examples of the best-fitting disk+dust decomposition model for a representative subsample of the 85 quasars. The residuals of the fit around ∼10% are consistent with the photometric uncertainties. We find a systematic excess emission in the 1.0-1.5µm range relative to the model. The median, mean, and standard deviation of the excess at restframe 1.2µm is 32%, 40% and 33%, respectively. This is comparable to the values found in Hernán-Caballero et al. [1] in spite of the stronger NIR continuum and inclusion of NIR nebular lines in the new disk template, supporting our previous claim that the excess originates in the dust component.
The distributions of A V for the extinction of the accretion disk and the relative luminosity of the dust and disk components (represented by the luminosity ratio between λ=3.0µm and λ=0.5µm) are shown in Figure 2. Negative values of A V are needed to model the disk of the bluest ∼25% of quasars because the disk template is an average of observed spectra that have not been de-reddened. Therefore the extinction values obtained in the fit are relative to the (unknown) average extinction of the sample in Shen [3]. Moderate values of A V (between -0.2 and 0.2 mag) fit all but the reddest 15% of quasars. The νL ν (3µm)/νL ν (0.5µm) ratio peaks at ∼1, indicating that the peak luminosities of the disk and dust components are typically similar, albeit the full range of variation being a whole order of magnitude. Figure 3 shows the infrared SEDs of the individual quasars, normalized at rest-frame 3µm, after subtraction of the disk component. The SEDs show little dispersion (<0.3 dex) between ∼1.5µm and ∼6µm, suggesting a largely invariable spectrum for the hot dust, that is consistent with black-body emission at a temperature close to the grain sublimation limit. The dispersion increases at longer wavelengths at least in part due to the onset of the broad silicate feature (which may be in emission or absorption). On the other hand, the larger dispersion at λ 1.5µm is a consequence of the uncertainty introduced by the subtraction of a disk component that is increasingly dominant at shorter wavelengths.
We obtain a template for the AGN-heated dust by averaging the AKARI+IRS spectra of the individual quasars. At λ<2µm only a few quasars have spectroscopic coverage, therefore we rely on the disksubtracted broadband photometry to extend the template to λ<2µm. We find that a 1400 K black-body is consistent with most of the broadband data-points.

APPLICATION TO LOCAL SEYFERT 1 GALAXIES
To test whether the templates for disk and dust emission obtained for luminous quasars are also representative of the nuclear emission in less luminous type 1 AGN, we have performed spectral decomposition of a sample of 13 local (z<0.07) Seyferts and quasars with rest-frame 3µm luminosities in the range 10 42.8−44.4 erg s −1 . The sample is a subset of the 23 broad-emission line AGN observed by Landt et al. [11]. We chose this sample because they obtained nearly simultaneous spectra in the optical (0.37-0.75µm) and NIR (0.8-2.4µm) with the FAST and SpeX spectrographs on the Tillinghast 1.5 m telescope and the NASA Infrared Telescope Facility, respectively. Out of their 23 sources, we select only the 13 that also have AKARI 2.5-5.0µm spectra from [12] and Spitzer/IRS 5.2-39µm spectra in the CASSIS [13,14] or ATLAS-IRS [15] databases. Therefore we have continuous spectroscopic coverage from ∼0.35µm to ∼35µm in the restframe. However, since our templates are not defined beyond 10µm, we perform the spectral decomposition only in the ∼0.35-10µm range.
Spectral decomposition on local type 1 AGN is more challenging compared to luminous quasars because the emission from the host galaxy is no longer negligible. The contribution from the stars to the total flux is important only at optical and NIR wavelengths, while the emission of the interstellar medium (ISM), in particular dust grains and aromatic hydrocarbons, is significant only in the mid-infrared. Since the integrated spectrum of the stars in a galaxy peaks at 1.6µm and its NIR shape does not vary substantially among the different spectral types, we can use a single template to represent the stellar emission of the host galaxy in all the sources. We choose the S0 galaxy template from Polletta et al. [16]. In the optical, the actual spectra of the host galaxies may diverge from the template at shorter wavelengths depending on the age and extinction of the stellar population, but the difference only becomes important at λ 0.4µm.
Another difficulty in decomposing local galaxies is that while the AGN emission is always spatially unresolved, the host galaxy is extended and its contribution to the observed spectrum depends on the spatial resolution of the observations and the size of the extraction aperture (see [17] for a discussion). The extraction aperture of the slit-less AKARI observations (60"×7.5") is wider than the IRS slit (3.6" for SL1 and SL2 modules), and both are larger than those of the FAST (3") and SpeX (0.8") observations. To compensate for this, we compute the best-fitting disk+dust+stellar model for the local AGN as follows: for each of the four spectral segments (FAST, SpeX, AKARI, IRS) we obtain the values of the parameters that minimize the residuals while fitting by least squares the model: where f disk (λ), f dust (λ), and f star (λ) are the templates and e −A V τ (λ) is the extinction correction for the disk template. The coefficients a, b, and A V must be the same for all four spectral segments, but coefficient c may take different values in each segment to compensate for the different amounts of host galaxy light in each aperture. We also apply a scaling factor between 0.7 and 1.3 to the fluxes in the FAST, AKARI and IRS segments to correct for any potential biases in the absolute flux calibration of up to ∼30% relative to the SpeX spectrum. The values of these scaling factors are treated as free parameters in the fitting algorithm and they are computed independently for each source.
We use the coefficients in the best fitting model to obtain a stitched spectrum that merges the FAST, SpeX, AKARI, and IRS segments. For this we first multiply each segment (except the SpeX one, which is taken as reference) by its corresponding scaling factor, and then subtract the excess stellar component given by: where the superscript i indicates the spectral segment. The resulting stitched spectrum corresponds to our prediction of the spectrum that would be observed if all the segments had the same PSF, were extracted in the same aperture, and had perfectly consistent absolute calibrations.
We evaluate the quality of the fits using the normalized root mean square error (RMSE), defined as: where f i and F (λ i ) are the flux densities at wavelength λ i in the stitched spectrum and the model, respectively. This statistic represents the typical relative residual between the stitched spectrum and the model, and is more informative than χ 2 in situations where residuals are not dominated by noise but differences between the model and the intrinsic spectrum of the source [17,18]. Figure 4 shows the original and stitched spectra for the 13 sources in the sample, as well as the best fitting decomposition into disk, dust, and stellar components. The contribution from the stellar component to the total emission varies from negligible (e.g. Mrk 79, Mrk 509, and Mrk 335) to dominant (e.g. Mrk 590), while the interstellar emission is negligible in all but NGC 7469. There is no obvious correlation between the spectral class (Seyfert types 1.0, 1.2, and 1.5) and model parameters such as A V of the disk or the relative luminosities of the components, albeit the sample size is too small to draw any strong conclusions. The fits are remarkably good for most sources (typical RMSE ∼10%), with significant discrepancies arising only at very short (λ<0.4µm) and long (λ>7-8µm) wavelengths. This suggests that the optical and NIR spectrum of the disk and dust emissions are relatively uniform among type 1 AGN regardless of luminosity. The discrepancies at short wavelength may be caused by uncertainty in the extinction affecting the disk component due to lack of restframe-UV observations, as well as increasing discrepancy between the assumed template for the stellar population and the actual spectrum of the host at short wavelengths. Discrepancies at λ>7-8µm are caused by source-to-source variation in the strength of the mid-infrared aromatic features relative to the stellar continuum, and in some cases like Mrk 590, NGC 4151, and PG 0844+349, an unusually strong silicate emission feature at ∼10µm. Hatziminaoglou et al. [19] showed that accurate modelling of the mid-infrared spectrum of Seyfert 1s requires a different approach, splitting the stellar and ISM emissions into separate spectral components and using multiple dust templates with different spectral indices and silicate strengths.    The grey lines represent the (AKARI+)IRS spectrum, while the red dots represent the broadband data points at wavelengths with no spectroscopic coverage. All the individual SEDs are normalized at restframe 3µm. The solid black line is the composite AKARI+IRS spectrum obtained for the same sample in Hernán-Caballero et al. [1]. The dashed black line is an extrapolation of this composite using a T=1400 K black-body spectrum.