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

The origin of carbon-enhanced metal-poor (CEMP) stars and their possible connections with the chemical elements produced by the first stellar generations is still highly debated. We briefly review observations of CEMP stars in different environments (Galactic stellar halo, ultra-faint and classical dwarf galaxies) and interpret their properties using cosmological chemical-evolution models for the formation of the Local Group. We discuss the implications of current observations for the properties of the first stars, clarify why the fraction of carbon-enhanced to carbon-normal stars varies in dwarf galaxies with different luminosity, and discuss the origin of the first CEMP(-no) star found in the Sculptor dwarf galaxy.

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.

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