We present asymptotic giant branch (AGB) models of metallicity $Z=10^{-4}$ and $Z=3\times 10^{-4}$, with the aim of understanding how the gas enrichment and the dust production change in very metal-poor environments and to assess the general contribution of AGB stars to the cosmic dust yield. The stellar yields and the dust produced are determined by the change in the surface chemical composition, with a transition occurring at $\sim 2.5~M_{\odot}$. Stars of mass $M < 2.5~M_{\odot}$ reach the carbon stage and produce carbon dust, whereas their higher mass counterparts produce mainly silicates and alumina dust; in both cases the amount of dust manufactured decreases towards lower metallicities. The $Z=10^{-4}$ models show a complex and interesting behaviour on this side, because the efficient destruction of the surface oxygen favours the achievement of the C-star stage, independently of the initial mass. The present results might indicate that the contribution from this class of stars to the overall dust enrichment in metal-poor environments is negligible at redshifts $z>5$.

Next-generation observations will revolutionize our understanding of binary black holes and will detect new sources, such as intermediate-mass black holes. Primary science goals include: Discover binary black holes throughout the observable Universe; Reveal the fundamental properties of black holes; Uncover the seeds of supermassive black holes.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

We present a 155ks NuSTAR observation of the $z\sim2$ hot dust-obscured galaxy (hot DOG) W1835+4355. We extracted spectra from the two NuSTAR detectors and analyzed them jointly with the archival XMM PN and MOS spectra. We performed a spectroscopic analysis based on both phenomenological and physically motivated models employing toroidal and spherical geometry for the obscurer. In all the modelings, the source exhibits a Compton-thick column density $N_{\rm H} \gtrsim 10^{24}$ cm$^{-2}$, a 2-10 keV luminosity $L_{2-10}\approx2\times10^{45}$ erg s$^{-1}$ , and a prominent soft excess ($\sim5-10$ % of the primary radiative output), which translates into a luminosity $\sim10^{44}$ erg s$^{-1}$. We modeled the spectral energy distribution from 1.6 to 850 $\mu m$ using a clumpy two-phase dusty torus model plus a modified blackbody to account for emission powered by star formation in the far-infrared. We employed several geometrical configurations consistent with those applied in the X-ray analysis. In all cases we obtained a bolometric luminosity $L_{\rm bol}\approx3-5\times10^{47}$ erg s$^{-1}$, which confirms the hyperluminous nature of this active galactic nucleus. Finally, we estimate a prodigious star formation rate of $\sim$3000 $M_{\odot}\,yr^{-1}$, which is consistent with the rates inferred for $z\approx2-4$ hyperluminous type I quasars. The heavily obscured nature, together with $L_{\rm bol}$, the ratio of X-ray to mid-infrared luminosity, the rest-frame optical morphology, and the host star formation rate are indicative of its evolutionary stage. We can interpret this as a late-stage merger event in the transitional, dust-enshrouded, evolutionary phase eventually leading to an optically bright AGN.

We present an ALMA high-resolution observation of the 840 um continuum and [CII] line emission in the WISE-SDSS selected hyper-luminous (WISSH) QSO J1015+0020 at z~4.4. Our analysis reveals an exceptional overdensity of [CII]-emitting companions with a very small (<150 km/s) velocity shift with respect to the QSO redshift. We report the discovery of the closest companion observed so far in submillimetre observations of high-z QSOs. It is only 2.2 kpc distant and merging with J1015+0020, while two other [CII] emitters are found at 8 and 17 kpc. Two strong continuum emitters are also detected at <3.5 arcsec. They are likely associated to the same overdense structure of J1015+0020, as they exceed by a factor of 100 the number of expected sources, considering the Log(N)-Log(S) at 850 um. The host galaxy of J1015+0020 shows a SFR of about 100 Msun/yr while the total SFR of the QSO and its companion galaxies is a factor of 10 higher, indicating that substantial stellar mass assembly at early epochs may have taken place in the QSO satellites. For J1015+0020 we compute a SMBH mass MBH~6E9 Msun and a dynamical mass Mdyn~4E10 Msun . This translates into an extreme ratio Mdyn/MBH~7. The total stellar mass of the QSO host galaxy plus the [CII] emitters already exceeds 1E11 Msun at z~4.4. These sources will likely merge and develop into a giant galaxy of 1.3E12 Msun. Under the assumption of constant mass or Eddington accretion rate equal to the observed values, we find that the growth timescale of the host galaxy is comparable or even shorter than that inferred for the SMBH.

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.

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.

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.

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.

Studying the coupling between the energy output produced by the central quasar and the host galaxy is fundamental to fully understand galaxy evolution. Quasar feedback is indeed supposed to dramatically affect the galaxy properties by depositing large amounts of energy and momentum into the ISM. In order to gain further insights on this process, we study the SEDs of sources at the brightest end of the quasar luminosity function, for which the feedback mechanism is supposed to be at its maximum. We model the rest-frame UV-to-FIR SEDs of 16 WISE-SDSS Selected Hyper-luminous (WISSH) quasars at 1.8 < z < 4.6 disentangling the different emission components and deriving physical parameters of both the nuclear component and the host galaxy. We also use a radiative transfer code to account for the contribution of the quasar-related emission to the FIR fluxes. Most SEDs are well described by a standard combination of accretion disk+torus and cold dust emission. However, about 30% of them require an additional emission component in the NIR, with temperatures peaking at 750K, which indicates the presence of a hotter dust component in these powerful quasars. We measure extreme values of both AGN bolometric luminosity (LBOL > 10^47 erg/s) and SFR (up to 2000 Msun/yr). A new relation between quasar and star-formation luminosity is derived (LSF propto LQSO^(0.73)) by combining several Herschel-detected quasar samples from z=0 to 4. Future observations will be crucial to measure the molecular gas content in these systems, probe the impact between quasar-driven outflows and on-going star-formation, and reveal the presence of merger signatures in their host galaxies.

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.

We present dust yields for asymptotic giant branch (AGB) and super--asymptotic giant branch (SAGB) stars of solar metallicity. Stars with initial mass $1.5~M_{\odot} \leq M_{\rm ini} \leq 3~M_{\odot}$ reach the carbon star stage during the AGB phase and produce mainly solid carbon and SiC. The size and the amount of the carbon particles formed follows a positive trend with themass of the star; the carbon grains with the largest size ($a_{\rm C} \sim 0.2\mu$m) are produced by AGB stars with $M_{\rm ini} = 2.5-3~M_{\odot}$, as these stars are those achieving the largest enrichment of carbon in the surface regions. The size of SiC grains, being sensitive to the surface silicon abundance, keeps around $a_{\rm SiC} \sim 0.1\mu$m. The mass of carbonaceous dust formed is in the range $10^{-4} - 5\times 10^{-3}~M_{\odot}$, whereas the amount of SiC produced is $2\times 10^{-4} - 10^{-3}~M_{\odot}$. Massive AGB/SAGB stars with $M_{\rm ini} > 3~M_{\odot}$ experience HBB, that inhibits formation of carbon stars. The most relevant dust species formed in these stars are silicates and alumina dust, with grain sizes in the range $0.1\mu m < a_{\rm ol} < 0.15\mu$m and $a_{\rm Al_2O_3} \sim 0.07\mu$m, respectively. The mass of silicates produced spans the interval $3.4\times 10^{-3}~M_{\odot} \leq M_{\rm dust} \leq 1.1\times 10^{-2}~M_{\odot}$ and increases with the initial mass of the star.

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.

We present the serendipitous ALMA detection of a faint submillimeter galaxy (SMG) lensed by a foreground z~1 galaxy. By optimizing the source detection to deblend the system, we accurately build the full spectral energy distribution of the distant galaxy from the I814 band to radio wavelengths. It is extremely red, with a I-K colour larger than 2.5. We estimate a photometric redshift of 3.28 and determine the physical parameters. The distant galaxy turns out to be magnified by the foreground lens by a factor of ~1.5, which implies an intrinsic Ks-band magnitude of ~24.5, a submillimeter flux at 870um of ~2.5 mJy and a SFR of ~150-300Msun/yr, depending on the adopted tracer. These values place our source towards the faint end of the distribution of observed SMGs, and in particular among the still few faint SMGs with a fully characterized spectral energy distribution, which allows us not only to accurately estimate its redshift but also to measure its stellar mass and other physical properties. The galaxy studied in this work is a representative of the population of faint SMGs, of which only few objects are known to date, that are undetected in optical and therefore are not typically accounted for when measuring the cosmic star formation history (SFH). This faint galaxy population thus likely represents an important and missing piece in our understanding of the cosmic SFH. Its observation and characterization is of major importance to achieve a solid picture of galaxy evolution.

We study the impact of star-forming mini-haloes, and the Initial Mass Function (IMF) of Population III (Pop III) stars, on the Galactic halo Metallicity Distribution Function (MDF) and on the properties of C-enhanced and C-normal stars at [Fe/H]<-3. For our investigation we use a data-constrained merger tree model for the Milky Way formation, which has been improved to self-consistently describe the physical processes regulating star-formation in mini-haloes, including the poor sampling of the Pop III IMF. We find that only when star-forming mini-haloes are included the low-Fe tail of the MDF is correctly reproduced, showing a plateau that is built up by C-enhanced metal-poor (CEMP) stars imprinted by primordial faint supernovae. The incomplete sampling of the Pop III IMF in inefficiently star-forming mini-haloes (< $10^{-3}$ $M_\odot$/yr) strongly limits the formation of Pair Instability Supernovae (PISNe), with progenitor masses $m_{\rm popIII}$=[140-260] $M_\odot$, even when a flat Pop III IMF is assumed. Second-generation stars formed in environments polluted at >50% level by PISNe are thus extremely rare, corresponding to $\approx$ 0.25% of the total stellar population at [Fe/H]<-2, which is consistent with recent observations. The low-Fe tail of the MDF strongly depends on the Pop III IMF shape and mass range. Given the current statistics, we find that a flat Pop III IMF model with $m_{\rm popIII}$=[10-300] $M_\odot$ is disfavoured by observations. We present testable predictions for Pop III stars extending down to lower masses, with $m_{\rm popIII}$=[0.1-300] $M_\odot$.

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $\beta$, of galaxies at $5 \leq z \leq 8$, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess - UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-$z$ galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-$z$ galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at $z \lesssim 7$.

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.

The growth of the first super massive black holes (SMBHs) at z > 6 is still a major challenge for theoretical models. If it starts from black hole (BH) remnants of Population III stars (light seeds with mass ~ 100 Msun) it requires super-Eddington accretion. An alternative route is to start from heavy seeds formed by the direct collapse of gas onto a ~ 10^5 Msun BH. Here we investigate the relative role of light and heavy seeds as BH progenitors of the first SMBHs. We use the cosmological, data constrained semi-analytic model GAMETE/QSOdust to simulate several independent merger histories of z > 6 quasars. Using physically motivated prescriptions to form light and heavy seeds in the progenitor galaxies, we find that the formation of a few heavy seeds (between 3 and 30 in our reference model) enables the Eddington-limited growth of SMBHs at z > 6. This conclusion depends sensitively on the interplay between chemical, radiative and mechanical feedback effects, which easily erase the conditions that allow the suppression of gas cooling in the low metallicity gas (Z < Zcr and JLW > Jcr). We find that heavy seeds can not form if dust cooling triggers gas fragmentation above a critical dust-to-gas mass ratio (D > Dcr). In addition, the relative importance of light and heavy seeds depends on the adopted mass range for light seeds, as this dramatically affects the history of cold gas along the merger tree, by both SN and AGN-driven winds.

Although dust content is usually assumed to depend uniquely on metallicity, recent observations of two extremely metal-poor dwarf galaxies have suggested that this may not always be true. At a similar oxygen abundance of ~ 3% Zsun, the dust-to-gas and dust-to-stellar mass ratios in SBS 0335-052 and IZw 18 differ by a factor 40-70 according to including molecular gas or excluding it. Here we investigate a possible reason for this dramatic difference through models based on a semi-analytical formulation of chemical evolution including dust. Results suggest that the greater dust mass in SBS 0335-052 is due to the more efficient grain growth allowed by the high density in the cold interstellar medium (ISM), observationally inferred to be almost 20 times higher than in IZw 18. Our models are able to explain the difference in dust masses, suggesting that efficient dust formation and dust content in galaxies, including those with the highest measured redshifts, depend sensitively on the ISM density, rather than only on metallicity.

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 population

We investigate the origin of the FIR continuum of SDSS J1148+5251, using it as a prototype for the more general class of high-luminosity high-redshift QSOs. We run the radiative transfer code TRADING to follow the transfer of radiation from the central source and from stellar sources through the dusty environment of the host galaxy. The model is based on the output of the semi-analytical merger tree code, GAMETE/QSOdust, which enables to predict the evolution of the host galaxy and of its nuclear black hole, following the star formation history and chemical evolution -- including dust -- in all the progenitor galaxies of SDSS J1148+5251. We find that the radiation emitted by the central source can also provide an important source of heating for the dust distributed in the host galaxy, powering at least 30% and up to 70% of the observed far infrared emission at rest-frame wavelengths [20 - 1000]micron. The remaining fraction is contributed by stellar sources and can only be achieved if the host galaxy is able to sustain a star formation rate of ~ 900 Msun/yr at z=6.4. This points to a co-evolution scenario where, during their hierarchical assembly, the first SMBHs and their host galaxies first grow at the same pace until the black hole reaches a mass of ~ 2 10^8 Msun and starts growing faster than its host, reaching the bright quasar phase when the black hole and stellar mass fall within the scatter of the scaling relation observed in local galaxies. This same evolutionary scenario has been recently shown to explain the properties of a larger sample of 5 < z <6.4 QSOs, and imply that current dynamical mass measurements may have missed an important fraction of the host galaxy stellar mass. We conclude that the FIR luminosity of high-z quasars is a sensitive tracer of the rapidly changing physical conditions in the host galaxy.

We investigate the metallicity distribution function (MDF) in the Galactic halo and the relative fraction of Carbon-normal and Carbon-rich stars. To this aim, we use an improved version of the semi-analytical code GAlaxy MErger Tree and Evolution (GAMETE), that reconstructs the hierarchical merger tree of the MW, following the star formation history and the metal and dust evolution in individual progenitors. The predicted scaling relations between the dust, metal and gas masses for MW progenitors show a good agreement with observational data of local galaxies and of Gamma Ray Burst (GRB) host galaxies at 0.1 < z < 6.3. We find that in order to reproduce the observed tail of the MDF at [Fe/H] < -4, faint SN explosions have to dominate the metal yields produced by Pop III stars, disfavoring a Pop III IMF that extends to stellar masses > 140 M_{sun}, into the Pair-Instability SN progenitor mass range. The relative contribution of C-normal and C-enhanced stars to the MDF and its dependence on [Fe/H] points to a scenario where the Pop III/II transition is driven by dust-cooling and the first low-mass stars form when the dust-to-gas ratio in their parent clouds exceeds a critical value of D_crit = 4.4 x 10^{-9}.

We investigate the evolutionary properties of a sample of quasars at 5<z<6.4 using the semi-analytical hierarchical model GAMETE/QSOdust. We find that the observed properties of these quasars are well reproduced by a common formation scenario in which stars form according to a standard IMF, via quiescent star formation and efficient merger-driven bursts, while the central BH grows via gas accretion and BH-BH mergers. Eventually, a strong AGN driven wind starts to clear up the ISM of dust and gas, damping the star formation and un-obscuring the line of sight toward the QSO. In this scenario, all the QSOs hosts have final stellar masses in the range $(4-6)\times 10^{11} M_{sun}$, a factor 3-30 larger than the upper limits allowed by the observations. We discuss alternative scenarios to alleviate this apparent tension: the most likely explanation resides in the large uncertainties that still affect dynamical mass measurements in these high-z galaxies. In addition, during the transition between the starburst-dominated and the active QSO phase, we predict that about 40% of the progenitor galaxies can be classified as Sub Millimeter Galaxies, although their number rapidly decreases with redshift.

We compare theoretical dust yields for stars with mass 1 Msun < mstar < 8 Msun, and metallicities 0.001 < Z < 0.008 with observed dust production rates (DPR) by carbon- rich and oxygen-rich Asymptotic Giant Branch (C-AGB and O-AGB) stars in the Large and Small Magellanic Clouds (LMC, SMC). The measured DPR of C-AGB in the LMC are reproduced only if the mass loss from AGB stars is very efficient during the carbon-star stage. The same yields over-predict the observed DPR in the SMC, suggesting a stronger metallicity dependence of the mass-loss rates during the carbon- star stage. DPR of O-AGB stars suggest that rapid silicate dust enrichment occurs due to efficient hot-bottom-burning if mstar > 3 Msun and Z > 0.001. When compared to the most recent observations, our models support a stellar origin for the existing dust mass, if no significant destruction in the ISM occurs, with a contribution from AGB stars of 70% in the LMC and 15% in the SMC.

The dust formation process in the winds of Asymptotic Giant Branch stars is discussed, based on full evolutionary models of stars with mass in the range $1$M$_{\odot} \leq$M$\leq 8$M$_{\odot}$, and metallicities $0.001 < Z <0.008$. Dust grains are assumed to form in an isotropically expanding wind, by growth of pre--existing seed nuclei. Convection, for what concerns the treatment of convective borders and the efficiency of the schematization adopted, turns out to be the physical ingredient used to calculate the evolutionary sequences with the highest impact on the results obtained. Low--mass stars with M$\leq 3$M$_{\odot}$ produce carbon type dust with also traces of silicon carbide. The mass of solid carbon formed, fairly independently of metallicity, ranges from a few $10^{-4}$M$_{\odot}$, for stars of initial mass $1-1.5$M$_{\odot}$, to $\sim 10^{-2}$M$_{\odot}$ for M$\sim 2-2.5$M$_{\odot}$; the size of dust particles is in the range $0.1 \mu$m$\leq a_C \leq 0.2\mu$m. On the contrary, the production of silicon carbide (SiC) depends on metallicity. For $10^{-3} \leq Z \leq 8\times 10^{-3}$ the size of SiC grains varies in the range $0.05 \mu {\rm m} < {\rm a_{SiC}} < 0.1 \mu$m, while the mass of SiC formed is $10^{-5}{\rm M}_{\odot} < {\rm M_{SiC}} < 10^{-3}{\rm M}_{\odot}$. Models of higher mass experience Hot Bottom Burning, which prevents the formation of carbon stars, and favours the formation of silicates and corundum. In this case the results scale with metallicity, owing to the larger silicon and aluminium contained in higher--Z models. At Z=$8\times 10^{-3}$ we find that the most massive stars produce dust masses $m_d \sim 0.01$M$_{\odot}$, whereas models of smaller mass produce a dust mass ten times smaller. The main component of dust are silicates, although corundum is also formed, in not negligible quantities ($\sim 10-20\%$).

We present new models for the evolution of stars with mass in the range 1Msun < M < 7.5Msun, followed from the pre-main-sequence through the asymptotic giant branch phase. The metallicity adopted is $Z=3*10^{-4} (which, with an alpha-enhancement of +0.4, corresponds to [Fe/H]=-2). Dust formation is described by following the growth of dust grains of various types as the wind expands from the stellar surface. Models with mass M>3Msun experience Hot Bottom Burning, thus maintaining the surface C/O below unity. Unlike higher Z models, the scarcity of silicon available in the envelope prevents the formation of silicates in meaningful quantities, sufficient to trigger the acceleration of the wind via radiation pressure on the dust grains formed. No silicate formation occurs below a threshold metallicity of Z=10^{-3}. Low--mass stars, with M< 2.5Msun become carbon stars, forming solid carbon dust in their surroundings. The total dust mass formed depends on the uncertain extent of the inwards penetration of the convective envelope during the Third Dredge--Up episodes following the Thermal Pulses. Carbon grains have sizes 0.08 micron < a_C < 0.12 micron and the total amount of dust formed (increasing with the mass of the star) is M_C=(2-6)*10^{-4}Msun. Our results imply that AGB stars with Z=3*10^{-4} can only contribute to carbon dust enrichment of the interstellar medium on relatively long timescales, > 300 Myr, comparable to the evolutionary time of a 3Msun star. At lower metallicities the scarcity of silicon available and the presence of Hot Bottom Burning even in M< 2Msun, prevents the formation of silicate and carbon grains. We extrapolate our conclusion to more metal--poor environments, and deduce that at Z < 10^{-4} dust enrichment is mostly due to metal condensation in supernova ejecta.

Most galaxies follow well-defined scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass. However, low-metallicity starbursts, rare in the Local Universe but more common at high redshift, deviate significantly from these scaling relations. On the "main sequence" of star formation, these galaxies have high SFR for a given M*; and on the mass-metallicity relation, they have excess M* for their low metallicity. In this paper, we characterize O/H, M*, and SFR for these deviant "low-metallicity starbursts", selected from a sample of ~1100 galaxies, spanning almost two orders of magnitude in metal abundance, a factor of ~10^6 in SFR, and of ~10^5 in stellar mass. Our sample includes quiescent star-forming galaxies and blue compact dwarfs at redshift 0, luminous compact galaxies at redshift 0.3, and Lyman Break galaxies at redshifts 1-3.4. Applying a Principal Component Analysis (PCA) to the galaxies in our sample with M*<10^{10} Msun gives a Fundamental Plane (FP) of scaling relations; SFR and stellar mass define the plane itself, and O/H its thickness. The dispersion for our sample in the edge-on view of the plane is 0.17 dex, independently of redshift and including the metal-poor starbursts. The same FP is followed by 55100 galaxies selected from the Sloan Digital Sky Survey, with a dispersion of 0.06dex. In a companion paper, we develop multi-phase chemical evolution models that successfully predict the observed scaling relations and the FP; the deviations from the main scaling relations are caused by a different (starburst or "active") mode of star formation. These scaling relations do not truly evolve, but rather are defined by the different galaxy populations dominant at different cosmological epochs.

Scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass give important insight on galaxy evolution. They are obeyed by most galaxies in the Local Universe and also at high redshift. In a companion paper, we compiled a sample of ~1100 galaxies from redshift 0 to ~3, spanning almost two orders of magnitude in metal abundance, a factor of $\sim10^6$ in SFR, and of ~10^5 in stellar mass. We have characterized empirically the star-formation "main sequence" (SFMS) and the mass-metallicity relation (MZR) for this sample, and also identified a class of low-metallicity starbursts, rare locally but more common in the distant universe. These galaxies deviate significantly from the main scaling relations, with high SFR and low metal content for a given M*. In this paper, we model the scaling relations and explain these deviations from them with a set of multi-phase chemical evolution models based on the idea that, independently of redshift, initial physical conditions in a galaxy's evolutionary history can dictate its location in the scaling relations. Our models are able to successfully reproduce the O/H, M*, and SFR scaling relations up to z~3, and also successfully predict the molecular cloud fraction as a function of stellar mass. These results suggest that the scaling relations are defined by different modes of star formation: an "active" starburst mode, more common at high redshift, and a quiescent "passive" mode that is predominant locally and governs the main trends.

Galaxy-scale gas outflows triggered by active galactic nuclei have been proposed as a key physical process to regulate the co-evolution of nuclear black holes and their host galaxies. The recent detection of a massive gas outflow in one of the most distant quasar, SDSS J1148+5251 at z = 6.4, presented by Maiolino et al. (2012) strongly supports this idea and suggests that strong quasar feedback is already at work at very early times. In a previous work, Valiante et al. (2011), we have presented a hierarchical semi-analytical model, GAMETE/ QSOdust, for the formation and evolution of high-redshift quasars, and we have applied it to the quasar SDSS J1148+5251, with the aim of investigating the star formation history, the nature of the dominant stellar populations and the origin and properties of the large dust mass observed in the host galaxy. A robust prediction of the model is that the evolution of the nuclear black hole and of the host galaxy are tightly coupled by quasar feedback in the form of strong galaxy-scale winds. In the present letter, we show that the gas outflow rate predicted by GAMETE/QSOdust is in good agreement with the lower limit of 3500 Msun/yr inferred by the observations. According to the model, the observed outflow at z = 6.4 is dominated by quasar feedback, as the outflow rate has already considerably depleted the gas content of the host galaxy, leading to a down-turn in the star formation rate at z < 7 - 8. Hence, supernova explosions give a negligible contribution to the observed winds at z = 6.4.

We calculate the dust formed around AGB and SAGB stars of metallicity Z=0.008 by following the evolution of models with masses in the range 1M<M<8M throughthe thermal pulses phase, and assuming that dust forms via condensation of molecules within a wind expanding isotropically from the stellar surface. We find that, because of the strong Hot Bottom Burning (HBB) experienced, high mass models produce silicates, whereas lower mass objects are predicted to be surrounded by carbonaceous grains; the transition between the two regimes occurs at a threshold mass of 3.5M. These fndings are consistent with the results presented in a previous investigation, for Z=0.001. However, in the present higher metallicity case, the production of silicates in the more massive stars continues for the whole AGB phase, because the HBB experienced is softer at Z=0.008 than at Z=0.001, thus the oxygen in the envelope, essential for the formation of water molecules, is never consumed completely. The total amount of dust formed for a given mass experiencing HBB increases with metallicity, because of the higher abundance of silicon, and the softer HBB, both factors favouring a higher rate of silicates production. This behaviour is not found in low mass stars,because the carbon enrichment of the stellar surface layers, due to repeated Third Drege Up episodes, is almost independent of the metallicity. Regarding cosmic dust enrichment by intermediate mass stars, we find that the cosmic yield at Z=0.008 is a factor 5 larger than at Z=0.001. In the lower metallicity case carbon dust dominates after about 300 Myr, but at Z=0.008 the dust mass is dominated by silicates at all times,with a prompt enrichment occurring after about 40 Myr, associated with the evolution of stars with masses M =7.5 -8M.

We compute the mass and composition of dust produced by stars with masses in the range 1Msun<M<8 Msun and with a metallicity of Z=0.001 during their AGB and Super AGB phases. Stellar evolution is followed from the pre-main sequence phase using the code ATON which provides, at each timestep, the thermodynamics and the chemical stucture of the wind. We use a simple model to describe the growth of the dust grains under the hypothesis of a time-independent, spherically symmetric stellar wind. We find that the total mass of dust injected by AGB stars in the interstellar medium does not increase monotonically with stellar mass and ranges between a minimum of 10^{-6}Msun for the 1.5Msun stellar model, up to 2x10^{-4} Msun, for the 6Msun case. Dust composition depends on the stellar mass: low-mass stars (M < 3Msun) produce carbon-rich dust, whereas more massive stars, experiencing Hot Bottom Burning, never reach the carbon-star stage, and produce silicates and iron. This is in partial disagreement with previous investigations in the literature, which are based on synthetic AGB models and predict that, when the initial metallicity is Z=0.001, C-rich dust is formed at all stellar masses. The differences are due to the different modelling of turbulent convection in the super-adiabaticity regime. Also in this case the treatment of super-adiabatic convection shows up as the most relevant issue affecting the dust-formation process. We also investigate Super AGB stars with masses 6.5Msun<M<8 Msun that evolve over a ONe core.Due to a favourable combination of mass loss and Hot Bottom Burning, these stars are predicted to be the most efficient silicate-dust producers, releasing [2 - 7]x 10^{-4} Msun masses of dust. We discuss the robustness of these predictions and their relevance for the nature and evolution of dust at early cosmic times.

We explore the minimal conditions which enable the formation of metal-enriched solar and sub-solar mass stars. We find that in the absence of dust grains, gas fragmentation occurs at densities nH ~ [10^4-10^5]cm^{-3} when the metallicity exceeds Z ~ 10^{-4} Zsun. The resulting fragmentation masses are > 10 Msun. The inclusion of Fe and Si cooling does not affect the thermal evolution as this is dominated by molecular cooling even for metallicities as large as Z = 10^{-2} Zsun. The presence of dust is the key driver for the formation of low-mass stars. We focus on three representative core-collapse supernova (SN) progenitors, and consider the effects of reverse shocks of increasing strength: these reduce the depletion factors, fdep = Mdust/(Mdust+Mmet), alter the shape of the grain size distribution function and modify the relative abundances of grain species and of metal species in the gas phase. We find that the lowest metallicity at which fragmentation occurs is Z=10^{-6} Zsun for gas pre-enriched by the explosion of a 20 Msun primordial SN (fdep > 0.22) and/or by a 35 Msun, Z=10^{-4} Zsun SN (fdep > 0.26); it is ~ 1 dex larger, when the gas is pre-enriched by a Z = 10^{-4} Zsun, 20 Msun SN (fdep > 0.04). Cloud fragmentation depends on the depletion factor and it is suppressed when the reverse shock leads to a too large destruction of dust grains. These features are all consistent with the existence of a minimum dust-to-gas ratio, Dcr, above which fragmentation is activated. We derive a simple analytic expression for Dcr which, for grain composition and properties explored in the present study, reads Dcr = [2.6 - 6.3] x 10^{-9}. When the dust-to-gas ratio of star forming clouds exceeds this value, the fragmentation masses range between 0.01 Msun and 1 Msun, thus enabling the formation of the first low-mass stars.

We present a semi-analytical model for the formation and evolution of a high redshift quasar (QSO). We reconstruct a set of hierarchical merger histories of a 10^13 Msun dark matter halo and model the evolution of the corresponding galaxy and of its central super massive black hole. The code GAMETE/QSODUST consistently follows: the black hole assembly via both coalescence with other black holes and gas accretion; the build up and star formation history of the quasar host galaxy, driven by binary mergers and mass accretion; the evolution of gas, stars, metals in the interstellar medium (ISM), accounting for mass exchanges with the external medium (infall and outflow processes); dust formation in Supernova (SN) ejecta and in the stellar atmosphere of Asymptotic Giant Branch (AGB) stars, dust destruction by interstellar shocks and grain growth in molecular clouds; the AGN feedback which powers a galactic-scale wind, self-regulating the black hole growth and eventually halting star formation.We use this model to study the case of SDSS J1148+5251 observed at redshift 6.4, exploring different star formation histories for the QSO host galaxy.

With the aim of investigating whether stellar sources can account for the >10^8 Msun dust masses inferred from mm/sub-mm observations of samples of 5<z<6.4 quasars,we develop a chemical evolution model which follows the evolution of metals and dust on the stellar characteristic lifetimes, taking into account dust destruction mechanisms.Using a grid of stellar dust yields as a function of the initial mass and metallicity over the range 1-40 Msun and 0-1 Zsun,we show that the role of AGB stars in cosmic dust evolution at high redshift might have been over-looked.We apply the chemical evolution model with dust to the host galaxy of the most distant quasar at z=6.4, SDSS J1148+5251.Given the current uncertainties on the star formation history of the host galaxy, we have considered two models: (i) a star formation history obtained in a numerical simulation by Li et al.(2007) which predicts that a large stellar bulge is already formed at z=6.4,and (ii) a constant star formation rate of 1000 Msun/yr, as suggested by the observations if most of the FIR luminosity is due to young stars.The total mass of dust predicted at z=6.4 by the first model is 2x10^8Msun,within the range of values inferred by observations,with a substantial contribution (80%) of AGB-dust.When a constant star formation rate is adopted,the contribution of AGB-dust decreases to 50% but the total mass of dust formed is a factor 2 smaller.Both models predict a rapid enrichment of the ISM with metals and a relatively mild evolution of the carbon abundance,in agreement with observational constraints. This supports the idea that stellar sources can account for the dust observed but show that the contribution of AGB stars to dust production cannot be neglected, even at the most extreme redshifts currently accessible to observations.

Our aim is to show how different hypotheses about Type Ia supernova progenitors can affect Galactic chemical evolution. We include different Type Ia SN progenitor models, identified by their distribution of time delays, in a very detailed chemical evolution model for the Milky Way which follows the evolution of several chemical species. We test the single degenerate and the double degenerate models for supernova Ia progenitors, as well as other more empirical models based on differences in the time delay distributions. We find that assuming the single degenerate or the double degenerate scenario produces negligible differences in the predicted [O/Fe] vs. [Fe/H] relation. On the other hand, assuming a percentage of prompt (exploding in the first 100 Myr) Type Ia supernovae of 50%, or that the maximum Type Ia rate is reached after 3-4 Gyr from the beginning of star formation, as suggested by several authors, produces more noticeable effects on the [O/Fe] trend. However, given the spread still existing in the observational data no model can be firmly excluded on the basis of only the [O/Fe] ratios. On the other hand, when the predictions of the different models are compared with the G-dwarf metallicity distribution, the scenarios with very few prompt Type Ia supernovae can be excluded. Models including the single degenerate or double degenerate scenario with a percentage of 10-13% of prompt Type Ia supernovae produce results in very good agreement with the observations. A fraction of prompt Type Ia supernovae larger than 30% worsens the agreement with observations and the same occurs if no prompt Type Ia supernovae are allowed. In particular, two empirical models for the Type Ia SN progenitors can be excluded: the one without prompt Type Ia supernovae and the one assuming delay time distribution going like t^{-0.5}.

The aim of this work is the computation of the cosmic Type Ia supernova rates at very high redshifts (z>2). We adopt various progenitor models in order to predict the number of explosions in different scenarios for galaxy formation and to check whether it is possible to select the best delay time distribution model, on the basis of the available observations of Type Ia supernovae. We also computed the Type Ia supernova rate in typical elliptical galaxies of different initial luminous masses and the total amount of iron produced by Type Ia supernovae in each case. It emerges that: it is not easy to select the best delay time distribution scenario from the observational data and this is because the cosmic star formation rate dominates over the distribution function of the delay times; the monolithic collapse scenario predicts an increasing trend of the SN Ia rate at high redshifts whereas the predicted rate in the hierarchical scheme drops dramatically at high redshift; for the elliptical galaxies we note that the predicted maximum of the Type Ia supernova rate depends on the initial galactic mass. The maximum occurs earlier (at about 0.3 Gyr) in the most massive ellipticals, as a consequence of downsizing in star formation. We find that different delay time distributions predict different relations between the Type Ia supernova rate per unit mass at the present time and the color of the parent galaxies and that bluer ellipticals present higher supernova Type Ia rates at the present time.