We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored.

We aim to provide the first data-constrained estimate of the metal-free (Population III; Pop III) star formation rate density $\dot{\rho}_{*}^{III}$ required at high-redshifts ($z \gtrsim 16$) in order to reproduce both the amplitude and the redshift of the EDGES 21-cm global signal. Our model accounts for the Lyman Alpha (Ly$\alpha$), radio and X-ray backgrounds from both Pop III and metal-enriched Population II (Pop II) stars. For the latter, we use the star formation rate density estimates (and the Ly$\alpha$ background) from the {\it Delphi} semi-analytic model that has been shown to reproduce all key observables for galaxies at $z \gtrsim 5$; the radio and X-ray backgrounds are fixed using low-$z$ values. The constraints on the free parameters characterizing the properties of the Pop III stars are obtained using a Markov Chain Monte Carlo analysis. Our results yield a $\dot{\rho}_{*}^{III}$ that whilst increasing from $z \sim 21-16$ thereafter shows a sharp decline which is in excellent agreement with the results found by \citet{valiante2016} to simulate the growth of $z \sim 6 - 7$ quasars and their host galaxies, suggesting that the bulk of Pop III star formation occurs in the rarest and most massive metal-poor halos at $z \lesssim 20$. This allows Pop III stars to produce a rapidly growing Ly$\alpha$ background between $z \sim 21-15$. Further, Pop III stars are required to provide a radio background that is about $3-4$ orders of magnitude higher than that provided by Pop II stars although Pop II stars dominate the X-ray background.

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

Metallicity and gas content are intimately related in the baryonic exchange cycle of galaxies, and galaxy evolution scenarios can be constrained by quantifying this relation. To this end, we have compiled a sample of ~400 galaxies in the Local Universe, dubbed "MAGMA" (Metallicity And Gas for Mass Assembly), which covers an unprecedented range in parameter space, spanning more than 5 orders of magnitude in stellar mass (Mstar), star-formation rate (SFR), and gas mass (Mgas), and a factor of ~60 in metallicity [Z, 12+log(O/H)]. Stellar masses and SFRs have been recalculated for all the galaxies using IRAC, WISE and GALEX photometry, and 12+log(O/H) has been transformed, where necessary, to a common metallicity calibration. To assess the true dimensionality of the data, we have applied multi-dimensional principal component analyses (PCAs) to our sample. In confirmation of previous work, we find that even with the vast parameter space covered by MAGMA, the relations between Mstar, SFR, Z and Mgas (MHI+MH2) require only two dimensions to describe the hypersurface. To accommodate the curvature in the Mstar-Z relation, we have applied a piecewise 3D PCA that successfully predicts observed 12+log(O/H) to an accuracy of ~0.1dex. MAGMA is a representative sample of isolated star-forming galaxies in the Local Universe, and can be used as a benchmark for cosmological simulations and to calibrate evolutionary trends with redshift.

The existence of massive stellar black hole binaries (MBHBs), with primary black hole (BH) masses $\ge 31 \, M_\odot$, was proven by the detection of the gravitational wave (GW) event GW150914 during the first LIGO/Virgo observing run (O1), and successively confirmed by seven additional GW signals discovered in the O1 and O2 data. By adopting the galaxy formation model \texttt{GAMESH} coupled with binary population synthesis (BPS) calculations, here we investigate the origin of these MBHBs by selecting simulated binaries compatible in mass and coalescence redshifts. We find that their cosmic birth rates peak in the redshift range $6.5 \leq z \leq 10$, regardless of the adopted BPS. These MBHBs are then old systems forming in low-metallicity ($Z \sim [0.01-0.1] \, Z_{\odot}$), low-stellar-mass galaxies, before the end of cosmic reionization, i.e. significantly beyond the peak of cosmic star formation. GW signals generated by coalescing MBHBs open up new possibilities to probe the nature of stellar populations in remote galaxies, at present too faint to be detected by available electromagnetic facilities.

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $\lambda 1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

The recent discovery of high redshift dusty galaxies implies a rapid dust enrichment of their interstellar medium (ISM). To interpret these observations, we run a cosmological simulation in a 30$h^{-1}$ cMpc/size volume down to $z \approx 4$. We use the hydrodynamical code dustyGadget, which accounts for the production of dust by stellar populations and its evolution in the ISM. We find that the cosmic dust density parameter ($\Omega_{\rm d}$) is mainly driven by stellar dust at $z \gtrsim 10$, so that mass- and metallicity-dependent yields are required to assess the dust content in the first galaxies. At $z \lesssim 9$ the growth of grains in the ISM of evolved systems (Log$(M_{\star}/M_{\odot})>8.5$) significantly increases their dust mass, in agreement with observations in the redshift range $4 \lesssim z < 8$. Our simulation shows that the variety of high redshift galaxies observed with ALMA can naturally be accounted for by modeling the grain-growth timescale as a function of the physical conditions in the gas cold phase. In addition, the trends of dust-to-metal (DTM) and dust-to-gas (${\cal D}$) ratios are compatible with the available data. A qualitative investigation of the inhomogeneous dust distribution in a representative massive halo at $z \approx 4$ shows that dust is found from the central galaxy up to the closest satellites along polluted filaments with $\rm Log({\cal D}) \leq -2.4$, but sharply declines at distances $d \gtrsim 30$ kpc along many lines of sight, where $\rm Log({\cal D}) \lesssim -4.0$.

The AGN bolometric correction is a key element to understand BH demographics and compute accurate BH accretion histories from AGN luminosities. However, current estimates still differ from each other by up to a factor of two to three, and rely on extrapolations at the lowest and highest luminosities. Here we revisit this fundamental issue presenting general hard X-ray ($K_{X}$) and optical ($K_{O}$) bolometric corrections, computed combining several AGN samples spanning the widest (about 7 dex) luminosity range ever used for this kind of studies. We analysed a total of $\sim 1000$ type 1 and type 2 AGN for which a dedicated SED-fitting has been carried out. We provide a bolometric correction separately for type 1 and type 2 AGN; the two bolometric corrections results to be in agreement in the overlapping luminosity range and therefore, for the first time, a universal bolometric correction for the whole AGN sample (both type 1 and type 2) has been computed. We found that $K_{X}$ is fairly constant at $log(L_{BOL}/L_{\odot}) < 11$, while it increases up to about one order of magnitude at $log(L_{BOL}/L_{\odot}) \sim 14.5$. A similar increasing trend has been observed when its dependence on either the Eddington ratio or the BH mass is considered, while no dependence on redshift up to $z\sim3.5$ has been found. On the contrary, the optical bolometric correction appears to be fairly constant (i.e. $K_{O} \sim 5$) whatever is the independent variable. We also verified that our bolometric corrections correctly predict the AGN bolometric luminosity functions. According to this analysis, our bolometric corrections can be applied to the whole AGN population in a wide range of luminosity and redshift.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

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 cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Here we present a detailed analysis of the properties and evolution of different dwarf galaxies, candidate to host the coalescence of black hole binary systems (BHB) generating GW150914-like events. By adopting a novel theoretical framework coupling the binary population synthesis code \texttt{SeBa} with the Galaxy formation model \texttt{GAMESH}, we can investigate the detailed evolution of these objects in a well resolved cosmological volume of 4~cMpc, having a Milky Way-like (MW) galaxy forming at its center. We identify three classes of interesting candidate galaxies: MW progenitors, dwarf satellites and dwarf galaxies evolving in isolation. We find that: (i) despite differences in individual histories and specific environments the candidates reduce to only nine representative galaxies; (ii) among them, $\sim44\%$ merges into the MW halo progenitors by the redshift of the expected signal, while the remaining dwarfs are found as isolated or as satellites of the MW and their evolution is strongly shaped by both peculiar dynamical history and environmental feedback; (iii) a stringent condition for the environments where GW150914-like binaries can form comes from a combination of the accretion history of their DM halos and the radiative feedback in the high redshift universe; (iv) by comparing with the observed catalogues from DGS and ALLSMOG surveys we find two observed dwarfs respecting the properties predicted by our model. We finally note how the present analysis opens the possibility to build future strategies for host galaxy identification.

Supernovae (SNe) are considered to have a major role in dust enrichment of high redshift galaxies and, due to the short lifetimes of interstellar grains, in dust replenishment of local galaxies. Here we explore how SN dust yields depend on the mass, metallicity, and rotation rate of the progenitor stars, and on the properties of the explosion. To this aim, assuming uniform mixing inside the ejecta, we quantify the dust mass produced by a sample of SN models with progenitor masses $13~M_{\odot} \leq M \leq 120~M_{\odot}$, metallicity $\rm -3 \leq [Fe/H] \leq 0$, rotation rate $\rm v_{\rm rot} = 0$ and $300$~km/s, that explode with a fixed energy of $1.2 \times 10^{51}$~erg (FE models) or with explosion properties calibrated to reproduce the $\rm ^{56}Ni$ - $M$ relation inferred from SN observations (CE models). We find that rotation favours more efficient dust production, particularly for more massive, low metallicity stars, but that metallicity and explosion properties have the largest effects on the dust mass and its composition. In FE models, SNe with $M \leq 20 - 25 ~M_{\odot}$ are more efficient at forming dust: between 0.1 and 1 $M_\odot$ is formed in a single explosion, with a composition dominated by silicates, carbon and magnetite grains when $\rm [Fe/H] = 0$, and by carbon and magnetite grains when $\rm [Fe/H] < 0$. In CE models, the ejecta are massive and metal-rich and dust production is more efficient. The dust mass increases with $M$ and it is dominated by silicates, at all [Fe/H].

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.

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI.

We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $\Lambda$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.

Shocks may have been prevalent in the early Universe, associated with virialization and supernova explosions, etc. Here, we study thermal evolution and fragmentation of shock-compressed clouds, by using a one-zone model with detailed thermal and chemical processes. We explore a large range of initial density (1-1e5 /cm^3), metallicity (0-1e-2 Z_sun), UV strength (0-500 times Galactic value), and cosmic microwave background temperature (10 and 30 K). Shock-compressed clouds contract isobarically via atomic and molecular line cooling, until self-gravitating clumps are formed by fragmentation. If the metals are only in the gas-phase, the clump mass is higher than ~ 3 M_sun in any conditions we studied. Although in some cases with a metallicity higher than ~ 1e-3 Z_sun, re-fragmentation of a clump is caused by metal-line cooling, this fragment mass is higher than ~ 30 M_sun. On the other hand, if about half the mass of metals is condensed in dust grains, as in the Galactic interstellar medium, dust cooling triggers re-fragmentation of a clump into sub-solar mass pieces, for metallicities higher than ~ 1e-5 Z_sun. Therefore, the presence of dust is essential in low-mass (< M_sun) star formation from a shock-compressed cloud.

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.

In this work we report deep MUSE observations of a Broad Absorption Line (BAL) quasar at z ~ 5, revealing a Lyman alpha nebula with a maximum projected linear size of ~ 60 kpc around the quasar (down to our 2-sigma SB limit per layer of ~ 9e-19 erg/s/cm^2/arcsec^2 for a 1 arcsec^2 aperture). After correcting for the cosmological surface brightness dimming, we find that our nebula, at z ~ 5, has an intrinsically less extended Lyman alpha emission than nebulae at lower redshift. However, such a discrepancy is greatly reduced when referring to comoving distances, which take into account the cosmological growth of dark matter (DM) haloes, suggesting a positive correlation between the size of Lyman alpha nebulae and the sizes of DM haloes/structures around quasars. Differently from the typical nebulae around radio-quiet non-BAL quasars, in the inner regions (~ 10 kpc) of the circumgalactic medium (CGM) of our source, the velocity dispersion of the Lyman alpha emission is very high (FWHM > 1000 km/s), suggesting that in our case we may be probing outflowing material associated with the quasar.

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.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

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.

We present a novel theoretical model to characterize the formation and coalescence sites of compact binaries in a cosmological context. This is based on the coupling between the binary population synthesis code SeBa with a simulation following the formation of a Milky Way-like halo in a well resolved cosmic volume of 4 cMpc, performed with the GAMESH pipeline. We have applied this technique to investigate when and where systems with properties similar to the recently observed LIGO/VIRGO events are more likely to form and where they are more likely to reside when they coalesce. We find that more than 70% of GW151226 and LVT151012-like systems form in galaxies with stellar mass M* > 10^8 Msun in the redshift range [0.06 - 3] and [0.14 - 11.3], respectively. All GW150914-like systems form in low-metallicity dwarfs with M* < 5 10^6 Msun at 2.4 < z < 4.2. Despite these initial differences, by the time they reach coalescence the observed events are most likely hosted by star forming galaxies with M* > 10^{10} Msun. Due to tidal stripping and radiative feedback, a non negligible fraction of GW150914-like candidates end-up in galaxies with properties similar to dwarf spheroidals and ultra-faint satellites.

We present new ALMA observations aimed at mapping molecular gas reservoirs through the CO(3-2) transition in three quasars at $z\simeq2.4$, LBQS 0109+0213, 2QZ J002830.4-281706, and [HB89] 0329-385. Previous [OIII]5007 observations of these quasars showed evidence for ionised outflows quenching star formation in their host galaxies. Systemic CO(3-2) emission has been detected only in one quasar, LBQS 0109+0213, where the CO(3-2) emission is spatially anti-correlated with the ionised outflow, suggesting that most of the molecular gas may have been dispersed or heated in the region swept by the outflow. In all three sources, including the one detected in CO, our constraints on the molecular gas mass indicate a significantly reduced reservoir compared to main-sequence galaxies at the same redshift, supporting a negative feedback scenario. In the quasar 2QZ J002830.4-281706, we tentatively detect an emission line blob blue-shifted by $v\sim-2000$ km/s with respect to the galaxy systemic velocity and spatially offset by 0.2 arcsec (1.7 kpc) with respect to the ALMA continuum peak. Interestingly, such emission feature is coincident in both velocity and space with the ionised outflow as seen in [OIII]5007. This tentative detection must be confirmed with deeper observations but, if real, it could represent the molecular counterpart of the ionised gas outflow driven by the AGN. Finally, in all ALMA maps we detect the presence of serendipitous line emitters within a projected distance $\sim 160$ kpc from the quasars. By identifying these features with the CO(3-2) transition, the serendipitous line emitters would be located within |$\Delta v$|$<$500 km/s from the quasars, hence suggesting an overdensity of galaxies in two out of three quasars.

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.

Type 2 active galactic nuclei (AGN) represent the majority of the AGN population. However, due to the difficulties in measuring their black hole (BH) masses, it is still unknown whether they follow the same BH mass-host galaxy scaling relations valid for quiescent galaxies and type 1 AGN. Here we present the locus of type 2 AGN having virial BH mass estimates in the $M_{BH}-\sigma_\star$ plane. Our analysis shows that the BH masses of type 2 AGN are $\sim0.9$ dex smaller than type 1 AGN at $\sigma_\star\sim 185$ km s$^{-1}$, regardless of the (early/late) AGN host galaxy morphology. Equivalently, type 2 AGN host galaxies have stellar velocity dispersions $\sim 0.2$ dex higher than type 1 AGN hosts at $M_{BH}\sim10^7$ M$_\odot$.

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.

Here we investigate the evolution of a Milky Way (MW) -like galaxy with the aim of predicting the properties of its progenitors all the way from $z \sim 20$ to $z = 0$. We apply GAMESH (Graziani et al. 2015) to a high resolution N-Body simulation following the formation of a MW-type halo and we investigate its properties at $z \sim 0$ and its progenitors in $0 < z < 4$. Our model predicts the observed galaxy main sequence, the mass-metallicity and the fundamental plane of metallicity relations in $0 < z < 4$. It also reproduces the stellar mass evolution of candidate MW progenitors in $0 \lesssim z \lesssim 2.5$, although the star formation rate and gas fraction of the simulated galaxies follow a shallower redshift dependence. We find that while the MW star formation and chemical enrichment are dominated by the contribution of galaxies hosted in Lyman $\alpha$-cooling halos, at z > 6 the contribution of star forming mini-halos is comparable to the star formation rate along the MW merger tree. These systems might then provide an important contribution in the early phases of reionization. A large number of mini-halos with old stellar populations, possibly Population~III stars, are dragged into the MW or survive in the Local Group. At low redshift dynamical effects, such as halo mergers, tidal stripping and halo disruption redistribute the baryonic properties among halo families. These results are critically discussed in light of future improvements including a more sophisticated treatment of radiative feedback and inhomogeneous metal enrichment.

We report ALMA observations of the most massive (star forming) galaxy in the redshift range 3<z<4 within the whole GOODS-S field. We detect a large elongated structure of molecular gas around the massive primeval galaxy, traced by the CO(4-3) emission, and extended over 40 kpc. We infer a mass of the large gaseous structure of Mgas~2-6x10^11 Msun. About 60% of this mass is not directly associated with either the central galaxy or its two lower mass satellites. The CO extended structure is also detected in continuum thermal emission. The kinematics of the molecular gas shows the presence of different components, which cannot be ascribed to simple rotation. Furthermore, on even larger scales, we detect nine additional CO systems within a radius of 250 kpc from the massive galaxy and mostly distributed in the same direction as the CO elongated structure found in the central 40 kpc. The stacked images of these CO systems show detections in the thermal continuum and in the X-rays, suggesting that these systems are forming stars at a rate of 30-120 Msun/yr. We suggest that the extended gas structure, combined with its kinematic properties, and the gas rich star forming systems detected on larger scales, are tracing the inner and densest regions of large scale accreting streams, feeding the central massive galaxy. These results corroborate models of galaxy formation, in which accreting streams are clumpy and undergo some star formation (hence enriching the streams with metals) even before accreting onto the central galaxy.

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.

Models and observations suggest that both power and effects of AGN feedback should be maximised in hyper-luminous (L_Bol>10^47 erg/s) quasars, i.e. objects at the brightest end of the AGN luminosity function. We present the first results of a multi-wavelength observing program, focusing on a sample of WISE/SDSS selected hyper-luminous (WISSH) broad-line quasars at z~1.5-5. The WISSH quasars project has been designed to reveal the most energetic AGN-driven outflows, estimate their occurrence at the peak of quasar activity and extend the study of correlations between outflows and nuclear properties up to poorly-investigated extreme AGN luminosities (L_Bol~10^47 -10^48 erg/s). We present NIR, long-slit LBT/LUCI1 spectroscopy of five WISSH quasars at z~2.3-3.5 showing prominent [OIII] emission lines with broad (FWHM~1200-2200 km/s) and skewed profiles. The luminosities of the broad [OIII] wings are the highest measured so far (L_[OIII]^broad >~5x10^44 erg/s) and reveal the presence of powerful ionised outflows with mass outflow rates Mdot >~1700 M_Sun/yr and kinetic powers Edot >~10^45 erg/s. Although these estimates are affected by large uncertainties, due to the use of [OIII] as tracer of ionized outflows and the very basic outflow model we assume, these results suggest that the AGN is highly efficient in pushing outwards large amounts of ionised gas in our targets. The mechanical outflow luminosities for WISSH quasars correspond to higher fractions (~1-3%) of L_Bol than those derived for lower L_Bol AGN. Our targets host very massive (M_BH>~2x10^9 M_Sun) black holes which are still accreting at a high rate (i.e. a factor of ~0.4-3 of the Eddington limit). These findings demonstrate that WISSH quasars offer the opportunity of probing the extreme end of both luminosity and SMBH mass functions and revealing powerful ionised outflows able to affect the evolution of their host galaxies.

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$.

Dust grains are classically thought to form in the winds of asymptotic giant branch (AGB) stars. However, there is increasing evidence today for dust formation in supernovae (SNe). To establish the relative importance of these two classes of stellar sources of dust, it is important to know the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. We have developed a new code (GRASH\_Rev) which follows the newly-formed dust evolution throughout the supernova explosion until the merging of the forward shock with the circumstellar ISM. We have considered four well studied SNe in the Milky Way and Large Magellanic Cloud: SN1987A, CasA, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations and estimate that between 1 and 8$\%$ of the observed mass will survive, leading to a SN dust production rate of $(3.9 \pm 3.7) \times 10^{-4}$ M$_{\odot}$yr$^{-1}$ in the Milky Way. This value is one order of magnitude larger than the dust production rate by AGB stars but insufficient to counterbalance the dust destruction by SNe, therefore requiring dust accretion in the gas phase.

The direct discovery of gravitational waves has finally opened a new observational window on our Universe, suggesting that the population of coalescing binary black holes is larger than previously expected. These sources produce an unresolved background of gravitational waves, potentially observables by ground-based interferometers. In this paper we investigate how modified theories of gravity, modeled using the ppE formalism, affect the expected signal, and analyze the detectability of the resulting stochastic background by current and future ground-based interferometers. We find the constraints that AdLIGO would be able to set on modified theories, showing that they may significantly improve the current bounds obtained from astrophysical observations of binary pulsars.

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$.

Negative feedback from active galactic nuclei (AGN) is considered a key mechanism in shaping galaxy evolution. Fast, extended outflows are frequently detected in the AGN host galaxies at all redshifts and luminosities, both in ionised and molecular gas. However, these outflows are only "potentially" able to quench star formation and we are still missing a decisive evidence of negative feedback in action. Here we present Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) H- and K-band integral-field spectroscopic observations of two quasars at $z\sim$2.4 characterised by fast, extended outflows detected through the [OIII]$\lambda$5007 line (Carniani et al. 2015). The high signal-to-noise ratio of our observations allows us to identify faint narrow (FWHM $< 500$ km/s), and spatially extended components in [OIII]$\lambda$5007 and H$\alpha$ emission associated with star formation in the host galaxy. Such star-formation powered emission is spatially anti-correlated with the fast outflows. The ionised outflows therefore appear to be able to suppress star formation in the region where the outflow is expanding. However the detection of narrow, spatially-extended H$\alpha$ emission indicates star formation rates of at least $\sim 50-90$ M$_{\odot}$/yr, suggesting either that AGN feedback does not affect the whole galaxy or that many feedback episodes are required before star formation is completely quenched. On the other hand, the narrow H$\alpha$ emission extending along the edges of the outflow cone may lead also to a positive feedback interpretation. Our results highlight the possible twofold role of galaxy-wide outflows in host galaxy evolution.

Dust grains are classically thought to form in the winds of AGB stars. However, nowadays there is increasing evidence for dust formation in SNe. In order to establish the relative importance of these two classes of stellar sources of dust it is important to know what is the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows to follow the dynamics of dust grains in the shocked SN ejecta and to compute the time evolution of the mass, composition and size distribution of the grains. We consider four well studied SNe in the Milky Way and LMC: SN 1987a, Cas A, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Conversely, in the other three SNe, the reverse shock has already destroyed between 10 and 40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 10^3 and 10^5 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1 and 8% of the currently observed mass will survive. This is in good agreement with the values adopted in chemical evolution models which consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.

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 study the population of asymptotic giant branch (AGB) stars in the Small Magellanic Cloud (SMC) by means of full evolutionary models of stars of mass 1Msun < M < 8Msun, evolved through the thermally pulsing phase. The models also account for dust production in the circumstellar envelope. We compare Spitzer infrared colours with results from theoretical modelling. We show that ~75% of the AGB population of the SMC is composed by scarcely obscured objects, mainly stars of mass M < 2.5Msun at various metallicity, formed between 700 Myr and 5 Gyr ago; ~ 70% of these sources are oxygen--rich stars, while ~ 30% are C-stars. The sample of the most obscured AGB stars, accounting for ~ 25% of the total sample, is composed almost entirely by carbon stars. The distribution in the colour-colour ([3.6]-[4.5], [5.8]-[8.0]) and colour-magnitude ([3.6]-[8.0], [8.0]) diagrams of these C-rich objects, with a large infrared emission, traces an obscuration sequence, according to the amount of carbonaceous dust in their surroundings. The overall population of C-rich AGB stars descends from 1.5-2Msun stars of metallicity Z=0.004, formed between 700 Myr and 2 Gyr ago, and from lower metallicity objects, of mass below 1.5Msun, 2-5 Gyr old. We also identify obscured oxygen-rich stars (M ~ 4-6Msun) experiencing hot bottom burning. The differences between the AGB populations of the SMC and LMC are also commented.

We quantify the role of Population (Pop) III core-collapse supernovae (SNe) as the first cosmic dust polluters. Starting from a homogeneous set of stellar progenitors with masses in the range [13 - 80] Msun, we find that the mass and composition of newly formed dust depend on the mixing efficiency of the ejecta and the degree of fallback experienced during the explosion. For standard Pop III SNe, whose explosions are calibrated to reproduce the average elemental abundances of Galactic halo stars with [Fe/H] < -2.5, between 0.18 and 3.1 Msun (0.39 - 1.76 Msun) of dust can form in uniformly mixed (unmixed) ejecta, and the dominant grain species are silicates. We also investigate dust formation in the ejecta of faint Pop III SN, where the ejecta experience a strong fallback. By examining a set of models, tailored to minimize the scatter with the abundances of carbon-enhanced Galactic halo stars with [Fe/H ] < -4, we find that amorphous carbon is the only grain species that forms, with masses in the range 2.7 10^{-3} - 0.27 Msun (7.5 10^{-4} - 0.11 Msun) for uniformly mixed (unmixed) ejecta models. Finally, for all the models we estimate the amount and composition of dust that survives the passage of the reverse shock, and find that, depending on circumstellar medium densities, between 3 and 50% (10 - 80%) of dust produced by standard (faint) Pop III SNe can contribute to early dust enrichment.

Observations of edge-on galaxies allow us to investigate the vertical extent and properties of dust, gas and stellar distributions. NGC 891 has been studied for decades and represents one of the best studied cases of an edge-on galaxy. We use deep PACS data together with IRAC, MIPS and SPIRE data to study the vertical extent of dust emission around NGC 891. We also test the presence of a more extended, thick dust component. By performing a convolution of an intrinsic vertical profile emission with each instrument PSF and comparing it with observations we derived the scaleheight of a thin and thick dust disc component. For all wavelengths considered the emission is best fit with the sum of a thin and a thick dust component. The scaleheight of both dust components shows a gradient passing from 70 $\mu$m to 250 $\mu$m. This could be due to a drop in dust heating (and thus dust temperature) with the distance from the plane, or to a sizable contribution ($\sim 15 - 80%$) of an unresolved thin disc of hotter dust to the observed surface brightness at shorter wavelengths. The scaleheight of the thick dust component, using observations from 70 $\mu$m to 250 $\mu$m has been estimated to be $(1.44\pm 0.12)$ kpc, consistent with previous estimates (extinction and scattering in optical bands and MIR emission). The amount of dust mass at distances larger than $\sim 2$ kpc from the midplane represents $2 - 3.3$ % of the total galactic dust mass and the relative abundance of small grains with respect to large grains is almost halved comparing to that in the midplane. The paucity of small grains high above the midplane might indicate that dust is hit by interstellar shocks or galactic fountains and entrained together with gas. The halo dust component is likely to be embedded in an atomic / molecular gas and heated by a thick stellar disc.

CR7 is the brightest $z=6.6 \, {\rm Ly}\alpha$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}\alpha$/HeII1640 line luminosities ($L_{\rm \alpha/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}\alpha$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_{\alpha} >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_{\alpha}$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.

AGN-driven outflows are invoked by galaxy evolutionary models to quench star formation and to explain the origin of the relations observed locally between super massive black holes and their host galaxies. This work aims to detect the presence of extended ionised outflows in luminous quasars where we expect the maximum activity both in star formation and in black hole accretion. Currently, there are only a few studies based on spatially resolved observations of outflows at high redshift, $z>2$. We analyse a sample of six luminous (${\rm L>10^{47} \ erg/s}$) quasars at $z\sim2.4$, observed in H-band using the near-IR integral field spectrometer SINFONI at VLT. We perform a kinematic analysis of the [OIII] emission line at $\lambda = 5007\AA$. [OIII] has a complex gas kinematic, with blue-shifted velocities of a few hundreds of km/s and line widths up to 1500 km/s. Using the spectroastrometric method we infer size of the ionised outflows of up to $\sim$2 kpc. The properties of the ionised outflows, mass outflow rate, momentum rate and kinetic power, are correlated with the AGN luminosity. The increase in outflow rate with increasing AGN luminosity is consistent with the idea that a luminous AGN pushes away the surrounding gas through fast outflows driven by radiation pressure, which depends on the emitted luminosity. We derive mass outflow rates of about 6-700 M$_{\odot}$/yr for our sample, which are lower than those observed in molecular outflows. Indeed physical properties of ionised outflows show dependences on AGN luminosity which are similar to those of molecular outflows but indicating that the mass of ionised gas is smaller than that of the molecular one. Alternatively, this discrepancy between ionised and molecular outflows could be explained with different acceleration mechanisms.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We derive the modified Friedmann equations for a generalization of the Dvali-Gabadadze-Porrati (DGP) model in which the brane has one additional compact dimension. The main new feature is the emission of gravitational waves into the bulk. We study two classes of solutions: First, if the compact dimension is stabilized, the waves vanish and one exactly recovers DGP cosmology. However, a stabilization by means of physical matter is not possible for a tension-dominated brane, thus implying a late time modification of 4D cosmology different from DGP. Second, for a freely expanding compact direction, we find exact attractor solutions with zero 4D Hubble parameter despite the presence of a 4D cosmological constant. The model hence constitutes an explicit example of dynamical degravitation at the full non-linear level. Without stabilization, however, there is no 4D regime and the model is ruled out observationally, as we demonstrate explicitly by comparing to supernova data.

In General Relativity, cosmic strings are well known to produce a static, locally flat spacetime with a wedge removed. If the tension exceeds a critical value, the deficit angle becomes larger than $ 2\pi $, leading to a compact exterior that ends in a conical singularity. In this work, we investigate dynamical solutions for cosmic strings with super-critical tensions. To this end, we model the string as a cylindrical shell of finite and stabilized transverse width and show that there is a marginally super-critical regime in which the stabilization can be achieved by physically reasonable matter. We show numerically that the static deficit angle solution is unstable for super-critical string tensions. Instead, the geometry starts expanding in axial direction at an asymptotically constant rate, and a horizon is formed in the exterior, which has the shape of a growing cigar. We are able to find the analytic form of the attractor solution describing the interior of the cosmic string. In particular, this enables us to analytically derive the relation between the string tension and the axial expansion rate. Furthermore, we show that the exterior conical singularity can be avoided for dynamical solutions. Our results might be relevant for theories with two extra dimensions, modeling our universe as a cosmic string with a three-dimensional axis. We derive the corresponding Friedmann equation, relating the on-brane Hubble parameter to the string tension or, equivalently, brane cosmological constant.

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 use Spitzer observations of the rich population of Asymptotic Giant Branch stars in the Large Magellanic Cloud (LMC) to test models describing the internal structure and nucleosynthesis of the most massive of these stars, i.e. those with initial mass above $\sim 4M_{\odot}$. To this aim, we compare Spitzer observations of LMC stars with the theoretical tracks of Asymptotic Giant Branch models, calculated with two of the most popular evolution codes, that are known to differ in particular for the treatment of convection. Although the physical evolution of the two models are significantly different, the properties of dust formed in their winds are surprisingly similar, as is their position in the colour-colour (CCD) and colour-magnitude (CMD) diagrams obtained with the Spitzer bands. This model independent result allows us to select a well defined region in the ($[3.6]-[4.5], [5.8]-[8.0]$) plane, populated by AGB stars experiencing Hot Bottom Burning, the progeny of stars with mass $M\sim 5.5M_{\odot}$. This result opens up an important test of the strength hot bottom burning using detailed near-IR (H and K bands) spectroscopic analysis of the oxygen-rich, high luminosity candidates found in the well defined region of the colour-colour plane. This test is possible because the two stellar evolution codes we use predict very different results for the surface chemistry, and the C/O ratio in particular, owing to their treatment of convection in the envelope and of convective boundaries during third dredge-up. The differences in surface chemistry are most apparent when the model stars reach the phase with the largest infrared emission.

Here we introduce GAMESH, a novel pipeline which implements self-consistent radiative and chemical feedback in a computational model of galaxy formation. By combining the cosmological chemical-evolution model GAMETE with the radiative transfer code CRASH, GAMESH can post process realistic outputs of a N-body simulation describing the redshift evolution of the forming galaxy. After introducing the GAMESH implementation and its features, we apply the code to a low-resolution N-body simulation of the Milky Way formation and we investigate the combined effects of self-consistent radiative and chemical feedback. Many physical properties, which can be directly compared with observations in the Galaxy and its surrounding satellites, are predicted by the code along the merger-tree assembly. The resulting redshift evolution of the Local Group star formation rates, reionisation and metal enrichment along with the predicted Metallicity Distribution Function of halo stars are critically compared with observations. We discuss the merits and limitations of the first release of GAMESH, also opening new directions to a full implementation of feedback processes in galaxy formation models by combining semi-analytic and numerical methods.

The cosmology of brane induced gravity in six infinite dimensions is investigated. It is shown that a brane with Friedmann-Robertson-Walker symmetries necessarily acts as a source of cylindrically symmetric gravitational waves, so called Einstein-Rosen waves. Their existence essentially distinguishes this model from its codimension-one counterpart and necessitates solving the nonlinear system of bulk and brane-matching equations. A numerical analysis is performed and two qualitatively different and dynamically separated classes of cosmologies are derived: degravitating solutions for which the Hubble parameter settles to zero despite the presence of a non-vanishing energy density on the brane and super-accelerating solutions for which Hubble grows unbounded. The parameter space of both the stable and unstable regime is derived and observational consequences are discussed: It is argued that the degravitating regime does not allow for a phenomenologically viable cosmology. On the other hand, the super-accelerating solutions are potentially viable, however, their unstable behavior questions their physical relevance.

We calculated theoretical evolutionary sequences of asymptotic giant branch (AGB) stars, including formation and evolution of dust grains in their circumstellar envelope. By considering stellar populations of the Large Magellanic Cloud (LMC), we calculate synthetic colour-colour and colour-magnitude diagrams, which are compared with those obtained by the Spitzer Space Telescope. The comparison between observations and theoretical predictions outlines that extremely obscured carbon-stars and oxygen-rich sources experiencing hot bottom burning (HBB) occupy well defined, distinct regions in the colour-colour ($[3.6]-[4.5]$, $[5.8]-[8.0]$) diagram. The C-rich stars are distributed along a diagonal strip that we interpret as an evolutionary sequence, becoming progressively more obscured as the stellar surface layers enrich in carbon. Their circumstellar envelopes host solid carbon dust grains with size in the range $0.05 < a < 0.2 \mu m$. The presence of SiC particles is expected only in the more metal-rich stars. The reddest sources, with $[3.6]-[4.5] > 2$, are the descendants of stars with initial mass $M_{in} \sim 2.5 - 3 M_{\odot}$ in the very latest phases of the AGB life. The oxygen-rich stars with the reddest colours ($[5.8]-[8.0] > 0.6$) are those experiencing HBB, the descendants of $\sim 5 M_{\odot}$ objects formed $10^{8}$ yr ago; alumina and silicates dust start forming at different distances from the central star. The overall dust production rate in the LMC is $\sim 4.5 \times 10^{-5} M_{\odot}/yr$, the relative percentages due to C- and M- star being respectively 85$%$ and 15 $%$.

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

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 origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

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.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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 O-rich Asymptotic Giant Branch (AGB) stars experience strong mass loss with efficient dust condensation and they are major sources of dust in the interstellar medium. Alumina dust (Al$_2$O$_3$) is an important dust component in O-rich circumstellar shells and it is expected to be fairly abundant in the winds of the more massive and O-rich AGB stars. By coupling AGB stellar nucleosynthesis and dust formation, we present a self-consistent exploration on the Al$_2$O$_3$ production in the winds of AGB stars with progenitor masses between $\sim$3 and 7 M$_{\odot}$ and metallicities in the range 0.0003 $\le$ Z $\le$ 0.018. We find that Al$_2$O$_3$ particles form at radial distances from the centre between $\sim2$ and 4 R$_*$ (depending on metallicity), which is in agreement with recent interferometric observations of Galactic O-rich AGB stars. The mass of Al$_2$O$_3$ dust is found to scale almost linearly with metallicity, with solar metallicity AGBs producing the highest amount (about 10$^{-3}$ M$_{\odot}$) of alumina dust. The Al$_2$O$_3$ grain size decreases with decreasing metallicity (and initial stellar mass) and the maximum size of the Al$_2$O$_3$ grains is $\sim$0.075 $\mu m$ for the solar metallicity models. Interestingly, the strong depletion of gaseous Al observed in the low-metallicity HBB AGB star HV 2576 seems to be consistent with the formation of Al$_2$O$_3$ dust as predicted by our models. We suggest that the content of Al may be used as a mass (and evolutionary stage) indicator in AGB stars experiencing HBB.

We trace the full evolution of low- and intermediate-mass stars ($1 M_{\odot} \leq M \leq 8M_{\odot}$) during the Asymptotic Giant Branch (AGB) phase in the {\it Spitzer} two-color and color-magnitude diagrams. We follow the formation and growth of dust particles in the circumstellar envelope with an isotropically expanding wind, in which gas molecules impinge upon pre--existing seed nuclei, favour their growth. These models are the first able to identify the main regions in the {\it Spitzer} data occupied by AGB stars in the Large Magellanic Cloud (LMC). The main diagonal sequence traced by LMC extreme stars in the [3.6]-[4.5] vs. [5.8]-[8.0] and [3.6]-[8.0] vs. [8.0] planes are nicely fit by carbon stars models; it results to be an evolutionary sequence with the reddest objects being at the final stages of their AGB evolution. The most extreme stars, with [3.6]-[4.5] $>$ 1.5 and [3.6]-[8.0] $>$ 3, are 2.5-3 $M_{\odot}$ stars surrounded by solid carbon grains. In higher mass ($>3 M_{\odot}$) models dust formation is driven by the extent of Hot Bottom Burning (HBB) - most of the dust formed is in the form of silicates and the maximum obscuration phase by dust particles occurs when the HBB experienced is strongest, before the mass of the envelope is considerably reduced.

Dust grains in low-metallicity star-forming regions may be responsible for the formation of the first low-mass stars. The minimal conditions to activate dust-induced fragmentation require the gas to be pre-enriched above a critical dust-to-gas mass ratio Dcr=[2.6--6.3]x10^-9 with the spread reflecting the dependence on the grain properties. The recently discovered Galactic halo star SDSS J102915+172927 has a stellar mass of 0.8 Msun and a metallicity of Z=4.5x10^-5 Zsun and represents an optimal candidate for the dust-induced low-mass star formation. Indeed, for the two most plausible Population III supernova progenitors, with 20 Msun and 35 Msun, the critical dust-to-gas mass ratio can be overcome provided that at least 0.4 Msun of dust condenses in the ejecta, allowing for moderate destruction by the reverse shock. Here we show that even if dust formation in the first supernovae is less efficient or strong dust destruction does occur, grain growth during the collapse of the parent gas cloud is sufficiently rapid to activate dust cooling and likely fragmentation into low-mass and long-lived stars. Silicates and magnetite grains can experience significant grain growth in the density range 10^9 /cc < nH<10^12 /cc by accreting gas-phase species (SiO, SiO2, and Fe) until their gas-phase abundance drops to zero, reaching condensation efficiencies =1. The corresponding increase in the dust-to-gas mass ratio allows dust-induced cooling and fragmentation to be activated at 10^12 /cc < nH < 10^14 /cc, before the collapsing cloud becomes optically thick to continuum radiation. We show that for all the initial conditions that apply to the parent cloud of SDSS J102915+172927, dust-driven fragmentation is able to account for the formation of the star.

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\%$).

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

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.

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.

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.

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.

In this paper we perform a detailed analysis of the effect of various approximations which have been used in the literature to compute the detection rates of compact binary coalescences for interferometric gravitational wave detectors. We evaluate the detection rates for the coalescence of BH-BH, NS-NS, and BH-NS binaries taking into account their specific statistical properties obtained from population synthesis models, the cosmic star formation rate history, and the effects of redshift on the emitted gravitational wave signals. The results are compared with those obtained with procedures that are based on different levels of approximations, such as adopting averaged values for the total mass and symmetric mass ratio for all the systems of a binary population, using these to compute the horizon distance for individual detectors, or estimating the coalescence rate density within this distance by its local value. We show that these approximations introduce a bias which depends not only on the approximation which is used, but also on the statistical properties of the considered binary population.

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.

The relative importance of metals and dust grains in the formation of the first low-mass stars has been a subject of debate. The recently discovered Galactic halo star SDSS J102915+172927 (Caffau et al. 2011) has a mass less than 0.8 Msun and a metallicity of Z = 4.5 10^{-5} Zsun. We investigate the origin and properties of this star by reconstructing the physical conditions in its birth cloud. We show that the observed elemental abundance trend of SDSS J102915+172927 can be well fitted by the yields of core-collapse supernovae with metal-free progenitors of 20 Msun and 35 Msun. Using these selected supernova explosion models, we compute the corresponding dust yields and the resulting dust depletion factor taking into account the partial destruction by the supernova reverse shock. We then follow the collapse and fragmentation of a star forming cloud enriched by the products of these SN explosions at the observed metallicity of SDSS J102915+172927. We find that [0.05 - 0.1] Msun mass fragments, which then lead to the formation of low-mass stars, can occur provided that the mass fraction of dust grains in the birth cloud exceeds 0.01 of the total mass of metals and dust. This, in turn, requires that at least 0.4 Msun of dust condense in the first supernovae, allowing for moderate destruction by the reverse shock. If dust formation in the first supernovae is less efficient or strong dust destruction does occur, the thermal evolution of the SDSS J102915+172927 birth cloud is dominated by molecular cooling, and only > 8 Msun fragments can form. We conclude that the observed properties of SDSS J102915+172927 support the suggestion that dust must have condensed in the ejecta of the first supernovae and played a fundamental role in the formation of the first low-mass stars.

We compute the gravitational wave background (GWB) generated by a cosmological population of (BH-BH) binaries using hybrid waveforms recently produced by numerical simulations of (BH-BH) coalescence, which include the inspiral, merger and ring-down contributions. A large sample of binary systems is simulated using the population synthesis code SeBa, and we extract fundamental statistical information on (BH-BH) physical parameters (primary and secondary BH masses, orbital separations and eccentricities, formation and merger timescales). We then derive the binary birth and merger rates using the theoretical cosmic star formation history obtained from a numerical study which reproduces the available observational data at redshifts $z < 8$. We evaluate the contributions of the inspiral, merger and ring-down signals to the GWB, and discuss how these depend on the parameters which critically affect the number of coalescing (BH-BH) systems. We find that Advanced LIGO/Virgo have a chance to detect the GWB signal from the inspiral phase with a $(S/N)=10$ only for the most optimistic model, which predicts the highest local merger rate of 0.85 Mpc$^{-3}$ Myr$^{-1}$. Third generation detectors, such as ET, could reveal the GWB from the inspiral phase predicted by any of the considered models. In addition, ET could sample the merger phase of the evolution at least for models which predict local merger rates between $[0.053 - 0.85]$ Mpc$^{-3}$ Myr$^{-1}$, which are more than a factor 2 lower the the upper limit inferred from the analysis of the LIGO S5 run\cite{Abadieetal2011}.

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.

This paper has been withdrawn by the authors because it had major revisions.

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.

We include an energy term based on Dark Matter (DM) self-annihilation during the cooling and subsequent collapse of the metal-free gas, in halos hosting the formation of the first stars in the Universe. We have found that the feedback induced on the chemistry of the cloud does modify the properties of the gas throughout the collapse. However, the modifications are not dramatic, and the typical Jeans mass within the halo is conserved throughout the collapse, for all the DM parameters we have considered. This result implies that the presence of Dark Matter annihilations does not substantially modify the Initial Mass Function of the First Stars, with respect to the standard case in which such additional energy term is not taken into account. We have also found that when the rate of energy produced by the DM annihilations and absorbed by the gas equals the chemical cooling (at densities yet far from the actual formation of a proto-stellar core) the structure does not halt its collapse, although that proceeds more slowly by a factor smaller than few per cent of the total collapse time.

We analyze the optical-near infrared spectra of 33 quasars with redshifts 3.9<z<6.4 to investigate the properties of dust extinction at these cosmic epochs. The SMC extinction curve has been shown to reproduce the dust reddening of most quasars at z<2.2; we investigate whether this curve also provides a good description of dust extinction at higher redshifts. We fit the observed spectra with synthetic absorbed quasar templates obtained by varying the intrinsic slope (alpha), the absolute extinction (A3000), and by using a grid of empirical and theoretical extinction curves. We find that seven quasars in our sample are affected by substantial extinction (A3000 >0.8), and characterized by very steep intrinsic slopes (alpha<-2.3). All of the individual quasars require extinction curve deviating from that of the SMC, with a tendency to flatten at lambda<2000 A (in the rest frame of the source). We obtain a mean extinction curve at z>4, both by performing a simultaneous fit of all quasars and by averaging the extinction curves inferred for individual quasars. In the case of broad absorption line quasars the mean extinction curve deviates from the SMC at a confidence level >95%. The difference between extinction curves in quasars at z>4 and in quasars at lower redshift is indicative of either a different dust production mechanism at high redshift, or a different mechanism for processing dust into the ISM. We suggest that the same transitions may also apply to normal, star-forming galaxies at z>4. In particular, the observed change in the average spectral slope of galaxies at z>4 may be partially ascribed to a variation in the extinction curve, rather than a lower dust content at high redshift. In this scenario, the extinction curve inferred at z>4 would imply a cosmic star-formation corrected for dust attenuation a factor of ~2 higher than estimated in the past.

Two classes of high energy sources in our galaxy are believed to host magnetars, neutron stars whose emission results from the dissipation of their magnetic field. The extremely high magnetic field of magnetars distorts their shape, and causes the emission of a conspicuous gravitational waves signal if rotation is fast and takes place around a different axis than the symmetry axis of the magnetic distortion. Based on a numerical model of the cosmic star formation history, we derive the cosmological background of gravitational waves produced by magnetars, when they are very young and fast spinning. We adopt different models for the configuration and strength of the internal magnetic field (which determines the distortion) as well as different values of the external dipole field strength (which governs the spin evolution of magnetars over a wide range of parameters). We find that the expected gravitational wave background differs considerably from one model to another. The strongest signals are generated for magnetars with very intense toroidal internal fields ($\sim 10^{16}$ G range) and external dipole fields of $\sim 10^{14}$, as envisaged in models aimed at explaining the properties of the Dec 2004 giant flare from SGR 1806-20. Such signals should be easily detectable with third generation ground based interferometers such as the Einstein Telescope.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

Astrophysical sources emit gravitational waves in a large variety of processes occurred since the beginning of star and galaxy formation. These waves permeate our high redshift Universe, and form a background which is the result of the superposition of different components, each associated to a specific astrophysical process. Each component has different spectral properties and features that it is important to investigate in view of a possible, future detection. In this contribution, we will review recent theoretical predictions for backgrounds produced by extragalactic sources and discuss their detectability with current and future gravitational wave observatories.

We study the age and metallicity distribution function (MDF) of metal-poor stars in the Milky Way halo as a function of galactocentric radius by combining N-body simulations and semi-analytical methods. We find that the oldest stars populate the innermost region, while extremely metal-poor stars are more concentrated within r < 60 kpc. The MDF of [Fe/H] < -2 stars varies only very weakly within the central 50 kpc, while the relative contribution of [Fe/H] < -2 stars strongly increases with r, varying from 16% within 7 kpc < r < 20 kpc up to > 40% for r > 20 kpc. This is due to the faster descent of the spatial distribution (as seen from Earth) of the more enriched population. This implies that the outer halo < 40 kpc is the best region to search for very metal-poor stars. Beyond ~ 60 kpc the density of [Fe/H] < -2 stars is maximum within dwarf galaxies. All these features are imprinted by a combination of (i) the virialization epoch of the star-forming haloes, and (ii) the metal enrichment history of the Milky Way environment.

We discuss expectations for the absorption of high-energy gamma-rays by gamma-gamma pair production with intergalactic radiation fields (IRFs) at very high redshifts (z~5-20), and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, a semi-analytical model incorporating both Population II and Population III stars is employed, which is consistent with a wide variety of existing high-z observations including QSO spectral measurements, WMAP Thomson depth constraints, near-IR source count limits, etc. We find that the UV IRF below the Lyman edge energy with intensities in the range of a few times 10^{-19} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1} can cause appreciable attenuation above ~12 GeV at z~5, down to ~6-8 GeV at z>~8-10. This may be observable in the spectra of blazars or gamma-ray bursts by the Fermi Gamma-ray Space Telescope or next generation facilities such as the Cherenkov Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable insight into early star formation and cosmic reionization.

We investigate the effects of the Cosmic Microwave Background (CMB) radiation field on the collapse of prestellar clouds. Using a semi-analytic model to follow the thermal evolution of clouds with varying initial metallicities and dust contents at different redshifts, we study self-consistently the response of the mean Jeans mass at cloud fragmentation to metal line-cooling, dust-cooling and the CMB. In the absence of dust grains, at redshifts z < 10 moderate characteristic masses (of 10s of Msun) are formed when the metallicity is 10^{-4} Zsun < Z < 10^{-2.5} Zsun; at higher metallicities, the CMB inhibits fragmentation and only very large masses (of ~ 100s of Msun) are formed. These effects become even more dramatic at z > 10 and the fragmentation mass scales are always > 100s of Msun, independent of the initial metallicity. When dust grains are present, sub-solar mass fragments are formed at any redshift for metallicities Z > 10^{-6} Zsun because dust-cooling remains relatively insensitive to the presence of the CMB. When Z > 10^{-3} Zsun, heating of dust grains by the CMB at z > 5 favors the formation of larger masses, which become super-solar when Z > 10^{-2} Zsun and z > 10. Finally, we discuss the implications of our result for the interpretation of the observed abundance patterns of very metal-poor stars in the galactic halo.

We explore the effects of magnetic energy dissipation on the formation of the first stars. For this purpose, we follow the evolution of primordial chemistry in the presence of magnetic fields in the post-recombination universe until the formation of the first virialized halos. From the point of virialization, we follow the protostellar collapse up to densities of $\sim10^{12}$ cm$^{-3}$ in a one-zone model. In the intergalactic medium (IGM), comoving field strengths of $\gtrsim0.1$ nG lead to Jeans masses of $10^8 M_\odot$ or more and thus delay gravitational collapse in the first halos until they are sufficiently massive. During protostellar collapse, we find that the temperature minimum at densities of $\sim10^3$ cm$^{-3}$ does not change significantly, such that the characteristic mass scale for fragmentation is not affected. However, we find a significant temperature increase at higher densities for comoving field strengths of $\gtrsim0.1$ nG. This may delay gravitational collapse, in particular at densities of $\sim10^9$ cm$^{-3}$, where the proton abundance drops rapidly and the main contribution to the ambipolar diffusion resistivity is due to collisions with Li$^+$. After the formation of the protostar, the increased gas temperatures may enhance the protostellar accretion rate. Our model confirms that initial weak magnetic fields may be amplified considerably during gravitational collapse and become dynamically relevant. For instance, a comoving field strength above $10^{-5}$ nG will be amplified above the critical value for the onset of jets which can magnetize the IGM.

We present the first results of a project, LSD, aimed at obtaining spatially-resolved, near-infrared spectroscopy of a complete sample of Lyman-Break Galaxies at z~3. Deep observations with adaptive optics resulted in the detection of the main optical lines, such as [OII], Hbeta and [OIII], which are used to study sizes, SFRs, morphologies, gas-phase metallicities, gas fractions and effective yields. Optical, near-IR and Spitzer/IRAC photometry is used to measure stellar mass. We obtain that morphologies are usually complex, with the presence of several peaks of emissions and companions that are not detected in broad-band images. Typical metallicities are 10-50% solar, with a strong evolution of the mass-metallicity relation from lower redshifts. Stellar masses, gas fraction, and evolutionary stages vary significantly among the galaxies, with less massive galaxies showing larger fractions of gas. In contrast with observations in the local universe, effective yields decrease with stellar mass and reach solar values at the low-mass end of the sample. This effect can be reproduced by gas infall with rates of the order of the SFRs. Outflows are present but are not needed to explain the mass-metallicity relation. We conclude that a large fraction of these galaxies are actively creating stars after major episodes of gas infall or merging.

Using the results of a numerical simulation which follows the evolution, metal enrichment and energy deposition of both Population III and Population II stars, we predict the redshift dependence of the formation rate of black hole remnants of Population III stars with masses 100- 500Msun and of neutron stars(black holes) remnants of Population II stars with masses 8-20Msun (20-40Msun). We describe the gravitational wave spectrum produced by Population III and Population II sources adopting the most appropriate signals available in the literature and we compute the stochastic backgrounds resulting from the cumulative emission of these sources throughout the history of the Universe. With the aim of assessing whether these backgrounds might act as foregrounds for signals generated in the Inflationary epoch, we compare their amplitudes with the sensitivity of currently planned and future ground/space-based interferometers.

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.

We study the effects of WIMP dark matter (DM) annihilations on the thermal and chemical evolution of the gaseous clouds where the first generation of stars in the Universe is formed. We follow the collapse of the gas inside a typical halo virializing at very high redshift, from well before virialization until a stage where the heating from DM annihilations exceeds the gas cooling rate. The DM energy input is estimated by inserting the energy released by DM annihilations (as predicted by an adiabatic contraction of the original DM profile) in a spherically symmetric radiative transfer scheme. In addition to the heating effects of the energy absorbed, we include its feedback upon the chemical properties of the gas, which is critical to determine the cooling rate in the halo, and hence the fragmentation scale and Jeans mass of the first stars. We find that DM annihilation does alter the free electron and especially the H2 fraction when the gas density is n>~ 10^4 cm^-3, for our fiducial parameter values. However, even if the change in the H2 abundance and the cooling efficiency of the gas is large (sometimes exceeding a factor 100), the effects on the temperature of the collapsing gas are far smaller (a reduction by a factor <~1.5), since the gas cooling rate depends very strongly on temperature: then, the fragmentation mass scale is reduced only slightly, hinting towards no dramatic change in the initial mass function of the first stars.

We study the evolution of the first stars in the universe (Population III) from the early pre-Main Sequence until the end of helium burning in the presence of WIMP dark matter annihilation inside the stellar structure. The two different mechanisms that can provide this energy source are the contemporary contraction of baryons and dark matter, and the capture of WIMPs by scattering off the gas with subsequent accumulation inside the star. We find that the first mechanism can generate an equilibrium phase, previously known as a "dark star", which is transient and present in the very early stages of pre-MS evolution. The mechanism of scattering and capture acts later, and can support the star virtually forever, depending on environmental characteristic of the dark matter halo and on the specific WIMP model.

Primordial gas in protogalactic dark matter (DM) halos with virial temperatures Tvir > 10^4 K begins to cool and condense via atomic hydrogen. Provided this gas is irradiated by a strong ultraviolet (UV) flux and remains free of H2 and other molecules, it has been proposed that the halo with Tvir ~10^4 K may avoid fragmentation, and lead to the rapid formation of a supermassive black hole (SMBH) as massive as M=10^5-10^6 Msun. This ``head--start'' would help explain the presence of SMBHs with inferred masses of several x 10^9 Msun, powering the bright quasars discovered in the Sloan Digital Sky Survey at redshift z>~6. However, high-redshift DM halos with Tvir~10^4K are likely already enriched with at least trace amounts of metals and dust produced by prior star-formation in their progenitors. Here we study the thermal and chemical evolution of low-metallicity gas exposed to extremely strong UV radiation fields. Our results, obtained in one-zone models, suggest that gas fragmentation is inevitable above a critical metallicity, whose value is between Zcr~3x10^{-4} Zsun (in the absence of dust) and as low as Zcr~ 5 x 10^{-6} Zsun (with a dust-to-gas mass ratio of about 0.01 Z/Zsun). We propose that when the metallicity exceeds these critical values, dense clusters of low--mass stars may form at the halo nucleus. Relatively massive stars in such a cluster can then rapidly coalesce into a single more massive object, which may produce an intermediate-mass BH remnant with a mass up to M <~10^2-10^3 Msun.

We explore the effect of cosmic radiative feedback from the sources of reionization on the thermal evolution of the intergalactic medium. We find that different prescriptions for this feedback predict quite different thermal and reionization histories. In spite of this, current data can not discriminate among different reionization scenarios. We find that future observations both from 21-cm and CMB experiments can be used to break the degeneracy among model parameters provided that we will be able to remove the foreground signal at the percent (or better) level.

We study the effects of WIMP dark matter (DM) on the collapse and evolution of the first stars in the Universe. Using a stellar evolution code, we follow the pre-Main Sequence (MS) phase of a grid of metal-free stars with masses in the range 5-600 solar mass forming in the centre of a 1e6 solar mass halo at redhisft z=20. DM particles of the parent halo are accreted in the proto-stellar interior by adiabatic contraction and scattering/capture processes, reaching central densities of order 1e12 GeV/cm3 at radii of the order of 10 AU. Energy release from annihilation reactions can effectively counteract the gravitational collapse, in agreement with results from other groups. We find this stalling phase (known as "dark" star) is transients and lasts from 2.1e3 yr (M=600 solar mass) to 1.8e4 yr (M=9 solar mass). Later in the evolution, DM scattering/capture rate becomes high enough that energy deposition from annihilations significantly alters the pre-MS evolution of the star in a way that depends on DM (i) velocity dispersion, (ii) density, (iii) elastic scattering cross section with baryons. For our fiducial set of parameters (10 km/s, 1e11 GeV/cm3, 1e-38 cm2) we find that the evolution of stars of mass lower than 40 solar masses "freezes" on the HR diagram before reaching the ZAMS. Stars with bigger masses manage to ignite nuclear reactions; however, DM "burning" prolonges their lifetimes by a factor 2 (5) for a 600 (40) solar mass star.

We propose a cosmological scenario for the formation and evolution of dwarf spheroidal galaxies (dSphs), satellites of the Milky Way (MW). An improved version of the semi-analytical code GAMETE (GAlaxy Merger Tree & Evolution) is used to follow the dSphs evolution simultaneously with the MW formation, matching the observed properties of both. In this scenario dSph galaxies represent fossil objects virializing at z = 7.2 +/- 0.7 (i.e. in the pre-reionization era z > z_rei = 6) in the MW environment, which at that epoch has already been pre-enriched up to [Fe/H] ~ -3; their dynamical masses are in the narrow range M = (1.6 +/- 0.7) x 10^8 M_sun, although a larger spread might be introduced by a more refined treatment of reionization. Mechanical feedback effects are dramatic in such low-mass objects, causing the complete blow-away of the gas ~100 Myr after the formation epoch: 99% of the present-day stellar mass, M_* = (3 +/- 0.7) x 10^6 M_sun, forms during this evolutionary phase, i.e. their age is >13 Gyr. Later on, star formation is re-ignited by returned gas from evolved stars and a second blow-away occurs. The cycle continues for about 1 Gyr during which star formation is intermittent. At z = 0 the dSph gas content is M_g = (2.68 +/- 0.97) x 10^4 M_sun. Our results match several observed properties of Sculptor, used as a template of dSphs: (i) the Metallicity Distribution Function; (ii) the Color Magnitude Diagram; (iii) the decrement of the stellar [O/Fe] abundance ratio for [Fe/H] > -1.5; (iv) the dark matter content and the light-to-mass ratio; (v) the HI gas mass content.

We compute the imprints left on the CMB by two cosmic reionization models consistent with current observations but characterized by alternative radiative feedback prescriptions (suppression and filtering) resulting in a different suppression of star formation in low-mass halos. The models imply different ionization and thermal histories and 21 cm background signals. The derived Comptonization, u, and free-free distortion, y_B, parameters are below current observational limits for both models. However, the value of u = 1.69 * 10^-7 (9.65 * 10^-8) for the suppression (filtering) model is in the detectability range of the next generation of CMB spectrum experiments. Through the dedicated Boltzmann code CMBFAST, modified to include the above ionization histories, we compute the CMB angular power spectrum (APS) of the TT, TE, and EE modes. For the EE mode the differences between these models are significantly larger than the cosmic and sampling variance over the multipole range l = 5-15, leaving a good chance of discriminating between these feedback mechanisms with forthcoming/future CMB polarization experiments. The main limitations come from foreground contamination: it should be subtracted at per cent level in terms of APS, a result potentially achievable by novel component separation techniques and mapping of Galactic foreground.

We use a semi-analytical model to study the impact of reionization, and the associated radiative feedback, on galaxy formation. Two feedback models have been considered: (i) a standard prescription, according to which star formation is totally suppressed in galaxies with circular velocity below a critical threshold (model CF06) and (ii) a characterization based on the filtering scale (model G00), allowing for a gradual reduction of the gas available for star formation in low-mass galaxies. In model CF06 reionization starts at z ~ 15-20, is 85% complete by z ~ 10; at the same z, the ionized fraction is 16% in model G00. The models match SDSS constraints on the evolution of the neutral hydrogen fraction at z < 7, but predict different Thomson optical depths, tau_e = 0.1017 (CF06), and 0.0631 (G00); such values are within 1 sigma of the WMAP 3-yr determination. Both models are in remarkable good agreement with additional existing data (evolution of Lyman-limit systems, cosmic star formation history, high-z galaxy counts, IGM thermal history), which therefore cannot be used to discriminate among different feedback models. Deviations among radiative feedback prescriptions emerge when considering the expected HI 21 cm background signal, where a ~ 15 mK absorption feature in the range 75-100 MHz is present in model G00 and a global shift of the emission feature preceding reionization towards larger frequencies occurs in the same model. Single dish observations with existing or forthcoming low-frequency radio telescopes can achieve mK sensitivity, allowing the identification of these features provided that foregrounds can be accurately subtracted.

A PopIII/Pop II transition from massive to normal stars is predicted to occur when the metallicity of the star forming gas crosses the critical range Z_cr = 10^(-5 +/- 1) Z_sun. To investigate the cosmic implications of such process we use numerical simulations which follow the evolution, metal enrichment and energy deposition of both Pop III and Pop II stars. We find that: (i) due to inefficient heavy element transport by outflows and slow "genetic" transmission during hierarchical growth, large fluctuations around the average metallicity arise; as a result Pop III star formation continues down to z=2.5, but at a low peak rate of 10^-5 M_sun yr^-1 Mpc^-3 occurring at z~6 (about 10^-4 of the PopII one); (ii) Pop III star formation proceeds in a "inside-out" mode in which formation sites are progressively confined at the periphery of collapsed structures, where the low gas density and correspondingly long free-fall timescales result in a very inefficient astration. These conclusions strongly encourage deep searches for pristine star formation sites at moderate (2<z<5) redshifts where metal free stars are likely to be hidden.

The presence of dust at high redshift requires efficient condensation of grains in SN ejecta, in accordance with current theoretical models. Yet, observations of the few well studied SNe and SN remnants imply condensation efficiencies which are about two orders of magnitude smaller. Motivated by this tension, we have (i) revisited the model of Todini & Ferrara (2001) for dust formation in the ejecta of core collapse SNe and (ii) followed, for the first time, the evolution of newly condensed grains from the time of formation to their survival - through the passage of the reverse shock - in the SN remnant. We find that 0.1 - 0.6 M_sun of dust form in the ejecta of 12 - 40 M_sun stellar progenitors. Depending on the density of the surrounding ISM, between 2-20% of the initial dust mass survives the passage of the reverse shock, on time-scales of about 4-8 x 10^4 yr from the stellar explosion. Sputtering by the hot gas induces a shift of the dust size distribution towards smaller grains. The resulting dust extinction curve shows a good agreement with that derived by observations of a reddened QSO at z =6.2. Stochastic heating of small grains leads to a wide distribution of dust temperatures. This supports the idea that large amounts (~ 0.1 M_sun) of cold dust (T ~ 40K) can be present in SN remnants, without being in conflict with the observed IR emission.

We study the stellar population history and chemical evolution of the Milky Way (MW) in a hierarchical LCDM model for structure formation. Using a Monte Carlo method based on the semi-analytical EPS formalism, we reconstruct the merger tree of our Galaxy and follow the evolution of gas and stars along the hierarchy. Our approach allows us to compare the observational properties of the MW with model results, exploring different properties of primordial stars, such as their IMF and the critical metallicity for low-mass star formation, Zcr. By matching our predictions to the Metallicity Distribution Function (MDF) of metal-poor stars in the Galactic halo we find that: (i) supernova feedback is required to reproduce the observed properties of the MW; (ii) stars with [Fe/H]<-2.5 form in halos accreting Galactic Medium (GM) enriched by earlier supernova explosions; (iii) the fiducial model (Zcr=10^-4Zsun,m_PopIII=200Msun) provides an overall good fit to the MDF but cannot account for the two stars with [Fe/H]<-5; the latter can be accommodated if Zcr<10^-6Zsun but such model overpopulates the range -5.3<[Fe/H]<-4 in which no stars have been detected; (iv) the current non-detection of metal-free stars robustly constrains either Zcr>0 or the masses of the first stars m_PopIII>0.9Msun; (v) the statistical impact of second generation stars, i.e stars forming out of gas polluted only by metal-free stars, is negligible in current samples; (vi) independently of Zcr, 60% of metals in the GM are ejected through winds by halos with masses M<6x10^9 Msun, showing that low-mass halos are the dominant population contributing to cosmic metal enrichment.

We study the spatial distribution of Galactic metal-free stars by combining an extremely high-resolution (7.8 X 10^5 solar masses per particle) Cold Dark Matter N-body simulation of the Milky-Way with a semi-analytic model of metal enrichment. This approach allows us to resolve halos with virial temperatures down to the 10^4K atomic cooling limit, and it is sufficiently flexible to make a number of robust conclusions, despite the extremely uncertain properties of the first stars. Galactic metal-free stars are formed over a large redshift range, which peaks at z~10, but continues down to z~5, contributing stars at wide range of Galactocentric radii. Stars containing only metals generated by primordial stars are similarly widespread. Neither changing the efficiency of metal dispersal by two orders of magnitude, nor drastically changing the approximations in our semi-analytical model can affect these result. Thus, if they have sufficiently long lifetimes, a significant number of stars formed in initially primordial star clusters should be found in the nearby Galactic halo regardless of the specifics of metal-free star formation. Observations of metal abundances in Galactic halo stars should be taken as directly constraining the properties of primordial stars, and the lack of metal-free halo stars today should be taken as strongly suggesting a 0.8 solar mass lower limit on the primordial initial mass function.

The thermal and fragmentation properties of star-forming clouds have important consequences on the corresponding characteristic stellar mass. The initial composition of the gas within these clouds is a record of the nucleosynthetic products of previous stellar generations. In this paper we present a model for the evolution of star-forming clouds enriched by metals and dust from the first supernovae, resulting from the explosions of metal-free progenitors with masses in the range 12 - 30 Msun and 140 - 260 Msun. Using a self-consistent approach, we show that: (i) metals depleted onto dust grains play a fundamental role, enabling fragmentation to solar or sub-solar mass scales already at metallicities Zcr = 10^{-6} Zsun; (ii) even at metallicities as high as 10^{-2} Zsun, metals diffused in the gas-phase lead to fragment mass scales which are ~ 100 Msun; (iii) C atoms are strongly depleted onto amorphous carbon grains and CO molecules so that CII plays a minor role in gas cooling, leaving OI as the main gas-phase cooling agent in low-metallicity clouds. These conclusions hold independently of the assumed supernova progenitors and suggest that the onset of low-mass star formation is conditioned to the presence of dust in the parent clouds.

We summarize some recent results on the metallicity and dust properties of Active Galactic Nuclei (AGN) at high redshift (1<z<6.4). By using the spectra of more than 5000 QSOs from the SDSS we find no evidence for any metallicity evolution in the redshift range 2<z<4.5, while there is a significant luminosity-metallicity dependence. These results are confirmed by the spectra of a smaller sample of narrow line AGNs at high-z (QSO2s and radio galaxies). The lack of metallicity evolution is interpreted both as a consequence of the cosmic downsizing and as a selection effect resulting from the joint QSO-galaxy evolution. The luminosity-metallicity relation is interpreted as a consequence of the mass-metallicity relation in the host galaxies of QSOs, but a relationship with the accretion rate is also possible. The lack of metallicity evolution is observed even in the spectra of the most distant QSOs known (z~6). This result is particularly surprising for elements such as Fe, C and Si, which are subject to a delayed enrichment, and requires that the hosts of these QSOs formed in short bursts and at very high redshift (z>10). The properties of dust in high-z QSOs are discussed within the context of the dust production mechanisms in the early universe. The dust extinction curve is observed to evolve beyond z>4, and by z~6 it is well described by the properties expected for dust produced by SNe, suggesting that the latter is the main mechanism of dust production in the early universe. We also show that the huge dust masses observed in distant QSOs can be accounted for by SN dust within the observational constraints currently available. Finally, we show that QSO winds, which have been proposed as an alternative mechanism of dust production, may also contribute significantly to the total dust budget at high redshift.

Kashlinsky et al. (2005) find a significant cosmic infrared background fluctuation excess on angular scales >50 arcsec that cannot be explained by instrumental noise or local foregrounds. The excess has been tentatively attributed to emission from primordial very massive (PopIII) stars formed <200 Myr after the Big Bang. Using an evolutionary model motivated by independent observations and including various feedback processes, we find that PopIII stars can contribute <40% of the total background intensity (\nu J_\nu ~ 1-2 nW m^-2 sr^-1 in the 0.8-8 \mum range) produced by all galaxies (hosting both PopIII and PopII stars) at z>5. The infrared fluctuation excess is instead very precisely accounted by the clustering signal of galaxies at z>5, predominantly hosting PopII stars with masses and properties similar to the present ones.

Motivated by theoretical predictions that first stars were predominantly very massive, we investigate the physics of the transition from an early epoch dominated by massive Pop III stars to a later epoch dominated by familiar low-mass Pop II/I stars by means of a numerically-generated catalogue of dark matter halos coupled with a self-consistent treatment of chemical and radiative feedback. Depending on the strength of the chemical feedback, Pop III stars can contribute a substantial fraction (several percent) of the cosmic star formation activity even at moderate redshifts, z = 5. We find that the three z = 10 sources tentatively detected in NICMOS UDFs should be powered by Pop III stars, if these are massive; however, this scenario fails to reproduce the derived WMAP electron scattering optical depth. Instead, both the UDFs and WMAP constraints can be fulfilled if stars at any time form with a more standard, slightly top-heavy, Larson IMF in the range 1 Msun < M < 100 Msun.

We performed R-band time series observations of the young, metal rich open cluster NGC 2301 for 12 nights in Feb. 2004. B images were also obtained and color magnitude diagrams, having limits of R=19.5 and B=21.5, were produced. Only asmall effort was made to determine cluster membership as our magnitude limits are far deeper than previously published values. Our photometric precision, for the brightest 5 magnitudes of sources, is 1-2 mmag. We determine that for the $\sim$4000 stars which have time-series data, 56% are variable and of these, approximately 13% are observed to exhibit periodic light curves ranging from tens of minutes to days. We present some examples of the light curves obtained and produce cuts in variability space based on parameters such as color and amplitude. The percentage of variability is approximately equal across all colors with the majority of variables having amplitudes of 0.15 magnitudes or less. In general, redder stars show larger variability amplitudes. We find a smooth decline in the number of periodic variables toward longer period. This decline is probably due to a transition from intrinsic to extrinsic variability and, in part, to our limited observing period of 12 nights. Essentially all the A and F main sequence stars in our sample are variable ($\sim$2 mmag and larger) and most present complex light curves containing multiple periods suggestive of their inclusion in the $\delta$Sct and $\gamma$Dor classes. A variable non-cluster member giant and two variable white dwarf candidates are discussed. Our equational description of variability is shown to be an excellent predictive tool for determining the cumulative fraction of variables that will be observed in a photometric survey. Our entire dataset is available electronically.

We study the detectability of primordial metal-free stars. Cosmological enrichment is a local process that takes place over an extended redshift range. While the duration of this transition depends on several unknown factors, in all cases late-forming metal-free stars are found in ~ 10^8 solar mass objects, which are just large enough to cool, but small enough to not be clustered near areas of previous star formation. We discuss the observational properties of these objects, some of which may have already been detected in ongoing Lyman-alpha surveys. If metal-free stars have masses between 140 and 260 solar masses, they are expected to end their lives as pair-production supernovae (PPSNe). Using the KEPLER code, we calculate a set of light curves for these objects. The peak luminosities of typical PPSNe are only slightly greater than those of Type Ia supernovae, but they remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova searches are able to place stringent limits on the fraction of very massive stars out to redshift 2. The planned Joint Dark Energy Mission satellite will be able to extend these constraints out to redshift 6

The thermal and chemical evolution of star-forming clouds is studied for different gas metallicities, Z, using the model of Omukai (2000), updated to include deuterium chemistry and the effects of cosmic microwave background (CMB) radiation. HD-line cooling dominates the thermal balance of clouds when Z \~ 10^{-5}-10^{-3} Z_sun and density ~10^{5} cm^{-3}. Early on, CMB radiation prevents the gas temperature to fall below T_CMB, although this hardly alters the cloud thermal evolution in low-metallicity gas. From the derived temperature evolution, we assess cloud/core fragmentation as a function of metallicity from linear perturbation theory, which requires that the core elongation E := (b-a)/a > E_NL ~ 1, where a (b) is the short (long) core axis length. The fragment mass is given by the thermal Jeans mass at E = E_NL. Given these assumptions and the initial (gaussian) distribution of E we compute the fragment mass distribution as a function of metallicity. We find that: (i) For Z=0, all fragments are very massive, > 10^{3}M_sun, consistently with previous studies; (ii) for Z>10^{-6} Z_sun a few clumps go through an additional high density (> 10^{10} cm^{-3}) fragmentation phase driven by dust-cooling, leading to low-mass fragments; (iii) The mass fraction in low-mass fragments is initially very small, but at Z ~ 10^{-5}Z_sun it becomes dominant and continues to grow as Z is increased; (iv) as a result of the two fragmentation modes, a bimodal mass distribution emerges in 0.01 < Z/Z_sun < 0.1. (v) For > 0.1Z_sun, the two peaks merge into a singly-peaked mass function which might be regarded as the precursor of the ordinary Salpeter-like IMF.

Interstellar dust plays a crucial role in the evolution of the Universe by assisting the formation of molecules, by triggering the formation of the first low-mass stars, and by absorbing stellar ultraviolet-optical light and subsequently re-emitting it at infrared/millimetre wavelengths. Dust is thought to be produced predominantly in the envelopes of evolved (age >1 Gyr), low-mass stars. This picture has, however, recently been brought into question by the discovery of large masses of dust in the host galaxies of quasars at redshift z>6, when the age of the Universe was less than 1 Gyr. Theoretical studies, corroborated by observations of nearby supernova remnants, have suggested that supernovae provide a fast and efficient dust formation environment in the early Universe. Here we report infrared observations of a quasar at redshift 6.2, which are used to obtain directly its dust extinction curve. We then show that such a curve is in excellent agreement with supernova dust models. This result demonstrates a supernova origin for dust in this high-redshift quasar, from which we infer that most of the dust at high redshifts has probably the same origin.

We show that the explosion of the first supernovae can trigger low-mass star formation via gravitational fragmentation of the supernova-driven gas shell. If the shell mass does not exceed the host galaxy gas mass, all explosions with energies E_SN > 10^{51} erg can lead to shell fragmentation. However, the minimum ambient density required to induce such fragmentation is much larger, n_0 > 300 cm^{-3}, for Type II supernovae than for pair-instability ones, which can induce star formation even at lower ambient densities. The typical mass of the unstable fragments is 10^{4-7} Msun; their density is in the range 110-6x10^7 cm{-3}. Fragments have a metallicity strictly lower than 10^{-2.6} Zsun and large values of the gravitational-to-pressure force ratio f ~ 8. Based on these findings, we conclude that the second generation of stars produced by such self-propagating star formation is predominantly constituted by low-mass, long-living, extremely metal-poor (or even metal-free, if mixing is suppressed) stars. We discuss the implications of such results for Pop III star formation scenarios and for the most iron-poor halo star HE0107-5240.

At redshift z>5 Type II supernovae (SNII) are the only known dust sources with evolutionary timescales shorter than the Hubble time. We extend the model of dust formation in the ejecta of SNII by Todini & Ferrara (2001) to investigate the same process in pair-instability supernovae (PISN), which are though to arise from the explosion of the first, metal free, very massive (140-260 Msun) cosmic stars. We find that 15%-30% of the PISN progenitor mass is converted into dust, a value >10 times higher than for SNII; PISN dust depletion factors (fraction of produced metals locked into dust grains) range between 0.3 and 0.7. These conclusions depend very weakly on the mass of the PISN stellar progenitor, which instead affects considerably the composition and size distribution. For the assumed temperature evolution, grain condensation starts 150-200 days after the explosion; the dominant compounds for all progenitor masses are SiO2 and Mg2SiO4 while the contribution of amorphous carbon and magnetite grains grows with progenitor mass; typical grain sizes range between 0.001 and a few 0.1 micron and are always smaller than 1 micron. We give a brief discussion of the implications of dust formation for the IMF evolution of the first stars, cosmic reionization and the intergalactic medium.

Studies of the broad emission-line regions (BLRs) in quasars have revealed solar or higher enrichment levels up to the highest redshifts. In combination with the presence of large amounts of dust in QSOs at $z$ about 6, this implies that substantial amounts of star formation and nucleosynthesis took place at significantly earlier epochs. Here, we examine whether a top-heavy stellar initial mass function (IMF) is indicated by current data, by modelling the contributions from different regions of the IMF, including Type Ia/II and pair instability supernovae, to the metal synthesis in BLRs. We find that, in order to reproduce the observations of roughly solar values of N/C and Fe/Mg in these objects, (i) stars with a present-day IMF are sufficient, regardless of their metallicity, (ii) zero-metallicity stars with a top-heavy IMF severely underproduce N/C, and (iii) the contribution of Type Ia SNe is not strongly required by the data. Therefore, stars of mass about 1-40 $M_\odot$ must have existed at $z$ about 10-20, possibly coeval with any hypothesized stars of masses exceeding about 100 $M_\odot$ at these epochs. This is in agreement with the nucleosynthetic abundance pattern detected in extremely metal-poor stars in the galactic halo.

Recent observations have shown that the intergalactic medium (IGM) is more transparent to Lyalpha photons close to Lyman Break Galaxies (LBGs) than at large distance from them, ie a proximity effect. Cosmological simulations including winds from LBGs have been so far unable to explain this trend. By coupling such simulations with the radiative transfer code CRASH, we investigate whether the addition of the ionizing radiation emitted by LBGs can increase the transmissivity by decreasing the neutral hydrogen fraction in the inner Mpc of the galaxy halo. The transmissivity as a function of distance is roughly reproduced only if LBGs are identified with dwarf galaxies (with masses < 10^9 solar masses), which are undergoing a vigorous (50 solar masses/yr) burst of star formation. Similar star formation rates in larger galaxies are not sufficient to overwhelm the large recombination rates associated with their denser environment. If so, photoionization partly reconciles theory with observations, although we discuss a number of uncertainties affecting both approaches.

The earliest stars to form in the Universe were the first sources of light, heat and metals after the Big Bang. The products of their evolution will have had a profound impact on subsequent generations of stars. Recent studies of primordial star formation have shown that, in the absence of metals (elements heavier than helium), the formation of stars with masses 100 times that of the Sun would have been strongly favoured, and that low-mass stars could not have formed before a minimum level of metal enrichment had been reached. The value of this minimum level is very uncertain, but is likely to be between 10^{-6} and 10^{-4} that of the Sun. Here we show that the recent discovery of the most iron-poor star known indicates the presence of dust in extremely low-metallicity gas, and that this dust is crucial for the formation of lower-mass second-generation stars that could survive until today. The dust provides a pathway for cooling the gas that leads to fragmentation of the precursor molecular cloud into smaller clumps, which become the lower-mass stars.

We conduct a comprehensive investigation of the detectability of the first stars and their enrichment signatures in galaxy clusters. We show that the mean metallicity of outflows from objects containing these Population III (PopIII) stars is well above the critical transition metallicity (Z_cr \sim 10^-4) that marks the formation of normal stars. Thus the fraction of PopIII objects formed as a function of redshift is heavily dependent on the distribution of metals and fairly independent of the precise value of Z_cr. Using an analytical model of inhomogenous structure formation, we study the evolution of PopIII objects as a function of the star formation efficiency, IMF, and efficiency of outflow generation. For all models, PopIII objects tend to be in the 10^6.5-10^7.0 solar mass range, just large enough to cool within a Hubble time, but small enough that they are not clustered near areas of previous star formation. Although the mean metallicity exceeds Z_cr at a redshift of 15 in all models, the peak of PopIII star formation occurs at z \sim 10, and such stars continue to form well into the observable range. We discuss the observational properties of these objects, some of which may have already been detected in ongoing surveys of high-redshift Lyman-alpha emitters. Finally, we combine our PopIII distributions with the yield models of Heger and Woosley (2002) to study their impact on the intracluster medium (ICM) in galaxy clusters. We find that PopIII stars can contribute no more than 20% of the iron observed in the ICM, although their peculiar elemental yields help to reconcile theoretical models with the observed Fe and Si/Fe abundances. However, these stars tend to overproduce S/Fe and their associated SN heating falls far short of the observed level of 1 keV per ICM gas particle.

Motivated by the relative lack of neutral hydrogen around Lyman Break Galaxies deduced from recent observations, we investigate the properties of the Lyalpha forest around high redshift galaxies. The study is based on improved numerical SPH simulations implementing, in addition to standard processes, a new scheme for multiphase and outflow physics description. Although on large scales our simulations reproduce a number of statistical properties of the IGM (because of the small filling factor of shock-heated gas), they underpredict the Lyalpha optical depth decrease inside 1 Mpc/h of the galaxies by a factor of ~2. We interpret this result as due to the combined effect of infall occurring along the filaments, which prevents efficient halo gas clearing by the outflow, and the insufficient increase of (collisional) hydrogen ionization produced by the temperature increase inside the hot, outflow-carved bubble. Unless an observational selection bias is present, we speculate that local photoionization could be the only viable explanation to solve the puzzle.

Recent studies suggest that the initial mass function (IMF) of the first stars was likely to be extremely top-heavy, unlike what is observed at present. We propose a scenario to generate fragmentation to lower masses once the first massive stars have formed and derive constraints on the primordial IMF. We estimate the mass fraction of pair-unstable supernovae, shown to be the dominant sources of the first heavy elements. These metals enrich the gas up to about $10^{-5}$ solar metallicity, when a transition to efficient cooling-driven fragmentation occurs producing 1 solar mass clumps. We argue that the remaining fraction of the first stars ends up in 100 solar mass VMBHs (Very Massive Black Holes). We obtain constraints on the fraction of first stars that contribute to the initial metal enrichment and the transition redshift for primordial IMF away from a top-heavy one, by making various assumptions about the fate of these VMBHs at late times. We conclude with a discussion of several hitherto unexplored implications of a high-mass dominated star formation mode in the early Universe.

If the first (PopIII) stars were very massive, their final fate is to collapse into very massive black holes. Once a proto-black hole has formed into the stellar core, accretion continues through a disk. It is widely accepted, although not confirmed, that magnetic fields drive an energetic jet which produces a burst of TeV neutrinos by photon-meson interaction, and eventually breaks out of the stellar envelope appearing as a Gamma Ray Burst (GRB). Based on recent numerical simulations and neutrino emission models, we predict the expected neutrino diffuse flux from these PopIII GRBs and compare it with the capabilities of present and planned detectors as AMANDA and IceCube. If beamed into 1% of the sky, we find that the rate of PopIII GRBs is $\le 4 \times 10^6$ yr$^{-1}$. High energy neutrinos from PopIII GRBs could dominate the overall flux in two energy bands [$10^4 - 10^5$] GeV and [$10^5 - 10^6$] GeV of neutrino telescopes. The enhanced sensitivities of forthcoming detectors in the high-energy band (AMANDA-II, IceCube) will provide a fundamental insight on the characteristic explosion energies of PopIII GRBs and will constitute a unique probe of the the Initial Mass Function (IMF) of the first stars and of the redshift $z_f$ marking the metallicity-driven transition from a top-heavy to a normal IMF. The current upper limit set by AMANDA-B10 implies that such transition must have occurred not later than $z_f =9.2$ for the most plausible jet energies. Based on such results, we speculate that PopIII GRBs, if not chocked, could be associated with a new class of events detected by BeppoSax, the Fast X-ray Transient (FXTs), which are bright X-ray sources, with peak energies in the 2-10 keV band and durations between 10-200 s.

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries as well as white dwarf-neutron star, white dwarf-black hole and black hole-neutron star systems. We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution. The gravitational signals emitted by each source during its early-inspiral phase add randomly to produce a stochastic background in the low frequency band with spectral strain amplitude between 10^{-18} Hz^{-1/2} and 5 10^{-17} Hz^{-1/2} at frequencies in the interval [5 10^{-6}-5 10^{-5}] Hz. The overall signal which, at frequencies above 10^{-4}Hz, is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range [1-10] mHz and acts like a confusion limited noise component which might limit the LISA sensitivity at frequencies above 1 mHz.

We study the gravitational wave emission from the first stars which are assumed to be Very Massive Objects (VMOs). We take into account various feedback (both radiative and stellar) effects regulating the collapse of objects in the early universe and thus derive the VMO initial mass function and formation rate. If the final fate of VMOs is to collapse, leaving very massive black hole remnants, then the gravitational waves emitted during each collapse would be seen as a stochastic background. The predicted spectral strain amplitude in a critical density Cold Dark Matter universe peaks in the frequency range \approx 5 \times 10^{-4}-5 \times 10^{-3} Hz where it has a value in the range \approx 10^{-20}-10^{-19} Hz^{-1/2}, and might be detected by LISA. The expected emission rate is roughly 4000 events/yr, resulting in a stationary, discrete sequence of bursts, i.e. a shot--noise signal.

Astrophysical sources of gravitational radiation are likely to have been formed since the beginning of star formation. Realistic source rates of formation throughout the Universe have been estimated from an observation-based determination of the star formation rate density evolution. Both the radiation emitted during the collapse to black holes and the spin-down radiation, induced by the r-mode instability, emitted by hot, young rapidly rotating neutron stars have been considered. We have investigated the overall signal produced by the ensemble of sources exploring the parameter space and discussing its possible detectability.

We estimate the spectral properties of the stochastic background of gravitational radiation emitted by a cosmological population of hot, young, rapidly rotating neutron stars. Their formation rate as a function of redshift is deduced from an observation-based determination of the star formation history in the Universe, and the gravitational energy is assumed to be radiated during the spin-down phase associated to the newly discovered r-mode instability. We calculate the overall signal produced by the ensemble of such neutron stars, assuming various cosmological backgrounds. We find that the spectral strain amplitude has a maximum $\approx (2-4)\times 10^{-26} {Hz}^{-1/2}$, at frequencies $\approx (30-60)$ Hz, while the corresponding closure density, $h^2 \Omega_{GW}$, has a maximum amplitude plateau of $\approx (2.2-3.3) \times 10^{-8}$ in the frequency range $(500-1700)$ Hz. We compare our results with a preliminary analysis done by Owen et al. (1998), and discuss the detectability of this background.

We analyse the stochastic background of gravitational radiation emitted by a cosmological population of core-collapse supernovae. The supernova rate as a function of redshift is deduced from an observation-based determination of the star formation rate density evolution. We then restrict our analysis to the range of progenitor masses leading to black hole collapse. In this case, the main features of the gravitational-wave emission spectra have been shown to be, to some extent, independent of the initial conditions and of the equation of state of the collapsing star, and to depend only on the black hole mass and angular momentum. We calculate the overall signal produced by the ensemble of black-hole collapses throughout the Universe, assuming a flat cosmology with vanishing cosmological constant. Within a wide range of parameter values, we find that the spectral strain amplitude has a maximum at a few hundred Hz with an amplitude between $10^{-28}$ and $10^{-27} Hz^{-1/2}$; the corresponding closure density, $\Omega_{GW}$, has a maximum amplitude ranging between $10^{-11}$ and $10^{-10}$ in the frequency interval $\sim 1.5-2.5$ kHz. Contrary to previous claims, our observation-based determination leads to a duty cycle of order 0.01, making our stochastic backgound a non-continuous one. Although the amplitude of our background is comparable to the sensitivity that can be reached by a pair of advanced LIGO detectors, the characteristic shot-noise structure of the predicted signal might be in principle exploited to design specific detection strategies.

The CMB anisotropy on arc second range is examined to test the power spectrum of perturbations in the small scale region and, in particular, to estimate the mass of dominant dark matter particles. It is shown that for the simplest evolutionary history with standard recombination, three and four beam observations could discriminate the mass of dark matter particles in the interval $0.5 KeV\leq M_{DM}\leq 4 KeV$ with an antenna beam (0.5 - 0.25) arc minute with amplitude $\approx 10^{-7}$.