Search

We present a 19.7 - 214 $\mu$m imaging atlas of local (4 - 181 Mpc; median 43 Mpc) active galactic nuclei (AGN) observed with FORCAST and HAWC+ on board the SOFIA telescope with angular resolutions ~ 3"- 20". This atlas comprises 22 Seyferts (17 Type 2 and 5 Type 1) with a total of 69 images, 41 of which have not been previously published. The AGN span a range of luminosities of log$_{10}$ ($L_{bol}$ [erg/s]) = [42, 46] with a median of log$_{10}$ ($L_{bol}$ [erg/s]) = 44.1 $\pm$ 1.0. We provide total fluxes of our sample using aperture photometry for point source objects and a 2-D Gaussian fitting for objects with extended host galaxy emission, which was used to estimate the unresolved nuclear component. Most galaxies in our sample are point-like sources, however, four sources (Centaurus A, Circinus, NGC 1068, and NGC 4388) show extended emission in all wavelengths. The 30 - 40 $\mu$m extended emission in NGC 4388 is coincident with the narrow line region at PA ~ 50$^{\circ}$, while the dusty extension at longer wavelengths arises from the host galaxy at PA ~ 90$^{\circ}$. Our new observations allow us to construct the best sampled spectral energy distributions (SEDs) available between 30 - 500 $\mu$m for a sample of nearby AGN. We estimate that the average peak wavelength of the nuclear SEDs is ~ 40 $\mu$m in $\nu$F$_{\nu}$ , which we associate with an unresolved extended dusty region heated by the AGN., Comment: 27 pages, 15 figures

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: ApJL accepted on September 26th, title slightly modified in accordance with ApJL standards, Fig. 2 updated with additional labels

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)

Cold gas and dust reprocess the X-ray emission of active galactic nuclei (AGN), producing spectro-polarimetric features in the X-ray band. The recent launch of IXPE allows for observations of this X-ray polarisation signal, encoding unique information on the circumnuclear medium of AGN. However, the models for interpreting these polarimetric data are under-explored and do not reach the same level of sophistication as the corresponding spectral models. We aim at closing the gap between the spectral and spectro-polarimetric modelling of AGN circumnuclear media by providing the tools for simulating X-ray polarisation in complex 3D transfer media alongside X-ray spectra. We lay out the framework for X-ray polarisation in 3D radiative transfer simulations and provide an implementation to the 3D radiative transfer code SKIRT, focussing on (de)polarisation due to scattering and fluorescent re-emission. As a first application, we studied a 2D toroidal reprocessor of cold gas, modelling the AGN circumnuclear medium. For the 2D torus model, we find a complex behaviour of the polarisation angle with photon energy, which we interpret as a balance between the reprocessed photon flux originating from different sky regions, with a direct link to the torus geometry. We calculated a large grid of AGN torus models and demonstrated how spatially resolved polarisation maps could form a useful tool for interpreting the geometrical information that is encoded in IXPE observations. With this work, we release high-resolution AGN torus templates that simultaneously describe X-ray spectra and spectro-polarimetry, for observational data fitting with XSPEC. The SKIRT code can now model X-ray polarisation simultaneously with X-ray spectra and provide synthetic spectro-polarimetric observations for complex 3D circumnuclear media, with all features of the established SKIRT framework available., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical properties, and the grain size distribution of dust, all of which can have a substantial impact on the theoretical radiative transfer calculations. In this work, we explore the role of the dust chemical composition in the AGN dusty torus through SED fitting to \emph{Spitzer}/IRS spectra of a sample of 49 nearby AGN with silicate features in emission. We implement a mineralogy model including the popular astronomical silicates and a set of oxides and amorphous silicates with different grain sizes. We find that best fits use principally porous alumina, periclase, and olivine. In terms of mass fractions, $\sim99\%$ of the dust is composed of dust grains of size $\rm{0.1 \mu m}$, with a $<1\%$ contribution from $\rm{3 \mu m}$ grains. Moreover, the astronomical silicates have a very low occurrence in the best fits, suggesting that they are not the most suited dust species to reproduce the silicate features in our sample., Comment: 16 pages, 13 figures, Accepted by MNRAS

Polar dust has been discovered in a number of local Active Galactic Nuclei (AGN), with radiation-driven torus models predicting a wind to be its main driver. However, little is known about its characteristics, spatial extent, or connection to the larger scale outflows. We present the first JWST/MIRI study aimed at imaging polar dust by zooming onto the centre of ESO 428-G14, part of the Galaxy Activity, Torus, and Outflow Survey (GATOS) survey of local AGN. We detect extended mid-infrared (MIR) emission within 200 pc from the nucleus. This polar structure is co-linear with a radio jet and lies perpendicular to a molecular gas lane that feeds and obscures the nucleus. Its morphology bears a striking resemblance to that of gas ionised by the AGN in the narrow-line region. We demonstrate that part of this spatial correspondence is due to contamination within the JWST filter bands from strong emission lines. Correcting for the contamination, we find the morphology of the dust continuum to be more compact, though still clearly extended out to ~ 100 pc. We estimate the emitting dust has a temperature of ~ 120 K. Using simple models, we find that the heating of small dust grains by the radiation from the central AGN and/or radiative jet-induced shocks is responsible for the extended MIR emission. Radiation-driven dusty winds from the torus is unlikely to be important. This has important implications for scales to which AGN winds can carry dust and dense gas out into their host galaxies., Comment: 17 pages, 8 figures, 2 tables. Published in MNRAS

In this paper we present the first-ever $L$- and $M$-band interferometric observations of Circinus, building upon a recent $N$-band analysis. We used these observations to reconstruct images and fit Gaussian models to the $L$ and $M$ bands. Our findings reveal a thin edge-on disk whose width is marginally resolved and is the spectral continuation of the disk imaged in the $N$ band to shorter wavelengths. Additionally, we find a point-like source in the $L$ and $M$ bands that, based on the $LMN$-band spectral energy distribution fit, corresponds to the $N$-band point source. We also demonstrate that there is no trace of direct sightlines to hot dust surfaces in the circumnuclear dust structure of Circinus. By assuming the dust is present, we find that obscuration of A$_{\rm V} \gtrsim 250$ mag is necessary to reproduce the measured fluxes. Hence, the imaged disk could play the role of the obscuring "torus" in the unified scheme of active galactic nuclei. Furthermore, we explored the parameter space of the disk + hyperbolic cone radiative transfer models and identify a simple modification at the base of the cone. Adding a cluster of clumps just above the disk and inside the base of the hyperbolic cone provides a much better match to the observed temperature distribution in the central aperture. This aligns well with the radiation-driven fountain models that have recently emerged. Only the unique combination of sensitivity and spatial resolution of the VLTI allows such models to be scrutinized and constrained in detail. We plan to test the applicability of this detailed dust structure to other MATISSE-observed active galactic nuclei in the future., Comment: Main article: 12 pages, 6 figures. Accepted for publication in A&A. arXiv admin note: text overlap with arXiv:2205.01575

Fits the infrared spectra from the nuclear regions of AGN can place constraints on the dust properties, distribution, and geometry by comparison with models. However, none of the currently available models fully describe the observations of AGN currently available. Among the aspects least explored, here we focus on the role of dust grain size. We offer the community a new spectral energy distribution (SED) library, hereinafter [GoMar23] model, which is based on the two-phase torus model developed before with the inclusion of the grain size as a model parameter, parameterized by the maximum grain size Psize or equivalently the mass-weighted average grain size < P >. We created 691,200 SEDs using the SKIRT code, where the maximum grain size can vary within the range Psize = 0.01 - 10.0um ( < P >= 0.007 - 3.41um). We fit this new and several existing libraries to a sample of 68 nearby and luminous AGNs with Spitzer/IRS spectra dominated by AGN-heated dust. We find that the [GoMar23] model can adequately reproduce up to 85-88% of the spectra. The dust grain size parameter significantly improves the final fit in up to 90% of these spectra. Statistical tests indicate that the grain size is the third most important parameter in the fitting procedure (after the size and half opening angle of the torus). The requirement of a foreground extinction by our model is lower compared to purely clumpy models. We find that 41% of our sample requires that the maximum dust grain size is as large as Psize =10um (< P >= 3.41um). Nonetheless, we also remark that disk+wind and clumpy torus models are still required to reproduce the spectra of a non-negligible fraction of objects, suggesting the need for several dust geometries to explain the infrared continuum of AGN. This work provides tentative evidence for dust grain growth in the proximity of the AGN., Comment: 26 pages, 14 figures, 4 tables, accepted for publication in A&A

Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT with the X-ray processes that govern the broadband X-ray spectra of obscured AGN, to study the structure of AGN circumnuclear media in full 3D, based on their reflected X-ray emission. We extended the SKIRT code with Compton scattering on free electrons, photo-absorption and fluorescence by cold atomic gas, scattering on bound electrons, and extinction by dust. This includes a novel treatment of extreme-forward scattering by dust, and a detailed description of anomalous Rayleigh scattering. To verify our X-ray implementation, we performed the first dedicated benchmark of X-ray torus models, comparing five X-ray radiative transfer codes. Finally, we illustrated the 3D nature of the code by producing synthetic X-ray images and spectra of clumpy torus models. SKIRT forms a powerful new tool to model AGN circumnuclear media in full 3D from X-ray to millimetre wavelengths, and is now publicly available. In the X-ray regime, we find an excellent agreement with the simulation results of the MYTorus and RefleX codes, which validates our X-ray implementation. We find some discrepancies with other codes, which motivates the need for a robust framework that can handle non-linear 3D radiative transfer effects. The new X-ray functionalities of the SKIRT code allow for uncomplicated access to a broad suite of 3D X-ray models for AGN that can easily be tested and modified. This will be particularly useful with the advent of X-ray microcalorimeter observations., Comment: 20 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

We present polarimetric maps of the Circinus galaxy nucleus in the $BVRI$ bands, obtained with VLT/FORS2. Circinus is the closest Seyfert 2 galaxy and harbours an archetypal obscured active galactic nucleus (AGN). Recent high angular resolution imaging revealed that a major fraction of its mid-infrared (MIR) emission is coming from the polar region. Previously, we demonstrated that these observations are consistent with a model of a compact dusty disc and a hyperboloid shell, resembling a hollow cone on larger scales. Here we focus on the AGN core, up to 40 pc from the central engine, and compare the observations to the radiative transfer models. Polarization maps reveal a conical structure, coinciding with the ionization cone. The wavelength-dependence of the polarization degree indicates that scattering on dust grains is producing polarization. The observed polarization degree ($\approx1-3\%$) is lower than predicted by the models; however, this is only a lower limit, since stellar emission dominates the total flux in the optical. The observed polarization angle ($\approx30$ degrees) is reproduced by the model of a dusty disc with a hollow cone that is illuminated by a tilted anisotropic central source. An accretion disc aligned with the ionization cone axis, and alternative dust geometries, such as a paraboloid shell, or a torus enveloped by ambient dust, are inconsistent with the data. We conclude that the optical polarimetric imaging supports earlier evidence for the presence of dust in the polar region, tentatively associated with dusty outflows., Comment: MNRAS (accepted)

Aims. Monte Carlo Radiative Transfer (MCRT) simulations are a powerful tool for understanding the role of dust in astrophysical systems and its influence on observations. However, due to the strong coupling of the radiation field and medium across the whole computational domain, the problem is non-local and non-linear and such simulations are computationally expensive in case of realistic 3D inhomogeneous dust distributions. We explore a novel technique for post-processing MCRT output to reduce the total computational run time by enhancing the output of computationally less expensive simulations of lower-quality. Methods. We combine principal component analysis (PCA) and non-negative matrix factorization (NMF) as dimensionality reduction techniques together with Gaussian Markov random fields and the Integrated nested Laplace approximation (INLA), an approximate method for Bayesian inference, to detect and reconstruct the non-random spatial structure in the images of lower signal-to-noise or with missing data. Results. We test our methodology using synthetic observations of a galaxy from the SKIRT Auriga project - a suite of high resolution magneto-hydrodynamic Milky Way-sized galaxies simulated in cosmological environment by 'zoom-in' technique. With this approach, we are able to reproduce high photon number reference images $\sim5$ times faster with median residuals below $\sim20\%$., Comment: To be published in Numerical methods and codes of Astronomy and Astrophysics

Dust is a major component of the interstellar medium. Through scattering, absorption and thermal re-emission, it can profoundly alter astrophysical observations. Models for dust composition and distribution are necessary to better understand and curb their impact on observations. A new approach for serial and computationally inexpensive production of such models is here presented. Traditionally these models are studied with the help of radiative transfer modelling, a critical tool to understand the impact of dust attenuation and reddening on the observed properties of galaxies and active galactic nuclei. Such simulations present, however, an approximately linear computational cost increase with the desired information resolution. Our new efficient model generator proposes a denoising variational autoencoder (or alternatively PCA), for spectral compression, combined with an approximate Bayesian method for spatial inference, to emulate high information radiative transfer models from low information models. For a simple spherical dust shell model with anisotropic illumination, our proposed approach successfully emulates the reference simulation starting from less than 1% of the information. Our emulations of the model at different viewing angles present median residuals below 15% across the spectral dimension, and below 48% across spatial and spectral dimensions. EmulART infers estimates for ~85% of information missing from the input, all within a total running time of around 20 minutes, estimated to be 6x faster than the present target high information resolution simulations, and up to 50x faster when applied to more complicated simulations., Comment: 85 pages, 33 figures, 9 tables

We present the BAT AGN Spectroscopic Survey (BASS) Near-infrared Data Release 2 (DR2), a study of 168 nearby ($\bar z$ = 0.04, $z$ < 0.6) active galactic nuclei (AGN) from the all-sky Swift Burst Array Telescope X-ray survey observed with Very Large Telescope (VLT)/X-shooter in the near-infrared (NIR; 0.8 - 2.4 $\mu$m). We find that 49/109 (45%) Seyfert 2 and 35/58 (60%) Seyfert 1 galaxies observed with VLT/X-shooter show at least one NIR high-ionization coronal line (CL, ionization potential $\chi$ > 100 eV). Comparing the emission of the [Si vi] $\lambda$1.9640 CL with the X-ray emission for the DR2 AGN, we find a significantly tighter correlation, with a lower scatter (0.37 dex) than for the optical [O iii] $\lambda$5007 line (0.71 dex). We do not find any correlation between CL emission and the X-ray photon index $\Gamma$. We find a clear trend of line blueshifts with increasing ionization potential in several CLs, such as [Si vi] $\lambda$1.9640, [Si x] $\lambda$1.4300, [S viii] $\lambda$0.9915, and [S ix] $\lambda$1.2520, indicating the radial structure of the CL region. Finally, we find a strong underestimation bias in black hole mass measurements of Sy 1.9 using broad H$\alpha$ due to the presence of significant dust obscuration. In contrast, the broad Pa$\alpha$ and Pa$\beta$ emission lines are in agreement with the $M$-$\sigma$ relation. Based on the combined DR1 and DR2 X-shooter sample, the NIR BASS sample now comprises 266 AGN with rest-frame NIR spectroscopic observations, the largest set assembled to date., Comment: 36 pages, 24 figures, 10 tables, published in ApJS as part of BASS DR2 special issue

Active galactic nuclei play a key role in the evolution of galaxies, but their inner workings and physical connection to the host are poorly understood due to a lack of angular resolution. Infrared interferometry makes it possible to resolve the circumnuclear dust in the nearby Seyfert 2 galaxy, Circinus. Previous observations have revealed complex structures and polar dust emission but interpretation was limited to simple models. MATISSE makes it possible to image these structures for the first time. We observed the Circinus Galaxy with VLTI/MATISSE, producing 150 correlated flux spectra and 100 closure phase spectra. We reconstructed images in the N-band at ~10 mas resolution. We fit blackbody functions with dust extinction to several aperture-extracted fluxes from the images to produce a temperature distribution of central dusty structures. We find significant substructure in the circumnuclear dust: central unresolved flux of ~0.5 Jy, a thin disk 1.9 pc in diameter oriented along ~45 deg,and a ~4x1.5 pc polar emission extending orthogonal to the disk. The polar emission exhibits patchiness, which we attribute to clumpy dust. Flux enhancements to the east and west of the disk are seen for the first time. We distinguish the temperature profiles of the disk and of the polar emission: the disk shows a steep temperature gradient indicative of denser material; the polar profile is flatter, indicating clumpiness and/or lower dust density. The unresolved flux is fitted with a high temperature, ~370 K. The polar dust remains warm (~200 K) out to 1.5 pc from the disk. The recovered morphology and temperature distribution resembles modeling of accretion disks with radiation-driven winds at large scales, but we placed new constraints on the subparsec dust. The subparsec features imaged here place new constraints on the physical modeling of circumnuclear dust in active galaxies., Comment: Accepted for publication in A&A

Reprocessed X-ray emission in Active Galactic Nuclei (AGN) can provide fundamental information about the circumnuclear environments of supermassive black holes. Recent mid-infrared studies have shown evidence of an extended dusty structure perpendicular to the torus plane. In this work, we build a self-consistent X-ray model for the Circinus Galaxy including the different physical components observed at different wavelengths and needed to reproduce both the morphological and spectral properties of this object in the mid-infrared. The model consists of four components: the accretion disk, the broad line region (BLR), a flared disk in the equatorial plane and a hollow cone in the polar direction. Our final model reproduces well the 3--70 keV Chandra and NuSTAR spectra of Circinus, including the complex Fe K$\alpha$ zone and the spectral curvature, although several additional Gaussian lines, associated to either ionized iron or to broadened Fe K$\alpha$/K$\beta$ lines, are needed. We find that the flared disk is Compton thick ($ N_{\rm H,d}= \rm 1.01^{+0.03}_{-0.24}\times 10^{25}\: cm^{-2}$) and geometrically thick ($CF=0.55^{+0.01}_{-0.05}$), and that the hollow cone has a Compton-thin column density ($ N_{\rm H,c}= \rm 2.18^{+0.47}_{-0.43}\times 10^{23}\: cm^{-2}$), which is consistent with the values inferred by mid-infrared studies. Including also the BLR, the effective line of sight column density is $ N_{\rm H}= \rm 1.47^{+0.03}_{-0.24}\times 10^{25}\: cm^{-2}$. This approach to X-ray modelling, i.e. including all the different reprocessing structures, will be very important to fully exploit data from future X-ray missions., Comment: Accepted for publication in MNRAS. 15 pages, 10 figures, 8 tables

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(\theta) \propto 1+\cos \theta$, where $\theta$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr, Comment: ApJ accepted. The code is available at https://cigale.lam.fr/news/

In this study, we use the SWIFT/BAT AGN sample, which has received extensive multiwavelength follow-up analysis as a result of the BAT AGN Spectroscopic Survey (BASS), to develop a diagnostic for nuclear obscuration by examining the relationship between the line-of-sight column densities ($N_{\rm{H}}$), the 2-10 keV-to-$12\,\rm{\mu m}$ luminosity ratio, and WISE mid-infrared colors. We demonstrate that heavily obscured AGNs tend to exhibit both preferentially ''redder'' mid-infrared colors and lower values of $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{\mu m}}$ than less obscured AGNs, and we derive expressions relating $N_{\rm{H}}$ to the $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{\mu m}}$ and $L_{22\,\rm{\mu m}}$/$L_{4.6\,\rm{\mu m}}$ luminosity ratios as well as develop diagnostic criteria using these ratios. Our diagnostic regions yield samples that are $\gtrsim80$% complete and $\gtrsim60$% pure for AGNs with log($N_{\rm{H}})\geq24$, as well as $\gtrsim85$% pure for AGNs with $\rm{log}(N_{\rm{H}})\gtrsim23.5$. We find that these diagnostics cannot be used to differentiate between optically star forming galaxies and active galaxies. Further, mid-IR contributions from host galaxies that dominate the observed $12~\rm{\mu m}$ emission can lead to larger apparent X-ray deficits and redder mid-IR colors than the AGNs would intrinsically exhibit, though this effect helps to better separate less obscured and more obscured AGNs. Finally, we test our diagnostics on two catalogs of AGNs and infrared galaxies, including the XMM-Newton XXL-N field, and we identify several known Compton-thick AGNs as well as a handful of candidate heavily obscured AGNs based upon our proposed obscuration diagnostics., Comment: 28 pages: 1-19 main text, 20-28 appendix and references. 15 figures. Accepted for publication in ApJS; part of BASS DR2 special issue

We present the largest currently existing subarcsecond 3-5 $\mu$m atlas of 119 local ($z < 0.3$) active galactic nuclei (AGN). This atlas includes AGN of 5 subtypes: 22 are Seyfert 1; 5 are intermediate Seyferts; 46 are Seyfert 2; 26 are LINERs; and 20 are composites/starbursts. Each AGN was observed with VLT ISAAC in the $L$- and/or $M$-bands between 2000 and 2013. We detect at 3$\sigma$ confidence 92 sources in the $L$-band and 83 sources in the $M$-band. We separate the flux into unresolved nuclear flux and resolved flux through two-Gaussian fitting. We report the nuclear flux, extended flux, apparent size, and position angle of each source, giving $3\sigma$ upper-limits for sources which are undetected. Using WISE W1- and W2-band photometry we derive relations predicting the nuclear $L$ and $M$ fluxes for Sy1 and Sy2 AGN based on their W1-W2 color and WISE fluxes. Lastly, we compare the measured mid-infrared colors to those predicted by dusty torus models SKIRTOR, CLUMPY, CAT3D, and CAT3D-WIND, finding best agreement with the latter. We find that models including polar winds best reproduce the 3-5$\mu$m colors, indicating that winds are an important component of dusty torus models. We find that several AGN are bluer than models predict. We discuss several explanations for this and find that it is most plausibly stellar light contamination within the ISAAC $L$-band nuclear fluxes., Comment: Main Text: 22 pages, 9 figures, 3 tables. Accepted by ApJ

The scatter (${\rm\sigma_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm \sigma_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm \sigma_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm \sigma_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$\beta$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation., Comment: 23 pages, 10 figures, Accepted to ApJ

Recent observations which resolved the mid-infrared (MIR) emission of nearby active galactic nuclei (AGN), surprisingly revealed that their dust emission appears prominently extended in the polar direction, at odds with the expectations from the canonical dusty torus. This polar dust, tentatively associated with dusty winds driven by radiation pressure, is found to have a major contribution to the MIR flux from scales of a few to hundreds of parsecs. When facing a potential change of paradigm, case studies of objects with the best intrinsic resolution are essential. One such source with a clear detection of polar dust is a nearby, well-known AGN in the Circinus galaxy. In the first paper, we successfully explained the peculiar MIR morphology of Circinus observed on large, tens of parsec scales with a model consisting of a compact dusty disc and an extended hollow dusty cone. In this work, we further refine the model on smaller, parsecs scales to test whether it can also explain the MIR interferometric data. We find that a model composed of a thin dusty disc seen almost edge-on and a polar outflow in the form of a hyperboloid shell can reproduce well the VLTI/MIDI observations at all wavelengths, baselines and position angles. In contrast, while providing a good fit to the integrated MIR spectrum, the dusty torus model fails to reproduce the spatially resolved interferometric data. We put forth the disc$+$hyperboloid wind model of Circinus AGN as a prototype for the dust structure in the AGN population with polar dust., Comment: MNRAS accepted

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

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Recent high angular resolution observations resolved for the first time the mid-infrared (MIR) structure of nearby active galactic nuclei (AGN). Surprisingly, they revealed that a major fraction of their MIR emission comes from the polar regions. This is at odds with the expectation based on AGN unification, which postulates a dusty torus in the equatorial region. The nearby, archetypical AGN in the Circinus galaxy offers one of the best opportunities to study the MIR emission in greater detail. New, high quality MIR images obtained with the upgraded VISIR instrument at the Very Large Telescope show that the previously detected bar-like structure extends up to at least 40 pc on both sides of the nucleus along the edges of the ionization cone. Motivated by observations across a wide wavelength range and on different spatial scales, we propose a phenomenological dust emission model for the AGN in the Circinus galaxy consisting of a compact dusty disk and a large-scale dusty cone shell, illuminated by a tilted accretion disk with an anisotropic emission pattern. Undertaking detailed radiative transfer simulations, we demonstrate that such a model is able to explain the peculiar MIR morphology and account for the entire IR spectral energy distribution. Our results call for caution when attributing dust emission of unresolved sources entirely to the torus and warrant further investigation of the MIR emission in the polar regions of AGN., Comment: Accepted to MNRAS. Version 2: typos corrected

Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy., Comment: 15 pages, 10 figures, accepted by Astronomy and Astrophysics

The primary source of emission of active galactic nuclei (AGN), the accretion disk, is surrounded by an optically and geometrically thick dusty structure ("the so-called dusty torus"). The infrared radiation emitted by the dust is nothing but a reprocessed fraction of the accretion disk emission, so the ratio of the torus to the AGN luminosity ($L_{\text{torus}}/L_{\text{AGN}}$) should correspond to the fraction of the sky obscured by dust, i.e. the covering factor. We undertook a critical investigation of the $L_{\text{torus}}/L_{\text{AGN}}$ as the dust covering factor proxy. Using state-of-the-art 3D Monte Carlo radiative transfer code, we calculated a grid of spectral energy distributions (SEDs) emitted by the clumpy two-phase dusty structure. With this grid of SEDs, we studied the relation between $L_{\text{torus}}/L_{\text{AGN}}$ and the dust covering factor for different parameters of the torus. We found that in the case of type 1 AGNs the torus anisotropy makes $L_{\text{torus}}/L_{\text{AGN}}$ underestimate low covering factors and overestimate high covering factors. In type 2 AGNs $L_{\text{torus}}/L_{\text{AGN}}$ always underestimates covering factors. Our results provide a novel easy-to-use method to account for anisotropy and obtain correct covering factors. Using two samples from the literature, we demonstrated the importance of our result for inferring the obscured AGN fraction. We found that after the anisotropy is properly accounted for, the dust covering factors show very weak dependence on $L_{\text{AGN}}$, with values in the range of $\approx0.6-0.7$. Our results also suggest a higher fraction of obscured AGNs at high luminosities than those found by X-ray surveys, in part owing to the presence of a Compton-thick AGN population predicted by population synthesis models., Comment: Published in MNRAS. V2: references updated and typos corrected, to match the published version

The narrow component of the iron K$\alpha$ is an almost ubiquitous feature in the X-ray spectra of active galactic nuclei (AGN) and is believed to originate in neutral material, possibly located in the molecular torus. This would imply a tight connection between the Fe K$\alpha$ equivalent width (EW) and the physical properties of the torus. In a recent work we have shown that the decrease of the covering factor of the torus with the luminosity, as expected by luminosity-dependent unification models, would be able to explain the decrease of Fe K$\alpha$ EW with the luminosity (i.e., the X-ray Baldwin effect). Recent developments in the study of the mid-IR (MIR) spectrum of AGN allow important parameters of the torus to be deduced, such as its covering factor ($f_{\rm\,obs}$) and equatorial column density ($N_{\rm\,H}^{\rm\,T}$), by applying clumpy torus models. Using XMM-Newton/EPIC observations of a sample of 24 type-I AGN, we investigate the relation between the physical parameters of the torus obtained by recent MIR works and the properties of the Fe K$\alpha$ line. We correct the values of the Fe K$\alpha$ EW by taking the inclination angle, the photon index, the equatorial column density, and half-opening angle of the torus into account using a physical torus model of X-ray reprocessed radiation. We find that the relation between Fe K$\alpha$ EW and $f_{\rm\,obs}$ shows a slope that is consistent with the expected value, albeit with a low statistical significance. A trend that is consistent with the theoretical prediction is also found when comparing the Fe K$\alpha$ EW to $N_{\rm\,H}^{\rm\,T}$. Our work seems to confirm that the bulk of the narrow Fe K$\alpha$ line is produced by the same material responsible for the MIR emission., Comment: A&A in press, 15 pages, 5 Figures, 3 tables - Few references updated

We modeled the AGN dusty torus as a clumpy two-phase medium, with high-density clumps embedded in a low-density interclump dust. To obtain spectral energy distributions and images of the torus at different wavelengths, we employed the 3D Monte Carlo radiative transfer code SKIRT. Apart from the grid of two-phase models, we calculated the corresponding sets of clumps-only models and models with a smooth dust distribution for comparison. We found that the most striking feature of the two-phase model is that it might offer a natural solution to the common issue reported in a number of papers -- the observed excess of the near-infrared emission., Comment: Proceedings of the Torus Workshop 2012 held at the University of Texas at San Antonio, 5-7 December 2012. C. Packham, R. Mason, and A. Alonso-Herrero (eds.). 8 pages, 5 figures. A grid of model SEDs available at https://sites.google.com/site/skirtorus/

We investigated the gravitational microlensing of active galactic nucleus dusty tori in the case of lensed quasars in the infrared domain. The dusty torus is modeled as a clumpy two-phase medium. To obtain spectral energy distributions and images of tori at different wavelengths, we used the 3D Monte Carlo radiative transfer code SKIRT. A ray-shooting technique has been used to calculate microlensing magnification maps. We simulated microlensing by the stars in the lens galaxy for different configurations of the lensed system and different values of the torus parameters, in order to estimate (a) amplitudes and timescales of high magnification events, and (b) the influence of geometrical and physical properties of dusty tori on light curves in the infrared domain. We found that, despite their large size, dusty tori could be significantly affected by microlensing in some cases, especially in the near-infrared domain (rest-frame). The very long time-scales of such events, in the range from several decades to hundreds of years, are limiting the practical use of this method to study the properties of dusty tori. However, our results indicate that, when studying flux ratios between the images in different wavebands of lensed quasars, one should not disregard the possibility that the near and mid-infrared flux ratios could be under the influence of microlensing., Comment: 10 pages, 6 figures, 1 table, MNRAS accepted. V3: corrected values in Table 1

We investigated the emission of active galactic nuclei dusty tori in the infrared domain, with a focus on the 10 micron silicate feature. We modeled the dusty torus as a clumpy two-phase medium with high-density clumps and a low-density medium filling the space between the clumps. We employed a three-dimensional radiative transfer code to obtain spectral energy distributions and images of tori at different wavelengths. We calculated a grid of models for different parameters and analyzed the influence of these parameters on the shape of the mid-infrared emission. A corresponding set of clumps-only models and models with a smooth dust distribution is calculated for comparison. We found that the dust distribution, the optical depth and a random arrangement of clumps in the innermost region, all have an impact on the shape and strength of the silicate feature. The 10 micron silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the silicate feature strength., Comment: 5 pages, 2 figures. Proceedings of the "8th Serbian Conference on Spectral Line Shapes in Astrophysics", Divcibare, Serbia, June 6-10 2011. Model SEDs available for download at https://sites.google.com/site/skirtorus/

We study the photocenter position variability due to variations in the quasar inner structure. We consider variability in the accretion disk emissivity and torus structure variability due to different illumination by the central source. We discuss possible detection of these effects by Gaia. Observations of the photocenter variability in two AGNs, SDSS J121855+020002 and SDSS J162011+1724327 have been reported and discussed. With investigation of the variations in the quasar inner structure we explore how much this effect can affect the position determination and whether it can be (or not) detected with Gaia mission. We used (a) a model of a relativistic disk, including the perturbation that can increase brightness of a part of the disk, and consequently offset the photocenter position, and (b) a model of a dusty torus which absorbs and re-emits the incoming radiation from accretion disk. We estimated the value of the photocenter offset due to these two effects. We found that perturbations in the inner structure can significantly offset the photocenter. It depends on the characteristics of perturbation and accretion disk and structure of the torus. In the case of two considered QSOs the observed photocenter offsets cannot be explained by variations in the accretion disk and other effects should be considered. We discussed possibility of exploding stars very close to the AGN source, and also possibility that there are two variable sources in the center of these two AGNs that may indicate a binary super-massive black hole system on a kpc (pc) scale. The Gaia mission seems to be very perspective, not only for astrometry, but also for exploring the inner structure of AGNs. We conclude that variations in the quasar inner structure can affect the observed photocenter (up to several mas). There is a chance to observe such effect in the case of bright and low-redshifted QSOs., Comment: 12 pages, 8 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics. Language improved, typos corrected

We investigate the emission of active galactic nuclei (AGN) dusty tori in the infrared domain. Following theoretical predictions coming from hydrodynamical simulations, we model the dusty torus as a 3D two-phase medium with high-density clumps and low-density medium filling the space between the clumps. Spectral energy distributions (SED) and images of the torus at different wavelengths are obtained using 3D Monte Carlo radiative transfer code SKIRT. Our approach of generating clumpy structure allows us to model tori with single clumps, complex structures of merged clumps or interconnected sponge-like structure. A corresponding set of clumps-only models and models with smooth dust distribution is calculated for comparison. We found that dust distribution, optical depth, clump size and their actual arrangement in the innermost region, all have an impact on the shape of near- and mid-infrared SED. The 10 micron silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the silicate feature strength. Finally, we find that having the dust distributed in a two-phase medium, might offer a natural solution to the lack of emission in the near-infrared, compared to observed data, which affects clumpy models currently available in the literature., Comment: 18 pages, 15 figures, 1 table. Accepted for publication in MNRAS. Model SEDs available for download at https://sites.google.com/site/skirtorus/

We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures., Comment: 15 pages, 7 figures, accepted for publication in ApJS

Interstellar dust grains efficiently absorb and scatter UV and optical radiation in galaxies, and therefore can significantly affect the apparent structure of spiral galaxies. We discuss the effect of dust attenuation on the observed structural properties of bulges and discs. We also present some first results on modelling the dust content of edge-on spiral galaxies using both optical and Herschel far-infrared data. Both of these results demonstrate the complex interplay of dust and star light in spiral galaxies., Comment: 6 pages, 2 figures, to appear in "Galaxies and their Masks", eds. D. L. Block, K. C. Freeman & I. Puerari, Springer (New York)