No Arabic abstract
We quantify the distribution of [OIII]+H$beta$ line strengths at z$simeq$7 using a sample of 20 bright (M$_{mathrm{UV}}$ $lesssim$ $-$21) galaxies. We select these systems over wide-area fields (2.3 deg$^2$ total) using a new colour-selection which precisely selects galaxies at z$simeq$6.63$-$6.83, a redshift range where blue Spitzer/IRAC [3.6]$-$[4.5] colours unambiguously indicate strong [OIII]$+$H$beta$ emission. These 20 galaxies suggest a log-normal [OIII]$+$H$beta$ EW distribution with median EW = 759$^{+112}_{-113}$ $mathrm{mathring{A}}$ and standard deviation = 0.26$^{+0.06}_{-0.05}$ dex. We find no evidence for strong variation in this EW distribution with UV luminosity. The typical [OIII]+H$beta$ EW at z$simeq$7 implied by our sample is considerably larger than that in massive star forming galaxies at z$simeq$2, consistent with a shift toward larger average sSFR (4.4 Gyr$^{-1}$) and lower metallicities (0.16 Z$_odot$). We also find evidence for the emergence of a population with yet more extreme nebular emission ([OIII]+H$beta$ EW$>$1200 $mathrm{mathring{A}}$) that is rarely seen at lower redshifts. These objects have extremely large sSFR ($>$30 Gyr$^{-1}$), as would be expected for systems undergoing a burst or upturn in star formation. While this may be a short-lived phase, our results suggest that 20% of the z$simeq$7 population has such extreme nebular emission, implying that galaxies likely undergo intense star formation episodes regularly at z$>$6. We argue that this population may be among the most effective ionizing agents in the reionization era, both in terms of photon production efficiency and escape fraction. We furthermore suggest that galaxies passing through this large sSFR phase are likely to be very efficient in forming bound star clusters.
Ly$alpha$ emission from galaxies can be utilized to characterize the ionization state in the intergalactic medium (IGM). We report our search for Ly$alpha$ emission at $z>7$ using a comprehensive Keck/MOSFIRE near-infrared spectroscopic dataset, as part of the Texas Spectroscopic Search for Ly$alpha$ Emission at the End of Reionization Survey. We analyze data from 10 nights of MOSFIRE observations which together target 72 high-$z$ candidate galaxies in the GOODS-N field, all with deep exposure times of 4.5-19 hr. Utilizing an improved automated emission-line search, we report 10 Ly$alpha$ emission lines detected ($>$4$sigma$) at $z>7$, significantly increasing the spectroscopically confirmed sample. Our sample includes large equivalent-width (EW) Ly$alpha$ emitters ($>$50r{A}), and additional tentative Ly$alpha$ emission lines detected at 3 - 4$sigma$ from five additional galaxies. We constrain the Ly$alpha$ EW distribution at $zsim7.6$, finding a significant drop from $zlesssim6$, suggesting an increasing fraction of neutral hydrogen (HI) in the IGM in this epoch. We estimate the Ly$alpha$ transmission through the IGM ($=$EW$_{zsimtext{7.6}}$/EW$_{zsimtext{2-6}}$), and infer an IGM HI fraction ($X_{text{HI}}$) of $49^{+19}_{-19}%$ at $zsim7.6$, which is lower in modest tension ($>$1$sigma$) with recent measurements at $z sim$ 7.6. The spatial distribution of the detected Ly$alpha$ emitters implies the presence of a potential highly ionized region at $zsim7.55$ which hosts four Ly$alpha$ emitters within a $sim$ 40 cMpc spatial separation. The prominence of this ionized region in our dataset could explain our lower inferred value of $X_{text{HI}}$, though our analysis is also sensitive to the chosen reference Ly$alpha$ EW distribution values and reionization models.
Reionization-era galaxies tend to exhibit weak Ly$alpha$ emission, likely reflecting attenuation from an increasingly neutral IGM. Recent observations have begun to reveal exceptions to this picture, with strong Ly$alpha$ emission now known in four of the most massive z=7$-$9 galaxies in the CANDELS fields, all of which also exhibit intense [OIII]$+$H$beta$ emission (EW$>$800 $mathrm{mathring{A}}$). To better understand why Ly$alpha$ is anonymously strong in a subset of massive z$simeq$7$-$9 galaxies, we have initiated an MMT/Binospec survey targeting a larger sample (N=22) of similarly luminous ($simeq$1$-$6 L$^{ast}_{mathrm{UV}}$) z$simeq$7 galaxies selected over very wide-area fields ($sim$3 deg$^2$). We confidently ($>$7$sigma$) detect Ly$alpha$ in 78% (7/9) of galaxies with strong [OIII]$+$H$beta$ emission (EW$>$800 $mathrm{mathring{A}}$) as opposed to only 8% (1/12) of galaxies with more moderate (EW=200$-$800 $mathrm{mathring{A}}$) [OIII]$+$H$beta$. We argue that the higher Ly$alpha$ EWs of the strong [OIII]$+$H$beta$ population likely reflect enhanced ionizing photon production efficiency owing to their large sSFRs ($gtrsim$30 Gyr$^{-1}$). We also find evidence that Ly$alpha$ transmission from massive galaxies declines less rapidly over $6<z<7$ than in low-mass lensed systems. In particular, our data suggest no strong evolution in Ly$alpha$ transmission, consistent with a picture wherein massive z$simeq$7 galaxies often reside in large ionized regions. We detect three closely-separated ($R$ = 1.7 physical Mpc) z$simeq$7 Ly$alpha$ emitters in our sample, conceivably tracing a large ionized structure that is consistent with this picture. We detect tentative evidence for an overdensity in this region, implying a large ionizing photon budget in the surrounding volume.
The distribution of Ly$alpha$ emission is an presently accessible method for studying the state of the intergalactic medium (IGM) into the reionization era. We carried out deep spectroscopic observations in order to search for Ly$alpha$ emission from galaxies with photometric redshifts $z$ = 5.5 - 8.3 selected from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Utilizing data from the Keck/DEIMOS spectrograph, we explore a wavelength coverage of Ly$alpha$ emission at $z$ ~ 5 - 7 with four nights of spectroscopic observations for 118 galaxies, detecting five emission lines with ~ 5$sigma$ significance: three in the GOODS-N and two in the GOODS-S field. We constrain the equivalent width (EW) distribution of Ly$alpha$ emission by comparing the number of detected objects with the expected number constructed from detailed simulations of mock emission lines that account for the observational conditions (e.g., exposure time, wavelength coverage, and sky emission) and galaxy photometric redshift probability distribution functions. The Ly$alpha$ EW distribution is well described by an exponential form, $text{dN/dEW}propto text{exp(-EW/}W_0)$, characterized by the $e$-folding scale ($W_0$) of ~ 60 - 100$AA$ at 0.3 < $z$ < 6. By contrast, our measure of the Ly$alpha$ EW distribution at 6.0 < $z$ < 7.0 rejects a Ly$alpha$ EW distribution with $W_0$ > 36.4$AA$ (125.3$AA$) at 1$sigma$ (2$sigma$) significance. This provides additional evidence that the EW distribution of Ly$alpha$ declines at $z$ > 6, suggesting an increasing fraction of neutral hydrogen in the IGM at that epoch.
Recent surveys have identified a seemingly ubiquitous population of galaxies with elevated [OIII]/H$beta$ emission line ratios at $z > 1$, though the nature of this phenomenon continues to be debated. The [OIII]/H$beta$ line ratio is of interest because it is a main component of the standard diagnostic tools used to differentiate between active galactic nuclei (AGN) and star-forming galaxies, as well as the gas-phase metallicity indicators $O_{23}$ and $R_{23}$. Here, we investigate the primary driver of increased [OIII]/H$beta$ ratios by median-stacking rest-frame optical spectra for a sample of star-forming galaxies in the 3D-HST survey in the redshift range $zsim1.4-2.2$. Using $N = 4220$ star-forming galaxies, we stack the data in bins of mass and specific star formation rates (sSFR) respectively. After accounting for stellar Balmer absorption, we measure [OIII]$lambda5007$AA/H$beta$ down to $mathrm{M} sim 10^{9.2} mathrm{M_odot}$ and sSFR $sim 10^{-9.6} mathrm{yr}^{-1}$, more than an order of magnitude lower than previous work at similar redshifts. We find an offset of $0.59pm0.05$ dex between the median ratios at $zsim2$ and $zsim0$ at fixed stellar mass, in agreement with existing studies. However, with respect to sSFR, the $z sim 2$ stacks all lie within 1$sigma$ of the median SDSS ratios, with an average offset of only $-0.06pm 0.05$. We find that the excitation properties of galaxies are tightly correlated with their sSFR at both $zsim2$ and $zsim0$, with a relation that appears to be roughly constant over the last 10 Gyr of cosmic time.
Using deep multi-wavelength photometry of galaxies from ZFOURGE, we group galaxies at $2.5<z<4.0$ by the shape of their spectral energy distributions (SEDs). We identify a population of galaxies with excess emission in the $K_s$-band, which corresponds to [OIII]+H$beta$ emission at $2.95<z<3.65$. This population includes 78% of the bluest galaxies with UV slopes steeper than $beta = -2$. We de-redshift and scale this photometry to build two composite SEDs, enabling us to measure equivalent widths of these Extreme [OIII]+H$beta$ Emission Line Galaxies (EELGs) at $zsim3.5$. We identify 60 galaxies that comprise a composite SED with [OIII]+H$beta$ rest-frame equivalent width of $803pm228$AA and another 218 galaxies in a composite SED with equivalent width of $230pm90$AA. These EELGs are analogous to the `green peas found in the SDSS, and are thought to be undergoing their first burst of star formation due to their blue colors ($beta < -1.6$), young ages ($log(rm{age}/yr)sim7.2$), and low dust attenuation values. Their strong nebular emission lines and compact sizes (typically $sim1.4$ kpc) are consistent with the properties of the star-forming galaxies possibly responsible for reionizing the universe at $z>6$. Many of the EELGs also exhibit Lyman-$alpha$ emission. Additionally, we find that many of these sources are clustered in an overdensity in the Chandra Deep Field South, with five spectroscopically confirmed members at $z=3.474 pm 0.004$. The spatial distribution and photometric redshifts of the ZFOURGE population further confirm the overdensity highlighted by the EELGs.