We present results from the first cosmological simulations which study the onset of primordial, metal-free (population III), cosmic star formation and the transition to the present-day, metal-rich star formation (population II-I), including molecular (H$_2$, HD, etc.) evolution, tracing the injection of metals by supernov{\ae} into the surrounding intergalactic medium and following the change in the initial stellar mass function (IMF) according to the metallicity of the corresponding stellar population. Our investigation addresses the role of a wide variety of parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and the associated transition in the star formation mode. All simulations present common trends. Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, inducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a result of the rapid pollution within high-density regions due to the first SN/pair-instability SN, local metallicity is quickly boosted above the critical metallicity for the transition. The population III regime lasts for a very short period during the first stages of star formation ($\sim 10^7\,\rm yr$), and its average contribution to the total star formation rate density drops rapidly below $\sim 10^{-3}-10^{-2}$

Cooling is the main process leading to the condensation of gas in the dark matter potential wells and consequently to star and structure formation. In a metal-free environment, the main available coolants are H, He, H2 and HD; once the gas is enriched with metals, these also become important in defining the cooling properties of the gas. We discuss the implementation in GADGET-2 of molecular and metal cooling at temperatures lower than 104 K, following the time-dependent properties of the gas and pollution from stellar evolution. We have checked the validity of our scheme by comparing the results of some test runs with previous calculations of cosmic abundance evolution and structure formation, finding excellent agreement. We have also investigated the relevance of molecule and metal cooling in some specific cases, finding that inclusion of HD cooling results in a higher clumping factor of the gas at high redshifts, while metal cooling at low temperatures can have a significant impact on the formation and evolution of cold objects

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 M star > 10 9 M ⊙ allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population.

We present results of high‐resolution hydrodynamical simulations of the formation and evolution of dwarf galaxies. Our simulations start from cosmological initial conditions at high redshift. They include metal‐dependent cooling, star formation, feedback from Type II and Ia supernovae and UV background radiation, with physical recipes identical to those applied in a previous study of Milky Way type galaxies. We find that a combination of feedback and the cosmic UV background results in the formation of galaxies with properties similar to the Local Group dwarf spheroidals, and that their effect is strongly moderated by the depth of the gravitational potential. Taking this into account, our models naturally reproduce the observed luminosities and metallicities. The final objects have halo masses between 2.3 × 108 and 1.1 × 109 M⊙, mean velocity dispersions between 6.5 and 9.7 km s−1, stellar masses ranging from 5 × 105 to 1.2 × 107 M⊙, median metallicities between [Fe/H]=−1.8 and −1.1 and half‐light radii of the order of 200–300 pc, all comparable with Local Group dwarf spheroidals. Our simulations also indicate that the dwarf spheroidal galaxies observed today lie near a halo mass threshold around 109 M⊙, in agreement with stellar kinematic data, where supernova feedback not only suffices to completely expel the interstellar medium and leave the residual gas free, but where the combination of feedback, UV radiation and self‐shielding establishes a dichotomy of age distributions similar to that observed in the Milky Way and M31 satellites.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) followup with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift ∼10, signatures of Pop III stars, sources and physics of reionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).