The first fireworks: A roadmap to Population III stars during the epoch of reionization through pair-instability supernovae

With the launch of JWST and other scheduled missions aimed at probing the distant universe, we are entering a new promising era for high-z astronomy. One of our main goals is the detection of the first population of stars (Population III or Pop III stars), and models suggest that Pop III star formation is allowed well into the Epoch of Reionization (EoR), rendering this an attainable achievement. In this paper, we focus on our chance of detecting massive Pop IIIs at the moment of their death as Pair-Instability Supernovae (PISNe). We estimate the probability of discovering PISNe during the EoR in galaxies with different stellar masses (7.5 ≤ Log(M/M) ≤ 10.5) from six dustyGadget simulations of 50 h-1 cMpc per side. We further assess the expected number of PISNe in surveys with JWST/NIRCam and Roman/WFI. On average, less than one PISN is expected in all examined JWST fields at z ≃ 8 with Δz = 1, and O(1) PISN may be found in a ~1 deg2 Roman field in the best-case scenario, although different assumptions on the Pop III IMF and/or Pop III star formation efficiency can decrease this number substantially. Including the contribution from unresolved low-mass haloes holds the potential for increased discoveries. JWST/NIRCam and Roman/WFI allow the detection of massive-progenitor (~250 M) PISNe throughout all the optimal F200W-F356W, F277W-F444W, and F158-F213 colours. PISNe are also pre-dominantly located at the outskirts of their hosting haloes, facilitating the disentangling of underlying stellar emission thanks to the spatial-resolution capabilities of the instruments.

    FULL ARTICLE

    A needle in a haystack? Catching Population III stars in the epoch of reionization: I. Population III star-forming environments

    Despite extensive search efforts, direct observations of the first [Population III (Pop III)] stars have not yet succeeded. Theoretical studies have suggested that late Pop III star formation is still possible in pristine clouds of high-mass galaxies, co-existing with Pop II stars, down to the epoch of reionization. Here, we re-assess this finding by exploring Pop III star formation in six 50 h-1 cMpc simulations performed with the hydrodynamical code dustyGadget. We find that Pop III star formation (~10-3.4-10-3.2 M yr-1 cMpc-3) is still occurring down to z ~ 6-8, i.e. well within the reach of deep JWST surveys. At these epochs, ≲ 10% of the rare massive galaxies with M ≳ 3 × 109 M are found to host Pop III stars, although with a Pop III/Pop II mass fraction ≲ 0.1%. Regardless of their mass, Pop III-hosting galaxies are mainly found on the main sequence, at high star-formation rates, probably induced by accretion of pristine gas. This scenario is also supported by their increasing star-formation histories and their preferential location in high-density regions of the cosmic web. Pop III stars are found both in the outskirts of metal-enriched regions and in isolated, pristine clouds. In the latter case, their signal may be less contaminated by Pop IIs, although its detectability will strongly depend on the specific line of sight to the source, due to the complex morphology of the host galaxy and its highly inhomogeneous dust distribution.

    FULL ARTICLE

    Peering into cosmic reionization: the Lyα visibility evolution from galaxies at z=4.5-8.5 with JWST

    The resonant scattering interaction between Ly α photons and neutral hydrogen implies that a partially neutral IGM can significantly impact the detectability of Ly α emission in galaxies. The redshift evolution of the Ly α equivalent width distribution of galaxies thus offers a key probe of the degree of ionization during the Epoch of Reionization (EoR). Previous in-depth investigations at z ≥ 7 were limited by ground-based instrument capabilities. We present an extensive study of Ly α emission from galaxies at 4 < z < 8.5, observed from the CEERS and JADES surveys in the JWST NIRSpec/PRISM configuration. The sample consists of 235 galaxies, among which we identify 65 as Ly α emitters. We first measure Ly α escape fractions from Balmer lines, and explore the correlations with the inferred galaxies’ physical properties, which are similar to those found at lower redshift. We also investigate the possible connection between the escape of Ly α photons and the inferred escape fractions of LyC photons obtained from indirect indicators. We then analyze the redshift evolution of the Ly α emitter fraction, finding lower average values at z = 5 and 6 compared to ground-based observations. At z = 7 we find a very large difference in Ly α visibility between the EGS and GOODS-South fields, possibly due to the presence of early reionized regions in the EGS. Such large variance is also expected in the Cosmic Dawn II radiation-hydrodynamical simulation. Our findings suggest a scenario in which the ending phase of the EoR is characterized by ∼ 1 pMpc ionized bubbles around a high fraction of moderately bright galaxies. Finally, we characterize such two ionized regions found in the EGS at z = 7.18 and z = 7.49 by estimating the radius of the ionized bubble that each of the spectroscopically-confirmed members could have created.

    FULL ARTICLE

    New insight on the nature of cosmic reionizers from the CEERS survey

    The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction ( fesc). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their fesc can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at 6≤z≤9 from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average 6<z<9 star-forming galaxies are compact in the rest-frame UV ( re∼0.4 kpc), are blue sources (UV- β slope ∼-2.17), and have a predicted fesc of about 0.13. A comparison of our results to models and predictions as well as an estimation of the ionizing budget suggests that low-mass galaxies with UV magnitudes fainter than M1500=−18 that we currently do not characterize with JWST observations probably played a key role in the process of reionization.

    FULL ARTICLE

    Identifying Lyα emitter candidates with Random Forest: Learning from galaxies in the CANDELS survey

    The physical processes that make a galaxy a Lyman alpha emitter have been extensively studied over the past 25 yr. However, the correlations between physical and morphological properties of galaxies and the strength of the Lyα emission line are still highly debated. Here, we investigate the correlations between the rest-frame Lyα equivalent width and stellar mass, star formation rate, dust reddening, metallicity, age, half-light semi-major axis, Sérsic index, and projected axis ratio in a sample of 1578 galaxies in the redshift range of 2 ≤ z ≤ 7.9 from the GOODS-S, UDS, and COSMOS fields. From the large sample of Lyα emitters (LAEs) in the dataset, we find that LAEs are typically common main sequence (MS) star-forming galaxies that show a stellar mass ≤109 M, star formation rate ≤ 100.5 M yr−1, E(B − V)≤0.2, and half-light semi-major axis ≤1 kpc. Building on these findings, we have developed a new method based on a random forest (RF) machine learning (ML) classifier to select galaxies with the highest probability of being Lyα emitters. When applied to a population in the redshift range z ∈ [2.5, 4.5], our classifier holds a (80 ± 2)% accuracy and (73 ± 4)% precision. At higher redshifts (z ∈ [4.5, 6]), we obtained an accuracy of 73% and precision of 80%. These results highlight the possibility of overcoming the current limitations in assembling large samples of LAEs by making informed predictions that can be used for planning future large-scale spectroscopic surveys.

    FULL ARTICLE

    Insights into the reionization epoch from cosmic-noon-C IV emitters in the VANDELS survey

    Recently, intense emission from nebular C III] and C IV emission lines have been observed in galaxies in the epoch of reionization (z > 6) and have been proposed as the prime way of measuring their redshift and studying their stellar populations. These galaxies might represent the best examples of cosmic reionizers, as suggested by recent low-z observations of Lyman continuum emitting galaxies, but it is hard to directly study the production and escape of ionizing photons at such high redshifts. The ESO spectroscopic public survey VANDELS offers the unique opportunity to find rare examples of such galaxies at cosmic noon (z ∼ 3), thanks to the ultra deep observations available. We have selected a sample of 39 galaxies showing C IV emission, whose origin (after a careful comparison to photoionization models) can be ascribed to star formation and not to active galactic nuclei. By using a multiwavelength approach, we determined their physical properties including metallicity and the ionization parameter and compared them to the properties of the parent population to understand what the ingredients are that could characterize the analogs of the cosmic reionizers. We find that C IV emitters are galaxies with high photon production efficiency and there are strong indications that they might also have a large escape fraction: given the visibility of C IV in the epoch of reionization, this could become the best tool to pinpoint the cosmic reioinzers.

    FULL ARTICLE

    Closing in on the sources of cosmic reionization: First results from the GLASS-JWST program

    The escape fraction of Lyman-continuum (LyC) photons (fesc) is a key parameter for determining the sources of cosmic reionization at z ≥ 6. At these redshifts, owing to the opacity of the intergalactic medium, the LyC emission cannot be measured directly. However, LyC leakers during the epoch of reionization could be identified using indirect indicators that have been extensively tested at low and intermediate redshifts. These include a high [O III]/[O II] flux ratio, high star-formation surface density, and compact sizes. In this work, we present observations of 29 4.5 ≤ z ≤ 8 gravitationally lensed galaxies in the Abell 2744 cluster field. From a combined analysis of JWST-NIRSpec and NIRCam data, we accurately derived their physical and spectroscopic properties: our galaxies have low masses (log(M)∼8.5), blue UV spectral slopes (β ∼ −2.1), compact sizes (re ∼ 0.3 − 0.5 kpc), and high [O III]/[O II] flux ratios. We confirm that these properties are similar to those characterizing low-redshift LyC leakers. Indirectly inferring the fraction of escaping ionizing photons, we find that more than 80% of our galaxies have predicted fesc values larger than 0.05, indicating that they would be considered leakers. The average predicted fesc value of our sample is 0.12, suggesting that similar galaxies at z ≥ 6 have provided a substantial contribution to cosmic reionization.

    FULL ARTICLE

    Reionization with galaxies and active galactic nuclei

    In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ∼1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( ⟨fsfesc⟩ ) and AGN ( fbhesc ) to find: (I) the ionizing budget is dominated by stellar radiation from low stellar mass ( M∗<10^9M⊙ ) galaxies at z > 6 with the AGN contribution (driven by Mbh>10^6M⊙ black holes in M∗≳10^9M⊙ galaxies) dominating at lower redshifts; (II) AGN only contribute 10−25 per cent to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (III) if the stellar mass dependence of ⟨fsfesc⟩ is shallower than fbhesc , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (IV) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ∼ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

    https://ui.adsabs.harvard.edu/abs/2020MNRAS.495.3065D/abstract

    Unveiling early black hole growth with multifrequency gravitational wave observations

    Third-generation ground-based gravitational wave interferometers, like the Einstein Telescope (ET), Cosmic Explorer, and the Laser Interferometer Space Antenna (LISA), will detect coalescing binary black holes over a wide mass spectrum and across all cosmic epochs. We track the cosmological growth of the earliest light and heavy seeds that swiftly transit into the supermassive domain using a semi-analytical model for the formation of quasars at z = 6.4, 2, and 0.2, in which we follow black hole coalescences driven by triple interactions. We find that light-seed binaries of several 102M⊙ are accessible to ET with a signal-to-noise ratio (S/N) of 10-20 at 6 < z < 15. They then enter the LISA domain with larger S/N as they grow to a few 10^4M⊙ . Detecting their gravitational signal would provide first time evidence that light seeds form, grow, and dynamically pair during galaxy mergers. The electromagnetic emission of accreting black holes of similar mass and redshift is too faint to be detected even for the deepest future facilities. ET will be our only chance to discover light seeds forming at cosmic dawn. At 2 < z < 8, we predict a population of ‘starved binaries’, long-lived marginally growing light-seed pairs, to be loud sources in the ET bandwidth (S/N > 20). Mergers involving heavy seeds ( ∼10^5−10^6M⊙ ) would be within reach up to z = 20 in the LISA frequency domain. The lower z model predicts 11.25 (18.7) ET (LISA) events per year, overall.

    https://ui.adsabs.harvard.edu/abs/2021MNRAS.500.4095V/abstract

    Cosmic archaeology with massive stellar black hole binaries

    The existence of massive stellar black hole binaries (MBHBs), with primary black hole (BH) masses ≥ 31 Mo, 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 02 data. By adopting the galaxy formation model 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 < 2 < 10, regardless of the adopted BPS.
    These MBHBs are then old systems forming in low-metallicity (Z~ [0.01-0.1] Zo), 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.

    fuoejfkd;wklp]