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 examine star formation activity in a distance- (z<0.1) and flux-limited sample of quasars (QSOs). Mid-infrared (MIR) spectral diagnostics at high spatial resolution (~0.4 arcsec) yield star formation rates (SFRs) in the inner regions (~300 pc to 1 kpc) for 13 of 20 of the sample members. We group these objects according to the size probed by the high angular resolution spectroscopy, with characteristic scales of <0.7 and ~0.7-1 kpc. Using the polycyclic aromatic hydrocarbon (PAH) feature at 11.3 um, we measure SFRs around 0.2 and 1.6 M_{\odot}yr^{-1}. We also measure the larger aperture PAH-derived SFRs in the individual IRS/Spitzer spectra of the sample and obtain a clear detection in ~58 percent of them. We compare smaller and larger aperture measurements and find that they are similar, suggesting that star formation activity in these QSOs is more centrally concentrated, with the inner region (>~1 kpc) accounting for the majority of star formation measured on these scales, and that PAH molecules can be present in most local MIR-bright QSOs within a few hundred pc from the central engine. By comparison with merger simulations, we find that our estimation of the SFR and black hole (BH) accretion rates are consistent with a scenario in which the star formation activity is centrally peaked as predicted by simulations.

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.

We present mid-infrared (MIR, 7.5-13.5 $\mu$m) imaging and spectroscopy observations obtained with the CanariCam (CC) instrument on the 10.4m Gran Telescopio CANARIAS for a sample of 20 nearby, MIR bright and X-ray luminous QSOs. We find that for the majority of QSOs the MIR emission is unresolved at angular scales nearly 0.3 arcsec, corresponding to physical scales $<600$ pc. We find that the higher-spatial resolution CC spectra have similar shapes to those obtained with Spitzer/IRS, and hence we can assume that the spectra are not heavily contaminated by extended emission in the host galaxy. We thus take advantage of the higher signal to noise Spitzer/IRS spectra, as a fair representation of the nuclear emission, to decompose it into a combination of active galactic nuclei (AGN), polycyclic aromatic hydrocarbon (PAH) and stellar components. In most cases the AGN is the dominant component, with a median contribution of 85 per cent of the continuum light at MIR (5-15 $\mu$m) within the IRS slit. This IR AGN emission is well reproduced by clumpy torus models. We find evidence for significant differences in the parameters that describe the dusty tori of QSOs when compared with the same parameters of Seyfert 1 and 2 nuclei. In particular, we find a lower number of clouds ($N_{0}<12$), steeper radial distribution of clouds ($q=1.5-3.0$), and clouds that are less optically thick ($\tau_{V}<100$) than in Seyfert 1, which could be attributed to dusty structures that have been partially evaporated and piled up by the higher radiation field in QSOs. We find that the combination of the angular width $\sigma_{torus}$, viewing angle $i$, and number of clouds along the equatorial line $N_{0}$, produces large escape probabilities ($P_{esc} > 2$ per cent) and low geometrical covering factors ($f_{2}<0.6$), as expected for AGN with broad lines in their optical spectra.

We present an analysis of the nuclear infrared (IR, 1.6 to 18 $\mu$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($\lambda_{C}=8.7\, \mu$m) and $7.5-13$ $\mu$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, \mu$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature.

We have carried out extensive and detailed photoionization modeling to successfully constrain the locations of different emission-line galaxies in optical and mid-infrared diagnostic diagrams. Our model grids cover a wide range in parameter space for the active galaxy continuum and starburst galaxies with different stellar population laws and metallicities. We compare the predicted AGN and star-formation mid-infrared line ratios [Ne III]15.56mm/[Ne II]12.81mm and [O IV]25.89mm/[Ne III]15.56mm to the observed values, and find that the best fit for the AGN is via a two-zone approximation. This two-zone approximation is a combination of a matter-bounded component, where [Ne III] and [O IV] are emitted efficiently, and a radiation-bounded component that maximizes [Ne II] emission. We overlay the predictions from this two-zone approximation onto the optical [O III]l5007/Hbeta and [N II]l6583/Halpha diagnostic diagram derived from the Sloan Digital Sky Survey, to find that the high-density and low-ionization radiation-bounded component in our two-zone AGN approximation model provides a good lower limit for [N II] emission. This establishes a new theoretical demarcation line for the minimum AGN contribution in this diagram. This new classification results by a factor of ~1.4 in a higher AGN population than predictions derived from previous divisions of star-forming galaxies. Similarly, we define a maximum AGN contribution in the [O III]/Hbeta and [N II]/Halpha diagram by using a two-zone approximation within a parameter range typical of the narrow-line region.

We present MIR sub-arcsec observations of Mrk 1383, a nearby QSO, part of a sample of nearby QSOs currently being observed with the 10.4 m Gran Telescopio de Canarias (GTC). The radial profile of the QSO is analyzed and we find that most of the MIR emission is unresolved (\theta ~ 0.4", 600 pc) and consistent with previous measurements from Spitzer (\theta ~ 4"). We derive the spectral energy distribution (SED) of the nuclear emission, relatively uncontaminated by starlight, combining measurements from the literature as a first step towards characterizing the torus parameters in this QSO.