No Arabic abstract
Recent observations of galaxies at $z gtrsim 7$, along with the low value of the electron scattering optical depth measured by the Planck mission, make galaxies plausible as dominant sources of ionizing photons during the epoch of reionization. However, scenarios of galaxy-driven reionization hinge on the assumption that the average escape fraction of ionizing photons is significantly higher for galaxies in the reionization epoch than in the local Universe. The NIRSpec instrument on the James Webb Space Telescope (JWST) will enable spectroscopic observations of large samples of reionization-epoch galaxies. While the leakage of ionizing photons will not be directly measurable from these spectra, the leakage is predicted to have an indirect effect on the spectral slope and the strength of nebular emission lines in the rest-frame ultraviolet and optical. Here, we apply a machine learning technique known as lasso regression on mock JWST/NIRSpec observations of simulated $z=7$ galaxies in order to obtain a model that can predict the escape fraction from JWST/NIRSpec data. Barring systematic biases in the simulated spectra, our method is able to retrieve the escape fraction with a mean absolute error of $Delta f_{mathrm{esc}} approx 0.12$ for spectra with $S/Napprox 5$ at a rest-frame wavelength of 1500 {AA} for our fiducial simulation. This prediction accuracy represents a significant improvement over previous similar approaches.
We describe a new method for simulating ionizing radiation and supernova feedback in the analogues of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3e11 Msun, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ~20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60 degree from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f_esc(i), and how it evolves as the particle ages. We discover that the average escape fraction f_esc is dominated by a small number of SFMC particles with high f_esc(i). On average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump, and from a galactic disk.
In this paper we calculate the escape fraction ($f_{rm esc}$) of ionizing photons from starburst galaxies. Using 2-D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the center of the disk along different angles through the superbubble and the gas disk. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of $sim 40 ^circ$, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scale heights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed $sim [1- cos (1 , {rm radian})] = 0.5$ from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time- and angle-averaged escape fraction on the mid-plane disk gas density (in the range $n_0=0.15-50$ cm $^{-3}$) and the disk scale height (between $z_0=10-600$ pc). We find that the escape fraction is related to the disk parameters (the mid-plane disk density and scale height) roughly so that $f_{rm esc}^alpha n_0^2 z_0^3$ (with $alphaapprox 2.2$) is a constant. For disks with a given WNM temperature, massive disks have lower escape fraction than low mass galaxies. For Milky Way ISM parameters, we find $f_{rm esc}sim 5%$, and it increases to $approx 10%$ for a galaxy ten times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the universe.
The HI gas content is a key ingredient in galaxy evolution, the study of which has been limited to moderate cosmological distances for individual galaxies due to the weakness of the hyperfine HI 21-cm transition. Here we present a new approach that allows us to infer the HI gas mass $M_{rm HI}$ of individual galaxies up to $zapprox 6$, based on a direct measurement of the [CII]-to-HI conversion factor in star-forming galaxies at $zgtrsim 2$ using $gamma$-ray burst afterglows. By compiling recent [CII]-158 $mu$m emission line measurements we quantify the evolution of the HI content in galaxies through cosmic time. We find that the HI mass starts to exceed the stellar mass $M_star$ at $zgtrsim 1$, and increases as a function of redshift. The HI fraction of the total baryonic mass increases from around $20%$ at $z = 0$ to about $60%$ at $zsim 6$. We further uncover a universal relation between the HI gas fraction $M_{rm HI}/M_star$ and the gas-phase metallicity, which seems to hold from $zapprox 6$ to $z=0$. The majority of galaxies at $z>2$ are observed to have HI depletion times, $t_{rm dep,HI} = M_{rm HI}/{rm SFR}$, less than $approx 2$ Gyr, substantially shorter than for $zsim 0$ galaxies. Finally, we use the [CII]-to-HI conversion factor to determine the cosmic mass density of HI in galaxies, $rho_{rm HI}$, at three distinct epochs: $zapprox 0$, $zapprox 2$, and $zsim 4-6$. These measurements are consistent with previous estimates based on 21-cm HI observations in the local Universe and with damped Lyman-$alpha$ absorbers (DLAs) at $zgtrsim 2$, suggesting an overall decrease by a factor of $approx 5$ in $rho_{rm HI}(z)$ from the end of the reionization epoch to the present.
Using results from high-resolution galaxy formation simulations in a standard Lambda-CDM cosmology and a fully conservative multi-resolution radiative transfer code around point sources, we compute the energy-dependent escape fraction of ionizing photons from a large number of star forming regions in two galaxies at five different redshifts from z=3.8 to 2.39. All escape fractions show a monotonic decline with time, from (at the Lyman-limit) ~6-10% at z=3.6 to ~1-2% at z=2.39, due to higher gas clumping at lower redshifts. It appears that increased feedback can lead to higher f_esc at z>3.4 via evacuation of gas from the vicinity of star forming regions and to lower f_esc at z<2.39 through accumulation of swept-up shells in denser environments. Our results agree well with the observational findings of citet{inoue..06} on redshift evolution of f_esc in the redshift interval z=2-3.6.
We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.