Massive black hole binaries (MBHBs) of $10^5 \, \rm M_\odot - 3 \times 10^7 \, \rm M_\odot $ merging in low redshift galaxies ($z\le4$) are sufficiently loud to be detected weeks before coalescence with the Laser Interferometer Space Antenna (LISA). This allows us to perform the parameter estimation $on$ $the$ $fly$, i.e. as a function of the time to coalescence during the inspiral phase, relevant for early warning of the planned LISA protected periods and for searches of electromagnetic signals. In this work, we study the evolution of the sky position, luminosity distance, chirp mass and mass ratio uncertainties as function of time left before merger. Overall, light systems with total intrinsic mass $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ are characterized by smaller uncertainties than heavy ones ($\rm M_{\rm tot} = 10^7 \, \rm M_\odot$) during the inspiral. Luminosity distance, chirp mass and mass ratio are well constrained at the end of the inspiral. Concerning sky position, at $z=1$, MBHBs with $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ can be localized with a median precision of $\simeq 10^2 \, \rm deg^2 (\simeq 1 \, \rm deg^2)$ at 1 month (1 hour) from merger, while the sky position of heavy MBHBs can be determined to $10 \, \rm deg^2$ only 1 hour before merger. However the uncertainty around the median values broadens with time, ranging in between 0.04 -- 20 $\rm deg^2$ (0.3 -- 3 $\times 10^3 \, \rm deg^2$) for light (heavy) systems at 1 hour before merger. At merger the sky localization improves down to $\simeq 10^{-1} \, \rm deg^2$ for all masses. For the benefit of the observer community, we provide the full set of data from our simulations and simple and ready-to-use analytical fits to describe the time evolution of uncertainties in the aforementioned parameters, valid for systems with total mass between $10^5$--$10^7 \, \rm M_\odot$ and redshift $0.3$--$3$.

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_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<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<10^10.5 Msun. 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 the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-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.

Accretion of stars on massive black holes (MBHs) can feed MBHs and generate tidal disruption events (TDEs). We introduce a new physically motivated model to self-consistently treat TDEs in cosmological simulations, and apply it to the assembly of a galaxy with final mass $3\times 10^{10}\mathrm{\, M}_\odot$ at $z=6$. This galaxy exhibits a TDE rate of $\sim 10^{-5}\mathrm{\, yr}^{-1}$, consistent with local observations but already in place when the Universe was one billion year old. A fraction of the disrupted stars participate in the growth of MBHs, dominating it until the MBH reaches mass $\sim 5 \times 10^5 \mathrm{\, M}_\odot$, but their contribution then becomes negligible compared to gas. TDEs could be a viable mechanism to grow light MBH seeds, but fewer TDEs are expected when the MBH becomes sufficiently massive to reach the luminosity of, and be detected as, an active galactic nucleus. Galaxy mergers bring multiple MBHs in the galaxy, resulting in an enhancement of the global TDE rate in the galaxy by $\sim 1$ order of magnitude during $100\mathrm{\, Myr}$ around mergers. This enhancement is not on the central MBH, but caused by the presence of MBHs in the infalling galaxies. This is the first self consistent study of TDEs in a cosmological environment and highlights that accretion of stars and TDEs are a natural process occurring in a Milky~Way-mass galaxy at early cosmic times.

We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on the effect of supernova feedback on the growth of black holes. By utilizing a semi-analytic galaxy formation model, not only do we account for the processes that drive the evolution of binaries at separations $\lesssim 1$ pc (gas-driven migration, stellar hardening and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by massive black hole pairs from kpc down to a few pc separation. We also include the effect of supernova feedback, which may eject a substantial amount of gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. In spite of the inclusion of these novel physical effects, we predict that the LISA detection rate should still be in excess of $\sim 2 \mbox{yr}^{-1}$, irrespective of the model for the seeds of the black hole population at high redshifts. However, scenarios in which black holes form from $\sim100 M_\odot$ seeds are more significantly impacted by the feedback from supernovae. We also present predictions for the mass ratio distribution of the merger population, and find that binaries typically have mass ratios between $\sim 0.1$ and $1$. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust, and show only a mild dependence on the model.

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. After the numerical MBH merger at the resolution limit, we post-process MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

Recent deep surveys have unravelled a population of faint active galactic nuclei (AGN) in the high redshift Universe, leading to various discussions on their nature and their role during the Epoch of Reionization. We use cosmological radiation-hydrodynamics simulations of a bright galaxy at z = 6 ($M_\star \gtrsim 10^{10} M_\odot$) hosting an actively growing super-massive black hole to study the properties of these objects. In particular, we study how the black hole and the galaxy co-evolve and what is the relative contribution of the AGN and of the stellar populations to the luminosity budget of the system. We find that the feedback from the AGN has no strong effect on the properties of the galaxy, and does not increase the total ionizing luminosity of the host. The average escape fraction of our galaxy is around $f_{\rm esc} \sim 5\%$. While our galaxy would be selected as an AGN in deep X-ray surveys, most of the UV luminosity is originating from stellar populations. This confirms that there is a transition in the galaxy population from star forming galaxies to quasar hosts, with bright Lyman-Break Galaxies (LBGs) with $M_{\rm UV}$ around -22 falling in the overlap region. Our results also suggest that faint AGN do not contribute significantly to reionizing the Universe.

We compute the tidal disruption event (TDE) rate around local massive black holes (MBHs) with masses as low as $2.5\times10^4\mathrm{ M}_\odot$, thus probing the dwarf regime for the first time. We select a sample of 37 galaxies for which we have the surface stellar density profile, a dynamical estimate of the mass of the MBH, and 6 of which, including our Milky Way, have a resolved nuclear star cluster (NSC). For the Milky Way, we find a total TDE rate of $\sim 10^{-4}\mathrm{ yr}^{-1}$ when taking the NSC in account, and $\sim 10^{-7} \mathrm{ yr}^{-1}$ otherwise. TDEs are mainly sourced from the NSC for light ($<3\times 10^{10}\mathrm{ M}_\odot$) galaxies, with a rate of few $10^{-5}\mathrm{ yr}^{-1}$, and an enhancement of up to 2 orders of magnitude compared to non-nucleated galaxies. We create a mock population of galaxies using different sets of scaling relations to explore trends with galaxy mass, taking into account the nucleated fraction of galaxies. Overall, we find a rate of few $10^{-5}\mathrm{ yr}^{-1}$ which drops when galaxies are more massive than $10^{11}\mathrm{ M}_\odot$ and contain MBHs swallowing stars whole and resulting in no observed TDE.

The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of $N$-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting a SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50% of Milky Way sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna -- LISA -- which will detect SMBH coalescence in the $10^4-10^7$$M_{\rm \odot}$ mass range.

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the lightcone of the Horizon-AGN simulation, for which we know the fate of each pair, to test three selection processes aiming at identifying true merging pairs. We find that the simplest one (selecting objects within two thresholds on projected distance $d$ and redshift difference $\Delta z$) gives similar results than the most complex one (based on a neural network analyzing $d$, $\Delta z$, redshift of the primary, masses/star formation rates/aspect ratio of both galaxies). Our best thresholds are $d_\mathrm{th}\sim100\mathrm{\, kpc}$ and $\Delta z_\mathrm{th} \sim 10^{-3}$, in agreement with recent results.

Supermassive black holes (BHs) of millions solar masses and above reside in the center of most local galaxies, but they also power active galactic nuclei and quasars, detected up to z=7. This observational evidence puts strong constraints on the BH growth and the mass of the first BH seeds. The scenario of "direct collapse" is very appealing as it leads to the formation of large mass BH seeds in the range 10^4-10^6 Msun, which eases explaining how quasars at z=6-7 are powered by BHs with masses >10^9 Msun. Direct collapse, however, appears to be rare, as the conditions required by the scenario are that gas is metal-free, the presence of a strong photo-dissociating Lyman-Werner flux, and large inflows of gas at the center of the halo, sustained for 10-100 Myr. We performed several cosmological hydrodynamical simulations that cover a large range of box sizes and resolutions, thus allowing us to understand the impact of several physical processes on the distribution of direct collapse BHs. We identify halos where direct collapse can happen, and derive the number density of BHs. We also investigate the discrepancies between hydrodynamical simulations, direct or post-processed, and semi-analytical studies. We find that for direct collapse to account for BHs in normal galaxies, the critical Lyman-Werner flux required for direct collapse must be much less than predicted by 3D simulations that include detailed chemical models. However, when supernova feedback is relatively weak, enough direct collapse BHs to explain z=6-7 quasars can be obtained for more realistic values of the critical Lyman Werner flux.

High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) 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 number counts of galaxies. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighborhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement of the galaxy number counts 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 un-observable 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.

Observations of high-redshift quasars provide information on the massive black holes (MBHs) powering them and the galaxies hosting them. Current observations of $z \gtrsim 6$ hosts, at sub-mm wavelengths, trace the properties of cold gas, and these are used to compare with the correlations between MBHs and galaxies characterising the $z=0$ population. The relations at $z=0$, however, rely on stellar-based tracers of the galaxy properties. We perform a very-high resolution cosmological zoom-in simulation of a $z=7$ quasar including state-of-the-art non-equilibrium chemistry, MBH formation, growth and feedback, to assess the evolution of the galaxy host and the central MBH, and compare the results with recent ALMA observations of high-redshift quasars. We measure both the stellar-based quantities used to establish the $z=0$ correlations, as well as the gas-based quantities available in $z \gtrsim 6$ observations, adopting the same assumptions and techniques used in observational studies. The high-redshift studies argued that MBHs at high redshift deviate from the local MBH-galaxy correlations. In our analysis of the single galaxy we evolve, we find that the high-redshift population sits on the same correlations as the local one, when using the same tracers used at $z=0$. When using the gas-based tracers, however, MBHs appear to be over-massive. The discrepancy between local and high-redshift MBHs seems caused by the different tracers employed, and necessary assumptions, and not by an intrinsic difference. Better calibration of the tracers, higher resolution data and availability of facilities that can probe the stellar population will be crucial to assess precisely and accurately high-redshift quasar hosts.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values.

In the local universe, black holes of $10^{5-6}$ 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}$ M$_{\odot}$ and up to $10^8$ 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.

While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$.

The dynamics of black hole seeds in high redshift galaxies is key to understand their ability to grow via accretion and to pair in close binaries during galactic mergers. To properly follow the dynamics of black holes we develop a physically motivated model to capture unresolved dynamical friction from stars, dark matter and gas. We first validate the model and then we use it to investigate the dynamics of seed black holes born at $z\sim9$ in dwarf proto-galaxies. We perform a suite of zoom cosmological simulations with spatial resolution as high as 10 pc and with a stellar and dark matter mass resolution of $2\times10^3 \, \textrm{M}_\odot$ and $2\times10^5 \, \textrm{M}_\odot$ respectively. We first explore the dynamics of a seed black hole in the galaxy where it is born and show that it is highly erratic if the seed mass is less than $10^5\, \textrm{M}_\odot$. The dynamics is dominated by the stellar component, whose distribution is irregular and patchy, thus inducing stochasticity in the orbits: the black hole may be anywhere in the proto-galaxy. When this dwarf merges into a larger galaxy, it is paramount to simulate the process with very high spatial and mass resolution in order to correctly account for the stripping of the stellar envelope of the satellite black hole. The outcome of the encounter could be either a tight binary or, at least temporary, a wandering black hole, leading to multiple black holes in a galaxy, each inherited from a different merger.

Distinguishing the photon output of an accreting supermassive black hole binary system from that of a single supermassive black hole accreting at the same rate is intrinsically difficult because the majority of the light emerges from near the ISCOs of the black holes. However, there are two possible signals that can distinctively mark binaries, both arising from the gap formed in circumbinary accretion flows inside approximately twice the binary separation. One of these is a "notch" cut into the thermal spectra of these systems in the IR/optical/UV, the other a periodically-varying excess hard X-ray luminosity whose period is of order the binary orbital period. Using data from detailed galaxy evolution simulations, we estimate the distribution function in mass, mass ratio, and accretion rate for accreting supermassive black hole binaries as a function of redshift and then transform this distribution function into predicted source counts for these two potential signals. At flux levels >~10^{-13}~erg~cm^{-2}~s^{-1}, there may be ~O(10^2) such systems in the sky, mostly in the redshift range ~0.5 < z < ~1. Roughly 10% should have periods short enough (<~5~yr) to detect the X-ray modulation; this is also the period range accessible to PTA observations.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

Super-Eddington accretion onto massive black hole seeds may be commonplace in the early Universe, where the conditions exist for rapid accretion. Direct collapse black holes are often invoked as a possible solution to the observation of super massive black holes (SMBHs) in the pre-reionisation Universe. We investigate here how feedback, mainly in the form of bipolar jets, from super-Eddington accreting seed black holes will affect their subsequent growth. We find that, nearly independent of the mass loading of the bipolar jets, the violent outflows generated by the jets evacuate a region of approximately 0.1 pc surrounding the black hole seed. However, the jet outflows are unable to break free of the halo and their impact is limited to the immediate vicinity of the black hole. The outflows suppress any accretion for approximately a dynamical time. The gas then cools, recombines and falls back to the centre where high accretion rates are again observed. The overall effect is to create an effective accretion rate with values of between 0.1 and 0.5 times the Eddington rate. If this episodic accretion rate is maintained for order 500 million years then the black hole will increase in mass by a factor of between 3 and 300 but far short of the factor of $10^4$ required for the seeds to become the SMBHs observed at $z>6$. Therefore, direct collapse black holes born into atomic cooling haloes and which experience strong negative mechanical feedback will require external influences (e.g. rapid major mergers with other haloes) to promote efficient accretion and reach SMBH masses within a few hundred million years.

In this work we include black hole (BH) seeding, growth and feedback into our semi-analytic galaxy formation model, Delphi. Our model now fully tracks the, accretion- and merger-driven, hierarchical assembly of the dark matter halo, baryonic and BH masses of high-redshift ($z>5$) galaxies. We use a minimal set of mass- and $z$-independent free parameters associated with star formation and BH growth (and feedback) and include suppressed BH growth in low-mass galaxies to explore a number of physical scenarios including: (i) two types of BH seeds (stellar and those from Direct Collapse BH; DCBH); (ii) the impact of reionization feedback; and (iii) the impact of instantaneous versus delayed galaxy mergers on the baryonic growth. While both reionization feedback and delayed galaxy mergers have no sensible impact on the evolving ultra-violet luminosity function, the latter limits the maximum BH masses achieved at these high-$z$. We then use this model, baselined against all available high-$z$ galaxy and BH data-sets, to predict the LISA detectability of merger events at $z > 5$. As expected, the merger rate is dominated by stellar BH mergers for all scenarios and our model predicts an expected upper limit of about 20 mergers in the case of instantaneous merging and no reionization feedback over the 4-year mission duration. Including the impact of delayed mergers and reionization feedback reduces this to about 12 events over the same observational time-scale.

Supermassive Black Holes (BHs) residing in brightest cluster galaxies (BCGs) are overly massive when considering the local relationships between the BH mass and stellar bulge mass or velocity dispersion. Due to the location of these BHs within the cluster, large-scale cluster processes may aid the growth of BHs in BCGs. In this work, we study a sample of 71 galaxy clusters to explore the relationship between the BH mass, stellar bulge mass of the BCG, and the total gravitating mass of the host clusters. Due to difficulties in obtaining dynamically measured BH masses in distant galaxies, we use the Fundamental Plane relationship of BHs to infer their masses. We utilize X-ray observations taken by $Chandra$ to measure the temperature of the intra-cluster medium (ICM), which is a proxy for the total mass of the cluster. We analyze the $\rm M_{BH}-kT$ and $\rm M_{BH}-M_{Bulge}$ relationships and establish the best-fitting power laws:$\log_{10}(M_{\rm BH} /10^9 M_{\odot})=-0.35+2.08 \log_{10}(kT / 1 \rm keV)$ and $\log_{10}(\rm M_{BH}/10^9M_{\odot})= -1.09+ 1.92 \log_{10}(M_{\rm bulge}/10^{11}M_{\odot})$. Both relations are comparable with that established earlier for a sample of brightest group/cluster galaxies with dynamically measured BH masses. Although both the $\rm M_{BH}-kT$ and the $\rm M_{BH}-M_{Bulge}$ relationships exhibit large intrinsic scatter, based on Monte Carlo simulations we conclude that dominant fraction of the scatter originates from the Fundamental Plane relationship. We split the sample into cool core and non-cool core resembling clusters, but do not find statistically significant differences in the $\rm M_{BH}-kT$ relation. We speculate that the overly massive BHs in BCGs may be due to frequent mergers and cool gas inflows onto the cluster center.

The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows that trigger nuclear starbursts that increase the central stellar density. Motivated by these numerical results, and by the observed over-representation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs) in the merging galaxies and in the final remnant. We post-process the simulation to obtain the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate, until the two BHs bind and are surrounded by the same stellar environment. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble time. We find that black hole binaries which are ejected and later merge ex-situ form in clusters of a characteristic mass $M_{\rm cl} \sim 10^{5.3}M_{\odot}$, whereas binaries which merge in-situ form in more massive clusters, $M_{\rm cl} \sim 10^{5.7}M_{\odot}$. The clusters which dominate the production of black hole binaries are similar in age and metallicity to the entire population. Finally, we estimate an approximate cosmic black hole merger rate of dynamically assembled binaries using the mean black hole mass for each cluster given its metallicity. We find an intrinsic rate of $\sim 6\,\mathrm{Gpc\,yr^{-1}}$ at $z=0$, a weakly increasing merger rate out to $z=1.5$, and then a decrease out to $z=4$. Our results can be used to provide a cosmological context and choose initial conditions in numerical studies of black hole binaries assembled in star clusters.

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Dense gas in the centre of galaxies feeds massive black holes, but can also become a source of obscuration and limit our ability to find faint Active Galactic Nuclei (AGN). We use a high resolution cosmological radiative hydrodynamics simulation to connect the properties of the gas in the central region (a few tens of parsecs) of a high redshift galaxy to the growth of a massive black hole during the first billion years of the Universe. We find that the feedback from the AGN efficiently controls the growth of the black hole and limits the duration of the high accretion episodes by emptying the gas reservoir. As the galaxy grows in mass, the production of metals results in the presence of dust-enriched gas in the galaxy centre that can obscure highly accreting black holes enough to strongly reduce their UV/optical visibility. We also find that the gas outside the very centre of the galaxy can contribute to the total column density and obscuration at a level at least comparable to the gas in the nuclear region. We suggest that this explains the different duty cycles required to explain the masses of high redshift quasars and the observed UV/optical luminosity functions: in our case, the AGN would be observed with an X-ray luminosity above $L_X = 10^{42}$ erg/s around 30% of the time, but with UV magnitude brighter than $M_{1450}$ = -23 only 4% of the time.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

We report the results of high-resolution subpixel imaging of the hard continuum and Fe K{\alpha} line of the Compton Thick (CT) Active Galactic Nucleus (AGN) ESO 428-G014, observed with Chandra ACIS. While the 3-4 keV emission is dominated by an extended component, a single nuclear point source is prominent in the 4-6 keV range. Instead, two peaks of similar intensity, separated by ~36 pc in projection on the plane of the sky are detected in the Fe K{\alpha} emission. The SE knot could be marginally associated with the heavily obscured hard continuum source. We discuss four possible interpretations of the nuclear morphology. (1) Given the bolometric luminosity and likely black hole (BH) mass of ESO 428-G014, we may be imaging two clumps of the CT obscuring torus in the Fe K{\alpha} line. (2) The Fe K{\alpha} knots may be connected with the fluorescent emission from the dusty bicone, or (3) with the light echo of a nuclear outburst. (4) We also explore the less likely possibility that we may be detecting the rare signature of merging nuclei. Considering the large-scale kpc-size extent of the hard continuum and Fe K{\alpha} emission (Papers I and II), we conclude that the AGN in ESO 428-G014 has been active for at least 104 yrs. Comparison with the models of Czerny et al (2009) suggests high accretion rates during this activity.

[abridged] Aims: We test the effects of re-orienting jets from an active galactic nucleus (AGN) on the intracluster medium in a galaxy cluster environment with short central cooling time. We investigate appearance and properties of the resulting cavities, and the efficiency of jets in providing near-isotropic heating to the cooling cluster core. Methods: We use numerical simulations to explore four models of jets over several active/inactive cycles. We keep the jet power and duration fixed, varying only the jet angle prescription. We track the total energy of the intracluster medium (ICM) in the cluster core over time, and the fraction of the jet energy transferred to the ICM, paying attention to where the energy is deposited. We also compare synthetic X-ray images of the simulated cluster to actual observations. Results: Jets whose re-orientation is minimal ($\lesssim 20^{\circ}$) typically produce conical structures of interconnected cavities, with the opening angle of the cones being $\sim 15-20^{\circ}$, extending to $\sim 300$ kpc from the cluster centre. Such jets transfer about $60\%$ of their energy to the ICM, yet they are not very efficient at heating the cluster core, as the jet energy is deposited further out. Jets that re-orient by $\gtrsim 20^{\circ}$ generally produce multiple pairs of detached cavities. Although smaller, these cavities are inflated within the central 50~kpc and are more isotropically distributed, resulting in more effective heating of the core. Such jets, over few hundreds Myr, can deposit up to $80\%$ of their energy where it is required. Consequently, these models come the closest to an heating/cooling balance and to mitigating runaway cooling of the core, even though all models have identical power/duration profiles. Additionally, the corresponding synthetic X-ray images exhibit structures closely resembling those seen in real cool-core clusters.

We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $\Lambda$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.

The gravitational interaction at parsec to sub-parsec scales between a circumbinary gas disc and a super massive black hole binary (SMBHB) is a promising mechanism to drive the migration of SMBHBs toward coalescence. The typical dynamical evolution can be separated in two main regimes: I) Slow migration ($T_{\rm mig}$ $\sim$ $10^{3-4}\times T_{\rm orb}$), where viscous torques are not efficient enough to redistribute the extracted angular momentum from the binary, leading to the formation of a low density cavity around the binary. II) Fast migration ($T_{\rm mig}$ $\sim$ $10^{1-2} \times T_{\rm orb}$), in which the redistribution of angular momentum is efficient and no low density cavity is formed in the circumbinary disc. Using N-Body/SPH simulations we study the effect of AGN feedback in this phase of a SMBHB evolution. We implement an AGN feedback model in the SPH code Gadget-3 that includes momentum feedback from winds, X-ray heating/radial-momentum and Eddington force. Using this implementation we run a set of simulations of SMBHB+disc in the two main shrinking regimes. From these simulations we conclude that the effect of the AGN mechanical feedback is negligible for SMBHBs in the slowly shrinking regime. However, in the fast shrinking regime the AGN wind excavate a "feedback cavity" leaving the SMBHB almost naked, thus stalling the orbital decay of the binary.

Understanding the growth of high redshift massive black holes (MBHs) is a problem of great astrophysical interest. The most luminous quasars at $z>6$ are frequently observed but they represent only the tip of the iceberg as the majority of the low luminosity AGN population remains undetected. In the present study, we perform a radiation hydrodynamics cosmological simulation to study the growth of "normal" black holes in the high redshift universe. In our simulation we model the formation of Pop III and Pop II stars along with their chemical, mechanical and radiative feedback. We consider both UV and X-ray emission from an accreting BH to simulate its radiative feedback. The selected halo has a mass of $\rm 3 \times 10^{10}~M_{\odot}$ at $z=7.5$ and we turn on radiative feedback from a MBH seed of $\rm 10^5~M_{\odot}$ along with in-situ star formation at $z=12$ when the halo mass reaches well above the atomic cooling limit. We find that the MBH accretes only about 2200 $\rm M_{\odot}$ during 320 Myr and the average mass accretion onto the MBH is a few times $\rm 10^{-6}~M_{\odot}/yr$. Our results suggest that the stunted growth of MBH is a consequence of supernovae in tandem with MBH feedback which drive large outflows and evacuate the gas from MBH vicinity. This may explain why a population of low luminosity AGN has not been detected so-far at $z>6$; the large contrast between the star formation rate and the MBH accretion rate may make then hard to detect even in upcoming deep surveys.

Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This study is an extension of Bellovary et. al. 2013, which analyzed the accretion of high mass, high redshift galaxies and their central black holes, and found that the gas content of the central black hole reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by merger gas. Moreover, we find that nearly 30$\%$ of the accreted mass budget of the SMBH enters the galaxy through the two major mergers in its history, which may account for the increase of merger-gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that merger gas enters the galaxy with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories.

We present a sample of 40 AGN in dwarf galaxies at redshifts $z \lesssim$ 2.4. The galaxies are drawn from the \textit{Chandra} COSMOS-Legacy survey as having stellar masses $10^{7}\leq M_{*}\leq3 \times 10^{9}$ M$_{\odot}$. Most of the dwarf galaxies are star-forming. After removing the contribution from star formation to the X-ray emission, the AGN luminosities of the 40 dwarf galaxies are in the range $L_\mathrm{0.5-10 keV} \sim10^{39} - 10^{44}$ erg s$^{-1}$. With 12 sources at $z > 0.5$, our sample constitutes the highest-redshift discovery of AGN in dwarf galaxies. The record-holder is cid\_1192, at $z = 2.39$ and with $L_\mathrm{0.5-10 keV} \sim 10^{44}$ erg s$^{-1}$. One of the dwarf galaxies has $M_\mathrm{*} = 6.6 \times 10^{7}$ M$_{\odot}$ and is the least massive galaxy found so far to host an AGN. All the AGN are of type 2 and consistent with hosting intermediate-mass black holes (BHs) with masses $\sim 10^{4} - 10^{5}$ M$_{\odot}$ and typical Eddington ratios $> 1\%$. We also study the evolution, corrected for completeness, of AGN fraction with stellar mass, X-ray luminosity, and redshift in dwarf galaxies out to $z$ = 0.7. We find that the AGN fraction for $10^{9}< M_{*}\leq3 \times 10^{9}$ M$_{\odot}$ and $L_\mathrm{X} \sim 10^{41}-10^{42}$ erg s$^{-1}$ is $\sim$0.4\% for $z \leq$ 0.3 and that it decreases with X-ray luminosity and decreasing stellar mass. Unlike massive galaxies, the AGN fraction seems to decrease with redshift, suggesting that AGN in dwarf galaxies evolve differently than those in high-mass galaxies. Mindful of potential caveats, the results seem to favor a direct collapse formation mechanism for the seed BHs in the early Universe.

We present a self-consistent prediction from a large-scale cosmological simulation for the population of `wandering' supermassive black holes (SMBHs) of mass greater than $10^6$ M$_{\odot}$ on long-lived, kpc-scale orbits within Milky Way (MW)-mass galaxies. We extract a sample of MW-mass halos from the Romulus25 cosmological simulation (Tremmel et al. 2017), which is uniquely able to capture the orbital evolution of SMBHs during and following galaxy mergers. We predict that such halos, regardless of recent merger history or morphology, host an average of $5.1 \pm 3.3$ SMBHs, including their central black hole, within 10 kpc from the galactic center and an average of $12.2 \pm 8.4$ SMBHs total within their virial radius, not counting those in satellite halos. Wandering SMBHs exist within their host galaxies for several Gyrs, often accreted by their host halo in the early Universe. We find, with $>4\sigma$ significance, that wandering SMBHs are preferentially found outside of galactic disks.

Galactic outflows, driven by active galactic nuclei (AGN), play a crucial role in galaxy formation and in the self-regulated growth of supermassive black holes (BHs). AGN feedback couples to and affects gas, rather than stars, and in many, if not most, gas-rich galaxies cold gas is rotationally supported and settles in a disc. We present a 2D analytical model for AGN-driven outflows in a gaseous disc and demonstrate the main improvements, compared to existing 1D solutions. We find significant differences for the outflow dynamics and wind efficiency. The outflow is energy-driven due to inefficient cooling up to a certain AGN luminosity ($\sim 10^{43}$erg/s in our fiducial model), above which the outflow remains momentum-driven in the disc up to galactic scales. We reproduce results of 3D simulations that gas is preferentially ejected perpendicular to the disc and find that the fraction of ejected interstellar medium is lower than in 1D models. The recovery time of gas in the disc, defined as the freefall time from the radius to which the AGN pushes the ISM at most, is remarkably short, of the order 1Myr. This indicates that AGN-driven winds cannot suppress BH growth for long. Without the inclusion of supernova feedback, we find a scaling of the black hole mass with the halo velocity dispersion of $M_\mathrm{BH} \propto \sigma ^{4.8}$.

Understanding the processes that drive the formation of black holes (BHs) is a key topic in observational cosmology. While the observed $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ correlation in bulge-dominated galaxies is thought to be produced by major mergers, the existence of a $M_{\mathrm{BH}}$--$M_{\star}$ relation, across all galaxy morphological types, suggests that BHs may be largely built by secular processes. Recent evidence that bulge-less galaxies, which are unlikely to have had significant mergers, are offset from the $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but lie on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, has strengthened this hypothesis. Nevertheless, the small size and heterogeneity of current datasets, coupled with the difficulty in measuring precise BH masses, makes it challenging to address this issue using empirical studies alone. Here, we use Horizon-AGN, a cosmological hydrodynamical simulation to probe the role of mergers in BH growth over cosmic time. We show that (1) as suggested by observations, simulated bulge-less galaxies lie offset from the main $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, (2) the positions of galaxies on the $M_{\mathrm{BH}}$--$M_{\star}$ relation are not affected by their merger histories and (3) only $\sim$35 per cent of the BH mass in today's massive galaxies is directly attributable to merging -- the majority ($\sim$65 per cent) of BH growth, therefore, takes place gradually, via secular processes, over cosmic time.

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

Observing the light emitted by the first accreting black holes (BHs) would dramatically improve our understanding of the formation of quasars at z > 6, possibly unveiling the nature of their supermassive black hole (SMBH) seeds. In previous works we explored the relative role of the two main competing BH seed formation channels, Population III remnants (low-mass seeds) and direct collapse BHs (high-mass seeds), investigating the properties of their host galaxies in a cosmological context. Building on this analysis, we predict here the spectral energy distribution and observational features of low- and high-mass BH seeds selected among the progenitors of a z~6 SMBH. We derive the processed emission from both accreting BHs and stars by using the photo-ionization code Cloudy, accounting for the evolution of metallicity and dust-to-gas mass ratio in the interstellar medium of the host galaxies, as predicted by the cosmological data- constrained model GAMETE/QSOdust. We show how future missions like JWST and ATHENA will be able to detect the light coming from SMBH progenitors already at z~16. We build upon previous complementary studies and propose a method based on the combined analysis of near infrared (NIR) colors, IR excess (IRX) and UV continuum slopes (i.e. color-color and IRX-Beta diagrams) to distinguish growing seed BH host galaxies from starburst-dominated systems in JWST surveys. Sources selected through this criterion would be the best target for follow-up X-ray observations.

Supermassive black holes (BHs) residing in the brightest cluster galaxies are over-massive relative to the stellar bulge mass or central stellar velocity dispersion of their host galaxies. As BHs residing at the bottom of the galaxy cluster's potential well may undergo physical processes that are driven by the large-scale characteristics of the galaxy clusters, it is possible that the growth of these BHs is (indirectly) governed by the properties of their host clusters. In this work, we explore the connection between the mass of BHs residing in the brightest group/cluster galaxies (BGGs/BCGs) and the virial temperature, and hence total gravitating mass, of galaxy groups/clusters. To this end, we investigate a sample of 17 BGGs/BCGs with dynamical BH mass measurements and utilize XMM-Newton X-ray observations to measure the virial temperatures and infer the $M_{\rm 500}$ mass of the galaxy groups/clusters. We find that the $M_{\rm BH} - kT$ relation is significantly tighter and exhibits smaller scatter than the $M_{\rm BH} - M_{\rm bulge}$ relations. The best-fitting power-law relations are $ \log_{10} (M_{\rm BH}/10^{9} \ \rm{M_{\odot}}) = 0.20 + 1.74 \log_{10} (kT/1 \ \rm{keV}) $ and $ \log_{10} (M_{\rm BH}/10^{9} \ \rm{M_{\odot}}) = -0.80 + 1.72 \log_{10} (M_{\rm bulge}/10^{11} \ M_{\odot})$. Thus, the BH mass of BGGs/BCGs may be set by physical processes that are governed by the properties of the host galaxy group/cluster. These results are confronted with the Horizon-AGN simulation, which reproduces the observed relations well, albeit the simulated relations exhibit notably smaller scatter.

While low mass, star forming galaxies are often considered as the primary driver of reionization, their actual contribution to the cosmic ultraviolet background is still uncertain, mostly because the escape fraction of ionizing photons is only poorly constrained. Theoretical studies have shown that efficient supernova feedback is a necessary condition to create paths through which ionizing radiation can escape into the intergalactic medium. We investigate the possibility that accreting supermassive black holes in early dwarf galaxies may provide additional feedback and enhance the leakage of ionizing radiation. We use a series of high resolution cosmological radiation hydrodynamics simulations where we isolate the different sources of feedback. We find that supernova feedback prevents the growth of the black hole, thus quenching its associated feedback. Even in cases where the black hole can grow, the structure of the interstellar medium is strongly dominated by supernova feedback. We conclude that, in the dwarf galaxy regime, supermassive black holes do not appear to play a significant role in enhancing the escape fraction and in contributing to the early UV background.

Feedback by Active Galactic Nuclei is often divided into quasar and radio mode, powered by radiation or radio jets, respectively. Both are fundamental in galaxy evolution, especially in late-type galaxies, as shown by cosmological simulations and observations of jet-ISM interactions in these systems. We compare AGN feedback by radiation and by collimated jets through a suite of simulations, in which a central AGN interacts with a clumpy, fractal galactic disc. We test AGN of $10^{43}$ and $10^{46}$ erg/s, considering jets perpendicular or parallel to the disc. Mechanical jets drive the more powerful outflows, exhibiting stronger mass and momentum coupling with the dense gas, while radiation heats and rarifies the gas more. Radiation and perpendicular jets evolve to be quite similar in outflow properties and effect on the cold ISM, while inclined jets interact more efficiently with all the disc gas, removing the densest $20\%$ in $20$ Myr, and thereby reducing the amount of cold gas available for star formation. All simulations show small-scale inflows of $0.01-0.1$ M$_\odot$/yr, which can easily reach down to the Bondi radius of the central supermassive black hole (especially for radiation and perpendicular jets), implying that AGN modulate their own duty cycle in a feedback/feeding cycle.

Super-massive black holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth, since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. Using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.

The observed massive end of the galaxy stellar mass function is steeper than its predicted dark matter halo counterpart in the standard $\Lambda $CDM paradigm. In this paper, we investigate the impact of active galactic nuclei (AGN) feedback on star formation in massive galaxies. We isolate the impact of AGNs by comparing two simulations from the HORIZON suite, which are identical except that one also includes super massive black holes (SMBH), and related feedback models. This allows us to cross-identify individual galaxies between simulations and quantify the effect of AGN feedback on their properties, including stellar mass and gas outflows. We find that massive galaxies ($ \rm M_{*} \geq 10^{11} M_\odot $) are quenched by AGN feedback to the extent that their stellar masses decrease by up to 80% at $z=0$. SMBHs affect their host halo through a combination of outflows that reduce their baryonic mass, particularly for galaxies in the mass range $ \rm 10^9 M_\odot \leq M_{*} \leq 10^{11} M_\odot $, and a disruption of central gas inflows, which limits in-situ star formation. As a result, net gas inflows onto massive galaxies, $ \rm M_{*} \geq 10^{11} M_\odot $, drop by up to 70%. We measure a redshift evolution in the stellar mass ratio of twin galaxies with and without AGN feedback, with galaxies of a given stellar mass showing stronger signs of quenching earlier on. This evolution is driven by a progressive flattening of the $\rm M_{\rm SMBH}-M_* $ relation with redshift, particularly for galaxies with $\rm M_{*} \leq 10^{10} M_\odot $. $\rm M_{\rm SMBH}/M_*$ ratios decrease over time, as falling average gas densities in galaxies curb SMBH growth.

The first billion years of the Universe is a pivotal time: stars, black holes (BHs) and galaxies form and assemble, sowing the seeds of galaxies as we know them today. Detecting, identifying and understand- ing the first galaxies and BHs is one of the current observational and theoretical challenges in galaxy formation. In this paper we present a population synthesis model aimed at galaxies, BHs and Active Galactic Nuclei (AGNs) at high redshift. The model builds a population based on empirical relations. Galaxies are characterized by a spectral energy distribution determined by age and metallicity, and AGNs by a spectral energy distribution determined by BH mass and accretion rate. We validate the model against observational constraints, and then predict properties of galaxies and AGN in other wavelength and/or luminosity ranges, estimating the contamination of stellar populations (normal stars and high-mass X-ray binaries) for AGN searches from the infrared to X-rays, and vice-versa for galaxy searches. For high-redshift galaxies, with stellar ages < 1 Gyr, we find that disentangling stellar and AGN emission is challenging at restframe UV/optical wavelengths, while high-mass X-ray binaries become more important sources of confusion in X-rays. We propose a color-color selection in JWST bands to separate AGN vs star-dominated galaxies in photometric observations. We also esti- mate the AGN contribution, with respect to massive, hot, metal-poor stars, at driving high ionization lines, such as C IV and He II. Finally, we test the influence of the minimum BH mass and occupa- tion fraction of BHs in low mass galaxies on the restframe UV/near-IR and X-ray AGN luminosity function.

In this study, we present a suite of high-resolution numerical simulations of an isolated galaxy to test a sub-grid framework to consistently follow the formation and dissociation of H$_2$ with non-equilibrium chemistry. The latter is solved via the package KROME, coupled to the mesh-less hydrodynamic code GIZMO. We include the effect of star formation (SF), modelled with a physically motivated prescription independent of H$_2$, supernova feedback and mass losses from low-mass stars, extragalactic and local stellar radiation, and dust and H$_2$ shielding, to investigate the emergence of the observed correlation between H$_2$ and SF rate surface densities. We present two different sub-grid models and compare them with on-the-fly radiative transfer (RT) calculations, to assess the main differences and limits of the different approaches. We also discuss a sub-grid clumping factor model to enhance the H$_2$ formation, consistent with our SF prescription, which is crucial, at the achieved resolution, to reproduce the correlation with H$_2$. We find that both sub-grid models perform very well relative to the RT simulation, giving comparable results, with moderate differences, but at much lower computational cost. We also find that, while the Kennicutt-Schmidt relation for the total gas is not strongly affected by the different ingredients included in the simulations, the H$_2$-based counterpart is much more sensitive, because of the crucial role played by the dissociating radiative flux and the gas shielding.

We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts.

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, H_AGN, H_noAGN and H_DM) significantly evolve with time. More specifically, at high redshift (z~5), the mean central density profiles of H_AGN and H_noAGN dark matter haloes tend to be much steeper than their H_DM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z~1.5, these mean halo density profiles in H_AGN have flattened, pummelled by powerful AGN activity ("quasar mode"): the integrated inner mass difference gaps with H_noAGN haloes have widened, and those with H_DM haloes have narrowed. Fast forward 9.5 billion years, down to z=0, and the trend reverses: H_AGN 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 H_noAGN and H_DM 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. As both dark matter and stellar inner density profiles respond quite sensitively to the presence of a central AGN, there is hope that future observational determinations of these quantities can be used constrain AGN feedback models.

Massive black holes (BHs) inhabit local galaxies, including the Milky Way and some dwarf galaxies. BH formation, occurring at early cosmic times, must account for the properties of BHs in today's galaxies, notably why some galaxies host a BH, and others do not. We investigate the formation, distribution and growth of BH `seeds' by using the adaptive mesh refinement code Ramses. We develop an implementation of BH formation in dense, low-metallicity environments, as advocated by models invoking the collapse of the first generation of stars, or of dense nuclear star clusters. The seed masses are computed one-by-one on-the-fly, based on the star formation rate and the stellar initial mass function. This self-consistent method to seed BHs allows us to study the distribution of BHs in a cosmological context and their evolution over cosmic time. We find that all high-mass galaxies tend to a host a BH, whereas low-mass counterparts have a lower probability of hosting a BH. After the end of the epoch of BH formation, this probability is modulated by the growth of the galaxy. The simulated BHs connect to low-redshift observational samples, and span a similar range in accretion properties as Lyman-Break Analogs. The growth of BHs in low-mass galaxies is stunted by strong supernova feedback. The properties of BHs in dwarf galaxies thus remain a testbed for BH formation. Simulations with strong supernova feedback, which is able to quench BH accretion in shallow potential wells, produce galaxies and BHs in better agreement with observational constraints.

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

Super-Eddington accretion has been suggested as a possible formation pathway of $10^9 \, M_\odot$ supermassive black holes (SMBHs) 800 Myr after the Big Bang. However, stellar feedback from BH seed progenitors and winds from BH accretion disks may decrease BH accretion rates. In this work, we study the impact of these physical processes on the formation of $z \sim 6$ quasar, including new physical prescriptions in the cosmological, data-constrained semi-analytic model GAMETE/QSOdust. We find that the feedback produced by the first stellar progenitors on the surrounding does not play a relevant role in preventing SMBHs formation. In order to grow the $z \gtrsim 6$ SMBHs, the accreted gas must efficiently lose angular momentum. Moreover disk winds, easily originated in super-Eddington accretion regime, can strongly reduce duty cycles. This produces a decrease in the active fraction among the progenitors of $z\sim6$ bright quasars, reducing the probability to observe them.

Numerical simulations have become a necessary tool to describe the complex interactions among the different processes involved in galaxy formation and evolution, unfeasible via an analytic approach. The last decade has seen a great effort by the scientific community in improving the sub-grid physics modelling and the numerical techniques used to make numerical simulations more predictive. Although the recently publicly available code GIZMO has proven to be successful in reproducing galaxy properties when coupled with the model of the MUFASA simulations and the more sophisticated prescriptions of the Feedback In Realistic Environments (FIRE) set-up, it has not been tested yet using delayed cooling supernova feedback, which still represent a reasonable approach for large cosmological simulations, for which detailed sub-grid models are prohibitive. In order to limit the computational cost and to be able to resolve the disc structure in the galaxies we perform a suite of zoom-in cosmological simulations with rather low resolution centred around a sub-L* galaxy with a halo mass of $3\times 10^{11}\,\rm M_\odot$ at $z=0$, to investigate the ability of this simple model, coupled with the new hydrodynamic method of GIZMO, to reproduce observed galaxy scaling relations (stellar to halo mass, stellar and baryonic Tully-Fisher, stellar mass-metallicity and mass-size). We find that the results are in good agreement with the main scaling relations, except for the total stellar mass, larger than that predicted by the abundance matching technique, and the effective sizes for the most massive galaxies in the sample, which are too small.

We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from $20\,\rm pc$ down to $1\,\rm pc$, in order to study the effects of the resolution on our results, remove numerical effects, and assess that resolving the influence radius of the orbiting black hole is a minimum condition to fully capture the formation of the binary. Our simulations include the relevant physical processes, namely star formation, supernova feedback, accretion onto the black holes and the ensuing feedback. We find that, in these mergers, dynamical friction from the smooth stellar component of the nucleus is the main process that drives black holes from kpc to pc scales. Gas does not play a crucial role and even clumps do not induce scattering or perturb the orbits. We compare the time needed for the formation of the binary to analytical predictions and suggest how to apply such analytical formalism to obtain estimates of binary formation times in lower resolution simulations.

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.

Black holes (BHs) are believed to be a key ingredient of galaxy formation. However, the galaxy-BH interplay is challenging to study due to the large dynamical range and complex physics involved. As a consequence, hydrodynamical cosmological simulations normally adopt sub-grid models to track the unresolved physical processes, in particular BH accretion; usually the spatial scale where the BH dominates the hydrodynamical processes (the Bondi radius) is unresolved, and an approximate Bondi-Hoyle accretion rate is used to estimate the growth of the BH. By comparing hydrodynamical simulations at different resolutions (300, 30, 3 pc) using a Bondi-Hoyle approximation to sub-parsec runs with non-parameterized accretion, our aim is to probe how well an approximated Bondi accretion is able to capture the BH accretion physics and the subsequent feedback on the galaxy. We analyse an isolated galaxy simulation that includes cooling, star formation, Type Ia and Type II supernovae, BH accretion and AGN feedback (radiation pressure, Compton heating/cooling) where mass, momentum, and energy are deposited in the interstellar medium through conical winds. We find that on average the approximated Bondi formalism can lead to both over- and under-estimations of the BH growth, depending on resolution and on how the variables entering into the Bondi-Hoyle formalism are calculated.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

The coeval AGN and galaxy evolution and the observed local relations between SMBHs and galaxy properties suggest some connection or feedback between SMBH growth and galaxy build-up. We looked for correlations between properties of X-ray detected AGN and their FIR detected host galaxies, to find quantitative evidences for this connection, highly debated in the latest years. We exploit the rich multi-wavelength data set available in the COSMOS field for a large sample (692 sources) of AGN and their hosts, in the redshift range $0.1<z<4$. We use X-ray data to select AGN and determine their properties (intrinsic luminosity and nuclear obscuration), and broad-band SED fitting to derive host galaxy properties (stellar mass $M_*$ and star formation rate SFR). We find that the AGN 2-10 keV luminosity ($L_{\rm X}$) and the host $8-1000~\mu m$ star formation luminosity ($L_{\rm IR}^{\rm SF}$) are significantly correlated. However, the average host $L_{\rm IR}^{\rm SF}$ has a flat distribution in bins of AGN $L_{\rm X}$, while the average AGN $L_{\rm X}$ increases in bins of host $L_{\rm IR}^{\rm SF}$, with logarithmic slope of $\sim0.7$, in the redshifts range $0.4<z<1.2$. We also discuss the comparison between the distribution of these two quantities and the predictions from hydro-dynamical simulations. Finally we find that the average column density ($N_H$) shows a positive correlation with the host $M_*$, at all redshifts, but not with the SFR (or $L_{\rm IR}^{\rm SF}$). This translates into a negative correlation with specific SFR. Our results are in agreement with the idea that BH accretion and SF rates are correlated, but occur with different variability time scales. The presence of a positive correlation between $N_H$ and host $M_*$ suggests that the X-ray $N_H$ is not entirely due to the circum-nuclear obscuring torus, but may also include a contribution from the host galaxy.

We estimate the accretion rates onto the supermassive black holes that power 20 of the highest-redshift quasars, at z>5.8, including the quasar with the highest redshift known to date -- ULAS J1120 at z=7.09. The analysis is based on the observed (rest-frame) optical luminosities and reliable "virial" estimates of the BH masses of the quasars, and utilizes scaling relations derived from thin accretion disk theory. The mass accretion rates through the postulated disks cover a wide range, dM_disk/dt~4-190 Msol/yr, with most of the objects (80\%) having dM_disk/dt~10-65 Msol/yr, confirming the Eddington-limited nature of the accretion flows. By combining our estimates of dM_disk/dt with conservative, lower limits on the bolometric luminosities of the quasars, we investigate which alternative values of \eta best account for all the available data. We find that the vast majority of quasars (~85\%) can be explained with radiative efficiencies in the range \eta~0.03-0.3, with a median value close to the commonly assumed \eta=0.1. Within this range, we obtain conservative estimates of \eta>0.14 for ULAS J1120 and SDSS J0100 (at z=6.3), and of >0.19 for SDSS J1148 (at z=6.41; assuming their BH masses are accurate). The implied accretion timescales are generally in the range t_acc = M_BH / dM_BH/dt ~0.1-1 Gyr, suggesting that most quasars could have had ~1-10 mass e-foldings since BH seed formation. Our analysis therefore demonstrates that the available luminosities and masses for the highest-redshift quasars can be explained self-consistently within the thin, radiatively efficient accretion disk paradigm. Episodes of radiatively inefficient, "super-critical" accretion may have occurred at significantly earlier epochs (i.e., z~10 or above).

The population of supermassive black holes (SMBHs) is composed by quiescent SMBHs, such as those seen in local galaxies including the Milky Way's, and active ones, resulting in quasars and active galactic nuclei (AGN). Outside our neighbourhood, all the information we have on SMBHs is derived from quasars and AGN, giving us a partial view. We study the evolution of the SMBH population, total and active, by the continuity equation, backwards in time from z=0 to z=4. Type-1 and type-2 AGN are differentiated in the model on the basis of the Eddington ratio distribution, chosen on the basis of observational estimates. The duty cycle is obtained by matching the luminosity function of quasars, and the average radiative efficiency is the only free parameter in the model. For higher radiative efficiencies (>~0.07) a large fraction of the SMBH population, most of them quiescent, must already be in place by z=4. For lower radiative efficiencies (~0.05), the duty cycle increases with the redshift and the SMBH population evolves dramatically since z=4. The mass function of active SMBHs does not depend on the choice of the radiative efficiency or of the local SMBH mass function, but it is mainly determined by the quasar luminosity function, once the Eddington ratio distribution is fixed. Only a direct measurement of the total BHMF at redshifts z>~2 could break these degeneracies giving important constraints on the average radiative efficiency. Focusing on type-1 AGN, for which observational estimates of the mass function and Eddington ratio distribution exist at various redshift, models with lower radiative efficiencies reproduce better the high-mass end of the mass function at high-z, but they tend to over-predict it at low-z, and vice-versa for models with higher radiative efficiencies.

By using cosmological hydrodynamical simulations we study the effect of supernova (SN) and active galactic nuclei (AGN) feedback on the mass transport of gas on to galactic nuclei and the black hole (BH) growth down to redshift z~6. We study the BH growth in relation with the mass transport processes associated with gravity and pressure torques, and how they are modified by feedback. Cosmological gas funelled through cold flows reaches the galactic outer region close to free-fall. Then torques associated to pressure triggered by gas turbulent motions produced in the circum-galactic medium by shocks and explosions from SNe are the main source of mass transport beyond the central ~ 100 pc. Due to high concentrations of mass in the central galactic region, gravitational torques tend to be more important at high redshift. The combined effect of almost free-falling material and both gravity and pressure torques produces a mass accretion rate of order ~ 1 M_sun/yr at ~ pc scales. In the absence of SN feedback, AGN feedback alone does not affect significantly either star formation or BH growth until the BH reaches a sufficiently high mass of $\sim 10^6$ M_sun to self-regulate. SN feedback alone, instead, decreases both stellar and BH growth. Finally, SN and AGN feedback in tandem efficiently quench the BH growth, while star formation remains at the levels set by SN feedback alone due to the small final BH mass, ~ few 10^5 M_sun. SNe create a more rarefied and hot environment where energy injection from the central AGN can accelerate the gas further.

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

The interplay between cosmic gas accretion on to galaxies and galaxy mergers drives the observed morphological diversity of galaxies. By comparing the state-of-the-art hydrodynamical cosmological simulations Horizon-AGN and Horizon-noAGN, we unambiguously identify the critical role of active galactic nuclei (AGN) in setting up the correct galaxy morphology for the massive end of the population. With AGN feedback, typical kinematic and morpho-metric properties of galaxy populations as well as the galaxy-halo mass relation are in much better agreement with observations. Only AGN feedback allows massive galaxies at the centre of groups and clusters to become ellipticals, while without AGN feedback those galaxies reform discs. It is the merger-enhanced AGN activity that is able to freeze the morphological type of the post-merger remnant by durably quenching its quiescent star formation. Hence morphology is shown to be driven not only by mass but also by the nature of cosmic accretion: at constant galaxy mass, ellipticals are galaxies that are mainly assembled through mergers, while discs are preferentially built from the in situ star formation fed by smooth cosmic gas infall.

The rate at which matter flows into a galactic nucleus during early phases of galaxy evolution can sometimes exceed the Eddington limit of the growing central black hole by several orders of magnitude. We discuss the necessary conditions for the black hole to actually accrete this matter at such a high rate, and consider the observational appearance and detectability of a hyperaccreting black hole. In order to be accreted at a hyper-Eddington rate, the infalling gas must have a sufficiently low angular momentum. Although most of the gas is accreted, a significant fraction accumulates in an optically thick envelope with luminosity $\sim L_{\rm Edd}$, probably pierced by jets of much higher power. If $\dot M > 10^3 \dot M_{\rm Edd}$, the envelope spectrum resembles a blackbody with a temperature of a few thousand K, but for lower (but still hyper-Eddington) accretion rates the spectrum becomes a very dilute and hard Wien spectrum. We consider the likelihood of various regimes of hyperaccretion, and discuss its possible observational signatures.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

Gravitational waves (GWs) provide a revolutionary tool to investigate yet unobserved astrophysical objects. Especially the first stars, which are believed to be more massive than present-day stars, might be indirectly observable via the merger of their compact remnants. We develop a self-consistent, cosmologically representative, semi-analytical model to simulate the formation of the first stars. By extrapolating binary stellar-evolution models at 10% solar metallicity to metal-free stars, we track the individual systems until the coalescence of the compact remnants. We estimate the contribution of primordial stars to the merger rate density and to the detection rate of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). Owing to their higher masses, the remnants of primordial stars produce strong GW signals, even if their contribution in number is relatively small. We find a probability of $\gtrsim1\%$ that the current detection GW150914 is of primordial origin. We estimate that aLIGO will detect roughly 1 primordial BH-BH merger per year for the final design sensitivity, although this rate depends sensitively on the primordial initial mass function (IMF). Turning this around, the detection of black hole mergers with a total binary mass of $\sim 300\,\mathrm{M}_\odot$ would enable us to constrain the primordial IMF.

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.

Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, we explore the impact of trace amounts of metals and dust enrichment. We perform three dimensional cosmological simulations for two halos of $\rm > 10^7~M_{\odot}$ with $\rm Z/Z_{\odot}= 10^{-4}-10^{-6}$ illuminated by an intense Lyman Werner flux of $\rm J_{21}=10^5$. Our results show that initially the collapse proceeds isothermally with $\rm T \sim 8000$ K but dust cooling becomes effective at densities of $\rm 10^{8}-10^{12} ~cm^{-3}$ and brings the gas temperature down to a few 100-1000 K for $\rm Z/Z_{\odot} \geq 10^{-6}$. No gravitationally bound clumps are found in $\rm Z/Z_{\odot} \leq 10^{-5}$ cases by the end of our simulations in contrast to the case with $\rm Z/Z_{\odot} = 10^{-4}$. Large inflow rates of $\rm \geq 0.1~M_{\odot}/yr$ are observed for $\rm Z/Z_{\odot} \leq 10^{-5}$ similar to a zero-metallicity case while for $\rm Z/Z_{\odot} = 10^{-4}$ the inflow rate starts to decline earlier due to the dust cooling and fragmentation. For given large inflow rates a central star of $\rm \sim 10^4~M_{\odot}$ may form for $\rm Z/Z_{\odot} \leq 10^{-5}$.

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 growth of the first super massive black holes (SMBHs) at z > 6 is still a major challenge for theoretical models. If it starts from black hole (BH) remnants of Population III stars (light seeds with mass ~ 100 Msun) it requires super-Eddington accretion. An alternative route is to start from heavy seeds formed by the direct collapse of gas onto a ~ 10^5 Msun BH. Here we investigate the relative role of light and heavy seeds as BH progenitors of the first SMBHs. We use the cosmological, data constrained semi-analytic model GAMETE/QSOdust to simulate several independent merger histories of z > 6 quasars. Using physically motivated prescriptions to form light and heavy seeds in the progenitor galaxies, we find that the formation of a few heavy seeds (between 3 and 30 in our reference model) enables the Eddington-limited growth of SMBHs at z > 6. This conclusion depends sensitively on the interplay between chemical, radiative and mechanical feedback effects, which easily erase the conditions that allow the suppression of gas cooling in the low metallicity gas (Z < Zcr and JLW > Jcr). We find that heavy seeds can not form if dust cooling triggers gas fragmentation above a critical dust-to-gas mass ratio (D > Dcr). In addition, the relative importance of light and heavy seeds depends on the adopted mass range for light seeds, as this dramatically affects the history of cold gas along the merger tree, by both SN and AGN-driven winds.

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

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.

The dynamics of massive black holes (BHs) in galaxy mergers is a rich field of research that has seen much progress in recent years. In this contribution we briefly review the processes describing the journey of BHs during mergers, from the cosmic context all the way to when BHs coalesce. If two galaxies each hosting a central BH merge, the BHs would be dragged towards the center of the newly formed galaxy. If/when the holes get sufficiently close, they coalesce via the emission of gravitational waves. How often two BHs are involved in galaxy mergers depends crucially on how many galaxies host BHs and on the galaxy merger history. It is therefore necessary to start with full cosmological models including BH physics and a careful dynamical treatment. After galaxies have merged, however, the BHs still have a long journey until they touch and coalesce. Their dynamical evolution is radically different in gas-rich and gas-poor galaxies, leading to a sort of "dichotomy" between high-redshift and low-redshift galaxies, and late-type and early-type, typically more massive galaxies.

Scaling relations between central black hole (BH) mass and host galaxy properties are of fundamental importance to studies of BH and galaxy evolution throughout cosmic time. Here we investigate the relationship between BH mass and host galaxy total stellar mass using a sample of 262 broad-line active galactic nuclei (AGN) in the nearby Universe (z < 0.055), as well as 79 galaxies with dynamical BH masses. The vast majority of our AGN sample is constructed using Sloan Digital Sky Survey spectroscopy and searching for Seyfert-like narrow-line ratios and broad H-alpha emission. BH masses are estimated using standard virial techniques. We also include a small number of dwarf galaxies with total stellar masses M_stellar < 10^9.5 Msun and a sub-sample of the reverberation-mapped AGNs. Total stellar masses of all 341 galaxies are calculated in the most consistent manner feasible using color-dependent mass-to-light ratios. We find a clear correlation between BH mass and total stellar mass for the AGN host galaxies, with M_BH proportional to M_stellar, similar to that of early-type galaxies with dynamically-detected BHs. However, the relation defined by the AGNs has a normalization that is lower by more than an order of magnitude, with a BH-to-total stellar mass fraction of M_BH/M_stellar ~ 0.025% across the stellar mass range 10^8 < M_stellar/Msun < 10^12. This result has significant implications for studies at high redshift and cosmological simulations in which stellar bulges cannot be resolved.

Massive Black Hole (MBH) seeds at redshift $z \gtrsim 10$ are now thought to be key ingredients to explain the presence of the super-massive ($10^{9-10} \, \mathrm{M_{\odot}}$) black holes in place $ < 1 \, \mathrm{Gyr}$ after the Big Bang. Once formed, massive seeds grow and emit copious amounts of radiation by accreting the left-over halo gas; their spectrum can then provide crucial information on their evolution. By combining radiation-hydrodynamic and spectral synthesis codes, we simulate the time-evolving spectrum emerging from the host halo of a MBH seed with initial mass $10^5 \, \mathrm{M_{\odot}}$, assuming both standard Eddington-limited accretion, or slim accretion disks, appropriate for super-Eddington flows. The emission occurs predominantly in the observed infrared-submm ($1-1000 \, \mathrm{\mu m}$) and X-ray ($0.1 - 100 \, \mathrm{keV}$) bands. Such signal should be easily detectable by JWST around $\sim 1 \, \mathrm{\mu m}$ up to $z \sim 25$, and by ATHENA (between $0.1$ and $10 \, \mathrm{keV}$, up to $z \sim 15$). Ultra-deep X-ray surveys like the Chandra Deep Field South could have already detected these systems up to $z \sim 15$. Based on this, we provide an upper limit for the $z \gtrsim 6$ MBH mass density of $\rho_{\bullet} \lesssim 2.5 \times 10^{2} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ assuming standard Eddington-limited accretion. If accretion occurs in the slim disk mode the limits are much weaker, $\rho_{\bullet} \lesssim 7.6 \times 10^{3} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ in the most constraining case.

Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics.

High redshift quasars at z>6 have masses up to ~$10^9$ M$_\odot$. One of the pathways to their formation includes direct collapse of gas, forming a supermassive star, precursor of the black hole seed. The conditions for direct collapse are more easily achievable in metal-free haloes, where atomic hydrogen cooling operates and molecular hydrogen (H2) formation is inhibited by a strong external UV flux. Above a certain value of UV flux (J_crit), the gas in a halo collapses isothermally at ~$10^4$ K and provides the conditions for supermassive star formation. However, H2 can self-shield, reducing the effect of photodissociation. So far, most numerical studies used the local Jeans length to calculate the column densities for self-shielding. We implement an improved method for the determination of column densities in 3D simulations and analyse its effect on the value of J_crit. This new method captures the gas geometry and velocity field and enables us to properly determine the direction-dependent self-shielding factor of H2 against photodissociating radiation. We find a value of J_crit that is a factor of two smaller than with the Jeans approach (~2000 J_21 vs. ~4000 J_21). The main reason for this difference is the strong directional dependence of the H2 column density. With this lower value of J_crit, the number of haloes exposed to a flux >J_crit is larger by more than an order of magnitude compared to previous studies. This may translate into a similar enhancement in the predicted number density of black hole seeds.

The growth of a supermassive black hole (BH) is determined by how much gas the host galaxy is able to feed it, which in turn is controlled by the cosmic environment, through galaxy mergers and accretion of cosmic flows that time how galaxies obtain their gas, but also by internal processes in the galaxy, such as star formation and feedback from stars and the BH itself. In this paper, we study the growth of a 10^12 Msun halo at z=2, which is the progenitor of al group of galaxies at z=0, and of its central BH by means of a high-resolution zoomed cosmological simulation, the Seth simulation. We study the evolution of the BH driven by the accretion of cold gas in the galaxy, and explore the efficiency of the feedback from supernovae (SNe). For a relatively inefficient energy input from SNe, the BH grows at the Eddington rate from early times, and reaches self-regulation once it is massive enough. We find that at early cosmic times z>3.5, efficient feedback from SNe forbids the formation of a settled disc as well as the accumulation of dense cold gas in the vicinity of the BH and starves the central compact object. As the galaxy and its halo accumulate mass, they become able to confine the nuclear inflows provided by major mergers and the BH grows at a sustained near-to-Eddington accretion rate. We argue that this mechanism should be ubiquitous amongst low-mass galaxies, corresponding to galaxies with a stellar mass below <10^9 Msun in our simulations.

The observational evidence that Super-Massive Black Holes ($M_{\bullet} \sim 10^{9-10} \, \mathrm{M_{\odot}}$) are already in place less than $1 \, \mathrm{Gyr}$ after the Big Bang poses stringent time constraints on the growth efficiency of their seeds. Among proposed possibilities, the formation of massive ($\sim 10^{3-6} \, \mathrm{M_{\odot}}$) seeds and/or the occurrence of super-Eddington ($\dot{M}>\dot{M}_{Edd}$) accretion episodes may contribute to the solution of this problem. In this work we analytically and numerically investigate the accretion flow onto high-redshift ($z \sim 10$) black holes to understand the physical requirements favoring rapid and efficient growth. Our model identifies a "feeding-dominated" accretion regime and a "feedback-limited" one, the latter being characterized by intermittent (duty cycles ${\cal D} \lesssim 0.5$) and inefficient growth, with recurring outflow episodes. We find that low-mass seeds ($\lesssim 10^{3-4} \, \mathrm{M_{\odot}}$) evolve in the feedback-limited regime, while more massive seeds ($\gtrsim 10^{5-6} \, \mathrm{M_{\odot}}$) grow very rapidly as they are found in the feeding-dominated regime. In addition to the standard accretion model with a fixed matter-energy conversion factor ($\epsilon = 0.1$), we have also explored slim disk models, appropriate for super-Eddington accretion, where radiation is trapped in the disk and the radiative efficiency is reduced ($\epsilon \lesssim 0.04$), which may ensure a continuous growth with $\dot{M} \gg \dot{M}_{Edd}$ (up to $\sim 300\dot{M}_{Edd}$ in our simulations). Under these conditions, outflows play a negligible role and a black hole can accrete $80\%-100\%$ of the gas mass of the host halo ($\sim 10^7 \, \mathrm{M_{\odot}}$) in $\sim 10 \, \mathrm{Myr}$, while in feedback-limited systems we predict that black holes can accrete only up to $\sim 15\%$ of the available mass.

Supermassive black holes are not only common in the present-day galaxies, but billion solar masses black holes also powered $z\geq 6$ quasars. One efficient way to form such black holes is the collapse of a massive primordial gas cloud into a so-called direct collapse black hole. The main requirement for this scenario is the presence of large accretion rates of $\rm \geq 0.1~M_{\odot}/yr$ to form a supermassive star. It is not yet clear how and under what conditions such accretion rates can be obtained. The prime aim of this work is to determine the mass accretion rates under non-isothermal collapse conditions. We perform high resolution cosmological simulations for three primordial halos of a few times $\rm 10^7~M_{\odot}$ illuminated by an external UV flux, $\rm J_{21}=100-1000$. We find that a rotationally supported structure of about parsec size is assembled, with an aspect ratio between $\rm 0.25 - 1$ depending upon the thermodynamical properties. Rotational support, however, does not halt collapse, and mass inflow rates of $\rm \sim 0.1~M_{\odot}/yr$ can be obtained in the presence of even a moderate UV background flux of strength $\rm J_{21} \geq 100$. To assess whether such large accretion rates can be maintained over longer time scales, we employed sink particles, confirming the persistence of accretion rates of $\rm \sim 0.1~M_{\odot}/yr$. We propose that complete isothermal collapse and molecular hydrogen suppression may not always be necessary to form supermassive stars, precursors of black hole seeds. Sufficiently high inflow rates can be obtained for UV flux $\rm J_{21}=500-1000$, at least for some cases. This value brings the estimate of the abundance of direct collapse black hole seeds closer to that high redshift quasars.

We introduce a force correction term to better model the dynamical friction (DF) experienced by a supermassive black hole (SMBH) as it orbits within its host galaxy. This new approach accurately follows the orbital decay of a SMBH and drastically improves over commonly used advection methods. The force correction introduced here naturally scales with the force resolution of the simulation and converges as resolution is increased. In controlled experiments we show how the orbital decay of the SMBH closely follows analytical predictions when particle masses are significantly smaller than that of the SMBH. In a cosmological simulation of the assembly of a small galaxy, we show how our method allows for realistic black hole orbits. This approach overcomes the limitations of the advection scheme, where black holes are rapidly and artificially pushed toward the halo center and then forced to merge, regardless of their orbits. We find that SMBHs from merging dwarf galaxies can spend significant time away from the center of the remnant galaxy. Improving the modeling of SMBH orbital decay will help in making robust predictions of the growth, detectability, and merger rates of SMBHs, especially at low galaxy masses or at high redshift.

We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.

Short-lived intermittent phases of super-critical (super-Eddington) growth, coupled with star formation via positive feedback, may account for early growth of massive black holes (MBH) and coevolution with their host spheroids. We estimate the possible growth rates and duty cycles of these episodes, both assuming slim accretion disk solutions, and adopting the results of recent numerical simulations. The angular momentum of gas joining the accretion disk determines the length of the accretion episodes, and the final mass a MBH can reach. The latter can be related to the gas velocity dispersion, and in galaxies with low-angular momentum gas the MBH can get to a higher mass. When the host galaxy is able to sustain inflow rates at 1-100 msunyr, replenishing and circulation lead to a sequence of short (~1e4-1e7 years), heavily obscured accretion episodes that increase the growth rates, with respect to an Eddington-limited case, by several orders of magnitude. Our model predicts that the ratio of MBH accretion rate to star formation rate is 1e2 or higher, leading, at early epochs, to a ratio of MBH to stellar mass higher than the "canonical" value of ~1e-3, in agreement with current observations. Our model makes specific predictions that long-lived super-critical accretion occurs only in galaxies with copious low-angular momentum gas, and in this case the MBH is more massive at fixed velocity dispersion.

We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (~10 pc) and temporal (~1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secular-related causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with $low$ angular momentum, but the gas that $loses$ angular momentum.

The formation of black holes inevitably affects the distribution of dark and baryonic matter in their vicinity, leading to an enhancement of the dark matter density, called spike, and if dark matter is made of WIMPs, to a strong enhancement of the dark matter annihilation rate. Spikes at the center of galaxies like the Milky Way are efficiently disrupted by baryonic processes, but mini-spikes can form and survive undisturbed at the center of dwarf spheroidal galaxies. We show that Fermi LAT satellite data allow to set very stringent limits on the existence of mini-spikes in dwarf galaxies: for thermal WIMPs with mass between 100 GeV and 1 TeV, we obtain a maximum black hole mass between 100 and 1000 $M_\odot$, ruling out black holes masses extrapolated from the M-{\sigma} relationship in a large region of the parameter space. We also performed Monte Carlo simulations of merger histories of black holes in dwarf spheroidals in a scenario where black holes form from the direct collapse of primordial gas in early halos, and found that this specific formation scenario is incompatible at the 84% CL with dark matter being in the form of thermal WIMPs.

The peculiar QSO J0927+2943 shows multiple sets of emission lines in its optical spectrum. This signature has been interpreted as the relative motion between a black hole, either recoiling or bound in a binary system, and its host galaxy, or as a superposition of two galaxies along the line of sight. In order to test these scenarios, we have collected 2mm CO(2-1) observations using the IRAM Plateau de Bure Interferometer, and optical images and spectroscopy at the Calar Alto observatory. Together with archival HST images, these data provide unique insights on the nature of this system. The recoiling/binary black hole scenarios are ruled out by the clear detection of a galactic-scale molecular gas reservoir at the same redshift of the QSO broad lines. The observations presented here also disfavour the superposition model, although with less constraints. Thus, the origin of the second, bright set of narrow emission lines in J0927+2943 is still unknown.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

In this paper we explore a possible route of black hole seed formation that appeal to a model by Davies, Miller & Bellovary who considered the case of the dynamical collapse of a dense cluster of stellar black holes subjected to an inflow of gas. Here, we explore this case in a broad cosmological context. The working hypotheses are that (i) nuclear star clusters form at high redshifts in pre-galactic discs hosted in dark matter halos, providing a suitable environment for the formation of stellar black holes in their cores, (ii) major central inflows of gas occur onto these clusters due to instabilities seeded in the growing discs and/or to mergers with other gas-rich halos, and that (iii) following the inflow, stellar black holes in the core avoid ejection due to the steepening to the potential well, leading to core collapse and the formation of a massive seed of $<~ 1000\, \rm M_\odot$. We simulate a cosmological box tracing the build up of the dark matter halos and there embedded baryons, and explore cluster evolution with a semi-analytical model. We show that this route is feasible, peaks at redshifts $z <~ 10$ and occurs in concomitance with the formation of seeds from other channels. The channel is competitive relative to others, and is independent of the metal content of the parent cluster. This mechanism of gas driven core collapse requires inflows with masses at least ten times larger than the mass of the parent star cluster, occurring on timescales shorter than the evaporation/ejection time of the stellar black holes from the core. In this respect, the results provide upper limit to the frequency of this process.

Supermassive black holes (BH) accrete gas from their surroundings and coalesce with companions during galaxy mergers, and both processes change the BH mass and spin. By means of high-resolution hydrodynamical simulations of galaxies, either idealised or embedded within the cosmic web, we explore the effects of interstellar gas dynamics and external perturbations on BH spin evolution. All these physical quantities were evolved on-the-fly in a self-consistent manner. We use a `maximal' model to describe the turbulence induced by stellar feedback to highlight its impact on the angular momentum of the gas accreted by the BH. Periods of intense star formation are followed by phases where stellar feedback drives large-scale outflows and hot bubbles. We find that BH accretion is synchronised with star formation, as only when gas is cold and dense do both processes take place. During such periods, gas motion is dominated by consistent rotation. On the other hand, when stellar feedback becomes substantial, turbulent motion randomises gas angular momentum. However BH accretion is strongly suppressed in that case, as cold and dense gas is lacking. In our cosmological simulation, at very early times (z>6), the galactic disc has not yet settled and no preferred direction exists for the angular momentum of the accreted gas, so the BH spin remains low. As the gas settles into a disc (6>z>3), the BH spin then rapidly reaches its maximal value. At lower redshifts (z<3), even when galaxy mergers flip the direction of the angular momentum of the accreted gas, causing it to counter-rotate, the BH spin magnitude only decreases modestly and temporarily. Should this be a typical evolution scenario for BH, it potentially has dramatic consequences regarding their origin and assembly, as accretion on maximally spinning BH embedded in thin Shakura-Sunyaev disc is significantly reduced.

Supermassive black holes (BH) at the centres of galaxies can rapidly change their mass and spin by gas accretion and mergers. Using hydrodynamical cosmological simulations, with prescriptions for BH growth and feedback from Active Galactic Nuclei, we study how the evolution of BH mass growth is driven by gas accretion and mergers. Using a semi-analytical approach to evolve spins, we also highlight the mechanisms responsible for driving the magnitude and the direction of spins as a function of cosmic time. We find that in the high-redshift universe galaxies maintain large values of gas accretion onto BHs, which therefore is the main driver of their mass and spin evolution. Sustained accretion of cold gas at high-redshift tends to align BH spins with the angular momentum of the surrounding gas and maximise their magnitude. Conversely, at low redshift, as BHs get more massive and galaxies more gas-poor, the contribution from binary coalescences to the total BH mass growth increases, especially at the high-mass end, and tends to decrease the magnitude of spins and change their direction.

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

Using a set of zoomed-in cosmological simulations of high-redshift progenitors of massive galaxies, we isolate and trace the history of gas that is accreted by central supermassive black holes. We determine the origins of the accreted gas, in terms of whether it entered the galaxy during a merger event or was smoothly accreted. Furthermore, we designate whether the smoothly accreted gas is accreted via a cold flow or is shocked upon entry into the halo. For moderate-mass (10^6 - 10^7 Msun) black holes at z ~ 4, there is a preference to accrete cold flow gas than gas of shocked or merger origin. However, this result is a consequence of the fact that the entire galaxy has a higher fraction of gas from cold flows. In general, each black hole tends to accrete the same fractions of smooth- and merger-accreted gas as is contained in its host galaxy, suggesting that once gas enters a halo it becomes well-mixed, and its origins are erased. We find that the angular momentum of the gas upon halo entry is a more important factor; black holes preferentially accrete gas that had low angular momentum when it entered the galaxy, regardless of whether it was accreted smoothly or through mergers.

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

We constrain the total accreted mass density in supermassive black holes at z>6, inferred via the upper limit derived from the integrated X-ray emission from a sample of photometrically selected galaxy candidates. Studying galaxies obtained from the deepest Hubble Space Telescope images combined with the Chandra 4 Msec observations of the Chandra Deep Field South, we achieve the most restrictive constraints on total black hole growth in the early Universe. We estimate an accreted mass density <1000Mo Mpc^-3 at z~6, significantly lower than the previous predictions from some existing models of early black hole growth and earlier prior observations. These results place interesting constraints on early black growth and mass assembly by accretion and imply one or more of the following: (1) only a fraction of the luminous galaxies at this epoch contain active black holes; (2) most black hole growth at early epochs happens in dusty and/or less massive - as yet undetected - host galaxies; (3) there is a significant fraction of low-z interlopers in the galaxy sample; (4) early black hole growth is radiatively inefficient, heavily obscured and/or is due to black hole mergers as opposed to accretion or (5) the bulk of the black hole growth occurs at late times. All of these possibilities have important implications for our understanding of high redshift seed formation models.

Massive black holes (MBHs) in contrast to stellar mass black holes are expected to substantially change their properties over their lifetime. MBH masses increase by several order of magnitude over the Hubble time, as illustrated by Soltan's argument. MBH spins also must evolve through the series of accretion and mergers events that grow the MBH's masses. We present a simple model that traces the joint evolution of MBH masses and spins across cosmic time. Our model includes MBH-MBH mergers, merger-driven gas accretion, stochastic fueling of MBHs through molecular cloud capture, and a basic implementation of accretion of recycled gas. This approach aims at improving the modeling of low-redshift MBHs and AGN, whose properties can be more easily estimated observationally. Despite the simplicity of the model, it captures well the global evolution of the MBH population from z\sim6 to today. Under our assumptions, we find that the typical spin and radiative efficiency of MBHs decrease with cosmic time because of the higher incidence of stochastic processes in gas-rich galaxies and MBH-MBH mergers in gas-poor galaxies. At z=0 the spin distribution in gas-poor galaxies peaks at spins 0.4-0.8, and it is not strongly mass dependent. MBHs in gas-rich galaxies have a more complex evolution, with low-mass MBHs at low redshift having low spins, and spins increasing at larger masses and redshifts. We also find that at z>1 MBH spins are on average highest in high luminosity AGN, while at lower redshifts these differences disappear.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

In powerful radio-quiet active galactic nuclei (AGN), black holes heavier than one billion solar masses form at a redshift ~1.5-2. Supermassive black holes in jetted radio-loud AGN seems to form earlier, at a redshift close to 4. The ratio of active radio-loud to radio-quiet AGN hosting heavy black holes is therefore a rather a strong function of redshift. We report on some recent evidence supporting this conclusion, gathered from the Burst Alert Telescope (BAT, onboard Swift) and by the Large Area Telescope (LAT, onboard Fermi). We suggest that the more frequent occurrence of relativistic jets in the most massive black holes at high redshifts, compared to later times, could be due to the average black hole spin being greater in the distant past, or else to the jet helping a fast accretion rate (or some combination of the two scenarios). We emphasize that the large total accretion efficiency of rapidly spinning black holes inhibits a fast growth, unless a large fraction of the available gravitational energy of the accreted mass is not converted into radiation, but used to form and maintain a powerful jet.

We explore how the co-evolution of massive black holes (MBHs) and galaxies is affected by environmental effects, addressing in particular MBHs hosted in the central galaxies of clusters (we will refer to these galaxies in general as 'CGs'). Recently the sample of MBHs in CGs with dynamically measured masses has increased, and it has been suggested that these MBH masses (M_BH) deviate from the expected correlations with velocity dispersion (sigma) and mass of the bulge (M_bulge) of the host galaxy: MBHs in CGs appear to be `over-massive'. This discrepancy is more pronounced when considering the M_BH-sigma relation than the M_BH-M_bulge one. We show that this behavior stems from a combination of two natural factors, (i) that CGs experience more mergers involving spheroidal galaxies and their MBHs, and (ii) that such mergers are preferentially gas-poor. We use a combination of analytical and semi-analytical models to investigate the MBH-galaxy co-evolution in different environments and find that the combination of these two factors explains the trends observed in current data-sets.

Mass and spin are the only two parameters needed to completely characterize black holes in General Relativity. However, the interaction between black holes and their environment is where complexity lies, as the relevant physical processes occur over a large range of scales. That is particularly relevant in the case of super-massive black holes (SMBHs), hosted in galaxy centers, and surrounded by swirling gas and various generations of stars. These compete with the SMBH for gas consumption and affect both dynamics and thermodynamics of the gas itself. How the behavior of such fiery environment influence the angular momentum of the gas accreted onto SMBHs, and, hence, black-hole spins is uncertain. We explore the interaction between SMBHs and their environment via first 3D sub-parsec resolution simulations (ranging from 0.1 pc to 1 kpc scales) that study the evolution of the SMBH spin by including the effects of star formation, stellar feedback, radiative transfer, and metal pollution according to the proper stellar yields and lifetimes. This approach is crucial to investigate the impact of star formation processes and feedback effects on the angular momentum of the material that could accrete on the central hole. We find that star formation and feedback mechanisms can locally inject significant amounts of entropy in the surrounding medium, and impact on the inflow inclination angles and Eddington fractions. As a consequence, the resulting trends show upper-intermediate equilibrium values for the spin parameter, a, of about 0.6 - 0.9, corresponding to radiative efficiencies \epsilon = 9% - 15%. These results suggest that star formation feedback taking place in the circum-nuclear disk during the in-fall cannot induce alone very strong chaotic trends in the gas flow, quite independently from the different numerical parameters.

We place firm upper limits on the global accretion history of massive black holes at z>5 from the recently measured unresolved fraction of the cosmic X-ray background. The maximum allowed unresolved intensity observed at 1.5 keV implies a maximum accreted-mass density onto massive black holes of rho_acc < 1.4E4 M_sun Mpc^{-3} for z>5. Considering the contribution of lower-z AGNs, the value reduces to rho_acc < 0.66E4 M_sun Mpc^{-3}. The tension between the need for the efficient and rapid accretion required by the observation of massive black holes already in place at z>7 and the strict upper limit on the accreted mass derived from the X-ray background may indicate that black holes are rare in high redshift galaxies, or that accretion is efficient only for black holes hosted by rare galaxies.

Massive black holes in galactic nuclei vary their mass M and spin vector J due to accretion. In this study we relax, for the first time, the assumption that accretion can be either chaotic, i.e. when the accretion episodes are randomly and isotropically oriented, or coherent, i.e. when they occur all in a preferred plane. Instead, we consider different degrees of anisotropy in the fueling, never confining to accretion events on a fixed direction. We follow the black hole growth evolving contemporarily mass, spin modulus a and spin direction. We discover the occurrence of two regimes. An early phase (M <~ 10 million solar masses) in which rapid alignment of the black hole spin direction to the disk angular momentum in each single episode leads to erratic changes in the black hole spin orientation and at the same time to large spins (a ~ 0.8). A second phase starts when the black hole mass increases above >~ 10 million solar masses and the accretion disks carry less mass and angular momentum relatively to the hole. In the absence of a preferential direction the black holes tend to spin-down in this phase. However, when a modest degree of anisotropy in the fueling process (still far from being coherent) is present, the black hole spin can increase up to a ~ 1 for very massive black holes (M >~ 100 million solar masses), and its direction is stable over the many accretion cycles. We discuss the implications that our results have in the realm of the observations of black hole spin and jet orientations.

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

The existence of massive black holes was postulated in the sixties, when the first quasars were discovered. In the late nineties their reality was proven beyond doubt, in the Milky way and a handful nearby galaxies. Since then, enormous theoretical and observational efforts have been made to understand the astrophysics of massive black holes. We have discovered that some of the most massive black holes known, weighing billions of solar masses, powered luminous quasars within the first billion years of the Universe. The first massive black holes must therefore have formed around the time the first stars and galaxies formed. Dynamical evidence also indicates that black holes with masses of millions to billions of solar masses ordinarily dwell in the centers of today's galaxies. Massive black holes populate galaxy centers today, and shone as quasars in the past; the quiescent black holes that we detect now in nearby bulges are the dormant remnants of this fiery past. In this review we report on basic, but critical, questions regarding the cosmological significance of massive black holes. What physical mechanisms lead to the formation of the first massive black holes? How massive were the initial massive black hole seeds? When and where did they form? How is the growth of black holes linked to that of their host galaxy? Answers to most of these questions are work in progress, in the spirit of these Reports on Progress in Physics.

I briefly outline recent theoretical developments on the formation of the first massive black holes (MBHs) that may grow into the population of MBHs powering quasars and inhabiting galactic centers today. I also touch upon possible observational tests that may give insights on what the properties of the first MBHs were.

The past 10 years have witnessed a change of perspective in the way astrophysicists think about massive black holes (MBHs), which are now considered to have a major role in the evolution of galaxies. This appreciation was driven by the realization that black holes of millions solar masses and above reside in the center of most galaxies, including the Milky Way. MBHs also powered active galactic nuclei known to exist just a few hundred million years after the Big Bang. Here, I summarize the current ideas on the evolution of MBHs through cosmic history, from their formation about 13 billion years ago to their growth within their host galaxies.

We explore the newly discovered "hangup-kick" effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.

Supermassive black holes (SMBHs) have been detected in the centers of most nearby massive galaxies. Galaxies today are the products of billions of years of galaxy mergers, but also billions of years of SMBH activity as active galactic nuclei (AGNs) that is connected to galaxy mergers. In this context, detection of AGN pairs should be relatively common. Observationally, however, dual AGN are scant, being just a few percent of all AGN. In this Letter we investigate the triggering of AGN activity in merging galaxies via a suite of high resolution hydrodynamical simulations. We follow the dynamics and accretion onto the SMBHs as they move from separations of tens of kiloparsecs to tens of parsecs. Our resolution, cooling and star formation implementation produce an inhomogeneous, multi-phase interstellar medium, allowing us to accurately trace star formation and accretion onto the SMBHs. We study the impact of gas content, morphology, and mass ratio, allowing us to study AGN activity and dynamics across a wide range of relevant conditions. We test when the two AGN are simultaneously detectable, for how long and at which separations. We find that strong dual AGN activity occurs during the late phases of the mergers, at small separations (<1-10 kpc) below the resolution limit of most surveys. Much of the SMBH accretion is not simultaneous, limiting the dual AGN fraction detectable through imaging and spectroscopy to a few percent, in agreement with observational samples.

In this paper, we compare the observationally derived black hole mass function (BHMF) of luminous broad-line quasars (BLQSOs) at 1<z<4.5 drawn from the Sloan Digital Sky Survey (SDSS) presented in Kelly et al. (2010), with models of merger driven BH growth in the context of standard hierarchical structure formation models. In the models, we explore two distinct black hole seeding prescriptions at the highest redshifts: "light seeds" - remnants of Population III stars and "massive seeds" that form from the direct collapse of pre-galactic disks. The subsequent merger triggered mass build-up of the black hole population is tracked over cosmic time under the assumption of a fixed accretion rate as well as rates drawn from the distribution derived by Merloni & Heinz. Our model snapshots are compared to the SDSS derived BHMFs of BLQSOs. Our key findings are that the duty cycle of SMBHs powering BLQSOs increases with increasing redshift for all models and models with Pop III remnants as black hole seeds are unable to fit the observationally derived BHMFs for BLQSOs, lending strong support for the massive seeding model (abridged).

We present results of simulations aimed at tracing the formation of nuclear star clusters (NCs) and black hole (BH) seeds, in a cosmological context. We focus on two mechanisms for the formation of BHs at high redshifts: as end-products of (1) Population III stars in metal free halos, and of (2) runaway stellar collisions in metal poor NCs. Our model tracks the chemical, radiative and mechanical feedback of stars on the baryonic component of the evolving halos. This procedure allows us to evaluate when and where the conditions for BH formation are met, and to trace the emergence of BH seeds arising from the dynamical channel, in a cosmological context. BHs start to appear already at z~30 as remnants of Population III stars. The efficiency of this mechanism begins decreasing once feedbacks become increasingly important. Around redshift z~15, BHs mostly form in the centre of mildly metal enriched halos inside dense NCs. The seed BHs that form along the two pathways have at birth a mass around 100-1000M\odot. The occupation fraction of BHs is a function of both halo mass and mass growth rate: at a given z, heavier and faster growing halos have a higher chance to form a native BH, or to acquire an inherited BH via merging of another system. With decreasing z, the probability of finding a BH shifts toward progressively higher mass halo intervals. This is due to the fact that, at later cosmic times, low mass systems rarely form a seed, and already formed BHs are deposited into larger mass systems due to hierarchical mergers. Our model predict that at z=0, all halos above 10^11M\odot should host a BH (in agreement with observational results), most probably inherited during their lifetime. Halos less massive then 10^9M\odot have a higher probability to host a native BH, but their occupation fraction decreases below 10%.

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal of providing insights to which galaxies do or do not host massive black holes. We adopt a novel approach to forming seed black holes in galaxy halos in cosmological SPH+N-body simulations. The formation of massive black hole seeds is dictated directly by the local gas density, temperature, and metallicity, and motivated by physical models of massive black hole formation. We explore seed black hole populations as a function of halo mass and redshift, and examine how varying the efficiency of massive black hole seed formation affects the relationship between black holes and their hosts. Seed black holes tend to form in halos with mass between 10^7 and 10^9 Msun, and the formation rate is suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find that the time of massive black hole formation and the occupation fraction of black holes are a function of the host halo mass. By z = 5, halos with mass M_halo > 3 x 10^9 Msun host massive black holes regardless of the efficiency of seed formation, while the occupation fraction for smaller halos increases with black hole formation efficiency. Our simulations explain why massive black holes are found in some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer in low-mass systems.

Motivated by recent observational results that focus on high redshift black holes, we explore the effect of scatter and observational biases on the ability to recover the intrinsic properties of the black hole population at high redshift. We find that scatter and selection biases can hide the intrinsic correlations between black holes and their hosts, with 'observable' subsamples of the whole population suggesting, on average, positive evolution even when the underlying population is characterized by no- or negative evolution. We create theoretical mass functions of black holes convolving the mass function of dark matter halos with standard relationships linking black holes with their hosts. Under these assumptions, we find that the local MBH - sigma correlation is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), overestimating the number density of all but the most massive black holes. Positive evolution or including scatter in the MBH - sigma correlation makes the discrepancy worse, as it further increases the number density of observable black holes. We notice that if the MBH - sigma correlation at z = 6 is steeper than today, then the mass function becomes shallower. This helps reproducing the mass function of z = 6 black holes proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of order 1/1000) host an active massive black hole at z = 6, or that most AGN are obscured, hindering their detection in optical surveys. Current measurements of the high redshift black hole mass function might be underestimating the density of low mass black holes if the active fraction or luminosity are a function of host or black hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper MBH - sigma relation at high redshift.

The formation of the first massive objects in the infant Universe remains impossible to observe directly and yet it sets the stage for the subsequent evolution of galaxies. While some black holes with masses > billion solar masses? have been detected in luminous quasars less than one billion years after the Big Bang, these individual extreme objects have limited utility in constraining the channels of formation of the earliest black holes. The initial conditions of black hole seed properties are quickly erased during the growth process. From deep, optimally stacked, archival X-ray observations, we measure the amount of black hole growth in z=6-8 galaxies (0.7-1 billion years after the Big Bang). Our results imply that black holes grow in tandem with their hosts throughout cosmic history, starting from the earliest times. We find that most copiously accreting black holes at these epochs are buried in significant amounts of gas and dust that absorb most radiation except for the highest energy X-rays. This suggests that black holes grow significantly more than previously thought during these early bursts, and due to obscuration they do not contribute to the re-ionization of the Universe with their ultraviolet emission.

In this paper, we examine if the properties of central black holes in galactic nuclei correlate with their host dark matter halo. We analyze the entire sample of galaxies where black hole mass, velocity dispersion, sigma, and asymptotic circular velocity, Vc, have all been measured. We fit M-sigma and M-Vc to a power law, and find that in both relationships the scatter and slope are similar. This model-independent analysis suggests that although the black hole masses are not uniquely determined by dark matter halo mass, when considered for the current sample as a whole, the M-Vc correlation may be as strong (or as weak) as M-sigma. Although the data are sparse, there appears to be more scatter in the correlation for both sigma and Vc at the low--mass end. This is not unexpected given our current understanding of galaxy and black hole assembly. In fact, there are several compelling reasons that account for this: (i) SMBH formation is likely less efficient in low-mass galaxies with large angular momentum content; (ii) SMBH growth is less efficient in low-mass disk galaxies that have not experienced major mergers; (iii) dynamical effects, such as gravitational recoil, increase scatter preferentially at the low-mass end. Therefore, the recent observational claim of the absence of central SMBHs in bulgeless, low mass galaxies, or deviations from the correlations defined by high-mass black holes in large galaxies today is, in fact, predicated by current models of black hole growth. We show how this arises as a direct consequence of the coupling between dark matter halos and central black holes at the earliest epochs.

We investigate the relative occurrence of radio--loud and radio-quiet quasars in the first billion years of the Universe, powered by black holes heavier than one billion solar masses. We consider the sample of high-redshfit blazars detected in the hard X-ray band in the 3-years all sky survey performed by the Burst Alert Telescope (BAT) onboard the Swift satellite. All the black holes powering these blazars exceed a billion solar mass, with accretion luminosities close to the Eddington limit. For each blazar pointing at us, there must be hundreds of similar sources (having black holes of similar masses) pointing elsewhere. This puts constraints on the density of billion solar masses black holes at high redshift (z>4), and on the relative importance of (jetted) radio-loud vs radio-quiet sources. We compare the expected number of high redshift radio--loud sources with the high luminosity radio-loud quasars detected in the Sloan Digital Sky Survey (SDSS), finding agreement up to z~3, but a serious deficit at z>3 of SDSS radio-loud quasars with respect to the expectations. We suggest that the most likely explanations for this disagreement are: i) the ratio of blazar to misaligned radio-sources decreases by an order of magnitude above z=3, possibly as a result of a decrease of the average bulk Lorentz factor; ii) the SDSS misses a large fraction of radio-loud sources at high redshifts, iii) the SDSS misses both radio-loud and radio-quiet quasars at high redshift, possibly because of obscuration or because of collimation of the optical-UV continuum in systems accreting near Eddington. These explanations imply very different number density of heavy black holes at high redshifts, that we discuss in the framework of the current ideas about the relations of dark matter haloes at high redshifts and the black hole they host.

Supermassive black holes presumably grow through numerous mergers throughout cosmic time. During each merger, supermassive black hole binaries are surrounded by a circumbinary accretion disk that imposes a significant (~1e4 G for a binary of 1e8 Msun) magnetic field. The motion of the binary through that field will convert the field energy to Poynting flux, with a luminosity ~1e43 erg/s (B/1e4 G)^2 (M/1e8 Msun)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. We find that the short timescales of many mergers will limit their detectability with most planned blind surveys to <1 per year over the whole sky, independent of the details of the emission process and flux distribution. Including an optimistic estimate for the radio flux makes detection even less likely, with <1 mergers per year over the whole sky. However, wide-field radio instruments may be able to localize systems identified in advance of merger by gravitational waves. Further, radio surveys may be able to detect the weaker emission produced by the binary's motion as it is modulated by spin-orbit precession and inspiral well in advance of merger.

The Suzaku AGN Spin Survey is designed to determine the supermassive black hole spin in six nearby active galactic nuclei (AGN) via deep Suzaku stares, thereby giving us our first glimpse of the local black hole spin distribution. Here, we present an analysis of the first target to be studied under the auspices of this Key Project, the Seyfert galaxy NGC 3783. Despite complexity in the spectrum arising from a multi-component warm absorber, we detect and study relativistic reflection from the inner accretion disk. Assuming that the X-ray reflection is from the surface of a flat disk around a Kerr black hole, and that no X-ray reflection occurs within the general relativistic radius of marginal stability, we determine a lower limit on the black hole spin of a > 0.88 (99% confidence). We examine the robustness of this result to the assumption of the analysis, and present a brief discussion of spin-related selection biases that might affect flux-limited samples of AGN.

LISA might detect gravitational waves from mergers of massive black hole binaries strongly lensed by intervening galaxies (Sereno et al. 2010). The detection of multiple gravitational lensing events would provide a new tool for cosmography. Constraints on cosmological parameters could be placed by exploiting either lensing statistics of strongly lensed sources or time delay measurements of lensed gravitational wave signals. These lensing methods do not need the measurement of the redshifts of the sources and the identification of their electromagnetic counterparts. They would extend cosmological probes to redshift z <= 10 and are then complementary to other lower or higher redshift tests, such as type Ia supernovae or cosmic microwave background. The accuracy of lensing tests strongly depends on the formation history of the merging binaries, and the related number of total detectable multiple images. Lensing amplification might also help to find the host galaxies. Any measurement of the source redshifts would allow to exploit the distance-redshift test in combination with lensing methods. Time-delay analyses might measure the Hubble parameter H_0 with accuracy of >= 10 km s^{-1}Mpc^{-1}. With prior knowledge of H_0, lensing statistics and time delays might constrain the dark matter density (delta Omega_M >= 0.08, due to parameter degeneracy). Inclusion of our methods with other available orthogonal techniques might significantly reduce the uncertainty contours for Omega_M and the dark energy equation of state.

Only a small fraction of local galaxies harbor an accreting black hole, classified as an active galactic nucleus (AGN). However, many stellar systems are plausibly expected to host black holes, from globular clusters to nuclear star clusters, to massive galaxies. The mere presence of stars in the vicinity of a black hole provides a source of fuel via mass loss of evolved stars. In this paper we assess the expected luminosities of black holes embedded in stellar systems of different sizes and properties, spanning a large range of masses. We model the distribution of stars and derive the amount of gas available to a central black hole through a geometrical model. We estimate the luminosity of the black holes under simple, but physically grounded, assumptions on the accretion flow. Finally we discuss the detectability of 'quiescent' black holes in the local Universe.

We perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism for black hole seed formation and growth, and find that satellite galaxies containing supermassive black hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like the Milky Way may host 5 -- 15 of these objects within their halos. The satellites that harbor black hole seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give a physical explanation for the origin and observed properties of some recently discovered off-nuclear ultraluminous X-ray sources such as HLX-1.

LISA should detect gravitational waves from tens to hundreds of systems containing black holes with mass in the range from 10 thousand to 10 million solar masses. Black holes in this mass range are not well constrained by current electromagnetic observations, so LISA could significantly enhance our understanding of the astrophysics of such systems. In this paper, we describe a framework for combining LISA observations to make statements about massive black hole populations. We summarise the constraints that LISA observations of extreme-mass-ratio inspirals might be able to place on the mass function of black holes in the LISA range. We also describe how LISA observations can be used to choose between different models for the hierarchical growth of structure in the early Universe. We consider four models that differ in their prescription for the initial mass distribution of black hole seeds, and in the efficiency of accretion onto the black holes. We show that with as little as 3 months of LISA data we can clearly distinguish between these models, even under relatively pessimistic assumptions about the performance of the detector and our knowledge of the gravitational waveforms.

The massive black holes we observe in galaxies today are the natural end-product of a complex evolutionary path, in which black holes seeded in proto-galaxies at high redshift grow through cosmic history via a sequence of mergers and accretion episodes. Electromagnetic observations probe a small subset of the population of massive black holes (namely, those that are active or those that are very close to us), but planned space-based gravitational-wave observatories such as the Laser Interferometer Space Antenna (LISA) can measure the parameters of ``electromagnetically invisible'' massive black holes out to high redshift. In this paper we introduce a Bayesian framework to analyze the information that can be gathered from a set of such measurements. Our goal is to connect a set of massive black hole binary merger observations to the underlying model of massive black hole formation. In other words, given a set of observed massive black hole coalescences, we assess what information can be extracted about the underlying massive black hole population model. For concreteness we consider ten specific models of massive black hole formation, chosen to probe four important (and largely unconstrained) aspects of the input physics used in structure formation simulations: seed formation, metallicity ``feedback'', accretion efficiency and accretion geometry. For the first time we allow for the possibility of ``model mixing'', by drawing the observed population from some combination of the ``pure'' models that have been simulated. A Bayesian analysis allows us to recover a posterior probability distribution for the ``mixing parameters'' that characterize the fractions of each model represented in the observed distribution. Our work shows that LISA has enormous potential to probe the underlying physics of structure formation.

We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries in the context of the LISA mission. Detection of multiple events would provide invaluable information on competing theories of gravity, evolution and formation of structures and, with complementary observations, constraints on H_0 and other cosmological parameters. Most of the optical depth for lensing is provided by intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple events are sizable for a broad range of formation histories. For the most optimistic models, up to 4 multiple events with a signal to noise ratio >= 8 are expected in a 5-year mission. Chances are significant even for conservative models with either light (<= 60%) or heavy (<= 40%) seeds. Due to lensing amplification, some intrinsically too faint signals are brought over threshold (<= 2 per year).

In this proceeding we explore a pathway to radio-loudness under the hypothesis that retrograde accretion onto giant spinning black holes leads to the launch of powerful jets, as seen in radio loud QSOs and recently in LAT/Fermi and BAT/Swift Blazars. Counter-rotation of the accretion disc relative to the BH spin is here associated to gas-poor galaxy mergers progenitors of giant (missing-light) ellipticals. The occurrence of retrograde accretion enters as unifying element that may account for the radio-loudness/galaxy morphology dichotomy observed in AGN.

Black hole spins affect the efficiency of the "classical" accretion processes, hence the radiative output from quasars. Spins also determine how much energy is extractable from the hole itself. Recently it became clear that massive black hole spins also affect the retention of black holes in galaxies, be cause of the impulsive "gravitational recoil", up to thousands km/s, due to anisotropic emission of gravitational waves at merger. I discuss here the evolution of massive black hole spins along the cosmic history, due to the combination of mergers and accretion events. I describe recent simulations of accreting black holes in merger remnants, and discuss the implication for the spins of black holes in quasars.

With the release of the first year Fermi catalogue, the number of blazars detected above 100 MeV lying at high redshift has been largely increased. There are 28 blazars at z>2 in the "clean" sample. All of them are Flat Spectrum Radio Quasars (FSRQs). We study and model their overall spectral energy distribution in order to find the physical parameters of the jet emitting region, and for all of them we estimate their black hole masses and accretion rates. We then compare the jet with the accretion disk properties, setting these sources in the broader context of all the other bright gamma-ray or hard X-ray blazars. We confirm that the jet power correlates with the accretion luminosity. We find that the high energy emission peak shifts to smaller frequencies as the observed luminosity increases, according to the blazar sequence, making the hard X-ray band the most suitable for searching the most luminous and distant blazars.

We explore the cosmic evolution of massive black hole (MBH) seeds forming within 'quasistars' (QSs), accreting black holes embedded within massive hydrostatic gaseous envelopes. These structures could form if the infall of gas into the center of a halo exceeds about 1 solar mass per year. We use a merger-tree approach to estimate the rate at which QSs might form as a function of redshift, and the statistical properties of the resulting QS and seed black hole populations. We relate the triggering of runaway infall to major mergers of gas-rich galaxies, and to a threshold for global gravitational instability, which we link to the angular momentum of the host. This is the main parameter of our models. Once infall is triggered, its rate is determined by the halo potential; the properties of the resulting QS and seed black hole depend on this rate. After the epoch of QSs, we model the growth of MBHs within their hosts in a merger-driven accretion scenario. We compare MBH seeds grown inside quasistars to a seed model that derives from the remnants of the first metal-free stars, and also study the case in which both channels of MBH formation operate simultaneously. We find that a limited range of QS/MBH formation efficiencies exists that allows one to reproduce observational constraints. Our models match the density of z = 6 quasars, the cumulative mass density accreted onto MBHs (according to Soltan's argument), and the current mass density of MBHs. The mass function of QSs peaks at mass ~ 1e6 solar masses, and we calculate the number counts for the JWST in the 2-10 micron band. We find that JWST could detect up to several QSs per field at z ~ 5 - 10.

As massive black holes (MBHs) grow from lower-mass seeds, it is natural to expect that a leftover population of progenitor MBHs should also exist in the present universe. Dwarf galaxies undergo a quiet merger history, and as a result, we expect that dwarfs observed in the local Universe retain some `memory' of the original seed mass distribution. Consequently, the properties of MBHs in nearby dwarf galaxies may provide clean indicators of the efficiency of MBH formation. In order to examine the properties of MBHs in dwarf galaxies, we evolve different MBH populations within a Milky Way halo from high-redshift to today. We consider two plausible MBH formation mechanisms: `massive seeds' formed via gas-dynamical instabilities and a Population III remnant seed model. `Massive seeds' have larger masses than PopIII remnants, but form in rarer hosts. We dynamically evolve all halos merging with the central system, taking into consideration how the interaction modifies the satellites, stripping their outer mass layers. We compute different properties of the MBH population hosted in these satellites. We find that for the most part MBHs retain the original mass, thus providing a clear indication of what the properties of the seeds were. We derive the black hole occupation fraction (BHOF) of the satellite population at z=0. MBHs generated as `massive seeds' have large masses that would favour their identification, but their typical BHOF is always below 40 per cent and decreases to less than per cent for observed dwarf galaxy sizes. In contrast, Population III remnants have a higher BHOF, but their masses have not grown much since formation, inhibiting their detection.

One of the sources of gravitational waves for the proposed space-based gravitational wave detector, the Laser Interferometer Space Antenna (LISA), are the inspirals of compact objects into supermassive black holes in the centres of galaxies - extreme-mass-ratio inspirals (EMRIs). Using LISA observations, we will be able to measure the parameters of each EMRI system detected to very high precision. However, the statistics of the set of EMRI events observed by LISA will be more important in constraining astrophysical models than extremely precise measurements for individual systems. The black holes to which LISA is most sensitive are in a mass range that is difficult to probe using other techniques, so LISA provides an almost unique window onto these objects. In this paper we explore, using Bayesian techniques, the constraints that LISA EMRI observations can place on the mass function of black holes at low redshift. We describe a general framework for approaching inference of this type --- using multiple observations in combination to constrain a parameterised source population. Assuming that the scaling of EMRI rate with black hole mass is known and taking a black hole distribution given by a simple power law, dn/d(ln M) = A (M/M_*)^b, we find that LISA could measure the parameters to a precision of D(ln A) ~ 0.08, and D(b) ~ 0.03 for a reference model that predicts ~1000 events. Even with as few as 10 events, LISA should constrain the slope to a precision ~0.3, which is the current level of observational uncertainty in the low-mass slope of the black hole mass function. We also consider a model in which A and b evolve with redshift, but find that EMRI observations alone do not have much power to probe such an evolution.

Evidence shows that massive black holes reside in most local galaxies. Studies have also established a number of relations between the MBH mass and properties of the host galaxy such as bulge mass and velocity dispersion. These results suggest that central MBHs, while much less massive than the host (~ 0.1%), are linked to the evolution of galactic structure. In hierarchical cosmologies, a single big galaxy today can be traced back to the stage when it was split up in hundreds of smaller components. Did MBH seeds form with the same efficiency in small proto-galaxies, or did their formation had to await the buildup of substantial galaxies with deeper potential wells? I briefly review here some of the physical processes that are conducive to the evolution of the massive black hole population. I will discuss black hole formation processes for `seed' black holes that are likely to place at early cosmic epochs, and possible observational tests of these scenarios.

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the ~ 1-10 Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, ~ 100-1000 solar masses. We primarily discuss two different source types -- mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-10 solar masses into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.

We investigate the physical properties of the 10 blazars at redshift greater than 2 detected in the 3-years all sky survey performed by the Burst Alert Telescope (BAT) onboard the Swift satellite. We find that the jets of these blazars are among the most powerful known. Furthermore, the mass of their central black hole, inferred from the optical-UV bump, exceeds a few billions of solar masses, with accretion luminosities being a large fraction of the Eddington one. We compare their properties with those of the brightest blazars of the 3-months survey performed by the Large Area Telescope (LAT) onboard the Fermi satellite. We find that the BAT blazars have more powerful jets, more luminous accretion disks and larger black hole masses than LAT blazars. These findings can be simply understood on the basis of the blazar sequence, that suggests that the most powerful blazars have a spectral energy distribution with a high energy peak at MeV (or even sub-MeV) energies. This implies that the most extreme blazars can be found more efficiently in hard X-rays, rather than in the high energy gamma-ray band. We then discuss the implications of our findings for future missions, such as the New Hard X-ray Mission (NHXM) and especially the Energetic X-ray Imaging Survey Telescope (EXIST) mission which, during its planned 2 years all sky survey, is expected to detect thousands of blazars, with a few of them at z greater than 6.

We assess the influence of massive black hole (MBH) ejections from galaxy centres, due to the gravitational radiation recoil, along the cosmic merger history of the MBH population. We discuss the 'danger' of the recoil for MBHs as a function of different MBH spin/orbit configurations and of the host halo cosmic bias, and on how that reflects on the 'occupation fraction' of MBHs. We assess ejection probabilities for mergers occurring in a gas-poor environment, where the MBH binary coalescence is driven by stellar dynamical processes, and the spin/orbit configuration is expected to be isotropically distributed. We contrast this case with the 'aligned' case. The latter is the most realistic situation for 'wet', gas-rich mergers, which are the expectation for high-redshift galaxies. We find that if all halos at z>5-7 host a MBH, the probability of the Milky Way (or similar size galaxy) to host a MBH today is less than 50%, unless MBHs form continuously in galaxies. The 'occupation fraction' of MBHs, intimately related to halo bias and MBH formation efficiency, plays a crucial role in increasing the retention fraction. Small halos, with shallow potential wells and low escape velocities, have a high ejection probability, but the MBH merger rate is very low along their galaxy formation merger hierarchy: MBH formation processes are likely inefficient in such shallow potential wells. Recoils can decrease the overall frequency of MBHs in small galaxies to ~60%, while they have little effect on the frequency of MBHs in large galaxies (at most a 20% effect).

Galactic nuclei host central massive objects either in the form of supermassive black holes or nuclear stellar clusters. Recent investigations have shown that both components co-exist in at least a few galaxies. In this paper we explore the possibility of a connection between nuclear star clusters and black holes that establishes at the moment of their formation. We here model the evolution of high redshift discs, hosted in dark matter halos with virial temperatures 10^4 K, whose gas has been polluted with metals just above the critical metallicity for fragmentation. A nuclear cluster forms as a result of a central starburst from gas inflowing from the unstable disc. The nuclear stellar cluster provides a suitable environment for the formation of a black hole seed, ensuing from runaway collisions among the most massive stars. Typical masses for the nuclear stellar clusters at the time of black hole formation (z~10) are inthe range 10^4-10^6 solar masses and have half mass radii < 0.5 pc. The black holes forming in these dense, high redshift clusters can have masses in the range ~300-2000 solar masses.

Massive black holes are nowadays believed to reside in most local galaxies. Studies have also established a number of relations between the MBH mass and properties of the host galaxy such as bulge mass and velocity dispersion. These results suggest that central MBHs, while much less massive than the host (~0.1%), are linked to the evolution of galactic structure. When did it all start? In hierarchical cosmologies, a single big galaxy today can be traced back to the stage when it was split up in hundreds of smaller components. Did MBH seeds form with the same efficiency in small proto-galaxies, or did their formation had to await the buildup of substantial galaxies with deeper potential wells? I briefly review here some of the physical processes that are conducive to the evolution of the massive black hole population. I will discuss black hole formation processes for `seed' black holes that are likely to place at early cosmic epochs, and possible observational tests of these scenarios.

Using high resolution hydrodynamical simulations, we explore the spin evolution of massive dual black holes orbiting inside a circumnuclear disc, relic of a gas-rich galaxy merger. The black holes spiral inwards from initially eccentric co or counter-rotating coplanar orbits relative to the disc's rotation, and accrete gas that is carrying a net angular momentum. As the black hole mass grows, its spin changes in strength and direction due to its gravito-magnetic coupling with the small-scale accretion disc. We find that the black hole spins loose memory of their initial orientation, as accretion torques suffice to align the spins with the angular momentum of their orbit on a short timescale (<1-2 Myr). A residual off-set in the spin direction relative to the orbital angular momentum remains, at the level of <10 degrees for the case of a cold disc, and <30 degrees for a warmer disc. Alignment in a cooler disc is more effective due to the higher coherence of the accretion flow near each black hole that reflects the large-scale coherence of the disc's rotation. If the massive black holes coalesce preserving the spin directions set after formation of a Keplerian binary, the relic black hole resulting from their coalescence receives a relatively small gravitational recoil. The distribution of recoil velocities inferred from a simulated sample of massive black hole binaries has median <70 km/s much smaller than the median resulting from an isotropic distribution of spins.

In this paper, we explore the establishment and evolution of the empirical correlation between black hole mass and velocity dispersion with redshift. We track the growth and accretion history of massive black holes starting from high redshift using two seeding models:(i) Population III remnants, and (ii) massive seeds from direct gas collapse. Although the seeds do not initially satisfy the $M_{\rm BH} - \sigma$ relation, the correlation is established and maintained at all times if self-regulating accretion episodes are associated with major mergers. The massive end of the $M_{\rm BH} - \sigma$ relation is established early, and lower mass MBHs migrate over time. How MBHs migrate toward the relation, the slope and the scatter of the relation all depend critically on the seeding model as well as the adopted self-regulation prescription. We expect flux limited AGN surveys and LISA to select accreting and merging MBHs respectively that have already migrated onto the $\msigma$ relation. This is a consequence of major mergers being more common at high redshift for the most massive, biased, galaxies that anchor the $\msigma$ relation early. We also predict the existence of a large population of low mass `hidden' MBHs at high redshift which can easily escape detection. Additionally, we find that if MBH seeds are massive, $\sim 10^5 M_{\odot}$, the low-mass end of the $\msigma$ flattens towards this asymptotic value, creating a characteristic `plume'.

We investigate the theoretical expectations for detections of supermassive binary black holes that can be identified as sub-parsec luminous quasars. To-date, only two candidates have been selected in a sample comprising 17,500 sources selected from the Sloan Digital Sky Survey (SDSS) Quasar Catalog at z<0.70 (Boroson & Lauer 2009) In this Letter, we use models of assembly and growth of supermassive black holes (SMBHs) in hierarchical cosmologies to study the statistics and observability of binary quasars at sub-parsec separations. Our goal is twofold: (1) test if such a scarce number of binaries is consistent with theoretical prediction of SMBH merger rates, and (2) provide additional predictions at higher redshifts, and at lower flux levels. We determine the cumulative number of expected binaries in a complete, volume limited sample. Motivated by Boroson & Lauer (2009), we apply the SDSS Quasar luminosity cut (M_i<-22) to our theoretical sample, deriving an upper limit to the observable binary fraction. We find that sub-parsec quasar binaries are intrinsically rare. Our best models predict ~0.01 deg^-2 sub-parsec binary quasars with separations below ~10^4 Schwarzschild radii (v_orb>2000 km/s) at z<0.7, which represent a fraction ~6x10^-4 of unabsorbed quasars in our theoretical sample. In a complete sample of ~10,000 sources, we therefore predict an upper limit of ~10 sub-parsec binary quasars. The number of binaries increases rapidly with increasing redshift. The decreasing lifetime with SMBH binary mass suggests that lowering the luminosity threshold does not lead to a significant increase in the number of detectable sub-parsec binary quasars.

We revisit the classical view that quasar contribution to the reionization of hydrogen is unimportant. When secondary ionization are taken into account, in many plausible scenarios for the formation and growth of supermassive black holes quasars contribute substantially or even dominantly at z>8, although their contribution generally falls below that of star-forming galaxies by z=6. Theoretical models that guide the design of the first generation of redshifted 21 cm experiments must, therefore, substantially account for the quasar contribution in order to be even qualitatively accurate.

We report on the results of spectral fits made to data obtained from a 168 ksec Suzaku observation of the Seyfert-1 galaxy Fairall 9. The source is clearly detected out to 30 keV. The observed spectrum is fairly simple; it is well-described by a power-law with a soft excess and disk reflection. A broad iron line is detected, and easily separated from distinct narrow components owing to the resolution of the CCDs in the X-ray Imaging Spectrometer (XIS). The broad line is revealed to be asymmetric, consistent with a disk origin. We fit the XIS and Hard X-ray Detector (HXD) spectra with relativistically-blurred disk reflection models. With the assumption that the inner disk extends to the innermost stable circular orbit, the best-fit model implies a black hole spin parameter of a = 0.60(7) and excludes extremal values at a high level of confidence. We discuss this result in the context of Seyfert observations and models of the cosmic distribution of black hole spin.

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass and spin parameter evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction. We develop an iterative scheme based on the adiabatic approximation to study the BH-disc coevolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen-Petterson effect) and interacts with the BH until the BH spin aligns with the outer angular momentum direction. For a BH aligning with a co-rotating disc, the fractional increase in mass is typically less than a few percent, while the spin modulus can increase up to a few tens of percent. The alignment timescale is between ~ 100 thousands and ~ 1 millions years for a maximally rotating BH accreting at the Eddington rate. BH-disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular co-rotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.

A strong hard X-ray luminosity from a blazar flags the presence of a very powerful jet. If the jet power is in turn related to the mass accretion rate, the most luminous hard X-ray blazars should pinpoint the largest accretion rates, and therefore the largest black hole masses. These ideas are confirmed by the Swift satellite observations of the blazar S5 0014+813, at the redshift z=3.366. Swift detected this source with all its three instruments, from the optical to the hard X-rays. Through the construction of its spectral energy distribution we are confident that its optical-UV emission is thermal in origin. Associating it to the emission of a standard optically thick geometrically thin accretion disk, we find a black hole mass of 40 billion solar masses, radiating at 40% the Eddington value. The derived mass is among the largest ever found. Super-Eddington slim disks or thick disks with the presence of a collimating funnel can in principle reduce the black hole mass estimate, but tends to produce spectra bluer than observed.

If the cosmological dark matter has a component made of small primordial black holes, they may have a significant impact on the physics of the first stars and on the subsequent formation of massive black holes. Primordial black holes would be adiabatically contracted into these stars and then would sink to the stellar center by dynamical friction, creating a larger black hole which may quickly swallow the whole star. If these primordial black holes are heavier than $\sim 10^{22} {\rm g}$, the first stars would likely live only for a very short time and would not contribute much to the reionization of the universe. They would instead become $10 - 10^3 M_\odot$ black holes which (depending on subsequent accretion) could serve as seeds for the super--massive black holes seen at high redshifts as well as those inside galaxies today.

In this Letter we explore the hypothesis that the quasar SDSSJ092712.65+294344.0 is hosting a massive black hole binary embedded in a circumbinary disc. The lightest, secondary black hole is active, and gas orbiting around it is responsible for the blue-shifted broad emission lines with velocity off-set of 2650 km/s, relative to the galaxy rest frame. As the tidal interaction of the binary with the outer disc is expected to excavate a gap, the blue-shifted narrow emission lines are consistent with being emitted from the low-density inhomogeneous gas of the hollow region. From the observations we infer a binary mass ratio q ~ 0.3, a mass for the primary of M1 ~ 2 billion Msun and a semi-major axis of 0.34 pc, corresponding to an orbital period of 370 years. We use the results of cosmological merger trees to estimate the likely-hood of observing SDSSJ092712.65+294344.0 as recoiling black hole or as a binary. We find that the binary hypothesis is preferred being one hundred times more probable than the ejection hypothesis. If SDSSJ092712.65+294344.0 hosts a binary, it would be the one closest massive black hole binary system ever discovered.

We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices implying obscuring column densities of N_H > 10^23 cm^-2. The AGN fraction should be considered a lower limit, since several more LABs not detected with Chandra show AGN signatures in their mid-infrared emission. We show that the UV luminosities of the AGN are easily capable of powering the extended Ly-alpha emission via photo-ionization alone. When combined with the UV flux from a starburst component, and energy deposited by mechanical feedback, we demonstrate that 'heating' by a central source, rather than gravitational cooling is the most likely power source of LABs. We argue that all LABs could be powered in this manner, but that the luminous host galaxies are often just below the sensitivity limits of current instrumentation, or are heavily obscured. No individual LABs show evidence for extended X-ray emission, and a stack equivalent to a >9Ms exposure of an average LAB also yields no statistical detection of a diffuse X-ray component. The resulting diffuse X-ray/Ly-alpha luminosity limit implies there is no hot (T>10^7 K) gas component in these halos, and also rules out inverse Compton scattering of cosmic microwave background photons, or local far-infrared photons, as a viable power source for LABs.

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or co-rotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an "orbital angular momentum flip" due to the gas-dynamical friction, i.e., it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi-Hoyle-Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly "lights up" at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ~10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is co-rotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double AGNs at separations of <100 pc.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line.

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

Massive black holes are key components of the assembly and evolution of cosmic structures and a number of surveys are currently on-going or planned to probe the demographics of these objects and to gain insight into the relevant physical processes. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses >10^7 solar masses. The whole cosmic population produces a stochastic background that could be detectable with upcoming Pulsar Timing Arrays. Sources sufficiently close and/or massive generate gravitational radiation that significantly exceeds the level of the background and could be individually resolved. We consider a wide range of massive black hole binary assembly scenarios, we investigate the distribution of the main physical parameters of the sources, such as masses and redshift, and explore the consequences for Pulsar Timing Arrays observations. Depending on the specific massive black hole population model, we estimate that on average at least one resolvable source produces timing residuals in the range ~5-50 ns. Pulsar Timing Arrays, and in particular the future Square Kilometre Array (SKA), can plausibly detect these unique systems, although the events are likely to be rare. These observations would naturally complement on the high-mass end of the massive black hole distribution function future surveys carried out by the Laser Interferometer Space Antenna (LISA)

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

Since the seminal work of Penrose (1969) and Blandford & Znajek (1977), it has been realized that black hole spin may be an important energy source in astrophysics. The radio-loud/radio-quiet dichotomy in the AGN population is usually attributed to differences in black hole spin, with correlations between black hole spin and host galaxy morphology being hypothesized in order to explain why radio-loud AGN occur in early-type galaxies. X-ray observations are uniquely able to answer: Does black hole spin play a crucial role in powering relativistic jets such as those seen from radio-loud active galactic nuclei (AGN), Galactic microquasars, and Gamma-Ray Bursts? Indeed, the importance of black hole spin goes beyond its role as a possible power source: the spin of a supermassive black hole is a fossil record of its formation and subsequent growth history.

We present a model for the formation of massive black holes ($\sim 1000 \msun$) due to stellar-dynamical processes in the first stellar clusters formed at early cosmic times ($z\sim10-20$). The high redshift black hole seeds form as a result of multiple successive instabilities that occur in low metallicity $Z\sim 10^{-5}Z_\odot$) protogalaxies. We focus on relatively massive halos at high redshift ($T_{\rm vir} > 10^4$ K, $z\gsim 10$) after the very first stars in the Universe have completed their evolution. This set of assumptions ensures that (i) atomic hydrogen cooling can contribute to the gas cooling process, (ii) a UV field has been created by the first stars, and (iii) the gas inside the halo has been mildly polluted by the first metals. The second condition implies that at low density $H_2$ is dissociated and does not contribute to cooling. The third condition sets a minimum threshold density for fragmentation, so that stars form efficiently only in the very inner core of the protogalaxy. Within this core, very compact stellar clusters form. A large fraction of these very dense clusters undergo core collapse before stars are able to complete stellar evolution. Runaway star-star collisions eventually lead to the formation of a very massive star, leaving behind a massive black hole remnant. Clusters unstable to runaway collisions are always the first, less massive ones that form. Typically a fraction $\sim 0.05$ of protogalaxies at $z\sim10-20$ form black hole seeds, with masses $\sim 1000-2000 \msun$, leading to a mass density in seeds of a few $\simeq10^2\msun/{\rm Mpc}^{-3}$. This density allows enough room for black hole growth by accretion during the quasar epoch.

Anisotropic gravitational radiation from a coalescing black hole binary is known to impart recoil velocities of up to ~1000 km/s to the remnant black hole. In this context, we study the motion of a recoiling black hole inside a galaxy modelled as an Hernquist sphere, and the signature that the hole imprints on the hot gas, using N-body/SPH simulations. Ejection of the black hole results in a sudden expansion of the gas ending with the formation of a gaseous core, similarly to what is seen for the stars. A cometary tail of particles bound to the black hole is initially released along its trail. As the black hole moves on a return orbit, a nearly spherical swarm of hot gaseous particles forms at every apocentre: this feature can live up to ~ 100 Myr. If the recoil velocity exceeds the sound speed initially, the black hole shocks the gas in the form of a Mach cone in density near each super-sonic pericentric passage. We find that the X-ray fingerprint of a recoiling black hole can be detected in Chandra X-ray maps out to a distance of Virgo. For exceptionally massive black holes the Mach cone and the wakes could be observed out to a few hundred of Mpc. Detection of the Mach cone is found to become of twofold importance: i) as a probe of high-velocity recoils and ii) as an assessment of the scatter of the mass-sigma relation at large black hole masses.

We analyze the available sample of double quasars, and investigate their physical properties. Our sample comprises 85 pairs, selected from the Sloan Digital Sky Survey (SDSS). We derive physical parameters for the engine and the host, and model the dynamical evolution of the pair. First, we compare different scaling relationships between massive black holes and their hosts (bulge mass, velocity dispersion, and their possible redshift dependences), and discuss their consistency. We then compute dynamical friction timescales for the double quasar systems to investigate their frequency and their agreement with scenarios for quasar triggering. Comparing typical merging timescales to expected quasar lifetimes, the fraction of double quasars should be roughly a factor of 10 larger than observed. Additionally, we find that, depending on the correlations between black holes and their hosts, the occurrence of double quasars could be redshift-dependent. Comparison of our models to the SDSS quasar catalog suggests that double quasars should be more common at high redshift. We compare the typical separations at which double quasars are observed to the predictions of merger simulations. We find that the distribution of physical separations peaks at ~30 kpc, with a tail at larger separations (~100-200 kpc). These large separation pairs are instead consistent with unequal mass mergers where gas is dynamically perturbed during the first pericentric passage, but the gas reaches the black hole only at the next apocenter, where the pair is observed.

All the physical processes involved in the formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures are likely to leave an imprint on the gravitational waves detectable by future space--borne missions, such as LISA. We report here the results of recent studies, carried out by means of dedicated simulations of black hole build--up, aiming at understanding the impact on LISA observations of two ingredients that are crucial in every massive black hole formation scenario, namely: (i) the nature and abundance of the first black hole seeds and (ii) the large gravitational recoils following the merger of highly spinning black holes. We predict LISA detection rates spanning two order of magnitude, in the range 3-300 events per year, depending on the detail of the assumed massive black hole seed model. On the other hand, large recoil velocities do not dramatically compromise the efficiency of LISA observations. The number of detections may drop substantially (by ~60%), in scenarios characterized byabundant light seeds, but if seeds are already massive and/or relatively rare, the detection rate is basically unaffected.

We present results from a new ultra-deep 400 ks Chandra observation of the SSA22 protocluster at z = 3.09. We have studied the X-ray properties of 234 z ~ 3 Lyman break galaxies (LBGs; protocluster and field) and 158 z = 3.09 Ly-alpha emitters (LAEs) in SSA22 to measure the influence of the high-density protocluster environment on the accretion activity of supermassive black holes (SMBHs) in these UV-selected star forming populations. We detect individually X-ray emission from active galactic nuclei (AGNs) in six LBGs and five LAEs; due to small overlap between the LBG and LAE source population, ten of these sources are unique. At least six and potentially eight of these sources are members of the protocluster. These sources have rest-frame 8-32 keV luminosities in the range of L_8-32 keV = (3-50) X 10^{43} ergs/s and an average observed-frame 2-8 keV to 0.5-2 keV band-ratio of ~0.8 (mean effective photon index of Gamma_eff = 1.1), suggesting significant absorption columns of N_H > 10^{22}-10^{24} cm^{-2}. We find that the fraction of LBGs and LAEs in the z = 3.09 protocluster harboring an AGN with L_8-32 keV > 3 X 10^{43} ergs/s is 9.5^{+12.7}_{-6.1}% and 5.1^{+6.8}_{-3.3}%, respectively. These AGN fractions are somewhat larger (by a mean factor of 6.1^{+10.3}_{-3.6}; significant at the 95% confidence level) than z ~ 3 sources found in lower-density "field" environments. Theoretical models imply that these results may be due to the presence of more actively growing and/or massive SMBHs in LBGs and LAEs within the protocluster compared to the field. Such a result is expected in a scenario where enhanced merger activity in the protocluster drives accelerated galaxy and SMBH growth at z > 2-3. (abridged)

During the final phases of inspiral, a massive black hole (MBH) binary experiences a recoil due to the asymmetric emission of gravitational waves. We use recent results from numerical relativity simulations together with models of the assembly and growth of MBHs in hierarchical cosmologies, to study the dynamics, statistics, and observability of recoling MBHs. We find that, at redshift z<3, kicked non-rotating holes are typically found between 1 and 30 kpc from their galaxy centers, while rapidly rotating ones are typically between 10 and a few hundred kpc. A recoiling hole that carries an accretion disk may shine as an "off-nuclear AGN" while it moves away from the center of its host galaxy. We predict that, depending on the hole spin distribution and the duration of their active phase, a population of off-nuclear AGN may already be detectable at low and intermediate redshifts in present deep Hubble Space Telescope observations. The James Webb Space Telescope may discover tens of wandering AGN per square degree, most of them moving within their host halos on unbound trajectories.

Using recent results from numerical relativity simulations of black hole mergers, we revisit previous studies of cosmological black hole spin evolution. We show that mergers are very unlikely to yield large spins, unless alignment of the spins of the merging holes with the orbital angular momentum is very efficient. We analyze the spin evolution in three specific scenarios: (1) spin evolves only through mergers, (2) spin evolves through mergers and prolonged accretion episodes, (3) spin evolves through mergers and short-lived (chaotic) accretion episodes. We study how different diagnostics can distinguish between these evolutionary scenarios, assessing the discriminating power of gravitational-wave measurements and X-ray spectroscopy. Gravitational radiation can produce three different types of spin measurements, yielding respectively the spins of the two black holes in a binary inspiral prior to merger, the spin of the merger remnant (as encoded in the ringdown waves), and the spin of ``single'' black holes during the extreme mass-ratio inspiral (EMRI) of compact objects. The latter spin population is also accessible to iron-line measurements. We compute and compare the spin distributions relevant for these different observations. If iron-line measurements and gravitational-wave observations of EMRIs only yield dimensionless spins j=J/M^2>0.9, then prolonged accretion should be responsible for spin-up, and chaotic accretion scenarios would be very unlikely. If only a fraction of the whole population of low-redshift black holes spins rapidly, spin-alignment during binary mergers (rather than prolonged accretion) could be responsible for spin-ups.

The relationship between galaxies and supermassive black holes (SMBH) found in their cores plays a key role in the formation and evolution of both of these major constituents of the universe, as well as the evolution of the intergalactic medium. Neither can be fully understood on their own, and studies of galaxy and SMBH co-formation and co-evolution are now among the central topics of research in cosmology. Yet the very origins, and the early growth phases of the SMBH are still not firmly established. We review our current understanding of the relevant processes and their astrophysical and cosmological context, with an emphasis on the observability of the SMBH growth mechanisms at high redshifts, and their leftover progeny at low redshifts.

We investigate the distribution of massive black holes (MBHs) in the Virgo cluster. Observations suggest that AGN activity is widespread in massive galaxies (M>1e10 solar masses), while at lower galaxy masses star clusters are more abundant, which might imply a limited presence of central black holes in these galaxy-mass regimes. We explore if this possible threshold in MBH hosting, is linked to nature, nurture, or a mixture of both. The nature scenario arises naturally in hierarchical cosmologies, as MBH formation mechanisms typically are efficient in biased systems, which would later evolve into massive galaxies. Nurture, in the guise of MBH ejections following MBH mergers, provides an additional mechanism that is more effective for low mass, satellite galaxies. The combination of inefficient formation, and lower retention of MBHs, leads to the natural explanation of the distribution of compact massive ob jects in Virgo galaxies. If MBHs arrive to the correlation with the host mass and velocity dispersion during merger-triggered accretion episodes, sustained tidal stripping of the host galaxies creates a population of MBHs which lie above the expected scaling between the holes and their host mass, suggesting a possible environmental dependence.

We investigate the evolution of high redshift seed black hole masses at late times and their observational signatures. The massive black hole seeds studied here form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this way, we follow the mass assembly of these black holes to the present time using a Monte-Carlo merger tree. Using this machinery we predict the black hole mass function at high redshifts and at the present time; the integrated mass density of black holes and the luminosity function of accreting black holes as a function of redshift. These predictions are made for a set of three seed models with varying black hole formation efficiency. Given the accuracy of current observational constraints, all 3 models can be adequately fit. Discrimination between the models appears predominantly at the low mass end of the present day black hole mass function which is not observationally well constrained. However, all our models predict that low surface brightness, bulgeless galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at z=0. This effect is more pronounced for low mass galaxies. This is the key discriminant between the models studied here and the Population III remnant seed model. We find that there exists a population of low mass galaxies that do not host nuclear black holes. Our prediction of the shape of the black hole mass - velocity dispersion relation at the low mass end is in agreement with the recent observational determination from the census of low mass galaxies in the Virgo cluster.

Supermassive black holes are nowadays believed to reside in most local galaxies. Accretion of gas and black hole mergers play a fundamental role in determining the two parameters defining a black hole: mass and spin. I briefly review here some of the physical processes that are conducive to the evolution of the massive black hole population. I'll discuss black hole formation processes that are likely to place at early cosmic epochs, and how massive black hole evolve in a hierarchical Universe. The mass of the black holes that we detect today in nearby galaxy has mostly been accumulated by accretion of gas. While black hole--black hole mergers do not contribute substantially to the final mass of massive black holes, they influence the occupancy of galaxy centers by black hole, owing to the chance of merging black holes being kicked from their dwellings due to the gravitational recoil. Similarly, accretion leaves a deeper imprint on the distribution of black hole spins than black hole mergers do. The differences in accretion histories for black holes hosted in elliptical or disc galaxies may reflect on different spin distributions.

We investigate the conjecture by Sikora, Stawarz & Lasota (2007) that the observed AGN--radio-loudness bimodality can be explained by the morphology-related bimodality of black-hole spin distribution in the centers of galaxies: central black holes in giant elliptical galaxies may have (on average) much larger spins than black holes in spiral/disc galaxies. We study how accretion from a warped disc influences the evolution of black hole spins and conclude that within the cosmological framework, where the most massive BHs have grown in mass via merger driven accretion, one indeed expects most supermassive black holes in elliptical galaxies to have on average higher spin than black holes in spiral galaxies, where random, small accretion episodes (e.g. tidally disrupted stars, accretion of molecular clouds) might have played a more important role.

In the last stages of a black hole merger, the binary can experience a recoil due to asymmetric emission of gravitational radiation. Recent numerical relativity simulations suggest that the recoil velocity can be as high as a few thousands kilometers per second for particular configurations. We consider here the effect of a worst case scenario for orbital and phase configurations on the hierarchical evolution of the massive black hole (MBH) population. The orbital configuration and spin orientation in the plane is chosen to be the one yielding the highest possible kick. Masses and spin magnitudes are instead derived self-consistently from the MBH evolutionary models. If seeds form early, e.g. as remnants of the first stars, almost the totality of the first few generation of binaries are ejected. The fraction of lost binaries decreases at later times due to a combination of the binary mass ratio distribution becoming shallower, and the deepening of the hosts potential wells. If seeds form at later times, in more massive halos, the retention rate is much higher. We show that the gravitational recoil does not pose a threat to the evolution of the MBH population that we observe locally in either case, although high mass seeds seem to be favored. The gravitational recoil is instead a real hazard for (i) MBHs in biased halos at high-redshift, where mergers are more common, and the potential wells still relatively shallow. Similarly, it is very challenging to retain (ii) MBHs merging in star clusters.

The formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures leaves a unique imprint on the background of gravitational waves at mHz frequencies. We study here, by means of dedicated simulations of black hole build--up, the possibility of constraining different models of black hole cosmic evolution using future gravitational wave space--borne missions, such as LISA. We consider two main scenarios for black hole formation, namely, one where seeds are light (~10^2 \msun, remnant of Population III stars), and one where seeds are heavy (>~10^4 \msun, direct collapse). In all the models we have investigated, massive black hole binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in 3 year LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high redshift halos. On the other extreme, a low event rate, about few tens in 3 years, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries (m<10^4\msun) are predicted in a fairly large number in Population III remnants models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1-1.

We compute the probability that intermediate mass black holes (IMBHs) capture companions due to dynamical interactions and become accreting sources, and explore the possibility that the accreting IMBHs would appear as ultra-luminous X-ray sources (ULXs). We focus on IMBHs originating from low-metallicity Population III stars. Two channels of IMBH formation are considered: from primordial halos in the framework of hierarchical clustering, and from non-mixed, zero-metallicity primeval gas in galactic discs. IMBHs can form binary systems due to tidal captures of single stars and exchange interactions with existing binary systems in galactic discs. We find that neither formation mechanism of the accreting IMBH binary is able to provide enough sources to explain the observed population of ULXs. Even at sub-ULX luminosity, the total number of accreting IMBHs with $L> 10^{36}$ erg s$^{-1}$ with dynamically captured companions is found to be $< 0.01$ per galaxy.

Supermassive black holes are nowadays believed to reside in most local galaxies. Observations have revealed us vast information on the population of local and distant black holes, but the detailed physical properties of these dark massive objects are still to be proven. Accretion of gas and black hole mergers play a fundamental role in determining the two parameters defining a black hole: mass and spin. We briefly review here the basic properties of the population of supermassive black holes, focusing on the still mysterious formation of the first massive black holes, and their evolution from early times to now.

We consider simple semi-analytic models that relate the active galactic nuclei (AGN) evolution to the merging history of their host dark matter haloes and quantify their ability of matching the AGN luminosity function and its spatial clustering at low and intermediate redshifts. In particular, we focus on the recent determinations of the AGN luminosity function in the hard X-ray band at z~0 which constitutes the most stringent observational test for our models. Indeed, while we find an acceptable agreement between the model bolometric luminosity function and the data at 1<z<2 and for luminosities larger than 1e10 solar bolometric luminosities, no semi-analytic model is capable of reproducing the number density of faint X-ray sources in the local universe. Some improvement can be obtained by advocating energy feedback that we model through a time-dependent Eddington ratio. Even in this case, however, the number density of faint AGNs is significantly below observations. This failure indicates that major mergers cannot constitute the only trigger to accretion episodes in the local AGN population.

With the aid of a simple yet robust approach we investigate the influence of dissipationless and dissipative merging on galaxy structure, and the consequent effects on the scaling laws followed by elliptical galaxies. Our results suggest that ellipticals cannot be originated by parabolic merging of low mass spheroids only, even in presence of substantial gas dissipation. However, we also found that scaling laws such as the Faber-Jackson, Kormendy, Fundamental Plane, and the Mbh -sigma relations, when considered over the whole mass range spanned by ellipticals in the local universe, are robust against merging. We conclude that galaxy scaling laws, possibly established at high redshift by the fast collapse in pre-existing dark matter halos of gas rich and clumpy stellar distributions, are compatible with a (small) number of galaxy mergers at lower redshift.

We study, by means of dedicated simulations of massive black hole build-up, the possibility to constraint the existence and nature of the AGN population at z>6 with available and planned X-ray and near infrared space telescopes. We find that X-ray deep-field observations can set important constraints to the faint-end of the AGN luminosity function at very high redshift. Planned X-ray telescopes should be able to detect AGN hosting black holes with masses down to >10^5 Msun (i.e., X-ray luminosities in excess of 10^42 erg s^-1), and can constrain the evolution of the population of massive black hole at early times (6<z<10). We find that this population of AGN should contribute substantially (~25%) to the unresolved fraction of the cosmic X-ray background in the 0.5-10 keV range, and that a significant fraction (~3-4%) of the total background intensity would remain unaccounted even after future X-ray observations. As byproduct, we compute the expected UV background from AGN at z>6 and we discuss the possible role of AGN in the reionization of the Universe at these early epochs, showing that AGN alone can provide enough ionizing photons only in the (improbable) case of an almost completely homogeneous inter-galactic medium. Finally, we show that super-Eddington accretion, suggested by the observed QSOs at z>6, must be a very rare event, confined to black holes living in the highest density peaks.

We investigate how hierarchical models for the co-evolution of the massive black hole (MBH) and AGN population can reproduce the observed faint X-ray counts. We find that the main variable influencing the theoretical predictions is the Eddington ratio of accreting sources. We compare three different models proposed for the evolution of AGN Eddington ratio, f_Edd: constant f_Edd=1, f_Edd decreasing with redshift, and f_Edd depending on the AGN luminosity, as suggested by simulations of galactic mergers including BHs and AGN feedback. We follow the full assembly of MBHs and host halos from early times to the present in a LambdaCDM cosmology. AGN activity is triggered by halo major mergers and MBHs accrete mass until they satisfy the observed correlation with velocity dispersion. We find that all three models can reproduce fairly well the total faint X-ray counts. The redshift distribution is however poorly matched in the first two models. The Eddington ratios suggested by merger simulations predicts no turn-off of the faint end of the AGN optical luminosity function at redshifts z>=1 down to very low luminosity.

Massive black holes (MBHs), with masses in the range 10^3-10^8 Msolar, which merge with a companion black hole of similar mass are expected to be the most powerful source of gravitational radiation in the frequency range probed by LISA. MBH binaries can be detected by LISA up to z~5-15. Gravitational waves from MBH mergers can serve as a powerful tool to study the early evolution of the MBH population, and possibly the role played by MBHs in the joint evolution with their hosts. I review scenarios for the co-evolution of MBHs and cosmological structures, where MBH seeds form in pre-galactic structures. These black holes evolve then in a hierarchical fashion, following the merger hierarchy of their hosts. Accretion of gas, traced by quasar activity, plays a fundamental role in determining the two parameters defining a black hole, mass and spin. Gravitational waves, together with observations in electromagnetic bands, can help constrain the evolution of both MBH mass and spin.

The Sloan Digital Sky survey detected luminous quasars at very high redshift, z>6. Follow-up observations indicated that at least some of these quasars are powered by supermassive black holes (SMBHs) with masses in excess of billion solar masses. SMBHs, therefore, seem to have already existed when the Universe was less than 1 Gyr old, and the bulk of galaxy formation still has to take place. We investigate in this paper to which extent accretion and dynamical processes influence the early growth of SMBHs. We assess the impact of (i) black hole mergers, (ii) the influence of the merging efficiency and (iii) the negative contribution due to dynamical effects which can kick black holes out of their host halos (gravitational recoil). We find that if accretion is always limited by the Eddington rate via a thin disc, the maximum radiative efficiency allowed to reproduce the LF at z=6 is of order 12%, when the adverse effect of the gravitational recoil is taken into consideration. Dynamical effects cannot be neglected in studies of high-redshift SMBHs. If black holes can accrete at super-critical rate during an early phase, reproducing the observed SMBH mass values is not an issue, even in the case that the recoil velocity is in the upper limits range, as the mass ratios of merging binaries are skewed towards low values, where the gravitational recoil effect is very mild. We propose that SMBH growth at early times is very selective, and efficient only for black holes hosted in high density peak halos.

We describe a mechanism by which supermassive black holes can form directly in the nuclei of protogalaxies, without the need for seed black holes left over from early star formation. Self-gravitating gas in dark matter halos can lose angular momentum rapidly via runaway, global dynamical instabilities, the so-called "bars within bars" mechanism. This leads to the rapid buildup of a dense, self-gravitating core supported by gas pressure - surrounded by a radiation pressure-dominated envelope - which gradually contracts and is compressed further by subsequent infall. These conditions lead to such high temperatures in the central region that the gas cools catastrophically by thermal neutrino emission, leading to the formation and rapid growth of a central black hole. We estimate the initial mass and growth rate of the black hole for typical conditions in metal-free halos with T_vir ~ 10^4 K, which are the most likely to be susceptible to runaway infall. The initial black hole should have a mass of <~20 solar masses, but in principle could grow at a super-Eddington rate until it reaches ~ 10^4-10^6 solar masses. Rapid growth may be limited by feedback from the accretion process and/or disruption of the mass supply by star formation or halo mergers. Even if super-Eddington growth stops at \~10^3-10^4 solar masses, this process would give black holes ample time to attain quasar-size masses by a redshift of 6, and could also provide the seeds for all supermassive black holes seen in the present universe.

Supermassive black holes (SMBHs) are nowadays believed to reside in most local galaxies, and the available data show an empirical correlation between bulge luminosity - or stellar velocity dispersion - and black hole mass, suggesting a single mechanism for assembling black holes and forming spheroids in galaxy halos. The evidence is therefore in favour of a co-evolution between galaxies, black holes and quasars. In cold dark matter cosmogonies, small-mass subgalactic systems form first to merge later into larger and larger structures. In this paradigm galaxy halos experience multiple mergers during their lifetime. If every galaxy with a bulge hosts a SMBH in its center, and a local galaxy has been made up by multiple mergers, then a black hole binary is a natural evolutionary stage. The evolution of the supermassive black hole population clearly has to be investigated taking into account both the cosmological framework and the dynamical evolution of SMBHs and their hosts. The seeds of SMBHs have to be looked for in the early Universe, as very luminous quasars are detected up to redshift higher than z=6. These black holes evolve then in a hierarchical fashion, following the merger hierarchy of their host halos. Accretion of gas, traced by quasar activity, plays a fundamental role in determining the two parameters defining a black hole: mass and spin. A particularly intriguing epoch is the initial phase of SMBH growth. It is very challenging to meet the observational constraints at z=6 if BHs are not fed at very high rates in their infancy.

We investigate the ability of hierarchical models of QSO formation and evolution to match the observed luminosity, number counts and spatial clustering of quasars at redshift z<2. These models assume that the QSO emission is triggered by galaxy mergers, that the mass of the central black hole correlates with halo properties and that quasars shine at their Eddington luminosity except, perhaps, during the very early stages of evolution. We find that models based on simple analytic approximations successfully reproduce the observed B-band QSO luminosity function at all redshifts, provided that some mechanisms is advocated to quench mass accretion within haloes larger than about 1e13 Msun that host bright quasars. These models also match the observed strength of QSO clustering at z~0.8. At larger redshifts, however, they underpredict the QSO biasing which, instead, is correctly reproduced by semi-analytic models in which the halo merger history and associated BHs are followed by Monte Carlo realizations of the merger hierarchy. We show that the disagreement between the luminosity function predicted by semi-analytic models and observations can be ascribed to the use of B-band data, which are a biased tracer of the quasar population, due to obscuration.

We discuss a model for the early assembly of supermassive black holes (SMBHs) at the center of galaxies that trace their hierarchical build-up far up in the dark halo `merger tree'. Motivated by the observations of luminous quasars around redshift z=6 with SMBH masses of billion solar masses, we assess the possibility of an early phase of stable super-critical quasi-spherical accretion in the BHs hosted by metal free halos with virial temperature larger than 10000 K. We assume that the first `seed' black holes formed with intermediate masses following the collapse of the first generation of stars, in mini-halos collapsing at z=20-30 from high peaks of density fluctuations. In high redshift halos with virial temperature larger than 10000 K, conditions exist for the formation of a fat disc of gas at T_gas=5000-10000 K. Cooling via hydrogen atomic lines is in fact effective in these comparatively massive halos. The cooling and collapse of an initially spherical configuration of gas leads to a rotationally supported disc at the center of the halo if baryons preserve their specific angular momentum during collapse. The conditions for the formation of the gas disc and accretion onto a central black holes out of this supply of gas are investigated, as well as the feedback of the emission onto the host and onto the intergalactic medium. We find that even a short phase of supercritical accretion eases the requirements set by the z=6 quasars.

We compute the expected low-frequency gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies. We follow the merging history of halos and associated holes via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. MBHs get incorporated through a series of mergers into larger and larger halos, sink to the centre owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become more massive, and form a binary system. Stellar dynamical processes dominates the orbital evolution of the binary at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the system. We discuss the observability of inspiraling MBH binaries by a low-frequency gravitational wave experiment such as the planned Laser Interferometer Space Antenna (LISA), discriminating between resolvable sources and unresolved confusion noise. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~90 individual events above a S/N=5 confidence level, while expected confusion noise is well below planned LISA capabilities.

We investigate the consequences of a model of the assembly and growth of massive black holes from primordial seeds, remnants of the first generation of stars in a hierarchical structure formation scenario. Our model traces the build-up of MBHs from an early epoch, and follows the merger history of dark matter halos and their associated holes via Monte Carlo realizations of the merger hierarchy from early times to the present time. The sequence of minor and major mergers experienced by galactic halos in their hierarchical growth affects the merger history of MBHs embedded in their nuclei. So, if the formation route for the assembly of SMBHs dates back to the early universe, a large number of BH interactions is inevitable. Binary black holes coalescence timescales can be long enough for a third BH to fall in and interact with the central binary. These BH triple interactions lead typically to the final expulsion of one of the three bodies and to the recoil of the binary. Also, asymmetric emission of gravitational waves in the last stages of the black hole merging can give a recoil velocity to the centre of mass of the coalescing binary. This scenario leads to the prediction of a population of intermediate mass BHs (IMBHs) wandering in galaxy halos at the present epoch. We compute the luminosity distribution produced by these IMBHs accreting from their circumstellar medium. We find that in a Milky Way-sized galaxy they are unable to account for sources with luminosities > 1e39 erg/s unless they carry a baryonic remnant from which they are able to accrete for a long time. We also find that, for typical spiral galaxies, the bright end of the point source distribution correlates with the mass of the galaxy, and the most luminous sources are expected to be found in the disk.

Massive black hole binaries, with masses in the range 1E3-1E8 Msun, are expected to be the most powerful sources of gravitational radiation at mHz frequencies, and hence are among the primary targets for the planned Laser Interferometer Space Antenna (LISA). We extend and refine our previous analysis (Sesana et al. 2004), detailing the gravitational wave signal expected from a cosmological population of massive black hole binaries. As done in our previous paper, we follow the merger history of dark matter halos, the dynamics of the massive black holes they host, and their growth via gas accretion and binary coalescences in a LCDM cosmology. Stellar dynamical processes dominates the orbital evolution of black hole binaries at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the pairs. We show that the GW signal from this population, in a 3 year LISA observation, will be resolved into approx 90 discrete events with S/N>5, among which approx. 35 will be observed above threshold until coalescence. These "merging events" involve relatively massive binaries, M=10E5 Msun, in the redshift range 2<z<6. The remaining approx. 55 events come from higher redshift, less massive binaries (M=5E3 Msun at z>6) and, although their S/N integrated over the duration of the observation can be substantial, the final coalescence phase is at too high frequency to be directly observable by space-based interferometers such as LISA. LISA will be able to detect a fraction approx. 90% of all the coalescences of massive black hole binaries occurring at z=5. The residual confusion noise from unresolved massive black hole binaries is expected to be at least an order of magnitude below the estimated stochastic noise.

We study the expected distribution of massive black hole (MBH) spins and its evolution with cosmic time in the context of hierarchical galaxy formation theories. Our model uses Monte Carlo realizations of the merger hierarchy in a LCDM cosmology, coupled to semi-analytical recipes, to follow the merger history of dark matter halos, the dynamics of the MBHs they host, and their growth via gas accretion and binary coalescences. The coalescence of comparable mass holes increases the spin of MBHs, while the capture of smaller companions in randomly-oriented orbits acts to spin holes down. We find that, given the distribution of MBH binary mass ratios in hierarchical models, binary coalescences alone do not lead to a systematic spin-up or spin-down of MBHs with time: the spin distribution retains memory of its initial conditions. By contrast, because of the Bardeen-Petterson effect, gas accretion via a thin disk tends to spin holes up even if the direction of the spin axis changes randomly in time. In our models, accretion dominates over black hole captures and efficiently spins holes up. The spin distribution is heavily skewed towards fast-rotating Kerr holes, is already in place at early epochs, and does not change much below redshift 5. If accretion is via a thin disk, about 70% of all MBHs are maximally rotating and have radiative efficiencies approaching 30% (assuming a "standard'' spin-efficiency conversion). Even in the conservative case where accretion is via a geometrically thick disk, about 80% of all MBHs have spin parameters a/m > 0.8 and accretion efficiencies > 12%. Rapidly spinning holes with high radiative efficiencies may satisfy constraints based on comparing the local MBH mass density with the mass density inferred from luminous quasars (Soltan's argument).

We compute the expected gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies in a hierarchical structure formation scenario in which seed holes of intermediate mass form far up in the dark halo merger tree. The merger history of DM halos and MBHs is followed from z=20 to the present in a LCDM cosmology. MBHs get incorporated through halo mergers into larger and larger structures, sink to the center owing to dynamical friction against the DM background, accrete cold material in the merger remnant, and form MBH binary systems. Stellar dynamical interactions cause the hardening of the binary at large separations, while gravitational wave emission takes over at small radii and leads to the final coalescence of the pair. The integrated emission from inspiraling MBH binaries results in a gravitational wave background (GWB). The characteristic strain spectrum has the standard h_c(f)\propto f^{-2/3} behavior only in the range 1E-9<f<1E-6 Hz. At lower frequencies the orbital decay of MBH binaries is driven by the ejection of background stars, and h_c(f) \propto f. At higher frequencies, f>1E-6 Hz, the strain amplitude is shaped by the convolution of last stable circular orbit emission. We discuss the observability of inspiraling MBH binaries by the planned LISA. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~60 (~250) individual events above a S/N=5 (S/N=1) confidence level. (Abridged)

Motivated by the recent detection by WMAP of a large optical depth to Thomson scattering -- implying a very early reionization epoch -- we assess a scenario where the universe was reionized by `miniquasars' powered by intermediate-mass black holes (IMBHs), the remnants of the first generation of massive stars. Pregalactic IMBHs form within minihalos above the cosmological Jeans mass collapsing at z=24, get incorporated through mergers into larger and larger systems, sink to the center owing to dynamical friction, and accrete cold material. The merger history of dark halos and associated IMBHs is followed by Monte Carlo realizations of the merger hierarchy in a LCDM cosmology. While seed IMBHs that are as rare as the 3.5-sigma peaks of the primordial density field evolve largely in isolation, a significant number of black hole binary systems will form if IMBHs populate the more numerous 3-sigma peaks instead. In the case of rapid binary coalescence a fraction of IMBHs will be displaced from galaxy centers and ejected into the IGM by the `gravitational rocket' effect, rather than accrete and shine as miniquasars. We show that, under a number of plausible assumptions for the amount of gas accreted onto IMBHs and their emission spectrum, miniquasars powered by IMBHs may be responsible for cosmological reionization at z~15. Reionization by miniquasars with a hard spectrum may be more `economical' than stellar reionization, as soft X-rays escape more easily from the dense sites of star formation and travel further than EUV radiation. Energetic photons will make the diffuse IGM warm and weakly ionized prior to the epoch of reionization breakthrough, set an entropy floor, and reduce gas clumping. Future 21 cm observations may detect a preheated, weakly-ionized IGM in emission against the CMB. (abridged)

We investigate a hierarchical structure formation scenario describing the evolution of a Super Massive Black Holes (SMBHs) population. The seeds of the local SMBHs are assumed to be 'pregalactic' black holes, remnants of the first POPIII stars. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. A simple model in which the damage done to a stellar cusps by decaying BH pairs is cumulative is able to reproduce the observed scaling relation between galaxy luminosity and core size. An accretion model connecting quasar activity with major mergers and the observed BH mass-velocity dispersion correlation reproduces remarkably well the observed luminosity function of optically-selected quasars in the redshift range 1<z<5. We finally asses the potential observability of the gravitational wave background generated by the cosmic evolution of SMBH binaries by the planned space-born interferometer LISA.

We investigate a hierarchical structure formation scenario in which galaxy stellar cores are created from the binding energy liberated by shrinking supermassive black hole (SMBH) binaries. The binary orbital decay heats the surrounding stars, eroding a preexisting 1/r^2 stellar cusp. We follow the merger history of dark matter halos and associated SMBHs via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. Massive black holes get incorporated through a series of mergers into larger and larger halos, sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, and form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce. A simple scheme is applied in which the loss cone is constantly refilled and a constant density core forms due to the ejection of stellar mass. We find that a model in which the effect of the hierarchy of SMBH interactions is cumulative and cores are preserved during galaxy mergers produces at the present epoch a correlation between the `mass deficit' (the mass needed to bring a flat inner density profile to a 1/r^2 cusp) and the mass of the nuclear SMBH, with a normalization and slope comparable to the observed relation. Models in which the mass displaced by the SMBH binary is replenished after every major galaxy merger appear instead to underestimate the mass deficit observed in `core' galaxies.

We assess models for the assembly of supermassive black holes (SMBHs) at the center of galaxies that trace their hierarchical build-up far up in the dark halo `merger tree'. We assume that the first `seed' black holes (BHs) formed in (mini)halos collapsing at z=20 from high-sigma density fluctuations. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. The merger history of dark matter halos and associated BHs is followed by cosmological Monte Carlo realizations of the merger hierarchy. A simple model, where quasar activity is driven by major mergers and SMBHs accrete at the Eddington rate a mass that scales with the fifth power of the velocity dispersion, is shown to reproduce the optical LF of quasars in the redshift range 1<z<4. Binary and triple BH interactions are followed in our merger tree. The assumptions underlying our scenario lead to the prediction of a population of massive BHs wandering in galaxy halos and the intergalactic medium at the present epoch, and contributing <10% to the total BH mass density. At all epochs the fraction of binary SMBHs in galaxy nuclei is of order 10%, while the fraction of binary quasars is less than 0.3%

We present an analysis of the optical data of the Hubble Deep Field South. We derive F300W(AB), F450W(AB), F606W(AB) and F814W(AB) number counts for galaxies in all four bands. The slope is steeper at shortest wavelengths: we estimated gamma(F300W(AB))=0.47, gamma(F450W(AB))=0.35, gamma(F606W(AB))=0.28 and gamma(F814W(AB))=0.28. Morphological number counts are actually dominated by late type galaxies, while early type galaxies show a decreasing slope at faint magnitudes. Combining this information with photometric redshifts, we notice that galaxies contributing with a steep slope to morphological number counts (i.e. spiral and irregular galaxies) have z>1, suggesting a moderate merging. However we emphasize that any cut in apparent magnitude at optical wavelengths results in samples biased against elliptical galaxies, affecting as a consequence the redshift distributions and the implications on the evolution of galaxies along the Hubble sequence.

We describe the construction of a catalogue of galaxies in the optical field of the Hubble Deep Field South. The HDF-S observations produced WFPC2 images in U, B, V, and I, the version 1 data have been made public on 23 November 1999. The effective field of view is 4.38 arcmin$^2$, and the 5$\sigma$ limiting magnitudes (in a FWHM aperture) are 28.87, 29.71, 30.19, 29.58 in the U, B, V and I bands, respectively. We created a catalogue for each pass-band (I$_{814}$, V$_{606}$, B$_{450}$, U$_{300}$), using simulations to account for incompleteness and spurious sources contamination. Along with photometry in all bands, we determined on the I$_{814}$-selected catalogue (I$_{814}<26$) structural parameters, such as a metric apparent size, derived from the petrosian radius, an asymmetry index, light concentration indexes and the mean surface brightness within the petrosian radius.