How HST∕WFC3 and JWST Can Measure Galaxy Assembly and AGN Growth

Recent observations and simulations have challenged the long-held paradigm that mergers are the dominant mechanism driving the growth of both galaxies and supermassive black holes (SMBH), in favour of non-merger (secular) processes. In this pilot study of merger-free SMBH and galaxy growth, we use Keck Cosmic Web Imager spectral observations to examine four low-redshift (0.043 < z < 0.073) disc-dominated ‘bulgeless’ galaxies hosting luminous active galactic nucleus (AGN), assumed to be merger-free. We detect blueshifted broadened [O iii] emission from outflows in all four sources, which the $\mathrm{[O\, \small {III}]}$/$\rm {H}\beta$ ratios reveal are ionized by the AGN. We calculate outflow rates in the range $0.12{-}0.7~\rm {M}_{\odot }~\rm {yr}^{-1}$, with velocities of $675{-}1710~\rm {km}~\rm {s}^{-1}$, large radial extents of $0.6{-}2.4~\rm {kpc}$, and SMBH accretion rates of $0.02{-}0.07~\rm {M}_{\odot }~\rm {yr}^{-1}$. We find that the outflow rates, kinematics, and energy injection rates are typical of the wider population of low-redshift AGN, and have velocities exceeding the galaxy escape velocity by a factor of ∼30, suggesting that these outflows will have a substantial impact through AGN feedback. Therefore, if both merger-driven and non-merger-driven SMBH growth lead to co-evolution, this suggests that co-evolution is regulated by feedback in both scenarios. Simulations find that bars and spiral arms can drive inflows to galactic centers at rates an order of magnitude larger than the combined SMBH accretion and outflow rates of our four targets. This work therefore provides further evidence that non-merger processes are sufficient to fuel SMBH growth and AGN outflows in disc galaxies.

Context. Merger events can trigger gas accretion onto supermassive black holes (SMBHs) located at the centre of galaxies and form close pairs of active galactic nuclei (AGNs). The fraction of AGNs in pairs offers critical insights into the dynamics of galaxy interactions, SMBH growth, and their co-evolution with host galaxies. However, the identification of dual AGNs is difficult, as it requires high-quality spatial and spectral data; hence, very few pairs have been found in the distant Universe so far. Aims. This study is aimed at providing a first observational estimate of the fraction of dual AGNs at 2 < z < 6 by analysing a sample of 16 AGNs observed with the JWST Near-InfraRed Spectrograph (NIRSpec) in integral field mode, as part of the GA-NIFS survey. For two AGNs in our sample, we also incorporated archival VLT/MUSE data to expand the search area. Methods. We searched for nearby companion galaxies and emission-line sources within the ∼20 × 20 kpc field of view of the NIRSpec data cubes, extending up to ∼50 kpc using the MUSE data cubes. We analysed the spectra of such emitters to determine their physical and kinematic properties. Results. We report the serendipitous discovery of a triple AGN system and four dual AGNs (two of which had been considered as candidates), with projected separations in the range 3−28 kpc. The results of this study more than double the number of known multiple AGNs at z > 3 at these separations. Their AGN classification is mainly based on standard optical emission line flux ratios, as observed with JWST/NIRSpec, and complemented with additional multi-wavelength diagnostics. The identification of these 3−5 multiple AGNs out of the 16 AGN systems in the GA-NIFS survey (i.e. ∼20−30%) suggests they might be more common than previously thought from other observational campaigns. Moreover, our inferred fraction of dual AGN moderately exceeds predictions from cosmological simulations that mimic our observational criteria (∼10%). Conclusions. This work highlights the exceptional capabilities of NIRSpec for detecting distant dual AGNs, prompting new investigations to constrain their fraction across cosmic time.

Galaxy pairs with separations of only a few kiloparsecs represent important stages in the merger-driven growth of supermassive black holes (SMBHs). However, such mergers are difficult to identify observationally due to the correspondingly small angular scales. In Paper I we presented a method of finding candidate kiloparsec-scale galaxy mergers that is leveraged on the selection of X-ray sources spatially offset from the centers of host galaxies. In this paper we analyze new Hubble Space Telescope (HST) WFC3 imaging for six of these sources to search for signatures of galaxy mergers. The HST imaging reveals that four of the six systems are on-going galaxy mergers with separations of 1.2−6.6 kpc (offset active galactic nuclei—AGNs). The nature of the remaining two spatially offset X-ray sources is ambiguous and may be associated with super-Eddington accretion in X-ray binaries. The ability of this sample to probe small galaxy separations and minor mergers makes it uniquely suited for testing the role of galaxy mergers for AGN triggering. We find that galaxy mergers with only one AGN are predominantly minor mergers with mass ratios similar to the overall population of galaxy mergers. By comparison, galaxy mergers with two AGNs are biased toward major mergers and larger nuclear gas masses. Finally, we find that the level of SMBH accretion increases toward smaller mass ratios (major mergers). This result suggests the mass ratio effects not only the frequency of AGN triggering but also the rate of SMBH growth in mergers.

We present a new technique for observationally identifying galaxy mergers spectroscopically rather than through host galaxy imaging. Our technique exploits the dynamics of supermassive black holes (SMBHs) powering active galactic nuclei (AGNs) in merger-remnant galaxies. Because structure in the universe is built up through galaxy mergers and nearly all galaxies host a central SMBH, some galaxies should possess two SMBHs near their centers as the result of a recent merger. These SMBHs spiral to the center of the resultant merger-remnant galaxy, and one or both of the SMBHs may power AGNs. Using the DEEP2 Galaxy Redshift Survey, we have examined 1881 red galaxies, of which 91 exhibit [O III] and Hbeta emission lines indicative of Seyfert 2 activity. Of these, 32 AGNs have [O III] emission-line redshifts significantly different from the redshifts of the host galaxies' stars, corresponding to velocity offsets of ~50 km/s to ~300 km/s. Two of these AGNs exhibit double-peaked [O III] emission lines, while the remaining 30 AGNs each exhibit a single set of velocity-offset [O III] emission lines. After exploring a variety of physical models for these velocity offsets, we argue that the most likely explanation is inspiralling SMBHs in merger-remnant galaxies. Based on this interpretation, we find that roughly half of the red galaxies hosting AGNs are also merger remnants, which implies that mergers may trigger AGN activity in red galaxies. The AGN velocity offsets we find imply a merger fraction of ~30% and a merger rate of ~3 mergers/Gyr for red galaxies at redshifts 0.34 < z < 0.82.

Context. The large-scale environments of active galactic nuclei (AGNs) reveal important information on the growth and evolution of supermassive black holes (SMBHs). Previous AGN clustering measurements using two-point correlation functions have hinted that AGNs with massive black holes preferentially reside in denser cosmic regions than AGNs with less massive SMBHs. At the same time, little to no dependence on the accretion rate has been found; however, the significance of such trends has been limited. Aims. Here, we present k th-nearest-neighbor ( k NN) statistics of 2MASS galaxies around AGNs from the Swift/BAT AGN Spectroscopic survey. These statistics have been shown to contribute additional higher order clustering information on the cosmic density field. Methods. By calculating the distances to the nearest seven galaxy neighbors in angular separation to each AGN within two redshift ranges (0.01 &lt; z &lt; 0.03 and 0.03 &lt; z &lt; 0.06), we compared their cumulative distribution functions to that of a randomly distributed sample to demonstrate the sensitivity of this method to the clustering of AGNs. We also split the AGNs into bins of bolometric luminosity, black hole mass, and Eddington ratio (while controlling for redshift) to search for trends between k NN statistics and fundamental AGN properties. Results. We find that AGNs with massive SMBHs have significantly closer neighbors than AGNs with less massive SMBHs (at the 99.98% confidence level), especially in our lower redshift range. We find less significant trends with luminosity or Eddington ratio. By comparing our results to empirical SMBH-galaxy-halo models implemented in N-body simulations, we show that small-scale k NN trends with black hole mass may go beyond stellar mass dependencies. Conclusions. This suggests that massive SMBHs in the local universe reside in more massive dark matter halos and denser regions of the cosmic web, which may indicate that environment is important for the growth of SMBHs, bolstering prior conclusions based on correlation functions.

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.

We present a comprehensive synthesis model for the active galactic nuclei (AGN) evolution and the growth of supermassive black hole (SMBH) in the Universe. We assume that black holes accrete in just three distinct physical states, or 'modes': at low Eddington ratio, only a radiatively inefficient, kinetically dominated mode is allowed [low kinetic (LK)]; at high Eddington ratio, instead, AGN may display both a purely radiative [radio quiet, high radiative (HR)] and a kinetic [radio loud, high kinetic (HK)] mode. We solve the continuity equation for the black hole mass function using the locally determined one as a boundary condition, and the hard X-ray luminosity function as tracer of the AGN growth rate distribution, supplemented with a luminosity-dependent bolometric correction and an absorbing column distribution. Differently from most previous semi-analytic and numerical models for black hole growth, we do not assume any specific distribution of Eddington ratios, rather we determine it empirically by coupling the mass and luminosity functions and a set of fundamental relations between observables in the three accretion modes. SMBH always show a very broad accretion rate distribution, and we discuss the profound consequences of this fact for our understanding of observed AGN fractions in galaxies, as well as for the empirical determination of SMBH mass functions with large surveys. We confirm previous results and clearly demonstrate that, at least for z≤ 1.5, SMBH mass function evolves antihierarchically, i.e. the most massive holes grew earlier and faster than less massive ones. For the first time, we find hints of a reversal of such a downsizing behaviour at redshifts above the peak of the black hole accretion rate density (z ≈2). We also derive tight constraints on the (mass-weighted) average radiative efficiency of AGN: under the simplifying assumption that the mass density of both high redshift (z ∼5) and 'wandering' black holes ejected from galactic nuclei after merger events are negligible compared to the local mass density, we find that 0.065 0.07, where ξ 0 is the local SMBH mass density in units of 4.3 x 10 5 M ⊙ Mpc -3 . We trace the cosmological evolution of the kinetic luminosity function of AGN, and find that the overall efficiency of SMBH in converting accreted rest mass energy into kinetic power, e kin , ranges between e kin ≃ (3-5) x 10 -3 , depending on the choice of the radio core luminosity function. Such a 'kinetic efficiency' varies however strongly with SMBH mass and redshift, being maximal for very massive holes at late times, as required for the AGN feedback by many galaxy formation models in cosmological contexts.

Context. Supermassive black holes (SMBHs) grow primarily through gas accretion, observed as active galactic nuclei (AGNs). While mergers can drive luminous AGN episodes, secular processes may fuel a substantial portion of cosmic black hole growth. Whether these mechanisms can sustain high black hole accretion rates remains uncertain. Aims. The aim of this study is to identify the secular mechanism driving high SMBH accretion rates, by targeting a galaxy with a moderately massive SMBH, high central gas densities, accretion rates of a few percent of the Eddington limit, and gas kinematics resolved close to black-hole-dominated scales. Methods. A blind search led to the identification of NGC 4593, which is representative of the AGN population driving BH mass density growth since z = 1. Combining HST imaging, VLT/MUSE spectroscopy, and ALMA imaging, we resolve molecular and ionised gas kinematics close to the sphere of influence of the SMBH. Results. A prominent single-arm (“ m = 1”) molecular gas spiral with log M mol /M ⊙ = 8.1 ± 0.3 extends from 1.3 kpc down to the SMBH’s sphere of influence (1.7 −0.2 +0.5 pc). Star formation in the spiral is inefficient (SFR = 4.9 × 10 −2 M ⊙ /yr, ⟨ t dep ⟩=3.9 ± 0.6 Gyr), whereas the gas inflow rate exceeds the SFR by two orders of magnitude and is sufficient to sustain the current SMBH accretion rate for ≥35 Myr, enabling ∼10% SMBH growth. A young, lopsided stellar component (log M ⋆ /M ⊙ = 7.5 − 9.3) exerts torques on the molecular gas, likely driving the gas inflow. This young stellar component may serve as both a cause and a product of sustained gas funnelling towards the SMBH. Conclusions. These findings provide evidence for the sustained secular m = 1 feeding mode at high SMBH accretion rates, linking kpc-scale gas dynamics to the black hole’s sphere of influence. This mechanism, consistent with simulation predictions, may represent a key contributor to SMBH growth in luminous AGNs since cosmic noon.

Presence of supermassive black holes (SMBHs) with mass (106 − 109) M⊙ at z = 10 has been recently revealed by James Webb Space Telescope (JWST) observations. In this study we generate seeds for the above range of SMBHs in various background cosmologies. We consider cosmic timescales required for black hole growth provided by three general relativistic cosmological models (ΛCDM, ωCDM and Dynamical Dark Energy(DDE) and the braneworld cosmology. The growth of SMBHs is studied through Eddington limited and super-Eddington accretion, where the accretion starts at z=30. It is found that growth of SMBHs by z=10 within Eddington limited accretion is possible through massive seeds (M ≥ 104 M⊙) in all cosmologies. Super Eddington accretion onto spinning black holes with mass of few tens of solar masses can result in SMBHs by z=10 in all cosmologies. The viable cosmologies considered here are found to be unable to strongly distinguish between the seed black hole masses. The seeds generated in this work are assumed to be of primordial origin in order to satisfy the criteria of formation of high redshift massive galaxies. The fraction of primordial black holes (PBHs) contributing to dark matter (fPBH) and their corresponding number densities for the mass range (105–108) M⊙ are calculated in both seed effect and Poisson effect. In seed effect, PBHs of mass ≥107 M⊙ contributes ≤10−2 to the dark matter fraction. The evolution of gas mass inside a PBH seeded dark matter halo is studied. The ratio of black hole to stellar mass is also evaluated for star formation efficiency in the range (0.1-1) and found to be (10−3 − 1) for MBH = 108 M⊙ and (10−2 − 10) for MBH = 109 M⊙.

Intermediate-mass black holes (IMBHs; M BH &lt; 2∗10 5 M ⊙ ) in galaxy centers are cruciel for painting a coherent picture of theformation and growth of supermassive black holes (SMBHs). Using Big Data analysis, we identified 305 IMBH candidatesfor IMBH and 1623 candidates of ‘light-weight’ SMBHs (2 ∗ 10 5 M odot &lt; M BH &lt; 10 6 M ⊙ ). For 35 host galaxies from thiscombined sample with the X-ray-confirmed active galactic nuclei (AGN) we collected and analyzed optical spectroscopicobservations. These data show that bulge stellar velocity dispersions (σ ∗ ) lie in the range of 24. . . 118 km/s and do notfollow the correlation with M BH established by larger SMBHs indicating that in the 10 5 − 10 6 M ⊙ range the accretion is theprevailing BH growth channel.

Context. Exploring supermassive black hole (SMBH) populations in protoclusters offers valuable insights into how environment affects SMBH growth. However, research on active galactic nuclei (AGN) within these areas is still limited by the small number of protoclusters known at high redshift and by the availability of associated deep X-ray observations. Aims. In order to understand how different environments affect AGN triggering and growth at high redshift, we investigated the X-ray AGN population in the field of the MUSE Quasar Nebula 01 (MQN01) protocluster at z ≃ 3.25. This field is known for hosting one of the largest Ly α nebulae and overdensities of UV continuum selected and sub-millimetre galaxies found so far at this redshift. Methods. We conducted an ultra-deep Chandra X-ray survey (634 ks) observation of the MQN01 field and produced a comparative analyses of the properties of the X-ray AGN detected in MQN01 against those observed in other selected protoclusters, such as Spiderweb and SSA22. Results. By combining the X-ray observations with deep MUSE and ALMA data of the same field, we identified six X-ray AGN within a volume of 16 comoving Mpc2 and ±1000 km s−1, corresponding to an X-ray AGN overdensity of δ ≈ 1000. This overdensity increases at the bright end (log(L2#x2212;10 keV/erg s−1) ≳ 44.5), exceeding what was observed in the Spiderweb and SSA22 within similar volumes. The AGN fraction measured in MQN01 is significantly higher (&gt; 20%) than in the field and increases with stellar masses, reaching a value of 100% for log(M*/M⊙) &gt; 10.5. Lastly, we observe that the average specific accretion rate (λsBHAR) for SMBH populations in MQN01 is higher than in the field and other protoclusters, generally increasing as one moves toward the centre of the overdense structures. Conclusions. Our results, especially the large fraction of highly accreting SMBHs in the inner regions of the MQN01 overdensity, suggest that protocluster environments offer ideal physical conditions for SMBH triggering and growth.

The presence of dual active galactic nuclei (AGNs) on scales of a few tens of kiloparsecs can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. This study uses the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (L bol > 1043 erg s−1) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at z ≤ 2. This is our underlying population. We then applied the luminosity threshold of L bol > 1043 erg s−1 to the neighboring SMBHs thereby identifying dual and multiple AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. We also find that the number of dual AGNs with separations of 0.5–4 kpc is twice the number of duals with separations of 4–30 kpc. All dual AGNs in our sample resulted from major mergers. Compared to single AGNs, duals have a lower black hole-to-halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate, and gas mass, are generally consistent with those of single AGN halos, albeit tending toward the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

We suggest a new picture of supermassive black hole (SMBH) growth in galaxy centers. Momentum-driven feedback from an accreting hole gives significant orbital energy but little angular momentum to the surrounding gas. Once central accretion drops, the feedback weakens and swept-up gas falls back towards the SMBH on near-parabolic orbits. These intersect near the black hole with partially opposed specific angular momenta, causing further infall and ultimately the formation of a small-scale accretion disk. The feeding rates into the disk typically exceed Eddington by factors of a few, growing the hole on the Salpeter timescale and stimulating further feedback. Natural consequences of this picture include (i) the formation and maintenance of a roughly toroidal distribution of obscuring matter near the hole; (ii) random orientations of successive accretion disk episodes; (iii) the possibility of rapid SMBH growth; (iv) tidal disruption of stars and close binaries formed from infalling gas, resulting in visible flares and ejection of hypervelocity stars; (v) super-solar abundances of the matter accreting on to the SMBH; and (vi) a lower central dark-matter density, and hence annihilation signal, than adiabatic SMBH growth implies. We also suggest a simple sub-grid recipe for implementing this process in numerical simulations.

Despite substantial effort, the precise physical processes that lead to the growth of super-massive black holes in the centers of galaxies are still not well understood. These phases of black hole growth are thought to be of key importance in understanding galaxy evolution. Forthcoming missions such as eROSITA, HETDEX, eBOSS, BigBOSS, LSST, and Pan-STARRS will compile by far the largest ever Active Galactic Nuclei (AGNs) catalogs which will allow us to measure the spatial distribution of AGNs in the universe with unprecedented accuracy. For the first time, AGN clustering measurements will reach a level of precision that will not only allow for an alternative approach to answering open questions in AGN and galaxy co-evolution but will open a new frontier, allowing us to precisely determine cosmological parameters. This paper reviews large-scale clustering measurements of broad line AGNs. We summarize how clustering is measured and which constraints can be derived from AGN clustering measurements, we discuss recent developments, and we briefly describe future projects that will deliver extremely large AGN samples which will enable AGN clustering measurements of unprecedented accuracy. In order to maximize the scientific return on the research fields of AGN and galaxy evolution and cosmology, we advise that the community develops a full understanding of the systematic uncertainties which will, in contrast to today’s measurement, be the dominant source of uncertainty.

We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent.