ABSTRACT In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $$L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$$. AGN feedback in some simulations prevents the existence of any bright AGN with $$L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($$L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$$) and fainter ($$L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simulations, at any redshift. 
                        more » 
                        « less   
                    
                            
                            Large-scale dual AGN in large-scale cosmological hydrodynamical simulations
                        
                    
    
            ABSTRACT Detecting dual active galactic nuclei (DAGNs) in observations and understanding theoretically which massive black holes (MBHs) compose them and in which galactic and large-scale environment they reside are becoming increasingly important questions as we enter the multimessenger era of MBH astronomy. This paper presents the abundance and properties of DAGN produced in nine large-scale cosmological hydrodynamical simulations. We focus on DAGN powered by AGN with $$L_{\rm bol}\geqslant 10^{43}\, \rm erg\, s^{-1}$$ and belonging to distinct galaxies, i.e. pairs that can be characterized with current and near-future electromagnetic observations. We find that the number density of DAGN separated by a few to 30 proper kpc varies from $$10^{-8}$$ (or none) to $$10^{-3} \, \rm comoving\, Mpc^{3}$$ in the redshift range $$z=0\!-\!7$$. At a given redshift, the densities of the DAGN numbers vary by up to two orders of magnitude from one simulation to another. However, for all simulations, the DAGN peak is in the range $$z=1\!-\!3$$, right before the peak of cosmic star formation or cosmic AGN activity. The corresponding fractions of DAGN (with respect to the total number of AGN) range from 0 per cent to 6 per cent. We find that simulations could produce too few DAGN at $z=0$ (or merge pairs too quickly) compared to current observational constraints while being consistent with preliminary constraints at high redshift ($$z\sim 3$$). Next-generation observatories (e.g. Advanced X-Ray Imaging Satellite [AXIS]) will be of paramount importance to detect DAGN across cosmic times. We predict the detectability of DAGN with future X-ray telescopes and discuss DAGN as progenitors for future Laser Interferometer Space Antenna (LISA) gravitational wave detections. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2108944
- PAR ID:
- 10563832
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 536
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3016-3040
- Size(s):
- p. 3016-3040
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            ABSTRACT We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model, the new model provides a better match to a large number of up-to-date empirical constraints, including: (1) new galaxy and AGN luminosity functions; (2) stellar spectra including binary stars; (3) obscured and unobscured AGN; (4) a measurement of the non-ionizing UV background; (5) measurements of the intergalactic H i and He ii photoionization rates at z ∼ 0−6; (6) the local X-ray background; and (7) improved measurements of the intergalactic opacity. In this model, AGN dominate the H i ionizing background at z ≲ 3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z ∼ 2, the slow evolution of the H i ionization rate inferred from the high-redshift H i Ly α forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of $$\approx 1{{\ \rm per\ cent}}$$ suffices to ionize the intergalactic medium at z = 3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the He ii Ly α forest in hydrodynamic simulations.more » « less
- 
            ABSTRACT We aim to determine the intrinsic far-Infrared (far-IR) emission of X-ray-luminous quasars over cosmic time. Using a 16 deg2 region of the Stripe 82 field surveyed by XMM-Newton and Herschel Space Observatory, we identify 2905 X-ray luminous (LX > 1042 erg/s) active galactic nuclei (AGN) in the range z ≈ 0–3. The IR is necessary to constrain host galaxy properties such as star formation rate (SFR) and gas mass. However, only 10 per cent of our AGN are detected both in the X-ray and IR. Because 90 per cent of the sample is undetected in the far-IR by Herschel, we explore the mean IR emission of these undetected sources by stacking their Herschel/SPIRE images in bins of X-ray luminosity and redshift. We create stacked spectral energy distributions from the optical to the far-IR, and estimate the median SFR, dust mass, stellar mass, and infrared luminosity using a fitting routine. We find that the stacked sources on average have similar SFR/Lbol ratios as IR detected sources. The majority of our sources fall on or above the main sequence line suggesting that X-ray selection alone does not predict the location of a galaxy on the main sequence. We also find that the gas depletion time scales of our AGN are similar to those of dusty star forming galaxies. This suggests that X-ray selected AGN host high star formation and that there are no signs of declining star formation.more » « less
- 
            ABSTRACT The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyse the population of faint AGNs ($$L_{\rm x, 2{-}10 \, keV}\leqslant 10^{42}\, \rm erg\,s^{ -1}$$) in the Illustris, TNG100, EAGLE, and SIMBA cosmological simulations, and find that the properties of their host galaxies vary from one simulation to another. In Illustris and EAGLE, faint AGNs are powered by low-mass BHs located in low-mass star-forming galaxies. In TNG100 and SIMBA, they are mostly associated with more massive BHs in quenched massive galaxies. We model the X-ray binary (XRB) populations of the simulated galaxies, and find that AGNs often dominate the galaxy AGN + XRB hard X-ray luminosity at z > 2, while XRBs dominate in some simulations at z < 2. Whether the AGN or XRB emission dominates in star-forming and quenched galaxies depends on the simulations. These differences in simulations can be used to discriminate between galaxy formation models with future high-resolution X-ray observations. We compare the luminosity of simulated faint AGN host galaxies to observations of stacked galaxies from Chandra. Our comparison indicates that the simulations post-processed with our X-ray modelling tend to overestimate the AGN + XRB X-ray luminosity; luminosity that can be strongly affected by AGN obscuration. Some simulations reveal clear AGN trends as a function of stellar mass (e.g. galaxy luminosity drop in massive galaxies), which are not apparent in the observations.more » « less
- 
            ABSTRACT The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of active galactic nuclei (AGNs), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 luminous and ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGNs that are Compton thick (CT; $$N_{\rm H}\ge 10^{24}\rm \, cm^{-2}$$) peaks at $$74_{-19}^{+14}{{\ \rm per\ cent}}$$ at a late merger stage, prior to coalescence, when the nuclei have projected separations (dsep) of 0.4–6 kpc. A similar peak is also observed in the median NH [$$(1.6\pm 0.5)\times 10^{24}\rm \, cm^{-2}$$]. The vast majority ($$85^{+7}_{-9}{{\ \rm per\ cent}}$$) of the AGNs in the final merger stages (dsep ≲ 10 kpc) are heavily obscured ($$N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$$), and the median NH of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray-selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the $$N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$$ gas almost completely covering the AGN in late mergers. CT AGNs tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity ($$L_{2-10}\lesssim 10^{43}\rm \, erg\, s^{-1}$$) AGNs in U/LIRGs.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
