Chasing the observational signatures of seed black holes at z > 7: candidate statistics

Supermassive black holes (SMBHs) of 10^9-10^10 Msun were already in place ~13 Gyr ago, at z>6. Super-Eddington growth of low-mass BH seeds (~100 Msun) or less extreme accretion onto ~10^5 Msun seeds have been recently considered as the main viable routes to these SMBHs. Here we study the statistics of these SMBH progenitors at z~6. The growth of low- and high-mass seeds and their host galaxies are consistently followed using the cosmological data constrained model GAMETE/QSOdust, which reproduces the observed properties of high-z quasars, like SDSS J1148+5251. We show that both seed formation channels can be in action over a similar redshift range, 15 < z < 18 and are found in dark matter halos with comparable mass, ~5x10^7 Msun. However, as long as the systems evolve in isolation (i.e. no mergers occur), noticeable differences in their properties emerge: at z>= 10 galaxies hosting high-mass seeds have smaller stellar mass and metallicity, the BHs accrete gas at higher rates and star formation proceeds less efficiently than in low-mass seeds hosts. At z<10 these differences are progressively erased, as the systems experience minor or major mergers and every trace of the BH origin gets lost.Comment: accepted for publication in MNRA

Clumping factors of HII, HeII and HeIII

Estimating the intergalactic medium ionization level of a region needs proper treatment of the reionization process for a large representative volume of the universe. The clumping factor, a parameter which accounts for the effect of recombinations in unresolved, small-scale structures, aids in achieving the required accuracy for the reionization history even in simulations with low spatial resolution. In this paper, we study for the first time the redshift evolution of clumping factors of different ionized species of H and He in a small but very high resolution simulation of the reionization process. We investigate the dependence of the value and redshift evolution of clumping factors on their definition, the ionization level of the gas, the grid resolution, box size and mean dimensionless density of the simulations.Comment: 12 pages, 10 figures, 1 table. Accepted by MNRA

Chasing the observational signatures of seed black holes at z > 7: Candidate statistics

Supermassive black holes (SMBHs) of 109–1010M were already in place ∼13 Gyr ago, at z > 6. Super-Eddington growth of low-mass BH seeds (∼100M) or less extreme accretion on to∼105M seeds have been recently considered as the main viable routes to these SMBHs. Here, we study the statistics of these SMBH progenitors at z ∼ 6. The growth of low- and high-mass seeds and their host galaxies are consistently followed using the cosmological data constrained model GAMETE/QSODUST, which reproduces the observed properties of high-z quasars, like SDSS J1148+5251.We show that both seed formation channels can be in action over a similar redshift range 15 < z < 18 and are found in dark matter haloes with comparable mass, ∼5 × 107M. However, as long as the systems evolve in isolation (i.e. no mergers occur), noticeable differences in their properties emerge: At z ≥ 10 galaxies hosting high-mass seeds have smaller stellar mass and metallicity, the BHs accrete gas at higher rates and star formation proceeds less efficiently than in low-mass seeds hosts. At z < 10 these differences are progressively erased, as the systems experience minor or major mergers and every trace of the BH origin gets lost

Radiative transfer of ionizing radiation through gas and dust: stellar source case

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H II regions and on small scale reionization. We find that H II regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions.Comment: 17 pages, 14 figures. Accepted for publication in MNRA

X-ray ionization of the intergalactic medium by quasars

We investigate the impact of quasars on the ionization of the surrounding intergalactic medium (IGM) with the radiative transfer code \texttt{CRASH4}, now accounting for X-rays and secondary electrons. After comparing with analytic solutions, we post-process a cosmic volume ($\approx 1.5\times 10^4$ Mpc$^3 h^{-3}$) containing a ULAS J1120+0641-like quasar (QSO) hosted by a $5 \times 10^{11} {\rm M}_\odot h^{-1}$ dark matter (DM) halo. We find that: (i) the average HII region ($R\sim3.2$~pMpc in a lifetime $t_f = 10^7$~yrs) is mainly set by UV flux, in agreement with semi-analytic scaling relations; (ii) a largely neutral ($x_{\textrm{HII}} < 0.001$), warm ($T\sim 10^3$~K) tail extends up to few Mpc beyond the ionization front, as a result of the X-ray flux; (iii) LyC-opaque inhomogeneities induce a line of sight (LOS) scatter in $R$ as high as few physical Mpc, consistent with the DLA scenario proposed to explain the anomalous size of the ULAS J1120+0641 ionized region. On the other hand, with an ionization rate $\dot{N}_{\gamma,0} \sim 10^{57}$~s$^{-1}$, the assumed DLA clustering and gas opacity, only one LOS shows an HII region compatible with the observed one. We deduce that either the ionization rate of the QSO is at least one order of magnitude lower or the ULAS J1120+0641 bright phase is shorter than $10^7$~yrs.Comment: Accepted for publication in MNRAS Main Journal, Accepted 2018 May 2

Chasing the observational signatures of seed black holes at z > 7: candidate observability

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

Where does galactic dust come from?

Here we investigate the origin of the dust mass (Mdust) observed in the Milky Way (MW) and of dust scaling relations found in a sample of local galaxies from the DGS and KINGFISH surveys. To this aim, we model dust production from Asymptotic Giant Branch (AGB) stars and supernovae (SNe) in simulated galaxies forming along the assembly of a Milky Way-like halo in a well resolved cosmic volume of 4cMpc using the GAMESH pipeline. We explore the impact of different sets of metallicity and mass-dependent AGB and SN dust yields on the predicted Mdust. Our results show that models accounting for grain destruction by the SN reverse shock predict a total dust mass in the MW that is a factor of ~4 lower than observed, and can not reproduce the observed galaxy-scale relations between dust and stellar masses, and dust-to-gas ratios and metallicity, with a smaller discrepancy in galaxies with low metallicity (12 + log(O/H) < 7.5) and low stellar masses (Mstar < 10^7 Msun). In agreement with previous studies, we suggest that competing processes in the interstellar medium must be at play to explain the observed trends. Our result reinforces this conclusion by showing that it holds independently of the adopted AGB and SN dust yields.Comment: 6 pages, 4 figures. Accepted version for publication in MNRA

The dust mass in z > 6 normal star forming galaxies

We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10^9 Msun allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy populationComment: 5 pages, 2 figures, accepted for publication as a Letter to MNRA

Radio recombination lines from obscured quasars with the SKA

We explore the possibility of detecting hydrogen radio recombination lines from 0 < z < 10 quasars. We compute the expected Hnalpha flux densities as a function of absolute magnitude and redshift by considering (i) the range of observed AGN spectral indices from UV to X-ray bands, (ii) secondary ionizations from X-ray photons, and (iii) stimulated emission due to nonthermal radiation. All these effects are important to determine the line fluxes. We find that the combination of slopes: alpha_X,hard = -1.11, alpha_X,soft = -0.7, alpha_EUV = -1.3, alpha_UV = -1.7, maximizes the expected flux, f_Hnalpha = 10 microJy for z = 7 quasars with M_AB = -27 in the n = 50 lines; allowed SED variations produce variations by a factor of 3 around this value. Secondaries boost the line intensity by a factor of 2 to 4, while stimulated emission in high-z quasars with M_AB = -26 provides an extra boost to RRL flux observed at nu = 1 GHz if recombinations arise in HII regions with T_e = 10^3-5 K, n_e = 10^3-5 cm^-3. We compute the sensitivity required for a 5sigma detection of Hnalpha lines using the SKA, finding that the SKA-MID could detect sources with M_AB < -27 (M_AB < -26) at z < 8 (z < 3) in less than 100 hrs of observing time. These observations could open new paths to searches for obscured SMBH progenitors, complementing X-ray, optical/IR and sub-mm surveys.Comment: 11 pages, 9 figures; to be published in Monthly Notices of the Royal Astronomical Society Main Journa