Around Active Galactic Nuclei Research Articles

Multi-epoch hard X-ray view of Compton-thick AGN Circinus Galaxy

ABSTRACT The circumnuclear material around active galactic nuclei (AGNs) is one of the essential components of the obscuration-based unification model. However, our understanding of the circumnuclear material in terms of its geometrical shape, structure, and its dependence on accretion rate is still debated. In this paper, we present the multi-epoch broad-band X-ray spectral modelling of a nearby Compton-thick AGN in Circinus galaxy. We utilize all the available hard X-ray (>10 keV) observations taken from different telescopes, i.e. BeppoSAX, Suzaku, NuSTAR, and AstroSat, at 10 different epochs across 22 yr from 1998 to 2020. The 3.0–79 keV broad-band X-ray spectral modelling using physically motivated models, namely MYTorus, borus02, and uxclumpy, infers the presence of a torus with a low covering factor of 0.28, an inclination angle of 77°–81° and Compton-thick line-of-sight column densities (NH,LOS = 4.13–9.26 × 1024 cm−2) in all the epochs. The joint multi-epoch spectral modelling suggests that the overall structure of the torus is likely to remain unchanged. However, we find tentative evidence for the variable line-of-sight column density on time-scales ranging from 1 d to 1 week to a few years, suggesting a clumpy circumnuclear material located at subparsec to tens of parsec scales.

GOALS-JWST: NIRCam and MIRI Imaging of the Circumnuclear Starburst Ring in NGC 7469

We present James Webb Space Telescope (JWST) imaging of NGC 7469 with the Near-Infrared Camera and the Mid-InfraRed Instrument. NGC 7469 is a nearby, z = 0.01627, luminous infrared galaxy that hosts both a Seyfert Type-1.5 nucleus and a circumnuclear starburst ring with a radius of ∼0.5 kpc. The new near-infrared (NIR) JWST imaging reveals 66 star-forming regions, 37 of which were not detected by Hubble Space Telescope (HST) observations. Twenty-eight of the 37 sources have very red NIR colors that indicate obscurations up to A v ∼ 7 and a contribution of at least 25% from hot dust emission to the 4.4 μm band. Their NIR colors are also consistent with young (<5 Myr) stellar populations and more than half of them are coincident with the mid-infrared (MIR) emission peaks. These younger, dusty star-forming regions account for ∼6% and ∼17% of the total 1.5 and 4.4 μm luminosity of the starburst ring, respectively. Thanks to JWST, we find a significant number of young dusty sources that were previously unseen due to dust extinction. The newly identified 28 young sources are a significant increase compared to the number of HST-detected young sources (4–5). This makes the total percentage of the young population rise from ∼15% to 48%. These results illustrate the effectiveness of JWST in identifying and characterizing previously hidden star formation in the densest star-forming environments around active galactic nuclei (AGN).

Torus and polar dust dependence on active galactic nucleus properties

We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14–195 keV; Swift/BAT), volume-limited (DL &lt; 40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy, and two-phase torus models and clumpy disc+wind models. We find that torus models (i.e. without including the polar dusty wind component) and disc+wind models provide the best fits for a comparable number of galaxies, 8 out of 24 (33.3%) and 9 out of 24 (37.5%), respectively. We find that the best-fit models depend on the hydrogen column density (NHX−ray), which is related to the X-ray (unobscured or obscured) and optical (Sy1/Sy2) classification. In particular, smooth, clumpy, and two-phase torus models are best at reproducing the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NHX−ray cm−2) = 23.5 ± 0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NHX−ray cm−2) = 21.0 ± 1.0) – specifically, in the near-infrared (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, namely, nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy, and two-phase torus (disc+wind) models, confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depends on the intrinsic AGN properties.

XMM-Newton and Swift observations of the Seyfert 1 AGN NGC 5940

Probing the physics of the accretion flow around active galactic nuclei (AGN) is crucial to understanding their emission mechanisms as well as being able to constrain the geometrical and variability properties of the different regions around them. The soft X-ray excess -- usually observed below $\sim2\,\mathrm{keV}$ in excess of the dominant X-ray powerlaw continuum -- is one prominent feature that is commonly seen in type 1 Seyfert AGN and therefore readily provides a useful diagnostic of the accretion flow mechanism around these systems. NGC 5940 is a Seyfert 1 AGN which reveals strong, prominent soft X-ray excess below $\sim2\,\mathrm{keV}$ as seen in both its XMM-Newton and Swift observations. Model fit to the data revealed that this feature could be equally well explained by the ionised partial covering, the thermal Comptonisation and the blurred reflection models. Although the other models cannot be decisively ruled out with the data at hand, the lack of significant broad iron $K_{\alpha}$ as well as any significant emission/absorption line features in the reflection grating spectrometer (RGS) data tend to favour the thermal Comptonisation origin for the soft X-ray excess in NGC 5940.

Hypercubes of AGN Tori (HYPERCAT). II. Resolving the Torus with Extremely Large Telescopes

Abstract Recent infrared interferometric observations revealed sub-parsec scale dust distributions around active galactic nuclei (AGNs). Using images of Clumpy torus models and NGC 1068 as an example, we demonstrate that the near- and mid-infrared nuclear emission of some nearby AGNs will be resolvable in direct imaging with the next generation of 30 m telescopes, potentially breaking degeneracies from previous studies that used integrated spectral energy distributions of unresolved AGN tori. To that effect we model wavelength-dependent point spread functions from the pupil images of various telescopes: James Webb Space Telescope, Keck, Giant Magellan Telescope, Thirty Meter Telescope, and Extremely Large Telescope. We take into account detector pixel scales and noise, and apply deconvolution techniques for image recovery. We also model 2D maps of the 10 μm silicate feature strength, S 10, of NGC 1068 and compare with observations. When the torus is resolved, we find S 10 variations across the image. However, to reproduce the S 10 measurements of an unresolved torus a dusty screen of A V &gt; 9 mag is required. We also fit the first resolved image of the K-band emission in NGC 1068 recently published by the GRAVITY Collaboration, deriving likely model parameters of the underlying dust distribution. We find that both (1) an elongated structure suggestive of a highly inclined emission ring, and (2) a geometrically thin but optically thick flared disk where the emission arises from a narrow strip of hot cloud surface layers on the far inner side of the torus funnel, can explain the observations.

Circumnuclear Multi-phase Gas in Circinus Galaxy IV: Absorption Owing to High-J CO Rotational Transitions

Abstract We studied the absorption features of CO lines against the continuum originating from the heated dust in the obscuring tori around active galactic nuclei (AGNs). We investigated the formation of absorption lines corresponding to the CO rotational transitions using three-dimensional non-LTE line transfer simulations considering the dust thermal emission. As in Papers I–III of this series, we performed post-processed radiative transfer calculations using the “radiation-driven fountain model” (Wada et al. 2016), which yields a geometrically thick obscuring structure around the nucleus. This model is consistent with the spectral energy distribution of the nearest type-2 Seyfert galaxy, the Circinus galaxy. We found that the continuum-subtracted channel maps of J = 4−3 and higher transitions show absorption regions along the disk midplane for an edge-on viewing angle. The spectra consist of multiple absorption and emission features, reflecting the internal inhomogeneous and turbulent structure of the torus. The deepest absorption feature is caused by the gas on the near side of the torus between r = 10 and 15 pc, which is located in front of the AGN-heated dust inside r ≃ 5 pc. We also found that a spatial resolution of 0.5–1.0 pc is necessary to resolve the absorption features. Moreover, the inclination angle must be close to the edge-on angle (i.e., ≳85°) to observe the absorption features. The findings of the present study imply that combining our radiation-hydrodynamic model with high-resolution observations of CO (7–6) by ALMA can provide new information about the internal structure of the molecular tori in nearby AGNs.

On the dark matter haloes of optical and IR-selected AGNs in the local universe

ABSTRACT We use the technique of total satellite luminosity, Lsat, to probe the dark matter haloes around active galactic nuclei (AGNs) in the SDSS Main Galaxy Sample. Our results focus on galaxies and AGNs that are the central galaxy of their halo. Our two AGN samples are constructed from optical emission-line diagnostics and from Wide-Field Infrared Survey Explorer (WISE) infrared colours. Both optically selected and WISE-selected AGN have Lsat values twice as high as non-active galaxy samples when controlling for stellar mass and mean stellar age. This implies that the haloes are twice as massive, but we cannot rule out that the increase in Lsat is due to these AGNs residing in younger haloes at the same mass. When only controlling for host galaxy stellar mass, WISE-selected AGNs also have higher Lsat values than optical AGNs at the factor of two level, consistent with previous results comparing the clustering of obscured and unobscured AGNs. However, controlling for stellar age in the two populations of host galaxies removes half of this difference, attenuating the statistical significance of the difference. We perform permutation tests to quantify the difference in the halo populations of each sample. The difference in star formation properties does not fully explain the difference in the two AGN populations, however. Although AGN luminosity correlates with mean stellar age, the difference in stellar age between the WISE and optical samples cannot be fully explained by differences in their AGN luminosity distributions.

Infrared lags in the light curves of AGNs measured using a deep survey

ABSTRACT Information on the structure around active galactic nuclei (AGNs) has long been derived from measuring lags in their varying light output at different wavelengths. In principle, infrared data would reach to larger radii, potentially even probing reprocessed radiation in any surrounding dusty torus. In practice, this has proved challenging because high quality data are required to detect such variability, and the observations must stretch over a long period to probe the likely month-scale lags in variability. In addition, large numbers of sources would need to be observed to start searching for any patterns in such lags. Here, we show that the UKIDSS Ultra Deep Survey, built up from repeated observations over almost a decade, provides an ideal data set for such a study. For 94 sources identified as strongly varying AGNs within its square-degree field, we find that the K-band light curves systematically lag the J-band light curves by an average of around a month. The lags become smaller at higher redshift, consistent with the band shift to optical rest-frame emission. The less luminous AGNs also display shorter lags, as would be expected if their physical size scales with luminosity.

Star formation in accretion discs and SMBH growth

ABSTRACT Accretion discs around active galactic nuclei (AGNs) are potentially unstable to star formation at large radii. We note that when the compact objects formed from some of these stars spiral into the central supermassive black hole (SMBH), there is no radiative feedback and therefore the accretion rate is not limited by radiation forces. Using a set of accretion disc models, we calculate the accretion rate on to the central SMBH in both gas and compact objects. We find that the time-scale for an SMBH to double in mass can decrease by factors ranging from ∼0.7 to as low as ∼0.1 in extreme cases, compared to gas accretion alone. Our results suggest that the formation of extremely massive black holes at high redshift may occur without prolonged super-Eddington gas accretion or very massive seed black holes. We comment on potential observational signatures as well as implications for other observations of AGNs.

Cosmological Simulation of Galaxy Groups and Clusters. I. Global Effect of Feedback from Active Galactic Nuclei

Abstract In this study we quantify the properties of the gas and dark matter around active galactic nuclei (AGNs) in simulated galaxy groups and clusters and analyze the effect of AGN feedback on the surrounding intracluster (group) medium. Our results suggest downsizing of AGN luminosity with host halo mass, supporting the results obtained from clustering studies of AGNs. By examining the temperature and density distribution of the gas in the vicinity of AGNs we show that due to feedback from the central engine, the gas gets displaced from the center of the group/cluster resulting in a reduction of the density but an enhancement of temperature. We show that these effects are pronounced at both high and low redshifts and propose new observables to study the effect of feedback in higher-redshift galaxies. We also show that the average stellar mass is decreased in halos in the presence of AGN feedback confirming claims from previous studies. Our work for the first time uses a fully cosmological hydrodynamic simulation to evaluate the global effects of AGN feedback on their host dark matter halos as well as galaxies at scales of galaxy groups and clusters.

A new tool to derive chemical abundances in type-2 active galactic nuclei

AbstarctWe present a new tool for the analysis of the optical emission lines of the gas in the Narrow Line Region (NLR) around Active Galactic Nuclei (AGNs). This new tool can be used in large samples of objects in a consistent way using different sets of optical emission-lines taking into the account possible variations from the O/H - N/O relation. The code compares certain observed emission-line ratios with the predictions from a large grid of photoionization models calculated under the most usual conditions in the NLR of AGNs to calculate the total oxygen abundance, nitrogen-to-oxygen ratio and ionization parameter. We applied our method to a sample of Seyfert 2 galaxies with optical emission-line fluxes from the literature. Our results confirm the high metallicity of the objects of the sample and provide consistent values with the direct method. The usage of models to calculate precise ICFs is mandatory when only optical emission lines are available to derive chemical abundances using the direct method in NLRs of AGN.

Morphology of AGN emission-line regions in SDSS-IV MaNGA survey

Extended narrow-line regions (NLRs) around active galactic nuclei (AGN) are shaped by the distribution of gas in the host galaxy and by the geometry of the circumnuclear obscuration, and thus they can be used to test the AGN unification model. In this work, we quantify the morphologies of the narrow-line regions in 308 nearby AGNs ($z=0-0.14$, \lbol $\sim 10^{42.4-44.1}$ \erg{}) from the MaNGA survey. Based on the narrow-line region maps, we find that a large fraction (81\%) of these AGN have bi-conical NLR morphology. The distribution of their measured opening angles suggests that the intrinsic opening angles of the ionization cones has a mean value of 85--98$^\circ$ with a finite spread of 39-44$^\circ$ (1-$\sigma$). Our inferred opening angle distribution implies a number ratio of type I to type II AGN of 1:1.6--2.3, consistent with other measurements of the type I / type II ratio at low AGN luminosities. Combining these measurements with the WISE photometry data, we find that redder mid-IR color (lower effective temperature of dust) corresponds to stronger and narrower photo-ionized bicones. This relation is in agreement with the unification model that suggests that the bi-conical narrow-line regions are shaped by a toroidal dusty structure within a few pc from the AGN. Furthermore, we find a significant alignment between the minor axis of host galaxy disks and AGN ionization cones. Together, these findings suggest that obscuration on both circumnuclear ($\sim $pc) and galactic ($\sim$ kpc) scales are important in shaping and orienting the AGN narrow-line regions.

Clustering of galaxies around AGNs in the HSC Wide survey

Abstract We have measured the clustering of galaxies around active galactic nuclei (AGNs) for which single-epoch virial masses of the super-massive black hole (SMBH) are available to investigate the relation between the large-scale environment of AGNs and the evolution of SMBHs. The AGN samples used in this work were derived from the Sloan Digital Sky Survey (SDSS) observations and the galaxy samples were from the 240 deg2 S15b data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). The investigated redshift range is 0.6–3.0, and the masses of the SMBHs lie in the range 107.5–1010 M⊙. The absolute magnitude of the galaxy samples reaches to Mλ310 ∼ −18 at rest-frame wavelength 310 nm for the low-redshift end of the samples. More than 70% of the galaxies in the analysis are blue. We found a significant dependence of the cross-correlation length on redshift, which primarily reflects the brightness-dependence of the galaxy clustering. At the lowest redshifts the cross-correlation length increases from 7 h−1 Mpc around Mλ310 = −19 mag to &amp;gt;10 h−1 Mpc beyond Mλ310 = −20 mag. No significant dependence of the cross-correlation length on BH mass was found for whole galaxy samples dominated by blue galaxies, while there was an indication of BH mass dependence in the cross-correlation with red galaxies. These results provides a picture of the environment of AGNs studied in this paper being enriched with blue star-forming galaxies, and a fraction of the galaxies are evolving into red galaxies along with the evolution of SMBHs in that system.

The mineralogy of newly formed dust in active galactic nuclei

The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 μm silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8±6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources.

Search for Magnetically Broadened Cascade Emission from Blazars with VERITAS

We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGN), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma rays from AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background (CMB) photons via inverse-Compton scattering to produce gamma rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around $10^{-14}$G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.

HIGH-RESOLUTION LINEAR POLARIMETRIC IMAGING FOR THE EVENT HORIZON TELESCOPE

ABSTRACT Images of the linear polarizations of synchrotron radiation around active galactic nuclei (AGNs) highlight their projected magnetic field lines and provide key data for understanding the physics of accretion and outflow from supermassive black holes. The highest-resolution polarimetric images of AGNs are produced with Very Long Baseline Interferometry (VLBI). Because VLBI incompletely samples the Fourier transform of the source image, any image reconstruction that fills in unmeasured spatial frequencies will not be unique and reconstruction algorithms are required. In this paper, we explore some extensions of the Maximum Entropy Method (MEM) to linear polarimetric VLBI imaging. In contrast to previous work, our polarimetric MEM algorithm combines a Stokes I imager that only uses bispectrum measurements that are immune to atmospheric phase corruption, with a joint Stokes Q and U imager that operates on robust polarimetric ratios. We demonstrate the effectiveness of our technique on 7 and 3 mm wavelength quasar observations from the VLBA and simulated 1.3 mm Event Horizon Telescope observations of Sgr A* and M87. Consistent with past studies, we find that polarimetric MEM can produce superior resolution compared to the standard CLEAN algorithm, when imaging smooth and compact source distributions. As an imaging framework, MEM is highly adaptable, allowing a range of constraints on polarization structure. Polarimetric MEM is thus an attractive choice for image reconstruction with the EHT.

A search for pair haloes around active galactic nuclei through a temporal analysis ofFermi-Large Area Telescope data

We develop a method to search for pair halos around active galactic nuclei (AGN) through a temporal analysis of gamma-ray data. The basis of our method is an analysis of the spatial distributions of photons coming from AGN flares and from AGN quiescent states and a further comparison of these two spatial distributions. This method can also be used for a reconstruction of a point spread function (PSF). We found no evidence for a pair halo component through this method by applying it to the Fermi-LAT data in the energy bands of 4.5-6, 6-10, and >10 GeV and set upper limits on the fraction of photons attributable to a pair halo component. An illustration of how to reconstruct the PSF of Fermi-LAT is given. We demonstrate that the PSF reconstructed by using this method is in good agreement with that which was obtained by using the gamma-ray data taken by LAT in the direction of the Crab pulsar and nebula.

SUBMILLIMETER-HCN DIAGRAM FOR ENERGY DIAGNOSTICS IN THE CENTERS OF GALAXIES

ABSTRACT Compiling data from literature and the Atacama Large Millimeter/submillimeter Array archive, we show enhanced HCN(4–3)/HCO+(4–3) and/or HCN(4–3)/CS(7–6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN enhancement). The number of sample galaxies is significantly increased from our previous work. We expect that this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect, were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN to HCO+ and HCN to CS in AGNs as compared to SB galaxies by a factor of a few to even ≳10 are a plausible explanation for the submillimeter HCN enhancement. However, a counterargument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we cannot fully distinguish these two scenarios at this moment owing to an insufficient amount of multi-transition, multi-species data, the former scenario is indicative of abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest that it is difficult with conventional gas-phase X-ray-dominated region ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral–neutral reactions that are efficiently activated in high-temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.

Modelling the broad Fe Kα reverberation in the AGN NGC 4151

The recent detection of X-ray reverberation lags, especially in the Fe Kα line region, around active galactic nuclei (AGN) has opened up the possibility of studying the time-resolved response (reflection) of hard X-rays from the accretion disc around supermassive black holes. Here, we use general relativistic transfer functions for reflection of X-rays from a point source located at some height above the black hole to study the time lags expected as a function of frequency and energy in the Fe Kα line region. We explore the models and the dependence of the lags on key parameters such as the height of the X-ray source, accretion disc inclination, black hole spin and black hole mass. We then compare these models with the observed frequency- and energy-dependence of the Fe Kα line lag in NGC 4151. Assuming the optical reverberation mapping mass of 4.6 × 107 M⊙, we get a best fit to the lag profile across the Fe Kα line in the frequency range (1-2) × 10− 5 Hz for an X-ray source located at a height |$h = 7^{+2.9}_{-2.6} R_G$| with a maximally spinning black hole and an inclination i < 30°.

MODELING THE MOLECULAR COMPOSITION IN AN ACTIVE GALACTIC NUCLEUS DISK

We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{\rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.