No Arabic abstract
The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $sim10^6 M_odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.
We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.
At the highest redshifts, z>6, several tens of luminous quasars have been detected. The search for fainter AGN, in deep X-ray surveys, has proven less successful, with few candidates to date. An extrapolation of the relationship between black hole (BH) and bulge mass would predict that the sample of z>6 galaxies host relatively massive BHs (>1e6 Msun), if one assumes that total stellar mass is a good proxy for bulge mass. At least a few of these BHs should be luminous enough to be detectable in the 4Ms CDFS. The relation between BH and stellar mass defined by local moderate-luminosity AGN in low-mass galaxies, however, has a normalization that is lower by approximately an order of magnitude compared to the BH-bulge mass relation. We explore how this scaling changes the interpretation of AGN in the high-z Universe. Despite large uncertainties, driven by those in the stellar mass function, and in the extrapolation of local relations, one can explain the current non-detection of moderate-luminosity AGN in Lyman Break Galaxies if galaxies below 1e11 Msun are characterized by the low-normalization scaling, and, even more so, if their Eddington ratio is also typical of moderate-luminosity AGN rather than luminous quasars. AGN being missed by X-ray searches due to obscuration or instrinsic X-ray weakness also remain a possibility.
The population of massive black holes (MBHs) in dwarf galaxies is elusive, but fundamentally important to understand the coevolution of black holes with their hosts and the formation of the first collapsed objects in the Universe. While some progress was made in determining the X-ray detected fraction of MBHs in dwarfs, with typical values ranging from $0%$ to $6%$, their overall active fraction, ${cal A}$, is still largely unconstrained. Here, we develop a theoretical model to predict the multiwavelength active fraction of MBHs in dwarf galaxies starting from first principles and based on the physical properties of the host, namely, its stellar mass and angular momentum content. We find multiwavelength active fractions for MBHs, accreting at typically low rates, ranging from $5%$ to $22%$, and increasing with the stellar mass of the host as ${cal A} sim(log_{10}M_{star})^{4.5}$. If dwarfs are characterized by low-metallicity environments, the active fraction may reach $sim 30%$ for the most massive hosts. For galaxies with stellar mass in the range $10^7<M_{star} [M_{odot}]<10^{10}$, our predictions are in agreement with occupation fractions derived from simulations and semi-analytical models. Additionally, we provide a fitting formula to predict the probability of finding an active MBH in a dwarf galaxy from observationally derived data. This model will be instrumental to guide future observational efforts to find MBHs in dwarfs. The James Webb Space Telescope, in particular, will play a crucial role in detecting MBHs in dwarfs, possibly uncovering active fractions $sim 3$ times larger than current X-ray surveys.
Massive black holes (MBHs) are nowadays recognized as integral parts of galaxy evolution. Both the approximate proportionality between MBH and galaxy mass, and the expected importance of feedback from active MBHs in regulating star formation in their host galaxies point to a strong interplay between MBHs and galaxies. MBHs must form in the first galaxies and be fed by gas in these galaxies, with continuous or intermittent inflows that, at times, can be larger than the Eddington rate. Feedback from supernovae and from the MBHs themselves modulates the growth of the first MBHs. While current observational data only probe the most massive and luminous MBHs, the tip of the iceberg, we will soon be able to test theoretical models of MBH evolution on more normal MBHs: the MBHs that are indeed relevant in building the population that we observe in local galaxies, including our own Milky Way.
We analyse the demographics of black holes (BHs) in the large-volume cosmological hydrodynamical simulation Horizon-AGN. This simulation statistically models how much gas is accreted onto BHs, traces the energy deposited into their environment and, consequently, the back-reaction of the ambient medium on BH growth. The synthetic BHs reproduce a variety of observational constraints such as the redshift evolution of the BH mass density and the mass function. Strong self-regulation via AGN feedback, weak supernova feedback, and unresolved internal processes result in a tight BH-galaxy mass correlation. Starting at z~2, tidal stripping creates a small population of BHs over-massive with respect to the halo. The fraction of galaxies hosting a central BH or an AGN increases with stellar mass. The AGN fraction agrees better with multi-wavelength studies, than single-wavelength ones, unless obscuration is taken into account. The most massive halos present BH multiplicity, with additional BHs gained by ongoing or past mergers. In some cases, both a central and an off-centre AGN shine concurrently, producing a dual AGN. This dual AGN population dwindles with decreasing redshift, as found in observations. Specific accretion rate and Eddington ratio distributions are in good agreement with observational estimates. The BH population is dominated in turn by fast, slow, and very slow accretors, with transitions occurring at z=3 and z=2 respectively.