Search

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

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

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures., Submitted in response to the 2018 NAS Exoplanet Science Strategy call. 5 pages, 2 figures

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

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time., Science White Paper submitted to the Astro2020 Decadal Survey. 5 pages plus references and cover page

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

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

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 $, 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.

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

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

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 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., 27 pages, 23 Figures, 7 tables. Accepted to ApJ

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

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

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

Lyman-alpha emitting galaxies can be used to study cosmological reionization, because a neutral intergalactic medium scatters Lyman-alpha photons into diffuse halos whose surface brightness falls below typical survey detection limits. Here we present the Lyman-alpha emitting galaxy LAE J095950.99+021219.1, identified at redshift z=6.944 in the COSMOS field using narrowband imaging and followup spectroscopy with the IMACS instrument on the Magellan I Baade telescope. With a single object spectroscopically confirmed so far, our survey remains consistent with a wide range of IGM neutral fraction at redshift seven, but further observations are planned and will help clarify the situation. Meantime, the object we present here is only the third Lyman-alpha selected galaxy to be spectroscopically confirmed at redshift seven, and is 2--3 times fainter than the previously confirmed redshift seven Lyman alpha galaxies., 15 pages including 3 figures. Accepted for publication in The Astrophysical Journal Letters

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 present a spectroscopically confirmed sample of Lyman alpha emitting galaxies (LAEs) at z ~ 4.5 in the Extended Chandra Deep Field South (ECDFS), which we combine with a sample of z ~ 4.5 LAEs from the Large Area Lyman Alpha (LALA) survey to build a unified Lya luminosity function (LF). We spectroscopically observed 64 candidate LAEs in the ECDFS, confirming 46 objects as z~4.5 LAEs. We did not detect significant flux from neither C_iv 1549�� nor the He_ii 1640�� emission in individual LAE spectra, even with a coadded spectrum. With the coadded line ratio of He_ii to Lya constraining the Population III star formation rate (SFR) to be 2$��$ significance, and the product $L^*��^*$ differs by a factor of >3. Similar LF differences have been used to infer evolution in the neutral gas fraction in the intergalactic medium at z>6, yet here the difference is likely due to cosmic variance, given that the two samples are from adjoining line-of-sight volumes. Combining our new sample of LAEs with those from previous LALA narrowband surveys at z = 4.5, we obtain one of the best measured Lya LFs to date of L* = 42.83 $\pm$ 0.06 and $��^*$ = -3.48 $\pm$ 0.09. We compare our new LF to others from the literature to study the evolution of the Lya luminosity density at 0 < z < 7. We find tentative evidence for evolution in the product $L^* ��^*$, which approximately tracks the cosmic SFR density, but since field-to-field and survey-to-survey variations are in some cases as large as the possible evolution, some caution is needed in interpreting this trend., 21 pages, 17 figures and 6 tables. Re-submitted to MNRAS to a new referee. Big changes with two referees' comments

Rapid mass assembly, likely from mergers or smooth accretion, has been predicted to play a vital role in star-formation in high-redshift Lyman-alpha (Lya) emitters. Here we predict the major merger, minor merger, and smooth accreting Lya emitter fraction from z~3 to z~7 using a large dark matter simulation, and a simple physical model that is successful in reproducing many observations over this large redshift range. The central tenet of this model, different from many of the earlier models, is that the star-formation in Lya emitters is proportional to the mass accretion rate rather than the total halo mass. We find that at z~3, nearly 35% of the Lya emitters accrete their mass through major (3:1) mergers, and this fraction increases to about 50% at z~7. This imply that the star-formation in a large fraction of high-redshift Lya emitters is driven by mergers. While there is discrepancy between the model predictions and observed merger fractions, some of this difference (~15%) can be attributed to the mass-ratio used to define a merger in the simulation. We predict that future, deeper observations which use a 3:1 definition of major mergers will find >30% major merger fraction of Lya emitters at redshifts >3., Accepted in MNRAS

We have coadded X-ray flux of all known Lyman Alpha Emitters (LAEs) in the 4 Msec Chandra Deep Field South (CDF-S) region, to place sensitive upper limits on the average unobscured star-formation rate (SFR_X) in these galaxies. A very small fraction of Lyman-Alpha galaxies in the field are individually detected in the X-rays, implying a low fraction of AGN activity. After excluding the few X-ray detected LAEs, we stack the undetected LAEs located in the 4 Ms CDF-S data and 250 ks Extended CDF-S (ECDFS) data, and compute a 1-��upper limit on SFR_X < 14, 28, 28, 140, 440, 880 M$_{\sun}$ yr$^{-1}$ for LAEs located at z = 2.1, 3.1, 3.2, 4.5, 5.7 and 6.5, respectively. The upper limit of SFR_X in LAEs can be then be compared to SFR$_{Ly��}$ derived from Lyman-Alpha line and thus can constrain on the Lyman-Alpha escape fraction ($f^{Esc}_{Ly��}$). We derive a lower limit on f(Lyman-Alpha Escape) > 14% (84 % confidence level, 1-��lower limit) for LAEs at redshift z ~ 2.1 and z ~ 3.1-3.2. At z > 4, the current LAE samples are not of sufficient size to constrain SFR_X well. By averaging all the LAEs at z> 2, the X-ray non-detection constrains f(Lyman-Alpha Escape) > 17% (84 % confidence level, 1-��lower limit), and rejects f(Lyman-Alpha Escape) < 5.7% at the 99.87% confidence level from 2.1 < z < 6.5., 11 pages, 3 tables, 3 figures, ApJ accepted

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)

Lyman alpha (Lya) emission lines should be attenuated in a neutral intergalactic medium (IGM). Therefore the visibility of Lya emitters at high redshifts can serve as a valuable probe of reionization at about the 50% level. We present an imaging search for z=7.7 Lya emitting galaxies using an ultra-narrowband filter (filter width= 9A) on the NEWFIRM imager at the Kitt Peak National Observatory. We found four candidate Lya emitters in a survey volume of 1.4 x 10^4 Mpc^3, with a line flux brighter than 6x10^-18 erg/cm^2/s (5 sigma in 2" aperture). We also performed a detailed Monte-Carlo simulation incorporating the instrumental effects to estimate the expected number of Lya emitters in our survey, and found that we should expect to detect one Lya emitter, assuming a non-evolving Lya luminosity function (LF) between z=6.5 and z=7.7. Even if one of the present candidates is spectroscopically confirmed as a z~8 Lya emitter, it would indicate that there is no significant evolution of the Lya LF from z=3.1 to z~8. While firm conclusions would need both spectroscopic confirmations and larger surveys to boost the number counts of galaxies, we successfully demonstrate the feasibility of sensitive near-infrared (1.06 um) narrow-band searches using custom filters designed to avoid the OH emission lines that make up most of the sky background., Published in ApJ, 3 figures

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies., Accepted by the Astrophysical Journal. 11 pages, 3 figures, 1 table (emulateapj)

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

We investigate the fraction of z ~ 0.3 Lyman alpha emitting galaxies (LAEs) which host active galactic nucleus activity, which is typically from 1 - 5% at 2 < z < 6. Using optical spectroscopy of 23 LAEs at 0.2 < z < 0.45 selected with GALEX UV data, we probe for AGN with a variety of methods, including line widths, diagnostic line ratios, high-ionization emission, infrared activity and X-ray luminosity. We found that our sample of low-redshift LAEs has an AGN fraction of 43 ^{+18}{-26}%, significantly higher than at high redshift. While previous results have shown that low-redshift LAEs have a lower space density than their high-redshift analogs, these results show that star-forming LAEs at low-redshift are rarer still. Accounting for differences in available AGN classification methods, we conclude that rest-frame optical spectroscopy is necessary to identify low-luminosity AGNs in LAEs at high redshift, and that limits on the X-ray flux are necessary to determine the extent to which the AGN contaminates the Lyman alpha fluxes., Replaced with ApJL accepted version, 5 pages, 2 figures

The Star Formation Camera (SFC) is a wide-field (~15'x19, >280 arcmin^2), high-resolution (18x18 mas pixels) UV/optical dichroic camera designed for the Theia 4-m space-borne space telescope concept. SFC will deliver diffraction-limited images at lambda > 300 nm in both a blue (190-517nm) and a red (517-1075nm) channel simultaneously. Our aim is to conduct a comprehensive and systematic study of the astrophysical processes and environments relevant for the births and life cycles of stars and their planetary systems, and to investigate and understand the range of environments, feedback mechanisms, and other factors that most affect the outcome of the star and planet formation process. This program addresses the origins and evolution of stars, galaxies, and cosmic structure and has direct relevance for the formation and survival of planetary systems like our Solar System and planets like Earth. We present the design and performance specifications resulting from the implementation study of the camera, conducted under NASA's Astrophysics Strategic Mission Concept Studies program, which is intended to assemble realistic options for mission development over the next decade. The result is an extraordinarily capable instrument that will provide deep, high-resolution imaging across a very wide field enabling a great variety of community science as well as completing the core survey science that drives the design of the camera. The technology associated with the camera is next generation but still relatively high TRL, allowing a low-risk solution with moderate technology development investment over the next 10 years. We estimate the cost of the instrument to be $390M FY08., 17 pages - Activity White Paper for the Astro2010 Decadal Survey Subcommittee on Programs

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