Search for: All records

ABSTRACT Population III stars are possible precursors to early supermassive black holes (BHs). The presence of soft UV Lyman–Werner (LW) background radiation can suppress Population III star formation in minihaloes and allow them to form in pristine atomic-cooling haloes. In the absence of molecular hydrogen ($$\rm H_2$$) cooling, atomic-cooling haloes enable rapid collapse with suppressed fragmentation. High background LW fluxes from preceding star-formation have been proposed to dissociate $$\rm H_2$$. This flux can be supplemented by LW radiation from one or more Population III star(s) in the same halo, reducing the necessary background level. Here, we consider atomic-cooling haloes in which multiple protostellar cores form close to one another nearly simultaneously. We assess whether the first star’s LW radiation can dissociate nearby $$\rm H_2$$, enabling rapid accretion on to a nearby protostellar core, and the prompt formation of a second, supermassive star (SMS) from warm, atomically-cooled gas. We use a set of hydrodynamical simulations with the code enzo, with identical LW backgrounds centred on a halo with two adjacent collapsing gas clumps. When an additional large local LW flux is introduced, we observe immediate reductions in both the accretion rates and the stellar masses that form within these clumps. While the LW flux reduces the $$\text{H}_2$$ fraction and increases the gas temperature, the halo core’s potential well is too shallow to promptly heat the gas to $$\gtrsim$$1000 K and increase the second protostar’s accretion rate. We conclude that this internal LW feedback scenario is unlikely to facilitate SMS or massive BH seed formation. 

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 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 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 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. 

The recent detections of a large number of candidate active galactic nuclei at high redshift (i.e.  $z\gtrsim 4$) has increased speculation that heavy seed massive black hole formation may be a required pathway. Here we re-implement the so-called Lyman-Werner (LW) channel model of Dijkstra et al. (2014) to calculate the expected number density of massive black holes formed through this channel. We further enhance this model by extracting information relevant to the model from the $\text{𝚁𝚎𝚗𝚊𝚒𝚜𝚜𝚊𝚗𝚌𝚎}$simulation suite. $\text{𝚁𝚎𝚗𝚊𝚒𝚜𝚜𝚊𝚗𝚌𝚎}$is a high-resolution suite of simulations ideally positioned to probe the high- $z$universe. Finally, we compare the LW-only channel against other models in the literature. We find that the LW-only channel results in a peak number density of massive black holes of approximately at $z\sim 10$. Given the growth requirements and the duty cycle of active galactic nuclei, this means that the LW-only is likely incompatible with recent JWST measurements and can, at most, be responsible for only a small subset of high- $z$active galactic nuclei. Other models from the literature (e.g. rapid assembly; relative velocities between baryons and dark matter) seem therefore better positioned, at present, to explain the high frequency of massive black holes at high $z$. 

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 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 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 We investigate the effect of dark stars (DSs) on the reionization history of the universe, and the interplay between them and feedback due to Lyman–Werner (LW) radiation in reducing the cosmic microwave background (CMB) optical depth to a value within the τ = 0.054 ± 0.007 range measured by Planck. We use a semianalytic approach to evaluate reionization histories and CMB optical depths, which includes Population II stars in atomic cooling halos and Population III stars in minihalos with LW feedback, preceded by a DS phase. We show that while LW feedback by itself can reduce the integrated optical depth to the last scattering surface to ∼0.05 only if the Population III star formation efficiency is less than ∼0.2%, the inclusion of a population of DSs can naturally lead to the measured CMB optical depth for much larger Population III star formation efficiencies ≳1%.