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

We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available.

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

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

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

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

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

Here we investigate how LyC-opaque systems present in the intergalactic medium at $z\approx3$ can distort the spectral shape of a uniform UV background (UVB) through radiative transfer (RT) effects. With this aim in mind, we perform a multi-frequency RT simulation through a cosmic volume of $10h^{-1}$~cMpc scale polluted by metals, and self-consistently derive the ions of all the species. The UVB spatial fluctuations are traced by the ratio of He$\, \rm \scriptstyle II\ $ and H$\, \rm \scriptstyle I\ $ column density, $\eta$, and the ratio of C$\,{\rm {\scriptstyle IV\ }}$ and Si$\,{\rm {\scriptstyle IV\ }}$ optical depths, $\zeta$. We find that: (i) $\eta$ spatially fluctuates through over-dense systems ($\Delta$) with statistically significant deviations $\delta\eta >25$\% in 18\% of the volume ; (ii) same fluctuations in $\zeta$ are also present in $34$\% of the enriched domain (only 8\% of the total volume) and derive from a combination of RT induced effects and in-homogeneous metal enrichment, both effective in systems with $\Delta > 1.5$.

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

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

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission.

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

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

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

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

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

We introduce CRASH-AMR, a new version of the cosmological Radiative Transfer (RT) code CRASH, enabled to use refined grids. This new feature allows us to attain higher resolution in our RT simulations and thus to describe more accurately ionisation and temperature patterns in high density regions. We have tested CRASH-AMR by simulating the evolution of an ionised region produced by a single source embedded in gas at constant density, as well as by a more realistic configuration of multiple sources in an inhomogeneous density field. While we find an excellent agreement with the previous version of CRASH when the AMR feature is disabled, showing that no numerical artifact has been introduced in CRASH-AMR, when additional refinement levels are used the code can simulate more accurately the physics of ionised gas in high density regions. This result has been attained at no computational loss, as RT simulations on AMR grids with maximum resolution equivalent to that of a uniform cartesian grid can be run with a gain of up to 60% in computational time.

We present a detailed analysis of the ionization and thermal structure of the intergalactic medium (IGM) around a high-redshift QSO using a large suite of cosmological, multi-frequency radiative transfer (RT) simulations, exploring the contribution from galaxies as well as the QSO, and the effect of X-rays and secondary ionization. We show that in high-z QSO environments both the central QSO and the surrounding galaxies concertedly control the reionization morphology of hydrogen and helium and have a non-linear impact on the thermal structure of the IGM. A QSO imprints a distinctive morphology on H II regions if its total ionizing photon budget exceeds that of the surrounding galaxies since the onset of hydrogen reionization; otherwise, the morphology shows little difference from that of H II regions produced only by galaxies. In addition, the spectral shape of the collective radiation field from galaxies and QSOs controls the thickness of the I-fronts. While a UV-obscured QSO can broaden the I-front, the contribution from other UV sources, either galaxies or unobscured QSO, is sufficient to maintain a sharp I-front. X-rays photons from the QSO are responsible for a prominent extended tail of partial ionization ahead of the I-front. QSOs leave a unique imprint on the morphology of He II / He III regions. We suggest that, while the physical state of the IGM is modified by QSOs, the most direct test to understand the role of galaxies and QSOs during reionization is to perform galaxy surveys in a region of sky imaged by 21 cm tomography.

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

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.

The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

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

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

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

Here we introduce CRASH3, the latest release of the 3D radiative transfer code CRASH. In its current implementation CRASH3 integrates into the reference algorithm the code Cloudy to evaluate the ionisation states of metals, self-consistently with the radiative transfer through H and He. The feedback of the heavy elements on the calculation of the gas temperature is also taken into account, making of CRASH3 the first 3D code for cosmological applications which treats self-consistently the radiative transfer through an inhomogeneous distribution of metal enriched gas with an arbitrary number of point sources and/or a background radiation. The code has been tested in idealized configurations, as well as in a more realistic case of multiple sources embedded in a polluted cosmic web. Through these validation tests the new method has been proven to be numerically stable and convergent. We have studied the dependence of the results on a number of physical quantities such as the source characteristics (spectral range and shape, intensity), the metal composition, the gas number density and metallicity.

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

We investigate the effect of helium on hydrogen reionisation using a hydrodynamical simulation combined with the cosmological radiative transfer code CRASH. The simulations are run in a 35.12/h comoving Mpc box using a variety of assumptions for the amplitude and power-law extreme-UV (EUV) spectral index, alpha, of the ionising emissivity. We use an empirically motivated prescription for ionising sources which ensures all of the models are consistent with constraints on the Thomson scattering optical depth and the hydrogen photo-ionisation rate at z=6. The inclusion of helium slightly delays reionisation due to the small number of ionising photons which reionise neutral helium instead of hydrogen. However, helium has a significant impact on the thermal state of the IGM. Models with alpha=3 produce IGM temperatures at the mean density at z=6 which are about 20 % higher compared to models without helium photo-heating. Harder EUV indices produce even larger IGM temperature boosts. A comparison to recent observational estimates of the IGM temperature at z=5 - 6 suggests that hydrogen reionisation was primarily driven by pop-II stellar sources with a soft EUV index, alpha<3. We also find that faint, as yet undetected galaxies, characterised by a luminosity function with a steepening faint-end slope and an increasing Lyman continuum escape fraction (fesc=0.5), are required to reproduce the ionising emissivity used in our simulations at z>6. Finally, we note there is some tension between recent observational constraints which indicate the IGM is > 10% neutral by volume z=7, and estimates of the ionising emissivity at z=6 which indicate only between 1 and 3 ionising photons are emitted per hydrogen atom over a Hubble time. This tension may be alleviated by either a lower neutral fraction at z=7 or an IGM which still remains a few % neutral by volume at z=6.