Unobscured Active Galactic Nuclei Research Articles

Euclid preparation

The Euclid space mission will cover over 14 000 deg2 with two optical and near-infrared spectro-photometric instruments, and is expected to detect around ten million active galactic nuclei (AGN). This unique data set will make a considerable impact on our understanding of galaxy evolution in general, and AGN in particular. For this work we identified the best colour selection criteria for AGN, based only on Euclid photometry or including ancillary photometric observations, such as the data that will be available with the Rubin Legacy Survey of Space and Time (LSST) and observations already available from Spitzer/IRAC. The analysis was performed for unobscured AGN, obscured AGN, and composite (AGN and star-forming) objects. We made use of the spectro-photometric realisations of infrared-selected targets at all-z (SPRITZ) to create mock catalogues mimicking both the Euclid Wide Survey (EWS) and the Euclid Deep Survey (EDS). Using these mock catalogues, we estimated the best colour selection, maximising the harmonic mean (F1) of: (a) completeness, that is, the fraction of AGN correctly selected with respect to the total AGN sample; and (b) purity, that is, the fraction of AGN inside the selection with respect to the selected sample. The selection of unobscured AGN in both Euclid surveys (Wide and Deep) is possible with Euclid photometry alone with F1 = 0.22–0.23 (Wide and Deep), which can increase to F1 = 0.43–0.38 (Wide and Deep) if we limit out study to objects at z > 0.7. Such a selection is improved once the Rubin/LSST filters, that is, a combination of the u, g, r, or z filters, are considered, reaching an F1 score of 0.84 and 0.86 for the EDS and EWS, respectively. The combination of a Euclid colour with the [3.6]−[4.5] colour, which is possible only in the EDS, results in an F1 score of 0.59, improving the results using only Euclid filters, but worse than the selection combining Euclid and LSST colours. The selection of composite (fAGN = 0.05–0.65 at 8–40 μm) and obscured AGN is challenging, with F1 ≤ 0.3 even when including Rubin/LSST or IRAC filters. This is unsurprising since it is driven by the similarities between the broad-band spectral energy distribution of these AGN and star-forming galaxies in the wavelength range 0.3–5 μm.

Systematic collapse of the accretion disc across the supermassive black hole population

ABSTRACT The structure of the accretion flow on to supermassive black holes is not well understood. Standard disc models match to zeroth-order in predicting substantial energy dissipation within optically thick material producing a characteristic strong blue/ultraviolet continuum. However, they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability. Based on stellar mass black holes within our Galaxy, accretion discs should undergo a transition into an X-ray hot, radiatively inefficient flow, below a (mass scaled) luminosity of $\sim 0.02\, L_{\rm {Edd}}$. While this has been seen in limited samples of nearby low-luminosity active galactic nuclei (AGN) and few rare changing-look AGN, it is not at all clear whether this transition is present in the wider AGN population across cosmic time. A key issue is the difficulty in disentangling a change in spectral state from increased dust obscuration and/or host galaxy contamination, effectively drowning out the AGN emission. Here, we use the new eROSITA eFEDS Survey to identify unobscured AGN from their X-ray emission, matched to excellent optical imaging from Subaru’s Hyper Suprime-Cam; allowing the subtraction of the host galaxy contamination. The resulting, uncontaminated, AGN spectra reveal a smooth transition from a strongly disc-dominated state in bright AGN, to the collapse of the disc into an inefficient X-ray plasma in the low-luminosity AGN, with the transition occurring at $\sim 0.02\, L_{\rm {Edd}}$; revealing fundamental aspects of accretion physics in AGN.

Birth of Rapidly Spinning, Overmassive Black Holes in the Early Universe

The James Webb Space Telescope (JWST) has unveiled numerous massive black holes (BHs) in faint, broad-line active galactic nuclei (AGNs). The discovery highlights the presence of dust-reddened AGN populations, referred to as “little red dots (LRDs),” more abundant than X-ray-selected AGNs, which are less influenced by obscuration. This finding indicates that the cosmic growth rate of BHs within this population does not decrease but rather increases at higher redshifts beyond z ∼ 6. The BH accretion rate density deduced from their luminosity function is remarkably higher than that from other AGN surveys in X-ray and infrared bands. To align the cumulative mass density accreted to BHs with the observed BH mass density at z ≃ 4–5, as derived from the integration of the BH mass function, the radiative efficiency must be doubled from the canonical 10% value, achieving significance beyond the >3σ confidence level. This suggests the presence of rapid spins with 96% of the maximum limit among these BHs under the thin-disk approximation, maintained by prolonged mass accretion instead of chaotic accretion with randomly oriented inflows. Moreover, we derive an upper bound for the stellar mass of galaxies hosting these LRDs, ensuring consistency with galaxy formation in the standard cosmological model, where the host stellar mass is limited by the available baryonic reservoir. Our analysis gives a lower bound for the BH-to-galaxy mass ratio that exceeds the typical value known in the nearby universe and aligns with that for JWST-detected unobscured AGNs. Accordingly, we propose a hypothesis that the dense, dust-rich environments within LRDs facilitate the emergence of rapidly spinning and overmassive BH populations during the epoch of reionization. This scenario predicts a potential association between relativistic jets and other high-energy phenomena with overmassive BHs in the early universe.

Accretion properties of X-ray AGN: evidence for radiation-regulated obscuration with redshift-dependent host galaxy contribution

ABSTRACT We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20\lt \log (N_{\rm H}/{\rm cm}^{-2}\lt 22$) and obscured ($22\lt \log (N_{\rm H}/{\rm cm}^{-2}\lt 24$) active galactic nuclei (AGNs) to shed light on the orientation versus evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS, and CDFS extragalactic surveys, we constrain their stellar masses, $M_\star$, intrinsic X-ray luminosities, $L_{\rm X}$, obscuring column densities, $N_{\rm H}$, and specific accretion rates $\lambda \propto L_{\rm X}/M_\star$. By combining these observables within a Bayesian non-parametric approach, we infer, for the first time, the specific accretion rate distribution (SARD) of obscured and unobscured AGN to $z\approx 3$, i.e. the probability of a galaxy with mass $M_\star$ at redshift z hosting an AGN with column density $N_{\rm H}$ and specific accretion rate $\lambda$. Our findings indicate that (1) both obscured and unobscured SARDs share similar shapes, shifting towards higher accretion rates with redshift, (2) unobscured SARDs exhibit a systematic offset towards higher $\lambda$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\log \lambda _{\rm break} \sim -2$ for $z \lt 0.5$, but shifts to higher $\lambda$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $\lambda -N_{\rm H}$ plane increases with redshift. These observations provide compelling evidence for AGN ‘downsizing’ and radiation-regulated nuclear-scale obscuration with an increasing host galaxy contribution towards higher redshifts.

Heavily obscured AGN detection: A radio versus X-ray challenge

Context. In the supermassive black hole (SMBH)-galaxy coevolution scenario, heavily obscured active galactic nuclei (AGN) represent a fundamental phase of SMBH growth during which most of the BH mass is accreted and the scaling relations with the host galaxy are set. Obscured nuclei are thought to constitute a major fraction of the whole AGN population, but their statistics and evolution across cosmic time are still highly uncertain. Therefore, it is pivotal to identify new ways to detect this vast and hidden population of growing SMBHs. A promising way to select heavily obscured AGN is through radio emission, which is largely unaffected by obscuration and can be used as a proxy for nuclear activity. Aims. In this work, we study the AGN radio detection effectiveness in the major deep extragalactic surveys, considering different AGN obscuration levels, redshift, and AGN bolometric luminosities. We particularly focus on comparing their radio and X-ray detectability, making predictions for present and future radio surveys. Methods. We extrapolated the predictions of the AGN population synthesis model of the cosmic X-ray background (CXB) to the radio band, by deriving the 1.4 GHz luminosity functions of unobscured (i.e., with hydrogen column densities log NH < 22), obscured (22 < log NH < 24), and Compton-thick (CTK, log NH > 24) AGN. We then used these functions to forecast the number of detectable AGN based on the area, flux limit, and completeness of a given radio survey and compare it with the AGN number resulting from X-ray predictions. Results. When applied to deep extragalactic fields covered both by radio and X-ray observations, we show that, while X-ray selection is generally more effective in detecting unobscured AGN, the surface density of CTK AGN radio detected is on average ten times larger than the X-ray one, and even greater at high redshifts, considering the current surveys and facilities. Our results suggest that thousands of CTK AGN are already present in current radio catalogs, but most of them escaped any detection in the corresponding X-ray observations. We also present expectations for the number of AGN to be detected by the Square Kilometer Array Observatory (SKAO) in its future deep and wide radio continuum surveys, finding that it will be able to detect more than 2000 AGN at z > 6 and tens of them at z > 10, more than half of which are expected to be CTK.

Euclidpreparation

The statistical distribution and evolution of key properties of active galactic nuclei (AGN), such as their accretion rate, mass, and spin, remains a subject of open debate in astrophysics. The ESAEuclidspace mission, launched on July 1 2023, promises a breakthrough in this field. We create detailed mock catalogues of AGN spectra from the rest-frame near-infrared down to the ultraviolet – including emission lines – to simulate whatEuclidwill observe for both obscured (type 2) and unobscured (type 1) AGN. We concentrate on the red grisms of the NISP instrument, which will be used for the wide-field survey, opening a new window for spectroscopic AGN studies in the near-infrared. We quantify the efficiency in the redshift determination as well as in retrieving the emission line flux of the Hα+[N II] complex, asEuclidis mainly focused on this emission line, given that it is expected to be the brightest one in the probed redshift range. Spectroscopic redshifts are measured for 83% of the simulated AGN in the interval where the Hαis visible (i.e. 0.89 < z < 1.83 at a line flux of > 2 × 10−16erg s−1cm−2, encompassing the peak of AGN activity atz ≃ 1 − 1.5) within the spectral coverage of the red grism. Outside this redshift range, the measurement efficiency decreases significantly. Overall, a spectroscopic redshift iscorrectly determined for about 90% of type 2 AGN down to an emission line flux of roughly 3 × 10−16erg s−1cm−2, and for type 1 AGN down to 8.5 × 10−16erg s−1cm−2. Recovered values for black hole mass show a small offset with respect to the input values by about 10%, but the agreement is good overall. With such a high spectroscopic coverage atz < 2, we will be able to measure AGN demography, scaling relations, and clustering from the epoch of the peak of AGN activity down to the present-day Universe for hundreds of thousands of AGN with homogeneous spectroscopic information.

Test for Echo: X-Ray Reflection Variability in the Seyfert-2 Active Galactic Nucleus NGC 4388

We report on a study of the narrow Fe Kα line and reflection spectrum in the well-known Seyfert-2 active galactic nucleus (AGN), NGC 4388. X-ray spectra summed from two extensive NICER monitoring campaigns, separated by years, show strong evidence of variation in the direct continuum and reflected emission, but only small variations in the obscuring gas. Fits to the spectra from individual NICER observations find a strong, positive correlation between the power-law photon index, Γ, and direct flux that is commonly observed in unobscured AGN. A search for a reverberation lag between the direct and reflected spectra—dominated by the narrow Fe Kα emission line—measures a timescale of t=16.37−0.38+0.46 days, or a characteristic radius of r=1.374−0.032+0.039×10−2 pc =3.4−0.1+0.1×104GM/c2 . Only one cycle of this tentative lag is observed, but it is driven by a particularly sharp drop in the direct continuum that leads to the subsequent disappearance of the otherwise prominent Fe Kα line. Physically motivated fits to high-resolution Chandra spectra of NGC 4388 measure a line production radius of r=2.9−0.7+1.2×104GM/c2 , formally consistent with the tentative lag. The line profile also prefers a Compton-thick reflector, indicating an origin in the disk and/or thick clumps within a wind. We discuss the strengths and weaknesses of our analysis and methods for testing our results in future observations, and we note the potential for X-ray reverberation lags to constrain black hole masses in obscured Seyferts wherein the optical broad line region is not visible.

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar’s brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.

Exploring the nature of UV-bright z ≳ 10 galaxies detected by JWST: star formation, black hole accretion, or a non-universal IMF?

ABSTRACT We use the Cosmic Archaeology Tool (CAT) semi-analytical model to explore the contribution of Population (Pop) III/II stars and active galactic nuclei (AGNs) to the galaxy ultraviolet (UV) luminosity function (LF) evolution at 4 ≤ z ≤ 20. We compare in particular with recent JWST data in order to explore the apparent tension between observations and theoretical models in the number density of bright galaxies at z ≳ 10. The model predicts a star formation history dominated by UV faint (MUV > −18) galaxies, with a Pop III contribution of $\lesssim 10~{{\ \rm per\ cent}}$ ($\lesssim 0.5~{{\ \rm per\ cent}}$) at z ≃ 20 (z ≃ 10). Stars are the primary sources of cosmic reionization, with $5~{{\ \rm per\ cent}}\!-\!10~{{\ \rm per\ cent}}$ of ionizing photons escaping into the intergalatic medium at 5 ≤ z ≤ 10, while the contribution of unobscured AGNs becomes dominant only at z ≲ 5. The predicted stellar and AGN UV LFs reproduce the observational data at 5 ≲ z ≲ 9–10. At higher redshift, CAT predicts a steeper evolution in the faint-end slope (MUV > –18), and a number density of bright galaxies (MUV ≃ −20) consistent with data at z ∼ 10–11, but smaller by 0.8 dex at z ∼ 12–13, and 1.2 dex at z ∼ 14–16, when compared to the values estimated by recent studies. Including the AGN emission does not affect the above findings, as AGNs contribute at most to $\lesssim 10~{{\ \rm per\ cent}}$ of the total UV luminosity at MUV < –19 and z ≳ 10. Interestingly, considering a gradual transition in the stellar initial mass function, modulated by metallicity and redshift as suggested by recent simulations, the model agrees with JWST data at z ∼ 12–13, and the disagreement at z ∼ 14–16 is reduced to 0.5 dex.

The high energy X-ray probe (HEX-P): probing the physics of the X-ray corona in active galactic nuclei

The hard X-ray emission in active galactic nuclei (AGN) and black hole X-ray binaries is thought to be produced by a hot cloud of electrons referred to as the corona. This emission, commonly described by a power law with a high-energy cutoff, is suggestive of Comptonization by thermal electrons. While several hypotheses have been proposed to explain the origin, geometry, and composition of the corona, we still lack a clear understanding of this fundamental component. NuSTAR has been playing a key role improving our knowledge of X-ray coronæ thanks to its unprecedented sensitivity above 10 keV. However, these constraints are limited to bright, nearby sources. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept combining high spatial resolution X-ray imaging and broad spectral coverage (0.2–80 keV) with a sensitivity superior to current facilities. In this paper, we highlight the major role that HEX-P will play in further advancing our insights of X-ray coronæ notably in AGN. We demonstrate how HEX-P will measure key properties and track the temporal evolution of coronæ in unobscured AGN. This will allow us to determine their electron distribution and test the dominant emission mechanisms. Furthermore, we show how HEX-P will accurately estimate the coronal properties of obscured AGN in the local Universe, helping address fundamental questions about AGN unification. In addition, HEX-P will characterize coronæ in a large sample of luminous quasars at cosmological redshifts for the first time and track the evolution of coronæ in transient systems in real time. We also demonstrate how HEX-P will enable estimating the coronal geometry using spectral-timing techniques. HEX-P will thus be essential to understand the evolution and growth of black holes over a broad range of mass, distance, and luminosity, and will help uncover the black holes’ role in shaping the Universe.

JWST Confirms the Nature of CID-42

The galaxy CID-42 (CXOC J100043.1+020637.2) at z = 0.359 has been proposed to contain a promising candidate for a gravitational-wave recoiling supermassive black hole (SMBH), a slingshot SMBH from a triple-SMBH interaction, or a kiloparsec-scale dual active galactic nucleus (AGN). These claims were primarily based on a pair of bright cores separated by ∼0.″5 resolved in optical Hubble Space Telescope (HST) imaging. Existing HST, Chandra, and ground-based imaging and spectroscopy are unable to confirm either scenario. With improved spatial resolution, depth, and IR wavelength coverage, NIRCam multiband imaging from the COSMOS-Web JWST treasury program well resolved the two cores in CID-42, revealing a significant stellar bulge for both cores (with stellar masses of ∼1010 M ⊙ for both). JWST imaging further revealed that only the SE core contains an unobscured AGN point source, based on both image decomposition and spectral energy distribution fitting. There is no evidence for AGN activity in the NW core based on the current data. These new observations unambiguously rule out the gravitational-wave-recoiling and slingshot-SMBH scenarios and establish CID-42 as a low-redshift merging pair of galaxies, with likely only one active AGN in the system. These results demonstrate the unparalleled capabilities of JWST (even with imaging alone) in studying the galactic-scale environment of merging galaxies and SMBHs.

BASS-XL: X-ray variability properties of unobscured active galactic nuclei

ABSTRACT We investigate the X-ray variability properties of Seyfert 1 Galaxies belonging to the BAT AGN Spectroscopic Survey (BASS). The sample includes 151 unobscured (NH < 1022 cm−2) AGNs observed with XMM–Newton for a total exposure time of ∼27 ms, representing the deepest variability study done so far with high signal-to-noise XMM–Newton observations, almost doubling the number of observations analysed in previous works. We constrain the relation between the normalized excess variance and the 2–10 keV AGN luminosities, black hole masses, and Eddington ratios. We find a highly significant correlation between $\sigma _{\rm NXS}^2$ and MBH , with a scatter of ∼0.85 dex. For sources with high L2–10 this correlation has a lower normalization, confirming that more luminous (higher mass) AGNs show less variability. We explored the $\sigma _{\rm NXS}^2$ versus MBH relation for the sub-sample of sources with MBH estimated via the ‘reverberation mapping’ technique, finding a tighter anticorrelation, with a scatter of ∼0.65 dex. We examine how the $\sigma _{\rm NXS}^2$ changes with energy by studying the relation between the variability in the hard (3–10 keV) and the soft (0.2–1 keV)/medium (1–3 keV) energy bands, finding that the spectral components dominating the hard energy band are more variable than the spectral components dominating in softer energy bands, on time-scales shorter than 10 ks.

AT 2021loi: A Bowen Fluorescence Flare with a Rebrightening Episode Occurring in a Previously Known AGN

The optical-ultraviolet transient AT 2021loi is located at the center of its host galaxy. Its spectral features identify it as a member of the Bowen fluorescence flare (BFF) class. The first member of this class was considered to be related to a tidal disruption event, but enhanced accretion onto an already active supermassive black hole was suggested as an alternative explanation. Having occurred in a previously known unobscured active galactic nucleus, AT 2021loi strengthens the latter interpretation. Its light curve is similar to those of previous BFFs, showing a rebrightening approximately 1 yr after the main peak (which was not explicitly identified but might be the case in all previous BFFs). An emission feature around 4680 Å, seen in the preflare spectrum, strengthens by a factor of ∼2 around the optical peak of the flare and is clearly seen as a double-peaked feature then, suggesting a blend of N iii λ4640 with He ii λ4686 as its origin. The appearance of O iii λ3133 and possible N iii λλ4097, 4103 (blended with Hδ) during the flare further support a Bowen fluorescence classification. Here we present ZTF, ATLAS, Keck, Las Cumbres Observatory, NEOWISE-R, Swift AMI, and Very Large Array observations of AT 2021loi, making it one of the best-observed BFFs to date. It thus provides some clarity on the nature of BFFs but also further demonstrates the diversity of nuclear transients.

A Census of WISE-selected Dual and Offset AGNs Across the Sky: New Constraints on Merger-driven Triggering of Obscured AGNs

Pairs of galaxies hosting active galactic nuclei (AGNs) are powerful probes of merger-driven supermassive black hole (SMBH) growth as they can resolve individual AGNs and trace mergers over a large range of physical separations. To exploit this on a large scale for the first time for both obscured and unobscured AGNs, we use the photometric redshifts of AGNs selected by the Wide-field Infrared Survey Explorer (WISE) to find probabilistic pairs (<100 kpc separations) across the sky, along with a comparison sample of inactive galaxy pairs. Our final sample of integrated pair probabilities yields 198 AGN–AGN pairs (dual AGNs) and 2767 AGN–galaxy pairs (offset AGNs) with uniformly measured AGN and host galaxy physical properties. We find the fraction of galaxy pairs hosting WISE AGNs is dominated by offset AGNs and is significantly elevated above that of inactive galaxies for large host stellar masses. We show how the AGN merger fraction directly increases with AGN extinction for both offset and dual AGNs, with up to ∼40% of heavily obscured AGNs found in galaxy pairs. Elevated AGN merger fractions coincide with increased host specific star formation rates that suggest merger-driven coevolution of galaxies and SMBHs. Among dual AGNs, the most rapid SMBH growth may occur within the less-massive galaxy. Relative to stochastic mechanisms, mergers produce an excess of AGNs at increasingly smaller separations, especially for obscured AGNs (up to a factor of ∼5), and are augmented by correlated triggering. Finally, this excess is stronger than for lower-luminosity optically selected AGNs, regardless of the AGN obscuration level.

A Candidate for the Least-massive Black Hole in the First 1.1 Billion Years of the Universe

We report a candidate of a low-luminosity active galactic nucleus (AGN) at z = 5 that was selected from the first near-infrared images of the JWST CEERS project. This source, named CEERS-AGN-z5-1 at absolute 1450 Å magnitude M 1450 = −19.5 ± 0.3, was found via a visual selection of compact sources from a catalog of Lyman break galaxies at z > 4, taking advantage of the superb spatial resolution of the JWST/NIRCam images. The 20 photometric data available from CFHT, Hubble Space Telescope, Spitzer, and JWST suggest that the continuum shape of this source is reminiscent of that for an unobscured AGN, and there is a clear color excess in the filters where the redshifted Hβ+[O iii] and Hα are covered. The estimated line luminosity is L Hβ+[O III] = 1043.0 erg s−1 and L Hα = 1042.9 erg s−1 with the corresponding rest-frame equivalent width EWHβ+[O III] = 1100 Å and EWHα = 1600 Å, respectively. Our spectral energy distribution fitting analysis favors the scenario that this object is either a strong broad-line emitter or even a super-Eddington accreting black hole (BH), although a possibility of an extremely young galaxy with moderate dust attenuation is not completely ruled out. The bolometric luminosity, L bol = 2.5 ± 0.3 × 1044 erg s−1, is consistent with those of z < 0.35 broad-line AGNs with M BH ∼ 106 M ⊙ accreting at the Eddington limit. This new AGN population in the first 1.1 billion years of the universe may close the gap between the observed BH mass range at high redshift and that of BH seeds. Spectroscopic confirmation is awaited to secure the redshift and its AGN nature.

NuSTAR Observations of a Heavily X-Ray-obscured AGN in the Dwarf Galaxy J144013+024744

We present a multiwavelength analysis of the dwarf Seyfert 2 galaxy J144013+024744, a candidate obscured active galactic nucleus (AGN) thought to be powered by an intermediate-mass black hole (IMBH, M • ≈ 104−106 M ⊙) of mass M • ∼ 105.2 M ⊙. To study its X-ray properties, we targeted J144013+024744 with NuSTAR for ≈100 ks. The X-ray spectrum was fitted with an absorbed power law, Pexmon, and a physical model (RXTorus). A Bayesian X-ray analysis was performed to estimate the posteriors. The phenomenological and the physical models suggest the AGN to be heavily obscured by a column density of N H = (3.4–7.0) × 1023 cm−2. In particular, the RXTorus model with a subsolar metallicity suggests the obscuring column to be almost Compton-thick. We compared the 2–10 keV intrinsic X-ray luminosity with the inferred X-ray luminosities based on empirical scaling relations for unobscured AGNs using L [Oiv] 25.89 μm, L [Oiii] λ5007, and L 6μm and found that the high-excitation [Oiv] line provides a better estimate of the intrinsic 2–10 keV X-ray luminosity ( erg s−1). Our results suggest that J144013+024744 is the first type 2 dwarf galaxy that shows X-ray spectroscopic evidence for obscuration. The column density that we estimated is among the highest measured to date for IMBH-powered AGNs, implying that a typical AGN torus geometry might extend to the low-mass end. This work has implications for constraining the BH occupation fraction in dwarf galaxies using X-ray observations.

Compton-thick AGN in the NuSTAR Era. IX. A Joint NuSTAR and XMM-Newton Analysis of Four Local AGN

We present the results of the broadband X-ray spectral analysis of simultaneous NuSTAR and XMM-Newton observations of four nearby Compton-thick active galactic nuclei (AGN) candidates selected from the Swift Burst Alert Telescope 150 month catalog. This work is part of a larger effort to identify and characterize all Compton-thick (N H ≥ 1024 cm−2) AGN in the local universe (z ≤ 0.05). We used three physically motivated models—MYTorus, borus02, and UXClumpy—to fit and characterize these sources. Of the four candidates analyzed, 2MASX J02051994-0233055 was found to be an unobscured (N H < 1022 cm−2) AGN, 2MASX J04075215-6116126 and IC 2227 to be Compton-thin (1022 cm−2 < N H < 1024 cm−2) AGN, and one, ESO 362−8, was confirmed to be a Compton-thick AGN. Additionally, every source was found to have a statistically significant difference between their line-of-sight and average torus hydrogen column density, further supporting the idea that the obscuring material in AGN is inhomogeneous. Furthermore, half of the sources in our sample (2MASX J02051994-0233055 and 2MASX J04075215-6116126) exhibited significant luminosity variation in the last decade, suggesting that this might be a common feature of AGN.

Investigating non-linear and stochastic hard X-ray variability of active galactic nuclei using recurrence analysis

ABSTRACT We present results of recurrence analysis of 46 active galactic nuclei (AGNs) using light curves from the 157-month catalogue of the Swift Burst Alert Telescope (BAT) in the 14–150 keV band. We generate recurrence plots and compute recurrence plot metrics for each object. We use the surrogate data method to compare all derived recurrence-based quantities to three sets of stochastic light curves with identical power spectrum, flux distribution, or both, in order to determine the presence of determinism, non-linearity, entropy, and non-stationarity. We compare these quantities with known physical characteristics of each system, such as black hole mass, Eddington ratio, and bolometric luminosity, radio loudness, obscuration, and spectroscopic type. We find that almost all AGNs in this sample exhibit substantial higher order modes of variability than is contained in the power spectrum, with approximately half exhibiting non-linear or non-stationary behaviour. We find that Type 2 AGNs are more likely to contain deterministic variability than Type 1 AGNs, while the same distinction is not found between obscured and unobscured AGNs. The complexity of variability among Type 1 AGNs is anticorrelated with Eddington ratio, while no relationship is found among Type 2 AGNs. The connections between the recurrence properties and AGN class suggest that hard X-ray emission is a probe of distinct accretion processes among classes of AGNs, which supports interpretations of changing-look AGNs and challenges the traditional unification model that classifies AGNs only on viewing angle.

Measurement of AGN dust extinction based on the near-infrared flux variability of WISE data

ABSTRACT We present the measurement of the line-of-sight extinction of the dusty torus for a large number of obscured active galactic nuclei (AGNs) based on the reddening of the colour of the variable flux component in near-infrared (NIR) wavelengths. We collected long-term monitoring data by Wide-field Infrared Survey Explorer (WISE) for 513 local AGNs catalogued by the Swift/BAT AGN Spectroscopic Survey (BASS) and found that the multi-epoch NIR flux data in two different bands (WISE W1 and W2) are tightly correlated for more than 90 per cent of the targets. The flux variation gradient (FVG) in the W1 and W2 bands was derived by applying linear regression analysis, and we reported that those for unobscured AGNs fall in a relatively narrow range, whereas those for obscured AGNs are distributed in a redder and broader range. The AGN’s line-of-sight dust extinction (AV) is calculated using the amount of the reddening in the FVG and is compared with the neutral hydrogen column density (NH) of the BASS catalogue. We found that the NH/AV ratios of obscured AGNs are greater than those of the Galactic diffuse interstellar medium (ISM) and are distributed with a large scatter by at most two orders of magnitude. Furthermore, we found that the lower envelope of the NH/AV of obscured AGNs is comparable to the Galactic diffuse ISM. These properties of the NH/AV can be explained by increase in the NH attributed to the dust-free gas clouds covering the line of sight in the broad-line region.

Extreme X-Ray Reflection in the Nucleus of the Seyfert Galaxy NGC 5033

NGC 5033 is an intriguing Seyfert galaxy because its subclassification may change with time, and because optical and submillimeter observations find that the massive black hole does not sit at the dynamical center of the galaxy, pointing to a past merger. We obtained a new optical spectrum of NGC 5033 using the 200″ Hale telescope at Palomar that clearly reveals a broad Hβ line (FWHM = 5400 ± 300 km s−1). This signals a clear view of the optical broad line region and requires Seyfert-1.5 designation. Some spectra obtained in the past suggest a Seyfert-1.9 classification, potentially signaling a variable or “changing-look” geometry. Our analysis of a 2019 Chandra spectrum of the massive black hole reveals very little obscuration, also suggesting a clean view of the central engine. However, the narrow Fe Kα emission line is measured to have an equivalent width of EW eV. This value is extremely high compared to typical values in unobscured active galactic nuclei . Indeed, the line is persistently strong in NGC 5033: the line equivalent width in a 2002 XMM-Newton snapshot is EW eV, similar to the EW eV equivalent width measured using ASCA in 1999. These results can likely be explained through a combination of an unusually high covering factor for reflection, and fluxes that are seen out of phase owing to light travel times. We examine the possibility that NGC 5033 may strengthen evidence for the X-ray Baldwin effect.