Search

IRAC data have long implied that early ($z\gtrsim6$) galaxies often have very high specific star formation rates (sSFR$\gtrsim$30 Gyr$^{-1}$), but JWST data have shown that at least some early galaxies are forming stars far less vigorously. Here, we systematically analyze the recent star formation histories (SFHs) of a large ($N=368$) sample of $z\sim6$ Lyman-break galaxies (LBGs) spanning $-22\lesssim M_\mathrm{UV}\lesssim-16$ assembled from ACS+NIRCam imaging in the GOODS and Abell 2744 fields. We find that very low H$\alpha$-to-UV luminosity ratios ($L_\mathrm{H\alpha}/L_\mathrm{UV}$) and strong recent downturns in star formation rate (SFR) are $\approx$5$\times$ more common among the UV-faintest subset of our sample ($\langle M_\mathrm{UV}\rangle=-17.4$) compared to the brightest subset ($\langle M_\mathrm{UV}\rangle=-20.0$). The frequency of high $L_\mathrm{H\alpha}/L_\mathrm{UV}$ and strong recent SFR upturns is approximately constant with UV luminosity. We discuss how bursty SFHs naturally reproduce this much greater diversity in recent SFHs among very UV-faint galaxies. Using public NIRSpec/prism data, we newly confirm recent strong SFR downturns among three LBGs in our sample, and validate our photometric inferences on key SFH signatures among $z\sim6$ LBGs in general. Our results imply that early galaxies frequently cycle through phases of rapid stellar mass assembly and other periods of much slower growth. This yields huge ($\gtrsim$1-2 mag) fluctuations in $M_\mathrm{UV}$ on rapid ($\sim$10-30 Myr) timescales, helping explain the surprising abundance of $z>10$ galaxies. Finally, we caution that this burstiness causes all existing high-redshift samples (particularly line-selected samples) to be far less complete to galaxies with long recent phases of low sSFR than those currently undergoing a burst., Comment: Submitted to ApJ. 29 pages, 13 figures, 1 appendix

The James Webb Space Telescope (JWST) has spectroscopically confirmed numerous galaxies at $z > 10$. While weak rest-ultraviolet emission lines have only been seen in a handful of sources, the stronger rest-optical emission lines are highly diagnostic and accessible at mid-infrared wavelengths with the Mid-Infrared Instrument (MIRI) of JWST. We report the photometric detection of the most distant spectroscopically confirmed galaxy JADES-GS-z14-0 at $z = 14.32^{+0.08}_{-0.20}$ with MIRI at $7.7\ \mu\mathrm{m}$. The most plausible solution for the stellar population properties is that this galaxy contains half a billion solar masses in stars with a strong burst of star formation in the most recent few million years. For this model, at least one-third of the flux at $7.7\ \mu\mathrm{m}$ comes from the rest-optical emission lines $\mathrm{H}\beta$ and/or $\mathrm{[OIII]}\lambda\lambda4959,5007$. The inferred properties of JADES-GS-z14-0 suggest rapid mass assembly and metal enrichment during the earliest phases of galaxy formation., Comment: Submitted; main text has 9 pages, 3 figures and 1 table; extended text has 15 pages, 5 figures, and 1 table

JWST has recently sparked a new era of Ly$\alpha$ spectroscopy, delivering the first measurements of the Ly$\alpha$ escape fraction and velocity profile in typical galaxies at $z\simeq6-10$. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Ly$\alpha$ properties in galaxies at $z\simeq5-6$ when the IGM is mostly ionized. Here we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Ly$\alpha$ velocity offsets, escape fractions, and EWs in $z\simeq5-6$ galaxies. We find that galaxies with large Ly$\alpha$ escape fractions (> 0.2) are common at $z\simeq5-6$, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Ly$\alpha$ becomes more prevalent in the galaxy population toward higher redshift from $z\sim3$ to $z\sim6$, although we find that this evolution slows considerably between $z\sim5$ and $z\sim6$, consistent with modest attenuation from residual HI in the mostly ionized IGM at $z\simeq5-6$. We find significant evolution in Ly$\alpha$ velocity profiles between $z\simeq2-3$ and $z\simeq5-6$, likely reflecting the influence of resonant scattering from residual intergalactic HI on the escape of Ly$\alpha$ emission near line center. This effect will make it challenging to use Ly$\alpha$ peak offsets as a probe of Lyman continuum leakage at $z\simeq5-6$. We use our $z\simeq5-6$ Ly$\alpha$ distributions to make predictions for typical Ly$\alpha$ properties at $z\gtrsim8$ and discuss implications of a recently-discovered Ly$\alpha$ emitter at $z\simeq8.5$ with a small peak velocity offset (156 km s$^{-1}$)., Comment: 31 pages, 12 figures, accepted by MNRAS

We characterize the earliest galaxy population in the JADES Origins Field (JOF), the deepest imaging field observed with JWST. We make use of the ancillary Hubble optical images (5 filters spanning $0.4-0.9\mu\mathrm{m}$) and novel JWST images with 14 filters spanning $0.8-5\mu\mathrm{m}$, including 7 medium-band filters, and reaching total exposure times of up to 46 hours per filter. We combine all our data at $>2.3\mu\mathrm{m}$ to construct an ultradeep image, reaching as deep as $\approx31.4$ AB mag in the stack and 30.3-31.0 AB mag ($5\sigma$, $r=0.1"$ circular aperture) in individual filters. We measure photometric redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts $z=11.5-15$. These objects show compact half-light radii of $R_{1/2}\sim50-200$pc, stellar masses of $M_{\star}\sim10^7-10^8 M_{\odot}$, and star-formation rates of $\mathrm{SFR}\sim0.1-1\,M_{\odot}\,\mathrm{yr}^{-1}$. Our search finds no candidates at $15

We describe {\it JWST}/NIRSpec prism measurements of Ly$\alpha$ emission in $z\gtrsim 5$ galaxies. We identify Ly$\alpha$ detections in 10 out of 69 galaxies with robust rest-optical emission line redshift measurements at $5\leq z<7$ in the CEERS and DDT-2750 observations of the EGS field. Galaxies at $z\simeq 6$ with faint continuum (F150W $=$ 27--29 mag) are found with extremely large rest-frame Ly$\alpha$ equivalent widths (ranging up to 286 A). Likely Ly$\alpha$ detections are also seen in two new $z>7$ galaxies ($z=$ 7.49 and 7.17) from the second epoch of CEERS observations, both showing large Ly$\alpha$ equivalent widths that likely indicate significant transmission through the IGM. We measure high Ly$\alpha$ escape fractions in the 12 Ly$\alpha$ emitters in our sample (median 0.28), two of which show $f_{\rm esc}^{ {\rm Ly}\alpha}$ near unity ($>0.80$). We find that $50_{-11}^{+11}$% of $z\simeq 6$ galaxies with [OIII]+H$\beta$ EW $>$ 1000 A have $f_{\rm esc}^{ {\rm Ly}\alpha}$ $>0.2$, consistent with the fractions found in lower-redshift samples with matched [OIII]+H$\beta$ EWs. While uncertainties are still significant, we find that only $10_{-5}^{+9}$% of $z>7$ galaxies with similarly strong rest optical emission lines show such large $f_{\rm esc}^{ {\rm Ly}\alpha}$, as may be expected if IGM attenuation of Ly$\alpha$ increases towards higher redshifts. We identify photometric galaxy overdensities near the $z\gtrsim 7$ Ly$\alpha$ emitters, potentially providing the ionizing flux necessary to create large ionized sightlines that facilitate Ly$\alpha$ transmission. Finally, we investigate the absence of Ly$\alpha$ emission in a comparable (and spectroscopically confirmed) galaxy overdensity at $z=7.88$ in the Abell 2744 field, discussing new prism spectra of the field obtained with the UNCOVER program., Comment: main text 20 pages, 13 figures; accepted for publication by MNRAS

The detection of numerous and relatively bright galaxies at redshifts z > 9 has prompted new investigations into the star-forming properties of high-redshift galaxies. Using local forms of the initial mass function (IMF) to estimate stellar masses of these galaxies from their light output leads to galaxy masses that are at the limit allowed for the state of the LambdaCDM Universe at their redshift. We explore how varying the IMF assumed in studies of galaxies in the early universe changes the inferred values for the stellar masses of these galaxies. We infer galaxy properties with the SED fitting code Prospector using varying IMF parameterizations for a sample of 102 galaxies from the JWST Advanced Deep Extragalactic Survey (JADES) spectroscopically confirmed to be at z > 6.7, with additional photometry from the JWST Extragalactic Medium Band Survey (JEMS) for twenty-one galaxies. We demonstrate that models with stellar masses reduced by a factor of three or more do not affect the modeled spectral energy distribution (SED)., Comment: The Significance statement is required for PNAS submission

The power-law slope of the rest-UV continuum ($f_{\lambda}\propto\lambda^{\beta}$) is a key metric of early star forming galaxies, providing one of our only windows into the stellar populations and physical conditions of $z>10$ galaxies. Expanding upon previous studies with limited sample sizes, we leverage deep imaging from JADES to investigate the UV slopes of 179 $z>9$ galaxies with apparent magnitudes of $m_{\rm F200W}=26-31$, which display a median UV slope of $\beta=-2.4$. We compare to a statistical sample of $z=5-9$ galaxies, finding a shift toward bluer rest-UV colors at all $\rm~M_{UV}$. The most UV-luminous $z>9$ galaxies are significantly bluer than their lower-redshift counterparts, representing a dearth of moderately-red galaxies in the first $500~$Myr. At yet earlier times, the $z>11$ galaxy population exhibits very blue UV slopes, implying very low attenuation from dust. We identify a robust sample of 44 galaxies with $\beta<-2.8$, which have SEDs requiring models of density-bounded HII regions and median ionizing photon escape fractions of $0.51$ to reproduce. Their rest-optical colors imply that this sample has weaker emission lines (median $m_{\rm F356W}-m_{\rm F444W}=0.19$ mag) than typical galaxies (median $m_{\rm F356W}-m_{\rm F444W}=0.39$ mag), consistent with the inferred escape fractions. This sample has relatively low stellar masses (median $\log(M/M_{\odot})=7.5$), and specific star-formation rates (median$=79\rm/Gyr$) nearly twice that of our full sample (median$=44\rm/Gyr$), suggesting they are more common among systems experiencing a recent upturn in star formation. We demonstrate that the shutoff of star formation provides an alternative solution for modelling of extremely blue UV colors, making distinct predictions for the rest-optical emission of these galaxies. Future spectroscopy will be required to distinguish between these physical pictures., Comment: 17 pages, 13 figures; submitted to MNRAS

Reionization is thought to be driven by faint star-forming galaxies, but characterizing this population has long remained very challenging. Here we utilize deep nine-band NIRCam imaging from JADES to study the star-forming and ionizing properties of 756 $z\sim6-9$ galaxies, including hundreds of very UV-faint objects ($M_\mathrm{UV}>-18$). The faintest ($m\sim30$) galaxies in our sample typically have stellar masses of $M_\ast\sim(1-3)\times10^7$ $M_\odot$ and young light-weighted ages ($\sim$50 Myr), though some show strong Balmer breaks implying much older ages ($\sim$500 Myr). We find no evidence for extremely massive galaxies ($>3\times10^{10}$ $M_\odot$). We infer a strong (factor $>$2) decline in the typical [OIII]$+$H$\beta$ EWs towards very faint $z\sim6-9$ galaxies, yet a weak UV luminosity dependence on the H$\alpha$ EWs at $z\sim6$. We demonstrate that these EW trends can be explained if fainter galaxies have systematically lower metallicities as well as more recently-declining star formation histories relative to the most UV-luminous galaxies in our sample. Our data provide evidence that the brightest galaxies are frequently experiencing a recent strong upturn in SFR. We also discuss how the EW trends may be influenced by a strong correlation between $M_\mathrm{UV}$ and Lyman continuum escape fraction. This alternative explanation has dramatically different implications for the contribution of galaxies along the luminosity function to cosmic reionization. Finally, we quantify the photometric overdensities around two $z>7$ strong Ly$\alpha$ emitters. One Ly$\alpha$ emitter lies close to a strong photometric overdensity while the other shows no significant nearby overdensity, perhaps implying that not all strong $z>7$ Ly$\alpha$ emitters reside in large ionized bubbles., Comment: Accepted for publication in MNRAS. Minor changes relative to previous version. Catalog of galaxy coordinates, inferred properties, and photometry available here: https://www.ryan-endsley.com/z6to9-jades-galaxy-catalog

We present a catalog of 717 candidate galaxies at $z > 8$ selected from 125 square arcminutes of NIRCam imaging as part of the JWST Advanced Deep Extragalactic Survey (JADES). We combine the full JADES imaging dataset with data from the JEMS and FRESCO JWST surveys along with extremely deep existing observations from HST/ACS for a final filter set that includes fifteen JWST/NIRCam filters and five HST/ACS filters. The high-redshift galaxy candidates were selected from their estimated photometric redshifts calculated using a template fitting approach, followed by visual inspection from seven independent reviewers. We explore these candidates in detail, highlighting interesting resolved or extended sources, sources with very red long-wavelength slopes, and our highest redshift candidates, which extend to $z_{phot} = 18$. We also investigate potential contamination by stellar objects, and do not find strong evidence from SED fitting that these faint high-redshift galaxy candidates are low-mass stars. Over 93\% of the sources are newly identified from our deep JADES imaging, including 31 new galaxy candidates at $z_{phot} > 12$. Using 42 sources in our sample with measured spectroscopic redshifts from NIRSpec and FRESCO, we find excellent agreement to our photometric redshift estimates, with no catastrophic outliers and an average difference of $\langle \Delta z = z_{phot}- z_{spec} \rangle= 0.26$. These sources comprise one of the most robust samples for probing the early buildup of galaxies within the first few hundred million years of the Universe's history., Comment: v2: 43 pages, 20 figures, accepted by The Astrophysical Journal, full online data catalog found at https://doi.org/10.5281/zenodo.7996499

JWST has revolutionized the field of extragalactic astronomy with its sensitive and high-resolution infrared view of the distant universe. Adding to the new legacy of JWST observations, we present the first NIRCam imaging data release from the JWST Advanced Deep Extragalactic Survey (JADES) providing 9 filters of infrared imaging of $\sim$25 arcmin$^2$ covering the Hubble Ultra Deep Field and portions of Great Observatories Origins Deep Survey (GOODS) South. Utilizing 87 on-sky dual-filter hours of exposure time, these images reveal the deepest ever near-infrared view of this iconic field. We supply carefully constructed 9-band mosaics of the JADES bands, as well as matching reductions of 5 additional bands from the JWST Extragalactic Medium-band Survey (JEMS). Combining with existing HST imaging, we provide 23-band space-based photometric catalogs and photometric redshifts for $\approx47,500$ sources. To promote broad engagement with the JADES survey, we have created an interactive {\tt FitsMap} website to provide an interface for professional researchers and the public to experience these JWST datasets. Combined with the first JADES NIRSpec data release, these public JADES imaging and spectroscopic datasets provide a new foundation for discoveries of the infrared universe by the worldwide scientific community., Comment: Several figures were modified to use better line styles. A brief comparison to IRAC Channel 1 photometry was added along with a few other clarifications. Paper has been accepted for publication in ApJS

We present an overview of the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), an ambitious program of infrared imaging and spectroscopy in the GOODS-S and GOODS-N deep fields, designed to study galaxy evolution from high redshift to cosmic noon. JADES uses about 770 hours of Cycle 1 guaranteed time largely from the Near-Infrared Camera (NIRCam) and Near-Infrared Spectrograph (NIRSpec) instrument teams. In GOODS-S, in and around the Hubble Ultra Deep Field and Chandra Deep Field South, JADES produces a deep imaging region of ~45 arcmin$^2$ with an average of 130 hrs of exposure time spread over 9 NIRCam filters. This is extended at medium depth in GOODS-S and GOODS-N with NIRCam imaging of ~175 arcmin$^2$ with an average exposure time of 20 hrs spread over 8-10 filters. In both fields, we conduct extensive NIRSpec multi-object spectroscopy, including 2 deep pointings of 55 hrs exposure time, 14 medium pointings of ~12 hrs, and 15 shallower pointings of ~4 hrs, targeting over 5000 HST and JWST-detected faint sources with 5 low, medium, and high-resolution dispersers covering 0.6-5.3 microns. Finally, JADES extends redward via coordinated parallels with the JWST Mid-Infrared Instrument (MIRI), featuring ~9 arcmin$^2$ with 43 hours of exposure at 7.7 microns and twice that area with 2-6.5 hours of exposure at 12.8 microns For nearly 30 years, the GOODS-S and GOODS-N fields have been developed as the premier deep fields on the sky; JADES is now providing a compelling start on the JWST legacy in these fields., Comment: 33 pages, submitted to ApJ Supplement. The JADES Collaboration web site is at https://jades-survey.github.io, and the initial data release is available at https://archive.stsci.edu/hlsp/jades with a viewer at http://jades.idies.jhu.edu

Several studies have detected Lyman-alpha (Ly$\alpha$) from bright ($M_\mathrm{UV}\lesssim-21.5$) galaxies during the early stages of reionization despite the significantly neutral intergalactic medium. To explain these detections, it has been suggested that $z>7$ Ly$\alpha$ emitters (LAEs) inhabit physical Mpc (pMpc)-scale ionized regions powered by overdensities of faint galaxies, but systematic searches for these overdensities near LAEs have been challenging. Here, we use CEERS JWST/NIRCam imaging to search for large-scale galaxy overdensities near two very UV-bright, $z=8.7$ LAEs in the EGS field. We colour select 27 $z=8.4-9.1$ candidates, including the one LAE in the footprint (EGSY8p7). From SED models, we infer moderately faint UV luminosities ($-21.2\lesssim{M_\mathrm{UV}}\lesssim -19.1$) and stellar masses of $M_*\approx10^{7.5-8.8}$ M$_\odot$. All are efficient ionizing agents ($\xi_{\mathrm{ion}}^{*}\approx10^{25.5-26.0}$ Hz erg$^{-1}$) and are generally morphologically simple with only one compact ($r_e\lesssim140$ to $\sim650$ pc) star-forming component. 13 candidates lie within 5 arcmin of EGSY8p7, leading to a factor-of-four galaxy overdensity at $\lesssim 5$ arcmin ($\sim 1.4$ projected pMpc at $z\sim8.7$) separations from EGSY8p7. Separations of $10-15$ arcmin ($\sim2.7-4.1$ projected pMpc) are consistent with an average field. The spatial distribution of our sample may qualitatively suggest an $R\geq2$ pMpc ionized bubble encompassing both LAEs in EGS, which is theoretically unexpected but may be possible for a galaxy population $4\times$ more numerous than the average to create with moderate escape fractions ($f_\mathrm{esc}\gtrsim0.15$) over long times ($\gtrsim200$ Myr). Upcoming spectroscopic follow-up will characterize the size of any ionized bubble that may exist and the properties of the galaxies powering such a bubble., Comment: 20 pages, 7 main figures, 3 appendices, accepted to MNRAS

We present JADES JWST/NIRSpec spectroscopy of GN-z11, the most luminous candidate $z>10$ Lyman break galaxy in the GOODS-North field with $M_{UV}=-21.5$. We derive a redshift of $z=10.603$ (lower than previous determinations) based on multiple emission lines in our low and medium resolution spectra over $0.8-5.3 \mu$m. We significantly detect the continuum and measure a blue rest-UV spectral slope of $\beta=-2.4$. Remarkably, we see spatially-extended Lyman-$\alpha$ in emission (despite the highly-neutral IGM expected at this early epoch), offset 555 km s$^{-1}$ redward of the systemic redshift. From our measurements of collisionally-excited lines of both low- and high-ionization (including [O II]$\lambda3727$, [Ne III]$\lambda 3869$ and C III]$\lambda1909$) we infer a high ionization parameter ($\log U\sim -2$). We detect the rarely-seen N IV]$\lambda1486$ and N III]$\lambda1748$ lines in both our low and medium resolution spectra, with other high ionization lines seen in the low resolution spectrum such as He II (blended with O III]) and C IV (with a possible P-Cygni profile). Based on the observed rest-UV line ratios, we cannot conclusively rule out photoionization from AGN, although the high C III]/He II and N III]/He II ratios are compatible with a star-formation explanation. If the observed emission lines are powered by star formation, then the strong N III]$\lambda1748$ observed may imply an unusually high $N/O$ abundance. Balmer emission lines (H$\gamma$, H$\delta$) are also detected, and if powered by star formation rather than an AGN we infer a star formation rate of $\sim 20-30 M_{\odot} yr^{-1}$ (depending on the IMF) and low dust attenuation. Our NIRSpec spectroscopy confirms that GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr after the Big Bang., Comment: Accepted for publication in Astronomy & Astrophysics, 17 pages, 12 figures. Line fluxes updated

We present JWST NIRCam 9-band near-infrared imaging of the luminous $z=10.6$ galaxy GN-z11 from the JWST Advanced Deep Extragalactic Survey (JADES) of the GOODS-N field. We find a spectral energy distribution (SED) entirely consistent with the expected form of a high-redshift galaxy: a clear blue continuum from 1.5 to 4 microns with a complete dropout in F115W. The core of GN-z11 is extremely compact in JWST imaging. We analyze the image with a two-component model, using a point source and a S\'{e}rsic profile that fits to a half-light radius of 200 pc and an index $n=0.9$. We find a low-surface brightness haze about $0.4''$ to the northeast of the galaxy, which is most likely a foreground object but might be a more extended component of GN-z11. At a spectroscopic redshift of 10.60 (Bunker et al. 2023), the comparison of the NIRCam F410M and F444W images spans the Balmer jump. From population synthesis modeling, here assuming no light from an active galactic nucleus, we reproduce the SED of GN-z11, finding a stellar mass of $\sim$$10^{9}~M_{\odot}$, a star-formation rate of $\sim$$20~M_{\odot}~\mathrm{yr}^{-1}$ and a young stellar age of $\sim$$20~\mathrm{Myr}$. As massive galaxies at high redshift are likely to be highly clustered, we search for faint neighbors of GN-z11, finding 9 galaxies out to $\sim$5 comoving Mpc transverse with photometric redshifts consistent with $z=10.6$, and a 10$^{\rm th}$ more tentative dropout only $3''$ away. This is consistent with GN-z11 being hosted by a massive dark-matter halo ($\approx8\times10^{10}~M_{\odot}$), though lower halo masses cannot be ruled out., Comment: 18 pages, 8 figures, accepted for publication in ApJ

We describe new JWST/NIRSpec observations of galaxies at $z\gtrsim7$ taken from the CEERS survey. Previous observations of this area have revealed associations of Ly$\alpha$ emitters at redshifts ($z=7.5$, $7.7$, $8.7$) where the intergalactic medium (IGM) is thought to be mostly neutral, leading to suggestions that these systems are situated in large ionized bubbles. We identify 21 $z\gtrsim7$ galaxies with robust redshifts in the CEERS dataset, including 10 in the Ly$\alpha$ associations. Their spectra are indicative of very highly ionized and metal poor gas, with line ratios (O32 $=17.84$ and Ne3O2 $=0.89$, linear scale) and metallicity ($12+\log{(\rm{O/H})}=7.84$) that are rarely seen at lower redshifts. We find that the most extreme spectral properties are found in the six $z\gtrsim7$ Ly$\alpha$ emitters in the sample. Each has a hard ionizing spectrum indicating that their visibility is likely enhanced by efficient ionizing photon production. Ly$\alpha$ velocity offsets are found to be very large ($\gtrsim300$ km s$^{-1}$), likely also contributing to their detectability. We find that Ly$\alpha$ in $z\gtrsim7$ galaxies is $6-12\times$ weaker than in lower redshift samples with matched rest-optical spectral properties. If the bubbles around the Ly$\alpha$ emitters are relatively small ($\lesssim0.5-1$ pMpc), we may expect such significant attenuation of Ly$\alpha$ in these ionized regions. We discuss several other effects that may contribute to weaker Ly$\alpha$ emission at $z\gtrsim7$. Deep spectroscopy of fainter galaxies in the vicinity of the Ly$\alpha$ emitters will better characterize the physical scale of the ionized bubbles in this field., Comment: 31 pages, 22 figures, accepted by MNRAS, references updated

Finding and characterising the first galaxies that illuminated the early Universe at cosmic dawn is pivotal to understand the physical conditions and the processes that led to the formation of the first stars. In the first few months of operations, imaging from the James Webb Space Telescope (JWST) have been used to identify tens of candidates of galaxies at redshift (z) greater than 10, less than 450 million years after the Big Bang. However, none of these candidates has yet been confirmed spectroscopically, leaving open the possibility that they are actually low-redshift interlopers. Here we present spectroscopic confirmation and analysis of four galaxies unambiguously detected at redshift 10.3

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

Spitzer/IRAC imaging has revealed that the brightest $z\sim7-8$ galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency among early galaxies. However, IRAC's limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here we use CEERS JWST/NIRCam data to characterize 116 UV-faint (median M$_{UV}=-19.5$) $z\sim6.5-8$ galaxies. The SEDs are typically dominated by young ($\sim$10-50 Myr), low-mass ($M_\ast\sim10^8\ M_\odot$) stellar populations, and we find no need for extremely high stellar masses ($\sim10^{11} M_\odot$). Considering previous studies of UV-bright (M$_{UV}\sim-22$) $z\sim7-8$ galaxies, we find evidence for a strong (5-10$\times$) increase in specific star formation rate toward lower luminosities (median sSFR=103 Gyr$^{-1}$ in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median $\xi_{ion}=10^{25.7}$ erg$^{-1}$ Hz). In spite of their much larger sSFRs, we find no significant increase in [OIII]$+$H$\beta$ EWs towards fainter M$_{UV}$ (median $\approx$780 $\mathring{A}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities ($\lesssim$3% $Z_\odot$) where [OIII] emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our CEERS galaxies. While the majority of our objects are very blue (median $\beta=-2.0$), we identify a significant tail of very dusty galaxies ($\beta\sim-1$) at $\approx$0.5$L_{UV}^\ast$ which may contribute significantly to the $z\sim7-8$ star formation rate density., Comment: Accepted in MNRAS. Updated to use the most recent NIRCam zeropoints. There are no significant changes to the conclusions relative to v1

The ultraviolet (UV) continuum slope ($\beta~$where$~f_\lambda\propto\lambda^\beta$) of galaxies is sensitive to a variety of properties, from the metallicity and age of the stellar population to dust attenuation throughout the galaxy. Considerable attention has focused on identifying reionization-era galaxies with very blue UV slopes ($\beta<-3$). Not only do such systems provide a signpost of low-metallicity stars, but they also identify galaxies likely to leak ionizing photons from their HII regions as such blue UV slopes require the reddening effect of nebular continuum to be diminished. In this paper we present a search for reionization-era galaxies with very blue UV colors in recent JWST/NIRCam imaging of the EGS field. We characterize UV slopes for a large sample of$~z\simeq7-11~$galaxies, finding a median of$~\beta=-2.0$. Two lower luminosity (M$_{\rm{UV}}\simeq-19.5$) and lower stellar mass (6-10$\times10^7$M$_\odot$) systems exhibit extremely blue UV slopes ($\beta=-2.9~$to$~-3.1$) and rest-optical photometry indicating weak nebular line emission. Each system is very compact (r$_e\lesssim$260pc) with very high star-formation-rate surface densities. We model the SEDs with a suite of BEAGLE models with varying levels of ionizing photon escape. The SEDs cannot be reproduced with our fiducial (f$_{\rm{esc,HII}}$=0) or alpha-enhanced (Z$_\star7$., Comment: 13 pages, 10 figures, 3 tables; Accepted to ApJ

With JWST, new opportunities to study the evolution of galaxies in the early Universe are emerging. Spitzer constraints on rest-optical properties of $z\gtrsim7$ galaxies demonstrated the power of using galaxy stellar masses and star formation histories (SFHs) to indirectly infer the cosmic star formation history. However, only the brightest individual $z\gtrsim8$ objects could be detected with Spitzer, making it difficult to robustly constrain activity at $z\gtrsim10$. Here, we leverage the greatly improved rest-optical sensitivity of JWST at $z\gtrsim8$ to constrain the ages of seven UV-bright ($M_{UV}\lesssim-19.5$) galaxies selected to lie at $z\sim8.5-11$, then investigate implications for $z\gtrsim15$ star formation. We infer the properties of individual objects with two spectral energy distribution modelling codes, then infer a distribution of ages for bright $z\sim8.5-11$ galaxies. We find a median age of $\sim20$ Myr, younger than that inferred at $z\sim7$ with a similar analysis, consistent with an evolution towards larger specific star formation rates at early times. The age distribution suggests that only $\sim3$ percent of bright $z\sim8.5-11$ galaxies would be similarly luminous at $z\gtrsim15$, implying that the number density of bright galaxies declines by at least an order of magnitude between $z\sim8.5-11$ and $z\sim15$. This evolution is challenging to reconcile with some early JWST results suggesting the abundance of bright galaxies does not significantly decrease towards very early times, but we suggest this tension may be eased if young stellar populations form on top of older stellar components, or if bright $z\sim15$ galaxies are observed during a burst of star formation., Comment: 16 pages, 10 figures, 1 appendix, accepted to MNRAS

The onset of the {\it JWST}-era provides a much-improved opportunity to characterize the resolved structure of early star forming systems. Previous {\it Spitzer} observations of $z\gtrsim 6$ galaxies revealed the presence of old stars and luminous HII regions (via [OIII]+H$\beta$ emission), but the poor resolution stunted our ability to map their locations with respect to the star forming regions identified in the rest-UV. In this paper, we investigate the internal structure of 12 of the most luminous $z\simeq 6-8$ galaxies in the EGS field observed with recent {\it JWST}/NIRCam imaging. The systems appear clumpy in the rest-UV, with more than half of the light coming from $\simeq 10^7$ to 10$^{9}$ M$_\odot$ star forming complexes that are $\simeq 150$ - 480 pc in size. The clumps tend to be dominated by young stars (median = 36 Myr), but we also find large variations in clump ages within individual galaxies. The [OIII]+H$\beta$ EW varies significantly across individual galaxies (reflecting differences in stellar and gas properties), but the HII regions largely track the UV-bright complexes. Perhaps surprisingly, the rest-optical continuum is just as clumpy as the UV, and we do not find older (and redder) nuclear stellar components that were previously undetected or faint in the UV. The majority of the stellar mass in bright $6

The variety of star formation histories (SFHs) of $z\gtrsim6$ galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations suggest that many $z\sim6-9$ galaxies are dominated by $\gtrsim200$ Myr stellar populations, implying significant star formation at $z\gtrsim9$, while others find that most reionization era galaxies are $\lesssim10$ Myr, consistent with little $z\gtrsim9$ star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{UV}\lesssim-21$) galaxies colour-selected to lie at $z\simeq6.6-6.9$, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of $\sim30-70$ Myr depending on SED model. We quantify the fractions of $\leq10$ Myr and $\geq250$ Myr galaxies, finding that these systems comprise $\sim15-30$ percent and $\sim20-25$ percent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for $z\gtrsim9$ stellar mass assembly and discuss improvements expected from JWST., Comment: 24 pages, 20 figures, 2 appendices, accepted to MNRAS

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates., Comment: 18 pages, 10 figures, submitted to MNRAS

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis., Comment: 32 pages, 20 figures. Submitted to the Astrophysical Journal

The reionization of hydrogen is closely linked to the first structures in the universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly$\alpha$) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly$\alpha$ observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity-halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly$\alpha$ visibility during reionization. We show that UV luminosity-halo mass scatter reduces Ly$\alpha$ visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{HI}}$, at $z \sim 7$ using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{HI}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{HI}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies., Comment: 14 pages, 7 figures, accepted for publication in MNRAS