Source author record

Marta Volonteri

Marta Volonteri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2025arXiv

Super-Eddington accretion in high-redshift quasar hosts: Black-hole driven outflows, galaxy quenching, and the nature of little red dots

The advent of the James Webb Space Telescope has revolutionised our understanding of the high-redshift Universe through its detection of bright, massive galaxies up to $z\gtrsim 10$ and its identification of peculiar sources called `little red dots' (LRDs). The origin of both classes of objects remains uncertain but is likely linked to the formation and early growth of the first massive black holes (MBHs), which may be more easily explained by invoking phases of super-Eddington accretion. In this study, we used a state-of-the-art zoom-in cosmological simulation of a quasar host to investigate whether these objects could resemble any of the peculiar sources observed with JWST during their assembly. We find that the impact of MBH feedback on star formation is typically moderate, with outflows preferentially escaping perpendicular to the galactic disc. However, for approximately ten percent of the galaxy's lifetime, the system enters a distinct quenched phase following rapid MBH growth driven by super-Eddington accretion. This phase culminates in a powerful feedback event, during which the MBH jet and disc-driven winds interact directly with the galactic disc and carve out a central cavity. We also find that, during the history of the quasar host progenitor, the spectral properties of the system can resemble both LRDs and quenched galaxies, depending on the specific evolutionary stage considered. These findings suggest that both conditions may represent transient phases in the life cycle of high-redshift galaxies.

preprint2022arXiv

Concordance between observations and simulations in the evolution of the mass relation between supermassive black holes and their host galaxies

We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars ($0.2 < z < 0.8$) and 32 broad-line active galactic nuclei (AGNs; $1.2<z<1.7$) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations which appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift ($z\sim$ 0.2 and 1.5, respectively) while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects.

preprint2022arXiv

Dual AGN in the Horizon-AGN simulation and their link to galaxy and massive black hole mergers, with an excursus on multiple AGN

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

preprint2022arXiv

Forming massive seed black holes in high-redshift quasar host progenitors

The presence of massive black holes (BHs) with masses of order $10^9\rm\, M_\odot$, powering bright quasars when the Universe was less than 1 Gyr old, poses strong constraints on their formation mechanism. Several scenarios have been proposed to date to explain massive BH formation, from the low-mass seed BH remnants of the first generation of stars to the massive seed BHs resulting from the rapid collapse of massive gas clouds. However, the plausibility of some of these scenarios to occur within the progenitors of high-z quasars has not yet been thoroughly explored. In this work, we investigate, by combining dark-matter only N-body simulations with a semi-analytic framework, whether the conditions for the formation of massive seed BHs from synchronised atomic-cooling halo pairs and/or dynamically-heated mini-haloes are fulfilled in the overdense regions where the progenitors of a typical high-redshift quasar host form and evolve. Our analysis shows that the peculiar conditions in such regions, i.e. strong halo clustering and high star formation rates, are crucial to produce a non-negligible number of massive seed BH host candidates: we find $\approx1400$ dynamically heated metal-free mini-haloes, including one of these which evolves to a synchronised pair and ends up in the massive quasar-host halo by $z=6$. This demonstrates that the progenitors of high-redshift quasar host haloes can harbour early massive seed BHs. Our results further suggest that multiple massive seed BHs may form in or near the quasar host's progenitors, potentially merging at lower redshifts and yielding gravitational wave events.

preprint2022arXiv

High-redshift quasars and their host galaxies II: multiphase gas and stellar kinematics

Observations of $z \gtrsim 6$ quasars provide information on the early phases of the most massive black holes (MBHs) and galaxies. Current observations at sub-mm wavelengths trace cold and warm gas, and future observations will extend information to other gas phases and the stellar properties. The goal of this study is to examine the gas life cycle in a $z \gtrsim 6$ quasar: from accretion from the halo to the galaxy and all the way into the MBH, to how star formation and the MBH itself affect the gas properties. Using 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, we investigate the distribution of the different gas phases in the interstellar medium across cosmic time. We assess the morphological evolution of the quasar host using different tracers (star- or gas-based) and the thermodynamic distribution of the MBH accretion-driven outflows, finding that obscuration in the disc is mainly due to molecular gas, with the atomic component contributing at larger scales and/or above/below the disc plane. Moreover, our results also show that molecular outflows, if present, are more likely the result of gas being lifted near the MBH than production within the wind because of thermal instabilities. Finally, we also discuss how different gas phases can be employed to dynamically constrain the MBH mass, and argue that resolutions below $\sim 100$ pc yield unreliable estimates because of the strong contribution of the nuclear stellar component to the potential at larger scales.

preprint2022arXiv

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

The semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Press-Schechter theory or that can be imported from dark matter simulations. On top of these merger trees, A-SLOTH applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. A-SLOTH samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star formation rate density at high redshift. A-SLOTH has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high-z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early Universe, such as supernova rates or gravitational wave events.

preprint2022arXiv

Seeking the growth of the first black hole seeds with JWST

In this paper we provide predictions for the BH population that would be observable with planned JWST surveys at $5 \le z \le 15$. We base our study on the recently developed Cosmic Archaeology Tool (CAT), which allows us to model BH seeds formation and growth, while being consistent with the general population of AGNs and galaxies observed at $4 \le z \le 7$. We find that JWST planned surveys will provide a complementary view on active BHs at $z > 5$, with JADES-Medium/-Deep being capable of detecting the numerous BHs that populate the faint-end of the distribution, COSMOS-Web sampling a large enough area to detect the rarest brightest systems, and CEERS/PRIMER bridging the gap between these two regimes. The relatively small field of view of the above surveys preferentially selects BHs with masses $6 \leq {\rm Log} (M_{\rm BH}/M_\odot) < 8$ at $7 \le z < 10$, residing in relatively metal poor (${\rm Log} (Z/Z_\odot) \ge -2$) and massive ($8\leq {\rm Log} (M_*/M_\odot) < 10$) galaxies. At $z \ge 10$, only JADES-Deep will have the sensitivity to detect growing BHs with masses $4 \leq {\rm Log} (M_{\rm BH}/M_\odot) < 6$, hosted by more metal poor ($-3 \leq {\rm Log} (Z/Z_\odot) < -2$) and less massive ($6 \leq {\rm Log} (M_*/M_\odot) < 8$) galaxies. In our model, the latter population corresponds to heavy BH seeds formed by the direct collapse of super-massive stars in their earliest phases of mass growth. Detecting these systems would provide invaluable insights on the nature and early growth of the first BH seeds.

preprint2022arXiv

The black hole population in low-mass galaxies in large-scale cosmological simulations

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

preprint2022arXiv

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

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$.

preprint2022arXiv

The origins of massive black holes

Massive black holes (MBHs) inhabit galaxy centers, power luminous quasars and Active Galactic Nuclei (AGN) and shape their cosmic environment with the energy they produce. The origins of MBHs remain a mystery and the recent detection by LIGO/Virgo of an almost 150 solar mass black hole has revitalized the question of whether there is a continuum between "stellar" and "massive" black holes and what the seeds of MBHs are. Seeds could have formed in the first galaxies, or could be also related to the collapse of horizon-sized regions in the early Universe. Understanding the origins of MBHs straddles fundamental physics, cosmology and astrophysics and it bridges the fields of gravitational wave physics and traditional astronomy. With several facilities in the next 10-15 years we foresee the possibility of discovering MBHs' avenues of formation. In this article we link three main topics: the channels of black hole seed formation, the journey from seeds to massive black holes, the diagnostics on the origins of MBHs. We highlight and critically discuss current unsolved problems and touch on recent developments that stirred the community.

preprint2020arXiv

Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters

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.

preprint2020arXiv

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

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. We use the same method to analyse the two simulations, post-processing 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.

preprint2020arXiv

Modelling a bright z = 6 galaxy at the faint end of the AGN luminosity function

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.

preprint2020arXiv

Real galaxy mergers from galaxy pair catalogs

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 $Δz$) gives similar results than the most complex one (based on a neural network analyzing $d$, $Δ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 $Δz_\mathrm{th} \sim 10^{-3}$, in agreement with recent results.

preprint2020arXiv

Reionization with galaxies and active galactic nuclei

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$.

preprint2020arXiv

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

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.

preprint2016arXiv

Exploring the nature of the Lyman-$α$ emitter CR7

CR7 is the brightest Lyman-$α$ emitter observed at $z>6$, which shows very strong Lyman-$α$ and HeII 1640Å 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.

preprint2016arXiv

From the first stars to the first black holes

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.

preprint2016arXiv

Gravitational Waves from the Remnants of the First Stars

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.

preprint2016arXiv

How to quench a galaxy

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 $Λ$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.

preprint2016arXiv

Inferences on the Relations Between Central Black Hole Mass and Total Galaxy Stellar Mass in the high-redshift Universe

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.

preprint2016arXiv

The cosmic evolution of massive black holes in the Horizon-AGN simulation

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.

preprint2016arXiv

The Horizon-AGN simulation: morphological diversity of galaxies promoted by AGN feedback

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.

preprint2015arXiv

A comparison of black hole growth in galaxy mergers with Gasoline and Ramses

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.

preprint2015arXiv

Black hole accretion versus star formation rate: theory confronts observations

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.

preprint2015arXiv

Black hole evolution: I. Supernova-regulated black hole growth

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.

preprint2015arXiv

Black hole formation and growth with non-Gaussian primordial density perturbations

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.

preprint2015arXiv

Growing black holes and galaxies: black hole accretion versus star formation rate

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.

preprint2015arXiv

How an improved implementation of H2 self-shielding influences the formation of massive stars and black holes

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.

preprint2015arXiv

Massive black holes in merging galaxies

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.

preprint2015arXiv

Off the Beaten Path: A New Approach to Realistically Model The Orbital Decay of Supermassive Black Holes in Galaxy Formation Simulations

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.

preprint2015arXiv

Relations Between Central Black Hole Mass and Total Galaxy Stellar Mass in the Local Universe

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.

preprint2015arXiv

Shining in the Dark: the Spectral Evolution of the First Black Holes

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{μm}$) and X-ray ($0.1 - 100 \, \mathrm{keV}$) bands. Such signal should be easily detectable by JWST around $\sim 1 \, \mathrm{μ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 $ρ_{\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, $ρ_{\bullet} \lesssim 7.6 \times 10^{3} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ in the most constraining case.

preprint2015arXiv

The case for super-critical accretion onto massive black holes at high redshift

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.

preprint2015arXiv

The evolution of high-redshift massive black holes

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.

preprint2015arXiv

The Growth Efficiency of High-Redshift Black Holes

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 ($ε= 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 ($ε\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.

preprint2014arXiv

Black hole evolution: II. Spinning black holes in a supernova-driven turbulent interstellar medium

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.

preprint2014arXiv

Black hole evolution: III. Statistical properties of mass growth and spin evolution using large-scale hydrodynamical cosmological simulations

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.

preprint2014arXiv

Constraining the high redshift formation of black hole seeds in nuclear star clusters with gas inflows

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.

preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

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.

preprint2014arXiv

Growth and activity of black holes in galaxy mergers with varying mass ratios

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.

preprint2014arXiv

New insights on the recoiling/binary black hole candidate J0927+2943 via molecular gas observations

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.

preprint2014arXiv

No WIMP Mini-Spikes in Dwarf Spheroidal Galaxies

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-σ 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.

preprint2014arXiv

Nuclear coups: dynamics of black holes in galaxy mergers

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.

preprint2013arXiv

Effects Of Circum-nuclear Disk Gas Evolution And The Spin Of Central Black Holes

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 ε= 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.

preprint2013arXiv

Massive black holes in central cluster galaxies

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.

preprint2013arXiv

The evolution of active galactic nuclei and their spins

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.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

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.

preprint2013arXiv

The Hot and Energetic Universe: The formation and growth of the earliest supermassive black holes

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.

preprint2012arXiv

Black Holes in the Early Universe

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.

preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

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.

preprint2012arXiv

Gravitational Recoil From Accretion-Aligned Black-Hole Binaries

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.

preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

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.

preprint2012arXiv

Observability of Dual Active Galactic Nuclei in Merging Galaxies

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.

preprint2012arXiv

The Formation and Evolution of Massive Black Holes

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.

preprint2012arXiv

The mass function of black holes at 1<z<4.5: comparison of models with observations

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).

preprint2011arXiv

Assessing the redshift evolution of massive black holes and their hosts

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.

preprint2011arXiv

Black hole growth in the early Universe is self-regulated and largely hidden from view

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.

preprint2011arXiv

Blazars in the early Universe

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.

preprint2011arXiv

Blindly Detecting Merging Supermassive Black Holes with Radio Surveys

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.

preprint2011arXiv

How important is the dark matter halo for black hole growth?

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.

preprint2011arXiv

Massive black holes in stellar systems: 'quiescent' accretion and luminosity

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.

preprint2011arXiv

The First Massive Black Hole Seeds and Their Hosts

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.

preprint2010arXiv

Constraining properties of the black hole population using LISA

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.

preprint2010arXiv

Evolution of massive black hole spins

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.

preprint2010arXiv

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

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.

preprint2010arXiv

Formation of Supermassive Black Holes

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.

preprint2010arXiv

Gravitational recoil: effects on massive black hole occupation fraction over cosmic time

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).

preprint2010arXiv

LISA extreme-mass-ratio inspiral events as probes of the black hole mass function

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.

preprint2010arXiv

Quasistars and the cosmic evolution of massive black holes

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.

preprint2010arXiv

Reconstructing the massive black hole cosmic history through gravitational waves

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.

preprint2010arXiv

Wandering Black Holes in Bright Disk Galaxy Halos

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.

preprint2009arXiv

Journey to the M_BH -sigma relation: the fate of low mass black holes in the Universe

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} - σ$ 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} - σ$ 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'.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

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.

preprint2009arXiv

The early evolution of massive black holes

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.

preprint2008arXiv

Cosmological black hole spin evolution by mergers and accretion

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.

preprint2008arXiv

Double quasars: probes of black hole scaling relationships and merger scenarios

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.

preprint2006arXiv

Evolution of supermassive black holes

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.

preprint2004arXiv

The distribution and cosmic evolution of massive black hole spins

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).