Bright Active Galactic Nuclei Research Articles

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.

PSF calibration of DAMPE for gamma-ray observations

The DArk Matter Particle Explorer (DAMPE) is dedicated to exploring critical scientific domains including the indirect detection of dark matter, cosmic ray physics, and gamma ray astronomy. This study introduces a novel method for calibrating the Point Spread Function (PSF) of DAMPE, specifically designed to enhance the accuracy of gamma-ray observations. By leveraging data from regions near pulsars and bright Active Galactic Nuclei (AGNs), we have refined the PSF calibration process, resulting in an improved angular resolution that closely matches our observational data. This advancement significantly boosts the precision of gamma-ray detection by DAMPE, thereby contributing to its mission objectives in dark matter detection and gamma ray astronomy.

X-Ray Bright Active Galactic Nuclei in Local Dwarf Galaxies: Insights from eROSITA

Although supermassive black holes (SMBHs) reside in the heart of virtually every massive galaxy, it remains debated whether dwarf galaxies commonly host SMBHs. Because low-mass galaxies may retain memory of the assembly history of their black holes (BHs), probing the BH occupation fraction of local dwarf galaxies might offer insights into the growth and seeding mechanisms of the first BHs. In this work, we exploit the Western half of the eROSITA all-sky survey (covering 20,000 deg2) and compile a catalog of accreting SMBHs in local (D < 200 Mpc) dwarf galaxies. Cleaning our sample from X-ray background sources, X-ray binaries, and ultraluminous X-ray sources, we identify 74 active galactic nucleus (AGN)–dwarf galaxy pairs. Using this large and uniform sample, we derive the luminosity function of the dwarf galaxy AGN, fitting it with a power-law function and obtaining dN/dLX=(15.9±2.2)×LX−1.63±0.05 . Measuring the offset between the dwarf galaxies' centroids and the X-ray sources, we find that ≈50% of the AGN are likely off-nuclear, in agreement with theoretical predictions. We compare the BH-to-stellar mass relation of our sample with the local and high-redshift relations, finding that our sources better adhere to the former, which suggests that local AGN across different mass scales undergo similar growth histories. Finally, we compare our sources with semianalytical models: while our sample’s shallowness prevents distinguishing between different seeding models, we find that the data favor models that keep SMBHs in dwarf galaxies active at a moderate rate, motivating model improvement by comparison to AGN in the dwarf galaxy regime.

Spectroastrometry and Reverberation Mapping of Active Galactic Nuclei. I. The Hβ Broad-line Region Structure and Black Hole Masses of Five Quasars

We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric (SA) observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among which Mrk 1239 and PDS 456 are for the first time spectroscopically monitored. We obtain multiyear monitoring data and perform multicomponent spectral decomposition to extract the broad Hβ profiles. We detect significant time lags between the Hβ and continuum variations, generally obeying the previously established BLR size–luminosity relation. Velocity-resolved Hβ time lags illustrate diverse, possibly evolving, BLR kinematics. We further measure the Hβ line widths from mean and rms spectra and the resulting virial products show good consistency among different observing seasons. Adopting a unity virial factor and the FWHM of the broad Hβ line from the mean spectra as a measure of velocity, the obtained black hole masses averaged over seasons are logM•/M⊙=8.02−0.14+0.09 , 6.92−0.12+0.12 , 8.01−0.25+0.16 , 7.44−0.14+0.13 , and 8.59−0.11+0.07 for the five objects, respectively. Black hole mass estimations using other line width measures are also reported. For objects with previous RM campaigns, our mass estimates are in agreement with earlier results. In a companion paper, we will employ BLR dynamical modeling to directly infer the black hole mass and thereby determine the virial factors.

Unified properties of supermassive black hole winds in radio-quiet and radio-loud AGN

ABSTRACT Powerful supermassive black hole (SMBH) winds in the form of ultra-fast outflows (UFOs) are detected in the X-ray spectra of several active galactic nuclei (AGNs) seemingly independently of their radio classification between radio quiet (RQ) and radio loud (RL). In this work, we explore the physical parameters of SMBH winds through a uniform analysis of a sample of X-ray bright RQ and RL AGN. We explored several correlations between different wind parameters and with respect to the AGN bolometric and Eddington luminosities. Our analysis shows that SMBH winds are not only a common trait of both AGN classes but also that they are most likely produced by the same physical mechanism. Consequently, we find that SMBH winds do not follow the radio-loudness dichotomy seen in jets. On average, a comparable amount of material accreted by the SMBH is ejected through such winds. The average wind power corresponds to about 3 per cent of the Eddington luminosity, confirming that they can drive AGN feedback. Moreover, the most energetic outflows are found in the most luminous sources. We find a possible positive correlation of the wind energetics, renormalized to the Eddington limit, with respect to $\lambda _{\mathrm{ Edd}}$, consistent with the correlation found with bolometric luminosity. We also observe a possible positive correlation between the energetics of the outflow and the X-ray radio-loudness parameter. In general, these results suggest an underlying relation between the acceleration mechanisms of accretion disc winds and jets.

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ∼0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

Proper motion of the radio jets in two blazars at redshift above 3

ABSTRACT There is still a limited number of high-redshift (z &amp;gt; 3) active galactic nuclei (AGNs) whose jet kinematics have been studied with very long baseline interferometry (VLBI). Without a dedicated proper motion survey, regularly conducted astrometric VLBI observations of bright radio-emitting AGN with sensitive arrays can be utilized to follow changes in the jets, by means of high-resolution imaging and brightness distribution modelling. Here, we present a first-time VLBI jet kinematic study of NVSS J080518 + 614423 (z = 3.033) and NVSS J165844 − 073918 (z = 3.742), two flat-spectrum radio quasars that display milliarcsecond-scale jet morphology. Archival astrometric observations carried out mainly with the Very Long Baseline Array, supplemented by recent data taken with the European VLBI Network, allowed us to monitor changes in their radio structure in the 7.6−8.6 GHz frequency band, covering almost two decades. By identifying individual jet components at each epoch, we were able to determine the apparent proper motion for multiple features in both sources. Apparent superluminal motions range $(1-14)\, c$, and are found to be consistent with studies of other high-redshift AGN targets. Using the physical parameters derived from the brightness distribution modelling, we estimate the Doppler-boosting factors (δ ≈ 11.2 and δ ≈ 2.7), the Lorentz factors (Γ ≈ 7.4 and Γ ≈ 36.6), and the jet viewing angles (θ ≈ 4.4° and θ ≈ 8.0°), for NVSS J080518 + 614423 and NVSS J165844 − 073918, respectively. The data revealed a stationary jet component with negligible apparent proper motion in NVSS J165844 − 073918.

On joint analysing &lt;i&gt;XMM-NuSTAR&lt;/i&gt; spectra of active galactic nuclei

A recently released &lt;i&gt;XMM-Newton&lt;/i&gt; note revealed a significant calibration issue between nuclear spectroscopic telescope array (&lt;i&gt;NuSTAR)&lt;/i&gt; and &lt;i&gt;XMM-Newton&lt;/i&gt; European Photon Imaging Camera (EPIC) and provided an empirical correction to the EPIC effective area. To quantify the bias caused by the calibration issue in the joint analysis of &lt;i&gt;XMM-NuSTAR&lt;/i&gt; spectra and verify the effectiveness of the correction, in this work, we perform joint-fitting of the &lt;i&gt;NuSTAR&lt;/i&gt; and EPIC-pn spectra for a large sample of 104 observation pairs of 44 X-ray bright active galactic nuclei (AGN). The spectra were extracted after requiring perfect simultaneity between the &lt;i&gt;XMM-Newton&lt;/i&gt; and &lt;i&gt;NuSTAR&lt;/i&gt; exposures (merging Good Time Intervals, GTIs from two missions) to avoid bias due to the rapid spectral variability of the AGN. Before the correction, the EPIC-pn spectra are systematically harder than the corresponding &lt;i&gt;NuSTAR&lt;/i&gt; spectra by &lt;inline-formula&gt;&lt;tex-math id="M1"&gt;\begin{document}$\Delta \varGamma \sim 0.1 $\end{document}&lt;/tex-math&gt;&lt;alternatives&gt;&lt;graphic specific-use="online" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="JUSTC-2023-0160_M1.jpg"/&gt;&lt;graphic specific-use="print" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="JUSTC-2023-0160_M1.png"/&gt;&lt;/alternatives&gt;&lt;/inline-formula&gt;, subsequently yielding significantly underestimated cutoff energy &lt;i&gt;E&lt;/i&gt;&lt;sub&gt;cut&lt;/sub&gt; and the strength of reflection component &lt;i&gt;R&lt;/i&gt; when performing joint-fitting. We confirm that the correction is highly effective and can commendably erase the discrepancy in best-fit &lt;i&gt;Γ&lt;/i&gt;, &lt;i&gt;E&lt;/i&gt;&lt;sub&gt;cut&lt;/sub&gt;, and &lt;i&gt;R&lt;/i&gt;. We thus urge the community to apply the correction when joint-fitting &lt;i&gt;XMM&lt;/i&gt;-&lt;i&gt;NuSTAR&lt;/i&gt; spectra, but note that the correction is limited to 3–12 keV and therefore not applicable when the soft X-ray band data are included. Besides, we show that as merging GTIs from two missions would cause severe loss of &lt;i&gt;NuSTAR&lt;/i&gt; net exposure time, in many cases, joint-fitting yields no advantage compared with utilizing &lt;i&gt;NuSTAR&lt;/i&gt; data alone. We finally present a technical note on filtering periods of high background flares for &lt;i&gt;XMM-Newton&lt;/i&gt; EPIC-pn exposures in the Small Window mode.

Uchuu–ν2GC galaxies and AGN: cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST

ABSTRACT Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift (z ≳ 6) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and JWST, will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological N-body simulation with 2.1 trillion particles in a volume of 25.7 Gpc3, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu–ν2GC catalogue, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu–ν2GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at z ∼ 6 − 7, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST).

Physics Beyond the Standard Model with Future X-Ray Observatories: Projected Constraints on Very-light Axion-like Particles with Athena and AXIS

Axion-like particles (ALPs) are well-motivated extensions of the Standard Model of Particle Physics and a generic prediction of some string theories. X-ray observations of bright active galactic nuclei (AGNs) hosted by rich clusters of galaxies are excellent probes of very-light ALPs, with masses . We evaluate the potential of future X-ray observatories, particularly Athena and the proposed AXIS, to constrain ALPs via observations of cluster-hosted AGNs, taking NGC 1275 in the Perseus cluster as our exemplar. Assuming perfect knowledge of the instrument calibration, we show that a modest exposure (200 ks) of NGC 1275 by Athena permits us to exclude all photon–ALP couplings g aγ > 6.3 × 10−14 GeV−1 at the 95% confidence level, as previously shown by Conlon et al., representing a factor of 10 improvement over current limits. We then proceed to assess the impact of realistic calibration uncertainties on the Athena projection by applying a standard Cash likelihood procedure, showing the projected constraints on g aγ weaken by a factor of 10 (back to the current most sensitive constraints). However, we show how the use of a deep neural network can disentangle the energy-dependent features induced by instrumental miscalibration and those induced by photon–ALP mixing, allowing us to recover most of the sensitivity to the ALP physics. In our explicit demonstration, the machine learning applied allows us to exclude g aγ > 2.0 × 10−13 GeV−1, complementing the projected constraints of next-generation ALP dark matter birefringent cavity searches for very-light ALPs. Finally, we show that a 200 ks AXIS/on-axis observation of NGC 1275 will tighten the current best constraints on very-light ALPs by a factor of 3.

X-Ray Spectroscopy of Interstellar Carbon: Evidence for Scattering by Carbon-bearing Material in the Spectrum of 1ES 1553+113

Molecules and particles make up ∼40%–70% of carbon in the interstellar medium, yet the exact chemical structure of these constituents remains unknown. We present carbon K-shell absorption spectroscopy of the Galactic interstellar medium obtained with the Low Energy Transmission Grating Spectrometer on board the Chandra Observatory that directly addresses this question. We probe several lines of sight, using bright active galactic nuclei as backlighters. We make our measurements differentially with respect to the bright source Mrk 421, in order to take the significant carbon K absorption in the instrument into account. In the spectrum of the blazar 1ES 1553+113 we find evidence for a novel feature: strong extinction on the low-energy side of the neutral C 1s−2p resonance, which is indicative of scattering by graphite particles. We find evidence for characteristic particle radii of order 0.1–0.15 μm. If this explanation for the feature is correct, limits on the mass of the available carbon along the line of sight may imply that the grains are partially aligned, and the X-rays from the source may have intrinsic polarization.

Molecular gas content and high excitation of a massive main-sequence galaxy at z = 3

We present new CO (J = 5 − 4 and 7 − 6) and [C I] (3P2−3P1 and 3P1−3P0) emission line observations of the star-forming galaxy D49 at the massive end of the main sequence at z = 3. We incorporate previous CO (J = 3 − 2) and optical-to-millimetre continuum observations to fit its spectral energy distribution. Our results hint at high-J CO luminosities exceeding the expected location on the empirical correlations with the infrared luminosity. [CI] emission fully consistent with the literature trends is found. We do not retrieve any signatures of a bright active galactic nucleus that could boost the J = 5 − 4, 7 − 6 lines in either the infrared or X-ray bands, but warm photon-dominated regions, shocks, or turbulence could in principle do so. We suggest that mechanical heating could be a favourable mechanism able to enhance the gas emission at fixed infrared luminosity in D49 and other main-sequence star-forming galaxies at high redshift, but further investigation is necessary to confirm this explanation. We derive molecular gas masses from dust, CO, and [C I] that all agree within the uncertainties. Given its high star formation rate ∼500 M⊙ yr−1 and stellar mass &gt; 1011.5 M⊙, the short depletion timescale of &lt; 0.3 Gyr might indicate that D49 is experiencing its last growth spurt and will soon transit to quiescence.

Quasi-periodic oscillations in the γ-ray light curves of bright active galactic nuclei

Context. The detection of quasi-periodic oscillations (QPOs) in the light curves of active galactic nuclei (AGNs) can provide insights into the physics of the super-massive black holes (SMBHs) powering these systems and could represent a signature of the existence of SMBH binaries, setting fundamental constraints on SMBH evolution in the Universe. Aims. The identification of long-term QPOs, characterized by periods on the order of several months to years, is particularly challenging and can only be achieved via all-sky monitoring instruments that provide unbiased, continuous light curves of astrophysical objects. The Fermi-LAT satellite, thanks to its monitoring observing strategy, is an ideal instrument for such a goal. Here, we aim to identify QPOs in the γ-ray light curves of the brightest AGNs within the Fermi-LAT catalog. Methods. We analyzed the light curves of the 35 brightest Fermi-LAT AGNs, including data from the beginning of the Fermi mission (August 2008) to April 2021, with energies from 100 MeV to 300 GeV. Two time binnings were investigated: 7 and 30 days. The search for quasi-periodic features was then performed using the continuous wavelet transform. The significance of the result was tested via Monte Carlo simulations of artificial light curves with the same power spectral density and probability distribution function as the original light curves. The significances were then corrected for the look-elsewhere effect and provided as post-trials. Results. We identified 24 quasars with candidate QPOs. Several of our candidates coincide with previous claims in the literature, namely: PKS 0537−441, S5 0716+714, Mrk 421, B2 1520+31, and PKS 2247−131. All our candidates are transient. The most significant multi-year QPO, with a period of about 1100 days, was observed in the quasar S5 1044+71. It is reported here for the first time.

Eddington accreting black holes in the epoch of reionization

ABSTRACT The evolution of the luminosity function (LF) of active galactic nuclei (AGNs) at redshift $z {\,\, \gtrsim \,\,}5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at z &amp;gt; 4 have to be revised upwards, and spark new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central supermassive black holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGNs shine at their Eddington limit: We derive the expected evolution as a function of the ‘duty cycle’ (fdc), i.e. the fraction of lifetime that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values (fdc ≃ 0.1) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} {\,\, \lesssim \,\,}0.05$) are required for AGNs to be the dominant population responsible for hydrogen reionization in the early Universe. We then show that this low-fdc evolution can be reconciled with the current constraints on helium reionization, although it implies a relatively large number of inactive SMBHs at $z{\,\, \gtrsim \,\,}5$, in tension with SMBH growth models based on heavy seeding.

Spectroastrometric Reverberation Mapping of Broad-line Regions

Spectroastrometry measures source astrometry as a function of wavelength/velocity. Reverberations of spectroastrometric signals naturally arise in broad-line regions (BLRs) of active galactic nuclei (AGNs) as a result of the continuum variations that drive responses of the broad emission lines with time delays. Such signals provide a new diagnostic for mapping BLR kinematics and geometry, complementary to the traditional intensity reverberation mapping (RM) technique. We present a generic mathematical formalism for spectroastrometric RM and show that under realistic parameters of a phenomenological BLR model, the spectroastrometric reverberation signals vary on a level of several to tens of microarcseconds, depending on the BLR size, continuum variability, and angular-size distance. We also derive the analytical expressions of spectroastrometric RM for an inclined ring-like BLR. We develop a Bayesian framework with a sophisticated Monte Carlo sampling technique to analyze spectroastrometric data and infer the BLR properties, including the central black hole mass and angular-size distance. We demonstrate the potential of spectroastrometric RM in spatially resolving BLR kinematics and geometry through a suite of simulation tests. The application to realistic observation data of 3C 273 obtains tentative, but enlightening results, reinforcing the practical feasibility of conducting spectroastrometric RM experiments on bright AGNs with the operating Very Large Telescope Interferometer as well as possibly with the planned next-generation 30 m class telescopes.

Polar Dust Emission in Quasar IR SEDs and Its Correlation with Narrow-line Regions

Polar dust has been found to play an important role in the mid-infrared emission of nearby Seyfert nuclei. If and how often polar dust exists among the quasar population is unknown due to the lack of spatially resolved observations. In this Letter, we report correlations between the prominence of active galactic nucleus (AGN) forbidden line emission (commonly associated with the narrow-line region) and the dust mid-IR energy output among the archetypal Palomar-Green quasar sample and other bright type-1 AGNs drawn from the Sloan Digital Sky Survey, Spitzer, and Wide-field Infrared Survey Explorer (WISE) archives. The AGN mid-IR color differences traced by WISE W2 (∼4.6 μm)−W3 (∼12 μm) and W2 (∼4.6 μm)−W4 (∼22 μm), and near-IR to mid-IR spectral energy distributions (SEDs) constrained with the Two Micron All Sky Survey, WISE, and Spitzer data have clear trends with the relative strength of the forbidden line regions traced by the optical [O iii] and mid-IR [O iv] emission lines. These observations indicate that, where the lines are strong, a large fraction of the AGN emission at λ ≳ 5 μm comes from dust in the forbidden line regions. We find that the widely quoted universal AGN template is a result of averaging quasar SEDs with different levels of polar dust emission above the torus output and that the typical intrinsic IR SED of compact torus dust emission alone falls with increasing wavelength past 5 μm (in νF ν ). In addition, the association of polar dust with the forbidden lines suggests an alternative to the receding torus hypothesis for the decrease in infrared output with increasing AGN luminosity.

On the Origin of the Strong Optical Variability of Emission-line Galaxies

Emission-line galaxies (ELGs) are crucial for understanding the formation and evolution of galaxies, while little is known about their variability. Here we report on the optical variability of a sample of ELGs selected in the COSMOS field, which has narrowband observations in two epochs separated by ≳12 yr. This sample was observed with the Suprime-Cam (SC) and Hyper Suprime-Cam (HSC) on the Subaru telescope in NB816 and bands, respectively. After carefully removing the wing effect of a narrowband filter, we check the optical variability in a sample of 181 spectroscopically confirmed ELGs. We find that 0 (0/68) Hα emitters, 11.9% (5/42) [O iii] emitters, and 0 (0/71) [O ii] emitters show significant variability () in the two-epoch narrowband observations. We investigate the presence of active galactic nuclei (AGN) in this variable ELG (var-ELG) sample with three methods, including X-ray luminosity, mid-infrared activity, and radio excess. We find zero bright AGN in this var-ELG sample but cannot rule out the contribution from faint AGN. We find that supernovae explosions (SNe) could also dominate the variability of the var-ELG sample. The merger morphology shown in the HST/F814W images of the entire var-ELG sample is in agreement with the enhancement of star formation, i.e., the SNe activity.

On the Coevolution of the AGN and Star-forming Galaxy Ultraviolet Luminosity Functions at 3 < z < 9

Studies of the rest-frame ultraviolet (UV) luminosity functions (LFs) typically treat star-forming galaxies and active galactic nuclei (AGNs) separately. However, modern ground-based surveys now probe volumes large enough to discover AGNs at depths sensitive enough for fainter galaxies, bridging these two populations. Using these observations as constraints, we present a methodology to empirically jointly model the evolution of the rest-UV LFs at z = 3–9. Our critical assumptions are that both populations have LFs well described by double power laws modified to allow for a flattening at the faint-end, and that all LF parameters evolve smoothly with redshift. This methodology provides a good fit to the observations and makes predictions to volume densities not yet observed, finding that the volume density of bright (M UV = −28) AGNs rises by ∼105 from z = 9 to z = 3, while bright (M UV = −21) star-forming galaxies rise by only ∼102 across the same epoch. The observed bright-end flattening of the z = 9 LF is unlikely to be due to AGN activity, and rather is due to a shallowing of the bright-end slope, implying a reduction of feedback in bright galaxies at early times. The intrinsic ionizing emissivity is dominated by star-forming galaxies at z > 3, even after applying a notional escape fraction. We find decent agreement between our AGN LFs and predictions based on different black hole seeding models, though all models underpredict the observed abundance of bright AGNs. We show that the wide-area surveys of the upcoming Euclid and Roman observatories should be capable of discovering AGNs to z ∼ 8.

Detecting clusters of galaxies and active galactic nuclei in an eROSITA all-sky survey digital twin

Context. The extended ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spectrum-Roentgen-Gamma (SRG) observatory is revolutionizing X-ray astronomy. The mission provides unprecedented samples of active galactic nuclei (AGN) and clusters of galaxies, with the potential of studying astrophysical properties of X-ray sources and measuring cosmological parameters using X-ray-selected samples with higher precision than ever before. Aims. We aim to study the detection, and the selection of AGN and clusters of galaxies in the first eROSITA all-sky survey, and to characterize the properties of the source catalog. Methods. We produced a half-sky simulation at the depth of the first eROSITA survey (eRASS1), by combining models that truthfully represent the population of clusters and AGN. In total, we simulated 1 116 758 clusters and 225 583 320 AGN. We ran the standard eROSITA detection algorithm, optimized for extragalactic sources. We matched the input and the source catalogs with a photon-based matching algorithm. Results. We perfectly recovered the bright AGN and clusters. We detected half of the simulated AGN with flux larger than 2 × 10−14 erg s−1 cm−2 as point sources and half of the simulated clusters with flux larger than 3 × 10−13 erg s−1 cm−2 as extended sources in the 0.5–2.0 keV band. We quantified the detection performance in terms of completeness, false detection rate, and contamination. We studied the population in the source catalog according to multiple cuts of source detection and extension likelihood. We find that the latter is suitable for removing contamination, and the former is very efficient in minimizing the false detection rate. We find that the detection of clusters of galaxies is mainly driven by flux and exposure time. It additionally depends on secondary effects, such as the size of the clusters on the sky plane and their dynamical state. The cool core bias mostly affects faint clusters classified as point sources, while its impact on the extent-selected sample is small. We measured the fraction of the area covered by our simulation as a function of limiting flux. We measured the X-ray luminosity of the detected clusters and find that it is compatible with the simulated values. Conclusions. We discuss how to best build samples of galaxy clusters for cosmological purposes, accounting for the nonuniform depth of eROSITA. This simulation provides a digital twin of the real eRASS1.

Localizing narrow Fe Kαemission within bright AGN

Context.The 6.4 keV Fe Kαemission line is a ubiquitous feature in X-ray spectra of active galactic nuclei (AGN), and its properties track the interaction between the variable primary X-ray continuum and the surrounding structure from which it arises.Aims.We clarify the nature and origin of the narrow Fe Kαemission using X-ray spectral, timing, and imaging constraints, plus possible correlations to AGN and host galaxy properties, for 38 bright nearby AGN (z &lt; 0.5) from the Burst Alert Telescope AGN Spectroscopic Survey.Methods.ModelingChandraandXMM-Newtonspectra, we computed line full-width half-maxima (FWHMs) and constructed Fe Kαline and 2–10 keV continuum light curves. The FWHM provides one estimate of the Fe Kαemitting region size,RFeKα, assuming virial motion. A second estimate comes from comparing the degree of correlation between the variability of the continuum and line-only light curves, compared to simulated light curves. Finally, we extractedChandraradial profiles to place upper limits onRFeKα.Results.For 90% (21/24) of AGN with FWHM measurements,RFeKαis smaller than the fiducial dust sublimation radius,Rsub. From timing analysis, 37 and 18 AGN show significant continuum and Fe Kαvariability, respectively. Despite a wide range of variability properties, the constraints on the Fe Kαphoton reprocessor size independently confirm thatRFeKαis smaller thanRsubin 83% of AGN. Finally, the imaging analysis yields loose upper limits for all but two sources; notably, the Circinus Galaxy and NGC 1068 show significant but subdominant extended Fe Kαemission out to ∼100 and ∼800 pc, respectively.Conclusions.Based on independent constraints, we conclude that the majority of the narrow Fe Kαemission in typical AGN predominantly arises from regions smaller than and presumably insideRsub, and thus it is associated either with the outer broad line region or outer accretion disk. However, the large diversity of continuum and narrow Fe Kαvariability properties are not easily accommodated by a universal scenario.