Search

We analyze JWST MIRI/MRS IFU observations of three Seyferts and showcase the intriguing polycyclic aromatic hydrocarbon (PAH) emission characteristics in regions of $\sim 500\,\rm pc$ scales over or around their active galactic nuclei (AGN). Combining the model predictions and the measurements of PAH features and other infrared emission lines, we find that the central regions containing a high fraction of neutral PAHs with small sizes, e.g., those in ESO137-G034, are in highly heated environments, due to collisional shock heating, with hard and moderately intense radiation fields. Such environments are proposed to be associated with inhibited growth or preferential erosion of PAHs, decreasing the average PAH size and the overall abundance of PAHs. We additionally find that the central regions containing a high fraction of ionized PAHs with large sizes, e.g., those in MCG-05-23-016, are likely experiencing severe photo-ionization because of the radiative effects from the radiative shock precursor besides the AGN. The severe photo-ionization can contribute to the ionization of all PAHs and further destruction of small PAHs. Overall, different Seyferts, even different regions in the same galaxy, e.g., those in NGC\,3081, can contain PAH populations of different properties. Specifically, Seyferts that exhibit similar PAH characteristics to ESO137-G034 and MCG-05-23-016 also tend to have similar emission line properties to them, suggesting that the explanations for PAH characteristics of ESO137-G034 and MCG-05-23-016 may also apply generally. These results have promising application in the era of JWST, especially in diagnosing different (i.e., radiative, and kinetic) AGN feedback modes., Comment: 20 pages (6 pages in the appendix), 6 figures in the main text, ApJL resubmitted after addressing referee's comments

Utilizing JWST MIRI/MRS IFU observations of the kiloparsec scale central regions, we showcase the diversity of ionized gas distributions and kinematics in six nearby Seyfert galaxies included in the GATOS survey. Specifically, we present spatially resolved flux distribution and velocity field maps of six ionized emission lines covering a large range of ionization potentials ($15.8-97.1$ eV). Based on these maps, we showcase the evidence of ionized gas outflows in the six targets, and find some highly disturbed regions in NGC\,5728, NGC\,5506, and ESO137-G034. We propose AGN-driven radio jets plausibly play an important role in triggering these highly disturbed regions. With the outflow rates estimated based on [Ne~{\footnotesize V}] emission, we find the six targets tend to have ionized outflow rates converged to a narrower range than previous finding. These results have important implication for the outflow properties in AGN of comparable luminosity., Comment: 34 pages (11 pages in the appendix), 18 figures in the main text, ApJ in press (accepted on July 26th)

We present arcsecond-scale observations of the active galactic nuclei (AGNs) of seven nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the 31.5 and 37.1 um filters of the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We isolate unresolved emission from the torus and find extended diffuse emission in six 37.1 um residual images in our sample. Using Spitzer/IRS spectra, we determine the dominant mid-infrared (MIR) extended emission source and attribute it to dust in the narrow line region (NLR) or star formation. We compare the optical NLR and radio jet axes to the extended 37.1 um emission and find coincident axes for three sources. For those AGNs with extended emission coincident with the optical axis, we find that spatial scales of the residual images are consistent with 0.1 - 1 kpc scale distances to which dust can be heated by the AGN. Using previously published subarcsecond 1 - 20 um imaging and spectroscopic data along with our new observations, we construct broadband spectral energy distributions (SEDs) of the AGNs at wavelengths 1 - 40 um. We find that three AGNs in our sample tentatively show a turnover in the SED between 30 - 40 um. Using results from Clumpy torus models and the Bayesian inference tool BayesClumpy, we find that the posterior outputs for AGNs with MIR turnover revealed by SOFIA/FORCAST have smaller uncertainties than AGNs that do not show a turnover., Comment: 16 pages, 9 figures

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission., Comment: 21 pages, 15 figures, accepted for publication in ApJ. The full list of Table 1 is available at http://www.kusastro.kyoto-u.ac.jp/~ichikawa/Table1_MR_20181107.txt

We present observations of NGC 1068 covering the $19.7-53.0$ $\mu$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $\mu$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $\mu$m within the central 700 pc. We compute $2.2-432$ $\mu$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $\mu$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $\mu$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $\mu$m ALMA observations. This result indicates that the $1-20$ $\mu$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $\mu$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $\mu$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $\mu$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $\mu$m can produce comparable results to those observed using IR interferometry., Comment: 12 pages, 7 Figures, Accepted for publication in ApJ

We present 31.5 micron imaging photometry of 11 nearby Seyfert galaxies observed from the Stratospheric Observatory For Infrared Astronomy (SOFIA) using the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST). We tentatively detect extended 31 micron emission for the first time in our sample. In combination with this new data set, subarcsecond resolution 1-18 micron imaging and 7.5-13 micron spectroscopic observations were used to compute the nuclear spectral energy distribution (SED) of each galaxy. We found that the turnover of the torus emission does not occur at wavelengths <31.5 micron, which we interpret as a lower-limit for the wavelength of peak emission. We used CLUMPY torus models to fit the nuclear infrared (IR) SED and infer trends in the physical parameters of the AGN torus for the galaxies in the sample. Including the 31.5 micron nuclear flux in the SED 1) reduces the number of clumpy torus models compatible with the data, and 2) modifies the model output for the outer radial extent of the torus for 10 of the 11 objects. Specifically, six (60%) objects show a decrease in radial extent while four (40%) show an increase. We find torus outer radii ranging from <1pc to 8.4 pc, Comment: 15 pages, 7 figures. Accepted for publication in MNRAS