Search

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$α$ 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$α$ 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$α$ 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$α$ line depths comparable to the deepest $z=7$ narrow-band surveys ($L_{\rm{Ly}α}\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$α$ 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., 11 pages, 7 figures, Submitted to ApJ

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-$α$ 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., 14 pages, 5 figures

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-$��$ 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., 9 pages, 5 figures, and 1 table. Accepted by Nature Astronomy

We study the dependence of Ly$��$ 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$��$-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$��$) 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$��$ equivalent width (EW(Ly$��$)), and the Ly$��$ escape faction ($f^{Ly��}_{esc}$), suggesting that small UV-continuum radii are crucial for Ly$��$ 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$��$ photons, and that Ly$��$-emitters show constant characteristic size independent of their redshift., 10 pages, 4 figures, 2 tables, accepted by ApJ; doi:10.3847/1538-4357/abf833

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$α$ 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$α$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$σ$ upper limit for the ratio of NV/Ly$α$, $f_{NV}/f_{Lyα} \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$., 8 pages, 4 figures; Submitted to ApJ

We present a new measurement of the Ly$��$ 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$��$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$��$ 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$��$ 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��$., 19 pages, 11 figures; accepted by ApJ

The H$��$ 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$��$ EWs correlate with galaxy properties and how selection biases could alter such correlations using a $z = 0.47$ narrowband-selected sample of 1572 H$��$ emitters from the Ly$��$ 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$��$ 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$��compared to high-mass HAEs. Combining the intrinsic EW -- stellar mass correlation with an observed SMF correctly reproduces the observed H$��$ 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., 21 Pages, 13 Figures, 2 Tables, accepted to MNRAS

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$��$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$��$ 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$��$ 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}��}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$��$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $��^*=10^{-2.85}$ Mpc$^{-3}$, $��= -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $��^*=10^{-2.8}$ Mpc$^{-3}$, $��=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $��^*=10^{-2.91}$ Mpc$^{-3}$, $��= -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$��$ 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 $��_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $��_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $��_{\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$��$ surveys., 16 pages, 8 figures, Accepted for publication in ApJ

We present our measurements of the H$��$, [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$��$-, 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$��$ 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 $��^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$��$, $��^\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 $��^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $��$ 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}��} = 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$��$ 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}��_{\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$��$, [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$��$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD., 21 pages, 12 figures, 6 Tables. Accepted by MNRAS for publication on 6 Jan 2020

We have studied ultraviolet images of 40 Green Pea galaxies and 15 local Lyman Break Galaxy Analogs to understand the relation between Ly$��$ 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$��$ equivalent width (EW(Ly$��$)) and the Ly$��$ escape fraction ($f^{Ly��}_{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$��$) and $f^{Ly��}_{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$��$ photons., 8 pages, 3 figures, 2 tables. Accepted for publication in ApJ

Ly$��$ emission from galaxies can be utilized to characterize the ionization state in the intergalactic medium (IGM). We report our search for Ly$��$ emission at $z>7$ using a comprehensive Keck/MOSFIRE near-infrared spectroscopic dataset, as part of the Texas Spectroscopic Search for Ly$��$ 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$��$ emission lines detected ($>$4$��$) at $z>7$, significantly increasing the spectroscopically confirmed sample. Our sample includes large equivalent-width (EW) Ly$��$ emitters ($>$50��), and additional tentative Ly$��$ emission lines detected at 3 - 4$��$ from five additional galaxies. We constrain the Ly$��$ 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$��$ 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$��$) with recent measurements at $z \sim$ 7.6. The spatial distribution of the detected Ly$��$ emitters implies the presence of a potential highly ionized region at $z\sim7.55$ which hosts four Ly$��$ 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$��$ EW distribution values and reionization models., 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 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 present the largest-ever sample of 79 Ly$��$ 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$��$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly��}$ $>$ 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$��$ luminosity function., 23 pages, 11 figures, 2 tables, accepted for publication in ApJ

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

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., 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

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., 6 pages, 3 figures

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap.", Comment: White paper submitted to the National Academy Committee on an Exoplanet Science Strategy; 6 pages (typo fixed)

We present an empirical parameterization of the [NII]/H$��$ 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$��$ in current low-resolution grism and narrow-band observations to derive intrinsic H$��$ fluxes, and (iii) to reliably forecast the number counts of H$��$ 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$��$ on stellar mass and [OIII]/H$��$ observed in local galaxy samples. We find large variations in the [NII]/H$��$ 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$��$ luminosities. Specifically, measurements of the intrinsic H$��$ 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$��$ emitters for Euclid and WFIRST. We also study the impact of blended H$��$ and [NII] on the accuracy of measured spectroscopic redshifts., 18 pages, 11 figures, 1 table, accepted for publication in ApJ

We present the results of an unbiased search for Ly�� emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Ly�� emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Ly�� emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Ly�� emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4��) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028��m. We identify this lines as Ly�� at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Ly�� rest-frame equivalent width (EW Ly��) yet published at z > 7 (140.3 +/- 19.0��)., 12 pages, 6 figures, ApJ Submitted

We studied Lyman-$��$ (Ly$��$) escape in a statistical sample of 43 Green Peas with HST/COS Ly$��$ spectra. Green Peas are nearby star-forming galaxies with strong [OIII]$��$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$��$ line emitters with rest-frame Ly$��$ equivalent width $>20$ ��. The Ly$��$ profiles of Green Peas are diverse. The Ly$��$ escape fraction, defined as the ratio of observed Ly$��$ flux to intrinsic Ly$��$ flux, shows anti-correlations with a few Ly$��$ kinematic features -- both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Ly$��$ profile. Using properties measured from SDSS optical spectra, we found many correlations -- Ly$��$ escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Ly$��$ profiles with the HI shell radiative transfer model and found Ly$��$ escape fraction anti-correlates with the best-fit $N_{HI}$. Finally, we fit an empirical linear relation to predict Ly$��$ escape fraction from the dust extinction and Ly$��$ 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$��$ escape and to probe the IGM optical depth along the line of sight of each $z>7$ Ly$��$ emission line galaxy in the JWST era., 15 pages, 11 figures, 3 tables, machine-readable tables included. ApJ in-press

Green Peas are nearby analogs of high-redshift Ly$��$-emitting galaxies. To probe their Ly$��$ escape, we study the spatial profiles of Ly$��$ 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$��$ emission from their 2D COS spectra, and of UV continuum from both the 2D spectra and NUV images. The Ly$��$ emission shows more extended spatial profiles than the UV continuum in most Green Peas. The deconvolved Full Width Half Maximum (FWHM) of the Ly$��$ 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$��$ sizes than UV sizes. We also compare the spatial profiles of Ly$��$ photons at blueshifted and redshifted velocities in eight Green Peas with sufficient data quality, and find the blue wing of the Ly$��$ line has a larger spatial extent than the red wing in four Green Peas with comparatively weak blue Ly$��$ line wings. We show that Green Peas and MUSE $z=3-6$ LAEs have similar Ly$��$ and UV continuum sizes, which probably suggests starbursts in both low-z and high-z LAEs drive similar gas outflows illuminated by Ly$��$ light. Five Lyman continuum (LyC) leakers in this sample have similar Ly$��$ 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., 8 pages, 6 figures. ApJ accepted

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

Using Ly�� emission line as a tracer of high redshift star forming galaxies, hundreds of Ly�� 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�� 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�� radiative transfer predict asymmetric Ly�� line profiles with broad red wing in LAEs with outflows. Here we report a z ~ 5.7 Ly�� emission line with a broad red wing extending to > 1000 km/s relative to the peak of Ly�� 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., Accepted to ApJ

We present a sample of 33 spectroscopically confirmed z ~ 3.1 Ly$��$-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 �� 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$��$ lines in the galaxy with the new [O III] detection, finding that the Ly$��$ line is shifted 52 km s$^{-1}$ redwards of the [O III] line, which defines the systemic velocity of the galaxy., 38 pages, 27 figures, 4 tables, Accepted for publication in MNRAS

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

Ly�� line EWs provide important clues to the physical nature of high redshift LAEs. However, measuring the Ly�� 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�� 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�� luminosity function and the Ly�� 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�� EW distribution dN/dEW = N exp(-EW/W$_0$), we obtain W$_0$ = 167+/-44�� (uncorrected for IGM absorption of Ly��). We discuss the likely range of IGM absorption effects in light of recent measurements of Ly�� line profiles and velocity offsets. Our data are consistent with Ly�� 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�� EW distribution with redshift which goes as a power-law form of W$_0$ $\prop$(1+z)$^$��, with �� = 1.1$\pm$0.1 (0.6$\pm$0.1) if no IGM corrections are applied to the Ly�� line; or �� = 1.7$\pm$0.1 (1.2$\pm$0.1) after applying a maximal IGM-absorption correction to Ly�� line) for an exponential (a gaussian) EW distribution from z = 0.3 to 6.5., published on MNRAS

In the absence of spectra, fitting template model spectra to observed photometric fluxes, known as Spectral Energy Distribution (SED) fitting, has become the workhorse for identifying high-z galaxies. In this paper, we present an analysis of the most recent and possibly most distant galaxies discovered in the Hubble Ultra Deep Field using a more robust method of redshift estimation based on Markov Chain Monte Carlo fitting (MCMC), rather than relying on the redshift of "best fit" models obtained using common chi^2 minimization techniques. The advantage of MCMC fitting is the ability to accurately estimate the probability density function of the redshift, as well as any other input model parameters, allowing us to derive accurate credible intervals by properly marginalizing over all other input model parameters. We apply our method to 13 recently identified sources and show that, despite claims based on chi^2 minimization, none of these sources can be securely ruled out as low redshift interlopers given the low signal-to-noise of currently available observations. We estimate that there is an average probability of 21% that these sources are low redshift interlopers., 18 Pages, 2 Figures, 1 Table. Accepted to ApJ

Galaxies selected on the basis of their emission line strength show low metallicities, regardless of their redshifts. We conclude this from a sample of faint galaxies at redshifts between 0.6 < z < 2.4, selected by their prominent emission lines in low-resolution grism spectra in the optical with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) and in the near-infrared using Wide-Field Camera 3 (WFC3). Using a sample of 11 emission line galaxies (ELGs) at 0.6 < z < 2.4 with luminosities of -22 < M_B < -19, which have [OII], H��, and [OIII] line flux measurements from the combination of two grism spectral surveys, we use the R23 method to derive the gas-phase oxygen abundances: 7.5 < 12+log(O/H) < 8.5. The galaxy stellar masses are derived using Bayesian based Markov Chain Monte Carlo (\piMC^2) fitting of their Spectral Energy Distribution (SED), and span the mass range 8.1 < log(M_*/M_\sun) < 10.1. These galaxies show a mass-metallicity (M-L) and Luminosity-Metallicity (L-Z) relation, which is offset by --0.6 dex in metallicity at given absolute magnitude and stellar mass relative to the local SDSS galaxies, as well as continuum selected DEEP2 samples at similar redshifts. The emission-line selected galaxies most resemble the local "green peas" galaxies and Lyman-alpha galaxies at z~0.3 and z~2.3 in the M-Z and L-Z relations and their morphologies. The G-M_{20} morphology analysis shows that 10 out of 11 show disturbed morphology, even as the star-forming regions are compact. These galaxies may be intrinsically metal poor, being at early stages of formation, or the low metallicities may be due to gas infall and accretion due to mergers., 24 pages with 7 figures

We show that populations of high redshift galaxies grow more luminous as they grow in linear size. This is because the brightness per unit area has a distinct upper limit due to the self-regulating nature of star-formation. As a corollary, we show that the observed increase in characteristic luminosity of Lyman Break Galaxies (L*) with time can be explained by their increase in size, which scales as 1/H(z). In contrast, Lyman-alpha selected galaxies have a characteristic, constant, small size between redshift z=2.25 and 6.5. Coupled with a characteristic surface brightness, this can explain their non-evolving ultraviolet continuum luminosity function. This compact physical size seems to be a critical determining factor in whether a galaxy will show Lyman-alpha emission. We base these conclusions on new size measurements for a sample of about 170 Lyman-alpha selected galaxies with Hubble Space Telescope broad-band imaging, over the redshift range 2.25 < z < 6. We combine these with a similar number of Lyman-alpha selected galaxies with half-light radii from the literature, and derive surface brightnesses for the entire combined sample., Submitted to ApJ letters, comments welcome

The study of Ly-alpha emission in the high-redshift universe is a useful probe of the epoch of reionization, as the Ly-alpha line should be attenuated by the intergalactic medium (IGM) at low to moderate neutral hydrogen fractions. Here we present the results of a deep and wide imaging search for Ly-alpha emitters in the COSMOS field. We have used two ultra-narrowband filters (filter width of ~8-9 {\deg}A) on the NEWFIRM camera, installed on the Mayall 4m telescope at Kitt Peak National Observatory, in order to isolate Ly-alpha emitters at z = 7.7; such ultra-narrowband imaging searches have proved to be excellent at detecting Ly-alpha emitters. We found 5-sigma detections of four candidate Ly-alpha emitters in a survey volume of 2.8 x 10^4 Mpc^3 (total survey area ~760 arcmin^2). Each candidate has a line flux greater than 8 x 10^-18 erg s^-1 cm^-2. Using these results to construct a luminosity function and comparing to previously established Ly-alpha luminosity functions at z = 5.7 and z = 6.5, we find no conclusive evidence for evolution of the luminosity function between z = 5.7 and z = 7.7. Statistical Monte Carlo simulations suggest that half of these candidates are real z = 7.7 targets, and spectroscopic follow-up will be required to verify the redshift of these candidates. However, our results are consistent with no strong evolution in the neutral hydrogen fraction of the IGM between z = 5.7 and z = 7.7, even if only one or two of the z = 7.7 candidates are spectroscopically confirmed., Comment: 29 pages, 5 figures, accepted to ApJ (12/11)

We present the highest redshift detections of resolved Lyman alpha emission, using Hubble Space Telescope/ACS F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science program in the GOODS CDF-S. We detect Lyman alpha emission from three spectroscopically confirmed z = 4.4 Lyman alpha emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lyman alpha emission. Comparing the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the escape of Lyman alpha photons at high redshift. While our data do not support a positional offset between the Lyman alpha and rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are significantly larger in Lyman alpha than in the rest-frame UV continuum. This result is confirmed when comparing object sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is significantly less than observed in similar filters from the ground. These results together imply that the Lyman alpha emission is not strictly confined to its indigenous star-forming regions. Rather, the Lyman alpha emission is more extended, with the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lyman alpha photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing a significant role in galaxies at these redshifts., Comment: Submitted to the Astrophysical Journal. 11 pages, 10 figures

[abridged] We propose a tiered, UV--near-IR, cosmological broad- and medium-band imaging and grism survey that covers ~10 deg^2 in two epochs to m_AB=28, ~3 deg^2 in seven epochs to m_AB=28, and ~1 deg^2 in 20 epochs to m_AB=30 mag. Such a survey is an essential complement to JWST surveys (1100nm and z>~8). We aim to: (1) understand in how galaxies formed from perturbations in the primordial density field by studying faint Ly\alpha-emitting and Lyman-break galaxies at 5.5~5x10^6 galaxies; (5) trace the strongly epoch-dependent galaxy merger rate and constrain how Dark Energy affected galaxy assembly and the growth of SMBHs; (6) study >~10^5 weak AGN/feeding SMBHs in the faint-end of the QSO LF, over 10 deg^2 and measure how the growth of SMBHs kept pace with galaxy assembly and spheroid growth, and how this process was shaped by various feedback processes over cosmic times since z~8. The proposed study is not feasible with current instrumentation but argues for a wide-field (>~250 arcmin^2), high-resolution (, Comment: 8 pages, LaTeX2e requires 'fancyheadings_b.sty' (included), 4 postscript figures. Science White Paper submitted to the Astro2010 Decadal Survey. Compiled high-resolution PDF versions are available from either: http://www8.nationalacademies.org/astro2010/DetailFileDisplay.aspx?id=378 or: http://www.public.asu.edu/~rjansen/sfc/

We present a simple physical model for populating dark matter halos with Lyman Alpha Emiiters(LAEs) and predict the physical properties of LAEs at z~3-7. The central tenet of this model is that the Ly-alpha luminosity is proportional to the star formation rate (SFR) which is directly related to the halo mass accretion rate. The only free parameter in our model is then the star-formation efficiency (SFE). An efficiency of 2.5% provides the best-fit to the Ly-alpha luminosity function (LF) at redshift z=3.1, and we use this SFE to construct Ly-alpha LFs at other redshifts. Our model reproduce the Ly-alpha LFs, stellar ages, SFR ~1-10; Msun/yr, stellar masses ~ 10^7-10^8 Msun and the clustering properties of LAEs at z~3-7. We find the spatial correlation lengths ro ~ 3-6 Mpc/h, in agreement with the observations. Finally, we estimate the field-to-field variation ~ 30% for current volume and flux limited surveys, again consistent with observations. Our results suggest that the star formation, and hence Ly-alpha emission in LAEs is powered by the accretion of new material, and that the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating the accreted mass, rather than the total mass of halos, to the Ly-alpha luminosity of LAEs naturally gives rise to the duty cycle of LAEs., Comment: Published in ApJ

We present the faintest spectroscopically confirmed sample of redshift four and five Lyman break galaxies to date. The sample is based on slitless grism spectra of the Hubble Ultra Deep Field region from the GRAPES (Grism ACS Program for Extragalactic Science) and PEARS (Probing Evolution and Reionization Spectroscopically) projects, using the G800L grism on the HST Advanced Camera for Surveys. We report here confirmations of 39 galaxies, pre-selected as candidate Lyman break galaxies using photometric selection criteria. We compare a "traditional" V-dropout selection to a more liberal one (with V-i > 0.9), and find that the traditional criteria are about 64% complete and 81% reliable. We also study the Lyman alpha emission properties of our sample. We find that Lyman alpha emission is detected in about 1/4 of the sample, and that our broad-band color selected sample includes 55% of previously published line-selected Lyman alpha sources. Finally, we examine our stacked 2D spectra. We demonstrate that strong, spatially extended (arcsecond scale) Lyman alpha emission is not a generic property of these Lyman break galaxies, but that a modest extension of the Lyman alpha photosphere (compared to the starlight) may be present in those galaxies with prominent Lyman alpha emission., Submitted to The Astrophysical Journal. Reduced spectra from both GRAPES and PEARS are available from STScI, at http://www.stsci.edu/science/grapes/ and at http://archive.stsci.edu/prepds/pears/

We examine the stellar populations of a sample of 98 redshift 4.5 Lyman alpha emitting galaxies using their broadband colors derived from deep photometry at the MMT. These galaxies were selected by narrowband excess from the Large Area Lyman Alpha survey. Twenty-two galaxies are detected in two or more of our MMT filters (g', r', i' and z'). By comparing broad and narrowband colors of these galaxies to synthetic colors from stellar population models, we determine their ages and stellar masses. The highest equivalent width objects have an average age of 4 Myr, consistent with ongoing star formation. The lowest EW objects show an age of 40 - 200 Myr, consistent with the expectation that larger numbers of stars are causing low EWs. We found masses ranging from 2e7 solar masses for the youngest objects in the sample to 2e9 solar masses for the oldest. It is possible that dust effects could produce large equivalent widths even in older populations by allowing the Lyman alpha photons to escape, even while the continuum is extinguished, and we present models for this scenario also., Comment: 8 pages. Submitted to The Astrophysical Journal. 3 figures. Uses "emulateapj.sty". Figure 1 has been replaced with a corrected version

High redshift galaxies selected on the basis of their strong Lyman-alpha emission tend to be young ages and small physical sizes. We show this by analyzing the spectral energy distribution (SED) of 9 Lyman-alpha emitting (LAE) galaxies at 4.0 < z < 5.7 n the Hubble Ultra Deep Field (HUDF). Rest-frame UV to optical 700A < lambda < 7500A luminosities, or upper limits, are used to constrain old stellar populations. We derive best fit, as well as maximally massive and maximally old, properties of all 9 objects. We show that these faint and distant objects are all very young, being most likely only a few millions years old, and not massive, the mass in stars being ~10^6-10^8 M_sun. Deep Spitzer Infrared Array Camera (IRAC) observations of these objects, even in cases where objects were not detected, were crucial in constraining the masses of these objects. The space density of these objects, ~1.25x10^-4 Mpc^-3 is comparable to previously reported space density of LAEs at moderate to high redshifts. These Lyman-alpha galaxies show modest star formation rates of ~8 M_sun yr^-1, which is nevertheless strong enough to have allowed these galaxies to assemble their stellar mass in less than a few x10^6 years. These sources appear to have small physical sizes, usually smaller than 1 Kpc, and are also rather concentrated. They are likely to be some of the least massive and youngest high redshift galaxies observed to date., 27 pages. Accepted. To appear in The Astrophysical Journal

We present a redshift z=6.535 galaxy discovered by its Lyman alpha emission in a 9180A narrowband image from the Large Area Lyman Alpha (LALA) survey. The Ly-a line luminosity (1.1e43 erg/s) is among the largest known for star forming galaxies at z=6.5. The line shows the distinct asymmetry that is characteristic of high-redshift Ly-a. The 2 sigma lower bound on the observer-frame equivalent width is > 530A. This is hard to reconcile with a neutral intergalactic medium unless the Ly-a line is intrinsically strong and is emitted from its host galaxy with an intrinsic Doppler shift of several hundred km/s. If the IGM is ionized, it corresponds to a rest frame equivalent width > 40A after correcting for Ly-a forest absorption. We also present complete spectroscopic followup of the remaining candidates with line flux > 2e-17 erg/cm2/s in our 1200 square arcminute narrowband image. These include another galaxy with a strong emission line at 9136A and no detected continuum flux, which however is most likely an [OIII] source at z=0.824 based on a weak detection of the [OIII] 4959A line., Revised version, accepted to The Astrophysical Journal, 22 April 2004 (original version submitted 27 February 2004)

We present Keck/LRIS optical spectra of 17 Lya-emitting galaxies and one Lyman break galaxy at z=4.5 discovered in the Large Area Lyman Alpha (LALA) survey. The survey has identified a sample of ~350 candidate Lya-emitting galaxies at z=4.5 in a search volume of 1.5 x 10^6 comoving Mpc^3. We targeted 25 candidates for spectroscopy; hence, the 18 confirmations presented herein suggest a selection reliability of 72%. The large equivalent widths (median W(rest)~80 A) but narrow physical widths (v < 500 km/s) of the Lya emission lines, along with the lack of accompanying high-ionization state emission lines, suggest that these galaxies are young systems powered by star formation rather than by AGN activity. Theoretical models of galaxy formation in the primordial Universe suggest that a small fraction of Lya-emitting galaxies at z=4.5 may still be nascent, metal-free objects. Indeed, we find with 90% confidence that 3 to 5 of the confirmed sources show W(rest) > 240 A, exceeding the maximum Lya equivalent width predicted for normal stellar populations. Nonetheless, we find no evidence for HeII 1640 emission in either individual or composite spectra, indicating that though these galaxies are young, they are not truly primitive, Population III objects., 12 pages, Accepted to ApJ

We report spectroscopic detection of three z~6 Lyman-alpha emitting galaxies, in the vicinity of the Hubble Ultra Deep Field, from the early data of the Gemini Lyman-$��$ at Reionisation Era (GLARE) project. Two objects, GLARE#3001 (z =5.79) and GLARE#3011 (z =5.94), are new detections and are fainter in $z'$ (z'_{AB} =26.37 and 27.15) than any Lyman break galaxy previously detected in Lyman-alpha. A third object, GLARE#1042 (z =5.83) has previously been detected in line emission from the ground; we report here a new spectroscopic continuum detection. Gemini/GMOS-S spectra of these objects, obtained using nod & shuffle, are presented together with a discussion of their photometric properties. All three objects were selected for spectroscopy via the i-drop Lyman Break technique, the two new detections from the GOODS v1.0 imaging data. The red i'-z' colors and high equivalent widths of these objects suggest a high-confidence z>5 Lyman-alpha identification of the emission lines. This brings the total number of known z>5 galaxies within 9 arcmin of the Hubble Ultra Deep Field to four, of which three are at the same redshift (z=5.8 within 2000 km/s suggesting the existence of a large-scale structure at this redshift., 5 pages, 2 figures. Revised to match accepted version

The Lyman-$\alpha$ line is expected to be strong in the presence of active star formation and the absence of dust, making it a good tool for finding chemically primitive galaxies in the early universe. We report on a new survey for high redshift Lyman-$\alpha$ sources, the Large Area Lyman Alpha (LALA) survey. Our survey achieves an unprecedented combination of volume and sensitivity by using narrow-band filters on the new $8192^2$ pixel CCD Mosaic Camera at the 4 meter Mayall telescope of Kitt Peak National Observatory. Well-detected sources with flux and equivalent width matching known high redshift Lyman-$\alpha$ galaxies have an observed surface density corresponding to $11000 \pm 700$ per square degree per unit redshift at $z=4.5$. Early spectroscopic followup from the Keck telescope suggests that $\sim 1/3$ of these are actually at $z\approx 4.5$, and has confirmed five $z > 4$ Lyman-$\alpha$ emitters so far. Combining our photometric survey with spectroscopic results, we estimate a net density of $\sim 4000 $ Lyman-$\alpha$ emitters per square degree per unit redshift at $z\approx 4.5$. The star formation rate density (estimated both from UV continuum and from line emission) is comparable to that of the Lyman break galaxy population within present uncertainties. The most extreme Lyman-$\alpha$ emitters in our sample have rest frame equivalent widths $> 100\AA$, consistent with the expectations for the first burst of star formation in a primitive, dust-free galaxy., Comment: 6 pages; proceedings from "The New Era of Wide-Field Astronomy," Preston, UK

ISO has enabled far-infrared spectroscopy of a variety of galaxies. Using the [CII] (158 micron) and [OI](63 micron) lines, we can characterize the physical conditions in the star-forming ISM. These observations also form the basis of our expectations for what the Herschel Space Observatory will observe for high redshift galaxies. While [CII] is suppressed in ULIRGs and normal galaxies with high dust temperatures, it is stronger than expected in metal poor galaxies by factors of a few. Young galaxies at high redshifts might be expected to be both metal poor and actively star-forming, leading to contrary expectations for the [CII] line strength. The best prediction for [CII] detection is derived by using the observed proportionality between [CII] and mid-IR emission from PAHs. Using the observed [CII]/7 micron ratio and number counts from ISO deep surveys we predict that HSO will be able to detect 100 sources/square-degree in the [CII] line., Comment: 4 Pages, uses toledoStyle.cls, to appear in `The Promise of FIRST', ESA SP-460, eds. G.L. Pilbratt et al

The Large Area Lyman Alpha survey has found ~ 150 Lyman-alpha emitters at z=4.5. While stellar models predict a maximum Lyman-alpha equivalent width (EW) of 240 angstrom, 60% of the Lyman-alpha emitters have EWs exceeding this value. We attempt to model the observed EW distribution by combining stellar population models with an extrapolation of Lyman break galaxy luminosity function at z=4, incorporating observational selection effects and Malmquist bias. To reproduce the high EWs seen in the sample we need to postulate a stellar initial mass function (IMF) with extreme slope alpha = 0.5 (instead of 2.35); zero metallicity stars; or narrow-lined active galactic nuclei. The models also reveal that only 7.5-15% of galaxies need show Lyman-alpha emission to explain the observed number counts. This raises the possibility that either star-formation in high redshift galaxies is episodic or the Lyman-alpha galaxies we are seeing are the youngest 7.5-15% and that Lyman-alpha is strongly quenched by dust at about 10 Mega-years of age., Comment: submitted to ApJ letters

Strong Lyman-alpha emission is a signpost of young stars and the absence of dust and thus indicates young galaxies. To find such a population of young galaxies at z=4.5 we started the Large Area Lyman Alpha survey (LALA). This survey achieves an unprecedented combination of volume and sensitivity by using narrow-band filters on a large format (36x36 square-arcminutes) camera on the 4 meter telescope at KPNO. The volume density and star-formation contribution of the Lyman-alpha emitters at z=4.5 is comparable to that of Lyman break galaxies. With many candidates and a few spectroscopic confirmations in hand we discuss what the properties of Ly-alpha emitters imply for galaxy and star formation in the early universe., Comment: 5 pages, 1 figure, To appear in ``Gas and Galaxy Evolution,'' eds. J.E. Hibbard, M.P. Rupen, & J.H. van Gorkom, ASP Conference Series (VLA 20th Anniversary Conference)

Lyman alpha galaxies at high redshifts offer a powerful probe of both the formation of galaxies and the reionization of the intergalactic medium. Lyman alpha line emission is an efficient tool for identifying young galaxies at high redshift, because it is strong in systems with young stars and little or no dust-- properties expected in galaxies undergoing their first burst of star-formation. Lyman alpha galaxies also provide a robust test of the reionization epoch that is independent of Gunn-Peterson trough observations in quasar spectra and is better able to distinguish line center optical depths tau=5 from tau=10^5. This is because neutral gas scatters Lyman alpha photons, dramatically ``blurring'' images of Lyman alpha galaxies embedded in a neutral intergalactic medium and rendering them undetectable. We present a photometrically selected sample of z=5.7 Lyman alpha emitters derived from the Large Area Lyman Alpha survey. The presence of these low-luminosity Lyman alpha sources at z=5.7 immediately implies that the reionization redshift was > 5.7. Comparing these objects to our earlier z=4.5 sample, we find that the number of z=5.7 emitters at fixed line luminosity marginally exceeds the no-evolution expectation, but falls well short of published model predictions. The equivalent width distribution is similar at the two redshifts. The large equivalent widths of the Lyman alpha line indicate young galaxies undergoing their first star formation., 10 pages, uses AASTeX macros; revised version, accepted by ApJ Letters

We report a strong correlation between the [CII] 158 micron cooling line and the mid-infrared flux in the 5-10 micron range in a wide variety of star-forming galaxies. The mid-infrared flux is dominated by Aromatic Feature Emission (AFE), which is thought to arise from large polycyclic aromatic hydrocarbon molecules or `PAHs' and generally associated with the smallest interstellar grains. The [CII] line is the dominant gas coolant in most regions of atomic interstellar gas, and therefore reflects the heating input to the gas. The ratio of these two quantities, [CII]/AFE, remains nearly constant with the ratio of the IRAS 60 micron band flux to the 100 micron band flux, R(60/100). This is in contrast to the drop in the [CII]/FIR ratio with increasing R(60/100), which signal higher dust temperatures and more intense radiation fields. We interpret the stable [CII]/AFE ratio as evidence that gas heating is dominated by the PAHs or small grains which are also AFE carriers over a wide range of conditions. The trend of decreasing [CII]/FIR with increasing radiation field suggests a decrease in the importance of PAHs or small grains relative to large grains both in gas heating and in dust cooling. We summarize the observed trends and suggest two plausible scenarios., Comment: 10 pages including 2 figures; accepted for publication in ApJ Letters

Using the near infrared fluxes of local galaxies derived from COBE/DIRBE J(1.25 $��$m) K (2.2 $��$m) \& L (3.5 $��$m) band maps and published Cepheid distances, we construct Tully-Fisher (TF) diagrams for the nearby galaxies. The measured dispersions in these luminosity-linewidth diagrams are remarkably small: $��_J = 0.09$ magnitudes, $��_K = 0.13$ magnitudes, and $��_L = 0.20$ magnitudes. For the J and K bands, Monte Carlo simulations give a 95\% confidence interval upper limit on the true scatter in the T-F diagram of $��_J \le 0.35$ and $��_K \le 0.45$. We determine Milky Way's luminosity and place it in the T-F diagram by fitting a bar plus ``standard'' values of its size and circular speed (Sun-Galactic center distance $R_0 = 8.5 \kpc$ and $��_0 =220 \kms$), the Milky Way lies within $1.5 ��$ of the TF relations. We can use the TF relation and the Cepheid distances to nearby bright galaxies to constrain $R_0$ and $��_0$. Alternatively, we can fix the parameters of the Galaxy to their standard values, ignore the Cepheid zero-point, and use the Tully-Fisher relation to determine the Hubble Constant directly: $H_0 = 66 \pm 12$ km/s/Mpc., 11 pages, aaspp macro, 3 figures, submitted to the Astrophysical Journal. Also available at http://www.astro.princeton.edu/~san/gal.html

Photometric monitoring of the gravitational lens system 0957+561A,B in the g and r bands with the Apache Point Observatory (APO) 3.5 m telescope during 1996 shows a sharp g band event in the trailing (B) image light curve at the precise time predicted from the observation of an event during 1995 in the leading (A) image with a delay of 415 days. This success confirms the "short delay," and the lack of any feature at a delay near 540 days rejects the "long delay" for this system, resolving a long-standing controversy. A series of statistical analyses of our light curve data yield a best fit delay of 417 +/- 3 days (95% confidence interval). Recent improvements in the modeling of the lens system (consisting of a galaxy and cluster) allow us to derive a value of the global (at z = 0.36) value of Hubble's constant H_0 using Refsdal's method, a simple and direct distance determination based on securely understood physics and geometry. The result is H_0 = 63 +/- 12 km/s/Mpc (for Omega = 1) where this 95% confidence interval is dominated by remaining lens model uncertainties., accepted by ApJ, AASTeX 4.0 preprint, 4 PostScript figures