Search

Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint active galactic nuclei (AGNs). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction (OF) at z  = 0 is not inherited from the simulation seeding modelling. The ability of BHs and their host galaxies to power an AGN depends on BH and galaxy subgrid modelling. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus can be used to differentiate galaxy formation models. AGN fractions at z  = 0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in the masses of simulated and observed BHs, BH OF affected by numerical choices, and an unknown fraction of obscured AGN. Our work advocates for more thorough comparisons with observations to improve the modelling of cosmological simulations, and our understanding of BH and galaxy physics in the low-mass regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time. [ABSTRACT FROM AUTHOR]

The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosted in one galaxy, after the galaxy merger has happened. We find that the relation between MBH and galaxy mass is similar to that of general AGN population and we compare the properties of dual AGN also with a control sample, discussing differences and similarities in masses and Eddington ratios. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. Between 30 and 80 per cent of dual AGN with separations between 4 and 30 kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50 kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive haloes, and mass loss of satellites shapes the galaxy–halo relation. [ABSTRACT FROM AUTHOR]

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M BH − M ⋆ relation and the BH mass function. All simulations predict tight M BH − M ⋆ relations, and struggle to produce BHs of |$M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$| in galaxies of |$M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$|⁠. While the time evolution of the mean M BH − M ⋆ relation is mild (⁠|$\rm \Delta M_{\rm BH}\leqslant 1\, dex$| for 0 |$\leqslant z \leqslant$| 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for |$M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$|⁠. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with |$M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$| in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (⁠|$\log _{10}\, f_{\rm Edd}\geqslant -2$|⁠) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations. [ABSTRACT FROM AUTHOR]

High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environments do not lead to a consensus on a possible excess of galaxies, the future missions JWST, WFIRST , and Euclid will provide new insights on quasar environments, and will substantially increase the number of study-cases. We are in a crucial period, where we need to both understand the current observations and predict how upcoming missions will improve our understanding of BH environments. Using the large-scale simulation Horizon-AGN, we find that statistically the most massive BHs reside in environments with the largest galaxy number counts. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighbourhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement at larger scales, in agreement with recent observations. We find that AGN feedback in the surrounding galaxies is able to decrease their luminosity and stellar mass, and therefore to make them unobservable when using restrictive galaxy selection criteria. Radiation from the quasars can spread over large distances, which could affect the formation history of surrounding galaxies, but a careful analysis of these processes requires radiative transfer simulations. [ABSTRACT FROM AUTHOR]

In the local Universe, black holes of |$10^{5-6}\, {\rm M_\odot }$| are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of |$10^{5-6}\, {\rm M_\odot }$| and up to |$10^8\, {\rm M_\odot }$|⁠ , forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel & Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population. [ABSTRACT FROM AUTHOR]

Using a suite of three large cosmological hydrodynamical simulations, HORIZON-AGN, HORIZON-NOAGN (no AGN feedback) and HORIZON-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density slope differences between objects formed in these three simulations (hereafter, HAGN, HnoAGN and HDM) significantly evolve with time. More specifically, at high redshift (z ∼ 5), the mean central density profiles of HAGN and HnoAGN dark matter haloes tend to be much steeper than their HDM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z ∼ 1.5, these mean halo density profiles in HAGN have flattened, pummelled by powerful AGN activity ('quasar mode'): the integrated inner mass difference gaps with HnoAGN haloes have widened, and those with HDM haloes have narrowed. Fast forward 9.5 billion years, down to z = 0, and the trend reverses: HAGN halo mean density profiles drift back to a more cusped shape as AGN feedback efficiency dwindles ('radio mode'), and the gaps in integrated central mass difference with HnoAGN and HDM close and broaden, respectively. On the galaxy side, the story differs noticeably. Averaged stellar profile central densities and inner slopes are monotonically reduced by AGN activity as a function of cosmic time, resulting in better agreement with local observations. [ABSTRACT FROM AUTHOR]