We investigate the relation between the detection of the $11.3\,\mu$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,\mu$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,\mu$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN.

We measure the 10 and $18\mu$m silicate features in a sample of 67 local ($z<0.1$) type 1 active galactic nuclei (AGN) with available {\it Spitzer} spectra dominated by non-stellar processes. We find that the $10\mu$m silicate feature peaks at $10.3^{+0.7}_{-0.9}\mu$m with a strength (Si$_{p}$ = ln f$_{p}$(spectrum)/f$_{p}$(continuum)) of $0.11^{+0.15}_{-0.36}$, while the $18\mu$m one peaks at $17.3^{+0.4}_{-0.7}\mu$m with a strength of $0.14^{+0.06}_{-0.06}$. We select from this sample sources with the strongest 10$\mu$m silicate strength ($\sigma_{Si_{10\mu m}}>0.28$, 10 objects). We carry out a detailed modeling of the IRS/{\it Spitzer} spectra by comparing several models that assume different geometries and dust composition: a smooth torus model, two clumpy torus models, a two-phase medium torus model, and a disk+outflow clumpy model. We find that the silicate features are well modeled by the clumpy model of Nenkova et al. 2008, and among all models those including outflows and complex dust composition are the best (Hoenig et al. 2017). We note that even in AGN-dominated galaxies it is usually necessary to add stellar contributions to reproduce the emission at the shortest wavelengths.

In order to understand the diversity of classes observed in active galactic nuclei (AGN), a geometrically and optically thick torus of gas and dust is required to obscure the central engine depending on the line of sight to the observer. We perform a simultaneous fitting of X-ray and mid-infrared (mid-IR) spectra to investigate if the same structure could produce both emissions and, if this the case, to obtain better constraints for the physical parameters of the torus. In this case, we take advantage of the fact that both emissions show important signatures of obscuration. We used the nearby type-2 active nucleus IC 5063 as a test object. This object is ideal because of the wealth of archival data including some high resolution data. It also has a relatively high AGN luminosity that dominates at both X-ray and mid-IR frequencies. We use high spectral resolution NuSTAR and IRS/Spitzer spectra. The AGN dusty models used several physically motivated models. We found that the combination of the smooth torus models at mid-IR by Fritz et al. (2006) and at X-rays by Balokovi\'c et al. (2018), with the viewing and half-opening angles linked to the same value, is the best choice to fit the spectra at both wavelengths. This allows us to determine all the parameters of its torus. This result suggests that the structure producing the continuum emission at mid-IR and the reflection component at X-ray is the same. Therefore, we prove that this technique can be used to infer the physical properties of the torus, at least when AGN dust dominates the mid-IR emission and the reflection component is significant at X-rays.

At distances from the active galaxy nucleus (AGN) where the ambient temperature falls below ~1500-1800 K, dust is able to survive. It is thus possible to have a large dusty structure present which surrounds the AGN. This is the first of two papers aiming at comparing six dusty torus models with available SEDs, namely Fritz et al. (2006), Nenkova et al. (2008B), Hoenig & Kishimoto (2010), Siebenmorgen et al. (2015), Stalevski et al. (2016), and Hoenig & Kishimoto (2017). In this first paper we use synthetic spectra to explore the discrimination between these models and under which circumstances they allow to restrict the torus parameters, while our second paper analyzes the best model to describe the mid-infrared spectroscopic data. We have produced synthetic spectra from current instruments: GTC/CanariCam and Spitzer /IRS and future JWST/MIRI and JWST/NIRSpec instruments. We find that for a reasonable brightness (F12um > 100mJy) we can actually distinguish among models except for the two pair of parent models. We show that these models can be distinguished based on the continuum slopes and the strength of the silicate features. Moreover, their parameters can be constrained within 15% of error, irrespective of the instrument used, for all the models but Hoenig & Kishimoto (2017). However, the parameter estimates are ruined when more than 50% of circumnuclear contributors are included. Therefore, future high spatial resolution spectra as those expected from JWST will provide enough coverage and spatial resolution to tackle this topic.

This is the second in a series of papers devoted to explore a set of six dusty models of active galactic nuclei (AGN) with available spectral energy distributions (SEDs). These models are the smooth torus by Fritz et al. (2006), the clumpy torus by Nenkova et al. (2008B), the clumpy torus by Hoenig & Kishimoto (2010), the two phase torus by Siebenmorgen et al. (2015), the two phase torus by Stalevski et al. (2016), and the wind model by Hoenig & Kishimoto (2017). The first paper explores discrimination among models and the parameter restriction using synthetic spectra (Gonzalez-Martin et al. 2019A). Here we perform spectral fitting of a sample of 110 AGN drawn from the Swift/BAT survey with Spitzer/IRS spectroscopic data. The aim is to explore which is the model that describes better the data and the resulting parameters. The clumpy wind-disk model by Hoenig & Kishimoto (2017) provides good fits for ~50% of the sample, and the clumpy torus model by Nenkova et al. (2008B) is good at describing ~30% of the objects. The wind-disk model by Hoenig & Kishimoto (2017) is better for reproducing the mid-infrared spectra of Type-1 Seyferts while Type-2 Seyferts are equally fitted by both models. Large residuals are found irrespective of the model used, indicating that the AGN dust continuum emission is more complex than predicted by the models or that the parameter space is not well sampled. We found that all the resulting parameters for our AGN sample are roughly constrained to 10-20% of the parameter space. The derived outer radius of the torus is smaller for the smooth torus by Fritz et al. (2006) and the two phase torus by Stalevski et al. (2016) than the one derived from the clumpy torus by (Nenkova et al. 2008B). Covering factors and line-of-sight viewing angles strongly depend on the model used. The total dust mass is the most robust derived quantity.

We study the detectability of the emission associated with the AGN dusty structure at sub-mm wavelengths using ALMA, in a theoretical and observational way. Theoretically, we use the Clumpy models from Nenkova et al. together with the mid-infrared to X-ray and the radio fundamental plane scaling relations. We find that it is more likely to detect bigger and denser dusty tori at the highest ALMA frequency (666 GHz/450 micron). We also find that with 1h at 353 GHz/850 micron and 10h at 666 GHz/450 micron we can detect, with a high detection limit, a 1 mJy torus (characteristic of bright AGN). This means, an object for which the unresolved SED at 12 micron has a flux ~1mJy. Observationally, we use four prototypical AGN: NGC\,1052 (low-luminosity AGN), NGC\,1068 (Type-2), NGC\,3516 (Type 1.5), and IZw1 (QSO), with radio, sub-millimeter, and mid-IR data available. All the mid-infrared spectra are best fit with the smooth model reported by Fritz et al. A power law and a single, or a composition of, synchrotron component/s reproduce the cm radio wavelengths. We combined and extrapolated both fits to compare the extrapolation of both torus and jet contributors at sub-mm wavelengths with data at these wavelengths. Our observational results are consistent with our theoretical results. The most promising candidate to detect the torus is the QSO IZw1 (therefore, highly accreting sources in general), although it cannot be resolved due to the distance of this source. We suggest that to explore the detection of a torus at sub-mm wavelengths, it is necessary to perform an SED analysis including radio data, with particular attention to the angular resolution.

We studied the circumnuclear MIR emission in a sample of 19 local active galactic nuclei (AGN) with high spatial resolution spectra using T-ReCS (Gemini) and CanariCam (GTC), together with IRS/Spitzer observations. We measured the flux and the equivalent width for the 11.3 micron PAH feature and the [SIV] line emission as a function of galactocentric distance. This allowed to study the star formation (SF) at sub-kpc scales from the nucleus for a large sample of nearby AGN. The [SIV] line emission could be tracing the AGN radiation field within a few thousand times the sublimation radius (R_sub), but it often peaks at distances greater than 1000 R_sub. One possibility is that the SF is contributing to the [SIV] total flux. We found an 11.3 micron PAH emission deficit within the inner few tens of parsecs from the AGN. This deficit might be due to the destruction of the molecules responsible for this feature or the lack of SF at these distances. We found a sensible agreement in the expected shift of the relation of the AGN bolometric luminosity and the SF rate. This indicates that numerical models attributing the link between AGN activity and host galaxy growth to mergers are in agreement with our data, for most inner galaxy parts.

Several authors have claimed that the less luminous active galactic nuclei (AGN) are not capable of sustaining the dusty torus structure. Thus, a gradual re-sizing of the torus is expected when the AGN luminosity decreases. Our aim is to confront mid-infrared observations of local AGN of different luminosities with this scenario. We decomposed about ~100 IRS/Spitzer spectra of LLAGN and powerful Seyferts in order to decontaminate the torus component from other contributors. We have used the affinity propagation (AP) method to cluster the data into five groups within the sample according to torus contribution to the 5-15 um range (Ctorus) and bolometric luminosity. The AP groups show a progressively higher torus contribution and an increase of the bolometric luminosity, from Group 1 (Ctorus~ 0% and logLbol ~ 41) and up to Group 5 (Ctorus ~80% and log(Lbol) ~44). We have fitted the average spectra of each of the AP groups to clumpy models. The torus is no longer present in Group 1, supporting the disappearance at low-luminosities. We were able to fit the average spectra for the torus component in Groups 3 (Ctorus~ 40% and log(Lbol)~ 42.6), 4 (Ctorus~ 60% and log(Lbol)~ 43.7), and 5 to Clumpy torus models. We did not find a good fitting to Clumpy torus models for Group 2 (Ctorus~ 18% and log(Lbol)~ 42). This might suggest a different configuration and/or composition of the clouds for Group 2, which is consistent with a different gas content seen in Groups 1, 2, and 3, according to the detections of H2 molecular lines. Groups 3, 4, and 5 show a trend to decrease of the width of the torus (which yields to a likely decrease of the geometrical covering factor), although we cannot confirm it with the present data. Finally, Groups 3, 4, and 5 show an increase on the outer radius of the torus for higher luminosities, consistent with a re-sizing of the torus according to the AGN luminosity.