Here we explore the infrared (IR) properties of the progenitors of high-z quasar host galaxies. Adopting the cosmological, data constrained semi-analytic model GAMETE/QSOdust, we simulate several independent merger histories of a luminous quasar at z ~ 6, following black hole growth and baryonic evolution in all its progenitor galaxies. We find that a fraction of progenitor galaxies (about 0.4 objects per single luminous quasar) at 6.5 < z < 8 has an IR luminosity of L_IR > 10^13 Lsun (hyper-luminous IR galaxies; HyLIRGs). HyLIRGs progenitors reside in the most massive halos, with dark matter (DM) masses of M_DM ~ 10^12.5 - 10^13 Msun. These systems can be easily observed in their ~ 1 mm-continuum emission in a few seconds of integration time with the Atacama Large Millimeter/submillimeter Array (ALMA), and at least 40% of them host nuclear BH activity that is potentially observable in the soft and hard X-ray band. Our findings are in line with recent observations of exceptional massive DM halos hosting HyLIRGs at z ~ 7, suggesting that z ~ 6 luminous quasars are indeed the signposts of these observed rare peaks in the high-z cosmic density field, and that massive IR-luminous galaxies at higher z are their natural ancestors.

Star-forming galaxies at high redshift show anomalous values of infrared excess, which can be described only by extremizing the existing relations between the shape of their ultraviolet continuum emission and their infrared-to-ultraviolet luminosity ratio, or by constructing \textit{ad-hoc} models of star formation and dust distribution. We present an alternative explanation, based on unveiled AGN activity, to the existence of such galaxies. In fact, the presence of a weak AGN configures as a natural scenario in order to explain the observed spectral properties of such high-$z$ objects in terms of a continuum slope distribution rather than altered infrared excesses, due to the different shape of the AGN continuum emission with respect to quiescent galaxies. To this aim, we directly compare the infrared-to-ultraviolet properties of high-redshift galaxies to those of known categories of AGN (quasars and Seyferts). We also infer the characteristics of their possible X-ray emission. We find a strong similarity between the spectral shapes and luminosity ratios of AGN with the corresponding properties of such galaxies. In addition, we derive expected X-ray fluxes that are compatible with energetics from AGN activity. We conclude that a moderate AGN contribution to the UV emission of such high-$z$ objects is a valid alternative to explain their spectral properties. Even the presence of an active nucleus in each source would not violate the expected quasar statistics. Furthermore, we suggest that the observed similarities between anomalous star-forming galaxies and quasars may provide a benchmark for future theoretical and observational studies on the galaxy population in the early Universe.

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.

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

Observations of the most luminous quasars at redshift z>6 reveal the existence of numerous supermasssive black holes (>10^9 Msun) already in place about twelve billion years ago. In addition, the interstellar medium of the galaxies hosting these black holes are observed to be chemically mature systems, with metallicities (Z>Zsun) and dust masses (>10^8 Msun) similar to that of more evolved, local galaxies. The connection between the rapid growth of the first supermassive black holes and the fast chemical evolution of the host galaxy is one of the most puzzling issues for theoretical models. Here we review state-of-the-art theoretical models that focus on this problem with particular emphasis on the conditions that lead to the formation of quasar seeds and their subsequent evolution at z>6.

Observational searches for faint active nuclei at $z > 6$ have been extremely elusive, with a few candidates whose high-$z$ nature is still to be confirmed. Interpreting this lack of detections is crucial to improve our understanding of high-$z$ supermassive black holes (SMBHs) formation and growth. In this work, we present a model for the emission of accreting BHs in the X-ray band, taking into account super-Eddington accretion, which can be very common in gas-rich systems at high-$z$. We compute the spectral energy distribution for a sample of active galaxies simulated in a cosmological context, which represent the progenitors of a $z \sim 6$ SMBH with $M_{\rm BH} \sim 10^9 \, M_\odot$. We find an average Compton thick fraction of $\sim 45\%$ and large typical column densities ($N_H \gtrsim 10^{23} \rm \, cm^2$). However, faint progenitors are still luminous enough to be detected in the X-ray band of current surveys. Even accounting for a maximum obscuration effect, the number of detectable BHs is reduced at most by a factor 2. In our simulated sample, observations of faint quasars are mainly limited by their very low active fraction ($f_{\rm act} \sim 1 \%$), which is the result of short, super-critical growth episodes. We suggest that to detect high-$z$ SMBHs progenitors, large area surveys with shallower sensitivities, such as Cosmos Legacy and XMM-LSS+XXL, are to be preferred with respect to deep surveys probing smaller fields, such as CDF-S.

The assembly of the first super massive black holes (SMBHs) at $z \gtrsim 6$ is still a subject of intense debate. If black holes (BHs) grow at their Eddington rate, they must start from $\gtrsim 10^4 \, M_\odot$ seeds formed by the direct collapse of gas. Here we explore the alternative scenario where $\sim 100 \, M_\odot$ BH remnants of the first stars grow at super-Eddington rate via radiatively inefficient slim accretion disks. We use an improved version of the cosmological, data-constrained semi-analytic model GAMETE/QSOdust, where we follow the evolution of nuclear BHs and gas cooling, disk and bulge formation of their host galaxies. Adopting SDSS J1148+5251 (J1148) at $z = 6.4$ as a prototype of luminous $z \gtrsim 6$ quasars, we find that $\sim$ 80% of its SMBH mass is grown by super-Eddington accretion, which can be sustained down to $z \sim 10$ in dense, gas-rich environments. The average BH mass at $z \sim 20$ is $M_{\rm BH} \gtrsim 10^4 \,M_\odot$, comparable to that of direct collapse BHs. At $z = 6.4$ the AGN-driven mass outflow rate is consistent with the observations and the BH-to-bulge mass ratio is compatible with the local scaling relation. However, the stellar mass in the central 2.5 kpc is closer to the value inferred from CO observations. Finally, $\sim 20 %$ of J1148 progenitors at $z=7.1$ have BH luminosities and masses comparable to ULAS J1120+0641, suggesting that this quasar may be one of the progenitors of J1148.