Search

The majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust1,2,3. The location and evolution of this obscuring material have been the subject of intense research in the past decades4,5, and are still debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies across the electromagnetic spectrum1,6,7,8. The origin of this trend may be driven by the increase in the inner radius of the obscuring material with incident luminosity, which arises from the sublimation of dust9; by the gravitational potential of the black hole10; by radiative feedback11,12,13,14; or by the interplay between outflows and inflows15. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism that regulates obscuration unclear. Here we report a systematic multi-wavelength survey of hard-X-ray-selected black holes that reveals that radiative feedback on dusty gas is the main physical mechanism that regulates the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas is located within a few to tens of parsecs of the accreting supermassive black hole (within the sphere of influence of the black hole), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.

Changing-look (CL) AGN are unique probes of accretion onto supermassive black holes (SMBHs), especially when simultaneous observations in complementary wavebands allow investigations into the properties of their accretion flows. We present the results of a search for CL behaviour in 412 Swift-BAT detected AGN with multiple epochs of optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS). 125 of these AGN also have 14-195 keV ultra-hard X-ray light-curves from Swift-BAT which are contemporaneous with the epochs of optical spectroscopy. Eight CL events are presented for the first time, where the appearance or disappearance of broad Balmer line emission leads to a change in the observed Seyfert type classification. Combining with known events from the literature, 21 AGN from BASS are now known to display CL behaviour. Nine CL events have 14-195 keV data available, and five of these CL events can be associated with significant changes in their 14-195 keV flux from BAT. The ultra-hard X-ray flux is less affected by obscuration and so these changes in the 14-195 keV band suggest that the majority of our CL events are not due to changes in line-of-sight obscuration. We derive a CL rate of 0.7-6.2 per cent on 10-25 year time-scales, and show that many transitions happen within at most a few years. Our results motivate further multi-wavelength observations with higher cadence to better understand the variability physics of accretion onto SMBHs., 14 pages, 7 figures, plus appendix. v2: updated references. Published in MNRAS

1ES 1927+654 is a paradigm-defying AGN and one of the most peculiar X-ray nuclear transients. In early 2018, this well-known AGN underwent a changing-look event, in which broad optical emission lines appeared and the optical flux increased. Yet, by July 2018, the X-ray flux had dropped by over two orders of magnitude, indicating a dramatic change to the inner accretion flow. With three years of observations with NICER, XMM-Newton, and NuSTAR, we present the X-ray evolution of 1ES 1927+654, which can be broken into three phases-(1) an early super-Eddington phase with rapid variability in X-ray luminosity and spectral parameters, (2) a stable super-Eddington phase at the peak X-ray luminosity, and (3) a steady decline back to the pre-outburst luminosity and spectral parameters. For the first time, we witnessed the formation of the X-ray corona, as the X-ray spectrum transitioned from thermally-dominated to primarily Comptonized. We also track the evolution of the prominent, broad 1 keV feature in the early X-ray spectra and show that this feature can be modeled with blueshifted reflection (z = -0.33) from a single-temperature blackbody irradiating spectrum using xillverTDE, a new flavor of the xillver models. Thus, we propose that the 1 keV feature could arise from reflected emission off the base of an optically thick outflow from a geometrically thick, super-Eddington inner accretion flow, connecting the inner accretion flow with outflows launched during extreme accretion events (e.g. tidal disruption events). Lastly, we compare 1ES 1927+654 to other nuclear transients and discuss applications of xillverTDE to super-Eddington accretors., Comment: 25 pages, 10 figures, accepted for publication in ApJ

We determine the low-redshift X-ray luminosity function (XLF), active black hole mass function (BHMF), and Eddington-ratio distribution function (ERDF) for both unobscured (Type 1) and obscured (Type 2) active galactic nuclei (AGN) using the unprecedented spectroscopic completeness of the BAT AGN Spectroscopic Survey (BASS) data release 2. In addition to a straightforward 1/Vmax approach, we also compute the intrinsic distributions, accounting for sample truncation by employing a forward modeling approach to recover the observed BHMF and ERDF. As previous BHMFs and ERDFs have been robustly determined only for samples of bright, broad-line (Type 1) AGNs and/or quasars, ours is the first directly observationally constrained BHMF and ERDF of Type 2 AGN. We find that after accounting for all observational biases, the intrinsic ERDF of Type 2 AGN is significantly skewed towards lower Eddington ratios than the intrinsic ERDF of Type 1 AGN. This result supports the radiation-regulated unification scenario, in which radiation pressure dictates the geometry of the dusty obscuring structure around an AGN. Calculating the ERDFs in two separate mass bins, we verify that the derived shape is consistent, validating the assumption that the ERDF (shape) is mass independent. We report the local AGN duty cycle as a function of mass and Eddington ratio, by comparing the BASS active BHMF with the local mass function for all SMBH. We also present the log N-log S of Swift-BAT 70-month sources., Comment: Accepted by APJS

Virial black hole (BH) mass (M BH) determination directly involves knowing the broad-line region (BLR) clouds’ velocity distribution, their distance from the central supermassive BH (R BLR), and the virial factor (f). Understanding whether biases arise in M BH estimation with increasing obscuration is possible only by studying a large (N > 100) statistical sample of obscuration-unbiased (hard) X-ray-selected active galactic nuclei (AGNs) in the rest-frame near-infrared (0.8–2.5 μm) since it penetrates deeper into the BLR than the optical. We present a detailed analysis of 65 local Burst Alert Telescope (BAT) selected Seyfert galaxies observed with Magellan/FIRE. Adding these to the near-infrared BAT AGN spectroscopic survey database, we study a total of 314 unique near-infrared spectra. While the FWHMs of Hα and near-infrared broad lines (He i, Paβ, Paα) remain unbiased to either BLR extinction or X-ray obscuration, the Hα broad-line luminosity is suppressed when N H ≳ 1021 cm−2, systematically underestimating M BH by 0.23–0.46 dex. Near-infrared line luminosities should be preferred to Hα until N H < 1022 cm−2, while at higher obscuration a less-biased R BLR proxy should be adopted. We estimate f for Seyfert 1 and 2 using two obscuration-unbiased M BH measurements, i.e., the stellar velocity dispersion and a BH mass prescription based on near-infrared and X-ray, and find that the virial factors do not depend on the redshift or obscuration, but some broad lines show a mild anticorrelation with M BH. Our results show the critical impact obscuration can have on BLR characterization and the importance of the near-infrared and X-rays for a less-biased view of the BLR.

We present the second catalog and data release of optical spectral line measurements and active galactic nucleus (AGN) demographics of the BAT AGN Spectroscopic Survey, which focuses on the Swift-BAT hard X-ray detected AGNs. We use spectra from dedicated campaigns and publicly available archives to investigate spectral properties of most of the AGNs listed in the 70 month Swift-BAT all-sky catalog; specifically, 743 of the 746 unbeamed and unlensed AGNs (99.6%). We find a good correspondence between the optical emission line widths and the hydrogen column density distributions using the X-ray spectra, with a clear dichotomy of AGN types for N (H) = 10(22) cm(-2). Based on optical emission-line diagnostics, we show that 48%-75% of BAT AGNs are classified as Seyfert, depending on the choice of emission lines used in the diagnostics. The fraction of objects with upper limits on line emission varies from 6% to 20%. Roughly 4% of the BAT AGNs have lines too weak to be placed on the most commonly used diagnostic diagram, [O iii]lambda 5007/H beta versus [N ii]lambda 6584/H alpha, despite the high signal-to-noise ratio of their spectra. This value increases to 35% in the [O iii]lambda 5007/[O ii]lambda 3727 diagram, owing to difficulties in line detection. Compared to optically selected narrow-line AGNs in the Sloan Digital Sky Survey, the BAT narrow-line AGNs have a higher rate of reddening/extinction, with H alpha/H beta > 5 (similar to 36%), indicating that hard X-ray selection more effectively detects obscured AGNs from the underlying AGN population. Finally, we present a subpopulation of AGNs that feature complex broad lines (34%, 250/743) or double-peaked narrow emission lines (2%, 17/743)., The Astrophysical Journal. Supplement Series, 261 (1), ISSN:1538-4365, ISSN:0067-0049

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, to develop a diagnostic for nuclear obscuration by examining the relationship between the line-of-sight column densities (N H), the 2–10 keV to 12 μ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 X,Obs./L 12 μm than less obscured AGNs, and we derive expressions relating N H to the L X,Obs./L 12 μm and L 22 μm/L 4.6 μm luminosity ratios, as well as develop diagnostic criteria using these ratios. Our diagnostic regions yield samples that are ≳80% complete and ≳60% pure for AGNs with log(N H/cm−2) ≥ 24, as well as ≳85% pure for AGNs with log ( N H / cm − 2 ) ≳ 23.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 μ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 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.

We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies., Comment: 22 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys., 22 pages, 14 figures, 9 tables. Accepted for publication in ApJ

We investigate the star formation properties of a large sample of ~2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to z~3 in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analysis of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range 9 < log M_{stellar}/M_{\odot} < 12 with uncertainties of ~0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_{stellar} and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift (z>3), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs., 12 pages, 6 figures, Accepted for publication in ApJ

To understand the origin of nuclear (≲100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed subarcsecond resolution Band-6 (211–275 GHz) Atacama Large Millimeter/submillimeter Array data of 98 nearby AGNs (z < 0.05) from the 70 month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling (∼1–200 pc), and spans broad ranges of 14–150 keV luminosity {$40\lt \mathrm{log}[{L}_{14-150}/(\mathrm{erg}\,{{\rm{s}}}^{-1})]\lt 45$}, black hole mass $[5\lt \mathrm{log}({M}_{\mathrm{BH}}/{M}_{\odot })\lt 10$], and Eddington ratio ($-4\lt \mathrm{log}{\lambda }_{\mathrm{Edd}}\lt 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14–150 keV luminosities. Its scatter is ≈0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to NH ∼ 1026 cm−2. While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission: AGN-related dust emission, outflow-driven shocks, and a small-scale (, The Astrophysical Journal, 938 (1), ISSN:0004-637X, ISSN:2041-8213