Search

Substantial populations of massive quiescent galaxies at $z\ge3$ challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at $z>3$. As a result, studies report that as much as 70\% of quiescent galaxies at $z>3$ may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at $3\lesssim z \lesssim 6$ using JWST NIRCam imaging data. We use empirically-constrained galaxy SED templates to define a region in the $F277W-F444W$ vs. $F150W-F277W$ color plane that captures quiescent galaxies at $z>3$. We apply this color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and identify 44 candidate $z\gtrsim3$ quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of post-starburst galaxies. We derive volume density estimates of $n\sim1-4\times10^{-5}$\,Mpc$^{-3}$ at $3< z, 12 pages, 4 figures, 2 tables, submitted to ApJ

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$α$-break by Hubble and with a Ly$α$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$α$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$β$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $μ$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., 30 pages, 13 figures, 6 tables, submitted to ApJ

We present JWST NIRSpec spectroscopy for 11 galaxy candidates with photometric redshifts of $z\simeq9-13$ and $M_{\rm\,UV} \in[-21,-18]$ newly identified in NIRCam images in the Cosmic Evolution Early Release Science (CEERS) Survey. We confirm emission line redshifts for 7 galaxies at $z=7.762-8.998$ using spectra at $\sim1-5μ$m either with the NIRSpec prism or its three medium resolution gratings. For $z\simeq9$ photometric candidates, we achieve a high confirmation rate of $\simeq$90\%, which validates the classical dropout selection from NIRCam photometry. No robust emission lines are identified in three galaxy candidates at $z>10$, where the strong [OIII] and H$β$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$α$ continuum break is not detected with the current sensitivity. Compared with HST-selected bright galaxies ($M_{\rm\,UV}\simeq-22$) that are similarly spectroscopically confirmed at $z\gtrsim8$, these NIRCam-selected galaxies are characterized by lower star formation rates (SFR$\simeq4\,M_{\odot}$~yr$^{-1}$) and lower stellar masses ($\simeq10^{8}\,M_{\odot}$), but with higher [OIII]+H$β$ equivalent widths ($\simeq$1100$Å$), and elevated production efficiency of ionizing photons ($\log(ξ_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($, 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

We present a spectroscopic survey of Ly$α$ emitters in the Extended Groth Strip (EGS) field, targeting the regime near the Epoch of Reionization. Using Keck/DEIMOS, we observed 947 high-$z$ candidates with photometric redshifts from 3 $< z_\text{phot}, 8 pages, 6 figures, 2 tables, submitted to MNRAS

We investigate the stellar mass - black hole mass ($\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$) relation with type 1 AGN down to $\mathcal{M}_\mathrm{BH}=10^7 M_\odot$, corresponding to a $\simeq -21$ absolute magnitude in rest-frame ultraviolet (UV), at $z = 2-2.5$. Exploiting the deep and large-area spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we identify 66 type 1 AGN with $\mathcal{M}_\mathrm{BH}$ ranging from $10^7$ to $10^{10} M_\odot$ that are measured with single-epoch virial method using C{\sc iv} emission lines detected in the HETDEX spectra. $\mathcal{M}_*$ of the host galaxies are estimated from optical to near-infrared photometric data taken with Spitzer, WISE, and ground-based 4-8m class telescopes by CIGALE SED fitting. We further assess the validity of SED fitting in two cases by host-nuclear decomposition performed through surface brightness profile fitting on spatially-resolved host galaxies with JWST/NIRCam CEERS data. We obtain the $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation covering the unexplored low-mass ranges of $\mathcal{M}_\mathrm{BH}~\sim~10^7-10^8~M_\odot$, and conduct forward modelling to fully account for the selection biases and observational uncertainties. The intrinsic $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation at $z\sim 2$ has a moderate positive offset of $0.52\pm0.14$~dex from the local relation, suggestive of more efficient black hole growth at higher redshift even in the low-mass regime of $\mathcal{M}_\mathrm{BH}~\sim~10^7-10^8~M_\odot$. Our $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation is inconsistent with the $\mathcal{M}_\mathrm{BH}$ suppression at the low-$\mathcal{M}_*$ regime predicted by recent hydrodynamic simulations at a $98\%$ confidence level, suggesting that feedback in the low-mass systems may be weaker than those produced in hydrodynamic simulations., Comment: 16 pages, 8 figures, accepted for publication in ApJ

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV8 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8., Comment: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of $\sim 20 \ M_\odot \textrm{yr}^{-1}$ at $z\sim4$ at a $1\sigma$ level, allowing us to quantify the amount of contamination from dusty star-forming sources in our quiescent sample. Starting with a parent sample of candidate massive quiescent galaxies from the Stevans et al. 2021 v1 SHELA catalog, we use the Bayesian \textsc{Bagpipes} spectral energy distribution fitting code to derive robust stellar masses ($M_*$) and star-formation rates (SFRs) for these sources, and select a conservative sample of 36 candidate massive ($M_* > 10^{11}M_\odot$) quiescent galaxies, with specific SFRs at $>2\sigma$ below the star-forming main sequence at $z\sim4$. Based on ALMA imaging, six of these candidate quiescent galaxies have the presence of significant dust-obscured star-formation, thus were removed from our final sample. This implies a $\sim 17\%$ contamination rate from dusty star-forming galaxies with our selection criteria using the v1 SHELA catalog. This conservatively-selected quiescent galaxy sample at $z=3-5$ will provide excellent targets for future observations to better constrain how massive galaxies can both grow and shut-down their star-formation in a relatively short time period., Comment: 16 pages, 13 figures

We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy-bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 square degree (H20) survey and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust, to the Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization.

We present analysis using a citizen science campaign to improve the cosmological measures from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the Hubble expansion rate, $H(z)$, and angular diameter distance, $D_A(z)$, at $z =$ 2.4, each to percent-level accuracy. This accuracy is determined primarily from the total number of detected Lyman-$\alpha$ emitters (LAEs), the false positive rate due to noise, and the contamination due to [O II] emitting galaxies. This paper presents the citizen science project, Dark Energy Explorers, with the goal of increasing the number of LAEs, decreasing the number of false positives due to noise and the [O II] galaxies. Initial analysis shows that citizen science is an efficient and effective tool for classification most accurately done by the human eye, especially in combination with unsupervised machine learning. Three aspects from the citizen science campaign that have the most impact are 1) identifying individual problems with detections, 2) providing a clean sample with 100% visual identification above a signal-to-noise cut, and 3) providing labels for machine learning efforts. Since the end of 2022, Dark Energy Explorers has collected over three and a half million classifications by 11,000 volunteers in over 85 different countries around the world. By incorporating the results of the Dark Energy Explorers we expect to improve the accuracy on the $D_A(z)$ and $H(z)$ parameters at $z =$ 2.4 by 10 - 30%. While the primary goal is to improve on HETDEX, Dark Energy Explorers has already proven to be a uniquely powerful tool for science advancement and increasing accessibility to science worldwide., Comment: 14 pages, 6 figures, accepted for publication in The Astrophysical Journal

A substantial number of ultra-high redshift (8 < z < 17) galaxy candidates have been detected with JWST, posing the question: are these observational results surprising in the context of current galaxy formation models? We address this question using the well-established Santa Cruz semi-analytic models, implemented within merger trees from the new suite of cosmological N-body simulations GUREFT, which were carefully designed for ultra-high redshift studies. Using our fiducial models calibrated at z=0, we present predictions for stellar mass functions, rest-frame UV luminosity functions, and various scaling relations. We find that our (dust-free) models predict galaxy number densities at z~11 (z~ 13) that are an order of magnitude (a factor of ~30) lower than the observational estimates. We estimate the uncertainty in the observed number densities due to cosmic variance, and find that it leads to a fractional error of 30-70% at z=11 (25-150% at z=13) for the available observed fields. We explore which processes in our models are most likely to be rate-limiting for the formation of luminous galaxies at these early epochs, considering the halo formation rate, gas cooling, star formation, and stellar feedback, and conclude that it is mainly efficient stellar-driven winds. We find that a modest boost of a factor of ~4 to the UV luminosities, which could arise from a top-heavy stellar initial mass function characteristic of Pop III stars, would bring our current models into agreement with the observations., 20 pages, 15 figures, submitted to MNRAS

[Abridged] We combine HST images from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey with JWST images from the Cosmic Evolution Early Release Science (CEERS) survey to measure the stellar and dust-obscured star formation distributions of a mass-complete ($>10^{10}M_\odot$) sample of 95 star-forming galaxies (SFGs) at $0.175$; 38 galaxies). At lower $S/N$, simulations show that simultaneous measurements of the size and Sérsic index become unreliable. We extend our study to fainter sources ($S/N>10$; 95 galaxies) by fixing their Sérsic index to unity. The Sérsic index of bright galaxies ($S/N>75$) has a median value of 0.7, which, together with their axis ratio distribution, suggests a disk-like morphology in the rest-MIR. Galaxies above the main sequence (MS; i.e., starbursts) have rest-MIR sizes that are a factor 2 smaller than their rest-optical sizes. The median rest-optical to rest-MIR size ratio of MS galaxies increases with stellar mass, from 1.1 at $10^{9.8}M_\odot$ to 1.6 at $10^{11}M_\odot$. This mass-dependent trend resembles the one found in the literature between the rest-optical and rest-near-infrared sizes of SFGs, suggesting that it is due to dust attenuation gradients affecting rest-optical sizes and that the sizes of the stellar and star-forming components of SFGs are, on average, consistent at all masses. There is, however, a small population of SFGs (15%) with a compact star-forming component embedded in a larger stellar structure. This could be the missing link between galaxies with an extended stellar component and those with a compact stellar component; the so-called blue nuggets., Submitted to A&A. 18 pages, 14 figures

The onset of galaxy formation is thought to be initiated by the infall of neutral, pristine gas onto the first protogalactic halos. However, direct constraints on the abundance of neutral atomic hydrogen (HI) in galaxies have been difficult to obtain at early cosmic times. Here we present spectroscopic observations with JWST of three galaxies at redshifts $z=8.8 - 11.4$, about $400-600$ Myr after the Big Bang, that show strong damped Lyman-$α$ absorption ($N_{\rm HI} > 10^{22}$ cm$^{-2}$) from HI in their local surroundings, an order of magnitude in excess of the Lyman-$α$ absorption caused by the neutral intergalactic medium at these redshifts. Consequently, these early galaxies cannot be contributing significantly to reionization, at least at their current evolutionary stages. Simulations of galaxy formation show that such massive gas reservoirs surrounding young galaxies so early in the history of the universe is a signature of galaxy formation in progress., Submitted

We present the mid-IR (MIR) morphologies for 64 star-forming galaxies at $0.210^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21~$μ$m), enabling us to measure the effective radii ($R_{\rm{eff}}$) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 $μ$m, which contains strong emission from Polycyclic aromatic hydrocarbon (PAH) features, a well-established tracer of star formation in galaxies. We define a ``PAH-band'' as the MIRI bandpass that contains these features at the redshift of the galaxy. We then compare the galaxy morphologies in the PAH-bands to those in rest-frame Near-UV (NUV) using HST ACS/F435W or ACS/F606W and optical/near-IR using HST WFC3/F160W imaging from UVCANDELS and CANDELS, where the NUV-band and F160W trace the profile of (unobscured) massive stars and the stellar continuum, respectively. The $R_{\rm{eff}}$ of galaxies in the PAH-band are slightly smaller ($\sim$10\%) than those in F160W for galaxies with $\rm{M_*\gtrsim10^{9.5}~M_\odot}$ at $z\leq1.2$, but the PAH-band and F160W have a similar fractions of light within 1 kpc. In contrast, the $R_{\rm{eff}}$ of galaxies in the NUV-band are larger, with lower fractions of light within 1 kpc compared to F160W for galaxies at $z\leq1.2$. Using the MIRI data to estimate the $\rm{SFR_{\rm{IR}}}$ surface density, we find the correlation between the $\rm{SFR_{\rm{IR}}}$ surface density and stellar mass has a steeper slope than that of the $\rm{SFR_{\rm{UV}}}$ surface density and stellar mass, suggesting more massive galaxies having increasing amounts of obscured fraction of star formation in their inner regions. This paper demonstrates how the high-angular resolution data from JWST/MIRI can reveal new information about the morphology of obscured-star formation., 28 pages, 11 figures, Accepted by ApJ

We describe the ensemble properties of the $1.9 < z < 3.5$ Lyman Alpha Emitters (LAEs) found in the HETDEX survey's first public data release, HETDEX Public Source Catalog 1 (Mentuch Cooper et al. 2023). Stacking the low-resolution ($R \sim$ 800) spectra greatly increases the signal-to-noise ratio, revealing spectral features otherwise hidden by noise, and we show that the stacked spectrum is representative of an average member of the set. The flux limited, Ly$\alpha$ signal-to-noise ratio restricted stack of 50K HETDEX LAEs shows the ensemble biweight ``average" $z \sim 2.6$ LAE to be a blue (UV continuum slope $\sim -2.4$ and E(B-V) $< 0.1$), moderately bright (M$_{\text{UV}} \sim -19.7$) star forming galaxy with strong Ly$\alpha$ emission (log $L_{Ly\alpha}$ $\sim$ 42.8 and $W_{\lambda}$(Ly$\alpha$) $\sim$ 114\AA), and potentially significant leakage of ionizing radiation. The restframe UV light is dominated by a young, metal poor stellar population with an average age 5-15 Myr and metallicity of 0.2-0.3 Z$_{\odot}$., Comment: 17 pages, 11 figures, 2 data files (ApJ Accepted)

We analyze rest-frame ultraviolet to optical spectra of three $z\simeq7.47$ - $7.75$ galaxies whose Ly$\alpha$-emission lines were previously detected with Keck/MOSFIRE observations, using the JWST/NIRSpec observations from the Cosmic Evolution Early Release Science (CEERS) survey. From NIRSpec data, we confirm the systemic redshifts of these Ly$\alpha$ emitters, and emission-line ratio diagnostics indicate these galaxies were highly ionized and metal poor. We investigate Ly$\alpha$ line properties, including the line flux, velocity offset, and spatial extension. For the one galaxy where we have both NIRSpec and MOSFIRE measurements, we find a significant offset in their flux measurements ($\sim5\times$ greater in MOSFIRE) and a marginal difference in the velocity shifts. The simplest interpretation is that the Ly$\alpha$ emission is extended and not entirely encompassed by the NIRSpec slit. The cross-dispersion profiles in NIRSpec reveal that Ly$\alpha$ in one galaxy is significantly more extended than the non-resonant emission lines. We also compute the expected sizes of ionized bubbles that can be generated by the Ly$\alpha$ sources, discussing viable scenarios for the creation of sizable ionized bubbles ($>$1 physical Mpc). The source with the highest-ionization condition is possibly capable of ionizing its own bubble, while the other two do not appear to be capable of ionizing such a large region, requiring additional sources of ionizing photons. Therefore, the fact that we detect Ly$\alpha$ from these galaxies suggests diverse scenarios on escape of Ly$\alpha$ during the epoch of reionization. High spectral resolution spectra with JWST/NIRSpec will be extremely useful for constraining the physics of patchy reionization., Comment: Submitted to ApJ (18 pages, 7 figures, 2 tables)

With the successful launch and commissioning of JWST we are now able to routinely spectroscopically probe the rest-frame optical emission of galaxies at $z>6$ for the first time. Amongst the most useful spectral diagnostics used in the optical is the Balmer/4000~\AA\ break; this is, in principle, a diagnostic of the mean ages of composite stellar populations. However, the Balmer break is also sensitive to the shape of the star formation history, the stellar (and gas) metallicity, the presence of nebular continuum emission, and dust attenuation. In this work we explore the origin of the Balmer/4000~\AA\ break using the SYNTHESIZER synthetic observations package. We then make predictions of the Balmer/4000~\AA\ break using the First Light and Reionisation Epoch Simulations (FLARES) at $5, Comment: 9 pages, 9 figures, submitted to MNRAS

The most luminous quasars at $z > 6$ are suspected to be both highly clustered and reside in the most massive dark matter halos in the early Universe, making them prime targets to search for galaxy overdensities and/or protoclusters. We search for Lyman-break dropout-selected galaxies using HST WFC3/ACS broadband imaging in the fields of three $6 < z < 7$ quasars, as well as their simultaneously observed coordinated-parallel fields, and constrain their photometric redshifts using EAZY. One field, J0305-3150, shows a volume density 10$\times$ higher than the blank-field UV luminosity function (UVLF) at M$_{UV} < -20$, with tentative evidence of a 3$σ$ overdensity in its parallel field located 15 cMpc away. Another field, J2054-0005, shows an angular overdensity within 500 ckpc from the quasar but still consistent with UVLF predictions within 3$σ$, while the last field, J2348-3054, shows no enhancement. We discuss methods for reducing uncertainty in overdensity measurements when using photometric selection and show that we can robustly select LBGs consistent with being physically associated with the quasar, corroborated by existing JWST/NIRCam WFSS data in the J0305 field. Even accounting for incompleteness, the overdensities in J0305 and J2054 are higher for brighter galaxies at short angular separations, suggesting preferential enhancement of more massive galaxies in the immediate vicinity of the quasar. Finally, we compare the LBG population with previously-identified [CII] and mm-continuum companions; the LBG overdensities are not accompanied by an enhanced number of dusty galaxies, suggesting that the overdense quasar fields are not in the bursty star-forming phase sometimes seen in high-redshift protoclusters., 22 pages (main text), 12 figures, 10 tables, 2 appendices. Final version after addressing referee report, accepted to ApJ May 2023

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.512, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products., Comment: 21 pages, 9 figures, submitted to AAS Journals

We report on the discovery of two low-luminosity, broad-line AGN at $z>5$ identified using JWST NIRSpec spectroscopy from the CEERS Survey. We detect broad H$α$ emission from both sources, with FWHM of $2038\pm286$ and $1807\pm207$ km s$^{-1}$, resulting in black hole (BH) masses that are 1-2 dex below that of existing samples of luminous quasars at $z>5$. The first source, CEERS 1670 at $z=5.242$, is 2-3 dex fainter than known quasars at similar redshifts and was previously identified as a candidate low-luminosity AGN based on its rest-frame optical SED. We measure a BH mass of $M_{\rm BH}=1.3\pm0.4\times 10^{7}~M_{\odot}$, confirming that this AGN is powered by the least-massive BH known in the universe at the end of cosmic reionization. The second source, CEERS 3210 at $z=5.624$, is inferred to be a heavily obscured, broad-line AGN caught in a transition phase between a dust-obscured starburst and an unobscured quasar. We estimate its BH mass to be $M_{\rm BH}\simeq 0.9-4.7 \times 10^{7}~M_{\odot}$, depending on the level of dust obscuration assumed. We derive host stellar masses, $M_\star$, allowing us to place constraints on the BH-galaxy mass relationship in the lowest mass range yet probed in the early universe. The $M_{\rm BH}/M_\star$ ratio for CEERS 1670, in particular, is consistent with or higher than the empirical relationship seen in massive galaxies at $z=0$. We examine the emission-line ratios of both sources and find that their location on the BPT and OHNO diagrams is consistent with model predictions for low-metallicity AGN with $Z/Z_\odot \simeq 0.2-0.4$. The spectroscopic identification of low-luminosity, broad-line AGN at $z>5$ with $M_{\rm BH}\simeq 10^{7}~M_{\odot}$ demonstrates the capability of JWST to push BH masses closer to the range predicted for the BH seed population and provides a unique opportunity to study the early stages of BH-galaxy assembly., 16 pages, 8 figures, submitted to ApJL

Active galactic nuclei (AGN) in the early Universe are thought to be prominent sources of energy and ionizing photons that affected the growth of their host galaxy and their environment. However, it is still unclear how the supermassive black holes (SMBHs) that fuel these AGN grew to the observed high masses already at high redshifts. Observations of high-redshift SMBH progenitors or lower-luminosity AGN will thus help characterize the evolution of SMBHs and their impact on the surroundings. With the launch of the JWST, fainter objects at high redshifts can now be detected, including lower-luminosity AGN. We assess the observability of such low luminosity AGN, using the cosmological simulation code GIZMO to provide a realistic environment for black hole growth in the early Universe. Soon after the first stars are born in the simulation run, we insert stellar-remnant black hole seeds of various initial masses, between $300$ and $10^4 {\rm \ M}_{\odot}$, at the center of a dark matter halo and follow their growth until $z\sim6$. Such stellar black hole seeds placed in a typical high-$z$ environment do not significantly accrete and grow to reach masses that can be observed with the JWST under conditions of standard Bondi-Hoyle accretion, as energy input from stellar feedback and chaotic dynamics prevent efficient gas accretion onto the black holes. To be observed with the JWST, rarer but still physically feasible growth regimes, involving Eddington or super-Eddington accretion, would be required. Alternatively, AGN observability may be boosted under even rarer conditions of extreme gravitational lensing., Comment: 12 pages, 10 figures, Accepted for publication in MNRAS

We use James Webb Space Telescope Near-Infrared Camera Wide Field Slitless Spectroscopy (NIRCam WFSS) and Near-Infrared spectrograph (NIRSpec) in the Cosmic Evolution Early Release survey (CEERS) to measure rest-frame optical emission-line of 155 galaxies at z>2. The blind NIRCam grism observations include a sample of galaxies with bright emission lines that were not observed on the NIRSpec masks. We study the changes of the Ha, [OIII]/Hb, and [NeIII]/[OII] emission lines in terms of redshift by comparing to lower redshift SDSS and CLEAR samples. We find a significant (>3$σ$) correlation between [OIII]/Hb with redshift, while [NeIII]/[OII] has a marginal (2$σ$) correlation with redshift. We compare [OIII]/Hb and [NeIII]/[OII] to stellar mass and Hb SFR. We find that both emission-line ratios have a correlation with Hb SFR and an anti-correlation with stellar mass across the redshifts 0, 16 pages, 11 figures

We present a robust sample of very high-redshift galaxy candidates from the first epoch of {\it JWST}/NIRCam imaging from the Next Generation Extragalactic Exploratory Deep (NGDEEP) Survey. The NGDEEP NIRCam imaging in the Hubble Ultra Deep Field Parallel Field 2 (HUDF-Par2) reaches $m=30.4$ (5$\sigma$, point-source) in F277W, making it the deepest public {\it JWST} GO imaging dataset to date. We describe our detailed data reduction process of the six-filter broad-band {\it JWST}/NIRCam imaging, incorporating custom corrections for systematic effects to produce high-quality calibrated images. Using robust photometric redshift selection criteria, we identify a sample of 38 $z \gtrsim 9$ galaxy candidates. These objects span a redshift range of $z=8.5-15.8$, and apparent magnitudes of $m_\mathrm{F277W} = 27-30.5$ AB mag, reaching $\sim 1.5$ mag deeper than previous public {\it JWST} imaging surveys. We calculate the rest-frame ultraviolet (UV) luminosity function at $z \sim$ 9 and 11, and present a new measurement of the luminosity function faint-end slope at $z \sim 11$. There is no significant evolution in the faint-end slope and number density from $z=9$ to 11. Comparing our results with theoretical predictions, we find that some models produce better agreement at the faint end than the bright end. These results will help to constrain how stellar feedback impacts star formation at these early epochs., Comment: 21 pages, 10 figures. Submitted to ApJL

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4, Updated with accepted ApJ version. Part of the CEERS Focus Issue. 27 pages, many figures (4 Figure Sets, available upon reasonable request)

The selection of high-redshift galaxies often involves spectral energy distribution (SED) fitting to photometric data, an expectation for contamination levels, and measurement of sample completeness -- all vetted through comparison to spectroscopic redshift measurements of a sub-sample. The first JWST data is now being taken over several extragalactic fields, to different depths and across various areas, which will be ideal for the discovery and classification of galaxies out to distances previously uncharted. As spectroscopic redshift measurements for sources in this epoch will not be initially available to compare with the first photometric measurements of z > 8 galaxies, robust photometric redshifts are of the utmost importance. Galaxies at z > 8 are expected to have bluer rest-frame ultraviolet (UV) colors than typically-used model SED templates, which could lead to catastrophic photometric redshift failures. We use a combination of BPASS and Cloudy models to create a supporting set of templates that match the predicted rest-UV colors of z > 8 the simulated galaxies in a mock catalog (Yung et al. 2022), which mimics expected field depths and areas of the Cosmic Evolution Early Release Science Survey (CEERS: m$_{5σ}$ ~ 28.6 over ~100 arcmin$^2$; Finkelstein et al. 2022a, Bagley et al. 2022). We use EAZY to highlight the improvements in redshift recovery with the inclusion of our new template set and suggest criteria for selecting galaxies at 8 < z < 10 with JWST, providing an important test case for observers venturing into this new era of astronomy., 14 pages, 10 figures, submitted to ApJ. SED templates published at: ceers.github.io/LarsonSEDTemplates

We perform a ground-based near-infrared spectroscopic survey using the Keck/MOSFIRE spectrograph to target Ly$α$ emission at $7.0$4$σ$ level in eight $z>7$ galaxies, which include additional members of the known $z\sim7.7$ Ly$α$-emitter (LAE) cluster (Tilvi et al. 2020). With the addition of these newly-discovered $z\sim7.7$ LAEs, this is currently the largest measured LAE cluster at $z>7$. The unusually-high Ly$α$ detection rate at $z\sim7.7$ in this field suggests significantly stronger Ly$α$ emission from the clustered LAEs than from the rest of our targets. We estimate the ionized bubble sizes around these LAEs and conclude that the LAEs are clustered within an extended ionized structure created by overlapping ionized bubbles which allow the easier escape of Ly$α$ from galaxies. It is remarkable that the brightest object in the cluster has the lowest measured redshift of the Ly$α$ line, being placed in front of the other LAEs in the line-of-sight direction. This suggests that we are witnessing the enhanced IGM transmission of Ly$α$ from galaxies on the rear side of an ionized area. This could be a consequence of Ly$α$ radiative transfer: Ly$α$ close to the central velocity is substantially scattered by the IGM while Ly$α$ from the rear-side galaxies is significantly redshifted to where it has a clear path., Submitted to ApJ (17 pages, 8 figures, 4 tables)

Studies of the high-redshift rest-frame ultraviolet (UV) luminosity functions (LFs) have typically treated the star-forming galaxy and active galactic nuclei (AGN) populations separately, as they have different survey depth and area requirements. However, the recent advent of wide-area deep ground-based imaging surveys now probe volumes large enough to discover AGNs, at depths sensitive enough for fainter star-forming galaxies, bridging these two populations. Using results from such surveys as observational constraints, we present a methodology to jointly empirically model the evolution of the rest-UV luminosity functions at z=3-9. We assume both populations have a LF well-described by a double power law modified to allow a flattening at the faint-end, and that all LF parameters evolve smoothly with redshift. This provides a good fit to the observations, and makes predictions to volume densities and redshifts not yet observed. We find that the volume density of bright (M_UV = -28) AGNs rises by five orders of magnitude from z=9 to z=3, while modestly bright (M_UV = -21) galaxies rise by only two orders of magnitude across the same epoch. The observed bright-end flattening of the z=9 LF is unlikely to be due to AGN, and rather is due to a shallowing of the bright-end slope, implying reduced feedback in bright galaxies at early times. Integrating our LFs we find that the intrinsic ionizing emissivity is dominated by galaxies at all z > 3, and this result holds even after applying a notional escape fraction. We compare our AGN LFs to predictions based on different black-hole seeding models, finding decent agreement on average, but that all models are unable to predict the observed abundance of bright AGNs. We make predictions for the upcoming Euclid and Roman observatories, showing that their respective wide-area surveys should be capable of discovering AGNs to z ~ 8., 22 pages, 10 figures, 5 tables. Submitted to ApJ. Comments welcome

We report the ALMA observation of a $z\gtrsim10$ galaxy candidate (GHZ1) discovered from the GLASS-JWST Early Release Science Program. Our ALMA program aims to detect the [OIII] emission line at the rest-frame 3393.0062 GHz ($88.36\mu$m) and far-IR continuum emission with the spectral window setup seamlessly covering a 26.125 GHz frequency range ($10.10, Comment: Accepted to ApJ after revising the figures and the analysis

We analyze the rest-frame near-UV and optical nebular spectra of three $z > 7$ galaxies from the Early Release Observations taken with the Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST). These three high-z galaxies show the detection of several strong-emission nebular lines, including the temperature-sensitive [O III] $λ$4363 line, allowing us to directly determine the nebular conditions and abundances for O/H, C/O, and Ne/O. We derive O/H abundances and ionization parameters that are generally consistent with other recent analyses. We analyze the mass-metallicity relationship (i.e., slope) and its redshift evolution by comparing between the three z > 7 galaxies and local star-forming galaxies. We also detect the C III] $λλ$1907,1909 emission in a z > 8 galaxy from which we determine the most distant C/O abundance to date. This valuable detection of log(C/O) = $-0.83\pm0.38$ provides the first test of C/O redshift evolution out to high-redshift. For neon, we use the high-ionization [Ne III] $λ$3869 line to measure the first Ne/O abundances at z>7, finding no evolution in this $α$-element ratio. We explore the tentative detection of [Fe II] and [Fe III] lines in a z>8 galaxy, which would indicate a rapid build up of metals. Importantly, we demonstrate that properly flux-calibrated and higher S/N spectra are crucial to robustly determine the abundance pattern in z>7 galaxies with NIRSpec/JWST., Accepted for publication on ApJL

Using spatially resolved H-alpha emission line maps of star-forming galaxies, we study the evolution of gradients in galaxy assembly over a wide range in redshift ($0.5, 25 pages, 12 figures, 4 tables, Submitted to ApJ

We report on the host properties of five X-ray luminous Active Galactic Nuclei (AGN) identified at $3 < z < 5$ in the first epoch of imaging from the Cosmic Evolution Early Release Science Survey (CEERS). Each galaxy has been imaged with the \textit{James Webb Space Telescope} (\jwst) Near-Infrared Camera (NIRCam), which provides spatially resolved, rest-frame optical morphologies at these redshifts. We also derive stellar masses and star formation rates for each host galaxy by fitting its spectral energy distribution using a combination of galaxy and AGN templates. The AGN hosts have an average stellar mass of ${\rm log}(M_{*}/{\rm M_{\odot}} )= 11.0$, making them among the most massive galaxies detected at this redshift range in the current CEERS pointings, even after accounting for nuclear light from the AGN. We find that three of the AGN hosts have spheroidal morphologies, one is a bulge-dominated disk and one host is dominated by point-like emission. None are found to show strong morphological disturbances that might indicate a recent interaction or merger event. Notably, all four of the resolved hosts have rest-frame optical colors consistent with a quenched or post-starburst stellar population. The presence of AGN in passively evolving galaxies at $z>3$ is significant because a rapid feedback mechanism is required in most semi-analytic models and cosmological simulations to explain the growing population of massive quiescent galaxies observed at these redshifts. Our findings are in general agreement with this picture and show that AGN can continue to inject energy into these systems after their star formation is curtailed, possibly helping to maintain their quiescent state., Submitted to ApJ; 12 Pages, 6 Figures

We analyze a sample of 25 [Ne V] $λ$3426 emission-line galaxies at $1.4, 17 pages + 5 (appendix), 7 figures + 2(appendix)

We present Ly$��$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$��$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$��$ emission lines. We derive the Ly$��$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Ly��}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}, 27 pages, 19 figures, accepted for publication in ApJ

We measure the Ly$��$ escape fraction of 935 [OIII]-emitting galaxies between $1.9 < z < 2.35$ by comparing stacked spectra from the Hubble Space Telescope/WFC3's near-IR grism to corresponding stacks from the Hobby Eberly Telescope Dark Energy Experiment's Internal Data Release 2. By measuring the stacks' H$��$ to Ly$��$ ratios, we determine the Ly$��$ escape fraction as a function of stellar mass, star formation rate, internal reddening, size, and [OIII]/H$��$ ratio. We show that the escape fraction of Ly$��$ correlates with a number of parameters, such as galaxy size, star formation rate, and nebular excitation. However, we also demonstrate that most of these relations are indirect, and the primary variables that control the escape of Ly$��$ are likely stellar mass and internal extinction. Overall, the escape of Ly$��$ declines from $\gtrsim 18\%$ in galaxies with $\log M/M_{\odot} \lesssim 9$ to $\lesssim 1\%$ for systems with $\log M/M_{\odot} \gtrsim 10$, with the sample's mean escape fraction being $6.0^{+0.6\%}_{-0.5\%}$., 21 pages, 12 figures, accepted by ApJ

Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Ly$α$) studies during the epoch of reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Ly$α$ from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at $6.06$. Furthermore, we compare the redshift evolution of $W_0$ between galaxies at different UV luminosities. The UV-bright ($M_{\text{UV}}L^{*}$) galaxies show weaker evolution with a decrease of 0.4 ($\pm$0.2) dex in $W_0$ at $z>6$ while UV-faint ($M_{\text{UV}}>-21$, or $L_{\text{UV}}6$. Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Ly$α$ transmission in the IGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity., 21 pages, 9 figures, 2 tables, 1st revision submitted to ApJ

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

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 cosmological ($30-400$ Mpc) distributions of neutral hydrogen (HI) in the inter-galactic medium (IGM) traced by Ly$��$ Emitters (LAEs) and QSOs at $z=2.1-2.5$, selected with the data of the on-going Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the eBOSS survey. We investigate spatial correlations of LAEs and QSOs with HI tomography maps reconstructed from HI Ly$��$ forest absorption in the spectra of background galaxies and QSOs obtained by the CLAMATO survey and this study, respectively. In the cosmological volume far from QSOs, we find that LAEs reside in regions of strong HI absorption, i.e. HI rich, which is consistent with results of previous galaxy-background QSO pair studies. Moreover, there is an anisotropy in the HI-distribution plot of transverse and line-of-sight distances; on average the HI absorption peak is blueshifted by $\sim 200$ km s$^{-1}$ from the LAE Ly$��$ redshift, reproducing the known average velocity offset between the Ly$��$ emission redshift and the galaxy systemic redshift. We have identified a $\sim$ 40-Mpc scale volume of HI underdensity that is a candidate for a giant HII bubble, where six QSOs and an LAE overdensity exist at $\left < z \right > =2.16$. The coincidence of the QSO and LAE overdensities with the HI underdensity indicates that the ionizing photon radiation of the QSOs has created a highly ionized volume of multiple proximity zones in a matter overdensity. Our results suggest an evolutionary picture where HI gas in an overdensity of galaxies becomes highly photoionized when QSOs emerge in the galaxies., 26 pages, 19 figures. Accepted for Publication in ApJ

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

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

Using cosmological volume simulations and a custom built sub-grid model for Pop~III star formation, we examine the baseline dust extinction in the first galaxies due to Pop~III metal enrichment in the first billion years of cosmic history. We find that while the most enriched, high-density lines of sight in primordial galaxies can experience a measurable amount of extinction from Pop~III dust ($E(B-V)_{\rm max}=0.07,\ A_{\rm V,max}\approx0.28$), the average extinction is very low with $\left< E(B-V) \right> \lesssim 10^{-3}$. We derive a power-law relationship between dark matter halo mass and extinction of $E(B-V)\propto M_{\rm halo}^{0.80}$. Performing a Monte Carlo parameter study, we establish the baseline reddening of the UV spectra of dwarf galaxies at high redshift due to Pop~III enrichment only. With this method, we find $\left-2.51\pm0.07$, which is both nearly halo mass and redshift independent., Comment: Published 01/2018 pages, 5 figures

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 present a flux-limited sample of $z\sim0.3$ Ly$\alpha$ emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism spectroscopic data. The published GALEX $z\sim0.3$ LAE sample is pre-selected from continuum-bright objects and thus is biased against high equivalent width (EW) LAEs. We remove this continuum pre-selection and compute the EW distribution and the luminosity function of the Ly$\alpha$ emission line directly from our sample. We examine the evolution of these quantities from $z\sim0.3$ to $2.2$ and find that the EW distribution shows little evidence for evolution over this redshift range. As shown by previous studies, the Ly$\alpha$ luminosity density from star-forming galaxies declines rapidly with declining redshift. However, we find that the decline in Ly$\alpha$ luminosity density from $z=2.2$ to $z=0.3$ may simply mirror the decline seen in the H$\alpha$ luminosity density from $z=2.2$ to $z=0.4$, implying little change in the volumetric Ly$\alpha$ escape fraction. Finally, we show that the observed Ly$\alpha$ luminosity density from AGNs is comparable to the observed Ly$\alpha$ luminosity density from star-forming galaxies at $z=0.3$. We suggest that this significant contribution from AGNs to the total observed Ly$\alpha$ luminosity density persists out to $z\sim2.2$., Comment: 27 pages, 19 figures, accepted for publication in ApJ

We report the results from some of the deepest Keck/Multi-Object Spectrometer For Infra-Red Exploration data yet obtained for candidate $z \gtrsim 7$ galaxies. Our data show one significant line detection with 6.5$��$ significance in our combined 10 hr of integration which is independently detected on more than one night, thus ruling out the possibility that the detection is spurious. The asymmetric line profile and non-detection in the optical bands strongly imply that the detected line is Ly$��$ emission from a galaxy at $z$(Ly$��)=7.6637 \pm 0.0011$, making it the fourth spectroscopically confirmed galaxy via Ly$��$ at $z>7.5$. This galaxy is bright in the rest-frame ultraviolet (UV; $M_{\rm UV} \sim -21.2$) with a moderately blue UV slope ($��=-2.2^{+0.3}_{-0.2}$), and exhibits a rest-frame Ly$��$ equivalent width of EW(Ly$��$) $\sim 15.6^{+5.9}_{-3.6}$ ��. The non-detection of the 11 other $z \sim$ 7-8 galaxies in our long 10 hr integration, reaching a median 5$��$ sensitivity of 28 �� in the rest-frame EW(Ly$��$), implies a 1.3$��$ deviation from the null hypothesis of a non-evolving distribution in the rest-frame EW(Ly$��$) between $36.5$, which may signal the evolving neutral fraction in the intergalactic medium at the end of the reionization epoch, although our weak evidence suggests the need for a larger statistical sample to allow for a more robust conclusion., 10 pages, 4 figures, ApJ, in press

Although the universe at redshifts greater than six represents only the first one billion years ( 6, with spectroscopically confirmed galaxies out to nearly z = 9. Using these large samples, we have begun to gain a physical insight into the processes inherent in galaxy evolution at early times. In this review, I will discuss i) the selection techniques for finding distant galaxies, including a summary of previous and ongoing ground and space-based searches, and spectroscopic followup efforts, ii) insights into galaxy evolution gleaned from measures such as the rest-frame ultraviolet luminosity function, the stellar mass function, and galaxy star-formation rates, and iii) the effect of galaxies on their surrounding environment, including the chemical enrichment of the universe, and the reionization of the intergalactic medium. Finally, I conclude with prospects for future observational study of the distant universe, using a bevy of new state-of-the-art facilities coming online over the next decade and beyond., Invited review, Accepted for publication in the Publications of the Astronomical Society of Australia. 40 pages, 10 figures, 3 tables

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10��z$ for descendants and $0.12��z$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate., 11 pages, 8 figures, 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

The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLy��/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLy��) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with "positive" Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8., Published in ApJ

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $��(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation., 31 pages, 24 figures, accepted for publication in ApJ

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies., 7 pages, 4 figures, 1 table. ApJL in press

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