Researcher profile

Fabio Pacucci

Fabio Pacucci contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2026arXiv

Galaxies in the Epoch of Reionization Are All Bark and No Bite -- Plenty of Ionizing Photons, Low Escape Fractions

Early results from JWST suggest that epoch-of-reionization (EoR) galaxies produce copious ionizing photons, which, if they escape efficiently, could cause reionization to occur too early. We study this problem using \jwst\ imaging and prism spectroscopy for 412 galaxies at 4.5 < z < 9.0. We fit these data simultaneously with stellar-population and nebular-emission models that include a parameter for the fraction of ionizing photons that escape the galaxy, $f_\mathrm{esc}$. We find that the ionization production efficiency, $ξ_\mathrm{ion}$ = Q(H) / L(UV), increases with redshift and decreasing UV luminosity, but shows significant scatter, $σ( \log ξ_\mathrm{ion})$ = 0.3 dex. The inferred escape fractions averaged over the population are low, ranging from $\langle f_\mathrm{esc} \rangle$ = $2.6\pm 1.4$\% at 6 < z < 9 to $6.5\pm 2.2$\% at 4.5 < z < 6 with weak or no indication of evolution with redshift. This implies that in our models most of the ionizing photons need to be absorbed to account for the nebular emission. We compute the impact of our results on reionization, including the distributions for $ξ_\mathrm{ion}$ and $f_\mathrm{esc}$, and the evolution and uncertainty of the UV luminosity function. Considering galaxies brighter than M(UV) < -16 mag, we would produce an IGM hydrogen-ionized fraction of $x_e = 0.5$ at 5.3 < z < 5.8, possibly too late compared to constraints from from QSO sightlines. Including fainter galaxies, M(UV) < -14 mag, we obtain $x_e = 0.5$ at 6.0 < z < 8.1, fully consistent with QSO and CMB data. This implies that EoR galaxies produce plenty of ionizing photons, but these do not efficiently escape. This may be a result of high gas column densities combined with burstier star-formation histories, which limit the time massive stars are able to clear channels through the gas for ionizing photons to escape.

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA&#39;s Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5&#39;&#39; imaging over a wide 24&#39; field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2022arXiv

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of &#34;dark&#34; pixels identified in the Lyman $α$ (Ly$α$) and Lyman $β$ (Ly$β$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$σ$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$σ$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$σ$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018.

preprint2022arXiv

A Search for H-Dropout Lyman Break Galaxies at z~12-16

We present two bright galaxy candidates at z~12-13 identified in our H-dropout Lyman break selection with 2.3 deg2 near-infrared deep imaging data. These galaxy candidates, selected after careful screening of foreground interlopers, have spectral energy distributions showing a sharp discontinuity around 1.7 um, a flat continuum at 2-5 um, and non-detections at <1.2 um in the available photometric datasets, all of which are consistent with z>12 galaxy. An ALMA program targeting one of the candidates shows a tentative 4sigma [OIII]88um line at z=13.27, in agreement with its photometric redshift estimate. The number density of the z~12-13 candidates is comparable to that of bright z~10 galaxies, and is consistent with a recently proposed double power-law luminosity function rather than the Schechter function, indicating little evolution in the abundance of bright galaxies from z~4 to 13. Comparisons with theoretical models show that the models cannot reproduce the bright end of rest-frame ultraviolet luminosity functions at z~10-13. Combined with recent studies reporting similarly bright galaxies at z~9-11 and mature stellar populations at z~6-9, our results indicate the existence of a number of star-forming galaxies at z>10, which will be detected with upcoming space missions such as James Webb Space Telescope, Nancy Grace Roman Space Telescope, and GREX-PLUS.

preprint2022arXiv

Are the Newly-Discovered $z \sim 13$ Drop-out Sources Starburst Galaxies or Quasars?

The detection of two $z\sim 13$ galaxy candidates has opened a new window on galaxy formation at an era only $330$ Myr after the Big Bang. Here, we investigate the physical nature of these sources: are we witnessing star forming galaxies or quasars at such early epochs? If powered by star formation, the observed ultraviolet (UV) luminosities and number densities can be jointly explained if: (i) these galaxies are extreme star-formers with star formation rates $5-24\times$ higher than those expected from extrapolations of average lower-redshift relations; (ii) the star formation efficiency increases with halo mass and is countered by increasing dust attenuation from $z \sim 10-5$; (iii) they form stars with an extremely top-heavy initial mass function. The quasar hypothesis is also plausible, with the UV luminosity produced by black holes of $\sim 10^8 \, \mathrm{M_{\odot}}$ accreting at or slightly above the Eddington rate ($f_{\rm Edd}\sim 1.0$). This black hole mass at $z\sim13$ would require very challenging, but not implausible, growth parameters. If spectroscopically confirmed, these two sources will represent a remarkable laboratory to study the Universe at previously inaccessible redshifts.

preprint2022arXiv

BASS XXII: The BASS DR2 AGN Catalog and Data

We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Doublespec. Many of the spectra have both higher resolution (R>2500, N~450) and/or very wide wavelength coverage (3200-10000 A, N~600) that are important for a variety of AGN and host galaxy studies. We include newly revised AGN counterparts for the full sample and review important issues for population studies, with 44 AGN redshifts determined for the first time and 780 black hole mass and accretion rate estimates. This release is spectroscopically complete for all AGN (100%, 858/858) with 99.8% having redshift measurements (857/858) and 96% completion in black hole mass estimates of unbeamed AGN (outside the Galactic plane). This AGN sample represents a unique census of the brightest hard X-ray selected AGN in the sky, spanning many orders of magnitude in Eddington ratio (Ledd=10^-5-100), black hole mass (MBH=10^5-10^10 Msun), and AGN bolometric luminosity (Lbol=10^40-10^47 ergs/s).

preprint2022arXiv

BAT AGN Spectroscopic Survey XXI: The Data Release 2 Overview

The BAT AGN Spectroscopic Survey (BASS) is designed to provide a highly complete census of the key physical parameters of supermassive black holes (SMBHs) that power local active galactic nuclei (AGN) (z<0.3), including their bolometric luminosity, black hole mass, accretion rates, and line-of-sight gas obscuration, and the distinctive properties of their host galaxies (e.g., star formation rates, masses, and gas fractions). We present an overview of the BASS data release 2 (DR2), an unprecedented spectroscopic survey in spectral range, resolution, and sensitivity, including 1449 optical (3200-10000 A) and 233 NIR (1-2.5 um) spectra for the brightest 858 ultra-hard X-ray (14-195 keV) selected AGN across the entire sky and essentially all levels of obscuration. This release provides a highly complete set of key measurements (emission line measurements and central velocity dispersions), with 99.9% measured redshifts and 98% black hole masses estimated (for unbeamed AGN outside the Galactic plane). The BASS DR2 AGN sample represents a unique census of nearby powerful AGN, spanning over 5 orders of magnitude in AGN bolometric luminosity, black hole mass, Eddington ratio, and obscuration. The public BASS DR2 sample and measurements can thus be used to answer fundamental questions about SMBH growth and its links to host galaxy evolution and feedback in the local universe, as well as open questions concerning SMBH physics. Here we provide a brief overview of the survey strategy, the key BASS DR2 measurements, data sets and catalogs, and scientific highlights from a series of DR2-based works.

preprint2022arXiv

Detectability of Wandering Intermediate-Mass Black Holes in the Milky Way Galaxy from Radio to X-rays

Intermediate-mass black holes (IMBHs, $10^{3-6} \, \rm{M_\odot}$), are typically found at the center of dwarf galaxies and might be wandering, thus far undetected, in the Milky Way (MW). We use model spectra for advection-dominated accretion flows to compute the typical fluxes, in a range of frequencies spanning from radio to X-rays, emitted by a putative population of $10^5 \, \rm{M_\odot}$ IMBHs wandering in five realistic, volume-weighted, MW environments. We predict that $\sim 27\%$ of the wandering IMBHs can be detected in the X-ray with Chandra, $\sim 37\%$ in the near-infrared with the Roman Space Telescope, $\sim 49\%$ in the sub-mm with CMB-S4 and $\sim 57\%$ in the radio with ngVLA. We find that the brightest fluxes are emitted by IMBHs passing through molecular clouds or cold neutral medium, where they are always detectable. We propose criteria to facilitate the selection of candidates in multi-wavelength surveys. Specifically, we compute the X-ray to optical ratio ($α_{\rm ox}$) and the optical to sub-mm ratio, as a function of the accretion rate of the IMBH. We show that at low rates the sub-mm emission of IMBHs is significantly higher than the optical, UV and X-ray emission. Finally, we place upper limits on the number $N_\bullet$ of these objects in the MW: $N_\bullet<2000$ and $N_\bullet<100$, based on our detectability expectations and current lack of detections in molecular clouds and cold neutral medium, respectively. These predictions will guide future searches of IMBHs in the MW, which will be instrumental to understanding their demographics and evolution.

preprint2022arXiv

Dynamics of Intermediate-Mass Black Holes Wandering in the Milky Way Galaxy Using the Illustris TNG50 Simulation

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, \rm M_{\odot}$ to $\sim 5 \times 10^{10} \, \rm M_{\odot}$). IMBHs originally located at the center of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest-resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in a MW analog galaxy. We showed that $\sim 87\%$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \rm ckpc \, h^{-1} \, Gyr^{-1}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic center has a mean of $\sim 180 \, \rm km \, s^{-1}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \rm km \, s^{-1}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys.

preprint2022arXiv

Implications for the Hubble tension from the ages of the oldest astrophysical objects

We use the ages of old astrophysical objects (OAO) in the redshift range $0 \lesssim z \lesssim 8$ as stringent tests of the late-time cosmic expansion history. Since the age of the Universe at any redshift is inversely proportional to $H_0$, requiring that the Universe be older than the oldest objects it contains at any redshift, provides an upper limit on $H_0$. Using a combination of galaxies imaged from the CANDELS program and various high-$z$ quasars, we construct an age-redshift diagram of $\gtrsim 100$ OAO up to $z \sim 8$. Assuming the $Λ$CDM model at late times, we find the 95\% confidence level upper limit $H_0<73.2\,{\rm km}/{\rm s}/{\rm Mpc}$, in slight disagreement with a host of local $H_0$ measurements. Taken at face value, and assuming that the OAO ages are reliable, this suggests that ultimately a combination of pre- and post-recombination ($z \lesssim 10$) new physics might be required to reconcile cosmic ages with early-time and local $H_0$ measurements. In the context of the Hubble tension, our results motivate the study of either combined global pre- and post-recombination modifications to $Λ$CDM, or local new physics which only affects the local $H_0$ measurements.

preprint2022arXiv

Lensing in the Darkness: A Bayesian Analysis of 22 Chandra Sources at $z \gtrsim 6$ Shows No Evidence of Lensing

More than $200$ quasars have been detected so far at $z > 6$, with only one showing clear signs of strong gravitational lensing. Some studies call for a missing population of lensed high-$z$ quasars, but their existence is still in doubt. A large fraction of high-$z$ quasars being lensed would have a significant effect on the shape of the intrinsic quasar luminosity function (QLF). Here, we perform the first systematic search for lensed X-ray-detected quasars at $z \gtrsim 6$ employing a Bayesian analysis, with the code BAYMAX, to look for morphological evidence of multiple images that may escape a visual inspection. We analyzed a sample of 22 quasars at $z > 5.8$ imaged by the Chandra X-ray observatory and found none with statistically significant multiple images. In the sub-sample of the 8 sources with photon counts $>20$ we exclude multiple images with separations $r>1&#39;&#39;$ and count ratios $f>0.4$, or with separations as small as $0.&#39;&#39;7$ and $f>0.7$ at $95\%$ confidence level. Comparing this non-detection with predictions from theoretical models suggesting a high and a low lensed fraction, we placed upper limits on the bright-end slope, $β$, of the QLF. Using only the sub-sample with 8 sources, we obtain, in the high-lensing model, a limit $β< 3.38$. Assuming no multiple source is present in the full sample of 22 sources, we obtain $β< 2.89$ and $β< 3.53$ in the high and low lensing models, respectively. These constraints strongly disfavor steep QLF shapes previously proposed in the literature.

preprint2022arXiv

Prospects to Explore High-redshift Black Hole Formation with Multi-band Gravitational Waves Observatories

The assembly of massive black holes in the early universe remains a poorly constrained open question in astrophysics. The merger and accretion of light seeds (remnants of Population III stars with mass below $\sim 1000M_{\odot}$) or heavy seeds (in the mass range $10^4-10^6 M_{\odot}$) could both explain the formation of massive black holes, but the abundance of seeds and their merging mechanism are highly uncertain. In the next decades, the gravitational-wave observatories coming online are expected to observe very high-redshift mergers, shedding light on the seeding of the first black holes. In this Letter we explore the potential and limitations for LISA, Cosmic Explorer and Einstein Telescope to constrain the mixture ratio of light and heavy seeds as well as the probability that central black holes in merging galaxies merge as well. Since the third generation ground-based gravitational-wave detectors will only observe light seed mergers, we demonstrate two scenarios in which the inference of the seed mixture ratio and merging probability can be limited. The synergy of multi-band gravitational-wave observations and electromagnetic observations will likely be necessary in order to fully characterize the process of high-redshift black hole formation.

preprint2021arXiv

The Search for the Farthest Quasar: Consequences for Black Hole Growth and Seed Models

The quest for high-redshift quasars has led to a series of record-breaking sources, with the current record holder at $z=7.642$. Here, we show how future detections of $z>8$ quasars impact the constraints on the parameters for black hole growth and seed models. Using broad flat priors on the growth parameters (Eddington ratio $\,f_{\rm Edd}$, duty cycle ${\cal D}$, seed mass $M_{\rm \bullet, seed}$ and radiative efficiency $ε$), we show that the large uncertainties in their determination decrease by a factor $\sim 5$ when a quasar&#39;s detection redshift goes from $z=9$ to $z=12$. In this high-redshift regime, $ε$ tends to the lowest value allowed, and the distribution for $M_{\rm \bullet, seed}$ peaks well inside the heavy seed domain. Remarkably, two quasars detected at $z > 7$ with low accretion rates (J1243+0100 and J0313-1806) already tighten the available parameter space, requiring $M_{\rm \bullet, seed} > 10^{3.5} \,{\rm M_\odot}$ and $ε< 0.1$. The radiative efficiency is a crucial unknown, with factor $\sim 2$ changes able to modify the predicted mass by $\sim 3$ orders of magnitude already at $z\sim 9$. The competing roles of inefficient accretion (decreasing $ε$) and black hole spin-up (increasing $ε$) significantly impact growth models. Finally, we suggest that yields currently predicted by upcoming quasar surveys (e.g., Euclid) will be instrumental for determining the most-likely seed mass regime. For example, assuming thin-disk accretion, a detection of a quasar with $M_\bullet \sim 10^{10} \,{\rm M_\odot}$ by $z\sim 9-10$ would exclude the entire parameter space available for light seeds and dramatically reduce the one for heavy seeds.

preprint2020arXiv

Effects of the Hubble Parameter on the Cosmic Growth of the First Quasars

Supermassive black holes (SMBHs) play a crucial role in the evolution of galaxies and are currently detected up to $z\sim 7.5$. Theories describing black hole (BH) growth are challenged by how rapidly seeds with initial mass $M_{\bullet} \lesssim 10^5 \, \mathrm{M_{\odot}}$, formed at $z \sim 20-30$, grew to $M_{\bullet} \sim 10^9 \, \mathrm{M_{\odot}}$ by $z\sim 7$. Here we study the effects of the value of the Hubble parameter, $H_0$, on models describing the early growth of BHs. First, we note that the predicted mass of a quasar at $z=6$ changes by $> 300$% if the underlying Hubble parameter used in the model varies from $H_0 = 65$ to $H_0 = 74$ km s$^{-1}$Mpc$^{-1}$, a range encompassing current estimates. Employing an MCMC approach based on priors from $z \gtrsim 6.5$ quasars and on $H_0$, we study the interconnection between $H_0$ and the parameters describing BH growth: seed mass $M_i$ and Eddington ratio $f_{\rm Edd}$. Assuming an Eddington ratio of $f_{\rm Edd} = 0.7$, in agreement with previous estimates, we find $H_0 = 73.6^{+1.2}_{-3.3}$ km s$^{-1}$Mpc$^{-1}$. In a second analysis, allowing all the parameters to vary freely, we find $\log(M_{i}/\mathrm{M_{\odot}}) > 4.5$ (at 95% CL), $H_0 = 74^{+1.5}_{-1.4}$ km s$^{-1}$Mpc$^{-1}$ and $f_{\rm Edd}=0.77^{+0.035}_{-0.026}$ at 68% CL. Our results on the typical Eddington ratio are in agreement with previous estimates. Current values of the Hubble parameter strongly favour heavy seed formation scenarios, with $M_i \gtrsim 10^4 \, \mathrm{M_{\odot}}$. In our model, with the priors on BH masses of quasars used, light seed formation scenarios are rejected at $\sim 3σ$.

preprint2020arXiv

Finding Direct-Collapse Black Holes at Birth

Direct-collapse black holes (DCBHs) are currently one of the leading contenders for the origins of the first quasars in the universe, over 300 of which have now been found at $z >$ 6. But the birth of a DCBH in an atomically-cooling halo does not by itself guarantee it will become a quasar by $z \sim$ 7, the halo must also be located in cold accretion flows or later merge with a series of other gas-rich halos capable of fueling the BH&#39;s rapid growth. Here, we present near infrared luminosities for DCBHs born in cold accretion flows in which they are destined to grow to 10$^9$ M$_{\odot}$ by $z \sim$ 7. Our observables, which are derived from cosmological simulations with radiation hydrodynamics with Enzo, reveal that DCBHs could be found by the James Webb Space Telescope at $z \lesssim$ 20 and strongly-lensed DCBHs might be found in future wide-field surveys by Euclid and the Wide-Field Infrared Space Telescope at $z \lesssim$ 15.

preprint2020arXiv

Reality or Mirage? Observational Test and Implications for the Claimed Extremely Magnified Quasar at $z = 6.3$

In the last two decades, approximately 200 quasars have been discovered at $z>6$, hosting active super-massive black holes with masses $M_{\bullet} \gtrsim 10^9 M_{\odot}$. While these sources reflect only the tip of the iceberg of the black hole mass distribution, their detection challenges standard growth models. The most massive $z>6$ black hole that was inferred thus far (J0100+2802, $M_{\bullet} \approx 1.2\times 10^{10} M_{\odot}$) was recently claimed to be lensed, with a magnification factor $μ=450$. Here, we perform a consistency check of this claim, finding that the detection of such source requires a bright-end slope $β\geq 3.7$ for the intrinsic quasar luminosity function, $Φ(L) \propto L^{-β}$. Commonly used values of $β\sim 2.8$ are rejected at $>3σ$. If the claim is confirmed, it is very unlikely that all the remaining $51$ sources in the SDSS sample are not magnified. Furthermore, it suffices that $\gtrsim 25\%$ of the remaining sources are lensed for the intrinsic luminosity function to differ significantly (i.e., $>3σ$) from the observed one. The presence of additional extremely magnified sources in the sample would lower the requirement to $\sim 4\%$. Our results urge the community to perform more extended multi-wavelength searches targeting $z>6$ lensed quasars, also among known samples. This effort could vitally contribute to solving the open problem of the growth of the brightest $z\sim 7$ quasars.

preprint2020arXiv

Separating Accretion and Mergers in the Cosmic Growth of Black Holes with X-ray and Gravitational Wave Observations

Black holes across a broad range of masses play a key role in the evolution of galaxies. The initial seeds of black holes formed at $z \sim 30$ and grew over cosmic time by gas accretion and mergers. Using observational data for quasars and theoretical models for the hierarchical assembly of dark matter halos, we study the relative importance of gas accretion and mergers for black hole growth, as a function of redshift ($0<z<10$) and black hole mass ($ 10^3 \, \mathrm{M_{\odot}} < M_{\bullet} < 10^{10} \, \mathrm{M_{\odot}}$). We find that (i) growth by accretion is dominant in a large fraction of the parameter space, especially at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z>6$; and (ii) growth by mergers is dominant at $M_{\bullet} < 10^5 \, \mathrm{M_{\odot}}$ and $z>5.5$, and at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z<2$. As the growth channel has direct implications for the black hole spin (with gas accretion leading to higher spin values), we test our model against $\sim 20$ robust spin measurements available thus far. As expected, the spin tends to decline toward the merger-dominated regime, thereby supporting our model. The next generation of X-ray and gravitational-wave observatories (e.g. Lynx, AXIS, Athena and LISA) will map out populations of black holes up to very high redshift ($z\sim 20)$, covering the parameter space investigated here in almost its entirety. Their data will be instrumental to providing a clear picture of how black holes grew across cosmic time.