ALMA observations have revealed the presence of dust in galaxies in the Epoch of Reionization (redshift $z>6$). However, the dust temperature, $T_d$, remains unconstrained, and this introduces large uncertainties, particularly in the dust mass determinations. Using an analytical and physically-motivated model, we show that dust in high-$z$, star-forming giant molecular clouds (GMC), largely dominating the observed far-infrared luminosity, is warmer ($T_d > 60\ \mathrm{K}$) than locally. This is due to the more compact GMC structure induced by the higher gas pressure and turbulence characterizing early galaxies. The compactness also delays GMC dispersal by stellar feedback, thus $\sim 40\%$ of the total UV radiation emitted by newly born stars remains obscured. A higher $T_d$ has additional implications: it (a) reduces the tension between local and high-$z$ IRX-$\beta$ relation, (b) alleviates the problem of the uncomfortably large dust masses deduced from observations of some EoR galaxies.

ALMA observations have revealed that [CII] 158$\mu$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($\Delta^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations.

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

Recent stacked ALMA observations have revealed that normal, star-forming galaxies at $z\approx 6$ are surrounded by extended ($\approx 10\,\mathrm{kpc}$) [CII] emitting halos which are not predicted by the most advanced, zoom-in simulations. We present a model in which these halos are the result of supernova-driven cooling outflows. Our model contains two free parameters, the outflow mass loading factor, $\eta$, and the parent galaxy dark matter halo circular velocity, $v_c$. The outflow model successfully matches the observed [CII] surface brightness profile if $\eta = 3.20 \pm 0.10$ and $v_c = 170 \pm 10{\,\rm km\,s^{-1}}$, corresponding to a dynamical mass of $\approx 10^{11}\, \mathrm{M}_\odot$. The predicted outflow rate and velocity range are $128 \pm 5 \,\mathrm{M}_\odot {\rm yr}^{-1}$ and $300-500 {\,\rm km\,s^{-1}}$, respectively. We conclude that: (a) extended halos can be produced by cooling outflows; (b) the large $\eta$ value is marginally consistent with starburst-driven outflows, but it might indicate additional energy input from AGN; (c) the presence of [CII] halos requires an ionizing photon escape fraction from galaxies $f_{\rm esc} \ll 1$. The model can be readily applied also to individual high-$z$ galaxies, as those observed, e.g., by the ALMA ALPINE survey now becoming available.

We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$\mu$m (and its known relation with star formation rate) and CIII]$\lambda$1909{\AA} emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average density, $n$, gas metallicity, $Z$, and burstiness parameter, $\kappa_s$, quantifying deviations from the Kennicutt-Schmidt relation. As an application, we consider COS-3018 (z=6.854), the only EoR Lyman Break Galaxy so far detected in both [CII] and CIII]. We show that COS-3018 is a moderate starburst ($\kappa_s \approx 3$), with $Z\approx 0.4\, Z_{\odot}$, and $n \approx 500\, {\rm cm^{-3}}$. Our method will be optimally applied to joint ALMA and JWST targets.

We present Atacama Large Millimiter/submillimiter Array (ALMA) observations of eight highly excited CO (J$_{\rm up}>8$) lines and continuum emission in two $z\sim6$ quasars: SDSS J231038.88+185519.7 (hereafter J2310), for which CO(8-7), CO(9-8), and CO(17-16) lines have been observed, and ULAS J131911.29+095951.4 (J1319), observed in the CO(14-13), CO(17-16) and CO(19-18) lines. The continuum emission of both quasars arises from a compact region ($< 0.9$ kpc). By assuming a modified black-body law, we estimate dust masses of Log$(M_{\rm dust}/M_{\odot})=8.75\pm0.07$ and Log$(M_{\rm dust}/M_{\odot})=8.8\pm0.2$ and dust temperatures of $T_{\rm dust}=76\pm3~{\rm K}$ and $T_{\rm dust}=66^{+15}_{-10}~{\rm K}$, respectively for J2310 and J1319. Only CO(8-7) and CO(9-8) in J2310 are detected, while $3\sigma$ upper limits on luminosities are reported for the other lines of both quasars. The CO line luminosities and upper limits measured in J2310 and J1319 are consistent with those observed in local AGN and starburst galaxies, and other $z\sim 6$ quasars, except for SDSS J1148+5251 (J1148), the only quasar at $z=6.4$ with a previous CO(17-16) line detection. By computing the CO SLEDs normalised to the CO(6-5) line and FIR luminosities for J2310, J1319, and J1149, we conclude that different gas heating mechanisms (X-ray radiation and/or shocks) may explain the different CO luminosities observed in these $z\sim6$ quasar. Future J$_{\rm up}>8$ CO observations will be crucial to understand the processes responsible for molecular gas excitation in luminous high-$z$ quasars.

We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their distribution in the galaxy disk, and the torus obscuration of the AGN radiation. We study the joint effect on the CO SLED excitation of (i) the X-ray luminosity from the AGN ($L_{X}$), (ii) the size of the molecular disk, (iii) the inclination angle ($\Omega$) of the torus with respect to the molecular disk, and (iv) the GMC mean density. We also discuss the possibility of an enhanced Cosmic Ray Ionization Rate (CRIR). We find that the X-ray Dominated Region (XDR) generated by the AGN in every case enhances the CO SLED for $J>5$, with increasing excitation of high-$J$ CO lines for increasing X-ray luminosity. Because high-$z$ galaxies are compact, the XDR region typically encloses the whole disk, thus its effect can be more important with respect to lower redshift objects. The impact of the XDR can be disentangled from an enhanced CRIR either if $L_X>10^{44} \rm \,erg\, s^{-1}$, or if $\Omega \geq 60^{\circ}$. We finally provide predictions on the CO(7-6)/[CII] and CO(17-16)/[CII] ratios as a function of $L_X$, which can be relevant for ALMA follow up of galaxies and quasars previously detected in [CII].

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

A tight relation between the [CII]158$\mu$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($\Sigma_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $\Sigma_{\rm CII}$ with the surface star formation rate ($\Sigma_*$). However, for $\Sigma_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $\Sigma_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($\kappa_s \gg 1$), parameterized by the "burstiness" parameter $\kappa_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$\mu$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies.

We study the photoevaporation of Jeans-unstable molecular clumps by isotropic FUV (6 eV $< {\rm h}\nu$ < 13.6 eV) radiation, through 3D radiative transfer hydrodynamical simulations implementing a non-equilibrium chemical network that includes the formation and dissociation of H$_2$. We run a set of simulations considering different clump masses ($M=10-200$ M$_\odot$) and impinging fluxes ($G_0=2\times 10^3-8\times 10^4$ in Habing units). In the initial phase, the radiation sweeps the clump as an R-type dissociation front, reducing the H$_2$ mass by a factor 40-90%. Then, a weak ($\mathcal{M}\simeq 2$) shock develops and travels towards the centre of the clump, which collapses while loosing mass from its surface. All considered clumps remain gravitationally unstable even if radiation rips off most of the clump mass, showing that external FUV radiation is not able to stop clump collapse. However, the FUV intensity regulates the final H$_2$ mass available for star formation: for example, for $G_0 < 10^4$ more than 10% of the initial clump mass survives. Finally, for massive clumps ($\sim 100$ M$_\odot$) the H$_2$ mass increases by 25-50% during the collapse, mostly because of the rapid density growth that implies a more efficient H$_2$ self-shielding.

The intrinsic strength of the Ly$\alpha$ line in young, star-forming systems makes it a special tool for studying high-redshift galaxies. However, interpreting observations remains challenging due to the complex radiative transfer involved. Here, we combine state-of-the-art hydrodynamical simulations of 'Althaea', a prototypical Lyman Break Galaxy (LBG, stellar mass $M_{\star}$$\simeq$ $10^{10}{\rm M}_{\odot})$ at $z=7.2$, with detailed radiative transfer computations of dust/continuum, [CII] 158 $\mu$m, and Ly$\alpha$ to clarify the relation between the galaxy properties and its Ly$\alpha$ emission. Althaea exhibits low ($f_\alpha< 1\%$) Ly$\alpha$ escape fractions and Equivalent Widths, EW $\lesssim 6$ Angstrom for the simulated lines of sight, with a large scatter. The correlation between escape fraction and inclination is weak, as a result of the rather chaotic structure of high-redshift galaxies. Low $f_\alpha$ values persist even if we artificially remove neutral gas around star forming regions to mimick the presence of HII regions. The high attenuation is primarily caused by dust clumps co-located with young stellar clusters. We can turn Althaea into a Lyman Alpha Emitter (LAE) only if we artificially remove dust from the clumps, yielding EWs up to $22$ Angstrom. Our study suggests that the LBG-LAE duty-cycle required by recent clustering measurements poses the challenging problem of a dynamically changing dust attenuation. Finally, we find an anti-correlation between the magnitude of Ly$\alpha$-[CII] line velocity shift and Ly$\alpha$ luminosity.

We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$\mu$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5\sigma$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects.

We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $\Sigma_{\rm [OIII]}/\Sigma_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations.

We report the discovery of Serenity-18, a galaxy at z=5.939 for which we could measure the content of molecular gas, M(H_2)~ 5 x10^9 M_sun, traced by the CO(6-5) emission, together with the metal-poor ([Fe/H]=-3.08 +- 0.12, [Si/H]=-2.86 +- 0.14) gas clump/filament which is possibly feeding its growth. The galaxy has an estimated star formation rate of ~100 M_sun yr^{-1}, implying that it is a typical main sequence galaxy at these redshifts. The metal-poor gas is detected through a damped Lyman-alpha absorber (DLA) observed at a spatial separation of 40 kpc and at the same redshift of Serenity-18, along the line of sight to the quasar SDSS J2310+1855 (z_em = 6.0025). The chemical abundances measured for the damped Lyman-alpha system are in very good agreement with those measured for other DLAs discovered at similar redshifts, indicating an enrichment due to massive PopII stars. The galaxy/Damped system we discovered is a direct observational evidence of the assembly of a galaxy at the edge of the reionization epoch.

Dwarf galaxies are ideal laboratories to study the physics of the interstellar medium (ISM). Emission lines have been widely used to this aim. Retrieving the full information encoded in the spectra is therefore essential. This can be efficiently and reliably done using Machine Learning (ML) algorithms. Here, we apply the ML code GAME to MUSE (Multi Unit Spectroscopic Explorer) and PMAS (Potsdam Multi Aperture Spectrophotometer) Integral Field Unit (IFU) observations of two nearby Blue Compact Galaxies (BCGs): Henize 2-10 and IZw18. We derive spatially resolved maps of several key ISM physical properties. We find that both galaxies show a remarkably uniform metallicity distribution. Henize 2-10 is a star forming dominated galaxy, with a Star Formation Rate (SFR) of about 1.2 M$_{\odot}$ yr$^{-1}$. Henize 2-10 features dense and dusty ($A_V$ up to 5-7 mag) star forming central sites. We find IZw18 to be very metal-poor ($Z$= 1/20 $Z_{\odot}$). IZw18 has a strong interstellar radiation field, with a large ionization parameter. We also use models of PopIII stars spectral energy distribution as a possible ionizing source for the HeII $\lambda$4686 emission detected in the IZw18 NW component. We find that PopIII stars could provide a significant contribution to the line intensity. The upper limit to the PopIII star formation is 52% of the total IZw18 SFR.

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 present ALMA observations of the CO(6-5) and [CII] emission lines and the sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic beam that is ten times smaller in angular size, we have achieved ten times better sensitivity in the CO(6-5) line, and two and half times better sensitivity in the [CII] line, enabling us to resolve the molecular gas emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO} = 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of $\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the observed CO velocity gradient is due to an inclined rotating disk, we derive a dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~ M_{\odot}$, which is a factor of approximately two smaller than the previously reported estimate based on [CII]. Regarding the central black hole, we provide a new estimate of the black hole mass based on the C~IV emission line detected in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~ M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$, where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio $v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows and/or complex kinematics due to a merger event. We estimate a global Toomre parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare, at the same angular resolution, the CO(6-5) and [CII] distributions, finding that dense molecular gas is more centrally concentrated with respect to [CII]. We find that the current BH growth rate is similar to that of its host galaxy.

We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.

The nature of the seeds of the observed high-z super-massive black holes (SMBH) is unknown. Although different options have been proposed, involving e.g. intermediate mass direct collapse black holes, BH remnants of massive stars remain the most natural explanation. To identify the most favorable conditions (if any) for their rapid growth, we study the accretion rate of a M_BH = 100M_sun BH formed in a typical z = 10 galaxy under different conditions (e.g. galaxy structure, BH initial position and velocity). We model the galaxy baryonic content and follow the BH orbit and accretion history for 300 Myr (the time span in 10 > z > 7), assuming the radiation-regulated accretion model by Park & Ricotti (2013). We find that, within the limits of our model, BH seeds cannot grow by more than 30 percent, suggesting that accretion on light-seed models are inadequate to explain high-z SMBH. We also compute the X-ray emission from such accreting stellar BH population in the (0.5, 8) keV band and find it comparable to the one produced by high-mass X-ray binaries. This study suggests that early BHs, by X-ray pre-heating of the intergalactic medium at cosmic dawn, might leave a specific signature on the HI 21 cm line power spectrum potentially detectable with SKA.

We present the analysis of the rest-frame optical-to-UV spectrum of APM 08279+5255, a well-known lensed broad absorption line (BAL) quasar at $z = 3.911$. The spectroscopic data are taken with the optical DOLoRes and near-IR NICS instruments at TNG, and include the previously unexplored range between C III] $\lambda$1910 and [O III] $\lambda\lambda$4959,5007. We investigate the possible presence of multiple BALs by computing "balnicity" and absorption indexes (i.e. BI, BI$_0$ and AI) for the transitions Si IV $\lambda$1400, C IV $\lambda$1549, Al III $\lambda$1860 and Mg II $\lambda$2800. No clear evidence for the presence of absorption features is found in addition to the already known, prominent BAL associated to C IV, which supports a high-ionization BAL classification for APM 08279+5255. We also study the properties of the [O III], H$\beta$ and Mg II emission lines. We find that [O III] is intrinsically weak ($F_{\rm [OIII]}/F_{\rm H\beta} \lesssim 0.04$), as it is typically found in luminous quasars with a strongly blueshifted C IV emission line ($\sim$2500 km s$^{-1}$ for APM 08279+5255). We compute the single-epoch black hole mass based on Mg II and H$\beta$ broad emission lines, finding $M_{\rm BH} = (2 \div 3) \times 10^{10}\mu^{-1}$ M$_\odot$, with the magnification factor $\mu$ that can vary between 4 and 100 according to CO and rest-frame UV-to-mid-IR imaging respectively. Using a Mg II equivalent width (EW)-to-Eddington ratio relation, the EW$_{\rm MgII} \sim 27$ \AA\ measured for APM 08279+5255 translates into an Eddington ratio of $\sim$0.4, which is more consistent with $\mu=4$. This magnification factor also provides a value of $M_{\rm BH}$ that is consistent with recent reverberation-mapping measurements derived from C IV and Si IV.

We present an updated, optimized version of GAME (GAlaxy Machine learning for Emission lines), a code designed to infer key interstellar medium physical properties from emission line intensities of UV/optical/far infrared galaxy spectra. The improvements concern: (a) an enlarged spectral library including Pop III stars; (b) the inclusion of spectral noise in the training procedure, and (c) an accurate evaluation of uncertainties. We extensively validate the optimized code and compare its performance against empirical methods and other available emission line codes (PYQZ and HII-CHI-MISTRY) on a sample of 62 SDSS stacked galaxy spectra and 75 observed HII regions. Very good agreement is found for metallicity. However, ionization parameters derived by GAME tend to be higher. We show that this is due to the use of too limited libraries in the other codes. The main advantages of GAME are the simultaneous use of all the measured spectral lines, and the extremely short computational times. We finally discuss the code potential and limitations.

The recent discovery of dusty galaxies well into the Epoch of Reionization (redshift $z>6$) poses challenging questions about the properties of the interstellar medium in these pristine systems. By combining state-of-the-art hydrodynamic and dust radiative transfer simulations, we address these questions focusing on the recently discovered dusty galaxy A2744_YD4 ($z=8.38$, Laporte et al. 2017}). We show that we can reproduce the observed spectral energy distribution (SED) only using different physical values with respect to the inferred ones by Laporte et al(2017), i.e. a star formation rate of $\mathrm{SFR} = 78\;\rm M_\odot \rm yr^{-1}$, a factor $\approx 4$ higher than deduced from simple Spectral Energy Distribution fitting. In this case we find: (a) dust attenuation (corresponding to $\tau_V=1.4$) is consistent with a Milky Way extinction curve; (b) the dust-to-metal ratio is low, $f_\mathrm{d} \sim 0.08$, implying that early dust formation is rather inefficient; (c) the luminosity-weighted dust temperature is high, $T_d=91\pm 23\, \rm K$, as a result of the intense ($\approx 100\times$ MW) interstellar radiation field; (d) due to the high $T_d$, the ALMA Band 7 detection can be explained by a limited dust mass, $M_d=1.6\times 10^6 $M$_\odot$. Finally, the high dust temperatures might solve the puzzling low infrared excess recently deduced for high-$z$ galaxies from the IRX-$\beta$ relation.

We present the first attempt to detect outflows from galaxies approaching the Epoch of Reionization (EoR) using a sample of 9 star-forming ($\rm SFR=31\pm 20~M_{\odot}~yr^{-1}$) $z\sim 5.5$ galaxies for which the [CII]158$\mu$m line has been previously obtained with ALMA. We first fit each line with a Gaussian function and compute the residuals by subtracting the best fitting model from the data. We combine the residuals of all sample galaxies and find that the total signal is characterised by a flux excess of $\sim 0.5$ mJy extended over $\sim 1000$ km~s$^{-1}$. Although we cannot exclude that part of this signal is due to emission from faint satellite galaxies, we show that the most probable explanation for the detected flux excess is the presence of broad wings in the [CII] lines, signatures of starburst-driven outflows. We infer an average outflow rate of $\rm \dot{M}=54\pm23~ M_{\odot}~yr^{-1}$, providing a loading factor $\eta=\rm \dot{M}/SFR=1.7\pm1.3$ in agreement with observed local starbursts. Our interpretation is consistent with outcomes from zoomed hydro-simulations of {\it Dahlia}, a $z\sim 6$ galaxy ($\rm SFR\sim 100~\rm M_{\odot}~yr^{-1}$) whose feedback-regulated star formation results into an outflow rate $\rm \dot{M}\sim 30~ M_{\odot}~yr^{-1}$. The quality of the ALMA data is not sufficient for a detailed analysis of the [CII] line profile in individual galaxies. Nevertheless, our results suggest that starburst-driven outflows are in place in the EoR and provide useful indications for future ALMA campaigns. Deeper observations of the [CII] line in this sample are required to better characterise feedback at high-$z$ and to understand the role of outflows in shaping early galaxy formation.

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.

We have performed a high sensitivity observation of the UFO/BAL quasar APM 08279+5255 at z=3.912 with NOEMA at 3.2 mm, aimed at detecting fast moving molecular gas. We report the detection of blueshifted CO(4-3) with maximum velocity (v95\%) of $-1340$ km s$^{-1}$, with respect to the systemic peak emission, and a luminosity of $L' = 9.9\times 10^9 ~\mu^{-1}$ K km s$^{-1}$ pc$^{-2}$ (where $\mu$ is the lensing magnification factor). We discuss various scenarios for the nature of this emission, and conclude that this is the first detection of fast molecular gas at redshift $>3$. We derive a mass flow rate of molecular gas in the range $\rm \dot M=3-7.4\times 10^3$ M$_\odot$/yr, and momentum boost $\dot P_{OF} / \dot P_{AGN} \sim 2-6$, therefore consistent with a momentum conserving flow. For the largest $\dot P_{OF}$ the scaling is also consistent with a energy conserving flow with an efficiency of $\sim$10-20\%. The present data can hardly discriminate between the two expansion modes. The mass loading factor of the molecular outflow $\eta=\dot M_{OF}/SFR$ is $>>1$. We also detect a molecular emission line at a frequency of 94.83 GHz, corresponding to a rest frame frequency of 465.8 GHz, which we tentatively identified with the cation molecule $\rm N_2H^+$(5-4), which would be the first detection of this species at high redshift. We discuss the alternative possibility that this emission is due to a CO emission line from the, so far undetected, lens galaxy. Further observations of additional transitions of the same species with NOEMA can discriminate between the two scenarios.

We study the CO line luminosity ($L_{\rm CO}$), the shape of the CO Spectral Line Energy Distribution (SLED), and the value of the CO-to-$\rm H_2$ conversion factor in galaxies in the Epoch of Reionization (EoR). To this aim, we construct a model that simultaneously takes into account the radiative transfer and the clumpy structure of giant molecular clouds (GMCs) where the CO lines are excited. We then use it to post-process state-of-the-art zoomed, high resolution ($30\, \rm{pc}$), cosmological simulation of a main-sequence ($M_{*}\approx10^{10}\, \rm{M_{\odot}}$, $SFR\approx 100\,\rm{M_{\odot}\, yr^{-1}}$) galaxy, "Alth{\ae}a", at $z\approx6$. We find that the CO emission traces the inner molecular disk ($r\approx 0.5 \,\rm{kpc}$) of Alth{\ae}a with the peak of the CO surface brightness co-located with that of the [CII] 158$\rm \mu m$ emission. Its $L_{\rm CO(1-0)}=10^{4.85}\, \rm{L_{\odot}}$ is comparable to that observed in local galaxies with similar stellar mass. The high ($\Sigma_{gas} \approx 220\, \rm M_{\odot}\, pc^{-2}$) gas surface density in Alth{\ae}a, its large Mach number (\mach$\approx 30$), and the warm kinetic temperature ($T_{k}\approx 45 \, \rm K$) of GMCs yield a CO SLED peaked at the CO(7-6) transition, i.e. at relatively high-$J$, and a CO-to-$\rm H_2$ conversion factor $\alpha_{\rm CO}\approx 1.5 \, \rm M_{\odot} \rm (K\, km\, s^{-1}\, pc^2)^{-1} $ lower than that of the Milky Way. The ALMA observing time required to detect (resolve) at 5$\sigma$ the CO(7-6) line from galaxies similar to Alth{\ae}a is $\approx13$ h ($\approx 38$ h).

We study the photoevaporation of molecular clumps exposed to a UV radiation field including hydrogen-ionizing photons ($h\nu > 13.6$ eV) produced by massive stars or quasars. We follow the propagation and collision of shock waves inside clumps and take into account self-shielding effects, determining the evolution of clump size and density with time. The structure of the ionization-photodissociation region (iPDR) is obtained for different initial clump masses ($M=0.01 - 10^4\,{\rm M}_\odot$) and impinging fluxes ($G_0=10^2 - 10^5$ in units of the Habing flux). The cases of molecular clumps engulfed in the HII region of an OB star and clumps carried within quasar outflows are treated separately. We find that the clump undergoes in both cases an initial shock-contraction phase and a following expansion phase, which lets the radiation penetrate in until the clump is completely evaporated. Typical evaporation time-scales are $\simeq 0.01$ Myr in the stellar case and 0.1 Myr in the quasar case, where the clump mass is 0.1 ${\rm M}_\odot$ and $10^3\,{\rm M}_\odot$ respectively. We find that clump lifetimes in quasar outflows are compatible with their observed extension, suggesting that photoevaporation is the main mechanism regulating the size of molecular outflows.

To improve our understanding of high-z galaxies we study the impact of H$_{2}$ chemistry on their evolution, morphology and observed properties. We compare two zoom-in high-resolution (30 pc) simulations of prototypical $M_{\star}\sim 10^{10} {\rm M}_{\odot}$ galaxies at $z=6$. The first, "Dahlia", adopts an equilibrium model for H$_{2}$ formation, while the second, "Alth{\ae}a", features an improved non-equilibrium chemistry network. The star formation rate (SFR) of the two galaxies is similar (within 50\%), and increases with time reaching values close to 100 ${\rm M}_{\odot}/\rm yr$ at $z=6$. They both have SFR-stellar mass relation consistent with observations, and a specific SFR of $\simeq 5\, {\rm Gyr}^{-1}$. The main differences arise in the gas properties. The non-equilibrium chemistry determines the H$\rightarrow$ H$_{2}$~transition to occur at densities $> 300\,{cm}^{-3}$, i.e. about 10 times larger than predicted by the equilibrium model used for Dahlia. As a result, Alth{\ae}a features a more clumpy and fragmented morphology, in turn making SN feedback more effective. Also, because of the lower density and weaker feedback, Dahlia sits $3\sigma$ away from the Schmidt-Kennicutt relation; Alth{\ae}a, instead nicely agrees with observations. The different gas properties result in widely different observables. Alth{\ae}a outshines Dahlia by a factor of 7 (15) in [CII]~$157.74\,\mu{\rm m}$ (H$_{2}$~$17.03\,\mu{\rm m}$) line emission. Yet, Alth{\ae}a is under-luminous with respect to the locally observed [CII]-SFR relation. Whether this relation does not apply at high-z or the line luminosity is reduced by CMB and metallicity effects remains as an open question.

We investigate quasar outflows at $z \geq 6$ by performing zoom-in cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we zoom in the $2 R_{200}$ region around a $2 \times 10^{12} M_{\odot}$ halo at $z = 6$, inside a $(500 ~ {\rm Mpc})^3$ comoving volume. We compare the results of our AGN runs with a control simulation in which only stellar/SN feedback is considered. Seeding $10^5 M_{\odot}$ BHs at the centers of $10^{9} M_{\odot}$ halos, we find the following results. BHs accrete gas at the Eddington rate over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{\rm BH} = 4 \times 10^9 M_{\odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($\dot{M}_{\rm out} \sim 2000 M_{\odot}$/yr) outflows of shock-heated low-density gas form at $z \sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over $z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is reduced by a factor of $5$ ($1000$). We analyse the relative contribution of multiple physical process: (i) disrupting cosmic filamentary cold gas inflows, (ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and find that AGN feedback has the following effects at $z = 6$. The inflowing gas mass fraction is reduced by $\sim 12 \%$, the high-density gas fraction is lowered by $\sim 13 \%$, and $\sim 20 \%$ of the gas outflows at a speed larger than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at $z \geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN feedback quenches their star formation dominantly by powerful outflows ejecting gas out of the host galaxy halo.

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.

We present new ALMA observations of the [OIII]88$\mu$m line and high angular resolution observations of the [CII]158$\mu$m line in a normal star forming galaxy at z$=$7.1. Previous [CII] observations of this galaxy had detected [CII] emission consistent with the Ly$\alpha$ redshift but spatially slightly offset relative to the optical (UV-rest frame) emission. The new [CII] observations reveal that the [CII] emission is partly clumpy and partly diffuse on scales larger than about 1kpc. [OIII] emission is also detected at high significance, offset relative to the optical counterpart in the same direction as the [CII] clumps, but mostly not overlapping with the bulk of the [CII] emission. The offset between different emission components (optical/UV and different far-IR tracers) is similar to what observed in much more powerful starbursts at high redshift. We show that the [OIII] emitting clump cannot be explained in terms of diffuse gas excited by the UV radiation emitted by the optical galaxy, but it requires excitation by in-situ (slightly dust obscured) star formation, at a rate of about 7 M$_{\odot}$/yr. Within 20 kpc from the optical galaxy the ALMA data reveal two additional [OIII] emitting systems, which must be star forming companions. We discuss that the complex properties revealed by ALMA in the z$\sim$7.1 galaxy are consistent with expectations by recent models and cosmological simulations, in which differential dust extinction, differential excitation and different metal enrichment levels, associated with different subsystems assembling a galaxy, are responsible for the different appearance of the system when observed with different tracers.

Since the beginning of the new millennium, more than 100 $z\sim 6$ quasars have been discovered through several surveys and followed-up with multi-wavelength observations. These data provided a large amount of information on the growth of supermassive black holes at the early epochs, the properties of quasar host galaxies and the joint formation and evolution of these massive systems. We review the properties of the highest-$z$ quasars known so far, especially focusing on some of the most recent results obtained in (sub-)millimeter bands. We discuss key observational challenges and open issues in theoretical models and highlight possible new strategies to improve our understanding of the galaxy-black hole formation and evolution in the early Universe.

We present the 78-ks Chandra observations of the $z=6.4$ quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance $\gtrsim9\sigma$). The X-ray spectrum is best-fitted by a power-law with photon index $\Gamma=1.9$ absorbed by a gas column density of $\rm N_{\rm H}=2.0^{+2.0}_{-1.5}\times10^{23}\,\rm cm^{-2}$. We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to $\sim 1.5\times 10^{45}~\rm erg~s^{-1}$, comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value ($\alpha_{\rm OX}=-1.76\pm 0.14$) of J1148 is consistent with the one found in quasars with similar rest-frame 2500 \AA ~luminosity ($L_{\rm 2500}\sim 10^{32}~\rm erg~s^{-1}$\AA$^{-1}$). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich ($Z\geq \rm Z_{\odot}$) molecular clouds distributed on $\sim$kpc scales in the host galaxy. Our analysis favors a black hole mass $M_{\rm BH} \sim 3\times 10^9 \rm M_\odot$ and a molecular hydrogen mass $M_{\rm H_2}\sim 2\times 10^{10} \rm M_\odot$, in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis.

The peculiar emission properties of the $z \sim 6.6$ Ly$\alpha$ emitter CR7 have been initially interpreted with the presence of either a direct collapse black hole (DCBH) or a substantial mass of Pop III stars. Instead, updated photometric observations by Bowler et al. (2016) seem to suggest that CR7 is a more standard system. Here we confirm that the original DCBH hypothesis is consistent also with the new data. Using radiation-hydrodynamic simulations, we reproduce the new IR photometry with two models involving a Compton-thick DCBH of mass $\approx 7 \times 10^6 \, \mathrm{M_{\odot}}$ accreting (a) metal-free ($Z=0$) gas with column density $N_H = 8 \times 10^{25} \, \mathrm{cm^{-2}}$, or (b) low-metallicity gas ($Z = 5 \times 10^{-3} \, \mathrm{Z_{\odot}}$) with $N_H = 3 \times 10^{24} \, \mathrm{cm^{-2}}$. The best fit model reproduces the photometric data to within $1 \sigma$. Such metals can be produced by weak star-forming activity occurring after the formation of the DCBH. The main contribution to the Spitzer/IRAC $3.6 \, \mathrm{\mu m}$ photometric band in both models is due to HeI/HeII $\lambda 4714, 4687$ emission lines, while the contribution of [OIII] $\lambda 4959, 5007$ emission lines, if present, is sub-dominant. Spectroscopic observations with JWST will be required to ultimately clarify the nature of CR7.

With the aim of improving predictions on far infrared (FIR) line emission from Giant Molecular Clouds (GMC), we study the effects of photoevaporation (PE) produced by external far-ultraviolet (FUV) and ionizing (extreme-ultraviolet, EUV) radiation on GMC structure. We consider three different GMCs with mass in the range $M_{\rm GMC} = 10^{3-6}\,\rm{M_{\odot}}$. Our model includes: (i) an observationally-based inhomogeneous GMC density field, and (ii) its time evolution during the PE process. In the fiducial case ($M_{\rm GMC}\approx10^5 M_{\odot}$), the photoevaporation time ($t_{pe}$) increases from 1 Myr to 30 Myr for gas metallicity $Z=0.05-1\,\rm Z_{\odot}$, respectively. Next, we compute the time-dependent luminosity of key FIR lines tracing the neutral and ionized gas layers of the GMCs, ([CII] at $158\,\rm{\mu m}$, [OIII] at $88\,\rm \mu m$) as a function of $G_0$, and $Z$ until complete photoevaporation at $t_{pe}$. We find that the specific [CII] luminosity is almost independent on the GMC model within the survival time of the cloud. Stronger FUV fluxes produce higher [CII] and [OIII] luminosities, however lasting for progressively shorter times. At $Z=Z_{\odot}$ the [CII] emission is maximized ($L_{\rm CII}\approx 10^4\,\rm{L_{\odot}}$ for the fiducial model) for $t<1\,\rm{Myr}$ and $\log G_0\geq 3$. Noticeably, and consistently with the recent detection by Inoue et al. (2016) of a galaxy at redshift $z\approx 7.2$, for $Z\leq 0.2\,\rm{Z_{\odot}}$ the [OIII] line might outshine [CII] emission by up to $\approx 1000$ times. We conclude that the [OIII] line is a key diagnostic of low metallicity ISM, especially in galaxies with very young stellar populations.

We introduce an intrinsic Ly\alpha\ emission line profile reconstruction method for high-$z$ quasars (QSOs). This approach utilises a covariance matrix of emission line properties obtained from a large, moderate-$z$ ($2 \leq z \leq 2.5$), high signal to noise (S/N > 15) sample of BOSS QSOs. For each QSO, we complete a Monte Carlo Markov Chain fitting of the continuum and emission line properties and perform a visual quality assessment to construct a large database of robustly fit spectra. With this dataset, we construct a covariance matrix to describe the correlations between the high ionisation emission lines Ly\alpha, C IV, Si IV + O IV] and C III], and find it to be well approximated by an $N$-dimensional Gaussian distribution. This covariance matrix characterises the correlations between the line width, peak height and velocity offset from systemic while also allowing for the existence of broad and narrow line components for Ly\alpha\ and C IV. We illustrate how this covariance matrix allows us to statistically characterise the intrinsic Ly\alpha\ line solely from the observed spectrum redward of 1275\AA. This procedure can be used to reconstruct the intrinsic Ly\alpha\ line emission profile in cases where Ly\alpha\ may otherwise be obscured. Applying this reconstruction method to our sample of QSOs, we recovered the Ly\alpha\ line flux to within 15 per cent of the measured flux at 1205\AA\ (1220\AA) ~85 (90) per cent of the time.

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

We study the redshift evolution of the quasar UV Luminosity Function (LF) for 0.5 < z < 6.5, by collecting the most up to date observational data and, in particular, the recently discovered population of faint AGNs. We fit the QSO LF using either a double power-law or a Schechter function, finding that both forms provide good fits to the data. We derive empirical relations for the LF parameters as a function of redshift and, based on these results, predict the quasar UV LF at z=8. From the inferred LF evolution, we compute the redshift evolution of the QSO/AGN comoving ionizing emissivity and hydrogen photoionization rate. If faint AGNs are included, the contribution of quasars to reionization increases substantially. However, their level of contribution critically depends on the detailed shape of the QSO LF, which can be constrained by efficient searches of high-z quasars. To this aim, we predict the expected (i) number of z>6 quasars detectable by ongoing and future NIR surveys (as EUCLID and WFIRST), and (ii) number counts for a single radio-recombination line observation with SKA-MID (FoV = 0.49 deg^2) as a function of the Hnalpha flux density, at 0<z<8. These surveys (even at z<6) will be fundamental to better constrain the role of quasars as reionization sources.

We present a new approach based on Supervised Machine Learning (SML) algorithms to infer key physical properties of galaxies (density, metallicity, column density and ionization parameter) from their emission line spectra. We introduce a numerical code (called GAME, GAlaxy Machine learning for Emission lines) implementing this method and test it extensively. GAME delivers excellent predictive performances, especially for estimates of metallicity and column densities. We compare GAME with the most widely used diagnostics (e.g. R$_{23}$, [NII]$\lambda$6584 / H$\alpha$ indicators) showing that it provides much better accuracy and wider applicability range. GAME is particularly suitable for use in combination with Integral Field Unit (IFU) spectroscopy, both for rest-frame optical/UV nebular lines and far-infrared/sub-mm lines arising from Photo-Dissociation Regions. Finally, GAME can also be applied to the analysis of synthetic galaxy maps built from numerical simulations.

We present zoom-in, AMR, high-resolution ($\simeq 30$ pc) simulations of high-redshift ($z \simeq 6$) galaxies with the aim of characterizing their internal properties and interstellar medium. Among other features, we adopt a star formation model based on a physically-sound molecular hydrogen prescription, and introduce a novel scheme for supernova feedback, stellar winds and dust-mediated radiation pressure. In the zoom-in simulation the target halo hosts "Dahlia", a galaxy with a stellar mass $M_*=1.6\times 10^{10}$M$_\odot$, representative of a typical $z\sim 6$ Lyman Break Galaxy. Dahlia has a total H2 mass of $10^{8.5}$M$_\odot$, that is mainly concentrated in a disk-like structure of effective radius $\simeq 0.6$ kpc and scale height $\simeq 200$ pc. Frequent mergers drive fresh gas towards the center of the disk, sustaining a star formation rate per unit area of $\simeq 15 $M$_\odot$ yr$^{-1}$ kpc$^{-2}$. The disk is composed by dense ($n \gtrsim 25$ cm$^{-3}$), metal-rich ($Z \simeq 0.5 $ Z$_\odot$) gas, that is pressure-supported by radiation. We compute the $158\mu$m [CII] emission arising from {Dahlia}, and find that $\simeq 95\%$ of the total [CII] luminosity ($L_{[CII]}\simeq10^{7.5}$ L$_\odot$) arises from the H2 disk. Although $30\%$ of the CII mass is transported out of the disk by outflows, such gas negligibly contributes to [CII] emission, due to its low density ($n \lesssim 10$ cm$^{-3}$) and metallicity ($Z\lesssim 10^{-1}$Z$_\odot$). Dahlia is under-luminous with respect to the local [CII]-SFR relation; however, its luminosity is consistent with upper limits derived for most $z\sim6$ galaxies.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

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

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$\mu$Jy, to derive differential number counts down to 60$\mu$Jy and 100$\mu$Jy at $\lambda=$1.3 mm and $\lambda=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$\mu$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{\lambda = 1.1~mm} = 130~\mu$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{\lambda = 1.3~mm} = 60~\mu$Jy. At the faintest flux limits, i.e. 30$\mu$Jy and 40$\mu$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

After two ALMA observing cycles, only a handful of [CII] $158\,\mu m$ emission line searches in z>6 galaxies have reported a positive detection, questioning the applicability of the local [CII]-SFR relation to high-z systems. To investigate this issue we use the Vallini et al. 2013 (V13) model, based on high-resolution, radiative transfer cosmological simulations to predict the [CII] emission from the interstellar medium of a z~7 (halo mass $M_h=1.17\times10^{11}M_{\odot}$) galaxy. We improve the V13 model by including (a) a physically-motivated metallicity (Z) distribution of the gas, (b) the contribution of Photo-Dissociation Regions (PDRs), (c) the effects of Cosmic Microwave Background on the [CII] line luminosity. We study the relative contribution of diffuse neutral gas to the total [CII] emission ($F _{diff}/F_{tot}$) for different SFR and Z values. We find that the [CII] emission arises predominantly from PDRs: regardless of the galaxy properties, $F _{diff}/F_{tot}\leq 10$% since, at these early epochs, the CMB temperature approaches the spin temperature of the [CII] transition in the cold neutral medium ($T_{CMB}\sim T_s^{CNM}\sim 20$ K). Our model predicts a high-z [CII]-SFR relation consistent with observations of local dwarf galaxies ($0.02<Z/Z_{\odot}<0.5$). The [CII] deficit suggested by actual data ($L_{CII}<2.0\times 10^7 L_{\odot}$ in BDF3299 at z~7.1) if confirmed by deeper ALMA observations, can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. The deviation from the local [CII]-SFR would then imply a modified Kennicutt-Schmidt relation in z>6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities ($Z<0.1 Z_{\odot}$).

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $\mu$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($\Delta I/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($\nu_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/\mu{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4\sigma$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,\mu$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

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.

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

The intensity mapping of the [CII] 157.7 $\rm \mu$m fine-structure emission line represents an ideal experiment to probe star formation activity in galaxies, especially in those that are too faint to be individually detected. Here, we investigate the feasibility of such an experiment for $z > 5$ galaxies. We construct the $L_{\rm CII} - M_{\rm h}$ relation from observations and simulations, then generate mock [CII] intensity maps by applying this relation to halo catalogs built from large scale N-body simulations. Maps of the extragalactic far-infrared (FIR) continuum, referred to as "foreground", and CO rotational transition lines and [CI] fine-structure lines referred to as "contamination", are produced as well. We find that, at 316 GHz (corresponding to $z_{\rm CII} = 5$), the mean intensities of the extragalactic FIR continuum, [CII] signal, all CO lines from $J=1$ to 13 and two [CI] lines are $\sim 3\times10^5$ Jy sr$^{-1}$, $\sim 1200$ Jy sr$^{-1}$, $\sim 800$ Jy sr$^{-1}$ and $\sim 100$ Jy sr$^{-1}$, respectively. We discuss a method that allows us to subtract the FIR continuum foreground by removing a spectrally smooth component from each line of sight, and to suppress the CO/[CI] contamination by discarding pixels that are bright in contamination emission. The $z > 5$ [CII] signal comes mainly from halos in the mass range $10^{11-12} \,M_\odot$; as this mass range is narrow, intensity mapping is an ideal experiment to investigate these early galaxies. In principle such signal is accessible to a ground-based telescope with a 6 m aperture, 150 K system temperature, a $128\times128$ pixels FIR camera in 5000 hr total integration time, however it is difficult to perform such an experiment by using currently available telescopes.

Past observations of QSO host galaxies at z >6 have found cold gas and star formation on compact scales of a few kiloparsecs. We present new high sensitivity IRAM PdBI follow-up observations of the [CII] 158micron emission line and FIR continuum in the host galaxy of SDSS J1148+5152, a luminous QSO at redshift 6.4189. We find that a large fraction of the gas traced by [CII] is at high velocities, up to ~1400 km/s relative to the systemic velocity, confirming the presence of a major quasar-driven outflow indicated by previous observations. The outflow has a complex morphology and reaches a maximum projected radius of ~30 kpc. The extreme spatial extent of the outflow allows us, for the first time in an external galaxy, to estimate mass-loss rate, kinetic power and momentum rate of the outflow as a function of the projected distance from the nucleus and the dynamical time-scale. These trends reveal multiple outflow events during the past 100 Myr, although the bulk of the mass, energy and momentum appear to have been released more recently, within the past ~20 Myr. Surprisingly, we discover that also the quiescent gas at systemic velocity is extremely extended. More specifically, we find that, while 30% of the [CII] within v\in(-200, 200) km/s traces a compact component that is not resolved by our observations, 70% of the [CII] emission in this velocity range is extended, with a projected FWHM size of 17.4+-1.4 kpc. We detect FIR continuum emission associated with both the compact and the extended [CII] components, although the extended FIR emission has a FWHM of 11+-3 kpc, thus smaller than the extended [CII] source. Overall, our results indicate that the cold gas traced by [CII] is distributed up to r~30 kpc. A large fraction of extended [CII] is likely associated with star formation on large scales, but the [CII] source extends well beyond the FIR continuum.

We report the serendipitous detection of the CO(17-16) emission line toward the quasar SDSSJ114816.64+525150.3 (J1148) at redshift z = 6.4 obtained with the Plateau de Bure Interferometer. The CO(17-16) line is possibly contaminated by OH+ emission, that may account for ~ 35 - 60% of the total flux observed. Photo-Dissociation and X-ray Dominated Regions (PDRs and XDRs) models show that PDRs alone cannot reproduce the high luminosity of the CO(17-16) line relative to low-J CO transitions and that XDRs are required. By adopting a composite PDR+XDR model we derive molecular cloud and radiation field properties in the nuclear region of J1148. Our results show that highly excited CO lines represent a sensitive and possibly unique tool to infer the presence of X-ray faint or obscured supermassive black hole progenitors in high-z galaxies.

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

We study the properties of the circumgalactic medium (CGM) of high-$z$ galaxies in the metal enrichment simulations presented in Pallottini et al. 2014. At $z=4$, we find that the simulated CGM gas density profiles are self-similar, once scaled with the virial radius of the parent dark matter halo. We also find a simple analytical expression relating the neutral hydrogen equivalent width (${\rm EW}_{\rm HI}$) of CGM absorbers as a function of the line of sight impact parameter ($b$). We test our predictions against mock spectra extracted from the simulations, and show that the model reproduces the ${\rm EW}_{\rm HI}(b)$ profile extracted from the synthetic spectra analysis. When compared with available data, our CGM model nicely predicts the observed ${\rm EW}_{\rm HI}(b)$ in $z\lesssim2$ galaxies, and supports the idea that the CGM profile does not evolve with redshift.

We study cosmic metal enrichment via AMR hydrodynamical simulations in a (10 Mpc/h)$^3$ volume following the Pop III-Pop II transition and for different Pop III IMFs. We have analyzed the joint evolution of metal enrichment on galactic and intergalactic scales at z=6 and z=4. Galaxies account for <9% of the baryonic mass; the remaining gas resides in the diffuse phases: (a) voids, i.e. regions with extremely low density ($\Delta$<1), (b) the true intergalactic medium (IGM, 1<$\Delta$<10) and (c) the circumgalactic medium (CGM, 10<$\Delta<10^{2.5}$), the interface between the IGM and galaxies. By z=6 a galactic mass-metallicity relation is established. At z=4, galaxies with a stellar mass $M_*=10^{8.5}M_\odot$ show log(O/H)+12=8.19, consistent with observations. The total amount of heavy elements rises from $\Omega^{SFH}_Z=1.52\, 10^{-6}$ at z=6 to 8.05 $10^{-6}$ at z=4. Metals in galaxies make up to ~0.89 of such budget at z=6; this fraction increases to ~0.95 at z=4. At z=6 (z=4) the remaining metals are distributed in CGM/IGM/voids with the following mass fractions: 0.06/0.04/0.01 (0.03/0.02/0.01). Analogously to galaxies, at z=4 a density-metallicity ($\Delta$-Z) relation is in place for the diffuse phases: the IGM/voids have a spatially uniform metallicity, Z~$10^{-3.5}$Zsun; in the CGM Z steeply rises with density up to ~$10^{-2}$Zsun. In all diffuse phases a considerable fraction of metals is in a warm/hot (T>$10^{4.5}$K) state. Due to these physical conditions, CIV absorption line experiments can probe only ~2% of the total carbon present in the IGM/CGM; however, metal absorption line spectra are very effective tools to study reionization. Finally, the Pop III star formation history is almost insensitive to the chosen Pop III IMF. Pop III stars are preferentially formed in truly pristine (Z=0) gas pockets, well outside polluted regions created by previous star formation episodes.

We present the search for the [CII] emission line in three $z>6.5$ Lyman-alpha emitters (LAEs) and one J-Dropout galaxy using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Plateau de Bure Interferometer (PdBI). We observed three bright $z\sim6.5-7$ LAEs discovered in the SUBARU deep field (SDF) and the Multiple Imaged lensed $z\sim 11$ galaxy candidate found behind the galaxy cluster MACSJ0647.7+7015. For the LAEs IOK-1 ($z=6.965$), SDF J132415.7+273058 ($z=6.541$) and SDF J132408.3+271543 ($z=6.554$) we find upper limits for the [CII] line luminosity of $<2.05$, $<4.52$ and $<10.56\times10^{8}{\rm L}_{\odot}$ respectively. We find upper limits to the FIR luminosity of the galaxies using a spectral energy distribution template of the local galaxy NGC 6946 and taking into account the effects of the Cosmic Microwave Background on the mm observations. For IOK-1, SDF J132415.7+273058 and SDF J132408.3+271543 we find upper limits for the FIR luminosity of $<2.33$, $3.79$ and $7.72\times10^{11}{\rm L}_{\odot}$ respectively. For the lensed galaxy MACS0647-JD, one of the highest redshift galaxy candidate to date with $z_{\rm ph}=10.7^{+0.6}_{-0.4}$ we put an upper limit in the [CII] emission of $<1.36\times10^{8}\times(\mu/15)^{-1}{\rm L}_{\odot}$ and an upper limit in the FIR luminosity of $<6.1\times10^{10}\times(\mu/15)^{-1}{\rm L}_{\odot}$ (where $\mu$ is the magnification factor). We explore the different conditions relevant for the search for [CII] emission in high redshift galaxies as well as the difficulties for future observations with ALMA and CCAT.

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

We present the results of observations taken with the X-shooter spectrograph devoted to the study of quasars at z ~ 6. This paper focuses on the properties of metals at high redshift traced, in particular, by the C IV doublet absorption systems. Six objects were observed with resolutions ~27 and 34 km/s in the visual, and 37.5 and 53.5 km/s in the near infrared. We detected 102 C IV lines in the range: 4.35 < z < 6.2 of which 27 are above z ~ 5. Thanks to the characteristics of resolution and spectral coverage of X-shooter, we could also detect 25 Si IV doublets associated with the C IV at z>5. The column density distribution function of the C IV line sample is observed to evolve in redshift for z>~ 5.3, with respect to the normalization defined by low redshift (1.5 < z <4) C IV lines. This behaviour is reflected in the redshift evolution of the C IV cosmic mass density, Omega_CIV, of lines with column density in the range 13.4 < log N(C IV) < 15, which is consistent with a drop of a factor of ~2 for z>~ 5.3. Considering only the stronger C IV lines (13.8 < log N(CIV) < 15), Omega_CIV gently rises by a factor of ~10 between z ~ 6.2 and z~ 1.5 with a possible flattening toward z ~ 0. The increase is well fitted by a power law: Omega_CIV = (2+-1)x10^{-8} [(1+z)/4]^{-3.1+-0.1}. An insight into the properties of the C IV absorbers and their evolution with redshift is obtained by comparing the observed column densities of associated C IV, Si IV and C II absorptions with the output of a set of CLOUDY photo-ionization models. As already claimed by cosmological simulations, we find that C IV is a good tracer of the metallicity in the low-density intergalacitc medium gas at z ~ 5-6 while at z ~ 3 it arises in gas with over density delta ~100.

By combining high resolution, radiative transfer cosmological simulations of z~6 galaxies with a sub-grid multi-phase model of their interstellar medium we derive the expected intensity of several far infrared (FIR) emission lines ([C II] 158 micron, [O I] 63 micron, and [N II] 122 micron) for different values of the gas metallicity, Z. For Z = Z_sun the [C II] spectrum is very complex due to the presence of several emitting clumps of individual size < 3 kpc; the peak is displaced from the galaxy center by ~100 km/s. While the [O I] spectrum is also similarly displaced, the [N II] line comes predominantly from the central ionized regions of the galaxy. When integrated over ~500 km/s, the [C II] line flux is 185 mJy km/s; 95% of such flux originates from the cold (T~250 K) H I phase, and only 5% from the warm (T~5000 K) neutral medium. The [O I] and [N II] fluxes are ~6 and ~90 times lower than the [C II] one, respectively. By comparing our results with observations of Himiko, the most extended and luminous Lyman Alpha Emitter (LAE) at z = 6.6, we find that the gas metallicity in this source must be sub-solar. We conclude that the [C II] line from z~6 galaxies is detectable by the ALMA full array in 1.9 < t_ON < 7.7 hr observing time, depending on Z.

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.

We present Plateau de Bure Interferometer observations of far infra-red emission lines in BRI 0952-0115, a lensed quasar at z=4.4 powered by a super-massive black hole (M_BH=2x10^9 M_sun). In this source, the resolved map of the [CII] emission at 158 micron allows us to reveal the presence of a companion galaxy, located at \sim 10 kpc from the quasar, undetected in optical observations. From the CO(5-4) emission line properties we infer a stellar mass M*<2.2x10^10 M_sun, which is significantly smaller than the one found in local galaxies hosting black holes with similar masses (M* \sim 10^12 M_sun). The detection of the [NII] emission at 205 micron suggests that the metallicity in BRI 0952-0115 is consistent with solar, implying that the chemical evolution has progressed very rapidly in this system. We also present PdBI observations of the [CII] emission line in SDSSJ1148+5251, one of the most distant quasar known, at z=6.4. We detect broad wings in the [CII] emission line, indicative of gas which is outflowing from the host galaxy. In particular, the extent of the wings, and the size of the [CII] emitting region associated to them, are indicative of a quasar-driven massive outflow with the highest outflow rate ever found (dM/dt>3500 M_sun/yr).

We present one of the first resolved maps of the [CII] 158 micron line, a powerful tracer of the star forming inter-stellar medium, at high redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-0115, the host galaxy of a lensed quasar at z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two [CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over a region of ~ 1 kpc in size near the quasar nucleus, and we infer a star formation surface density >150 Msun/yr/kpc^2, similar to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~ 12 kpc, and located at ~ 10 kpc from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, in the vicinity of the accreting black hole, and in more extended star forming galaxies.

We present a novel Bayesian method for the joint reconstruction of cosmological matter density fields, peculiar velocities and power-spectra in the quasi-nonlinear regime. We study its applicability to the Ly-alpha forest based on multiple quasar absorption spectra. Our approach to this problem includes a multiscale, nonlinear, two-step scheme since the statistics describing the matter distribution depends on scale, being strongly non-Gaussian on small scales (< 0.1 h^{-1} Mpc) and closely lognormal on scales >~10 h^{-1} Mpc. The first step consists on performing 1D highly resolved matter density reconstructions along the line-of-sight towards z~2-3 quasars based on an arbitrary non-Gaussian univariate model for matter statistics. The second step consists on Gibbs-sampling based on conditional PDFs. The matter density field is sampled in real space with Hamiltonian-sampling using the Poisson/Gamma-lognormal model, while redshift distortions are corrected with linear Lagrangian perturbation theory. The power-spectrum of the lognormal transformed variable which is Gaussian distributed (and thus close to the linear regime) can consistently be sampled with the inverse Gamma distribution function. We test our method through numerical N-body simulations with a computational volume large enough (> 1 h^{-3} Gpc^3) to show that the linear power-spectra are nicely recovered over scales larger than >~20 h^{-1} Mpc, i.e. the relevant range where features imprinted by the baryon-acoustics oscillations (BAOs) appear.

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

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

We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to their selection and to the characterization of their distribution of UV slopes beta, using deep IR images from HST. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of beta. The same fit is used to determine their rest-frame UV magnitudes M1600 through a simple interpolation. We estimate the effect of observational uncertainties with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density function of UV slopes PDF(beta) as a function of rest-frame magnitude. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our newly defined criteria improve the selection of z~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 objects. We find that z~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: all L<L* galaxies have an average UV slope <beta>\simeq-2.1, while brighter objects only are slightly redder (<beta> -1.9). We find an intrinsic dispersion ~ 0.3 for PDF(beta) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower than previously found: log(SFRD)\simeq-1.09 M_sun/yr/Mpc^3. Finally, we discuss how the UV slope of z~4 galaxies changes as a function of the dust-corrected UV magnitude (i.e. SFR) showing that most galaxies with a high SFR (> 80 M_sun/yr) are highly extincted objects. [Abridged]

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.

Claim of dust extinction for this GRB has been debated in the past. We suggest that the discrepant results occur primarily because most of previous studies have not simultaneously investigated the X-ray to near-IR spectral energy distribution of this GRB. The difficulty with this burst is that the X-ray afterglow is dominated by strong flares at early times and is poorly monitored at late times. In addition, the Z band photometry, which is the most sensitive to dust extinction, has been found to be affected by strong systematics. In this paper we carefully re-analyze the Swift/XRT afterglow observations of this GRB, using extensive past studies of X-ray flare properties when computing the X-ray afterglow flux level and exploiting the recent reanalysis of the optical (UV rest frame) data of the same GRB. We extract the X-ray to optical/near-IR afterglow SED for the three epochs where the best spectral coverage is available: 0.47, 1.25, and 3.4 days after the trigger. A spectral power-law model has been fitted to the extracted SEDs. We discuss that no spectral breaks or chromatic temporal breaks are expected in the epochs of interest. To fit any UV rest-frame dust absorption, we tested the Small Magellanic Cloud (SMC) extinction curve, the mean extinction curve (MEC) found for a sample of QSO at $z>4$ and its corresponding attenuation curve, as well as a starburst attenuation curve, and the extinction curve consistent with a supernova dust origin (SN-type). The SMC extinction curve and the SN-type one provide good fit to the data at all epochs, with an average amount of dust absorption at $\lambda_{rest} = 3000 \AA$ of $A_{3000} = 0.25\pm 0.07$ mag. These results indicate that the primeval galaxy at $z = 6.3$ hosting this GRB has already enriched its ISM with dust.

We present a novel, fast method to recover the density field through the statistics of the transmitted flux in high redshift quasar absorption spectra. The proposed technique requires the computation of the probability distribution function of the transmitted flux (P_F) in the Ly-alpha forest region and, as a sole assumption, the knowledge of the probability distribution function of the matter density field (P_Delta). We show that the probability density conservation of the flux and matter density unveils a flux-density (F-Delta) relation which can be used to invert the Ly-alpha forest without any assumption on the physical properties of the intergalactic medium. We test our inversion method at z=3 through the following steps: [i] simulation of a sample of synthetic spectra for which P_Delta is known; [ii] computation of P_F; [iii] inversion of the Ly-alpha forest through the F-Delta relation. Our technique, when applied to only 10 observed spectra characterized by a signal-to noise ratio S/N >= 100 provides an exquisite (relative error epsilon_Delta <~ 12 % in >~ 50 % of the pixels) reconstruction of the density field in >~ 90 % of the line of sight. We finally discuss strengths and limitations of the method.

X-shooter, with its characteristics of resolution, spectral coverage and efficiency, provides a unique opportunity to obtain spectra of the highest-redshift quasars (z ~ 6) that will allow us to carry out successful investigations on key cosmological issues, from the details of the re-ionization process, to the evolution of the first galaxies and AGNs. In this paper, we present the spectra of three z ~ 6 quasars: one obtained during the commissioning of X-shooter and two in the context of our ongoing GTO programme. Combining this sample with data in the literature, we update the value of the C IV cosmic mass density in the range 4.5 < z < 5, confirming the constant trend with redshift between 2.5 and 5.

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.

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.

We investigate the cosmic reionization history by comparing semi-analytical models of the Lya forest with observations of high-z quasars and gamma ray bursts absorption spectra. In order to constrain the reionization epoch z_rei, we consider two physically motivated scenarios in which reionization ends either early (ERM, z_rei>= 7) or late (LRM, z_rei~6). We analyze the transmitted flux in a sample of 17 quasars spectra at 5.7<z<6.4 and in the spectrum of the gamma ray burst 050904 at z=6.3, studying the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that current observational data do not require any sudden change in the ionization state of the IGM at z~6, favouring indeed a highly ionized Universe at these epochs, as predicted by the ERM. Moreover, we test the predictions of this model through Lya emitters observations, finding that the ERM provide a good fit to the evolution of the luminosity function of Lya emitting galaxies in the redshift range z=5.7-6.5. The overall result points towards an extended reionization process which starts at z>=11 and completes at z_rei>=7, in agreement with the recent WMAP5 data.

During the last few decades scientists have been able to test the bases of the physics paradigms, where the quantum mechanics has to match the cosmological scales. Between the extremes of this scenario, biological phenomena and their complexity take place, challenging the laws we observe in the atomic and sub-atomic world. In order to explore the details of this world, new huge experimental facilities are under construction. These projects involve people coming from several countries and give physicists the opportunity to work together with chemists, biologists and other scientists. The Roman Young Researchers Meeting is a conference, organised by Ph. D. students and young postdocs connected to the Roman area. It is aimed primarily at graduate students and post-docs, working in physics. The 1st conference has been held on the 21st of July 2009 at the University of Roma Tor Vergata. It was organised in three sessions, devoted to Astrophysics and Cosmology, Soft and Condensed Matter Physics and Theoretical and Particle Physics. In this proceeding we collect the contributions which have been presented and discussed during the meeting, according to the specific topics treated.

We propose the Apparent Shrinking Criterion (ASC) to interpret the spatial extent, R_w, of transmitted flux windows in the absorption spectra of high-z quasars. The ASC can discriminate between the two regimes in which R_w corresponds either to the physical size, R_HII, of the quasar HII region, or to the distance, R^{max}_w, at which the transmitted flux drops to =0.1 and a Gunn-Peterson (GP) trough appears. In the first case (HR regime), one can determine the IGM mean HI fraction, x_HI; in the second (PR regime), the value of R_w allows to measure the local photoionization rate and the local enhancement of the photoionization rate, Gamma_G, due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+SPH numerical simulations, and applied to 15 high-z (z>5.8) quasars sample from Fan et al. (2006). All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of x_HI. Four sample quasars show evidence for a local enhancement (up to 50%) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.

We provide measurements of the neutral hydrogen fraction xHI at z~6, by comparing semi-analytical models of the Lyalpha forest with observations of high-z quasars and Gamma Ray Bursts absorption spectra. We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74<zem<6.42 studying separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with our models, we conclude that xHI evolves smoothly from 10^{-4.4} at z=5.3 to 10^{-4.2} at z=5.6, with a robust upper limit xHI<0.36 at z=6.3. We show the results of the first-ever detected transverse proximity effect in the HI Lyalpha forest, produced by the HII region of the faint quasar RD J1148+5253 at z=5.70 intervening along the LOS of SDSS J1148+5251 at z=6.42. Moreover, we propose a novel method to study cosmic reionization using absorption line spectra of high-redshift GRBs afterglows. We show that the time evolution and the statistics of gaps in the observed spectra represent exquisite tools to discriminate among different reionization models. By applying our methods to GRB 050904 detected at z=6.29, we show that the observation of this burst provides strong indications of a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI=6.4\pm 0.3\times 10^{-5} along that line of sight.

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) in three different scenarios: i) GRBs follow the cosmic star formation and their LF is constant in time; ii) GRBs follow the cosmic star formation but the LF varies with redshift; iii) GRBs form preferentially in low-metallicity environments. We then test model predictions against the Swift 3-year data, showing that scenario i) is robustly ruled out. Moreover, we show that the number of bright GRBs detected by Swift suggests that GRBs should have experienced some sort of luminosity evolution with redshift, being more luminous in the past. Finally we propose to use the observations of the afterglow spectrum of GRBs at z>5.5 to constrain the reionization history and we applied our method to the case of GRB 050904.

We use a semi-analytic model of Lyman alpha emitters (LAEs) to constrain the reionization history. By considering two physically motivated scenarios in which reionization ends either early (ERM, z_i ~ 7) or late (LRM, z_i ~ 6), we fix the global value of the IGM neutral fraction (e.g. chi_{HI}=3 times 10^{-4}, 0.15 at z=6.56 for the ERM and LRM, respectively) leaving only the star formation efficiency and the effective escape fraction of Lya photons as free parameters. The ERM fits the observed LAE luminosity function (LF) at z=5.7 and 6.56 requiring no redshift evolution or mass dependence of the star formation efficiency, and LAE star formation rates (SFR) of 3-103 solar masses/year, contributing approximately 8% of the cosmic SFR density at z=5.7. The LRM requires a physically uncomfortable drop of approximately 4.5 times in the SFR of the emitters from z=6.5 to 5.7. Thus, the data seem to imply that the Universe was already highly ionized at z=6.56. The mass-dependent Lya transmissivity is between 0.36-0.51 (ERM) and less than 0.26 (LRM) at z=6.56. The LF data at z=4.5 imply an extra Lya line damping factor of approximately 0.25 possibly due to dust; the presence of a (clumpy) dust component with E(B-V) ~ 0.28 is also required to reproduce the observed large Lya equivalent widths at the same redshift. Additional useful information can be extracted from the line profile (weighted) skewness, found to be S_W=10-17 Angstrom for the two reionization models, which shows an interesting L_alpha-chi_{HI} anti-correlation, holding under the model assumptions. The shortcomings of the model and strategies to overcome them are discussed.

We propose to study cosmic reionization using absorption line spectra of high-redshift Gamma Ray Burst (GRB) afterglows. We show that the statistics of the dark portions (gaps) in GRB absorption spectra represent exquisite tools to discriminate among different reionization models. We then compute the probability to find the largest gap in a given width range [Wmax, Wmax + dW] at a flux threshold Fth for burst afterglows at redshifts 6.3 < z < 6.7. We show that different reionization scenarios populate the (Wmax, Fth) plane in a very different way, allowing to distinguish among different reionization histories. We provide here useful plots that allow a very simple and direct comparison between observations and model results. Finally, we apply our methods to GRB 050904 detected at z = 6.29. We show that the observation of this burst strongly favors reionization models which predict a highly ionized intergalactic medium at z~6, with an estimated mean neutral hydrogen fraction xHI = 6.4 \pm 0.3 \times 10^-5 along the line of sight towards GRB 050904.

We analyze the transmitted flux in a sample of 17 QSOs spectra at 5.74<zem<6.42 to obtain tighter constraints on the volume-averaged neutral hydrogen fraction, xHI, at z~6. We study separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with Lyalpha forest simulations, we conclude that xHI evolves smoothly from 10^{-4.4} at z=5.3 to 10^{-4.2} at z=5.6, with a robust upper limit xHI<0.36 at z=6.3. The frequency and physical sizes of the peaks imply an origin in cosmic underdense regions and/or in HII regions around faint quasars or galaxies. In one case (the intervening HII region of the faint quasar RD J1148+5253 at z=5.70 along the LOS of SDSS J1148+5251 at z=6.42) the increase of the peak spectral density is explained by the first-ever detected transverse proximity effect in the HI Lyalpha forest; this indicates that at least some peaks result from a locally enhanced radiation field. We then obtain a strong lower limit on the foreground QSO lifetime of tQ>11 Myr. The observed widths of the peaks are found to be systematically larger than the simulated ones. Reasons for such discrepancy might reside either in the photoionization equilibrium assumption or in radiative transfer effects.

By means of carefully calibrated semi-analytical reionization models, we estimate the minimum mass of star-forming haloes required to match the current data. Models which do not include haloes of total mass M < 10^9 M_sun fail at reproducing the Gunn-Peterson and electron scattering optical depths simultaneously, as they contribute too few (many) photons at high (low, z \approx 6) redshift. Marginally acceptable solutions require haloes with M \approx 5 \times 10^7 M_sun at z \approx 10, corresponding to virial temperatures (\sim 10^4K) for which cooling can be ensured by atomic transitions. However, a much better match to the data is obtained if minihaloes (M \sim 10^6 M_sun) are included in the analysis. We have critically examined the assumptions made in our model and conclude that reionization in the large-galaxies-only scenario can remain viable only if metal-free stars and/or some other exotic sources at z > 6 are included.

We investigate the possibility of constraining the ionization state of the intergalactic medium (IGM) close to the end of reionization (z ~ 6) by measuring the size of the HII regions in high-z quasars spectra. We perform a combination of multiphase smoothed particle hydrodynamics (SPH) and 3D radiative transfer (RT) simulations to reliably predict the properties of typical high-z quasar HII regions, embedded in a partly neutral IGM (x_HI=0.1). In this work we assume a fixed configuration for the quasar lifetime and luminosity, i.e. t_Q=10^7 yr and N_gamma=5.2 10^56 s-1. From the analysis of mock spectra along lines of sight through the simulated QSO environment, we find that the HII region size derived from quasar spectra is on average 30 per cent smaller than the physical one. Additional maximum likelihood analysis shows that this offset induces an overestimate of the neutral hydrogen fraction, x_HI, by a factor of 3. By applying the same statistical method to a sample of observed QSOs our study favors a mostly ionized (x_HI < 0.06) universe at z=6.1.

We use a semi-analytical approach to simulate absorption spectra of QSOs at high redshifts with the aim of constraining the cosmic reionization history. We consider two physically motivated and detailed reionization histories: (i) an Early Reionization Model (ERM) in which the intergalactic medium is reionized by PopIII stars at $z\approx 14$, and (ii) a more standard Late Reionization Model (LRM) in which overlapping, induced by QSOs and normal galaxies, occurs at $z\approx 6$. From the analysis of current Ly$\alpha$ forest data at $z < 6$, we conclude that it is impossible to disentangle the two scenarios, which fit equally well the observed Gunn-Peterson optical depth, flux probability distribution function and dark gap width distribution. At $z>6$, however, clear differences start to emerge which are best quantified by the dark gap and peak width distributions. We find that 35 (zero) per cent of the lines of sight within $5.7< z <6.3$ show dark gaps widths $>50$ Angstrom in the rest frame of the QSO if reionization is not (is) complete at $z \gtrsim 6$. Similarly, the ERM predicts peaks of width $\sim 1$ Angstrom in 40 per cent of the lines of sight in the redshift range $6.0-6.6$; in the same range, LRM predicts no peaks of width $>0.8$ Angstrom. We conclude that the dark gap and peak width statistics represent superb probes of cosmic reionization if about ten QSOs can be found at $z > 6$. We finally discuss strengths and limitations of our method.