Search for: All records

Abstract Detecting the first generation of stars, Population III (Pop III), has been a long-standing goal in astrophysics, yet they remain elusive even in the JWST era. Here we present a novel NIRCam-based selection method for Pop III galaxies, and carefully validate it through completeness and contamination simulations. We systematically search ≃ 500 arcmin²across JWST legacy fields for Pop III candidates, including GLIMPSE, which, assisted by gravitational lensing, has produced JWST’s deepest NIRCam imaging thus far. We discover one promising Pop III galaxy candidate (GLIMPSE-16043) at $z=6.50_{-0.24}^{+0.03}$, a moderately lensed galaxy ( $\mu =2.9_{-0.2}^{+0.1}$) with an intrinsic UV magnitude of $M_{\mathrm{UV}}=-15.89_{-0.14}^{+0.12}$. It exhibits key Pop III features: strong Hαemission (rest-frame EW 2810 ± 550 Å); a Balmer jump; no dust (UV slopeβ = −2.34 ± 0.36); and undetectable metal lines (e.g., [Oiii]; [Oiii]/Hβ < 0.44), implying a gas-phase metallicity ofZ_gas/Z_⊙ < 0.5%. These properties indicate the presence of a nascent, metal-deficient young stellar population (<5 Myr) with a stellar mass of ≃10⁵M_⊙. Intriguingly, this source deviates significantly from the extrapolated UV–metallicity relation derived from recent JWST observations atz= 4–10, consistent with UV enhancement by a top-heavy Pop III initial mass function or the presence of an extremely metal-poor active galactic nucleus. We also derive the first observational constraints on the Pop III UV luminosity function atz ≃ 6–7. The volume density of GLIMPSE-16043 (≈10⁻⁴cMpc⁻³) is in excellent agreement with theoretical predictions, independently reinforcing its plausibility. This study demonstrates the power of our novel NIRCam method to finally reveal distant galaxies even more pristine than the Milky Way’s most metal-poor satellites, thereby promising to bring us closer to the first generation of stars than we have ever been before. 

Abstract A key obstacle to accurate models of the first stars and galaxies is the vast range of distance scales that must be considered. While star formation occurs on sub-parsec scales within dark matter (DM) minihalos, it is influenced by large-scale baryon-dark matter streaming velocities (v_bc) and Lyman-Werner (LW) radiative feedback which vary significantly on scales of ∼100 Mpc. We present a novel approach to this issue in which we utilize artificial neural networks (NNs) to emulate the Population III (PopIII) and Population II (PopII) star formation histories of many small-scale cells given by a more complex semi-analytic framework based on DM halo merger trees. Within each simulation cell, the NN takes a set of input parameters that depend on the surrounding large-scale environment, such as the cosmic overdensity,δ(x⃗), andv_bcof the cell, then outputs the resulting star formation far more efficiently than is possible with the semi-analytic model. This rapid emulation allows us to self-consistently determine the LW background intensity on ∼100 Mpc scales, while simultaneously including the detailed merger histories (and corresponding star formation histories) of the low-mass minihalos that host the first stars. Comparing with the full semi-analytic framework utilizing DM halo merger trees, our NN emulators yield star formation histories with redshift-averaged errors of ∼7.3% and ∼5.2% for PopII and PopIII, respectively. When compared to a simpler sub-grid star formation prescription reliant on halo mass function integration, we find that the diversity of halo merger histories in our simulation leads to enhanced spatial fluctuations, an earlier transition from PopIII to PopII dominated star formation, and more scatter in star formation histories overall. 

Abstract Long troughs observed in thez> 5.5 Lyαand Lyβforests are thought to be caused by the last remaining neutral patches during the end phases of reionization — termed neutral islands. If this is true, then the longest troughs mark locations where we are most likely to observe the reionizing intergalactic medium (IGM). A key feature of the neutral islands is that they are bounded by ionization fronts (I-fronts) which emit Lyman series lines. In this paper, we explore the possibility of directly imaging the outline of neutral islands with a narrowband survey targeting Lyα. In a companion paper, we quantified the intensity of I-front Lyαemissions during reionization and its dependence on the spectrum of incident ionizing radiation and I-front speed. Here we apply those results to reionization simulations to model the emissions from neutral islands. We find that neutral islands would appear as diffuse structures that are tens of comoving Mpc across, with surface brightnesses in the range ≈ 1 - 5× 10^-21erg s^-1cm^-2arcsec^-2. The islands are brighter if the spectrum of ionizing radiation driving the I-fronts is harder, and/or if the I-fronts are moving faster. We develop mock observations for current and futuristic observatories and find that, while extremely challenging, detecting neutral islands is potentially within reach of an ambitious observing program with wide-field narrowband imaging. Our results demonstrate the potentially high impact of low-surface brightness observations for studying reionization. 

Abstract During reionization, intergalactic ionization fronts (I-fronts) are sources of Lyαline radiation produced by collisional excitation of hydrogen atoms within the fronts. In principle, detecting this emission could provide direct evidence for a reionizing intergalactic medium (IGM). In this paper, we use a suite of high-resolution one-dimensional radiative transfer simulations run on cosmological density fields to quantify the parameter space of I-front Lyαemission. We find that the Lyαproduction efficiency — the ratio of emitted Lyαflux to incident ionizing flux driving the front — depends mainly on the I-front speed and the spectral index of the ionizing radiation. IGM density fluctuations on scales smaller than the typical I-front width produce scatter in the efficiency, but they do not significantly boost its mean value. The Lyαflux emitted by an I-front is largest if 3 conditions are met simultaneously: (1) the incident ionizing flux is large; (2) the incident spectrum is hard, consisting of more energetic photons; (3) the I-front is traveling through a cosmological over-density, which causes it to propagate more slowly. We present a convenient parameterization of the efficiency in terms of I-front speed and incident spectral index. We make these results publicly available as an interpolation table and we provide a simple fitting function for a representative ionizing background spectrum. Our results can be applied as a sub-grid model for I-front Lyα emissions in reionization simulations with spatial and/or temporal resolutions too coarse to resolve I-front structure. In a companion paper, we use our results to explore the possibility of directly imaging Lyαemission around neutral islands during the last phases of reionization. 

Abstract We present a method that calibrates a semianalytic model to the Renaissance Simulations, a suite of cosmological hydrodynamical simulations with high-redshift galaxy formation. This approach combines the strengths of semianalytic techniques and hydrodynamical simulations, enabling the extension to larger volumes and lower redshifts that are inaccessible to simulations due to computational expense. Using a sample of Renaissance star formation histories from an average density region of the Universe, we construct a four-parameter prescription for metal-enriched star formation characterized by an initial bursty stage followed by a steady stage where stars are formed at constant efficiencies. Our model also includes a treatment of Pop III star formation where a minimum halo mass and log-normal distribution of stellar mass are adopted to match the numerical simulations. Star formation is generally well reproduced for halos with masses ≲10⁹M_⊙. Between 11 <z< 25 our model produces metal-enriched star formation rate densities (SFRDs) that typically agree with Renaissance within a factor of ∼2 for the average density region. Additionally, the total metal-enriched stellar mass only differs from Renaissance by about 10% atz∼ 11. For regions that are either more overdense or rarefied and not included in the calibration, we produce metal-enriched SFRDs that agree with Renaissance within a factor of ∼2 at high-zbut eventually differ by higher factors for later times. This is likely due to environmental dependencies not included in the model. Our star formation prescriptions can easily be adopted in other analytic or semianalytic works to match our calibration to Renaissance. 

Abstract We present a new self-consistent semianalytic model of the first stars and galaxies to explore the high-redshift (z≥ 15) Population III (PopIII) and metal-enriched star formation histories. Our model includes the detailed merger history of dark matter halos generated with Monte Carlo merger trees. We calibrate the minimum halo mass for PopIII star formation from recent hydrodynamical cosmological simulations that simultaneously include the baryon–dark matter streaming velocity, Lyman–Werner (LW) feedback, and molecular hydrogen self-shielding. We find an overall increase in the resulting star formation rate density (SFRD) compared to calibrations based on previous simulations (e.g., the PopIII SFRD is over an order of magnitude higher atz= 35−15). We evaluate the effect of the halo-to-halo scatter in this critical mass and find that it increases the PopIII stellar mass density by a factor ∼1.5 atz≥ 15. Additionally, we assess the impact of various semianalytic/analytic prescriptions for halo assembly and star formation previously adopted in the literature. For example, we find that models assuming smooth halo growth computed via abundance matching predict SFRDs similar to the merger tree model for our fiducial model parameters, but that they may underestimate the PopIII SFRD in cases of strong LW feedback. Finally, we simulate subvolumes of the Universe with our model both to quantify the reduction in total star formation in numerical simulations due to a lack of density fluctuations on spatial scales larger than the simulation box, and to determine spatial fluctuations in SFRD due to the diversity in halo abundances and merger histories. 

Abstract Line-intensity mapping (IM) experiments seek to perform statistical measurements of large-scale structure with spectral lines such as 21 cm, CO, and Lyα. A challenge in these observations is to ensure that astrophysical foregrounds, such as galactic synchrotron emission in 21 cm measurements, are properly removed. One method that has the potential to reduce foreground contamination is to cross correlate with a galaxy survey that overlaps with the IM volume. However, telescopes sensitive to high-redshift galaxies typically have small field of views compared to IM surveys. Thus, a galaxy survey for cross correlation would necessarily consist of pencil beams that sparsely fill the IM volume. In this paper, we develop the formalism to forecast the sensitivity of cross correlations between IM experiments and pencil-beam galaxy surveys. We find that a random distribution of pencil beams leads to very similar overall sensitivity as a lattice spaced across the IM survey and derive a simple formula for random configurations that agrees with the Fisher matrix formalism. We explore examples of combining high-redshift James Webb Space Telescope (JWST) observations with both an SPHEREx-like LyαIM survey and a 21 cm experiment based on the Hydrogen Epoch of Reionization Array (HERA). We find that the JWST-SPHEREx case is promising, leading to a total signal-to-noise ratio of ∼5 after 100 total hours of JWST (atz= 7). We find that HERA is not well-suited for this approach owing to its drift-scan strategy, but that a similar experiment that can integrate down on one field could be. 

Abstract We discover analytic equations that can infer the value of Ω_mfrom the positions and velocity moduli of halo and galaxy catalogs. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from GadgetN-body simulations to perform field-level likelihood-free inference, and show that our model can infer Ω_mwith ∼6% accuracy from halo catalogs of thousands ofN-body simulations run with six different codes: Abacus, CUBEP³M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict Ω_mfrom halo catalogs of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that, by tuning a single free parameter, our equations can also infer the value of Ω_mfrom galaxy catalogs of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogs from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and Ω_m, one that is not affected by galaxy formation physics down to scales as small as 10h⁻¹kpc. 

Abstract Fuzzy dark matter (FDM) is a proposed modification for the standard cold dark matter (CDM) model motivated by small-scale discrepancies in low-mass galaxies. Composed of ultralight (mass ∼ 10²²eV) axions with kiloparsec-scale de Broglie wavelengths, this is one of a class of candidates that predicts that the first collapsed objects form in relatively massive dark matter halos. This implies that the formation history of the first stars and galaxies would be very different, potentially placing strong constraints on such models. Here we numerically simulate the formation of the first stars in an FDM cosmology, following the collapse in a representative volume all the way down to primordial protostar formation including a primordial nonequilibrium chemical network and cooling for the first time. We find two novel results: first, the large-scale collapse results in a very thin and flat gas “pancake”; second, despite the very different cosmology, this pancake fragments until it forms protostellar objects indistinguishable from those in CDM. Combined, these results indicate that the first generation of stars in this model are also likely to be massive and, because of the sheet morphology, do not self-regulate, resulting in a massive Population III starburst. We estimate the total number of first stars forming in this extended structure to be 10⁴over 20 Myr using a simple model to account for the ionizing feedback from the stars, and should be observable with the James Webb Space Telescope. These predictions provide a potential smoking gun signature of FDM and similar dark matter candidates. 

ABSTRACT Dark matter haloes that reach the H i-cooling mass without prior star formation or external metal pollution represent potential sites for the formation of small – extremely faint – Population III galaxies at high redshifts. Gravitational lensing may in rare cases boost their fluxes to detectable levels, but to find even a small number of such objects in randomly selected regions of the sky requires very large areas to be surveyed. Because of this, a small, wide-field telescope can in principle offer better detection prospects than a large telescope with a smaller field of view. Here, we derive the minimum comoving number density required to allow gravitational lensing to lift such objects at redshift z = 5−16 above the detection thresholds of blind surveys carried out with the James Webb space telescope (JWST), the Roman space telescope (RST) and Euclid. We find that the prospects for photometric detections of Pop III galaxies are promising, and that they are better for RST than for JWST and Euclid. However, the Pop III galaxies favoured by current simulations have number densities too low to allow spectroscopic detections based on the strength of the He ii1640 emission line in any of the considered surveys unless very high star formation efficiencies (ϵ ≳ 0.1) are evoked. We argue that targeting individual cluster lenses instead of the wide-field surveys considered in this paper results in better spectroscopic detection prospects, while for photometric detection, the wide-field surveys perform considerably better.