Fraction Of Active Galactic Nuclei Research Articles

Modelling the galaxy radio continuum from star formation and active galactic nuclei in the Shark semi-analytic model

Abstract We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGN) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission is modelled based on the jet production rate, which depends on the black hole mass, accretion rate and spin, and includes synchrotron self-absorption. Shark reproduces radio luminosity functions (RLFs) at $1.4\, \rm GHz$ and $150\, \rm MHz$ for 0 ≤ z ≤ 4, and scaling relations between radio luminosity, star formation rate and infrared luminosity of galaxies in the local and distant universe in good agreement with observations. The model also reproduces observed number counts of radio sources from 150 MHz to 8.4 GHz to within a factor of two on average, though larger discrepancies are seen at the very bright fluxes at higher frequencies. We use this model to understand how the radio continuum emission from radio-quiet AGNs can affect the measured RLFs of galaxies. We find current methods to exclude AGNs from observational samples result in large fractions of radio-quiet AGNs contaminating the ‘star-forming galaxies’ selection and a brighter end to the resulting RLFs.We investigate how this effects the infrared-radio correlation (IRRC) and show that AGN contamination can lead to evolution of the IRRC with redshift. Without this contamination our model predicts a redshift- and stellar mass-independent IRRC, except at the dwarf-galaxy regime.

Active Galactic Nuclei Selection and Demographics: A New Age with JWST/MIRI

Understanding the coevolution of supermassive black holes and their host systems requires a comprehensive census of active galactic nuclei (AGNs) behavior across a wide range of redshift, luminosity, obscuration level, and galaxy properties. We report significant progress with JWST toward this goal from the Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES). Based on comprehensive spectral energy distribution (SED) analysis of 3273 MIRI-detected sources, we identify 217 AGN candidates over a survey area of ∼34 arcmin2, including a primary sample of 111 AGNs in normal massive galaxies (M * > 109.5 M ☉) at z ∼ 0–4, an extended sample of 86 AGN candidates in low-mass galaxies (M * < 109.5 M ☉), and a high-z sample of 20 AGN candidates at z ∼ 4–8.4. Notably, about 80% of our MIRI-selected AGN candidates are new discoveries despite the extensive pre-JWST AGN searches. Even among the massive galaxies where the previous AGN search is believed to be thorough, 34% of the MIRI AGN identifications are new, highlighting the impact of obscuration on previous selections. By combining our results with the efforts at other wavelengths, we build the most complete AGN sample to date and examine the relative performance of different selection techniques. We find the obscured AGN fraction increases from L AGN,bol ∼ 1010 L ⊙ to 1011 L ⊙ and then drops toward higher luminosity. Additionally, the obscured AGN fraction gradually increases from z ∼ 0 to z ∼ 4 with most high-z AGNs obscured. We discuss how AGN obscuration, intrinsic SED variations, galaxy contamination, survey depth, and selection techniques complicate the construction of a complete AGN sample.

AGNs in massive galaxy clusters: Role of galaxy merging, infalling groups, cluster mass, and dynamical state

Context. There is compelling evidence that active galactic nuclei (AGNs) in high-density regions have undergone a different evolution than their counterparts in the field, indicating that they are strongly affected by their environment. However, we still lack a comprehensive understanding of the dominant mechanisms that trigger the nucleus and the processes that drive the evolution of AGNs in clusters. Aims. To investigate (and possibly disentangle) the various factors that may affect the prevalence of AGNs in cluster galaxies, we selected a sample of 19 thoroughly studied X-ray-selected galaxy clusters from the LoCuSS survey. All these clusters are considered massive, with M500 ≳ 2 × 1014 M⊙, and span a narrow redshift range between z ∼ 0.16 and 0.28. Methods. We divided the cluster surroundings into two concentric annuli with a width of R500 radius. We considered the first annulus as the central cluster region and the second as the outskirts. We further divided the cluster sample based on the presence of infalling X-ray-detected groups, cluster mass, or dynamical state. We determined the AGN fraction in cluster galaxies of the various sub-samples by correlating the X-ray point-like sources selected from the 4XMM DR10 catalogue with the highly complete spectroscopic catalogue of cluster members obtained with Hectospec. We subsequently used the optical spectra to determine the type of nuclear activity and we visually inspected the host morphology for indications of galaxy mergers or other interactions. Results. We found that the X-ray AGN fraction in the outskirts is consistent with the field, but it is significantly lower in cluster centres, in agreement with previous results for massive clusters. We show that these results do not depend on cluster mass, at least within our cluster mass range, nor on the presence of X-ray-detected infalling groups. Furthermore, we did not find any evidence of a spatial correlation between infalling groups and AGNs. Nevertheless, a significant excess of X-ray AGNs is found in the outskirts of relaxed clusters at the 2σ confidence level, compared both to non-relaxed clusters and to the field. Finally, according to the literature, the fraction of broad- to narrow-line AGNs in clusters is roughly consistent with the field. However, broad-line AGNs may be preferably located in cluster centres. In the outskirts, the optical spectra of X-ray AGNs present narrow emission lines or they are dominated by stellar emission. Conclusions. Our results suggest that the mechanisms that trigger AGN activity may vary between cluster centres and the outskirts. Ram pressure can efficiently remove the gas from infalling galaxies, thereby triggering AGN activity in some cases. However, the reduced availability of gas globally diminishes the fraction of AGNs in cluster centers. The surplus of X-ray AGNs identified in the outskirts of relaxed clusters may be attributed to an increased frequency of galaxy mergers, a notion that is further supported by the disturbed morphology observed in several galaxies.

O iv]- and [Ne v]-weak Active Galactic Nuclei Hidden by Compton-thick Material in Late Mergers

We study “buried” active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow-line regions. Employing an indicator of a [O iv] 25.89 or [Ne v] 14.32 μm line to 12 μm AGN luminosity ratio, we find 17 buried AGN candidates that are [O iv]-weak (L [O IV]/L 12,AGN ≤ −3.0) or [Ne v]-weak (L [Ne V]/L 12,AGN ≤ −3.4) among 30 AGNs in local U/LIRGs. For the [O iv]-weak AGNs, we estimate their covering fractions of Compton-thick (CT; N H ≥ 1024 cm−2) material with an X-ray clumpy torus model to be fCT(spec)=0.55±0.19 on average. This value is consistent with the fraction of CT AGNs ( fCT(stat)=53%±12% ) among the [O iv]-weak AGNs in U/LIRGs and much larger than that in Swift/Burst Alert Telescope (BAT) AGNs (23% ± 6%). The fraction of [O iv]-weak AGNs increases from 27−10+13% (early) to 66−12+10% (late mergers). Similar results are obtained with the [Ne v] line. The [O iv]- or [Ne v]-weak AGNs in late mergers show larger N H and Eddington ratios (λ Edd) than those of the Swift/BAT AGNs, and the largest N H is ≳1025 cm−2 at logλEdd∼−1 , close to the effective Eddington limit for CT material. These suggest that (1) the circumnuclear material in buried AGNs is regulated by the radiation force from high-λ Edd AGNs on the CT obscurers, and (2) their dense material with large fCT(spec) (∼0.5 ± 0.1) in U/LIRGs is a likely cause of a unique structure of buried AGNs, whose amount of material may be maintained through merger-induced supply from their host galaxies.

Constraints on supermassive black hole binaries from JWST and NANOGrav

We use the recent statistics of dual active galactic nuclei (AGN) in the James Webb Space Telescope (JWST) data at z ∼ 3.4 to address two aspects of the feedback and evolution scenarios of supermassive black hole binaries. We find that the JWST data provide evidence for the members of a binary BH being ‘lit’ at the same time, rather than independently – a scenario which is consistent with gas-rich mergers being responsible for concurrent AGN activity. This conclusion is supported by the recent NANOGrav Pulsar Timing Array (PTA) measurements, whose upper limits on the stochastic gravitational wave strain amplitude lie below those expected from extrapolating the dual AGN fraction. The results indicate either a ‘stalling’ of the binaries at the separations probed by NANOGrav, or rapid gas-driven inspirals.

The impact of AGN X-ray selection on the AGN halo occupation distribution

The connection between active galactic nuclei (AGN) and their host dark matter halos provides powerful insights into how supermassive black holes (SMBHs) grow and coevolve with their host galaxies. Here we investigate the impact of observational AGN selection on the AGN halo occupation distribution (HOD) by forward-modeling AGN activity into cosmological N-body simulations. By assuming straightforward relationships between the SMBH mass, galaxy mass, and (sub)halo mass, as well as a uniform broken power law distribution of Eddington ratios, we find that luminosity-limited AGN samples result in biased HOD shapes. While AGN defined by an Eddington ratio threshold produce AGN fractions that are flat across halo mass (unbiased by definition), luminosity-limited AGN fractions peak around galaxy-group-sized halo masses and then decrease with increasing halo mass. With higher luminosities, the rise of the AGN fraction starts at higher halo masses, the peak is shifted towards higher halo masses, and the decline at higher halo masses is more rapid. These results are consistent with recent HOD constraints from AGN clustering measurements, which find (1) characteristic halo mass scales of $ Vir $h^ M_ odot and (2) a shallower rise of the number of satellite AGN with increasing halo mass than for the overall galaxy population. Thus the observational biases due to AGN selection can naturally explain the constant, characteristic halo mass scale inferred from large-scale AGN clustering amplitudes over a range of redshifts, as well as the measured inconsistencies between AGN and galaxy HODs. We conclude that AGN selection biases can have significant impacts on the inferred AGN HOD, and can therefore lead to possible misinterpretations of how AGN populate dark matter halos and the AGN-host galaxy connection.

Unveiling Luminous Lyα Emitters at z ≈ 6 through JWST/NIRCam Imaging in the COSMOS Field

We study a sample of 14 spectroscopically confirmed Lyα emitters (LAEs) in the late era of reionization (at redshift z ≈ 6) based on the JWST/NIRCam imaging data set. These LAEs with high Lyα luminosity of L(Lyα) ∼1042.4–1043.4 erg s−1 have been covered by the (ongoing) COSMOS-Web survey over 0.28 deg2 in four NIRCam bands (F115W, F150W, F277W, and F444W). With JWST/NIRCam imaging, we determine the UV continua with M UV ranging from −20.5 to −18.5 mag. The UV slopes have a median value of β ≈ −2.35, and the steepest slopes can reach β < −3. Under the excellent spatial resolution of JWST, we identify three objects in the sample as potential merging/interacting systems. The 14 LAEs (and their components) are compact in morphology, residing substantially below the mass–size relation of high-z galaxies. We further investigate their physical properties, including the stellar mass (M *) and star formation rates (SFRs). Most of the LAEs lie on the SFR–M * main-sequence relation, while two of them, featured as “little red dots,” likely host active galactic nuclei (AGNs), implying a ∼10% AGN fraction. Moreover, we reveal that a new correlation may exist between Lyα equivalent width and the offset between Lyα and UV emission (Δd Lyα ), with a median Δd Lyα ∼ 1 kpc. This could be explained by the Lyα radiative transfer process in both the interstellar medium and circumgalactic medium. The results usher in a new era of detailed analysis on high-z LAEs with the JWST capability.

Active Galactic Nuclei in a Mid-infrared Selected Galaxy Sample at z > 0.13: [Ne v]λ3426 Line Emission as a Benchmark

We present a 24 μm-selected spectroscopic sample z > 0.13 (median 〈z〉 = 0.41) in the Lockman Hole field, consisting of 4035 spectra. Our aim is to identify active galactic nuclei (AGNs) and determine their fraction in this mid-infrared-selected sample. In this work, we use the [Ne v]λ3426 emission line to spectroscopically identify AGNs. Combined with broad-line Type I AGNs selected in our previous study, our sample consists of 887 (∼22%) spectroscopically confirmed AGNs. We perform a stacking analysis on the remaining spectra and find that in various mid-infrared (MIR) wedge-selected AGN candidates, the stacked spectra still show significant [Ne v]λ3426 emission. In contrast, no clear [Ne v]λ3426 signal is detected in non-AGN candidates falling outside the wedges. Assuming a range of AGN MIR spectral energy distribution (SED) slope of −0.3 <α < 0.7, and an average star-forming relation derived from 65 star-forming templates, we develop a robust method to separate the AGN and star-forming contributions to the MIR SEDs using the rest-frame L 12/L 1.6 versus L 4.5/L 1.6 diagram. We separate the objects into bins of L 12 and find that the AGN fraction increases with increasing L 12. We also find that the stacked [Ne v]λ3426 strength scales with L 12. The pure AGN luminosity at 12 μm exhibits a positive correlation with the star formation rates, indicating possible coevolution and common gas supply between the AGNs and their host galaxies. Varying population properties across the redshift range explored contribute to the observed correlation.

Interacting galaxies in the IllustrisTNG simulations – VII: the connection between the most luminous active galactic nuclei and galaxy interactions

ABSTRACT We investigate the connection between the most luminous active galactic nuclei (AGN), galaxy pairs, and post-mergers in the IllustrisTNG simulation. We select galaxy pairs and post-mergers with a mass ratio between 1:10 &amp;lt; μ &amp;lt; 1:1 and a redshift between 0 &amp;lt; z &amp;lt; 1. We compare the incidence of luminous AGN in pairs with matched non-pair controls, finding that AGN with luminosity Lbol &amp;gt; 1044 erg s−1 occur in ∼26 per cent of paired galaxies with a companion within 20 kpc, compared with ∼14 per cent in controls (a fractional excess of ∼2). The enhancement of AGN in galaxy pairs is luminosity dependent, with the highest excess (a factor of ∼6 ± 2 at the closest separations) for AGN with Lbol &amp;gt; 1045 erg s−1. Additionally, pairs exhibit a modest yet statistically significant excess of luminous AGN up to ∼150 kpc in separation. For pairs which merge between 0 &amp;lt; z &amp;lt; 1, AGN fractions are elevated between 1.5 Gyr before and after coalescence, with the highest excess closest to coalescence. Our results indicate that pre-coalescence interactions drive excesses of luminous AGN, but that luminous AGN in galaxy pairs are not ubiquitous. Finally, we investigate what fraction of AGN can be associated with an interaction (either having a companion within 100 kpc or a merger within the last 500 Myr). For AGN with Lbol &amp;gt; 1045 erg s−1, ∼55 per cent are interacting, compared with a 30 per cent interaction fraction in AGN with 1044 &amp;lt; Lbol &amp;lt; 1044.5 erg s−1. Our results support a picture in which interactions play a dominant role in (but are not the sole cause of) triggering the most luminous AGN.

WHaD diagram: Classifying the ionizing source with one single emission line

Context. The usual approach to classify the ionizing source using optical spectroscopy is based on the use of diagnostic diagrams that compare the relative strength of pairs of collisitional metallic lines (e.g., [O III] and [N II]) to recombination hydrogen lines (e.g., Hβ and Hα). Despite it having been accepted as the standard procedure, it presents known problems, including confusion regimes and/or limitations related to the required signal-to-noise (S/N) of the emission lines involved. These problems not only affect our intrinsic understanding of the interstellar medium and its properties, but also the fundamental galaxy properties, such as the star formation rate and the oxygen abundance. This raises key questions related to the fraction of active galactic nuclei and other essential parameters. Aims. We attempt to minimize the problems introduced by the use of these diagrams, in particular, their implementation when the available information is limited due to either the fact that not all lines are available or they do not have the required S/N value. Methods. We explored the existing alternatives in the literature to minimize the confusion among different ionizing sources. We have proposed a new, simple diagram that uses the equivalent width and the velocity dispersion from one single emission line, Hα, to classify the ionizing sources. Results. We used aperture-limited and spatially resolved spectroscopic data from the nearby Universe (z ∼ 0.01) to demonstrate that the new diagram, which we have named WHaD, segregates the different ionizing sources in a more efficient way than earlier procedures. A new set of regions have been defined in this diagram to select among different ionizing sources. Conclusions. The new proposed diagram is well positioned to assist in determining the ionizing source when only Hα is available or when the S/N of the emission lines is too low to obtain reliable fluxes for the weakest emission lines in classical diagnostic diagrams (e.g., Hβ).

Black Hole Mass and Eddington-ratio Distributions of Less-luminous Quasars at z ∼ 4 in the Subaru Hyper Suprime-Cam Wide Field

We investigate the black hole mass function (BHMF) and Eddington-ratio distribution function (ERDF) of broad-line active galactic nuclei (AGNs) at z = 4, based on a sample of 52 quasars with i < 23.2 at 3.50 ≤ z ≤ 4.25 from the Hyper Suprime-Cam Subaru Strategic Program S16A-Wide2 data set, and 1462 quasars with i < 20.2 in the same redshift range from the Sloan Digital Sky Survey data release 7 quasar catalog. Virial black hole (BH) masses of quasars are estimated using the width of the CIV 1549 Å line and the continuum luminosity at 1350 Å. To obtain the intrinsic broad-line AGN BHMF and ERDF, we correct for the incompleteness in the low-mass and/or low-Eddington-ratio ranges caused by the flux-limited selection. The resulting BHMF is constrained down to logMBH/M⊙∼7.5 . In comparison with broad-line AGN BHMFs at z ∼ 2 in the literature, we find that the number density of massive SMBHs peaks at higher redshifts, consistent with the downsizing evolutionary scenario. Additionally, the resulting ERDF shows a negative dependence on BH mass, suggesting more massive SMBHs tend to accrete at lower-Eddington ratios at z = 4. With the derived intrinsic broad-line AGN BHMF, we also evaluate the active fraction of broad-line AGNs among the entire SMBH population at z = 4. The resulting active fraction may suggest a positive dependence on BH mass. Finally, we examine the time evolution of broad-line AGN BHMF between z = 4 and 6 through solving the continuity equation. The results suggest that the broad-line AGN BHMFs at z = 4–6 only show evolution in their normalization, but with no significant changes in their shape.

MaNGA AGN dwarf galaxies (MAD) – I. A new sample of AGNs in dwarf galaxies with spatially-resolved spectroscopy

ABSTRACT The finding of active galactic nuclei (AGNs) in dwarf galaxies has important implications for galaxy evolution and supermassive black hole formation models. Yet, how AGNs in dwarf galaxies form is still debated, in part due to scant demographics. We use the MaNGA survey, comprising ∼10 000 galaxies at z &amp;lt; 0.15, to identify AGN dwarf galaxies using a spaxel-by-spaxel classification in three spatially-resolved emission-line diagnostic diagrams (the [N ii]-, [S ii]-, and [O i]-BPT) and the WHAN diagram. This yields a sample of 664 AGN dwarf galaxies, the largest to date, and an AGN fraction of $\sim 20~{{\ \rm per\ cent}}$ that is significantly larger than that of single-fibre-spectroscopy studies (i.e. $\sim 1~{{\ \rm per\ cent}}$). This can be explained by the lower bolometric luminosity (&amp;lt;1042 erg s−1) and accretion rate (sub-Eddington) of the MaNGA AGN dwarf galaxies. We additionally identify 1176 SF-AGNs (classified as star-forming in the [N ii]-BPT but as AGNs in the [S ii]- and [O i]-BPT), 122 Composite, and 173 LINER sources. The offset between the optical centre of the galaxy and the median position of the AGN spaxels is more than 3 arcsec for ∼62% of the AGNs, suggesting that some could be off-nuclear. We also identify seven new broad-line AGNs with log MBH = 5.0–5.9 M⊙. Our results show how integral-field spectroscopy is a powerful tool for uncovering faint and low-accretion AGNs and better constraining the demographics of AGNs in dwarf galaxies.

Distribution of Hα Emitters in Merging Galaxy Clusters

Studies of star formation in various galaxy cluster mergers have reached apparently contradictory conclusions regarding whether mergers stimulate star formation, quench it, or have no effect. Because the mergers studied span a range of time since pericenter (TSP), it is possible that the apparent effect on star formation is a function of the TSP. We use a sample of 12 bimodal mergers to assess the star formation as a function of TSP. We measure the equivalent width of the Hα emission line in ∼100 member galaxies in each merger, classify galaxies as emitters or nonemitters, and then classify emitters as star-forming galaxies (SFGs) or active galactic nucleus (AGN) based on the [N ii] λ6583 line. We quantify the distribution of SFG and AGN relative to nonemitters along the spatial axis defined by the subcluster separation. The SFG and AGN fractions vary from merger to merger but show no trend with TSP. The spatial distribution of SFG is consistent with that of nonemitters in eight mergers, but show significant avoidance of the system center in the remaining four mergers, including the three with the lowest TSP. If there is a connection between star formation activity and TSP, probing it further will require more precise TSP estimates and more mergers with TSP in the range of 0–400 Myr.

The AGN fraction in high-redshift protocluster candidates selected by Planck and Herschel

ABSTRACT A complete understanding of the mass assembly history of structures in the Universe requires the study of the growth of galaxies and their supermassive black holes (SMBHs) as a function of their local environment over cosmic time. In this context, it is important to quantify the effects that the early stages of galaxy cluster development have on the growth of SMBHs. We used a sample of Herschel/SPIRE sources of ∼ 228 red and compact Planck-selected protocluster (PC) candidates to estimate the active galactic nuclei (AGNs) fraction from a large sample of galaxies within these candidates. We estimate the AGN fraction by using the mid-infrared (mid-IR) photometry provided by the WISE/AllWISE data of ∼650 counterparts at high redshifts. We created an AllWISE mid-IR colour–colour selection using a clustering machine learning algorithm and two WISE colour cuts using the 3.4 $\mu$m (W1), 4.6 $\mu$m (W2), and 12 $\mu$m (W3) passbands, to classify sources as AGN. We also compare the AGN fraction in PCs with that in the field to better understand the influence of the environment on galaxy development. We found an AGN fraction of fAGN = 0.113 ± 0.03 in PC candidates and an AGN fraction of fAGN = 0.095 ± 0.013 in the field. We also selected a subsample of ‘red’ SPIRE subsample with a higher overdensity significance, obtaining fAGN = 0.186 ± 0.044, versus fAGN = 0.037 ± 0.010 of ‘non-red sources’, consistent with higher AGN fractions for denser environments. We conclude that our results point towards a higher AGN fraction in PCs, similar to other studies.

EPOCHS VII: discovery of high-redshift (6.5 &amp;lt; z &amp;lt; 12) AGN candidates in JWST ERO and PEARLS data

ABSTRACT We present an analysis of a sample of robust high-redshift galaxies selected from the ‘blank’ fields of the Prime Extragalactic Areas for Reionization Science (PEARLS) survey and Early Release Observations (ERO) data from JWST with the aim of selecting candidate high-redshift active galactic nuclei (AGN). Sources were identified from this parent sample using a threefold selection procedure, which includes spectral energy distribution (SED) fitting to identify sources that are best fitted by AGN SED templates, a further selection based on the relative performance of AGN and non-AGN models, and finally morphological fitting to identify compact sources of emission, resulting in a purity-oriented procedure. Using this procedure, we identify a sample of nine AGN candidates at 6.5 &amp;lt; z &amp;lt; 12, from which we constrain their physical properties as well as measure a lower bound on the AGN fraction in this redshift range of 5 ± 1 per cent. As this is an extreme lower limit due to our focus on purity and our SEDs being calibrated for unobscured Type 1 AGN, this demonstrates that AGN are perhaps quite common at this early epoch. The rest-frame UV colours of our candidate objects suggest that these systems are potentially candidate obese black hole galaxies (OBGs). We also investigate Chandra and VLA maps of these areas from which we calculate detection limits. Of note is a z = 11.9 candidate source exhibiting an abrupt morphological shift in the reddest band as compared to bluer bands, indicating a potential merger or an unusually strong outflow.

X-ray properties and obscured fraction of AGN in the J1030 Chandra field

The 500ks Chandra ACIS-I observation of the field around the z = 6.31 quasar SDSS J1030+0524 is currently the fifth deepest extragalactic X-ray survey. The rich multi-band coverage of the field allowed an effective identification and redshift determination of the X-ray source counterparts; to date, a catalog of 243 extragalactic X-ray sources with either a spectroscopic or photometric redshift estimate in the range z ≈ 0 − 6 is available over an area of 355 arcmin2. Given its depth and the multi-band information, this catalog is an excellent resource to investigate X-ray spectral properties of distant active galactic nuclei (AGN) and derive the redshift evolution of their obscuration. We performed a thorough X-ray spectral analysis for each object in the sample, and measured its nuclear column density NH and intrinsic (de-absorbed) 2–10 keV rest-frame luminosity, L2 − 10. Whenever possible, we also used the presence of the Fe Kα emission line to improve the photometric redshift estimates. We measured the fractions of AGN hidden by column densities in excess of 1022 and 1023 cm−2 (f22 and f23, respectively) as a function of L2 − 10 and redshift, and corrected for selection effects to recover the intrinsic obscured fractions. At z ∼ 1.2, we found f22 ∼ 0.7 − 0.8 and f23 ∼ 0.5 − 0.6, respectively, in broad agreement with the results from other X-ray surveys. No significant variations in X-ray luminosity were found within the limited luminosity range probed by our sample (log L2 − 10 ∼ 42.8 − 44.3). When focusing on luminous AGN with log L2 − 10 ∼ 44 to maximize the sample completeness up to large cosmological distances, we did not observe any significant change in f22 or f23 over the redshift range z ∼ 0.8 − 3. Nonetheless, the obscured fractions we measure are significantly higher than is seen in the local Universe for objects of comparable intrinsic luminosity, pointing toward an increase in the average AGN obscuration toward early cosmic epochs, as also observed in other X-ray surveys. We finally compared our results with recent analytic models that ascribe the greater obscuration observed in AGN at high redshifts to the dense interstellar medium (ISM) of their hosts. When combined with literature measurements, our results favor a scenario in which the total column density of the ISM and the characteristic surface density of its individual clouds both increase toward early cosmic epochs as NH, ISM∝(1 + z)δ, with δ ∼ 3.3 − 4 and Σc, * ∝ (1 + z)2, respectively.

The co-evolution of strong AGN and central galaxies in different environments

ABSTRACT We exploit a sample of 80 000 Sloan Digital Sky Survey central galaxies to investigate the effect of active galactic nucleus (AGN) feedback on their evolution. We trace the demographics of optically selected AGN (Seyferts) as a function of their internal properties and environment. We find that the preeminence of AGN as the dominant ionizing mechanism increases with stellar mass, overtaking star formation for galaxies with Mstellar ≥ 1010 M⊙. The AGN fraction changes systematically with the galaxies’ star formation activity. Within the blue cloud, this fraction increases as star formation activity declines, reaching a maximum near the green valley (${\sim} 17 \pm 4~{{\ \rm per\ cent}}$), followed by a decrease as the galaxies transition into the red sequence. This systematic trend provides evidence that AGN feedback plays a key role in regulating and suppressing star formation. In general, Seyfert central galaxies achieve an early-type morphology while they still host residual star formation. This suggests that, in all environments, the morphology of Seyfert galaxies evolves from late- to early-type before their star formation is fully quenched. Stellar mass plays an important role in this morphological transformation: while low-mass systems tend to emerge from the green valley with an elliptical morphology (T-Type ∼ −2.5 ± 0.7), their high-mass counterparts maintain a spiral morphology deeper into the red sequence. In high-stellar mass centrals, the fraction of Seyferts increases from early- to late-type galaxies, indicating that AGN feedback may be linked with the morphology and its transformation. Our analysis further suggests that AGN are fuelled by their own host halo gas reservoir, but when in group centrals can also increase their gas reservoir via interactions with satellite galaxies.

Automated detection of double nuclei galaxies using Gothic and the discovery of a large sample of dual AGN

ABSTRACT We present a novel algorithm to detect double nuclei galaxies (DNG) called Gothic (Graph-bOosTed iterated HIll Climbing) – that detects whether a given image of a galaxy has two or more closely separated nuclei. Our aim is to test for the presence of dual/multiple active galactic nuclei (AGN) in galaxies that visually represent a DNG. Although galaxy mergers are common, the detection of dual AGN is rare. Their detection is very important as they help us understand the formation of supermassive black hole (SMBH) binaries, SMBH growth and AGN feedback effects in multiple nuclei systems. There is thus a need for an algorithm to do a systematic survey of existing imaging data for the discovery of DNGs and dual AGNs. We have tested Gothic on an established sample of DNGs with a 100 per cent detection rate and subsequently conducted a blind search of 1 million SDSS DR16 galaxies (with spectroscopic data available) lying in the redshift range of z = 0 to 0.75. From the list of candidate DNGs found, we have detected 159 dual AGNs, of which 2 are triple AGN systems. Our results show that dual AGNs are not common, and triple AGN even rarer. The colour (u–r) magnitude plots of the DNGs indicate that star formation is quenched as the nuclei come closer and as the AGN fraction increases. The quenching is especially prominent for dual/triple AGN galaxies that lie at the extreme end of the red sequence.

A High Fraction of Heavily X-Ray-obscured Active Galactic Nuclei

We present new estimates on the fraction of heavily X-ray-obscured, Compton-thick (CT) active galactic nuclei (AGNs) out to a redshift of z ≤ 0.8. From a sample of 540 AGNs selected by mid-infrared (MIR) properties in observed X-ray survey fields, we forward model the observed-to-intrinsic X-ray luminosity ratio () with a Markov Chain Monte Carlo simulation to estimate the total fraction of CT AGNs (f CT), many of which are missed in typical X-ray observations. We create model N H distributions and convert these to using a set of X-ray spectral models. We probe the posterior distribution of our models to infer the population of X-ray-nondetected sources. From our simulation we estimate a CT fraction of f CT = . We perform an X-ray stacking analysis for sources in Chandra X-ray Observatory fields and find that the expected soft (0.5–2 keV) and hard (2–7 keV) observed fluxes drawn from our model to be within 0.48 and 0.12 dex of our stacked fluxes, respectively. Our results suggests at least 50% of all MIR-selected AGNs, possibly more, are CT (N H ≳ 1024 cm−2), which is in excellent agreement with other recent work using independent methods. This work indicates that the total number of AGNs is higher than can be identified using X-ray observations alone, highlighting the importance of a multiwavelength approach. A high f CT also has implications for black hole (BH) accretion physics and supports models of BH and galaxy coevolution that include periods of heavy obscuration.

Unraveling Joint Evolution of Bars, Star Formation, and Active Galactic Nuclei of Disk Galaxies

We aim to unravel the interplay between bars, star formation (SF), and active galactic nuclei (AGNs) in barred galaxies. To this end, we utilize the SDSS DR12 to select a sample of nearby (0.02 < z < 0.06) disk galaxies that are suitable for bar examination (M r < −20.12 and inclination ≲53°). We identify 3662 barred galaxies and measure the length and axis ratio of each bar. We invent new bar parameters that mitigate the stellar and bulge mass biases and show, for the first time, that the evolution of non-AGN and AGN-hosting barred galaxies should be tracked using different bar parameters: the bar length for non-AGN galaxies and the bar axis ratio for AGN-hosting galaxies. Our analysis confirms that barred galaxies have a higher specific SF rate than unbarred control galaxies. Moreover, we find a positive correlation of bar length with both the SF enhancement and the centrally star-forming galaxy fraction, indicating the interconnectivity of bars and SF through the bar-driven gas inflow. We also find that, while the AGN fraction of barred galaxies is the same as that of the unbarred control sample, galaxies hosting more massive black holes (BHs) have rounder (i.e., higher axis ratio) bars, implying that the bar is not a cause of AGN activity; rather, AGNs appear to regulate bars. Our findings corroborate theoretical predictions that bars in non-AGN galaxies grow in length, and bars in AGN-hosting galaxies become rounder as BHs grow and eventually get destroyed.