We study the formation of the first massive galaxies and their observational properties for comparisons with JWST or ALMA data by performing cosmological radiative-hydrodynamics simulations. We find that galaxies in overdense regions have high star formation rates larger than ∼10 M⊙ yr-1 at z=10 and their stellar masses reach 109 M⊙. The star formation rates of our model galaxies at z≲12 nicely match recent JWST data. In addition, we show that the morphology of galaxies drastically changes with time via major mergers and stellar feedback.

Local Group (LG) very metal-poor massive stars are the best proxy for the First Stars of the Universe, fundamental to model the early evolution of the first galaxies, and key to unravelling new evolutionary pathways restricted to low metallicities, such as chemically homogeneous evolution. Yet, due to the great leap in distance required to reach metal-poor galaxies of the LG and vicinity, no comprehensive spectroscopic studies have been carried out at sub-$$$$ metallicities so far.

We focus on the massive star population of the 1/10Ȯ galaxy Sextans A. After five observing campaigns at the 10.4-m Gran Telescopio Canarias (GTC), we have assembled a spectroscopic catalogue of more than 150 OB stars. This catalogue will be fundamental to test stellar evolution at very low metallicity, to detect the first binary systems at 1/10Ȯ and to unveil the most recent star formation history of this galaxy.

The first stars in the Universe have inherited their composition from primordial nucleosynthesis, so they have no metal. These stars, which are also named population III (pop III) stars, began the process of reionization in the Universe and contributed to the metal enrichment with heavy elements. Previous studies showed that they should have been rotating fast due to small or no angular momentum loss, reaching easily the critical velocity since they are massive and have very low stellar winds, thus their mass loss is very low or zero. Our aim is to study how the production of primary nitrogen is affected due to high rotation in the pop III stars. So, we compared grids of pop III stars with zero, average, and high rotation. All these models have been computed using Geneva code (GENEC) in the mass range of 9M⊙ ≤ Mini ≤ 120M⊙. Due to the rotational mixing, the carbon produced in the He-burning core is diffused towards the H-burning shell, triggering the CNO cycle and producing primary nitrogen. In some models the transition of the shell from a pp-chain H-burning to a CNO H-burning induces a strong energy release and a complete change of the stellar structure and the nucleosynthesis. The production of nitrogen is boosted for the high rotation models.

During the recombination of the universe, supersonic relative motion between baryons and dark matter (DM) generally existed. In the presence of such streaming motions, gas clumps can collapse outside of virial radii of their closest dark matter halos. Such baryon dominant objects are thought to be self-gravitating and are called supersonically induced gas objects; SIGOs. We perform three-dimensional hydrodynamical simulations by including H2 chemical reactions and stream velocity and follow SIGO’s formation from z = 200 to z = 25. SIGOs can be formed under the influence of stream velocity, and cooling is effective in contracting gas clouds. We follow its further evolution with higher resolution. We find that there are SIGOs which become Jeans unstable outside of the virial radius of the closest DM halos. Those SIGOs are gravitationally unstable and trigger star formation.

Getting a better understanding of the evolution and nucleosynthetic yields of the most metal-poor stars (Z ≲ 10−5) is critical because they are part of the big picture of the history of the primitive universe. Yet many of the remaining unknowns of stellar evolution lie in the birth, life, and death of these objects. We review stellar evolution of intermediate-mass Z ≤ 10−5 models existing in the literature, with a particular focus on the problem of their final fates. We emphasise the importance of the mixing episodes between the stellar envelope and the nuclearly processed core, which occur after stars exhaust their central He (second dredge-up and dredge-out episodes). The depth and efficiency of these episodes are critical to determine the mass limits for the formation of electron-capture SNe. Our knowledge of these phenomena is not complete because they are strongly affected by the choice of input physics. These uncertainties affect stars in all mass and metallicity ranges. However, difficulties in calibration pose additional challenges in the case of the most metal-poor stars. We also consider the alternative SN I1/2 channel to form SNe out of the most metal-poor intermediate-mass objects. In this case, it is critical to understand the thermally pulsing Asymptotic Giant Branch evolution until the late stages. Efficient second dredge-up and, later, third dredge-up episodes could be able to pollute stellar envelopes enough for the stars to undergo thermal pulses in a way very similar to that of higher initial Z objects. Inefficient second and/or third dredge-up may leave an almost pristine envelope, unable to sustain strong stellar winds. This may allow the H-exhausted core to grow to the Chandrasekhar mass before the envelope is completely lost, and thus let the star explode as an SN I1/2. After reviewing the information available on these two possible channels for the formation of SNe, we discuss existing nucleosynthetic yields of stars of metallicity Z ≤ 10−5 and present an example of nucleosynthetic calculations for a thermally pulsing Super-Asymptotic Giant Branch star of Z = 10−5. We compare theoretical predictions with observations of the lowest [Fe/H] objects detected. The review closes by discussing current open questions as well as possible fruitful avenues for future research.

The double distance and metallicity frontier marked by the SMC has been finally broken with the aid of powerful multi-object spectrographs installed at 8-10m class telescopes. VLT, GTC and Keck have enabled studies of massive stars in dwarf irregular galaxies of the Local Group with poorer metal-content than the SMC. The community is working to test the predictions of evolutionary models in the low-metallicity regime, set the new standard for the metal-poor high-redshift Universe, and test the extrapolation of the physics of massive stars to environments of decreasing metallicity. In this paper, we review current knowledge on this topic.

The properties of the first generation of stars and their supernova (SN) explosions remain unknown due to the lack of their actual observations. Pop III stars may have been very massive and predicted to be exploded as pair-instability SNe, but the observed metal-poor stars show the abundance patterns which are more consistent with yields of core-collapse SNe. We study the multicolor light curves for a metal-free core-collapse SN models (massive stars of 25-100 solar mass range) to determine the indicators for the detection and identification of first generation SNe. We use mixing-fallback supernova explosion models which explain the observed abundance patterns of metal poor stars. Numerical calculations of the multicolor light curves are performed using the multigroup radiation hydrodynamic code STELLA. The calculated light curves of metal-free SNe are compared with our calculations of non-zero metallicity models and observed SNe.

Gamma-ray bursts (GRBs) are unique probes of the first generation (Pop III) stars. We show that a relativistic gamma-ray burst (GRB) jet can potentially pierce the envelope of a very massive Pop III star even if the Pop III star has a supergiant hydrogen envelope without mass loss, thanks to the long-lived powerful accretion of the envelope itself. While the Pop III GRB is estimated to be energetic (Eγ,iso ~ 1055 erg), the supergiant envelope hides the initial bright phase in the cocoon component, leading to a GRB with a long duration ~1000 (1 + z) s and an ordinary isotropic luminosity ~ 1052 erg s−1 (~ 10−9 erg cm−2 s−1 at redshift z ~ 20), although these quantities are found to be sensitive to the core and envelope mass. We also show that Pop III.2 GRBs (which are primordial but affected by radiation from other stars) occur >100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future X-ray and radio missions. The radio transient surveys are already constraining the Pop III GRB rate and promising in the future.

We review the properties of massive Population III and very metal-poor stars, including briefly their formation, IMF, their main sequence evolution, possible mass loss mechanisms, atmosphere modeling etc. For detailed predictions concerning the properties of these stars we refer to Schaerer (2002) and Schaerer (2002) and references therein. Extending these calculations, Schaerer (2007) present new calculations concerning the ionizing power, Ly-α strength and related properties for different metallicities as well as for a range of power-law and log-normal IMFs. For illustrations from these studies see the Figures below. New detailed calibrations for solar metallicity O-type stars have recently been presented by Martins et al. (2005).