Source author record

Fabio Pacucci

Fabio Pacucci 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

28works
7topics
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

28 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's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' 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 "dark" 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''$ and count ratios $f>0.4$, or with separations as small as $0.''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'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'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.

preprint2016arXiv

First Identification of Direct Collapse Black Hole Candidates in the Early Universe in CANDELS/GOODS-S

The first black hole seeds, formed when the Universe was younger than 500 Myr, are recognized to play an important role for the growth of early (z ~ 7) super-massive black holes. While progresses have been made in understanding their formation and growth, their observational signatures remain largely unexplored. As a result, no detection of such sources has been confirmed so far. Supported by numerical simulations, we present a novel photometric method to identify black hole seed candidates in deep multi-wavelength surveys. We predict that these highly-obscured sources are characterized by a steep spectrum in the infrared (1.6-4.5 micron), i.e. by very red colors. The method selects the only 2 objects with a robust X-ray detection found in the CANDELS/GOODS-S survey with a photometric redshift z > 6. Fitting their infrared spectra only with a stellar component would require unrealistic star formation rates (>2000 solar masses per year). To date, the selected objects represent the most promising black hole seed candidates, possibly formed via the direct collapse black hole scenario, with predicted mass >10^5 solar masses. While this result is based on the best photometric observations of high-z sources available to date, additional progress is expected from spectroscopic and deeper X-ray data. Upcoming observatories, like the JWST, will greatly expand the scope of this work.

preprint2015arXiv

Early galaxy formation in warm dark matter cosmologies

We present a framework for high-redshift ($z \geq 7$) galaxy formation that traces their dark matter (DM) and baryonic assembly in four cosmologies: Cold Dark Matter (CDM) and Warm Dark Matter (WDM) with particle masses of $m_x =$ 1.5, 3 and 5 ${\rm keV}$. We use the same astrophysical parameters regulating star formation and feedback, chosen to match current observations of the evolving ultra violet luminosity function (UV LF). We find that the assembly of observable (with current and upcoming instruments) galaxies in CDM and $m_x \geq 3 {\rm keV}$ WDM results in similar halo mass to light ratios (M/L), stellar mass densities (SMDs) and UV LFs. However the suppression of small-scale structure leads to a notably delayed and subsequently more rapid stellar assembly in the $1.5 {\rm keV}$ WDM model. Thus galaxy assembly in $m_x \leq 2 {\rm keV}$ WDM cosmologies is characterized by: (i) a dearth of small-mass halos hosting faint galaxies; and (ii) a younger, more UV bright stellar population, for a given stellar mass. The higher M/L ratio (effect ii) partially compensates for the dearth of small-mass halos (effect i), making the resulting UV LFs closer to CDM than expected from simple estimates of halo abundances. We find that the redshift evolution of the SMD is a powerful probe of the nature of DM. Integrating down to a limit of $M_{UV} =-16.5$ for the James Webb Space Telescope (JWST), the SMD evolves as $\log$(SMD)$\propto -0.63 (1+z)$ in $m_x = 1.5 {\rm keV}$ WDM, as compared to $\log$(SMD)$\propto -0.44 (1+z)$ in CDM. Thus high-redshift stellar assembly provides a powerful testbed for WDM models, accessible with the upcoming JWST.

preprint2015arXiv

Gravitational Waves from Direct Collapse Black Holes Formation

The possible formation of Direct Collapse Black Holes (DCBHs) in the first metal-free atomic cooling halos at high redshifts ($z > 10$) is nowadays object of intense study and several methods to prove their existence are currently under development. The abrupt collapse of a massive ($\sim 10^4 - 10^5 \, \mathrm{M_{\odot}}$) and rotating object is a powerful source of gravitational waves emission. In this work, we employ modern waveforms and the improved knowledge on the DCBHs formation rate to estimate the gravitational signal emitted by these sources at cosmological distances. Their formation rate is very high ($\sim 10^4 \, \mathrm{yr^{-1}}$ up to $z\sim20$), but due to a short duration of the collapse event ($\sim 2-30\, \mathrm{s}$, depending on the DCBH mass) the integrated signal from these sources is characterized by a very low duty-cycle (${\cal D}\sim 10^{-3}$), i.e. a shot-noise signal. Our results show that the estimated signal lies above the foreseen sensitivity of the Ultimate-DECIGO observatory in the frequency range $(0.8-300) \, \mathrm{mHz}$, with a peak amplitude $Ω_{gw} = 1.1 \times 10^{-54}$ at $ν_{max} = 0.9 \, \mathrm{mHz}$ and a peak Signal-to-Noise Ratio $\mathrm{SNR}\sim 22$ at $ν= 20 \, \mathrm{mHz}$. This amplitude is lower than the Galactic confusion noise, generated by binary systems of compact objects in the same frequency band. For this reason, advanced techniques will be required to separate this signal from background and foreground noise components. As a proof-of-concept, we conclude by proposing a simple method, based on the auto-correlation function, to recognize the presence of a ${\cal D} \ll 1$ signal buried into the continuous noise. The aim of this work is to test the existence of a large population of high-z DCBHs, by observing the gravitational waves emitted during their infancy.

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

Simulating the growth of Intermediate Mass Black Holes

Theoretical models predict that a population of Intermediate Mass Black Holes (IMBHs) of mass $M_\bullet \approx 10^{4-5} \, \mathrm{M_{\odot}}$ might form at high ($z > 10$) redshift by different processes. Such objects would represent the seeds out of which $z > 6$ Super-Massive Black Holes (SMBHs) grow. We numerically investigate the radiation-hydrodynamic evolution governing the growth of such seeds via accretion of primordial gas within their parent dark matter halo of virial temperature $T_{vir} \sim 10^4 \, \mathrm{K}$. We find that the accretion onto a Direct Collapse Black Hole (DCBH) of initial mass $M_0=10^5 \, \mathrm{M_{\odot}}$ occurs at an average rate $\dot{M}_{\bullet} \simeq 1.35 \, \dot{M}_{Edd} \simeq 0.1 \, \mathrm{M_{\odot} \, yr^{-1}}$, is intermittent (duty-cycle $ < 50\%$) and lasts $\approx 142 \, \mathrm{Myr}$; the system emits on average at super-Eddington luminosities, progressively becoming more luminous as the density of the inner mass shells, directly feeding the central object, increases. Finally, when $\approx 90\%$ of the gas mass has been accreted (in spite of an average super-Eddington emission) onto the black hole, whose final mass is $\sim 7 \times 10^6 \, \mathrm{M_{\odot}}$, the remaining gas is ejected from the halo due to a powerful radiation burst releasing a peak luminosity $L_{peak}\sim 3\times 10^{45} \, \mathrm{erg \, s^{-1}}$. The IMBH is Compton-thick during most of the evolution, reaching a column density $N_H \sim 10^{25} \, \mathrm{cm^{-2}}$ in the late stages of the simulation. We briefly discuss the observational implications of the model.

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

Essential physics of early galaxy formation

We present a theoretical model embedding the essential physics of early galaxy formation (z = 5-12) based on the single premise that any galaxy can form stars with a maximal limiting efficiency that provides enough energy to expel all the remaining gas, quenching further star formation. This simple idea is implemented into a merger-tree based semi-analytical model that utilises two mass and redshift-independent parameters to capture the key physics of supernova feedback in ejecting gas from low-mass halos, and tracks the resulting impact on the subsequent growth of more massive systems via halo mergers and gas accretion. Our model shows that: (i) the smallest halos (halo mass $M_h \leq 10^{10} M_\odot$) build up their gas mass by accretion from the intergalactic medium; (ii) the bulk of the gas powering star formation in larger halos ($M_h \geq 10^{11.5} M_\odot$) is brought in by merging progenitors; (iii) the faint-end UV luminosity function slope evolves according to $α= -1.75 \log \,z -0.52$. In addition, (iv) the stellar mass-to-light ratio is well fit by the functional form $\log\, M_* = -0.38 M_{UV} -0.13\, z + 2.4$, which we use to build the evolving stellar mass function to compare to observations. We end with a census of the cosmic stellar mass density (SMD) across galaxies with UV magnitudes over the range $-23 \leq M_{UV} \leq -11$ spanning redshifts $5 < z < 12$: (v) while currently detected LBGs contain $\approx 50$% (10%) of the total SMD at $z=5$ (8), the JWST will detect up to 25% of the SMD at $z \simeq 9.5$.

preprint2014arXiv

The X-ray spectra of the first galaxies: 21cm signatures

The cosmological 21cm signal is a physics-rich probe of the early Universe, encoding information about both the ionization and the thermal history of the intergalactic medium (IGM). The latter is likely governed by X-rays from star-formation processes inside very high redshift (z > 15) galaxies. Due to the strong dependence of the mean free path on the photon energy, the X-ray SED can have a significant impact on the interferometric signal from the cosmic dawn. Recent Chandra observations of nearby, star-forming galaxies show that their SEDs are more complicated than is usually assumed in 21cm studies. In particular, these galaxies have ubiquitous, sub-keV thermal emission from the hot interstellar medium (ISM), which generally dominates the soft X-ray luminosity (with energies < 1 keV, sufficiently low to significantly interact with the IGM). Using illustrative soft and hard SEDs, we show that the IGM temperature fluctuations in the early Universe would be substantially increased if the X-ray spectra of the first galaxies were dominated by the hot ISM, compared with X-ray binaries with harder spectra. The associated large-scale power of the 21cm signal would be higher by roughly a factor of three. More generally, we show that the peak in the redshift evolution of the large-scale (k = 0.2 1/Mpc) 21cm power is a robust probe of the soft-band SED of the first galaxies, and importantly, is not degenerate with their bolometric luminosities. On the other hand, the redshift of the peak constrains the X-ray luminosity and halo masses which host the first galaxies.

preprint2013arXiv

Detectability of Free Floating Planets in Open Clusters with JWST

Recent observations have shown the presence of extra-solar planets in Galactic open stellar clusters, as in the Praesepe (M44). These systems provide a favorable environment for planetary formation due to the high heavy-element content exhibited by the majority of their population. The large stellar density, and corresponding high close-encounter event rate, may induce strong perturbations of planetary orbits with large semimajor axes. Here we present a set of N-body simulations implementing a novel scheme to treat the tidal effects of external stellar perturbers on planetary orbit eccentricity and inclination. By simulating five nearby open clusters we determine the rate of occurrence of bodies extracted from their parent stellar system by quasi-impulsive tidal interactions. We find that the specific free-floating planet production rate (total number of free-floating planets per unit of time, normalized by the total number of stars) is proportional to the stellar density of the cluster, with a constant of proportionality equal to (23 +/- 5)10^-6 pc^3 Myr^-1. For the Pleiades (M45) we predict that about 26% of stars should have lost their planets. This raises the exciting possibility of directly observing these wandering planets with the James Webb Space Telescope in the NIR band. Assuming a surface temperature of the planet of 500 K, a free-floating planet of Jupiter size inside the Pleiades would have a specific flux @4.4 micron of approximately 400 nJy, which would lead to a very clear detection (S/N of order 100) in only one hour of integration.

preprint2013arXiv

Focusing on Warm Dark Matter with Lensed High-redshift Galaxies

We propose a novel use of high-redshift galaxies, discovered in deep Hubble Space Telescope (HST) fields around strong lensing clusters. These fields probe small comoving volumes (about 1000 cubic Mpc) at high magnification (μ > 10), and can detect otherwise inaccessible ultra-faint galaxies. Even a few galaxies found in such small volumes require a very high number density of collapsed dark matter (DM) halos. This implies significant primordial power on small scales, allowing these observations to rule out popular alternatives to standard cold dark matter (CDM) models, such as warm dark matter (WDM). In this work, we analytically compute WDM halo mass functions at z = 10, including the effects of both particle free-streaming and residual velocity dispersion. The preliminary number density corresponding to the two galaxies at z about 10 already detected by the Cluster Lensing And Supernova survey with Hubble (CLASH) constrains the WDM particle mass to m_x > 1 (0.9) keV at 84% (99.9%) confidence limit (for a thermal relic relativistic at decoupling). This limit depends only on the WDM halo mass function and, unlike previous constraints on m_x, is independent of any astrophysical modeling. The forthcoming HST Frontier Fields can significantly tighten these constraints.