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

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 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 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

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

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

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 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 an empirical parameterization of the [NII]/H$\alpha$ flux ratio as a function of stellar mass and redshift valid at 0 < z < 2.7 and 8.5 < log(M) < 11.0. This description can easily be applied to (i) simulations for modeling [NII]$\lambda6584$ line emission, (ii) deblend [NII] and H$\alpha$ in current low-resolution grism and narrow-band observations to derive intrinsic H$\alpha$ fluxes, and (iii) to reliably forecast the number counts of H$\alpha$ emission-line galaxies for future surveys, such as those planned for Euclid and the Wide Field Infrared Survey Telescope (WFIRST). Our model combines the evolution of the locus on the Baldwin, Phillips & Terlevich (BPT) diagram measured in spectroscopic data out to z = 2.5 with the strong dependence of [NII]/H$\alpha$ on stellar mass and [OIII]/H$\beta$ observed in local galaxy samples. We find large variations in the [NII]/H$\alpha$ flux ratio at a fixed redshift due to its dependency on stellar mass; hence, the assumption of a constant [NII] flux contamination fraction can lead to a significant under- or overestimate of H$\alpha$ luminosities. Specifically, measurements of the intrinsic H$\alpha$ luminosity function derived from current low-resolution grism spectroscopy assuming a constant 29% contamination of [NII] can be overestimated by factors of ~8 at log(L) > 43.0 for galaxies at redshifts z = 1.5. This has implications for the prediction of H$\alpha$ emitters for Euclid and WFIRST. We also study the impact of blended H$\alpha$ and [NII] on the accuracy of measured spectroscopic redshifts., Comment: 18 pages, 11 figures, 1 table, accepted for publication in ApJ

Searching for extreme emission line galaxies allows us to find low-mass metal-poor galaxies that are good analogs of high redshift Ly$\alpha$ emitting galaxies. These low-mass extreme emission line galaxies are also potential Lyman-continuum leakers. Finding them at very low redshifts ($z\lesssim0.05$) allows us to be sensitive to even lower stellar masses and metallicities. We report on a sample of extreme emission line galaxies at $z\lesssim0.05$ (blueberry galaxies). We selected them from SDSS broadband images on the basis of their broad band colors, and studied their properties with MMT spectroscopy. From the whole SDSS DR12 photometric catalog, we found 51 photometric candidates. We spectroscopically confirm 40 as blueberry galaxies. (An additional 7 candidates are contaminants, and 4 remain without spectra.) These blueberries are dwarf starburst galaxies with very small sizes ($< 1\hbox{kpc}$), and very high ionization ([OIII]/[OII]$\sim10-60$). They also have some of the lowest stellar masses ($\log(\hbox{M}/\hbox{M}_{\odot})\sim6.5-7.5$) and lowest metallicities ($7.1<12+\log(\hbox{O/H})<7.8$) starburst galaxies. Thus they are small counterparts to green peas and high redshift Ly$\alpha$ emitting galaxies., Comment: 9 pages, 8 figures, 2 tables updated. Machine-readable table included in the source. ApJ accepted

We studied Lyman-$\alpha$ (Ly$\alpha$) escape in a statistical sample of 43 Green Peas with HST/COS Ly$\alpha$ spectra. Green Peas are nearby star-forming galaxies with strong [OIII]$\lambda$5007 emission lines. Our sample is four times larger than the previous sample and covers a much more complete range of Green Pea properties. We found that about 2/3 of Green Peas are strong Ly$\alpha$ line emitters with rest-frame Ly$\alpha$ equivalent width $>20$ \AA. The Ly$\alpha$ profiles of Green Peas are diverse. The Ly$\alpha$ escape fraction, defined as the ratio of observed Ly$\alpha$ flux to intrinsic Ly$\alpha$ flux, shows anti-correlations with a few Ly$\alpha$ kinematic features -- both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Ly$\alpha$ profile. Using properties measured from SDSS optical spectra, we found many correlations -- Ly$\alpha$ escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Ly$\alpha$ profiles with the HI shell radiative transfer model and found Ly$\alpha$ escape fraction anti-correlates with the best-fit $N_{HI}$. Finally, we fit an empirical linear relation to predict Ly$\alpha$ escape fraction from the dust extinction and Ly$\alpha$ red peak velocity. The standard deviation of this relation is about 0.3 dex. This relation can be used to isolate the effect of IGM scatterings from Ly$\alpha$ escape and to probe the IGM optical depth along the line of sight of each $z>7$ Ly$\alpha$ emission line galaxy in the JWST era., Comment: 15 pages, 11 figures, 3 tables, machine-readable tables included. ApJ in-press

Green Peas are nearby analogs of high-redshift Ly$\alpha$-emitting galaxies. To probe their Ly$\alpha$ escape, we study the spatial profiles of Ly$\alpha$ and UV continuum emission of 24 Green Pea galaxies using the Cosmic Origins Spectrograph (COS) on Hubble Space Telescope (HST). We extract the spatial profiles of Ly$\alpha$ emission from their 2D COS spectra, and of UV continuum from both the 2D spectra and NUV images. The Ly$\alpha$ emission shows more extended spatial profiles than the UV continuum in most Green Peas. The deconvolved Full Width Half Maximum (FWHM) of the Ly$\alpha$ spatial profile is about 2 to 4 times that of the UV continuum in most cases. Since Green Peas are analogs of high-z LAEs, it suggests that most high-z LAEs likely have larger Ly$\alpha$ sizes than UV sizes. We also compare the spatial profiles of Ly$\alpha$ photons at blueshifted and redshifted velocities in eight Green Peas with sufficient data quality, and find the blue wing of the Ly$\alpha$ line has a larger spatial extent than the red wing in four Green Peas with comparatively weak blue Ly$\alpha$ line wings. We show that Green Peas and MUSE $z=3-6$ LAEs have similar Ly$\alpha$ and UV continuum sizes, which probably suggests starbursts in both low-z and high-z LAEs drive similar gas outflows illuminated by Ly$\alpha$ light. Five Lyman continuum (LyC) leakers in this sample have similar Ly$\alpha$ to UV continuum size ratios (~1.4-4.3) to the other Green Peas, indicating their LyC emission escape through ionized holes in the interstellar medium., Comment: 8 pages, 6 figures. ApJ accepted

We present a narrowband survey with three adjacent filters for z=2.8--2.9 Lyman Alpha Emitter (LAE) galaxies in the Extended Chandra Deep Field South (ECDFS), along with spectroscopic followup. With a complete sample of 96 LAEs in the narrowband NB466, we confirm a large-scale structure at z~ 2.8. Compared to the blank field in NB470 and NB475, the LAE density excess in the NB466 field is ~6.0+/-0.8 times the standard deviation expected at z~2.8, assuming a linear bias of 2. The overdense large scale structure in NB466 can be decomposed into 4 protoclusters, whose overdensities are 4.6 - 6.6. These 4 protoclusters are expected to evolve into a Coma-like cluster at z~ 0. In the meanwhile, we investigate the average star-formation rates derived from Ly{\alpha}, rest-frame UV and X-ray, the Ly{\alpha} luminosity functions, the Ly{\alpha} photon densities and their dependence on the environment. We find that the Ly{\alpha} photon density in the overdense field (NB466) is ~50\% higher than that in the blank field (NB470+NB475). The 3 brightest LAEs, including a quasar at z=2.81, are all detected in X-ray and in NB466. These three LAE-AGNs contribute an extra 20--30\% Ly{\alpha} photon density. Furthermore, we find that LAEs in overdense regions are younger and less dusty. We conclude that the structure we found is a significant and rare density peak, and narrowband imaging is an efficient method to detect and study such structures in the high-z universe., Comment: 22 pages, 20 figures, 5 tables. Accepted to ApJ Supplement

Deep rest-frame optical spectroscopy is critical for characterizing and understanding the physical conditions and properties of the ionized gas in galaxies. Here, we present a new spectroscopic survey called "Metal Abundances across Cosmic Time" or $\mathcal{MACT}$, which will obtain rest-frame optical spectra for $\sim$3000 emission-line galaxies. This paper describes the optical spectroscopy that has been conducted with MMT/Hectospec and Keck/DEIMOS for $\approx$1900 $z=0.1-1$ emission-line galaxies selected from our narrowband and intermediate-band imaging in the Subaru Deep Field. In addition, we present a sample of 164 galaxies for which we have measured the weak [OIII]$\lambda$4363 line (66 with at least 3$\sigma$ detections and 98 with significant upper limits). This nebular emission line determines the gas-phase metallicity by measuring the electron temperature of the ionized gas. This paper presents the optical spectra, emission-line measurements, interstellar properties (e.g., metallicity, gas density), and stellar properties (e.g., star formation rates, stellar mass). Paper II of the $\mathcal{MACT}$ survey (Ly et al.) presents the first results on the stellar mass--gas metallicity relation at $z\lesssim1$ using the sample with [OIII]$\lambda$4363 measurements., Comment: 64 pages, 33 figures, 17 tables. Updated to match published version in the Astrophysical Journal Supplement Series

We present the results from a stellar population modeling analysis of a sample of 162 z=4.5, and 14 z=5.7 Lyman alpha emitting galaxies (LAEs) in the Bootes field, using deep Spitzer/IRAC data at 3.6 and 4.5 um from the Spitzer Lyman Alpha Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 um for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z=4.5 LAEs and 9/14 (64%) of the z=5.7 LAEs are detected at >3-sigma in at least one IRAC band. Individual z=4.5 IRAC-detected LAEs have a large range of stellar mass, from 5x10^8 to 10^11 Msol. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr - 1 Gyr, while the remainder have ages < 100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8 -- 20 x 10^8 Msol) and young (64 - 570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star-formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z=4.5 LAE correlation is elevated 4-5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyman alpha photon escape., Comment: 27 pages, 23 Figures, 7 tables. Accepted to ApJ

We analyze archival Ly$\alpha$ spectra of 12 "Green Pea" galaxies observed with the Hubble Space Telescope, model their Ly$\alpha$ profiles with radiative transfer models, and explore the dependence of Ly$\alpha$ escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies with [OIII]$\lambda$5007 equivalent widths of hundreds of \AA. All 12 Green Pea galaxies in our sample show Ly$\alpha$ lines in emission, with a Ly$\alpha$ equivalent width distribution similar to high redshift Ly$\alpha$ emitters. Combining the optical and UV spectra of Green Pea galaxies, we estimate their Ly$\alpha$ escape fractions and find correlations between Ly$\alpha$ escape fraction and kinematic features of Ly$\alpha$ profiles. The escape fraction of Ly$\alpha$ in these galaxies ranges from 1.4% to 67%. We also find that the Ly$\alpha$ escape fraction depends strongly on metallicity and moderately on dust extinction. We compare their high-quality Ly$\alpha$ profiles with single HI shell radiative transfer models and find that the Ly$\alpha$ escape fraction anti-correlates with the derived HI column densities. Single shell models fit most Ly$\alpha$ profiles well, but not the ones with highest escape fractions of Ly$\alpha$. Our results suggest that low HI column density and low metallicity are essential for Ly$\alpha$ escape, and make a galaxy a Ly$\alpha$ emitter., Comment: 9 figures, ApJ accepted

We report on two regularly rotating galaxies at redshift z=2, using high resolution spectra of the bright [CII] 158 micron emission line from the HIFI instrument on the Herschel Space Observatory. Both SDSS090122.37+181432.3 ("S0901") and SDSS J120602.09+514229.5 ("the Clone") are strongly lensed and show the double-horned line profile that is typical of rotating gas disks. Using a parametric disk model to fit the emission line profiles, we find that S0901 has a rotation speed v sin(i) = 120 +/- 7 km/s and gas velocity dispersion sigma < 23 km/s. The best fitting model for the Clone is a rotationally supported disk having v sin(i) = 79 +/- 11 km/s and sigma < 4km/s. However the Clone is also consistent with a family of dispersion-dominated models having sigma = 92 +/- 20 km/s. Our results showcase the potential of the [CII] line as a kinematic probe of high redshift galaxy dynamics: [CII] is bright; accessible to heterodyne receivers with exquisite velocity resolution; and traces dense star-forming interstellar gas. Future [CII] line observations with ALMA would offer the further advantage of spatial resolution, allowing a clearer separation between rotation and velocity dispersion., Comment: 20 pages, 4 figures; in press at The Astrophysical Journal

We present a sample of 33 spectroscopically confirmed z ~ 3.1 Ly$\alpha$-emitting galaxies (LAEs) in the Cosmological Evolution Survey (COSMOS) field. This paper details the narrow-band survey we conducted to detect the LAE sample, the optical spectroscopy we performed to confirm the nature of these LAEs, and a new near-infrared spectroscopic detection of the [O III] 5007 \AA\ line in one of these LAEs. This detection is in addition to two [O III] detections in two z ~ 3.1 LAEs we have reported on previously (McLinden et al 2011). The bulk of the paper then presents detailed constraints on the physical characteristics of the entire LAE sample from spectral energy distribution (SED) fitting. These characteristics include mass, age, star-formation history, dust content, and metallicity. We also detail an approach to account for nebular emission lines in the SED fitting process - wherein our models predict the strength of the [O III] line in an LAE spectrum. We are able to study the success of this prediction because we can compare the model predictions to our actual near-infrared observations both in galaxies that have [O III] detections and those that yielded non-detections. We find a median stellar mass of 6.9 $\times$ 10$^8$ M$_{\odot}$ and a median star formation rate weighted stellar population age of 4.5 $\times$ 10$^6$ yr. In addition to SED fitting, we quantify the velocity offset between the [O III] and Ly$\alpha$ lines in the galaxy with the new [O III] detection, finding that the Ly$\alpha$ line is shifted 52 km s$^{-1}$ redwards of the [O III] line, which defines the systemic velocity of the galaxy., Comment: 38 pages, 27 figures, 4 tables, Accepted for publication in MNRAS

Using Ly{\alpha} emission line as a tracer of high redshift star forming galaxies, hundreds of Ly{\alpha} emission line galaxies (LAEs) at z > 5 have been detected. These LAEs are considered to be low mass young galaxies, critical to the reionization of the universe and the metal enrichment of circumgalactic medium (CGM) and intergalactic medium (IGM). It is assumed that outflows in LAEs can help ionizing photons and Ly{\alpha} photons escape out of galaxies. However we still know little about the outflows in high redshifts LAEs due to observational difficulties, especially at redshift > 5. Models of Ly{\alpha} radiative transfer predict asymmetric Ly{\alpha} line profiles with broad red wing in LAEs with outflows. Here we report a z ~ 5.7 Ly{\alpha} emission line with a broad red wing extending to > 1000 km/s relative to the peak of Ly{\alpha} line, which has been detected in only a couple of z > 5 LAEs till now. If the broad red wing is ascribed to gas outflow instead of AGN activity, the outflow velocity could be larger than the escape velocity (~ 500 km/s) of typical halo mass of z ~ 5.7 LAEs, being consistent with the picture that outflows in LAEs disperse metals to CGM and IGM., Comment: Accepted to ApJ

Ly{\alpha} line EWs provide important clues to the physical nature of high redshift LAEs. However, measuring the Ly{\alpha} EW distribution of high-z narrowband selected LAEs can be hard because many sources do not have broadband photometry. We investigate the possible biases in measuring the intrinsic Ly{\alpha} EW distribution for a LAE sample at z $\sim$ 4.5 in the Extended CDFS. Only weak Malmquist-type bias in both the intrinsic Ly{\alpha} luminosity function and the Ly{\alpha} EW distribution were found. However, the observed EW distribution is severely biased if one only considers LAEs with detections in the continuum. Taking the broadband non-detections into account requires fitting the distribution of the broadband-to-narrowband ratio, which then gives a larger EW distribution scale length. Assuming an exponential form of the intrinsic Ly{\alpha} EW distribution dN/dEW = N exp(-EW/W$_0$), we obtain W$_0$ = 167+/-44{\AA} (uncorrected for IGM absorption of Ly{\alpha}). We discuss the likely range of IGM absorption effects in light of recent measurements of Ly{\alpha} line profiles and velocity offsets. Our data are consistent with Ly{\alpha} EW being independent of UV luminosity (i.e., we do not see evidence for the "Ando" effect). Our simulations also imply that broad-band images should be 0.5-1 magnitude deeper than narrowband images for an effective and reasonably complete LAE survey. Comparing with consistent measurements at other redshifts, we see a strong evolution in Ly{\alpha} EW distribution with redshift which goes as a power-law form of W$_0$ $\prop$(1+z)$^${\xi}, with {\xi} = 1.1$\pm$0.1 (0.6$\pm$0.1) if no IGM corrections are applied to the Ly{\alpha} line; or {\xi} = 1.7$\pm$0.1 (1.2$\pm$0.1) after applying a maximal IGM-absorption correction to Ly{\alpha} line) for an exponential (a gaussian) EW distribution from z = 0.3 to 6.5., Comment: published on MNRAS

Motivated by recent observations of galaxies dominated by emission lines, which show evidence of being metal poor with young stellar populations, we present calculations of multiple model grids with a range of abundances, ionization parameters, and stellar ages, finding that the predicted spectral line diagnostics are heavily dependent on all three parameters. These new model grids extend the ionization parameter to larger values than typically explored. We compare these model predictions with previous observations of such objects, including two new Lyman-$\alpha$ emitting galaxies (LAE) that we have observed. Our models give improved constraints on the metallicity and ionization parameter of these previously studied objects, as we are now able to consider high ionization parameter models. However, similar to previous work, these models have difficulty predicting large line diagnostics for high ionization potential species, requiring future work refining the modelling of FUV photons. Our model grids are also able to constrain the metallicity and ionization parameter of our LAEs, and give constraints on their Ly$\alpha$ escape fractions, all of which are consistent with recent lower redshift studies of LAEs., Comment: 20 pages, 13 figures, 1 table. Submitted to AJ. We encourage comments and suggestions

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground., Comment: 19 pages, four figures. Resubmitted to The Astrophysical Journal after revisions to address the referee's report