Search

We investigate the morphology and physical properties of a sample of 22 IR-selected dusty star-forming galaxies at Cosmic Noon (z ~ 2), using James Webb Space Telescope NIRCam images obtained in the EGS field for the Cosmic Evolution Early Release Science survey. The exceptional resolution of the NIRCam images allows us to spatially resolve these galaxies up to 4.4um and identify their bulge/core even when very extinguished by dust. Based on red-green-blue images using the F115W, F200W and F444W filters, we divide each galaxy in several uniformly colored regions, fit their respective Spectral Energy Distribution and measure dust attenuations, stellar masses, star formation rates and ages. After classifying each region as star-forming or quiescent, we assign galaxies to three classes, depending on whether active star-formation is located in the core, in the disk or in both. (i) ~70% of our DSFGs have a compact highly dust attenuated star-forming core that can contain up to 80% of the star-formation of the galaxy but only 20-30% of its stellar mass, and is always surrounded by a larger, less attenuated massive disk (no blue nuggets); (ii) 64% (27%) of disks are significantly (strongly) lopsided, likely due to asymmetric cold gas accretion, major mergers and/or large scale instabilities; (iii) 23% of galaxies have a star-forming core embedded in a quiescent disk, they are undergoing outside-in quenching, often facilitated by their strong lopsidedness inducing small and large scale instabilities; (iv) some galaxies host highly heterogeneous disks in term of RGB colors: these are driven by in-homogeneous dust attenuation; and (v) we find surprising evidence for clump-like substructures being quiescent and/or residing in quiescent regions. This work demonstrates the major impact JWST/NIRCam has on understanding the complexity of the evolution of distant massive galaxies., 22 pages, 21 figures, submitted to A&A

With the launch of JWST and other scheduled missions aimed at probing the distant Universe, we are entering a new promising era for high-$z$ astronomy. One of our main goals is the detection of the first population of stars (Population III or Pop III stars), and models suggest that Pop III star formation is allowed well into the Epoch of Reionization (EoR), rendering this an attainable achievement. In this paper, we focus on our chance of detecting massive Pop IIIs at the moment of their death as Pair-Instability Supernovae (PISNe). We estimate the probability of discovering PISNe during the EoR in galaxies with different stellar masses ($7.5 \leq \mathrm{Log}(M_\star/\mathrm{M_\odot}) \leq 10.5$) from six dustyGadget simulations of $50h^{-1}$ cMpc per side. We further assess the expected number of PISNe in surveys with JWST/NIRCam and Roman/WFI. On average, less than one PISN is expected in all examined JWST fields, while $\simeq 1.5 \, \eta_\mathrm{III}$ PISNe may be found in a $\sim 1$ deg$^2$ Roman field, with potential for increased discoveries when considering more optimistic choices for the Pop III initial mass function and formation efficiency $\eta_\mathrm{III}$, and when including the contribution of coarsely-resolved environments. JWST/NIRCam and Roman/WFI allow the detection of massive-progenitor ($\sim 250$ $\mathrm{M_\odot}$) PISNe throughout all the optimal F200W-F356W, F277W-F444W, and F158-F213 colors. PISNe are also predominantly located at the outskirts of their hosting haloes, facilitating the disentangling of underlying stellar emission thanks to the spatial-resolution capabilities of the instruments., Comment: 15 pages, 7 figures, 1 table. Submitted to MNRAS

A population of massive galaxies invisible or very faint in deep optical/near-infrared surveys, but brighter at longer wavelengths has been uncovered in the past years. However, the nature of these optically dark/faint galaxies (OFGs, among other nomenclatures) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. Particularly, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on understanding why galaxies like OFGs are so faint at optical/near-infrared wavelengths. We find that stellar mass is the primary proxy of dust attenuation among those studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$), exhibiting 30% smaller effective radius than the average SFGs at the same stellar mass and redshift. On the contrary, they do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show stellar mass as a primary proxy of dust attenuation and compact stellar light profiles behind thick dust columns obscuring typical massive SFGs., Submitted to A&A. 13 pages, 9 figures

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

During the first 500 million years of cosmic history, the first stars and galaxies formed and seeded the cosmos with heavy elements. These early galaxies illuminated the transition from the cosmic "dark ages" to the reionization of the intergalactic medium. This transitional period has been largely inaccessible to direct observation until the recent commissioning of JWST, which has extended our observational reach into that epoch. Excitingly, the first JWST science observations uncovered a surprisingly high abundance of early star-forming galaxies. However, the distances (redshifts) of these galaxies were, by necessity, estimated from multi-band photometry. Photometric redshifts, while generally robust, can suffer from uncertainties and/or degeneracies. Spectroscopic measurements of the precise redshifts are required to validate these sources and to reliably quantify their space densities, stellar masses, and star formation rates, which provide powerful constraints on galaxy formation models and cosmology. Here we present the results of JWST follow-up spectroscopy of a small sample of galaxies suspected to be amongst the most distant yet observed. We confirm redshifts z > 10 for two galaxies, including one of the first bright JWST-discovered candidates with z = 11.4, and show that another galaxy with suggested z ~ 16 instead has z = 4.9, with strong emission lines that mimic the expected colors of more distant objects. These results reinforce the evidence for the rapid production of luminous galaxies in the very young Universe, while also highlighting the necessity of spectroscopic verification for remarkable candidates., Submitted to Nature

We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging components extending ~3.6" (30kpc) across the sky. The BGG contributes 69% of the group's total galactic stellar mass, with one of the merging components containing 76% of the total mass of the BGG and a SFR>1810Msun/yr. Most importantly, we detect intra-halo light (IHL) in several HST and JWST/NIRCam bands, allowing us to construct a state-of-the-art rest-frame UV-NIR Spectral Energy Distribution of the IHL for the first time at this high redshift. This allows stellar population characterisation of both the IHL and member galaxies, as well as the morphology distribution of group galaxies vs. their star-formation activity when coupled with Herschel data. We create a stacked image of the IHL, giving us a sensitivity to extended emission of 28.5 mag/arcsec2 at rest-frame 1um. We find that the IHL is extremely dust poor (Av~0), containing an evolved stellar population of log10(t50/yr)=8.8, corresponding to a formation epoch for 50% of the stellar material 0.63Gyr before z=1.85. There is no evidence of ongoing star-formation in the IHL. The IHL in this group at z=1.85 contributes ~10% of the total stellar mass, comparable with what is observed in local clusters. This suggests that the evolution of the IHL fraction is more self-similar with redshift than predicted by some models, challenging our understanding of IHL formation during the assembly of high-redshift clusters. JWST is unveiling a new side of group formation at this redshift, which will evolve into Virgo-like structures in the local Universe., 14 pages + appendix, 11 figures, 4 tables. Accepted to A&A on 15th May 2023

Spectroscopic studies of extreme-ionization galaxies (EIGs) are critical to our understanding of exotic systems throughout cosmic time. These EIGs exhibit spectral features requiring >54.42 eV photons: the energy needed to fully ionize helium into He2+ and emit He II recombination lines. They are likely key contributors to reionization, and they can also probe exotic stellar populations or accretion onto massive black holes. To facilitate the use of EIGs as probes of high ionization, we focus on ratios constructed from strong rest-frame UV/optical emission lines, specifically [O III] 5008, H-beta, [Ne III] 3870, [O II] 3727,3729, and [Ne V] 3427. These lines probe the relative intensity at energies of 35.12, 13.62, 40.96, 13.62 eV, and 97.12, respectively, covering a wider range of ionization than traced by other common rest-frame UV/optical techniques. We use ratios of these lines ([Ne V]/[Ne III] = Ne53 and [Ne III]/[O II]), which are closely separated in wavelength, and mitigates effects of dust attenuation and uncertainties in flux calibration. We make predictions from photoionization models constructed from Cloudy that use a broad range of stellar populations and black hole accretion models to explore the sensitivity of these line ratios to changes in the ionizing spectrum. We compare our models to observations from the Hubble Space Telescope and James Webb Space Telescope of galaxies with strong high-ionization emission lines at z ~ 0, z ~ 2, and z ~ 7. We show that the Ne53 ratio can separate galaxies with ionization from 'normal' stellar populations from those with AGN and even 'exotic' Population III models. We introduce new selection methods to identify galaxies with photoionization driven by Population III stars or intermediate-mass black hole accretion disks that could be identified in upcoming high-redshift spectroscopic surveys., 16 pages, 5 figures, 1 table. Accepted in ApJ

We present the first publicly released catalog of sources obtained from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). HETDEX is an integral field spectroscopic survey designed to measure the Hubble expansion parameter and angular diameter distance at 1.88, Comment: 38 pages, 20 figures. Data access and details about the catalog can be found online at http://hetdex.org/. A copy of the catalogs presented in this work (Version 3.2) is available to download at Zenodo doi:10.5281/zenodo.7448504

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 fraction of ionizing photons escaping from galaxies, $f_{esc}$, is at the same time a crucial parameter in modelling reionization and a very poorly known quantity, especially at high redshift. Recent observations are starting to constrain the values of $f_{esc}$ in low-z star-forming galaxies, but the validity of this comparison remains to be verified. Applying at high-z the empirical relation between $f_{esc}$ and the UV slope trends derived from the Low-z Lyman Continuum Survey, we use the DELPHI semi-analytical galaxy formation model to estimate the global ionizing emissivity of high-z galaxies, which we use to compute the resulting reionization history. We find that both the global ionizing emissivity and reionization history match the observational constraints. Assuming that the low-z correlations hold during the epoch of reionization, we find that galaxies with $-16 \lesssim M_{UV} \lesssim -13.5$ are the main drivers of reionization. We derive a population-averaged $\langle f_{esc} \rangle \simeq 8\%, 10\%, 20\%$ at z=4.5, 6, 8., 5+1 page, 3 figures, submitted to A&A

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

We present ALMA Band 7 observations of a remarkably bright galaxy candidate at $z_{\rm phot}$=$16.7^{+1.9}_{-0.3}$ ($M_{\rm UV}$=$-21.6$), S5-z17-1, identified in JWST Early Release Observation data of Stephen's Quintet. We do not detect the dust continuum at 866 $\mu$m, ruling out the possibility that S5-z17-1 is a low-$z$ dusty starburst with a star-formation rate (SFR) of $\gtrsim 30\,M_{\odot}$ yr$^{-1}$. We detect a 5.1$\sigma$ line feature at $338.726\pm0.007$ GHz exactly coinciding with the JWST source position, with a 2% likelihood of the signal being spurious. The most likely line identification would be [OIII]52$\mu$m at $z=16.01$ or [CII]158$\mu$m at $z=4.61$, whose line luminosities do not violate the non-detection of the dust continuum in both cases. Together with three other $z\gtrsim$ 11-13 candidate galaxies recently observed with ALMA, we conduct a joint ALMA and JWST spectral energy distribution (SED) analysis and find that the high-$z$ solution at $z\sim$11-17 is favored in every candidate as a very blue (UV continuum slope of $\approx-2.3$) and luminous ($M_{\rm UV}$ $\approx$ [$-$24:$-$21]) system. Still, we find in some candidates that reasonable SED fits ($\Delta$ $\chi^{2}\lesssim4$) are reproduced by type-II quasar and/or quiescent galaxy templates with strong emission lines at $z\sim3$-5, where such populations predicted from their luminosity functions and EW([OIII]+H$\beta$) distributions are abundant in survey volumes used for the $z\sim$11-17 candidates. While these recent ALMA observation results have strengthened the likelihood of the high-$z$ solutions, lower-$z$ possibilities are not completely ruled out in some of the $z\sim$11-17 candidates., Comment: 23 pages, 8 figures, 7 tables, submitted to ApJ. A typo in Table4 is corrected

MACS0647$-$JD is a triply-lensed $z\sim11$ galaxy originally discovered with the Hubble Space Telescope. Here we report new JWST imaging, which clearly resolves MACS0647$-$JD as having two components that are either merging galaxies or stellar complexes within a single galaxy. Both are very small, with stellar masses $\sim10^8\,M_\odot$ and radii $r, Comment: 27 pages, 14 figures, accepted for publication in ApJL

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is a large-volume spectroscopic survey without pre-selection of sources, searching ~ 540 deg^2 for Lyman-alpha emitting galaxies (LAEs) at 1.9 < z < 3.5. Taking advantage of such a wide-volume survey, we perform a pilot study using early HETDEX data to search for lensed Lyman-alpha emitters. After performing a proof-of-concept using a previously known lensed LAE covered by HETDEX, we perform a search for previously unknown lensed LAEs in the HETDEX spectroscopic sample. We present a catalog of 26 potential LAEs lensed by foreground, red, non-star-forming galaxies at z ~ 0.4 - 0.7. We estimate the magnification for each candidate system, finding 12 candidates to be within the strong lensing regime (magnification $μ$ > 2). Follow-up observations of these potential lensed LAEs have the potential to confirm their lensed nature and explore these distant galaxies in more detail., 28 pages, 31 figures, 2 tables, accepted for publication in ApJ

The abundance of bright galaxies at z>8 can provide key constraints on models of galaxy formation and evolution, as the predicted abundance varies greatly when different physical prescriptions for gas cooling and star formation are implemented. We present the results of a search for bright z=9-10 galaxies selected from pure-parallel Hubble Space Telescope imaging programs. We include 132 fields observed as part of the Brightest of Reionizing Galaxies survey, the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, and the WFC3 Infrared Spectroscopic Parallel survey. These observations cover a total of 620 sq. arcmin, about 70% of which is also covered with Spitzer Space Telescope infrared imaging. We identify thirteen candidate galaxies in the range 8.3, Comment: Submitted to AAS Journals; key results in Fig 17; 45 pages, 19 figures

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