Search

We use JWST/NIRISS slitless spectroscopy from the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey to investigate the physical condition of star-forming galaxies at $1.7 < z < 3.4$. At these redshifts, the deep NGDEEP NIRISS slitless spectroscopy covers the [O II]$\lambda\lambda$3726,3729, [O III]$\lambda\lambda$4959,5007, H$\beta$ and H$\alpha$ emission features for galaxies with stellar masses $\log(\mathrm{M_\ast/M_\odot}) \gtrsim 7$, nearly a factor of a hundred lower than previous studies. We focus on the [O III]/[O II] (O$_{32}$) ratio which is primarily sensitive to the ionization state and with a secondary dependence on the gas-phase metallicity of the interstellar medium. We find significant ($\gtrsim5\sigma$) correlations between the O$_{32}$ ratio and galaxy properties as O$_{32}$ increases with decreasing stellar mass, decreasing star formation rate (SFR), increasing specific SFR (sSFR$\equiv \mathrm{SFR}/M_*$), and increasing equivalent width (EW) of H$\beta$ and H$\alpha$. These trends suggest a tight connection between the ionization parameter and these galaxy properties. Galaxies at $z\sim2-3$ exhibit a higher O$_{32}$ than local normal galaxies with the same stellar masses and SFRs, indicating that they have a higher ionization parameter and lower metallicity than local normal galaxies. In addition, we observe an evolutionary trend in the O$_{32}$ -- EW(H$\beta$) relation from $z\sim0$ and $z\gtrsim5$, such that higher redshift galaxies have higher EW(H$\beta$) and higher O$_{32}$ at fixed EW. We argue that both the enhanced recent star formation activity and the higher star formation surface density may contribute to the increase in O$_{32}$ and the ionization parameter., Comment: 27 pages, 14 figures

We investigate the Lyman-$\alpha$ (Ly$\alpha$) and Lyman continuum (LyC) properties of the Sunburst Arc, a $z=2.37$ gravitationally lensed galaxy with a multiply-imaged, compact region leaking LyC and a triple-peaked Ly$\alpha$ profile indicating direct Ly$\alpha$ escape. Non-LyC-leaking regions show a redshifted Ly$\alpha$ peak, a redshifted and central Ly$\alpha$ peak, or a triple-peaked Ly$\alpha$ profile. We measure the properties of the Ly$\alpha$ profile from different regions of the galaxy using $R\sim5000$ Magellan/MagE spectra. We compare the Ly$\alpha$ spectral properties to LyC and narrowband Ly$\alpha$ maps from Hubble Space Telescope (HST) imaging to explore the subgalactic Ly$\alpha-$LyC connection. We find strong correlations (Pearson correlation coefficient $r>0.6$) between the LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) and Ly$\alpha$ (1) peak separation $v_{\rm{sep}}$, (2) ratio of the minimum flux density between the redshifted and blueshifted Ly$\alpha$ peaks to continuum flux density $f_{\rm{min}}/f_{\rm{cont}}$, and (3) equivalent width. We favor a complex \ion{H}{1} geometry to explain the Ly$\alpha$ profiles from non-LyC-leaking regions and suggest two \ion{H}{1} geometries that could diffuse and/or rescatter the central Ly$\alpha$ peak from the LyC-leaking region into our sightline across transverse distances of several hundred parsecs. Our results emphasize the complexity of Ly$\alpha$ radiative transfer and its sensitivity to the anisotropies of \ion{H}{1} gas on subgalactic scales. Large differences in the physical scales on which we observe spatially variable direct escape Ly$\alpha$, blueshifted Ly$\alpha$, and escaping LyC photons in the Sunburst Arc underscore the importance of resolving the physical scales that govern Ly$\alpha$ and LyC escape., Comment: Submitted to The Astrophysical Journal with revisions from the first referee report. Comments welcome

A multitude of JWST studies reveal a surprising over-abundance of over-massive accreting super-massive blackholes (SMBHs) -- leading to a deepening tension between theory and observation in the first billion years of cosmic time. Across X-ray to infrared wavelengths, models built off of pre-JWST predictions fail to easily reproduce observed AGN signatures (or lack thereof), driving uncertainty around the true nature of these sources. Using a sample of JWST AGN identified via their broadened Halpha emission and covered by the deepest X-ray surveys, we find neither any measurable X-ray emission nor any detection of high-ionization emission lines frequently associated with accreting SMBHs. We propose that these sources are accreting at or beyond the Eddington limit, which reduces the need for efficient production of heavy SMBH seeds at cosmic dawn. Using a theoretical model of super-Eddington accretion, we can produce the observed relative dearth of both X-ray and ultraviolet emission, as well as the high Balmer decrements, without the need for significant dust attenuation. This work indicates that super-Eddington accretion is easily achieved through-out the early Universe, and further study is required to determine what environments are required to trigger this mode of black hole growth., Comment: Submitted

We present the BoRG-JWST survey, a combination of two JWST Cycle 1 programs aimed at obtaining NIRSpec spectroscopy of representative, UV-bright $75$, highlighting the large abundance of high-redshift galaxies even in individual WFC3 pointings. The primary sample span an absolute magnitude range $-20.47$ sources currently known and comparable to the brightest sources at $z>10$. Prominent [O III]+H$\beta$ lines are found across the full sample, while a stack of sources reveals a plethora of other rest-optical lines and additional rest-UV C III]1909 \r{A} emission. Despite their luminosities, none of the low-resolution spectra display evidence for Type 1 AGN activity based on a search for broad-line emission. Lastly, we present a spectroscopic data release of 188 confirmed $0.5\lesssim z\lesssim5.0$ sources from filler MSA observations, highlighting the legacy value of the survey and a representative benchmark for comparisons to deep field observations., Comment: Submitted to ApJ, comments welcome

With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe., Comment: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

We measure the gas-phase abundances of the elements He, N, O, Ne, S, Ar, and Fe in the Lyman-continuum emitting region of the Sunburst Arc, a highly magnified galaxy at redshift $z=2.37$. We detect the temperature-sensitive auroral lines [SII]$\lambda\lambda4069,4076$, [OII]$\lambda\lambda7320,7330$, [SIII]$\lambda6312$, [OIII]$\lambda4363$, and [NeIII]$\lambda3343$ in a stacked spectrum of 5 multiple images of the Lyman-continuum emitter (LCE), from which we directly measure the electron temperature in the low, intermediate, and high ionization zones. We also detect the density-sensitive doublets of [OII]$\lambda\lambda3727,3729$, [SII]$\lambda\lambda6717,6731$, and [ArIV]$\lambda\lambda4713,4741$, which constrain the density in both the low- and high-ionization gas. With these temperature and density measurements, we measure gas-phase abundances with similar rigor as studies of local galaxies. We measure a gas-phase metallicity for the LCE of $12+\log(\textrm{O}/\textrm{H}) = 7.97 \pm 0.05$, and find an enhanced nitrogen abundance $\log(\textrm{N}/\textrm{O}) = -0.65^{+0.16}_{-0.25}$. This nitrogen abundance is consistent with enrichment from a population of Wolf-Rayet stars, additional signatures of which are reported in a companion paper. Abundances of sulfur, argon, neon, and iron are consistent with local low-metallicity HII regions and low-redshift galaxies. This study represents the most complete chemical abundance analysis of a galaxy at Cosmic Noon to date, which enables direct comparisons between local HII regions and those in the distant universe., Comment: 15 pages, 4 figures, 3 tables. Submitted to ApJ

We present JWST/NIRSpec high-resolution spectroscopy G395H/F290LP of MACS0647-JD, a gravitationally lensed galaxy merger at $z=10.167$. The new spectroscopy, which is acquired for the two lensed images (JD1 and JD2), detects and resolves emission lines in the rest-frame ultraviolet (UV) and blue optical, including the resolved [OII]3726,3729 doublet, [NeIII]3870, [HeI]3890, H$\delta$, H$\gamma$, and [OIII]4363. This is the first observation of the resolved [OII]3726,3729 doublet for a galaxy at $z>8$. We measure a line flux ratio [OII]3729/3726 $= 0.9 \pm 0.3$, which corresponds to an estimated electron density of $\log(n_{e} / \rm{cm}^{-3}) = 2.9 \pm 0.5$. This is significantly higher than the electron densities of local galaxies reported in the literature. We compile the measurements from the literature and further analyze the redshift evolution of $n_{e}$. We find that the redshift evolution follows the power-law form of $n_{e} = A\times (1+z)^{p}$ with $A=54^{+31}_{-23}$ cm$^{-3}$ and $p=1.2^{+0.4}_{-0.4}$. This power-law form may be explained by a combination of metallicity and morphological evolution of galaxies, which become, on average, more metal-poor and more compact with increasing redshift., Comment: 20 pages, 11 figures. Accepted for publication in ApJ

JWST spectroscopy has revolutionized our understanding of galaxies in the early universe. Covering wavelengths up to $5.3\,{\rm \mu m}$, NIRSpec can detect rest-frame optical emission lines H$\alpha$ out to $z = 7$ and [O III] to $z = 9.5$. Observing these lines in more distant galaxies requires longer wavelength spectroscopy with MIRI. Here we present MIRI MRS IFU observations of the lensed galaxy merger MACS0647$-$JD at $z = 10.165$. With exposure times of 4.2 hours in each of two bands, we detect H$\alpha$ at $9\sigma$, [O III]$\,\lambda5008$ at $11\sigma$, and [O III]$\,\lambda4960$ at $3\sigma$. Combined with previously reported NIRSpec spectroscopy that yields seven emission lines including the auroral line [O III]$\,\lambda4363$, we present the first direct metallicity measurement of a $z > 10$ galaxy: $12+{\rm log(O/H)}= 7.79\pm0.09$, or $0.13^{+0.02}_{-0.03}\,Z_{\odot}$. This is similar to galaxies at $z \sim 4 - 9$ with direct metallicity measurements, though higher than expected given the high specific star formation rate ${\rm log(sSFR / yr^{-1})} = -7.4 \pm 0.3$. We further constrain the ionization parameter ${\rm log}(U)$ = $-1.9 \pm 0.1$, ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})$ = $25.3\pm0.1$, and star formation rate $5.0\pm0.6\,M_{\odot}/{\rm yr}$ within the past $10\,{\rm Myr}$. These observations demonstrate the combined power of JWST NIRSpec and MIRI for studying galaxies in the first $500$ million years., Comment: 16 pages, 8 figures, accepted for publication in ApJ

Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$\lambda\lambda$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies., Comment: Accepted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 16 pages, 6 figures

Ultraviolet (UV; rest-frame $\sim1200-2000$ A) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at $5.6

We describe a custom outlier rejection algorithm for JWST/NIRSpec integral field spectroscopy. This method uses a layered sigma clipping approach that adapts clipping thresholds based upon the spatial profile of the science target. We find that this algorithm produces a robust outlier rejection while simultaneously preserving the signal of the science target. Originally developed as a response to unsatisfactory initial performance of the jwst pipeline outlier detection step, this method works either as a standalone solution, or as a supplement to the current pipeline software. Comparing leftover (i.e., not flagged) artifacts with the current pipeline's outlier detection step, we find that our method results in one fifth as many residual artifacts as the jwst pipeline. However, we find a combination of both methods removes nearly all artifacts -- an approach that takes advantage of both our algorithm's robust outlier rejection and the pipeline's use of individual dithers. This combined approach is what the TEMPLATES Early Release Science team has converged upon for our NIRSpec observations. Finally, we publicly release the code and Jupyter notebooks for the custom outlier rejection algorithm., Comment: 10 pages, including 5 figures. Submitted to the Publications of the Astronomical Society of the Pacific (PASP). Code associated with paper released at https://github.com/aibhleog/baryon-sweep

This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [O III] 5007 Angstrom and H beta; the other two targets have have JWST coverage of PAH 3.3 micron and complementary ALMA data covering the [C II] 158 micron emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has technical goals to establish best practices for the Integral Field Units (IFU) within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES's observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks., Comment: 28 pages, including 8 figures and 3 tables. Submitted to The Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program, https://iopscience.iop.org/collections/apj-231204-01-templates

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$\sigma$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$\beta$), and dust content, of galaxies at $1

We present the first results from the Web Epoch of Reionization Lyman-$\alpha$ Survey (WERLS), a spectroscopic survey of Lyman-$\alpha$ emission using Keck I/MOSFIRE and LRIS. WERLS targets bright ($J<26$) galaxy candidates with photometric redshifts of $5.5\lesssim z \lesssim 8$ selected from pre-JWST imaging embedded in the Epoch of Reionization (EoR) within three JWST deep fields: CEERS, PRIMER, and COSMOS-Web. Here, we report 11 $z\sim7-8$ Lyman-$\alpha$ emitters (LAEs; 3 secure and 8 tentative candidates) detected in the first five nights of WERLS MOSFIRE data. We estimate our observed LAE yield is $\sim13$%, broadly consistent with expectations assuming some loss from redshift uncertainty, contamination from sky OH lines, and that the Universe is approximately half-ionized at this epoch, whereby observable Lyman-$\alpha$ emission is unlikely for galaxies embedded in a neutral intergalactic medium. Our targets are selected to be UV-bright, and span a range of absolute UV magnitudes with $-23.1 < M_{\text{UV}} < -19.8$. With two LAEs detected at $z=7.68$, we also consider the possibility of an ionized bubble at this redshift. Future synergistic Keck+JWST efforts will provide a powerful tool for pinpointing beacons of reionization and mapping the large scale distribution of mass relative to the ionization state of the Universe., Comment: 27 pages, 8 figures; ApJ submitted

In this letter, we report the discovery of the highest redshift, heavily obscured, radio-loud AGN candidate selected using JWST NIRCam/MIRI, mid-IR, sub-mm, and radio imaging in the COSMOS-Web field. Using multi-frequency radio observations and mid-IR photometry, we identify a powerful, radio-loud (RL), growing supermassive black hole (SMBH) with significant spectral steepening of the radio SED ($f_{1.28 \mathrm{GHz}} \sim 2$ mJy, $q_{24\mu m} = -1.1$, $\alpha_{1.28-3\mathrm{GHz}}=-1.2$, $\Delta \alpha = -0.4$). In conjunction with ALMA, deep ground-based observations, ancillary space-based data, and the unprecedented resolution and sensitivity of JWST, we find no evidence of AGN contribution to the UV/optical/NIR data and thus infer heavy amounts of obscuration (N$_{\mathrm{H}} > 10^{23}$ cm$^{-2}$). Using the wealth of deep UV to sub-mm photometric data, we report a singular solution photo-z of $z_\mathrm{phot}$ = 7.7$^{+0.4}_{-0.3}$ and estimate an extremely massive host-galaxy ($\log M_{\star} = 11.4 -12\,\mathrm{M}_{\odot}$) hosting a powerful, growing SMBH (L$_{\mathrm{Bol}} = 4-12 \times 10^{46}$ erg s$^{-1}$). This source represents the furthest known obscured RL AGN candidate, and its level of obscuration aligns with the most representative but observationally scarce population of AGN at these epochs., Comment: Accepted to ApJL

We present JWST and ALMA results for the lensing system SPT0418-47, which includes a strongly-lensed, dusty star-forming galaxy at redshift z=4.225 and an associated multiply-imaged companion. JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star Formation (TEMPLATES). This data set provides robust, mutiwavelength detection of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C II] emission confirms that SPT0418B lies at the same redshift as SPT0418A. From a source plane reconstruction, we infer that the projected physical separation of the two galaxies is $4.42\pm 0.05$ kpc. We derive total magnifications of $\mu=29.5\pm1.2$ and $\mu=4.2\pm 0.9$ for SPT0418A and SPT0418B, respectively. We use both CIGALE and PROSPECTOR to derive stellar masses. The stellar mass ratio of SPT0418A and SPT0418B is approximately 4 to 1 ($4.5\pm 1.0$ for CIGALE and $4.2^{+1.9}_{-1.6}$ for PROSPECTOR). We also see evidence of extended structure associated with SPT0418A in the lensing reconstruction that is suggestive of a tidal feature. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels., Comment: Submitted to AAS Journals,14 pages, 7 figures

We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$\alpha$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [NII]/H$\alpha$ ratios and broad H$\alpha$ emission consistent with the presence of an AGN in a $\gtrsim$1kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing ALMA rest-frame 380$\mu$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z\sim$ 4 and beyond., Comment: 14 pages, 6 figures, accepted for publication in ApJ

The measured ages of massive, quiescent galaxies at $z\sim 3-4$ imply that massive galaxies quench as early as $z\sim 6$. While the number of spectroscopic confirmations of quiescent galaxies at $z < 3$ has increased over the years, there are only a handful at $z > 3.5$. We report spectroscopic redshifts of one secure ($z=3.757$) and two tentative ($z = 3.336$, $z=4.673$) massive ($\log(M_\ast/M_\odot) > 10.3$) quiescent galaxies with 11 hours of Keck/MOSFIRE $K$-band observations. Our candidates were selected from the FENIKS survey, which uses deep Gemini/Flamingos-2 $K_b$ $K_r$ imaging optimized for increased sensitivity to the characteristic red colors of galaxies at $z > 3$ with strong Balmer/4000 \AA\ breaks. The rest-frame $UVJ$ and $(ugi)_s$ colors of 3/4 quiescent candidates are consistent with $1-2$ Gyr old stellar populations. This places these galaxies as the oldest objects at these redshifts, and challenges the notion that quiescent galaxies at $z > 3$ are all recently-quenched, "post-starburst" galaxies. Our spectroscopy shows that the other quiescent-galaxy candidate is a broad-line AGN ($z = 3.594$) with strong, redshifted $H\beta$+[O III] emission with a velocity offset $>1000$ km/s, indicative of a powerful outflow. The star-formation history of our highest redshift candidate suggests that its progenitor was already in place by $z \sim 7-11$, reaching $\sim$ 10$^{11} M_{\odot}$ by $z \simeq 8$. These observations reveal the limit of what is possible with deep near-infrared photometry and targeted spectroscopy from the ground and demonstrate that secure spectroscopic confirmation of quiescent galaxies at $z > 4$ is only feasible with JWST., Comment: 20 pages, 11 figures, submitted to ApJ

Extreme emission line galaxies (EELGs) are considered the primary contributor to cosmic reionization and are valuable laboratories to study the astrophysics of massive stars. It is strongly expected that Roman's High Latitude Wide Area Survey (HLWAS) will find many strongly gravitationally lensed [O III] emitters at Cosmic Noon (1 < z < 2.8). Roman imaging and grism spectroscopy alone will simultaneously confirm these strong lens systems and probe their interstellar medium (ISM) and stellar properties on small scales ($\lesssim$ 100 pc). Moreover, these observations will synergize with ground-based and space-based follow-up observations of the discovered lensed [O III] emitters in multi-wavelength analyses of their properties (e.g., massive stars and possible escape of ionizing radiation), spatially resolved on the scales of individual star cluster complexes. Only Roman can uniquely sample a large number of lensed [O III] emitters to study the small scale (~ 100 pc) ISM and stellar properties of these extreme emission line galaxies, detailing the key physics of massive stars and the ISM that govern cosmic reionization., Comment: Submitted in response to the call for Roman Telescope CCS white papers

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

Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies., Comment: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

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 <5$, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam's wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies., Comment: 12 pages, 4 figures, 2 tables, submitted to ApJ

We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $\lambda\lambda$1907,1909, [OII] $\lambda$3727, [NeIII] $\lambda$3869, [NeIII] $\lambda$3968, H$\delta$ $\lambda$4101, H$\gamma$ $\lambda$4340, and [OIII] $\lambda$4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity $Z =$ log(O/H) = $7.5 - 8.0$, or $(0.06 - 0.2)$ $Z_\odot$, ionization parameter log($U$) $\sim -1.9\pm0.2$, and an ionizing photon production efficiency ${\rm log}(\xi_{\rm ion})=25.2\pm0.2\,$erg$^{-1}$ Hz. The spectrum has a softened Lyman-$\alpha$ break, evidence for a strong Ly$\alpha$ damping wing, suggesting that MACS0647-JD was unable to ionize its surroundings beyond its immediate vicinity ($R_{\text{HII}} \ll 1$ pMpc). The Ly$\alpha$ damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift $z = 10.6 \pm 0.3$. MACS0647-JD has a stellar mass log($M/M_\odot$) = $8.1 \pm 0.3$, including $\sim$ 6$\times 10^7 M_\odot$ in component A, most of which formed recently (within $\sim$ 20 Myr) with a star formation rate $2\pm1 M_\odot$ / yr, all within an effective radius $70\pm24\,$pc. The smaller component B ($r \sim 20$) pc is likely older ($\sim$100 Myr) with more dust ($A_V \sim 0.1$ mag), as found previously. Spectroscopy of a fainter companion galaxy C separated by a distance of \about\ 3$\,$kpc reveals a Lyman break consistent with $z = 10.17$. MACS0647-JD is likely the most distant galaxy merger known., Comment: 22 pages, 6 figures, accepted for publication in ApJ

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)

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_UV<-20), with blue continua (-2.38 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

During the first 500 million years of cosmic history, the first stars and galaxies formed, seeding the Universe with heavy elements and eventually reionizing the intergalactic medium. Observations with JWST have uncovered a surprisingly high abundance of candidates for early star-forming galaxies, with distances (redshifts, $z$), estimated from multi-band photometry, as large as $z\approx 16$, far beyond pre-JWST limits. While generally robust, such photometric redshifts can suffer from degeneracies and occasionally catastrophic errors. Spectroscopic measurement is required to validate these sources and to reliably quantify physical properties that can constrain galaxy formation models and cosmology. Here we present JWST spectroscopy that confirms redshifts for two very luminous galaxies with $z > 11$, but also demonstrates that another candidate with suggested $z\approx 16$ instead has $z = 4.9$, with an unusual combination of nebular line emission and dust reddening that mimics the colors expected for much more distant objects. These results reinforce evidence for the early, rapid formation of remarkably luminous galaxies, while also highlighting the necessity of spectroscopic verification. The large abundance of bright, early galaxies may indicate shortcomings in current galaxy formation models, or deviation from physical properties (such as the stellar initial mass function) that are generally believed to hold at later times., Comment: Accelerated Article Preview (AAP) version published in Nature

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ 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$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ 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 $\mu$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., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

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 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\mu$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$\beta$ lines would be redshifted beyond the wavelength range observed by NIRSpec, and the Lyman-$\alpha$ 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$\beta$ equivalent widths ($\simeq$1100$\r{A}$), and elevated production efficiency of ionizing photons ($\log(\xi_{\rm\,ion}/{\rm\,Hz\,erg}^{-1})\simeq25.8$) induced by young stellar populations ($<10$~Myrs) accounting for $\simeq20\%$ of the galaxy mass, highlighting the key contribution of faint galaxies to cosmic reionization. Taking advantage of the homogeneous selection and sensitivity, we also investigate metallicity and ISM conditions with empirical calibrations using the [OIII]/H$\beta$ ratio. We find that galaxies at $z\sim8-9$ have higher SFRs and lower metallicities than galaxies at similar stellar masses at $z\sim2-6$, which is generally consistent with the current galaxy formation and evolution models., Comment: 24 pages, 12 figures, 2 tables. ApJL Focus Issue in press

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., Comment: 16 pages, 5 figures, 1 table. Accepted in ApJ

We present results from the Cosmic Evolution Early Release Survey (CEERS) on the stellar-population parameters for 28 galaxies with redshifts $4

We perform a ground-based near-infrared spectroscopic survey using the Keck/MOSFIRE spectrograph to target Ly$\alpha$ emission at $7.0$4$\sigma$ level in eight $z>7$ galaxies, which include additional members of the known $z\sim7.7$ Ly$\alpha$-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$\alpha$ detection rate at $z\sim7.7$ in this field suggests significantly stronger Ly$\alpha$ 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$\alpha$ from galaxies. It is remarkable that the brightest object in the cluster has the lowest measured redshift of the Ly$\alpha$ 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$\alpha$ from galaxies on the rear side of an ionized area. This could be a consequence of Ly$\alpha$ radiative transfer: Ly$\alpha$ close to the central velocity is substantially scattered by the IGM while Ly$\alpha$ from the rear-side galaxies is significantly redshifted to where it has a clear path., Comment: Submitted to ApJ (17 pages, 8 figures, 4 tables)

We present the first JWST observations of the $z=4.11$ luminous radio galaxy TN J1338-1942, obtained as part of the ``Prime Extragalactic Areas for Reionization and Lensing Science'' (``PEARLS'') project. Our NIRCam observations, designed to probe the key rest-frame optical continuum and emission line features at this redshift, enable resolved spectral energy distribution modelling that incorporates both a range of stellar population assumptions and radiative shock models. With an estimated stellar mass of $\log_{10}(M/\text{M}_{\odot}) \sim 10.9$, TN J1338--1942 is confirmed to be one of the most massive galaxies known at this epoch. Our observations also reveal extremely high equivalent-width nebular emission coincident with the luminous AGN jets that is best fit by radiative shocks surrounded by extensive recent star-formation. We estimate the total star-formation rate (SFR) could be as high as $\sim1600\,\text{M}_{\odot}\,\text{yr}^{-1}$, with the SFR that we attribute to the jet induced burst conservatively $\gtrsim500\,\text{M}_{\odot}\,\text{yr}^{-1}$. The mass-weighted age of the star-formation, $t_{\text{mass}} <4$ Myr, is consistent with the likely age of the jets responsible for the triggered activity and significantly younger than that measured in the core of the host galaxy. The extreme scale of the potential jet-triggered star-formation activity indicates the potential importance of positive AGN feedback in the earliest stages of massive galaxy formation, with our observations also illustrating the extraordinary prospects for detailed studies of high-redshift galaxies with JWST., Comment: 17 pages, 15 figures, Accepted for publication in MNRAS after minor revisions

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\sigma}$ ~ 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., Comment: 14 pages, 10 figures, submitted to ApJ. SED templates published at: ceers.github.io/LarsonSEDTemplates

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

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 \targb\ is a low-$z$ dusty starburst with a star-formation rate 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 $\simeq-2.3$) and luminous ($M_{\rm UV}\simeq[-$24:$-21]$) system. Still, we find in several 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 identification of 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 several of the $z\sim$11--17 candidates, indicating the need to consider the relative surface densities of the lower-$z$ contaminants in the ultra high-$z$ galaxy search., Comment: 25 pages, 8 figures, 7 tables. ApJ in press

We present the data release and data reduction process for the Epoch 1 NIRCam observations for the Cosmic Evolution Early Release Science Survey (CEERS). These data consist of NIRCam imaging in six broadband filters (F115W, F150W, F200W, F277W, F356W and F444W) and one medium band filter (F410M) over four pointings, obtained in parallel with primary CEERS MIRI observations (Yang et al. in prep). We reduced the NIRCam imaging with the JWST Calibration Pipeline, with custom modifications and reduction steps designed to address additional features and challenges with the data. Here we provide a detailed description of each step in our reduction and a discussion of future expected improvements. Our reduction process includes corrections for known pre-launch issues such as 1/f noise, as well as in-flight issues including snowballs, wisps, and astrometric alignment. Many of our custom reduction processes were first developed with pre-launch simulated NIRCam imaging over the full 10 CEERS NIRCam pointings. We present a description of the creation and reduction of this simulated dataset in the Appendix. We provide mosaics of the real images in a public release, as well as our reduction scripts with detailed explanations to allow users to reproduce our final data products. These represent one of the first official public datasets released from the Directors Discretionary Early Release Science (DD-ERS) program., Comment: 27 pages, 14 figures, submitted to ApJ. Accompanying CEERS public Data Release 0.5 available at ceers.github.io/releases.html

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<100\,\rm pc$. The brighter larger component "A" is intrinsically very blue ($\beta\sim-2.6$), likely due to very recent star formation and no dust, and is spatially extended with an effective radius $\sim70\,\rm pc$. The smaller component "B" appears redder ($\beta\sim-2$), likely because it is older ($100-200\,\rm Myr$) with mild dust extinction ($A_V\sim0.1\,\rm mag$), and a smaller radius $\sim20\,\rm pc$. We identify galaxies with similar colors in a high-redshift simulation, finding their star formation histories to be out of phase. With an estimated stellar mass ratio of roughly 2:1 and physical projected separation $\sim400\,\rm pc$, we may be witnessing a galaxy merger 400 million years after the Big Bang. We also identify a candidate companion galaxy C $\sim3\,{\rm kpc}$ away, likely destined to merge with galaxies A and B. The combined light from galaxies A+B is magnified by factors of $\sim$8, 5, and 2 in three lensed images JD1, 2, and 3 with F356W fluxes $\sim322$, $203$, $86\,\rm nJy$ (AB mag 25.1, 25.6, 26.6). MACS0647$-$JD is significantly brighter than other galaxies recently discovered at similar redshifts with JWST. Without magnification, it would have AB mag 27.3 ($M_{UV}=-20.4$). With a high confidence level, we obtain a photometric redshift of $z=10.6\pm0.3$ based on photometry measured in 6 NIRCam filters spanning $1-5\rm\mu m$, out to $4300\,\r{A}$ rest-frame. JWST NIRSpec observations planned for January 2023 will deliver a spectroscopic redshift and a more detailed study of the physical properties of MACS0647$-$JD., Comment: 27 pages, 14 figures, accepted for publication in ApJL

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0., Comment: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

We report the discovery of four galaxy candidates observed 450 - 600 Myr after the Big Bang with photometric redshifts between z ~ 8.3 - 10.2 measured using the JWST NIRCam imaging of the galaxy cluster WHL0137-08 observed in 8 filters spanning 0.8-5.0 micron, plus 9 Hubble filters spanning 0.4-1.7 micron. One candidate is gravitationally lensed with a magnification of ~8, while the other three are located in a nearby NIRCam module with expected magnifications of <~ 1.1. Using SED fitting, we estimate the stellar masses of these galaxies are typically in the range log M*/Msun = 8.3 - 8.7. All appear young with mass-weighted ages < 240 Myr, low dust content A_V < 0.15 mag, and specific star formation rates sSFR ~0.25 - 10 /Gyr for most. One z ~ 9 candidate is consistent with an age < 5 Myr and a sSFR ~10 /Gyr, as inferred from a strong F444W excess, implying [OIII]+H-beta rest-frame equivalent width ~2000 Angstrom, although an older z ~ 10 object is also allowed. Another z ~ 9 candidate is lensed into an arc 2.4 arcsec long with a magnification of ~8. This arc is the most spatially-resolved galaxy at z ~ 9 known to date, revealing structures ~30 pc across. Follow-up spectroscopy of WHL0137-08 with JWST/NIRSpec will be useful to spectroscopically confirm these high-redshift galaxy candidates and to study their physical properties in more detail., Comment: 20 pages, 5 figures, accepted for publication in ApJ

We analyze a sample of 25 [Ne V] $\lambda$3426 emission-line galaxies at $1.4

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., Comment: Submitted to ApJ; 12 Pages, 6 Figures

The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$\mu$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $\mu > 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe., Comment: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We present rest-frame optical emission-line flux ratio measurements for five $z>5$ galaxies observed by the JWST Near-Infared Spectrograph (NIRSpec) in the SMACS 0723 Early Release Observations. We add several quality-control and post-processing steps to the NIRSpec pipeline reduction products in order to ensure reliable relative flux calibration of emission lines that are closely separated in wavelength, despite the uncertain \textit{absolute} spectrophotometry of the current version of the reductions. Compared to $z\sim3$ galaxies in the literature, the $z>5$ galaxies have similar [OIII]$\lambda$5008/H$\beta$ ratios, similar [OIII]$\lambda$4364/H$\gamma$ ratios, and higher ($\sim$0.5 dex) [NeIII]$\lambda$3870/[OII]$\lambda$3728 ratios. We compare the observations to MAPPINGS V photoionization models and find that the measured [NeIII]$\lambda$3870/[OII]$\lambda$3728, [OIII]$\lambda$4364/H$\gamma$, and [OIII]$\lambda$5008/H$\beta$ emission-line ratios are consistent with an interstellar medium that has very high ionization ($\log(Q) \simeq 8-9$, units of cm~s$^{-1}$), low metallicity ($Z/Z_\odot \lesssim 0.2$), and very high pressure ($\log(P/k) \simeq 8-9$, units of cm$^{-3}$). The combination of [OIII]$\lambda$4364/H$\gamma$ and [OIII]$\lambda$(4960+5008)/H$\beta$ line ratios indicate very high electron temperatures of $4.1<\log(T_e/{\rm K})<4.4$, further implying metallicities of $Z/Z_\odot \lesssim 0.2$ with the application of low-redshift calibrations for ``$T_e$-based'' metallicities. These observations represent a tantalizing new view of the physical conditions of the interstellar medium in galaxies at cosmic dawn., Comment: Accepted for publication in AAS Journals. 14 pages, 6 figures, 3 tables

We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10m Keck telescope to search for Ly$\alpha$ emission from candidate galaxies at z ~9-10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Ly$\alpha$ from one object in ~8.1 hours of integration, at z = 8.665 +/- 0.001 with an integrated signal-to-noise > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field, and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Ly$\alpha$ detected at z=8.683 (Zitrin et al. 2015). The detection of Ly$\alpha$ suggests the existence of large (~1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV Luminosity Function down to M_UV = -13). The Ly$\alpha$ we detect would be able to escape the predominantly-neutral intergalactic medium at this epoch if our galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. 2021. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 which is powering one (or more) ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field., Comment: 19 pages, 12 figures, 8 tables, Published in ApJ