No Arabic abstract
Population III stars can regulate star formation in the primordial Universe in several ways. They can ionize nearby halos, and even if their ionizing photons are trapped by their own halos, their Lyman-Werner (LW) photons can still escape and destroy H$_2$ in other halos, preventing them from cooling and forming stars. LW escape fractions are thus a key parameter in cosmological simulations of early reionization and star formation but have not yet been parametrized for realistic halos by halo or stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV escape from 9--120 M$_{odot}$ Pop III stars in $10^5$ to $10^7$ M$_{odot}$ halos with ZEUS-MP. We find that photons in the LW lines (i.e. those responsible for destroying H$_{2}$ in nearby systems) have escape fractions ranging from 0% to 85%. No LW photons escape the most massive halo in our sample, even from the most massive star. Escape fractions for photons elsewhere in the 11.18--13.6~eV energy range, which can be redshifted into the LW lines at cosmological distances, are generally much higher, being above 60% for all but the least massive stars in the most massive halos. We find that shielding of H$_2$ by neutral hydrogen, which has been neglected in most studies to date, produces escape fractions that are up to a factor of three smaller than those predicted by H$_2$ self-shielding alone.
Direct collapse black holes forming in pristine, atomically-cooling haloes at $z approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.
Pristine, atomically-cooled haloes are leading contenders for the sites of primordial quasar formation because atomic cooling triggers rapid baryon collapse that can create 10$^4$ - 10$^5$ M$_{odot}$ black hole seeds. However, until now no numerical simulations with a wide range of halo spins and assembly histories have followed the collapse for the times required to form a black hole. We have now performed cosmological simulations of baryon collapse in atomically-cooled haloes for times that are sufficient for supermassive stars to form and die as direct-collapse black holes (DCBHs). Our simulations reveal that fragmentation of the accretion disk at the center of the halo after $sim$ 500 kyr is nearly ubiquitous and in most cases leads to the formation of binary or multiple supermassive stellar systems. They also confirm that rapid baryon collapse proceeds for the times required for these stars to form DCBHs. Our discovery raises the exciting possibility of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the James Webb Space Telescope and ground-based telescopes in the coming decade.
The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$alpha$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($Sigma$), we propose a model relating $f_{esc}$ to $Sigma$ and find $f_{esc}proptoSigma^{0.4pm0.1}$. Since $Sigma$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs) account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($alpha$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)
We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.
Identifying the mechanisms driving the escape of Lyman Continuum (LyC) photons is crucial to find Lyman Continuum Emitter (LCE) candidates. To understand the physical properties involved in the leakage of LyC photons, we investigate the connection between the HI covering fraction, HI velocity width, the Lyman alpha (LyA) properties and escape of LyC photons in a sample of 22 star-forming galaxies including 13 LCEs. We fit the stellar continua, dust attenuation, and absorption lines between 920 and 1300 A to extract the HI covering fractions and dust attenuation. Additionally, we measure the HI velocity widths of the optically thick Lyman series and derive the LyA equivalent widths (EW), escape fractions (fesc), peak velocities and fluxes at the minimum of the LyA profiles. Overall, we highlight strong correlations between the presence of low HI covering fractions and (1) low LyA peak velocities; (2) more flux at the profile minimum; and (3) larger EW(LyA), fesc(LyA), and fesc(LyC). Hence, low column density channels are crucial ISM ingredients for the leakage of LyC and LyA photons. Additionally, galaxies with narrower HI absorption velocity widths have higher LyA equivalent widths, larger LyA escape fractions, and lower LyA peak velocity separations. This suggests that these galaxies have low HI column density. Finally, we find that dust regulates the amount of LyA and LyC radiation that actually escapes the ISM. Overall, the ISM porosity is one origin of strong LyA emission and enables the escape of ionizing photons in low-z leakers. However, this is not enough to explain the largest fesc(LyC) observed, which indicates that the most extreme LCEs are likely density-bounded along all lines of sight to the observer. Overall, the neutral gas porosity constrains a lower limit to the escape fraction of LyC and LyA photons, providing a key estimator of the leakage of ionizing photons.