Search

ASPIRE (A SPectroscopic survey of bIased halos in the Reionization Era) is a quasar legacy survey primarily using JWST to target a sample of 25 $z>6$ quasars with NIRCam slitless spectroscopy and imaging. The first study in this series found evidence of a strong overdensity of galaxies around J0305-3150, a luminous quasar at $z=6.61$, within a single NIRCam pointing obtained in JWST Cycle 1. Here, we present the first results of a JWST Cycle 2 mosaic that covers 35 arcmin$^2$ with NIRCam imaging/WFSS of the same field to investigate the spatial extent of the putative protocluster. The F356W grism data targets [OIII]+H$\beta$ at $5.310$ cMpc., Comment: submitted to ApJ june 2024. 18 pages, 7 figures. companion paper II by Champagne+24b is also available on arXiv today

We present indirect constraints on the absolute escape fraction of ionizing photons ($f_{\rm esc}^{\rm LyC}$) of the system GN 42912 which comprises two luminous galaxies ($M_{\rm UV}$ magnitudes of -20.89 and -20.37) at $z\sim7.5$, GN 42912-NE and GN 42912-SW, to determine their contribution to the ionizing photon budget of the Epoch of Reionization (EoR). The high-resolution James Webb Space Telescope NIRSpec and NIRCam observations reveal the two galaxies are separated by only ~0.1$"$ (0.5 kpc) on the sky and have a 358 km s$^{-1}$ velocity separation. GN 42912-NE and GN 42912-SW are relatively massive for this redshift (log($M_\ast/M_\odot$) $\sim$ 8.4 and 8.9, respectively), with gas-phase metallicities of 18 per cent and 23 per cent solar, O$_{32}$ ratios of 5.3 and $>5.8$, and $\beta$ slopes of $-1.92$ and $-1.51$, respectively. We use the Mg II$\lambda\lambda$2796,2803 doublet to constrain $f_{\rm esc}^{\rm LyC}$. Mg II has an ionization potential close to that of neutral hydrogen and, in the optically thin regime, can be used as an indirect tracer of the LyC leakage. We establish realistic conservative upper limits on $f_{\rm esc}^{\rm LyC}$ of 8.5 per cent for GN 42912-NE and 14 per cent for GN 42912-SW. These estimates align with $f_{\rm esc}^{\rm LyC}$ trends observed with $\beta$, O$_{32}$, and the H$\beta$ equivalent width at $z<4$. The small inferred ionized region sizes ($<0.3$ pMpc) around both galaxies indicate they have not ionized a significant fraction of the surrounding neutral gas. While these $z>7$ $f_{\rm esc}^{\rm LyC}$ constraints do not decisively determine a specific reionization model, they support a minor contribution from these two relatively luminous galaxies to the EoR., Comment: 18 pages, submitted to MNRAS

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

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

The James Webb Space Telescope (JWST) has uncovered a ubiquitous population of dust-obscured compact sources at $z\gtrsim 4$. Many of these objects exhibit signs of active galactic nucleus (AGN) activity, making their study crucial for understanding the formation of supermassive black holes (SMBHs) and their growth with host galaxies. In this work, we examine low and medium resolution JWST/NIRSpec spectra from the JADES GTO public data release in the GOODS-N field of a red, luminous ($M_{\rm B}\sim-22.2$ mag) and compact ($<500$ pc) source at $z=4.13$. The rest-optical ($\lambda_{\rm rest} > 4000$ A) continuum of this source is strongly dominated by a massive (log$_{10}$[$M_*/M_\odot] \sim 10.6$), quenched (log$_{10}$[sSFR/yr$^{-1}$] $< -11$) galaxy, as indicated by the clear presence of a Balmer break and stellar absorption lines. Star-formation history modeling reveals a starburst episode followed by rapid quenching about 200 Myr ago. The spectrum shows extremely broad (FWHM $\sim 2500$ km/s) H$\alpha$ emission and elevated optical line ratios, indicating an actively accreting SMBH. Moreover, our work has potentially revealed clear AGN signatures in the rest-UV in LRDs for the first time, via a detection of a strong Ly$\alpha$ emission and a broad MgII, doublet. The derived black hole mass of log$_{10}$($M_{\rm BH}/M_\odot) \sim 7.3$ results in $M_{\rm BH}/M_*\sim 0.04$ %, consistent with the local relations, unlike the elevated ratios in other high-$z$ reddened AGN. Finally, we use JWST data from AGN at $z=4-10$ to explore an evolutionary link between high-$z$ reddened AGN, early quiescent galaxies, and local ellipticals., Comment: 20 pages, 7 figures, 2 tables. Submitted to ApJ

JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal poor in oxygen and carbon but heavily-enriched in nitrogen. This abundance pattern is almost never seen in lower redshift galaxies but is commonly observed in globular cluster stars. We have recently demonstrated that this peculiar abundance pattern appears in a compact ($\simeq 20$ pc) metal-poor galaxy undergoing a strong burst of star formation. This galaxy was originally selected based on strong CIV emission, indicating a hard radiation field rarely seen locally. In this paper, we present JWST/NIRSpec observations of another reionization-era galaxy known to power strong CIV emission, the $z=7.04$ gravitationally-lensed galaxy A1703-zd6. The emission line spectrum reveals this is a metal poor galaxy ($12+\log(\rm O/H) = 7.47\pm0.19$) dominated by a young stellar population ($1.6^{+0.5}_{-0.4}$ Myr) that powers a very hard ionizing spectrum (CIV EW = 19.4 $\unicode{x212B}$, He II EW = 2.2 $\unicode{x212B}$). The ISM is highly-enriched in nitrogen ($\log(\rm N/O)=-0.6$) with very high electron densities ($8-19\times10^4$ cm$^{-3}$) and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission (EW$\gtrsim20$ $\unicode{x212B}$) in two new $z\gtrsim 6$ metal poor galaxies. To put these results in context, we search for UV line emission in a sample of 737 $z\gtrsim 4$ galaxies with NIRSpec spectra, establishing that 40(30)% of systems with [OIII]+H$\beta$ EW $>2000\unicode{x212B}$ have NIV] (CIV) detections with EW$>5$ $\unicode{x212B}$ ($>10$ $\unicode{x212B}$). These results suggest high N/O ratios and hard ionizing sources appear in a brief phase following a burst of star formation in compact high density stellar complexes., Comment: 29 pages, 18 figures; submitted to ApJ

Low-luminosity AGNs with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their co-evolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we perform a systematic search for broad-line ${\rm H\alpha}$ emitters (BHAEs) at $z\approx 4-5$ in 25 fields of the ASPIRE (A SPectroscopic survey of biased halos In the Reionization Era) project, covering a total area of 275 arcmin$^2$. We identify 16 BHAEs with FWHM of the broad components spanning from $\sim$ 1000 km s$^{-1}$ to 3000 km s$^{-1}$. Assuming the broad linewidths arise due to Doppler broadening around BHs, the implied BH masses range from $10^7$ to $10^{8}~M_\odot$, with broad ${\rm H\alpha}$-converted bolometric luminosity of $10^{44.5}-10^{45.5}$ erg s$^{-1}$ and Eddington ratios of $0.07-0.47$. The spatially extended structure of the F200W stacked image may trace the stellar light from the host galaxies. The ${\rm H\alpha}$ luminosity function indicates an increasing AGN fraction towards the higher ${\rm H\alpha}$ luminosities. We find possible evidence for clustering of BHAEs: two sources are at the same redshift with a projected separation of 519 kpc; one BHAE appears as a composite system residing in an overdense region with three close companion ${\rm H\alpha}$ emitters. Three BHAEs exhibit blueshifted absorption troughs indicative of the presence of high-column-density gas. We find the broad-line and photometrically selected BHAE samples exhibit different distributions in the optical continuum slopes, which can be attributed to their different selection methods. The ASPIRE broad-line ${\rm H\alpha}$ sample provides a good database for future studies of faint AGN populations at high redshift., Comment: 19 pages, 13 figures, 4 tables. Accepted by the ApJ

This paper presents an overview of the high-redshift extragalactic science case for the Single Aperture Large Telescope for Universe Studies (SALTUS) far-infrared NASA probe-class mission concept. Enabled by its 14m primary reflector, SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions. SALTUS would be a versatile observatory capable of responding to the scientific needs of the extragalactic community in the 2030s, and a natural follow-on to the near- and mid-IR capabilities of JWST. Key early-universe science goals for SALTUS focus on understanding the role of galactic feedback processes in regulating galaxy growth across cosmic time, and charting the rise of metals and dust from the early universe to the present. This paper summarizes these science cases and the performance metrics most relevant for high-redshift observations., Comment: Submitted to JATIS as part of a special issue on SALTUS. Key performance metrics in Figs. 1 and 2

Extreme emission line galaxies (EELGs) exhibit large equivalent widths (EW) in their rest-optical emission lines ([OIII]$\lambda5007$ or H$\alpha$ rest-frame EW$ > 750\r{A}$) which can be tied to a recent upturn in star formation rate, due to the sensitivity of the nebular line emission and the rest-optical continuum to young ($<10$Myr) and evolved stellar populations, respectively. By studying a sample of 85 star forming galaxies (SFGs), spanning the redshift and magnitude interval $3 $ M$_{UV}>-21$, in the JWST Advanced Deep Extragalactic Survey (JADES) with NIRSpec/prism spectroscopy, we determine that SFGs initiate an EELG phase when entering a significant burst of star formation, with the highest EWs observed in EELGs with the youngest luminosity-weighted ages ($<5$ Myr old) and the highest burst intensity (those with the greatest excess between their current and long-term average SFR). We spectroscopically confirm that a greater proportion of SFGs are in an EELG phase at high redshift in our UV-selected sample ($61\pm4\%$ in our $z>5.7$ high-redshift bin, compared to $23^{+4}_{-1}\%$ in our lowest-redshift bin $35.7$ have observed Lyman-$\alpha$ emission, potentially lying within large ionised regions. The high detection rate of Lyman-$\alpha$ emitters in our EELG selection suggests that the physical conditions associated with entering an EELG phase also promote the escape of Lyman-$\alpha$ photons., Comment: 34 pages, 25 figures

Nearly a decade ago, we began to see indications that reionization-era galaxies power hard radiation fields rarely seen at lower redshift. Most striking were detections of nebular CIV emission in what appeared to be typical low mass galaxies, requiring an ample supply of 48 eV photons to triply ionize carbon. The nature of this population has long remained unclear owing to limitations of ground-based spectroscopy. We have obtained deep JWST/NIRSpec R=1000 spectroscopy of the two z>6 CIV-emitting galaxies known prior to JWST. Here we present a rest-UV to optical spectrum of one of these two systems, the multiply-imaged z=6.1 lensed galaxy RXCJ2248-ID. NIRCam imaging reveals two compact (<22pc) clumps separated by 220pc, with one comprising a dense concentration of massive stars ($>10,400M_{\odot}$/yr/kpc$^2$) formed in a recent burst. We stack spectra of 3 images of the galaxy (J=24.8-25.9), yielding a very deep spectrum providing a high S/N template of strong emission line sources at z>6. The spectrum reveals narrow high ionization lines (HeII, CIV, NIV]) with line ratios consistent with powering by massive stars. The rest-optical spectrum is dominated by very strong emission lines ([OIII] EW=2800\AA), albeit with weak emission from low-ionization transitions ([OIII]/[OII]=184). The electron density is found to be very high($6.4-31\times10^4$cm$^{-3}$) based on three UV transitions. The ionized gas is metal poor ($12+\log(\rm O/H)=7.43^{+0.17}_{-0.09}$), yet highly enriched in nitrogen ($\log(\rm N/O)=-0.39^{+0.11}_{-0.10}$). The spectrum appears broadly similar to that of GNz11 at z=10.6, without showing the same AGN signatures. We suggest that the hard radiation field and rapid nitrogen enrichment may be a short-lived phase that many z>6 galaxies go through as they undergo strong bursts of star formation. We comment on the potential link of such spectra to globular cluster formation., Comment: 22 pages, 13 figures, 4 tables; Accepted to MNRAS

We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $\sigma(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $\sigma(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$., Comment: 28 pages, 12 figures, 2 Appendix figures

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

Using deep JWST imaging from JADES, JEMS and SMILES, we characterize optically-faint and extremely red galaxies at $z>3$ that were previously missing from galaxy census estimates. The data indicate the existence of abundant, dusty and post-starburst-like galaxies down to $10^8$M$_\odot$, below the sensitivity limit of Spitzer and ALMA. Modeling the NIRCam and HST photometry of these red sources can result in extreme, high values for both stellar mass and star formation rate (SFR); however, including 7 MIRI filters out to 21$\mu$m results in decreased mass (median 0.6 dex for log$_{10}$M$^*$/M$_{\odot}>$10), and SFR (median 10$\times$ for SFR$>$100 M$_{\odot}$/yr). At $z>6$, our sample includes a high fraction of little red dots (LRDs; NIRCam-selected dust-reddened AGN candidates). We significantly measure older stellar populations in the LRDs out to rest-frame 3$\mu$m (the stellar bump) and rule out a dominant contribution from hot dust emission, a signature of AGN contamination to stellar population measurements. This allows us to measure their contribution to the cosmic census at $z>3$, below the typical detection limits of ALMA ($L_{\rm IR}<10^{12}L_\odot$). We find that these sources, which are overwhelmingly missed by HST and ALMA, could effectively double the obscured fraction of the star formation rate density at $4

With the remarkable sensitivity and resolution of JWST in the infrared, measuring rest-optical kinematics of galaxies at $z>5$ has become possible for the first time. This study pilots a new method for measuring galaxy dynamics for highly multiplexed, unbiased samples by combining FRESCO NIRCam grism spectroscopy and JADES medium-band imaging. Here we present one of the first JWST kinematic measurements for a galaxy at $z>5$. We find a significant velocity gradient, which, if interpreted as rotation yields $V_{rot} = 240\pm50$km/s and we hence refer to this galaxy as Twister-z5. With a rest-frame optical effective radius of $r_e=2.25$kpc, the high rotation velocity in this galaxy is not due to a compact size as may be expected in the early universe but rather a high total mass, ${\rm log(M}_{dyn}/{\rm M}_\odot)=11.0\pm0.2$. This is a factor of roughly 4x higher than the stellar mass within the effective radius. We also observe that the radial H$\alpha$ equivalent width profile and the specific star formation rate map from resolved stellar population modeling is centrally depressed by a factor of $\sim1.5$ from the center to $r_e$. Combined with the morphology of the line-emitting gas in comparison to the continuum, this centrally suppressed star formation is consistent with a star-forming disk surrounding a bulge growing inside-out. While large, rapidly rotating disks are common to z~2, the existence of one after only 1Gyr of cosmic time, shown for the first time in ionized gas, adds to the growing evidence that some galaxies matured earlier than expected in the history of the universe., Comment: Fig. 3 shows the main result

Recent JWST observations have revealed an unexpected abundance of massive galaxy candidates in the early Universe, extending further in redshift and to lower luminosity than what had previously been found by sub-millimeter surveys. These JWST candidates have been interpreted as challenging the $\Lambda$CDM cosmology, but, so far, they have mostly relied only on rest-frame ultraviolet data and lacked spectroscopic confirmation of their redshifts. Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between $z_{\rm spec}=5-9$ from the JWST FRESCO survey. We find no tension with the $\Lambda$CDM model in our sample. However, three ultra-massive galaxies (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) require an exceptional fraction of 50% of baryons converted into stars -- two to three times higher than even the most efficient galaxies at later epochs. The contribution from an active nucleus is unlikely because of their extended emission. Ultra-massive galaxies account for as much as 17% of the total cosmic star formation rate density at $z\sim5-6$., Comment: Nature in press. Updated to the accepted version. 24 pages, 4 main figures, 7 supplementary figures, 3 supplementary tables

Early JWST observations have uncovered a new population of red sources that might represent a previously overlooked phase of supermassive black hole growth (Kocevski et al. 2023; Matthee et al. 2023, Labb\'e et al. 2023). One of the most intriguing examples is an extremely red, point-like object that was found to be triply-imaged by the strong lensing (SL) cluster Abell 2744 (Furtak et al. 2023). Here we present deep JWST/NIRSpec observations of this object, Abell2744-QSO1. The spectroscopy confirms that the three images are of the same object, and that it is a highly reddened ($A_V\simeq3$) broad emission-line Active Galactic Nucleus (AGN) at a redshift of $z_{\mathrm{spec}}=7.0451\pm0.0005$. From the width of H$\beta$ ($\mathrm{FWHM}=2800\pm250\,\frac{\mathrm{km}}{\mathrm{s}}$) we derive a black hole mass of $M_{\mathrm{BH}}=4_{-1}^{+2}\times10^7\,\mathrm{M}_{\odot}$. We infer a very high ratio of black hole to galaxy mass of at least 3%, an order of magnitude more than is seen in local galaxies (Bennert et al. 2011), and possibly as high as 100%. The lack of strong metal lines in the spectrum together with the high bolometric luminosity ($L_{\mathrm{bol}}=(1.1\pm0.3)\times10^{45}\,\frac{\mathrm{erg}}{\mathrm{s}}$) indicate that we are seeing the black hole in a phase of rapid growth, accreting at 30% of the Eddington limit. The rapid growth and high black hole to galaxy mass ratio of A2744-QSO1 suggest that it may represent the missing link between black hole seeds (Volonteri et al. 2021) and the first luminous quasars (Fan et al. 2022)., Comment: Published in Nature. Updated to the accepted version

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 study of the environments of 17 Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$\alpha$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$\alpha$ velocity offsets ($\Delta v_\text{Ly$\alpha$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$\alpha$ escape fractions ($f_\mathrm{esc,\,Ly\alpha}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$\alpha$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$\alpha$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs., Comment: 15 pages, 8 figures. Accepted for publication in A&A

We present the properties of 17 faint Ly$\alpha$ emitting galaxies (LAEs) at $z>5.8$ from the JWST Advanced Deep Extragalactic Survey (JADES) in the Hubble Ultra Deep Field/GOODS-S. These LAEs span a redshift range $z\approx5.8-8.0$ and a UV magnitude range $M_{UV}\approx-17$ to $-20.6$, with the Ly$\alpha$ equivalent width (EW) in the range $\approx 25-350$ \AA. The detection of other rest-optical emission lines in the spectra of these LAEs enables the determination of accurate systemic redshifts and Ly{\alpha} velocity offsets, as well as the physical and chemical composition of their stars and interstellar media. These faint LAEs are consistent with metal-poor systems with high ionization parameters, similar to the general galaxy population at $z>6$. We measured an average ionizing photon production efficiency, log($\xi_\rm{ion}$/erg$^{-1}$ Hz) $\approx25.57$ across our LAEs, which does not evolve strongly with redshift. We report an anti-correlation between the Ly$\alpha$ escape fraction (f_\rm{esc}) and the velocity offset from systemic redshift, consistent with model expectations. We further find that the strength and velocity offset of Ly$\alpha$ are neither correlated with galaxy spectroscopic properties nor with $\xi_\rm{ion}$. We find a decrease in $f_\rm{esc}$(Ly$\alpha$) with redshift, indicative of decreasing sizes of ionized bubbles around LAEs at high redshifts. We used a range of galaxy properties to predict Lyman continuum $f_\rm{esc}$ for our LAEs, finding that the ionizing photon output into the intergalactic medium remains roughly constant across the observed Ly$\alpha$ EW, showing a mild increase at fainter M$_{UV}$ and at higher redshifts. We derived correlations between the ionizing photon output from LAEs and $M_{UV}$, Ly$\alpha$ EW and redshifts, which can be used to constrain the ionizing photon contribution of LAEs at $z > 6$ towards cosmic reionization., Comment: Accepted for publication in A&A. 24 pages, 13 figures, spectra presented in the appendix. This updated version includes the addition of one more LAE to the sample, ID 004296 at z=6.712

The physical processes that establish the morphological evolution and the structural diversity of galaxies are key unknowns in extragalactic astrophysics. Here we report the finding of the morphologically-mature galaxy JADES-GS+53.18343-27.79097, which existed within the first 700 million years of the Universe's history. This star-forming galaxy with a stellar mass of 400 million solar masses consists of three components, a highly-compact core with a half-light radius of less than 100 pc, an actively star-forming disc with a radius of about 400 pc, and a star-forming clump, which all show distinctive star-formation histories. The central stellar mass density of this galaxy is within a factor of two of the most massive present-day ellipticals, while being globally 1000 times less massive. The radial profile of the specific star-formation rate is rising toward the outskirts. This evidence suggests the first detection of inside-out growth of a galaxy as a proto-bulge and a star-forming disc in the Epoch of Reionization., Comment: 42 pages, 17 figures

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

We describe the NIRSpec component of the JWST Deep Extragalactic Survey (JADES), and provide deep spectroscopy of 253 sources targeted with the NIRSpec micro-shutter assembly in the Hubble Ultra Deep Field and surrounding GOODS-South. The multi-object spectra presented here are the deepest so far obtained with JWST, amounting to up to 28 hours in the low-dispersion ($R\sim 30-300$) prism, and up to 7 hours in each of the three medium-resolution $R\approx 1000$ gratings and one high-dispersion grating, G395H ($R\approx2700$). Our low-dispersion and medium-dispersion spectra cover the wavelength range $0.6-5.3\mu$m. We describe the selection of the spectroscopic targets, the strategy for the allocation of targets to micro-shutters, and the design of the observations. We present the public release of the reduced 2D and 1D spectra, and a description of the reduction and calibration process. We measure spectroscopic redshifts for 178 of the objects targeted extending up to $z=13.2$. We present a catalog of all emission lines detected at $S/N>5$, and our redshift determinations for the targets. Combined with the first JADES NIRCam data release, these public JADES spectroscopic and imaging datasets provide a new foundation for discoveries of the infrared universe by the worldwide scientific community., Comment: Accepted for publication in A&A. Data products available from https://archive.stsci.edu/hlsp/jades

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

Constraining when and how reionisation began is pivotal for understanding when the first galaxies formed. Lyman-alpha (Ly$\alpha$) emission from galaxies is currently our most promising probe of these early stages. At z>7 the majority of galaxies detected with Ly$\alpha$ are in candidate overdensities. Here we quantify the probability of these galaxies residing in large ionised bubbles. We create (1.6 Gpc)$^3$ reionising intergalactic medium (IGM) simulations, providing sufficient volume to robustly measure bubble size distributions around UV-bright galaxies and rare overdensities. We find $M_{\rm UV} \lesssim -16$ galaxies and overdensities are $\gtrsim$10-1000x more likely to trace ionised bubbles compared to randomly selected positions. The brightest galaxies and strongest overdensities have bubble size distributions with highest characteristic size and least scatter. We compare two models: gradual reionisation driven by numerous UV-faint galaxies versus more rapid reionisation by rarer brighter galaxies, producing larger bubbles at fixed neutral fraction. We demonstrate that recently observed z~7 overdensities are highly likely to trace large ionised bubbles, corroborated by their high Ly$\alpha$ detection rates. However, the z~8.7 association of Ly$\alpha$ emitters in EGS and GN-z11, with Ly$\alpha$ at z=10.6, are unlikely to trace large bubbles in our fiducial model -- 11% and 7% probability of >1 proper Mpc bubbles, respectively. Ly$\alpha$ detections at such high redshifts could be explained by: a less neutral IGM than previously expected; larger ionised regions at fixed neutral fraction; or if intrinsic Ly$\alpha$ flux is unusually strong in these galaxies. We discuss how to test these scenarios with JWST and the prospects for using upcoming wide-area surveys to distinguish between reionisation models., Comment: 15 pages (+ 3 pages in Appendix), 13 figures, submitted to MNRAS

We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $\mu$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $\delta_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3

Local and low-redshift ($z$<$3$) galaxies are known to broadly follow a bimodal distribution: actively star forming galaxies with relatively stable star-formation rates, and passive systems. These two populations are connected by galaxies in relatively slow transition. In contrast, theory predicts that star formation was stochastic at early cosmic times and in low-mass systems: these galaxies transitioned rapidly between starburst episodes and phases of suppressed star formation, potentially even causing temporary quiescence -- so-called mini-quenching events. However, the regime of star-formation burstiness is observationally highly unconstrained. Directly observing mini-quenched galaxies in the primordial Universe is therefore of utmost importance to constrain models of galaxy formation and transformation. Early quenched galaxies have been identified out to redshift $z<5$, and these are all found to be massive ($M_{*}>10^{10}~M_{\odot}$) and relatively old. Here we report a (mini-)quenched galaxy at z$=$7.3, when the Universe was only 700~Myr old. The JWST/NIRSpec spectrum is very blue ($U$-$V$$=$0.16$\pm0.03$~mag), but exhibits a Balmer break and no nebular emission lines. The galaxy experienced a short starburst followed by rapid quenching; its stellar mass (4-6$\times 10^8~M_\odot$) falls in a range that is sensitive to various feedback mechanisms, which can result in perhaps only temporary quenching., Comment: Accepted for publication in Nature

We report the discovery of a remarkable Ly$\alpha$ emitting galaxy at z = 7.2782, JADES-GS-z7-LA, with EW$_0$(Ly$\alpha$) $= 388.0 \pm 88.8$\AA and UV magnitude -17.0. The spectroscopic redshift is confirmed via rest-frame optical lines [O II], H$\beta$ and [O III] in its JWST/NIRSpec Micro-Shutter Assembly (MSA) spectrum. The Ly$\alpha$ line is detected in both lower resolution PRISM as well as medium resolution G140M grating spectra. The LSF-deconvolved Ly$\alpha$ FWHM in the grating is $383.9 \pm 56.2$ km/s and the Ly$\alpha$ velocity offset compared to the systemic redshift is $113.3 \pm 80.0$ km/s, indicative of very little neutral gas or dust within the galaxy. We estimate the Ly$\alpha$ escape fraction to be >70%. JADES-GS-z7-LA has a O32 ratio of $11.1 \pm 2.2$ and a R23 ratio of $11.2 \pm 2.6$, consistent with low metallicity and high ionization parameters. Deep NIRCam imaging also revealed a close companion source (separated by 0.23"), which exhibits similar photometry to that of JADES-GS-z7-LA, with a photometric excess in the F410M NIRCam image consistent with [O III]+H$\beta$ emission at the same redshift. The spectral energy distribution of JADES-GS-z7-LA indicates a "bursty" star formation history, with a low stellar mass of $\approx 10^7$ $M_\odot$. Assuming that the Ly$\alpha$ transmission through the intergalactic medium is the same as its measured escape fraction, an ionized region of size > 1.5 pMpc is needed to explain the high Ly$\alpha$ EW and low velocity offset compared to systemic seen in JADES-GS-z7-LA. Owing to its UV-faintness, we show that it is incapable of single-handedly ionizing a region large enough to explain its Ly$\alpha$ emission. Therefore, we suggest that JADES-GS-z7-LA (and possibly the companion source) may be a part of a larger overdensity, presenting direct evidence of overlapping ionized bubbles at $z>7$., Comment: 14 pages, 8 figures, Accepted for publication in A&A

We report the discovery of an extreme galaxy overdensity at $z = 5.4$ in the GOODS-S field using JWST/NIRCam imaging from JADES and JEMS alongside JWST/NIRCam wide field slitless spectroscopy from FRESCO. We identified potential members of the overdensity using HST+JWST photometry spanning $\lambda = 0.4-5.0\ \mu\mathrm{m}$. These data provide accurate and well-constrained photometric redshifts down to $m \approx 29-30\,\mathrm{mag}$. We subsequently confirmed $N = 81$ galaxies at $5.2 < z < 5.5$ using JWST slitless spectroscopy over $\lambda = 3.9-5.0\ \mu\mathrm{m}$ through a targeted line search for $\mathrm{H} \alpha$ around the best-fit photometric redshift. We verified that $N = 42$ of these galaxies reside in the field while $N = 39$ galaxies reside in a density around $\sim 10$ times that of a random volume. Stellar populations for these galaxies were inferred from the photometry and used to construct the star-forming main sequence, where protocluster members appeared more massive and exhibited earlier star formation (and thus older stellar populations) when compared to their field galaxy counterparts. We estimate the total halo mass of this large-scale structure to be $12.6 \lesssim \mathrm{log}_{10} \left( M_{\mathrm{halo}}/M_{\odot} \right) \lesssim 12.8$ using an empirical stellar mass to halo mass relation, which is likely an underestimate as a result of incompleteness. Our discovery demonstrates the power of JWST at constraining dark matter halo assembly and galaxy formation at very early cosmic times., Comment: Resubmitted to ApJ based on reviewer report; main text has 15 pages, 6 figures and 1 table; appendix has 1 page, 2 figure sets, and 2 tables

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

Large dust reservoirs (up to $\sim 10^8 \, \mathrm{M_\odot}$) have been detected in galaxies out to redshift $z \sim 8$, when the age of the universe was only about 600 Myr. Generating significant amounts of dust within such a short timescale has proven challenging for theories of dust formation and has prompted the revision of the modelling of potential sites of dust production such as the atmospheres of asymptotic giant branch (AGB) stars in low-metallicity environments, supernovae (SNe) ejecta, and the accelerated growth of grains in the interstellar medium (ISM). However, degeneracies between different evolutionary pathways remain when the total dust mass of galaxies is the only available observable. Here we report observations of the $2175 \, \mathring{\rm A}$ dust attenuation feature, well known in the Milky Way (MW) and galaxies at $z \lesssim 3$, in the near-infrared spectra of galaxies up to $z \sim 7$, corresponding to the first billion years of cosmic time. The relatively short timescale implied for the formation of carbonaceous grains giving rise to this feature suggests a rapid production process, likely in Wolf-Rayet (WR) stars or SN ejecta., Comment: 28 pages, 8 figures. Accepted for publication in Nature

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

Recent JWST/NIRCam imaging taken for the ultra-deep UNCOVER program reveals a very red dropout object at $z_{\mathrm{phot}}\simeq7.6$, triply imaged by the galaxy cluster Abell 2744 ($z_{\mathrm{d}}=0.308$). All three images are very compact, i.e. unresolved, with a de-lensed size upper-limit of $r_{e}\lesssim35$ pc. The images have apparent magnitudes of $m_{\mathrm{F444W}}\sim25-26$ AB, and the magnification-corrected absolute UV magnitude of the source is $M_{\mathrm{UV},1450}=-16.81\pm0.09$. From the sum of observed fluxes and from a spectral energy distribution (SED) analysis, we obtain estimates of the bolometric luminosities of the source of $L_{\mathrm{bol}}\gtrsim10^{43} \frac{\mathrm{erg}}{\mathrm{s}}$ and $L_{\mathrm{bol}}\sim10^{44}-10^{46} \frac{\mathrm{erg}}{\mathrm{s}}$, respectively. Based on its compact, point-like appearance, its position in color-color space and the SED analysis, we tentatively conclude that this object is a UV-faint dust-obscured quasar-like object, i.e. an active galactic nucleus (AGN) at high redshift. We also discuss other alternative origins for the object's emission features, including a massive star cluster, Population III, supermassive, or dark stars, or a direct-collapse black hole. Although populations of red galaxies at similar photometric redshifts have been detected with JWST, this object is unique in that its high-redshift nature is corroborated geometrically by lensing, that it is unresolved despite being magnified -- and thus intrinsically even more compact -- and that it occupies notably distinct regions in both size-luminosity and color-color space. Planned UNCOVER JWST/NIRSpec observations, scheduled in Cycle 1, will enable a more detailed analysis of this object., Comment: Accepted for publication in ApJ. Considerable updates to the SED-fitting section, including the discussion of a heavily dust-obscured type 1 AGN as the most likely scenario

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

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

Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization., Comment: MNRAS Accepted (01 Nov 2022)

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 band 6 ALMA observations of a heavily-obscured radio-loud ($L_{1.4\ \mathrm{GHz}}=10^{25.4}$ W Hz$^{-1}$) AGN candidate at $z_\mathrm{phot}=6.83\pm0.06$ found in the 1.5 deg$^2$ COSMOS field. The ALMA data reveal detections of exceptionally strong [CII]158$\mu$m ($z_\mathrm{[CII]}=6.8532$) and underlying dust continuum emission from this object (COS-87259), where the [CII] line luminosity, line width, and 158$\mu$m continuum luminosity are comparable to that seen from $z\sim7$ sub-mm galaxies and quasar hosts. The 158$\mu$m continuum detection suggests a total infrared luminosity of $9\times10^{12}$ $L_\odot$ with corresponding very large obscured star formation rate (1300 $M_\odot$/yr) and dust mass ($2\times10^9$ $M_\odot$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization era galaxy with $M_\ast\approx1.7\times10^{11}$ $M_\odot$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbors an AGN that is heavily obscured ($\tau_{_{\mathrm{9.7\mu m}}}=2.3$) with a bolometric luminosity of approximately $5\times10^{13}$ $L_\odot$. Such a very high AGN luminosity suggests this object is powered by an $\approx$1.6 $\times$ 10$^9$ $M_\odot$ black hole if accreting near the Eddington limit, and is effectively a highly-obscured version of an extremely UV-luminous ($M_{1450}\approx-27.3$) $z\sim7$ quasar. Notably, these $z\sim7$ quasars are an exceedingly rare population ($\sim$0.001 deg$^{-2}$) while COS-87259 was identified over a relatively small field. Future very wide-area surveys with, e.g., Roman and Euclid have the potential to identify many more extremely red yet UV-bright $z\gtrsim7$ objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population., Comment: 12 pages, 7 figures, 1 table. Updated to accepted version (MNRAS)

We report 18 dust continuum detections ($\geq 3.3\sigma$) at $\sim88{\rm \mu m}$ and $158{\rm \mu m}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm \mu m}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$., Comment: 17 pages, 9 figures, 4 tables, Accepted for publication in 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

Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$\alpha$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly$\alpha$ profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at $z\sim7$ remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [CII]-based redshifts from ALMA. Here we report MMT/Binospec Ly$\alpha$ spectroscopy of eight UV-bright ($\mathrm{M_{UV}}^{}\sim-22$) galaxies at $z\simeq7$ selected from the ALMA REBELS survey. We detect Ly$\alpha$ in 4 of 8 galaxies and use the [CII] systemic redshifts to investigate the Ly$\alpha$ velocity profiles. The Ly$\alpha$ lines are significantly redshifted from systemic (average velocity offset=223 km/s) and broad (FWHM$\approx$300$-$650 km/s), with two sources showing emission extending to $\approx$750 km/s. We find that the broadest Ly$\alpha$ profiles are associated with the largest [CII] line widths, suggesting a potential link between the Ly$\alpha$ FWHM and the dynamical mass. Since Ly$\alpha$ photons at high velocities transmit efficiently through the $z=7$ IGM, our data suggest that velocity profiles play a significant role in boosting the Ly$\alpha$ visibility of the most UV-luminous reionization-era galaxies., Comment: 17 pages, 7 figures, submitted to MNRAS

Our understanding of reionization has advanced considerably over the past decade, with several results now demonstrating that the IGM transitioned from substantially neutral at $z=7$ to largely reionized at $z=6$. However, little remains known about the sizes of ionized bubbles at $z\gtrsim7$ as well as the galaxy overdensities which drive their growth. Fortunately, rest-UV spectroscopic observations offer a pathway towards characterizing these ionized bubbles thanks to the resonant nature of Lyman-alpha photons. In a previous work, we presented Ly$\alpha$ detections from three closely-separated Lyman-break galaxies at $z\simeq6.8$, suggesting the presence of a large ($R>1$ physical Mpc) ionized bubble in the 1.5 deg$^2$ COSMOS field. Here, we present new deep Ly$\alpha$ spectra of ten UV-bright ($\mathrm{M}_{\mathrm{UV}}^{} \leq -20.4$) $z\simeq6.6-6.9$ galaxies in the surrounding area, enabling us to better characterize this potential ionized bubble. We confidently detect (S/N$>$7) Ly$\alpha$ emission at $z=6.701-6.882$ in nine of ten observed galaxies, revealing that the large-scale volume spanned by these sources (characteristic radius $R = 3.2$ physical Mpc) traces a strong galaxy overdensity ($N/\langle N\rangle \gtrsim 3$). Our data additionally confirm that the Ly$\alpha$ emission of UV-bright galaxies in this volume is significantly enhanced, with 40% (4/10) showing strong Ly$\alpha$ emission (equivalent width$>$25 $\mathrm{\mathring{A}}$) compared to the 8$-$9% found on average at $z\sim7$. The median Ly$\alpha$ equivalent width of our observed galaxies is also $\approx$2$\times$ that typical at $z\sim7$, consistent with expectations if a very large ($R\sim3$ physical Mpc) ionized bubble is allowing the Ly$\alpha$ photons to cosmologically redshift far into the damping wing before encountering HI., Comment: 13 pages, 3 figures, 3 tables. Accepted in MNRAS with minor revisions