No Arabic abstract
We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.
Recent suggestions of a photon underproduction crisis (Kollmeier etal 2014) have generated concern over the intensity and spectrum of ionizing photons in the metagalactic ultraviolet background (UVB). The balance of hydrogen photoionization and recombination determines the opacity of the low-redshift intergalactic medium (IGM). We calibrate the hydrogen photoionization rate ($Gamma_{rm H}$) by comparing {it Hubble Space Telescope} spectroscopic surveys of the low-redshift column density distribution of HI absorbers and the observed ($z < 0.4$) mean Lya flux decrement, $D_A = (0.014)(1+z)^{2.2}$, to new cosmological simulations. The distribution, $f(N_{rm HI}, z) equiv d^2 {cal N} / d(log N_{rm HI}) dz$, is consistent with an increased UVB that includes contributions from both quasars and galaxies. Our recommended fit, $Gamma_{rm H}(z) = (4.6 times 10^{-14}$ s$^{-1})(1+z)^{4.4}$ for $0 < z < 0.47$, corresponds to unidirectional LyC photon flux $Phi_0 approx 5700$~cm$^{-2}$~s$^{-1}$ at $z = 0$. This flux agrees with observed IGM metal ionization ratios (CIII/CIV and SiIII/SiIV) and suggests a 25-30% contribution of Lya absorbers to the cosmic baryon inventory. The primary uncertainties in the low-redshift UVB are the contribution from massive stars in galaxies and the LyC escape fraction ($f_{rm esc}$), a highly directional quantity that is difficult to constrain statistically. We suggest that both quasars and low-mass starburst galaxies are important contributors to the ionizing UVB at $z < 2$. Their additional ionizing flux would resolve any crisis in photon underproduction.
We demonstrate a new method for measuring the escape fraction of ionizing photons using Hubble Space Telescope imaging of resolved stars in NGC 4214, a local analog of high-redshift starburst galaxies that are thought to be responsible for cosmic reionization. Specifically, we forward model the UV through near-IR spectral energy distributions of $sim$83,000 resolved stars to infer their individual ionizing flux outputs. We constrain the local escape fraction by comparing the number of ionizing photons produced by stars to the number that are either absorbed by dust or consumed by ionizing the surrounding neutral hydrogen in individual star-forming regions. We find substantial spatial variation in the escape fraction (0-40%). Integrating over the entire galaxy yields a global escape fraction of 25% (+16%/-15%). This value is much higher than previous escape fractions of zero reported for this galaxy. We discuss sources of this apparent tension, and demonstrate that the viewing angle and the 3D ISM geometric effects are the cause. If we assume the NGC 4214 has no internal dust, like many high-redshift galaxies, we find an escape fraction of 59% (an upper-limit for NGC 4214). This is the first non-zero escape fraction measurement for UV-faint (M$_{rm FUV}$ = -15.9) galaxies at any redshift, and supports the idea that starburst UV-faint dwarf galaxies can provide a sufficient amount of ionizing photons to the intergalactic medium.
Recent observations show that the measured rates of star formation in the early universe are insufficient to produce re-ionization, and therefore, another source of ionizing photons is required. In this emph{Letter}, we examine the possibility that these can be supplied by the fast accretion shocks formed around the cores of the most massive haloes ($10.5< log M/M_{odot} < 12$) on spatial scales of order 1 kpc. We model the detailed physics of these fast accretion shocks, and apply these to a simple 1-D spherical hydrodynamic accretion model for baryonic infall in dark matter halos with an Einasto density distribution. The escape of UV photons from these halos is delayed by the time taken to reach the critical accretion shock velocity for escape of UV photons; 220 km s$^{-1}$, and by the time it takes for these photons to ionize the surrounding baryonic matter in the accretion flow. Assuming that in the universe at large the baryonic matter tracks the dark matter, we can estimate the epoch of re-ionization in the case that accretion shocks act alone as the source of UV photons. We find that 50% of the volume (and 5-8% of the mass) of the universe can be ionized by $z sim 7-8$. The UV production rate has an uncertainty of a factor of about 5 due to uncertainties in the cosmological parameters controlling the development of large scale structure. Because our mechanism is a steeply rising function of decreasing redshift, this uncertainty translates to a re-ionization redshift uncertainty of less than $pm0.5$. We also find that, even without including the UV photon production of stars, re-ionization is essentially complete by $z sim 5.8$. Thus, fast accretion shocks can provide an important additional source of ionizing photons in the early universe.
Using broadband photometry from the Hubble Space Telescope in combination with Very Large Telescope narrowband Halpha observations of the nearby spiral galaxy NGC 300, we explore a method for estimating the escape fractions of hydrogen-ionizing photons from HII regions within this galaxy. Our goal in this concept study is to evaluate the spectral types of the most massive stars using the broadband data and estimating their ionizing photon output with the help of stellar atmosphere models. A comparison with the Halpha flux that gives the amount of ionized gas in the HII region provides a measure of the escape fraction of ionizing photons from that region. We performed some tests with a number of synthetic young clusters with varying parameters to assess the reliability of the method. However, we found that the derived stellar spectral types and consequently the expected ionizing photon luminosity of a region is highly uncertain. The tests also show that on one hand we tended to overestimate the integrated photon output of a region for young ages and low numbers of stars, and on the other hand we mostly underestimated the combined ionizing luminosity for a large stellar number and older cluster ages. We conclude that the proposed method of using stellar broadband photometry to infer the leakage of ionizing photons from HII regions is highly uncertain and dominated by the errors of the resulting stellar spectral types. Therefore this method is not suitable. Stellar spectra are needed to reliably determine the stellar types and escape fractions. Studies to this end have been carried out for the Magellanic Clouds.
We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=5.7 and 6.6 in a 0.65 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field, using deep images taken with Subaru/Suprime-Cam, UKIRT/WFCAM, and Spitzer/IRAC. We produce stacked multiband images at each redshift from 165 (z=5.7) and 91 (z=6.6) IRAC-undetected objects, to derive typical spectral energy distributions (SEDs) of z~6-7 LAEs for the first time. The stacked LAEs have as blue UV continua as the HST/WFC3 z-dropout galaxies of similar Muv, with a spectral slope beta ~ -3, but at the same time they have red UV-to-optical colors with detection in the 3.6um band. Using SED fitting we find that the stacked LAEs have low stellar masses of ~(3-10)*10^7 Msun, very young ages of ~1-3 Myr, negligible dust extinction, and strong nebular emission from the ionized interstellar medium, although the z=6.6 object is fitted similarly well with high-mass models without nebular emission; inclusion of nebular emission reproduces the red UV-to-optical color while keeping the UV color sufficiently blue. We infer that typical LAEs at z~6-7 are building blocks of galaxies seen at lower redshifts. We find a tentative decrease in the Lyman alpha escape fraction from z=5.7 to 6.6, which may imply an increase in the intergalactic medium neutral fraction. From the minimum contribution of nebular emission required to fit the observed SEDs, we place an upper limit on the escape fraction of ionizing photons to be f_esc^ion~0.6 at z=5.7 and ~0.9 at z=6.6. We also compare the stellar populations of our LAEs with that of stacked HST/WFC3 z-dropout galaxies.