No Arabic abstract
Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly analogues of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $zsim3$ and $zsim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $zsim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{rm esc}$, lower-metallicity LyC leakers at $zsim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{rm 158mu m}$ and [OIII]$_{rm 88mu m}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.
An analysis of the dynamics of a star formation event is performed. It is shown that galaxies able to drive leftover gas to sufficient altitudes in a few million years are characterized by two basic properties: small sizes (<1kpc) and high star formation rate surface densities (Sigma_SFR > 10 Msun/yr/kpc2). For the parameter space of relevance, the outflow is primarily driven by supernovae with radiation pressure being significant but subdominant. Our analysis provides the unifying physical origin for a diverse set of observed LyC leakers, including the green-peas galaxies, [SII]-weak galaxies, Lyman-alpha emitters, with these two characteristics as the common denominator. Among verifiable physical properties of LyC leakers, we predict that (1) the newly formed stellar masses are are typically in the range of 1e8-1e10 Msun, except perhaps ULIRGs, (2) the outflow velocities are typically in the range typically of 100-600km/s, but may exceed 1e3 km/s in ULIRGs, with a strong positive correlation between the stellar masses formed and the outflow velocities, (3) the overall escape fraction of galaxies is expected to increase with increasing redshift, given the cosmological trend that galaxies become denser and more compact with increasing redshift. In addition, two interesting by-product predictions are also borne out. First, ULIRGs appear to be in a parameter region where they should be prodigious LyC leakers, unless there is a large ram-pressure. Second, Lyman break galaxies (LBGs) are not supposed to be prodigious LyC leakers in our model, given their claimed effective radii exceeding 1kpc.
Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.
We present $U_{336}V_{606}J_{125}H_{160}$ follow-up $HST$ observations of 16 $zsim3$ candidate LyC emitters in the HS1549+1919 field. With these data, we obtain high spatial-resolution photometric redshifts of all sub-arcsecond components of the LyC candidates in order to eliminate foreground contamination and identify robust candidates for leaking LyC emission. Of the 16 candidates, we find one object with a robust LyC detection that is not due to foreground contamination. This object (MD5) resolves into two components; we refer to the LyC-emitting component as MD5b. MD5b has an observed 1500AA to 900AA flux-density ratio of $(F_{UV}/F_{LyC})_{obs}=4.0pm2.0$, compatible with predictions from stellar population synthesis models. Assuming minimal IGM absorption, this ratio corresponds to a relative (absolute) escape fraction of $f_{esc,rel}^{MD5b}=75-100$% ($f_{esc,abs}^{MD5b}=14-19$%). The stellar population fit to MD5b indicates an age of $lesssim50$Myr, which is in the youngest 10% of the $HST$ sample and the youngest third of typical $zsim3$ Lyman break galaxies, and may be a contributing factor to its LyC detection. We obtain a revised, contamination-free estimate for the comoving specific ionizing emissivity at $z=2.85$, indicating (with large uncertainties) that star-forming galaxies provide roughly the same contribution as QSOs to the ionizing background at this redshift. Our results show that foreground contamination prevents ground-based LyC studies from obtaining a full understanding of LyC emission from $zsim3$ star-forming galaxies. Future progress in direct LyC searches is contingent upon the elimination of foreground contaminants through high spatial-resolution observations, and upon acquisition of sufficiently deep LyC imaging to probe ionizing radiation in high-redshift galaxies.
We have recently discovered five Lyman continuum leaking galaxies at z~0.3, selected for their compactness, intense star-formation, and high [OIII]/[OII] ratio (Izotov et al. 2016ab). Here we derive their ionizing photon production efficiency, a fundamental quantity for inferring the number of photons available to reionize the Universe, for the first time for galaxies with confirmed strong Lyman continuum escape (fesc~6-13%). We find an ionizing photon production per unit UV luminosity, which is a factor 2-6 times higher than the canonical value when reported to their observed UV luminosity. After correction for extinction this value is close to the canonical value. The properties of our five Lyman continuum leakers are found to be very similar to those of the confirmed z=3.218 leaker Ion2 from de Barros et al. (2016) and very similar to those of typical star-forming galaxies at z>~6. Our results suggest that UV bright galaxies at high-z such as Lyman break galaxies can be Lyman continuum leakers and that their contribution to cosmic reionization may be underestimated.
We calculate Lyman Alpha Emitter (LAE) angular correlation functions (ACFs) at $z simeq 6.6$ and the fraction of lifetime (for the 100 Myrs preceding $zsimeq6.6$) galaxies spend as Lyman Break Galaxies (LBGs) or as LBGs with Lyman Alpha (Ly$alpha$) emission using a model that combines SPH cosmological simulations (GADGET-2), dust attenuation and a radiative transfer code (pCRASH). The ACFs are a powerful tool that significantly narrows the 3D parameter space allowed by LAE Ly$alpha$ and UV luminosity functions (LFs) alone. With this work, we simultaneously constrain the escape fraction of ionizing photons $f_{esc}=0.05-0.5$, the mean fraction of neutral hydrogen in the intergalactic medium (IGM) $langle chi_{HI} rangle leq 0.01$ and the dust-dependent ratio of the escape fractions of Ly$alpha$ and UV continuum photons $f_{alpha}/f_c=0.6-1.2$. Our results show that reionization has the largest impact on the amplitude of the ACFs, and its imprints are clearly distinguishable from those of $f_{esc}$ and $f_alpha/f_c$. We also show that galaxies with a critical stellar mass of $M_* = 10^{8.5} (10^{9.5})M_{odot}$ produce enough luminosity to stay visible as LBGs (LAEs). Finally, the fraction of time during the past 100 Myrs prior to z=6.6 a galaxy spends as a LBG or as a LBG with Ly$alpha$ emission increases with the UV magnitude (and the stellar mass $M_*$): considering observed (dust and IGM attenuated) luminosities, the fraction of time a galaxy spends as a LBG (LAE) increases from 65% to 100% (0-100%) as $M_{UV}$ decreases from $M_{UV} = -18.0$ to $-23.5$ ($M_*$ increases from $10^8-10^{10.5} M_{odot}$). Thus in our model the brightest (most massive) LBGs most often show Ly$alpha$ emission.