Search

We investigate the `intrinsic' H$\alpha$ EW distributions of $z \sim 0.4 - 2.2$ narrowband-selected H$\alpha$ samples from HiZELS and DAWN using a forward modeling approach. We find an EW - stellar mass anti-correlation with steepening slopes $-0.18\pm0.03$ to $-0.24^{+0.06}_{-0.08}$ at $z \sim 0.4$ and $z\sim 2.2$, respectively. Typical EW increases as $(1+z)^{1.78^{+0.22}_{-0.23}}$ for a $10^{10}$ M$_\odot$ emitter from $15^{+2.4}_{-2.3}$\r{A} ($z \sim 0.4$) to $67.7^{+10.4}_{-10.0}$\r{A} ($z \sim 2.2$) and is steeper with decreasing stellar mass highlighting the high EW nature of low-mass high-$z$ systems. We model this redshift evolving anti-correlation, $W_0(M,z)$, and find it produces H$\alpha$ luminosity and SFR functions strongly consistent with observations validating the model and allowing us to use $W_0(M,z)$ to investigate the relative contribution of H$\alpha$ emitters towards cosmic SF. We find EW$_0 > 200$ \r{A} emitters contribute significantly to cosmic SF activity at $z \sim 1.5 - 2$ making up $\sim 40$% of total SF consistent with sSFR $> 10^{-8.5}$ yr$^{-1}$ ($\sim 45 - 55$%). Overall, this highlights the importance of high EW systems at high-$z$. Our $W_0(M,z)$ model also reproduces the cosmic sSFR evolution found in simulations and observations and show that tension between the two can simply arise from selection effects in observations. Lastly, we forecast Roman and Euclid grism surveys using $W_0(M,z)$ including observational efficiency and limiting resolution effects where we predict $\sim 24000$ and $\sim 30000$ $0.5 < z < 1.9$ H$\alpha$ emitters per deg$^{-2}$, respectively, down to $>5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ including $10^{7.2 - 8}$ M$_\odot$ galaxies at $z > 1$ with EW$_0 >1000$\r{A}. Both Roman and Euclid will enable us to observe with unprecedented detail some of the most bursty/high EW, low-mass star-forming galaxies near cosmic noon., Comment: 25 pages, 12 Figures, and 12 Tables. Submitted to MNRAS. Comments are welcomed. (Abridged Abstract)

Strong emission from doubly ionized oxygen is a beacon for some of the most intensely star forming galaxies known. JWST enables the search for this beacon in the early universe with unprecedented sensitivity. In this work, we extend the study of faint [OIII]$_{5008}$ selected galaxies by an order of magnitude in line luminosity. We use publicly available UNCOVER DR1 JWST/NIRCam and HST imaging data of the cluster lensing field, Abell 2744, to identify strong (rest-frame EW$>500$\AA) [OIII]$_{5008}$ emitters at $z\sim7$ based on excess F410M flux. We find $N=68$ $z\sim7$ [OIII] candidates, with a subset of $N=33$ that have deep HST coverage required to rule-out lower redshift interlopers (13.68 arcmin$^2$ with F814W $5\sigma$ depth $>28$ AB). Such strong emission lines can lead to very red colors that could be misinterpreted as evidence for old, massive stellar populations, but are shown to be due to emission lines where we have spectra. Using this deep HST sample and completeness simulations, which calculate the effective survey volume of the UNCOVER lensing field as a function of [OIII] luminosity, we derive a new [OIII] luminosity function (LF) extending to $41.09<\rm{log}_{10}(L/\rm{erg\,s}^{-1})<42.35$ which is an order of magnitude deeper than previous $z\sim6$ [OIII] LFs based on JWST slitless spectroscopy. This LF is well fit by a power law with a faint-end slope of $\alpha=-2.07^{+0.22}_{-0.23}$. There is little or no evolution between this LF and published [OIII] LFs at redshifts $3\lesssim z\lesssim7$, and no evidence of a turnover at faint luminosities. The sizes of these extreme [OIII] emitters are broadly similar to their low redshift counterparts, the green peas. The luminosity function of [OIII] emitters matches that of Lyman-$\alpha$ at the bright end, suggesting that many of them should be Lyman-$\alpha$ emitters., Comment: submitted to ApJ

We searched the Chandra and XMM archives for observations of 900 green pea galaxies to find AGN signatures. Green peas are low-mass galaxies with prominent emission lines, similar in size and star formation rate to high-redshift dwarf galaxies. Of the 29 observations found, 9 show X-ray detections with $S/N>3$. The 2-10 keV X-ray luminosity for these 9 sources exceeds $10^{40}~\mathrm{erg~s}^{-1}$, with 2 sources exceeding $10^{41}~\mathrm{erg~s}^{-1}$, suggesting the presence of intermediate-mass black holes (IMBH) or low-luminosity AGN (LLAGN) with BH masses between $100-10^6M_\mathrm{\odot}$. All X-ray detected sources (plus 6 additional sources) show He~II$\lambda4686$ emission and a broad component of the H$\alpha$ emission line, indicating winds. The line widths of the broad H$\alpha$ and He II$\lambda4686$ emitting gas clouds are weakly correlated ($R^{2}=0.15$), suggesting He II$\lambda4686$ emission is inconsistent with winds from super-Eddington accretors. However, the ratio of X-ray luminosity to star formation rate shows an anti-correlation with metallicity in 5 out of 9 X-ray detected sources, implying ultraluminous X-ray sources are key contributors to the observed X-ray luminosity. This could be due to super-Eddington accretors or IMBH. The X-ray emission is much higher than that produced by Wolf-Rayet stars and supernovae-driven winds. Thus, the X-ray luminosity in these 9 sources can only be explained by black holes with masses over $100~M_\mathrm{\odot}$. Our findings suggest the presence of LLAGN in these galaxies, with broad H$\alpha$ line widths implying BH masses of $10^4-10^6M_\mathrm{\odot}$. Given Green Peas' role as significant Lyman Continuum leakers, LLAGN in these galaxies could have contributed significantly to cosmic reionization., Comment: Submitted to ApJ. 17 pages, 11 figures and 3 tables. Comments welcome

This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [O III] 5007 Angstrom and H beta; the other two targets have have JWST coverage of PAH 3.3 micron and complementary ALMA data covering the [C II] 158 micron emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has technical goals to establish best practices for the Integral Field Units (IFU) within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES's observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks., Comment: 28 pages, including 8 figures and 3 tables. Submitted to The Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program, https://iopscience.iop.org/collections/apj-231204-01-templates

We present the first results of the Hubble Deep Hydrogen Alpha (HDH$\alpha$) project, which analyzes the space-borne deep H$\alpha$ narrowband imaging data in the GOODS-S region. The HDH$\alpha$ data comprises 72 orbits' images taken with the HST ACS/WFC F658N filter. The exposure time varies across a total area of $\sim$76.1 $\rm{arcmin}^2$, adding up to a total exposure time of 195.7 ks, among which 68.8 ks are spent in the deepest region. These images are aligned, reprojected, and combined to have the same pixel grid as the Hubble Legacy Fields (HLF). The scientific goals of the HDH$\alpha$ include establishing a sample of emission-line galaxies (ELGs) including [O III] emitters at $z\sim$ 0.3, [O II] emitters at $z\sim$ 0.8, and Lyman-$\alpha$ emitters (LAEs) at $z \sim 4.4$, studying the line morphology of ELGs with high resolution imaging data, and statistically analyzing the line luminosity functions and line equivalent-width distributions of ELGs selected with HST. Furthermore, the HDH$\alpha$ project enhances the legacy value of the GOODS-S field by contributing the first HST-based narrowband image to the existing data sets, which includes the HST broadband data and other ancillary data from X-ray to radio taken by other facilities. In this paper, we describe the data reduction process of the HDH$\alpha$, select ELGs based on HST's F658N and broadband data, validate the redshifts of the selected candidates by cross matching with the public spectroscopic catalogs in the GOODS-S, and present a final catalog of the confirmed [O III] emitters at $z\sim$ 0.3, [O II] emitters at $z\sim$ 0.8, and LAEs at $z \sim 4.4$., Comment: 27 pages, 14 figures, 9 tables, accepted by ApJS

Around 70 $z>6.5$ luminous quasars have been discovered, strongly biased toward the bright end, thus not providing a comprehensive view on quasar abundance beyond cosmic dawn. We present the predicted results of Roman/Rubin high-redshift quasar survey, yielding 3 times more, $2-4$ magnitudes deeper quasar samples, probing high-redshift quasars across broad range of luminosities, especially faint quasars at $L_\mathrm{bol}\sim 10^{10}\;L_{\odot}$ or $M_\mathrm{1450} \sim-22$ that are currently poorly explored. We include high-$z$ quasars, galactic dwarfs and low-$z$ compact galaxies with similar colors as quasar candidates. We create mock catalogs based on population models to evaluate selection completeness and efficiency. We utilize classical color dropout method in $z$ and $Y$ bands to select primary quasar candidates, followed up with Bayesian selection method to identify quasars. We show that overall selection completeness $> 80\%$ and efficiency $\sim 10\%$ at $6.56.5$, 20 at $z > 7.5$ and 2 at $z > 8.5$. The quasar yields depend sensitively on the assumed quasar luminosity shape and redshift evolution. Brown dwarf rejection through proper motion up to 50$\%$ can be made for stars brighter than 25 mag, low-$z$ galaxies dominate at fainter magnitude. Our results show that Roman/Rubin are able to discover a statistical sample of the earliest and faintest quasars in the Universe. The new valuable datasets worth follow up studies with James Webb Space Telescope and Extremely Large Telescopes, to determine quasar luminosity function faint end slope and constraint the supermassive black holes growth in the early Universe., Comment: Accepted for publication in ApJ

Extreme emission line galaxies (EELGs) are considered the primary contributor to cosmic reionization and are valuable laboratories to study the astrophysics of massive stars. It is strongly expected that Roman's High Latitude Wide Area Survey (HLWAS) will find many strongly gravitationally lensed [O III] emitters at Cosmic Noon (1 < z < 2.8). Roman imaging and grism spectroscopy alone will simultaneously confirm these strong lens systems and probe their interstellar medium (ISM) and stellar properties on small scales ($\lesssim$ 100 pc). Moreover, these observations will synergize with ground-based and space-based follow-up observations of the discovered lensed [O III] emitters in multi-wavelength analyses of their properties (e.g., massive stars and possible escape of ionizing radiation), spatially resolved on the scales of individual star cluster complexes. Only Roman can uniquely sample a large number of lensed [O III] emitters to study the small scale (~ 100 pc) ISM and stellar properties of these extreme emission line galaxies, detailing the key physics of massive stars and the ISM that govern cosmic reionization., Comment: Submitted in response to the call for Roman Telescope CCS white papers

Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies., Comment: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

The slitless grism on the Nancy Grace Roman Space Telescope will enable deep near-infrared spectroscopy over a wide field of view. We demonstrate Roman's capability to detect Ly$\alpha$ galaxies at $z>7$ using a multi-position-angle (PA) observational strategy. We simulate Roman grism data using a realistic foreground scene from the COSMOS field. We also input fake Ly$\alpha$ galaxies spanning redshift z=7.5-10.5 and a line-flux range of interest. We show how a novel data cube search technique -- CUBGRISM -- originally developed for GALEX can be applied to Roman grism data to produce a Ly$\alpha$ flux-limited sample without the need for continuum detections. We investigate the impact of altering the number of independent PAs and exposure time. A deep Roman grism survey with 25 PAs and a total exposure time of $70$hrs can achieve Ly$\alpha$ line depths comparable to the deepest $z=7$ narrow-band surveys ($L_{{\rm{Ly}}\alpha}\gtrsim10^{43}$erg s$^{-1}$). Assuming a null result, where the opacity of the intergalactic medium (IGM) remains unchanged from $z\sim7$, this level of sensitivity will detect $\sim400$ deg$^{-2}$ Ly$\alpha$ emitters from $z=7.25-8.75$. A decline from this expected number density is the signature of an increasing neutral hydrogen fraction and the onset of reionization. Our simulations indicate that a deep Roman grism survey has the ability to measure the timing and magnitude of this decline, allowing us to infer the ionization state of the IGM and helping us to distinguish between models of reionization., Comment: 11 pages, 7 figures, Submitted to ApJ

Although double-peaked narrow emission-line galaxies have been studied extensively in the past years, only a few are reported with the green pea galaxies (GPs). Here we present our discovery of five GPs with double-peaked narrow [OIII] emission lines, referred to as DPGPs, selected from the LAMOST and SDSS spectroscopic surveys. We find that these five DPGPs have blueshifted narrow components more prominent than the redshifted components, with velocity offsets of [OIII]$\lambda$5007 lines ranging from 306 to 518 $\rm km\, s^{-1}$ and full widths at half maximums (FWHMs) of individual components ranging from 263 to 441 $\rm km\, s^{-1}$. By analyzing the spectra and the spectral energy distributions (SEDs), we find that they have larger metallicities and stellar masses compared with other GPs. The H$\alpha$ line width, emission-line diagnostic, mid-infrared color, radio emission, and SED fitting provide evidence of the AGN activities in these DPGPs. They have the same spectral properties of Type 2 quasars. Furthermore, we discuss the possible nature of the double-peaked narrow emission-line profiles of these DPGPs and find that they are more likely to be dual AGN. These DPGP galaxies are ideal laboratories for exploring the growth mode of AGN in the extremely luminous emission-line galaxies, the co-evolution between AGN and host galaxies, and the evolution of high-redshift galaxies in the early Universe., Comment: 13 pages, 8 figures, 3 tables; Accepted for publication in MNRAS

We use large simulations of Lyman-Alpha Emitters with different fractions of ionized intergalactic medium to quantify the clustering of Ly$\alpha$ emitters as measured by the Void Probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace Roman Space Telescope. Since Roman will be able to carry out blind spectroscopic surveys of Ly$\alpha$ emitters continuously between $7 3-4\sigma$ at several redshifts between $75-8\sigma$) across the epoch of reionization, and would yield a detailed history of the reionization of the IGM and its effect on Lyman-$\alpha$ Emitter clustering., Comment: 19 total pages, 7 total figures, appendix with 3 additional figures

We probe what constraints for the global ionized hydrogen fraction the Void Probability Function (VPF) clustering can give for the Lyman-Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyman-alpha emission, and measuring the drop in the luminosity function of Lyman-$\alpha$ emitters (LAEs) has been used to constrain the ionization fraction in narrowband surveys. However, the clustering of LAEs is independent from the luminosity function's inherent evolution, and can offer additional constraints for reionization under different models. The VPF measures how likely a given circle is to be empty. It is a volume-averaged clustering statistic that traces the behavior of higher order correlations, and its simplicity offers helpful frameworks for planning surveys. Using the \citet{Jensen2014} simulations of LAEs within various amount of ionized intergalactic medium, we predict the behavior of the VPF in one (301x150.5x30 Mpc$^3$), four (5.44$\times 10^6$ Mpc$^3$), or eight (1.1$\times 10^7$ Mpc$^3$) fields of LAGER imaging. We examine the VPF at 5 and 13 arcminutes, corresponding to the minimum scale implied by the LAE density and the separation of the 2D VPF from random, and the maximum scale from the 8-field 15.5 deg$^2$ LAGER area. We find that even a single DECam field of LAGER (2-3 deg$^2$) could discriminate between mostly neutral vs. ionized. Additionally, we find four fields allows the distinction between 30, 50, and 95 percent ionized; and that eight fields could even distinguish between 30, 50, 73, and 95 percent ionized., Comment: 14 pages, 5 figures

The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizations of 2-deg^2 lightcones, containing a total of >25 million simulated galaxies with -16 < MUV < -25 spanning z ~ 0 to 10. This dataset enables a new set of experiments with the impacts of survey size on the derived galaxy formation and cosmological constraints. The intrinsic and observable galaxy properties are predicted using a well-established, physics-based semi-analytic modelling approach. We provide forecasts for number density, cosmic SFR, field-to-field variance, and angular two-point correlation functions, and demonstrate how the future wide-field surveys will be able to improve these measurements relative to current generation surveys. We also present a comparison between these lightcones and others that have been constructed with empirical models. The mock lightcones are designed to facilitate the exploration of multi-instrument synergies and connecting with current generation instruments and legacy surveys. In addition to Roman, we also provide photometry for a number of other instruments on upcoming facilities, including Euclid and Rubin, as well as the instruments that are part of many legacy surveys. Full object catalogues and data tables for the results presented in this work are made available through a web-based, interactive portal https://www.simonsfoundation.org/semi-analytic-forecasts., Comment: 25 pages, 16 figures. Accepted for publication in MNRAS

The Early Release Observations (EROs) of JWST beautifully demonstrate the promise of JWST in characterizing the universe at cosmic dawn. We analyze the ERO spectra of three $z \sim 8$ galaxies to determine their metallicities, gas temperatures and ionization. These galaxies offer the first opportunity to understand the physical properties of epoch-of-reionization galaxies through detailed rest-optical emission line spectroscopy. We show that these objects have metal abundances $12+\log[O/H] \approx 6.9 - 8.2$, based on both the $T_e$ method and on a recent calibration of the $R_{23}$ metallicity indicator. Since the spectra are some of the earliest science data from JWST, we compare several line ratios with values expected from robust physics, to validate our measurement procedures. We compare the abundances and emission line ratios to a nearby sample of Green Pea galaxies -- a population of nearby emission line galaxies whose UV properties resemble epoch-of-reionization galaxies, and which often have large Lyman continuum escape fractions. The JWST data show striking further similarities between these high redshift galaxies and nearby Green Peas. The $z\sim 8$ galaxies span the metallicity range covered by Green Peas. They also show the compact morphology that is typical of emission line dominated galaxies at all redshifts. Based on these similarities with Green Peas, it is likely that these are the first rest-optical spectra of galaxies that are actively driving cosmological reionization, Comment: 13 pages, 6 figures; Submitted to ApJ Letters

We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sensitive probes of cosmic reionization and narrow-band imaging is a robust and effective method for selecting a large number of LAEs. In this work, we present results from the spectroscopic follow-up of LAE candidates in two LAGER fields, COSMOS and WIDE-12, using observations from Keck/LRIS. We report the successful detection of Ly$\alpha$ emission in 15 candidates (11 in COSMOS and 4 in WIDE-12 fields). Three of these in COSMOS have matching confirmations from a previous LAGER spectroscopic follow-up and are part of the overdense region, LAGER-$z7$OD1. Additionally, two candidates that were not detected in the LRIS observations have prior spectroscopic confirmations from Magellan. Including these, we obtain a spectroscopic confirmation success rate of $\sim$$80$% for LAGER LAE candidates. Apart from Ly$\alpha$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$\sigma$ upper limit for the ratio of NV/Ly$\alpha$, $f_{NV}/f_{Ly\alpha} \lesssim 0.27$, which implies that ionizing emission from these sources is mostly dominated by star formation. Including confirmations from this work, a total of 33 LAE sources from LAGER are now spectroscopically confirmed. LAGER has more than doubled the sample of spectroscopically confirmed LAE sources at $z\sim7$., Comment: 8 pages, 4 figures; Submitted to ApJ

Dark matter halos that reach the HI-cooling mass without prior star formation or external metal pollution represent potential sites for the formation of small - extremely faint - Population III galaxies at high redshifts. Gravitational lensing may in rare cases boost their fluxes to detectable levels, but to find even a small number of such objects in randomly selected regions of the sky requires very large areas to be surveyed. Because of this, a small, wide-field telescope can in principle offer better detection prospects than a large telescope with a smaller field of view. Here, we derive the minimum comoving number density required to allow gravitational lensing to lift such objects at redshift $z=5-16$ above the detection thresholds of blind surveys carried out with the James Webb space telescope (JWST), the Roman space telescope (RST) and Euclid. We find that the prospects for photometric detections of Pop III galaxies is promising, and that they are better for RST than for JWST and Euclid. However, the Pop III galaxies favoured by current simulations have number densities too low to allow spectroscopic detections based on the strength of the HeII1640 emission line in any of the considered surveys unless very high star formation efficiencies ($\epsilon > 0.1$) are envoked. We argue that targeting individual cluster lenses instead of the wide field surveys considered in this paper results in better spectroscopic detection prospects, while for photometric detection, the wide field surveys perform considerably better., Comment: 15 pages, 4 figures

We explore the presence of active galactic nuclei (AGN)/black holes (BH) in Green Pea galaxies (GPs), motivated by the presence of high ionization emission lines such as HeII and [NeIII] in their optical spectra. In order to identify AGN candidates, we used mid-infrared (MIR) photometric observations from the all-sky Wide-field Infrared Survey Explorer (WISE) mission for a sample of 1004 GPs. Considering only $>5\sigma$ detections with no contamination from neighboring sources in AllWISE, we select 31 GPs out of 134 as candidate AGN based on a stringent 3-band WISE color diagnostic. Using multi-epoch photometry in W1 and W2 bands based on time-resolved unWISE coadd images, we find two sources exhibiting variability in both the WISE bands among 112 GPs with W1$\leqslant16$ mag and no contamination from neighboring sources in unWISE. These two variable sources were selected as AGN by the WISE 3-band color diagnostic as well. Compared to variable AGN fractions observed among low-mass galaxy samples in previous studies, we find a higher fraction ($\sim1.8\%$) of MIR variable sources among GPs, which demonstrates the uniqueness and importance of studying these extreme objects. Through this work, we demonstrate that MIR diagnostics are promising tools to select AGN that may be missed by other selection techniques (including optical emission-line ratios and X-ray emission) in star-formation dominated, low-mass, low-metallicity galaxies., Comment: 10 pages, 6 figures; Resubmitted to ApJ

We present a new measurement of the Ly$\alpha$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest $z\sim 7$ survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify $N=95$ new $z=6.9$ Ly$\alpha$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$\alpha$ luminosity functions. The best-fit Schechter luminosity function parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Ly$\alpha$ luminosity density declines at the same rate as the UV continuum LF from $z=5.7$ to $z=6.9$. This is consistent with an intergalactic medium that was fully ionized as early as redshift $z\sim 7$, or with a volume-averaged neutral hydrogen fraction of $x_{HI} < 0.33$ at $1\sigma$., Comment: 19 pages, 11 figures; accepted by ApJ

We study the dependence of Ly$\alpha$ escape from galaxies on UV continuum size and luminosity using a sample of 40 Green Pea (GP) galaxies, which are the best local analogs of high-redshift Ly$\alpha$-emitters (LAEs). We use the Cosmic Origins Spectrograph near-ultraviolet images from the \textit{Hubble Space Telescope} to measure the UV size and luminosity with $0.047''$ spatial resolution. Like most galaxies the GPs show a log-normal size distribution. They also show a positive correlation between size and UV-continuum luminosity. The slope of the size-continuum luminosity relation for GPs is consistent with those of continuum-selected star-forming galaxies at low and high redshifts. A distinctive feature of GPs is a very compact typical radius of 0.33 kpc with a population spread (1$\sigma$) of 0.19 kpc. The peak of the size distribution and the intercept of the size-luminosity relation of GPs are noticeably smaller than those of continuum-selected star-forming galaxies at similar redshifts. There are statistically significant anti-correlations found between the circularized half-light radius ($r_{\rm cir,50}$), the Ly$\alpha$ equivalent width (EW(Ly$\alpha$)), and the Ly$\alpha$ escape faction ($f^{Ly\alpha}_{esc}$), suggesting that small UV-continuum radii are crucial for Ly$\alpha$ emission. GPs and high-redshift LAEs have similar sizes, once spatial resolution effects are properly considered. Our results show that a compact/small size is crucial for escape of Ly$\alpha$ photons, and that Ly$\alpha$-emitters show constant characteristic size independent of their redshift., Comment: 10 pages, 4 figures, 2 tables, accepted by ApJ; doi:10.3847/1538-4357/abf833

The H$\alpha$ equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic star-formation activity. Previous assessments show that EW strongly anti-correlates with stellar mass as $M^{-0.25}$ similar to the sSFR -- stellar mass relation. However, such a correlation may be driven/formed by selection effects. In this study, we investigate how H$\alpha$ EWs correlate with galaxy properties and how selection biases could alter such correlations using a $z = 0.47$ narrowband-selected sample of 1572 H$\alpha$ emitters from the Ly$\alpha$ Galaxies in the Epoch of Reionization (LAGER) survey. The sample covers 3 deg$^2$ of COSMOS and $1.1\times10^5$ cMpc$^3$. We assume an intrinsic EW distribution to form mock samples of H$\alpha$ emitters (HAEs) and propagate the selection criteria to match observations, giving us control on how selection biases can affect the underlying results. We find EW intrinsically correlates with stellar mass as $W_0 \propto M^{-0.16\pm0.03}$ and decreases by a factor of $\sim 3$ from $10^{7}$ to $10^{10}$ M$_\odot$. We find low-mass HAEs to be $\sim 320$ times more likely to have rest-frame EW$ > 200$\AA compared to high-mass HAEs. Combining the intrinsic EW -- stellar mass correlation with an observed SMF correctly reproduces the observed H$\alpha$ LF, while not correcting for selection effects underestimates the number of bright HAEs. This suggests that the intrinsic EW -- stellar mass correlation is physically significant and reproduces three statistical distributions of galaxy populations (LF, SMF, EW distribution). At lower masses, we find there are more high-EW outliers compared to high masses, even after taking into account selection effects. Our results suggest that high sSFR outliers indicative of bursty SF activity are intrinsically more prevalent in low-mass HAEs and not a byproduct of selection effects., Comment: 21 Pages, 13 Figures, 2 Tables, accepted to MNRAS

Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million years old and could be experiencing rapid evolution of the neutral hydrogen fraction in the intergalactic medium. The protocluster is identified by an overdensity of 6 times the average galaxy density, and with 21 narrowband selected Lyman-$\alpha$ galaxies, among which 16 have been spectroscopically confirmed. At redshifts similar to or above this record, smaller protogroups with fewer members have been reported. LAGER-z7OD1 shows an elongated shape and consists of two subprotoclusters, which would have merged into one massive cluster with a present-day mass of $3.7 \times 10^{15}$ solar masses. The total volume of the ionized bubbles generated by its member galaxies is found to be comparable to the volume of the protocluster itself, indicating that we are witnessing the merging of the individual bubbles and that the intergalactic medium within the protocluster is almost fully ionized. LAGER-z7OD1 thus provides a unique natural laboratory to investigate the reionization process., Comment: 9 pages, 5 figures, and 1 table. Accepted by Nature Astronomy

We calculate the void probability function (VPF) in simulations of Lyman-$\alpha$ emitters (LAEs) across a wide redshift range ($z=3.1,\ 4.5,\ 5.7,\ 6.6$). The VPF measures the zero-point correlation function (i.e. places devoid of galaxies) and naturally connects to higher order correlation functions while being computationally simple to calculate. We explore the Poissonian and systematic errors on the VPF, specify its accuracy as a function of average source density and the volume probed, and provide the appropriate size scales to measure the VPF. At small radii the accuracy of the VPF is limited by galaxy density, while at large radii the VPF is limited by the number of independent volumes probed. We also offer guidelines for understanding and quantifying the error in the VPF. We approximate the error in the VPF by using independent sub-volumes of the catalogs, after finding that jackknife statistics underestimate the uncertainty. We use the VPF to probe the strength of higher order correlation functions by measuring and examining the hierarchical scaling between the correlation functions using count-in-cells. The negative binomial model (NBM) has been shown to best describe the scaling between the two point correlation function and VPF for low-redshift galaxy observations. We further test the fit of the NBM by directly deriving the volume averaged two-point correlation function from the VPF and vice versa. We find the NBM best describes the $z=3.1, 4.5, 5.7$ simulated LAEs, with a 1$\sigma$ deviation from the model in the $z=6.6$ catalog. This suggests that LAEs show higher order clustering terms similar to those of normal low redshift galaxies., Comment: 22 pages, 6 figures

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 ($>$50\r{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 $z\sim7.6$, finding a significant drop from $z\lesssim6$, suggesting an increasing fraction of neutral hydrogen (HI) in the IGM in this epoch. We estimate the Ly$\alpha$ transmission through the IGM ($=$EW$_{z\sim\text{7.6}}$/EW$_{z\sim\text{2-6}}$), and infer an IGM HI fraction ($X_{\text{HI}}$) of $49^{+19}_{-19}\%$ at $z\sim7.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 $z\sim7.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., Comment: Published 2020 November 27 in ApJ. Changed to match the published version (HI fraction & Figure 12 updated and other minor changes). 27 pages, 15 figures, 3 tables

We have studied ultraviolet images of 40 Green Pea galaxies and 15 local Lyman Break Galaxy Analogs to understand the relation between Ly$\alpha$ photon escape and central UV photometric properties. We measured star formation intensity (SFI, star formation rate per unit area) from the central 250 pc region ($S_{\rm 250pc}$) using COS/NUV images from the \textit{Hubble Space Telescope}. The measured $S_{\rm 250pc}$ of our sample Green Peas ranges from 2.3--46 $M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$, with a geometric mean of $15 M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$ and a standard deviation of 0.266 dex, forming a relatively narrow distribution. The Lyman Break Galaxy Analogs show a similarly narrow distribution of $S_{\rm 250pc}$ (0.271 dex), though with a larger mean of 28 $M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$. We show that while the Ly$\alpha$ equivalent width (EW(Ly$\alpha$)) and the Ly$\alpha$ escape fraction ($f^{Ly\alpha}_{esc}$) are not significantly correlated with the central SFI ($S_{\rm 250pc}$), both are positively correlated with the ratio of surface brightness to galaxy stellar mass ($S_{\rm 250pc}/M_{\rm star}$), with correlation coefficients ($p$-values) of 0.702 ($1\times 10^{-8}$) and 0.529 ($5\times 10^{-4}$) with EW(Ly$\alpha$) and $f^{Ly\alpha}_{esc}$, respectively. These correlations suggest a scenario where intense central star formation can drive a galactic wind in galaxies with relatively shallow gravitational potential wells, thus clearing channels for the escape of Ly$\alpha$ photons., Comment: 8 pages, 3 figures, 2 tables. Accepted for publication in ApJ

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$\alpha$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$\alpha$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$\alpha$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}\alpha}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$\alpha$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $\phi^*=10^{-2.85}$ Mpc$^{-3}$, $\alpha = -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $\phi^*=10^{-2.8}$ Mpc$^{-3}$, $\alpha=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $\phi^*=10^{-2.91}$ Mpc$^{-3}$, $\alpha = -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$\alpha$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $\rho_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $\rho_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $\rho_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$\alpha$ surveys., Comment: 16 pages, 8 figures, Accepted for publication in ApJ

We present our measurements of the H$\alpha$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$\alpha$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$\sigma$ depths of $8.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and comoving volumes of $(1-7)\times10^{5}$ Mpc$^3$ making our survey one of the deepest narrowband surveys. We measure the observed luminosity functions and find best-fits of $\phi^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$\alpha$, $\phi^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc$^{-3}$ and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s$^{-1}$ for [OIII], and $\phi^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $\alpha$ fixed to $-1.75$, $-1.6$, and $-1.3$, respectively. An excess of bright $> 10^{42}$ erg s$^{-1}$ [OIII] emitters is observed and may be due to AGN contamination. Dust corrections are applied assuming $A_{\rm{H}\alpha} = 1$ mag. We also design our own empirical rest-frame $g - r$ calibration using SDSS DR12 data, test it against our $z = 0.47$ H$\alpha$ emitters with $z$COSMOS $1$D spectra, and calibrate it for $(g - r)$ between $-0.8$ and $1.3$ mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as $\log_{10}\rho_{\rm{SFR}}/(\rm{M}_\odot\ \rm{yr}^{-1}\ \rm{Mpc}^{-3}) = -1.63\pm0.04$, $-1.07\pm0.06$, and $-0.90\pm0.10$ for H$\alpha$, [OIII], and [OII], respectively. We find our [OIII] and [OII] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multi-wavelength SFRDs, while the H$\alpha$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD., Comment: 21 pages, 12 figures, 6 Tables. Accepted by MNRAS for publication on 6 Jan 2020

We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$\alpha$ $\lambda6563$, [OIII]$\lambda\lambda4959,5007$, and [OII]$\lambda\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$\alpha$ emitters and for $0.8

Extragalactic studies have demonstrated there is a moderately tight ($\approx$0.3 dex) relationship between galaxy stellar mass ($M_{\star}$) and star formation rate (SFR) that holds for star-forming galaxies at $M_{\star} \sim 3 \times 10^8$-10$^{11}~M_{\odot}$, i.e., the "star formation main sequence." However, it has yet to be determined whether such a relationship extends to even lower mass galaxies, particularly at intermediate or higher redshifts. We present new results using observations for 714 narrowband H$\alpha$-selected galaxies with stellar masses between $10^6$ and $10^{10}~M_{\odot}$ (average of $10^{8.2}~M_{\odot}$) at $z \approx$ 0.07-0.5. These galaxies have sensitive UV to near-infrared photometric measurements and optical spectroscopy. The latter allows us to correct our H$\alpha$ SFRs for dust attenuation using Balmer decrements. Our study reveals: (1) for low-SFR galaxies, our H$\alpha$ SFRs systematically underpredict compared to FUV measurements, consistent with other studies; (2) at a given stellar mass ($\approx $10$ ^{8}~M_{\odot}$), log(specific SFR) evolves as $ A \log(1+z) $ with $ A = 5.26 \pm 0.75 $, and on average, specific SFR increases with decreasing stellar mass; (3) the SFR-$M_{\star}$ relation holds for galaxies down to $\sim$10$^6~M_{\odot}$ ($\sim$1.5 dex below previous studies), and over lookback times of up to 5 Gyr, follows a redshift-dependent relation of $\log{({\rm SFR})} \propto \alpha \log(M_{\star}/M_{\odot}) + \beta z$ with $\alpha = 0.60 \pm 0.01$ and $\beta = 1.86 \pm 0.07$; and (4) the observed dispersion in the SFR-$M_{\star}$ relation at low stellar masses is $\approx$0.3 dex. Accounting for survey selection effects using simulated galaxies, we estimate the true dispersion is $\approx$0.5 dex., Comment: Updated references to match published version

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1

We present Keck/MOSFIRE $H$-band spectroscopy targeting C$\,$III] $\lambda$1907,1909 in a $z=7.5056$ galaxy previously identified via Ly$\alpha$ emission. We detect strong line emission at $1.621\pm0.002\,\mu$m with a line flux of ($2.63\pm0.52$)$\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We tentatively identify this line as [C$\,$III] $\lambda$1907, but we are unable to detect C$\,$III] $\lambda$1909 owing to sky emission at the expected location. This gives a galaxy systemic redshift, $z_{sys}=7.5032\pm0.0003$, with a velocity offset to Ly$\alpha$ of $\Delta$v$_{Ly\alpha}$ = $88\pm27$ km s$^{-1}$. The ratio of combined C$\,$III]/Ly$\alpha$ is 0.30-0.45, one of the highest values measured for any $z>2$ galaxy. We do not detect Si$\,$III] $\lambda\lambda$1883, 1892, and place an upper limit on Si$\,$III]/C$\,$III] $<$ 0.35 ($2\sigma$). Comparing our results to photoionization models, the C$\,$III] equivalent width (W$_{\text{CIII]}} = 16.23\pm2.32\,$\AA), low Si$\,$III]/C$\,$III] ratio, and high implied [O$\,$III] equivalent width (from the $Spitzer$/IRAC [3.6]$-$[4.5]$\simeq$0.8 mag color) require sub-Solar metallicities ($Z\simeq0.1-0.2 Z_{\odot}$) and a high ionization parameter, log$\,$U $\gtrsim -1.5$. These results favor models that produce higher ionization, such as the BPASS models for the photospheres of high-mass stars, and that include both binary stellar populations and/or an IMF that extends to 300 M$_{\odot}$. The combined C$\,$III] equivalent width and [3.6]$-$[4.5] color are more consistent with ionization from young stars than AGN, however we cannot rule out ionization from a combination of an AGN and young stars. We make predictions for $James~Webb~Space~Telescope$ spectroscopy using these different models, which will ultimately test the nature of the ionizing radiation in this source., Comment: 19 pages, 12 figures, accepted for publication in ApJ

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts., Comment: 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function., Comment: 23 pages, 11 figures, 2 tables, accepted for publication in ApJ

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency., Comment: 12 pages, 12 figures, submitted to PASP

We spectroscopically confirmed two narrow-band selected redshift 7.0 Ly$\alpha$ galaxies and studied their restframe UV spectra. The Ly$\alpha$ and other UV nebular lines are very useful to confirm the galactic redshifts and diagnose the different mechanisms driving the ionizing emission. We observed two narrowband-selected $z$=7.0 Ly$\alpha$ candidates in the LAGER Chandra Deep Field South (CDFS) field with IMACS at Magellan telescope and confirmed they are Ly$\alpha$ emitters at $z$=6.924 and 6.931. In one galaxy, we also obtained deep NIR spectroscopy, which yields non-detections of the high-ionization UV nebular lines. We measured upper-limits of the ratios of CIV$\lambda$1548/Ly$\alpha$, HeII$\lambda$1640/Ly$\alpha$, OIII]$\lambda$1660/Ly$\alpha$, and CIII]$\lambda$1909/Ly$\alpha$ from the NIR spectra. These upper-limits imply that the ionizing emission in this galaxy is dominated by normal star formation instead of AGN., Comment: ApJ submitted

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade., Comment: 8 pages (including cover page and references), 3 figures. Science white paper submitted to Astro2020. arXiv admin note: substantial text overlap with arXiv:1802.01539

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization., Comment: 7 pages, 3 figures, Science white paper submitted to the Astro2020 Decadal Survey

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn., Comment: Science White paper submitted to Astro2020 Decadal Survey

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources., Comment: 6 pages, 3 figures

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production., Comment: 14 pages, 3 figures, 2 tables

Realizing the utility of Ly$\alpha$ emission to trace the evolution of the intergalactic medium (IGM) during the epoch of reionization requires deep spectroscopy across the boundary of optical and near-infrared (NIR) spectrographs at $z\sim7.2$ when Ly$\alpha$ emission is at $\sim$1$\mu$m. Our Texas Spectroscopic Search for Ly$\alpha$ Emission at the End of Reionization includes 18 nights of deep spectroscopic observations using the Keck DEIMOS (optical) and MOSFIRE (NIR) spectrographs. Within this dataset we observe Ly$\alpha$ emission from 183 photometric-redshift selected galaxies at $z =$ 5.5 - 8.3 from the Cosmic Assembly Near infrared Deep Extragalactic Legacy Survey (CANDELS). Our overlapping MOSFIRE observations, over 84 galaxies total, provide the deepest NIR spectroscopic data yet obtained for Ly$\alpha$ from galaxies $z > 7$, with $>16$ hr integration time for four observed galaxies. Here we analyze these four targets, and we report the discovery of a new $z = 7.60$ Ly$\alpha$ detection as well as provide an updated observation of the previously confirmed $z=7.51$ Ly$\alpha$ emission from Finkelstein et al. (2013) with a $\sim$3$\times$ longer exposure time. Our analysis of these Ly$\alpha$ emission line profiles reveal a significant asymmetric shape. The two detected Ly$\alpha$ emission lines from bright sources ($M_{\text{UV}}<-20.25$) could imply that these bright galaxies inhabit ionized bubbles in a partially neutral IGM, although deeper exposures may yet reveal Ly$\alpha$ emission in the fainter sources., Comment: 10 pages, 2 figures, 3 tables, Submitted to ApJ

The gas metallicity of galaxies is often estimated using strong emission lines such as the optical lines of [OIII] and [OII]. The most common measure is "R23", defined as ([OII]$\lambda$$\lambda$3726, 3729 + [OIII]$\lambda$$\lambda$4959,5007)/H$\beta$. Most calibrations for these strong-line metallicity indicators are for continuum selected galaxies. We report a new empirical calibration of R23 for extreme emission-line galaxies using a large sample of about 800 star-forming green pea galaxies with reliable T$_e$-based gas-phase metallicity measurements. This sample is assembled from Sloan Digital Sky Survey (SDSS) Data Release 13 with the equivalent width of the line [OIII]$\lambda$5007 $>$ 300 \AA\ or the equivalent width of the line H$\beta$ $>$ 100 \AA\ in the redshift range 0.011 $<$ z $<$ 0.411. For galaxies with strong emission lines and large ionization parameter (which manifests as log [OIII]$\lambda$$\lambda$4959,5007/[OII]$\lambda$$\lambda$3726,3729 $\geq$ 0.6), R23 monotonically increases with log(O/H) and the double-value degeneracy is broken. Our calibration provides metallicity estimates that are accurate to within $\sim$ 0.14 dex in this regime. Many previous R23 calibrations are found to have bias and large scatter for extreme emission-line galaxies. We give formulae and plots to directly convert R23 and [OIII]$\lambda$$\lambda$4959,5007/[OII]$\lambda$$\lambda$3726,3729 to log(O/H). Since green peas are best nearby analogs of high-redshift Lyman-$\alpha$ emitting galaxies, the new calibration offers a good way to estimate the metallicities of both extreme emission-line galaxies and high-redshift Lyman-$\alpha$ emitting galaxies. We also report on 15 galaxies with metallicities less than 1/12 solar, with the lowest metallicities being 12+log(O/H) = 7.25 and 7.26., Comment: 17 pages, 10 figures; Accepted for publication in ApJ

We explore the relation between the star formation rate surface density ($\Sigma$SFR) and the interstellar gas pressure for nearby compact starburst galaxies. The sample consists of 17 green peas and 19 Lyman break analogs. Green peas are nearby analogs of Ly$\alpha$ emitters at high redshift and Lyman break analogs are nearby analogs of Lyman break galaxies at high redshift. We measure the sizes for green peas using Hubble Space Telescope Cosmic Origins Spectrograph (COS) NUV images with a spatial resolution of $\sim$ 0.05$^{''}$. We estimate the gas thermal pressure in HII regions by $P = N_{total}Tk{_B} \simeq 2n_{e}Tk{_B}$. The electron density is derived using the [SII] doublet at 6716,6731 \AA, and the temperature is calculated from the [OIII] lines. The correlation is characterized by $\Sigma$ SFR = 2.40$\times$10$^{-3\,}$M$_{\odot\,}$yr$^{-1\,}$kpc$^{-2}$$\left(\frac{P/k_{B}}{10^{4}cm^{-3}K}\right)^{1.33}$. Green peas and Lyman break analogs have high $\Sigma$SFR up to 1.2 M$_{\odot\,}$yr$^{-1\,}$kpc$^{-2}$ and high thermal pressure in HII region up to P/k$_B$ $\sim$10$^{7.2}{\rm\, K\, cm}^{-3}$. These values are at the highest end of the range seen in nearby starburst galaxies. The high gas pressure and the correlation, are in agreement with those found in star-forming galaxies at z $\sim$ 2.5. These extreme pressures are shown to be responsible for driving galactic winds in nearby starbursts. These outflows may be a crucial in enabling Lyman-$\alpha$ and Lyman-continuum to escape., Comment: 16 pages, 6 figures, 1 table; Accepted for publication in ApJ

We derive direct measurement gas-phase metallicities of $7.4 < 12 + \log(O/H) < 8.4$ for 14 low-mass Emission Line Galaxies (ELGs) at $0.3 < z < 0.8$ identified in the Faint Infrared Grism Survey (FIGS). We use deep slitless G102 grism spectroscopy of the Hubble Ultra Deep Field (HUDF), dispersing light from all objects in the field at wavelengths between 0.85 and 1.15 microns. We run an automatic search routine on these spectra to robustly identify 71 emission line sources, using archival data from VLT/MUSE to measure additional lines and confirm redshifts. We identify 14 objects with $0.3 < z < 0.8$ with measurable O[III]$\lambda$4363 \AA\ emission lines in matching VLT/MUSE spectra. For these galaxies, we derive direct electron-temperature gas-phase metallicities with a range of $7.4 < 12 + \log(O/H) < 8.4$. With matching stellar masses in the range of $10^{7.9} M_{\odot} < M_{\star} < 10^{10.4} M_{\odot}$, we construct a mass-metallicity (MZ) relation and find that the relation is offset to lower metallicities compared to metallicities derived from alternative methods (e.g.,$R_{23}$, O3N2, N2O2) and continuum selected samples. Using star formation rates (SFR) derived from the $H\alpha$ emission line, we calculate our galaxies' position on the Fundamental Metallicity Relation (FMR), where we also find an offset toward lower metallicities. This demonstrates that this emission-line-selected sample probes objects of low stellar masses but even lower metallicities than many comparable surveys. We detect a trend suggesting galaxies with higher Specific Star Formation (SSFR) are more likely to have lower metallicity. This could be due to cold accretion of metal-poor gas that drives star formation, or could be because outflows of metal-rich stellar winds and SNe ejecta are more common in galaxies with higher SSFR., Comment: 14 pages, 11 figures, accepted in ApJ

We use both photometric and spectroscopic data from the {\it Hubble Space Telescope} to explore the relationships among 4000 \AA\ break (D4000) strength, colors, stellar masses, and morphology, in a sample of 352 galaxies with log$(M_{*}/M_{\odot}) > 9.44$ at 0.6 $\lesssim z \lesssim$ 1.2. We have identified authentically quiescent galaxies in the $UVJ$ diagram based on their D4000 strengths. This spectroscopic identification is in good agreement with their photometrically-derived specific star formation rates (sSFR). Morphologically, most (that is, 66 out of 68 galaxies, $\sim$ 97 \%) of these newly identified quiescent galaxies have a prominent bulge component. However, not all of the bulge-dominated galaxies are quenched. We found that bulge-dominated galaxies show positive correlations among the D4000 strength, stellar mass, and the S\'ersic index, while late-type disks do not show such strong positive correlations. Also, bulge-dominated galaxies are clearly separated into two main groups in the parameter space of sSFR vs. stellar mass and stellar surface density within the effective radius, $\Sigma_{\rm e}$, while late-type disks and irregulars only show high sSFR. This split is directly linked to the `blue cloud' and the `red sequence' populations, and correlates with the associated central compactness indicated by $\Sigma_{\rm e}$. While star-forming massive late-type disks and irregulars (with D4000 $<$ 1.5 and log$(M_{*}/M_{\odot}) \gtrsim 10.5$) span a stellar mass range comparable to bulge-dominated galaxies, most have systematically lower $\Sigma_{\rm e}$ $\lesssim$ $10^{9}M_{\odot}\rm{kpc^{-2}}$. This suggests that the presence of a bulge is a necessary but not sufficient requirement for quenching at intermediate redshifts., Comment: 23 pages, 11 figures, and 3 tables. Accepted for publication in ApJ

A population of so-called "blueberry" starbursting dwarf galaxies with extremely blue colors, low-metallicities, and enormous ionization ratios, has recently been found by Yang et al. (2017). Yet we still do not know their detailed properties, such as morphologies, AGN occupations, massive star contents, infrared emission, dust properties, etc. As a pilot study of the blueberries, we investigate the spectroscopic and photometric properties of two blueberry candidates, RGG B and RGG 5, for which Hubble Space Telescope high-resolution images are available. We find that RGG B and RGG 5 perhaps are likely to be two merging dwarf galaxy systems. RGG B may have a close merging companion; yet the current evidence still cannot exclude the possibility that RGG B is just disturbed by in-situ star formation through, e.g., outflows, rather than undergoing a merger. RGG 5 presents stellar shells in the outskirt which can be the powerful evidence of galaxy merging. We also find that, all of the blueberries, including RGG B and RGG 5, are located close to the theoretical maximum-starburst-line in the BPT diagram, have very high ionization parameters, and relatively low hardness ionizing radiation fields, exhibit nitrogen overabundances, and show extremely red mid-IR colors, and reside in the so-called "ULIRGs/LINERs/Obscured AGN" region. The blueberry galaxies may not harbor AGN.

We present the results from the application of a two-dimensional emission line detection method, EMission-line two-Dimensional (EM2D), to the near-infrared G102 grism observations obtained with the Wide-Field Camera 3 (WFC3) as part of the Cycle 22 {\em Hubble Space Telescope} Treasury Program: the Faint Infrared Grism Survey (FIGS). Using the EM2D method, we have assembled a catalog of emission line galaxies (ELGs) with resolved star formation from each of the four FIGS fields. Not only can one better assess the global properties of ELGs, but the EM2D method allows for the analysis and an improved study of the individual emission-line region {\it within} each galaxy. This paper includes a description of the methodology, advantages, and the first results of the EM2D method applied to ELGs in FIGS. The advantage of 2D emission line measurements includes significant improvement of galaxy redshift measurements, approaching the level of accuracy seen in high-spectral-resolution data, but with greater efficiency; and the ability to identify and measure the properties of multiple sites of star-formation and over scales of $\sim$ 1 kpc within individual galaxies out to z $\sim$ 4. The EM2D method also significantly improves the reliability of high-redshift ($z\sim7$) Lyman-$\alpha$ detections. Coupled with the wide field of view and high efficiency of space-based grism observations, EM2D provides a noteworthy improvement on the physical parameters that can be extracted from grism observations., Comment: Accepted by ApJ

We constrain the stellar population properties of a sample of 52 massive galaxies, with stellar mass log Ms>10.5, over the redshift range 0.5