Search

We present a resolved study of $>900$ progenitors of Milky Way Analogs (MWAs) at $0.3

The James Webb Space Telescope (JWST) has recently discovered a new population of objects at high redshift referred to as `Little Red Dots' (LRDs). Their nature currently remains elusive, despite their surprisingly high inferred number densities. This emerging population of red point-like sources is reshaping our view of the early Universe and may shed light on the formation of high-redshift supermassive black holes. Here we present a spectroscopically confirmed LRD CANUCS-LRD-z8.6 at $z_{\rm spec}=8.6319\pm 0.0005$ hosting an Active Galactic Nucleus (AGN), using JWST data. This source shows the typical spectral shape of an LRD (blue UV and red optical continuum, unresolved in JWST imaging), along with broad H$\beta$ line emission, detection of high-ionization emission lines (CIV, NIV]) and very high electron temperature indicative of the presence of AGN. This is also combined with a very low metallicity ($Z<0.1 Z_\odot$). The presence of all these diverse features in one source makes CANUCS-LRD-z8.6 unique. We show that the inferred black hole mass of CANUCS-LRD-z8.6 ($M_{\rm BH}=1.0^{+0.6}_{-0.4}\times 10^{8}\rm ~M_\odot$) strongly challenges current standard theoretical models and simulations of black hole formation, and forces us to adopt `ad hoc' prescriptions. Indeed if massive seeds, or light seeds with super-Eddington accretion, are considered, the observed BH mass of CANUCS-LRD-z8.6 at $z=8.6$ can be reproduced. Moreover, the black hole is over-massive compared to its host, relative to the local $M_{\rm BH}-M_*$ relations, pointing towards an earlier and faster evolution of the black hole compared to its host galaxy., Comment: 4 main figures; 8 supplementary figures; 5 supplementary tables

We examine the quiescent fractions of massive galaxies in six $z\gtrsim3$ spectroscopically-confirmed protoclusters in the COSMOS field, one of which is newly confirmed and presented here. We report the spectroscopic confirmation of MAGAZ3NE~J100143+023021 at $z=3.122^{+0.007}_{-0.004}$ by the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) survey. MAGAZ3NE~J100143+023021 contains a total of 79 protocluster members (28 spectroscopic and 51 photometric). Three spectroscopically-confirmed members are star-forming ultra-massive galaxies ($\log(M_{\star}/{\rm M}_\odot)>11$; UMGs), the most massive of which has $\log(M_{\star}/{\rm M}_\odot)=11.15^{+0.05}_{-0.06}$. Combining Keck/MOSFIRE spectroscopy and the COSMOS2020 photometric catalog, we use a weighted Gaussian kernel density estimator to map the protocluster and measure its total mass $2.25^{+1.55}_{-0.65}\times10^{14}~{\rm M}_{\odot}$ in the dense ``core'' region. For each of the six COSMOS protoclusters, we compare the quiescent fraction to the status of the central UMG as star-forming or quiescent. We observe that galaxies in these protoclusters appear to obey galactic conformity: elevated quiescent fractions are found in protoclusters with $UVJ$ quiescent UMGs and low quiescent fractions are found in protoclusters containing $UVJ$ star-forming UMGs. This correlation of star-formation/quiescence in UMGs and the massive galaxies nearby in these protoclusters is the first evidence for the existence of galactic conformity at $z>3$. Despite disagreements over mechanisms behind conformity at low redshifts, its presence at these early cosmic times would provide strong constraints on the physics proposed to drive galactic conformity., Comment: 6 figures, 5 tables. Accepted for publication in the Astrophysical Journal

We present a comprehensive analysis of the observed Stellar-to-Halo mass relationship (SHMR) spanning redshifts from 0.2 to 4.5. This was enabled through galaxy clustering and abundance measurements from two large (effective area ~ 1.61 deg^2) and homogeneously prepared photometric catalogs - UltraVISTA ultra-deep stripes DR3 (COSMOS) and FENIKS v1 (UDS). To translate these measurements into the SHMR, we introduce a novel halo occupation distribution (HOD) fitting approach (``smooth-$z$'') whereby HOD parameters between neighboring z-bins are connected via physically motivated continuity (smoothing) priors. As a result, the high constraining power at z <~ 2, due to a much wider dynamical range in stellar mass (~ 3 dex), helps constrain the SHMR at z >~ 2, where that range shrinks down to <~ 1 dex. We find that the halo mass is tightly coupled to star formation: the halo mass with peak integrated star-forming efficiency (SFE), M_h^peak remains constant within ~ 10^12.2 - 10^12.4 Msolar throughout the redshifts probed. Furthermore, we show that if we had relied on COSMOS alone (as opposed to COSMOS+UDS), as has been done by many preceding studies, M_h^peak would be systematically lower by up to ~0.15 dex at z < 1.5, highlighting the importance of mitigating cosmic variance. Finally, for the first time, we show how the SFE evolves with redshift as halos grow in mass along their progenitor merger trees, instead of at fixed halo masses., Comment: Submitted to ApJ. 36 pages, 18 figures. Comments welcome

We present a new analytical model for the attenuation to Epoch of Reionization (EoR) galaxies by proximate neutral hydrogen gas. Many galaxy spectra in the EoR taken by JWST have shown a flux deficit at wavelengths just redward of the Lyman break, and this has been regarded as resulting from Ly$\alpha$ damping wing absorption by the increasing amount of neutral hydrogen in the line-of-sight. However, previous attenuation models for the intergalactic medium (IGM) commonly used in photometric redshift template-fitting codes assume that the Lyman break is rather sharp, which leads to systematic overestimation of photometric redshifts at $z>7$. In this letter, we build and empirically calibrate a new attenuation model that takes the increased Ly$\alpha$ damping wing absorption into account. Our model consists of the canonical IGM absorption and an additional absorption component due to dense neutral hydrogen gas clouds proximate to the galaxy, and we derive the redshift evolution of HI column density of the proximate clouds by calibrating the model using CANUCS JWST observations. The resulting total transmission curve resolves the photometric redshift bias at $z>7$, an improvement that is robust to choice of template-fitting code, template set, and photometric catalog used. The new attenuation model can be easily implemented in existing template-fitting codes, and significantly improves the photometric redshift performance in the EoR., Comment: Submitted to ApJL. Comments are welcome. The new attenuation model has been implemented in the latest version EAzY (v0.8 or later) and Phosphoros (v2.0.2 or later)

Despite the ubiquity of clumpy star-forming galaxies at high-redshift, the origin of clumps are still largely unconstrained due to the limited observations that can validate the mechanisms for clump formation. We postulate that if clumps form due to the accretion of metal-poor gas that leads to violent disk instability, clumpy galaxies should have lower gas-phase metallicities compared to non-clumpy galaxies. In this work, we obtain the near-infrared spectrum for 42 clumpy and non-clumpy star-forming galaxies of similar masses, SFRs, and colors at $z\approx0.7$ using the Gemini Near-Infrared Spectrograph (GNIRS) and infer their gas-phase metallicity from the {\nii} and {\halpha} line ratio. We find that clumpy galaxies have lower metallicities compared to non-clumpy galaxies, with an offset in the weighted average metallicity of $0.07\pm0.02$ dex. We also find an offset of $0.06\pm0.02$ dex between clumpy and non-clumpy galaxies in a comparable sample of 23 star-forming galaxies at $z\approx1.5$ using existing data from the FMOS-COSMOS survey. Similarly, lower {\nii}/{\halpha} ratio are typically found in galaxies that have more of their $\mathrm{UV_{rest}}$ luminosity originating from clumps, suggesting that \enquote{clumpier} galaxies are more metal poor. We also derive the intrinsic velocity dispersion and line-of-sight rotational velocity for galaxies from the GNIRS sample. The majority of galaxies have $\sigma_0/v_c \approx 0.2$, with no significant difference between clumpy and non-clumpy galaxies. Our result indicates that clump formation may be related to the inflow of metal-poor gas; however, the process that forms them does not necessarily require significant, long-term kinematic instability in the disk., Comment: 19 pages, 10 figures, 1 table, 2nd revision submitted to ApJ following minor comments

We present galaxy MACS0416-Y1 at z$_{\rm{spec}} = 8.312$ as observed by the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). MACS0416-Y1 has been shown to have extreme dust properties, thus, we study the physical properties and star formation histories of its resolved components. Overall, we find that MACS0416-Y1 is undergoing a star formation burst in three resolved clumps. The central clump is less massive compared to the other clumps and possibly formed in the merging process of the two larger clumps. Although the star formation history indicates an ongoing star formation burst, this gas-rich galaxy shows comparable star formation efficiency to cosmic noon galaxies. Using NIRSpec prism spectroscopy, we measure metallicity, $12 +\log\rm{(O/H)} = 7.76\pm0.03$ , ionisation parameter, $\log U = -2.48\pm0.03$, and electron temperature $\rm{T}_e = 18000\pm 4000 K $. The emission line ratios of the galaxy indicate an evolved Interstellar medium (ISM) similar to $z\sim2$ star-forming galaxies. Further, we find possible presence of ionisation from an active galactic nuclei (AGN) using emission line diagnostics, however, we do not detect broad line component in H$\beta$ emission line. As this gas-rich galaxy is undergoing a major merger, we hypothesise that the high dust temperature in MACS0416-Y1 is caused by the star formation burst or a possible narrow-line AGN.

We present the design and observations of low resolution JWST/NIRSpec PRISM spectroscopy from the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 JWST Treasury program. Targets are selected using JWST/NIRCam photometry from UNCOVER and other programs, and cover a wide range of categories and redshifts to ensure the legacy value of the survey. These categories include the first galaxies at $z\gtrsim10$, faint galaxies during the Epoch of Reionization ($z\sim6-8$), high redshift AGN ($z\gtrsim6$), Population III star candidates, distant quiescent and dusty galaxies ($1\lesssim z \lesssim 6$), and filler galaxies sampling redshift--color--magnitude space from $z\sim 0.1-13$. Seven NIRSpec MSA masks across the extended Abell 2744 cluster were observed, along with NIRCam parallel imaging in 8 filters (F090W, F115W, F150W, F200W, F277W, F356W, F410M, F444W, F480M) over a total area of ~26 arcmin$^2$, overlapping existing HST coverage from programs including the Hubble Frontier Fields and BUFFALO. We successfully observed 553 objects down to $m_{\mathrm{F444W}}\sim30\mathrm{AB}$, and by leveraging mask overlaps, we reach total on-target exposure times ranging from 2.4-16.7h. We demonstrate the success rate and distribution of confirmed redshifts, and also highlight the rich information revealed by these ultradeep spectra for a subset of our targets. An updated lens model of Abell 2744 is also presented, including 14 additional spectroscopic redshifts and finding a total cluster mass of $M_{\mathrm{SL}}=(2.1\pm0.3)\times10^{15}\,\mathrm{M}_{\odot}$. We publicly release reduced 1D and 2D spectra for all objects observed in Summer 2023 along with a spectroscopic redshift catalog and the updated lens model of the cluster (https://jwst-uncover.github.io/DR4.html)., Comment: 19 pages, 9 figures, 4 tables, submitted to ApJ, comments welcome! Data available at: https://jwst-uncover.github.io/DR4.html (v2: figure format correction)

We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log $M/M_\odot=10.2-11.2$) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.26$\pm0.04$ dex in [C/H], 0.16$\pm0.03$ dex in [Fe/H], and 0.07$\pm0.04$ dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies forming at higher redshift consistently show higher [Mg/Fe] and lower [Fe/H] and [Mg/H], regardless of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily driven by lower-redshift samples naturally including galaxies with longer star-formation timescales. In contrast, the lower [Mg/H] likely reflects earlier-forming galaxies expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is somewhat lower at z=1-3 compared to the lower redshift relation. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies., Comment: Accepted to ApJ; 22 pages, 6 figures, 3 tables

The epoch of reionisation progressed through the emission of ionising photons from galaxies to their local intergalactic medium. In this work, we characterise the dwarf star-forming galaxies as candidates for the source of ionising photons that drove EoR. We investigate the ionising properties and star formation histories of star-forming dwarf galaxies at the last stages of EoR at $4.8<\rm{z}<7$ using observations from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The magnification due to gravitational lensing allows us to probe large dynamic ranges in stellar mass ($2\times 10^{6}\leq\rm{M}_*/\rm{M}_\odot\leq5\times 10^{9}$) and UV magnitudes ($-22.68\leq$M$_{UV}\leq=-15.95$).We find a median UV slope \buv of $-2. 56\pm0.23$ and the production efficiency of ionising photons $\log$ \xiion $=25.39\pm0.6$ for the full sample ($4.8<\rm{z}<7$) with a median stellar mass of $6.3\pm0.5\times10^{7} \rm{M}_\odot$. We find both \buv and \xiion are marginally correlated with the stellar mass of the galaxy, indicating a possible greater contribution of dwarf galaxies to the reionisation of the Universe. We find that on average, galaxies in our sample are experiencing a recent rise/burst of star formation which translates to a higher scatter in \xiion and a large scatter in H$\alpha$ equivalent widths. Finally, we investigate the trends of H$\alpha$ and [OIII]+H$\beta$ EWs with UV magnitude and find M$_{UV}$ is correlated between H$\alpha$ but not with [OIII]+H$\beta$ EWs indicating low metallicities and recent burst in the UV faint galaxies., Comment: Accepted For Publication in MNRAS

Spatially resolved studies are key to understanding when, where, and how stars form within galaxies. Using slitless grism spectra and broadband imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS) we study the spatially resolved properties of a strongly lensed ($\mu$ = 5.4$\pm$1.8) z = 0.8718 galaxy pair consisting of a blue face-on galaxy (10.2 $\pm$ 0.2 log($M/M_\odot$)) with multiple star-forming clumps and a dusty red edge-on galaxy (9.9 $\pm$ 0.3 log($M/M_\odot$)). We produce accurate H$\alpha$ maps from JWST/NIRISS grism data using a new methodology that accurately models spatially varying continuum and emission line strengths. With spatially resolved indicators, we probe star formation on timescales of $\sim$ 10 Myr (NIRISS H$\alpha$ emission line maps) and $\sim$ 100 Myr (UV imaging and broadband SED fits). Taking the ratio of the H$\alpha$ to UV flux ($\eta$), we measure spatially resolved star formation burstiness. We find that in the face-on galaxy both H$\alpha$ and broadband star formation rates (SFRs) drop at large galactocentric radii by a factor of $\sim$ 4.7 and 3.8 respectively, while SFR over the last $\sim$ 100 Myrs has increased by a factor of 1.6. Additionally, of the 20 clumps identified in the galaxy pair we find that 7 are experiencing bursty star formation, while 10 clumps are quenching and 3 are in equilibrium (either being in a state of steady star formation or post-burst). Our analysis reveals that the blue face-on galaxy disk is predominantly in a quenching or equilibrium phase. However, the most intense quenching within the galaxy is seen in the quenching clumps. This pilot study demonstrates what JWST/NIRISS data can reveal about spatially varying star formation in galaxies at Cosmic Noon., Comment: MNRAS Accepted, 15 pages, 12 figures

Using deep imaging from the CANDELS and HFF surveys, we present bulge+disc decompositions with GalfitM for $\sim$17,000 galaxies over $0.2 \leq z\leq 1.5$. We use various model parameters to select reliable samples of discs and bulges, and derive their stellar masses using an empirically calibrated relation between mass-to-light ratio and colour. Across our entire redshift range, we show that discs follow stellar mass-size relations that are consistent with those of star-forming galaxies, suggesting that discs primarily evolve via star formation. In contrast, the stellar mass-size relations of bulges are mass-independent. Our novel dataset further enables us to separate components into star-forming and quiescent based on their specific star formation rates. We find that both star-forming discs and star-forming bulges lie on stellar mass-size relations that are similar to those of star-forming galaxies, while quiescent discs are typically smaller than star-forming discs and lie on steeper relations, implying distinct evolutionary mechanisms. Similar to quiescent galaxies, quiescent bulges show a flattening in the stellar mass-size relation at $\sim$10$^{10}$M$_\odot$, below which they show little mass dependence. However, their best-fitting relations have lower normalisations, indicating that at a given mass, bulges are smaller than quiescent galaxies. Finally, we obtain rest-frame colours for individual components, showing that bulges typically have redder colours than discs, as expected. We visually derive UVJ criteria to separate star-forming and quiescent components and show that this separation agrees well with component colour. HFF bulge+disc decomposition catalogues used for these analyses are publicly released with this paper., Comment: 30 pages, 18 figures, and 6 tables. Resubmitted to MNRAS after addressing a thorough and constructive referee report

Strong gravitational lensing in galaxy clusters has become an essential tool in astrophysics, allowing us to directly probe the dark matter distribution and study magnified background sources. The precision and reliability of strong lensing models rely heavily on the number and quality of multiple images of background sources with spectroscopic redshifts. We present an updated strong lensing model of the galaxy cluster MACS J0416.1-2403 with the largest sample of multiple images with spectroscopic redshifts in a galaxy cluster field to date. Furthermore, we aim to demonstrate the effectiveness of JWST particularly its NIRISS camera, for strong lensing studies. We use the JWST 's NIRCam imaging and NIRSpec and NIRISS spectroscopy from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The cluster mass model is constrained using Lenstool software. Our new dataset, used for constraining the lens model, comprises 303 secure multiple images from 111 background sources and includes systems with previously known MUSE redshift and systems for which we obtained spectroscopic redshift for the first time using NIRISS and NIRSpec spectroscopy. The total number of secure spectroscopic systems is >20% higher than in the previous strong lensing studies of this cluster. The derived strong lensing model can reproduce multiple images with the root-mean-square distance of 0.53''. We also provide a full catalogue with 415 multiple images, including less reliable candidates. We furthermore demonstrate the effectiveness of JWST particularly NIRISS, for strong lensing studies. As NIRISS F115W, F150W, and F200W grism spectroscopy captures at least two of the [OII] {\lambda}3727, [OIII] {\lambda}{\lambda}4959, 5007, and H{\alpha} lines at 1

We present an analysis of the number density of galaxies as a function of stellar mass (i.e., the stellar mass function, SMF) in the COSMOS field at z~3.3, making a comparison between the SMF in overdense environments and the SMF in the coeval field. In particular, this region contains the Elent\'ari proto-supercluster, a system of 6 extended overdensities spanning ~70 cMpc on a side. A clear difference is seen in the high-mass slope of these SMFs, with overdense regions showing an increase in the ratio of high-mass galaxies to low-mass galaxies relative to the field, indicating a more rapid build-up of stellar mass in overdense environments. This result qualitatively agrees with analyses of clusters at z~1, though the differences between protocluster and field SMFs at z~3.3 are smaller. While this is consistent with overdensities enhancing the evolution of their member galaxies, potentially through increased merger rates, whether this enhancement begins in protocluster environments or even earlier in group environments is still unclear. Though the measured fractions of quiescent galaxies between the field and overdense environments do not vary significantly, implying that this stellar mass enhancement is ongoing and any starbursts triggered by merger activity have not yet quenched, we note that spectroscopic observations are biased towards star-forming populations, particularly for low-mass galaxies. If mergers are indeed responsible, high resolution imaging of Elent\'ari and similar structures at these early epochs should then reveal increased merger rates relative to the field. Larger samples of well-characterized overdensities are necessary to draw broader conclusions in these areas., Comment: 15 pages, 8 figures, plus references and appendices. Submitted to ApJ, posted after responding to referee comments. Further comments welcome

We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.03$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe., Comment: Accepted for publication in ApJ

We present spectroscopic confirmation of an ultra-massive galaxy (UMG) with $\log(M_\star/M_\odot) = 10.98 \pm 0.07$ at $z_\mathrm{spec} = 4.8947$ in the Extended Groth Strip (EGS), based on deep observations of Ly$\alpha$ emission with Keck/DEIMOS. The ultra-massive galaxy (UMG-28740) is the most massive member in one of the most significant overdensities in the EGS, with four additional photometric members with $\log(M_\star/M_\odot) > 10.5$ within $R_\mathrm{proj} \sim 1$ cMpc. The Ly$\alpha$ profile is highly asymmetric ($A_f = 3.56$), suggesting the presence of neutral gas within the interstellar medium, circumgalactic medium, or via AGN-driven outflows. Spectral energy distribution (SED) fitting using a large suite of star formation histories and two sets of high-quality photometry from ground- and space-based facilities consistently estimates the stellar mass of UMG-28740 to be $\log(M_\star/M_\odot) \sim 11$ with a small standard deviation between measurements ($\sigma = 0.07$). While the best-fit SED models agree on stellar mass, we find discrepancies in the estimated star formation rate for UMG-28740, resulting in either a star-forming or quiescent system. JWST/NIRCam photometry of UMG-28740 strongly favors a quiescent scenario, demonstrating the need for high-quality mid-IR observations. Assuming the galaxy to be quiescent, UMG-28740 formed the bulk of its stars at $z > 10$ and is quenching at $z \sim 8$, resulting in a high star formation efficiency at high redshift ($\epsilon \sim 0.2$ at $z \sim 5$ and $\epsilon \gtrsim 1$ at $z \gtrsim 8$). As the most massive galaxy in its protocluster environment, UMG-28740 is a unique example of the impossibly early galaxy problem., Comment: 11 pages, 6 figures, 2 tables, Accepted to OJA

In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ~17 arcmin^2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z>6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields., Comment: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ($z=1-3$, $H_{AB}\le23$), as well as 53 star-forming galaxies at $z=1-4$. With 16~hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log$M_*/M_{\odot}\sim10.2-11.5$) and detailed star-formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages ($0.8-3.3$~Gyr) and star formation timescales ($\sim0.5-4$~Gyr), and four out of 20 galaxies were already quiescent by $z=3$. On average, the $z>1.75$ $[z<1.75]$ galaxies formed 50\% of their stellar mass before $z=4$ $[z=3]$. Furthermore, the typical SFHs of galaxies in these two redshift bins ($z_{\text{mean}}=2.2~[1.3]$) indicate that galaxies at higher redshift formed earlier and over shorter star-formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, number density calculations imply that mergers and/or late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances, resolved stellar populations and kinematics., Comment: Accepted in ApJ; 25 pages, 14 figures, 2 tables (excluding appendices)

The changes in colors across a galaxy are intimately connected to the galaxy's formation, growth, quenching history, and dust content. A particularly important epoch in the growth of galaxies is near $z \sim 2$ often referred to as `cosmic noon', where galaxies on average reach the peak of their star formation. We study a population of 125 cluster galaxies at $z \sim 1.6$ in three Hubble Space Telescope (HST) filters, F475W, F625W, and F160W, roughly corresponding to the rest-frame FUV, NUV, and r band, respectively. By comparing to a control sample of 200 field galaxies at similar redshift, we reveal clear, statistically significant differences in the overall spatially resolved colors and color gradients in galaxies across these two different environments. On average, cluster galaxies have redder UV colors in both the inner and outer regions bounded by $r_{\mathrm{50}}$, as well as an overall wider dispersion of outside-in color gradients. The presence of these observed differences, along with evidence from ancillary data from previous studies, strongly suggests that the environment drives these population-level color differences, by affecting the stellar populations and/or dust content., Comment: Accepted for publication in ApJ

We report an updated mass and magnification model of galaxy cluster Abell 370 using new NIRCam and NIRISS data from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using Lenstool and a combination of archival HST and MUSE data with new JWST data as constraints, we derive an improved gravitational lensing model and extract magnifications of background galaxies with uncertainties. Using our best fit model, we perform a search for new multiply imaged systems via predicted positions. We report no new multiply imaged systems with identifiable redshifts, likely due to already very deep HST and Spitzer data, but confirm a $z\sim8$ multiply imaged system by measuring its redshift with NIRISS and NIRSpec spectra. We find that the overall shape of the critical curve for a source at $z = 9.0$ is similar to previous models of Abell 370, with small changes. We investigate the $z\sim8$ galaxy with two images observable with an apparent magnitude in the F125W band of $26.0\pm0.2$ and $25.6\pm0.1$. After correcting for the magnifications of the images, 7.2$^{+0.2}_{-1.2}$ and 8.7$^{+0.4}_{-0.4}$, we use SED fitting to find an intrinsic stellar mass of log($M^*/M_{\odot})$ = 7.35$^{+0.04}_{-0.05}$, intrinsic SFR of 3.5$^{+2.2}_{-1.4}$ M$_{\odot}$/yr, and $M_{UV}$ of -21.3$^{+0.2}_{-0.2}$, which is close to the knee of the luminosity function at that redshift. Our model, and corresponding magnification, shear, and convergence maps are available on request and will be made publicly available on MAST in a CANUCS data release (DOI: 10.17909/ph4n-6n76)., Comment: 15 pages

The resolved mass assembly of Milky-Way-mass galaxies has been previously studied in simulations, the local universe, and at higher redshifts using infrared (IR) light profiles. To better characterize the mass assembly of Milky Way Analogues (MWAs), as well as their changes in star-formation rate and color gradients, we construct resolved stellar mass and star-formation rate maps of MWA progenitors selected with abundance matching techniques up to z $\sim$ 2 using deep, multi-wavelength imaging data from the Hubble Frontier Fields. Our results using stellar mass profiles agree well with previous studies that utilize IR light profiles, showing that the inner 2 kpc of the galaxies and the regions beyond 2 kpc exhibit similar rates of stellar mass growth. This indicates the progenitors of MWAs from $z\sim 2$ to the present do not preferentially grow their bulges or their disks. The evolution of the star-formation rate (SFR) profiles indicate greater decrease in SFR density in the inner regions versus the outer regions. S\'ersic parameters indicate modest growth in the central regions at lower redshifts, perhaps indicating slight bulge growth. However, the S\'ersic index does not rise above $n \sim 2$ until $z < 0.5$, meaning these galaxies are still disk dominated systems. We find that the half-mass radii of the MWA progenitors increase between $1.5 < z < 2$, but remain constant at later epochs ($z < 1.5$). This implies mild bulge growth since $z\sim 2$ in MWA progenitors, in line with previous MWA mass assembly studies., Comment: 19 Pages, 12 figures, accepted by ApJ

The most distant galaxies detected by JWST are assembling in a Universe that is less than 5\% of its present age. At these times, the progenitors of galaxies like the Milky Way are expected to be about 10,000 times less massive than they are now, with masses quite comparable to that of massive globular clusters seen in the local Universe. Composed today primarily of old stars and correlating with the properties of their parent dark matter halos, the first globular clusters are thought to have formed during the earliest stages of galaxy assembly. In this article we explore the connection between star clusters and galaxy assembly by showing JWST observations of a strongly lensed galaxy at zspec = 8.304, exhibiting a network of massive star clusters (the 'Firefly Sparkle') cocooned in a diffuse arc. The Firefly Sparkle exhibits the hallmarks expected of a future Milky Way-type galaxy captured during its earliest and most gas-rich stage of formation. The mass distribution of the galaxy seems to be concentrated in ten distinct clusters, with individual cluster masses that straddle the boundary between low-mass galaxies and high-mass globular clusters. The cluster ages suggest that they are gravitationally bound with star formation histories showing a recent starburst possibly triggered by the interaction with a companion galaxy at the same redshift at a projected distance of $\sim$2 kpc away from the Firefly Sparkle. The central star cluster shows nebular-dominated spectra consistent with high temperatures and a top-heavy initial mass function, the product of formation in a very metal poor environment. Combined with abundance matching that suggests that this is likely to be a progenitor of galaxies like our own, the Firefly Sparkle provides an unprecedented case study of a Milky Way-like galaxy in the earliest stages of its assembly in only a 600 million year old Universe., Comment: submitted, comments are welcome! 36 pages, 9 figures

Many quiescent galaxies discovered in the early Universe by \textit{JWST} raise fundamental questions on when and how these galaxies became and stayed quenched. Making use of the latest version of the semi-analytic model GAEA that provides good agreement with the observed quenched fractions up to $z\sim 3$, we make predictions for the expected fractions of quiescent galaxies up to $z\sim 7$ and analyze the main quenching mechanism. We find that in a simulated box of $685~{\rm Mpc}$ on a side, the first quenched massive ($M_{\star} \sim 10^{11} {\rm M}_{\odot}$), Milky Way mass, and low mass ($M_{\star} \sim 10^{9.5} {\rm M}_{\odot}$ ) galaxies appear at $z\sim 4.5$, $z\sim 6.2$, and before $z = 7$. Most quenched galaxies identified at early redshifts remain quenched for more than 1 Gyr. Independently of galaxy stellar mass, the dominant quenching mechanism at high redshift is accretion disk feedback (quasar winds) from a central massive black hole, which is triggered by mergers in massive and MW-mass galaxies, and by disk instabilities in low-mass galaxies. Environmental stripping becomes increasingly more important at lower redshift., Comment: 8 pages, 3 figures, Accepted publication in APJL

We have identified 189 candidate $z > 1.3$ protoclusters and clusters in the LSST Deep Drilling Fields. This sample will enable the measurement of the metal enrichment and star formation history of clusters during their early assembly period through the direct measurement of the rate of supernovae identified through the LSST. The protocluster sample was selected from galaxy overdensities in a $Spitzer$/IRAC colour-selected sample using criteria that were optimised for protocluster purity using a realistic lightcone. Our tests reveal that $60-80\%$ of the identified candidates are likely to be genuine protoclusters or clusters, which is corroborated by a $\sim4\sigma$ stacked X-ray signal from these structures. We provide photometric redshift estimates for 47 candidates which exhibit strong peaks in the photo-$z$ distribution of their candidate members. However, the lack of a photo-$z$ peak does not mean a candidate is not genuine, since we find a stacked X-ray signal of similar significance from both the candidates that exhibit photo-$z$ peaks and those that do not. Tests on the lightcone reveal that our pursuit of a pure sample of protoclusters results in that sample being highly incomplete ($\sim4\%$) and heavily biased towards larger, richer, more massive, and more centrally concentrated protoclusters than the total protocluster population. Most ($\sim75\%$) of the selected protoclusters are likely to have a maximum collapsed halo mass of between $10^{13}-10^{14}$ M$_{\odot}$, with only $\sim25\%$ likely to be collapsed clusters above $10^{14}$ M$_{\odot}$. However, the aforementioned bias ensures our sample is $\sim50\%$ complete for structures that have already collapsed into clusters more massive than $10^{14}$ M$_{\odot}$., Comment: 18 pages, 14 figures

We present the construction of a deep multi-wavelength PSF-matched photometric catalog in the UDS field following the final UKIDSS UDS DR11 release. The catalog includes photometry in 24 filters, from the MegaCam-uS (0.38 microns) band to the Spitzer-IRAC (8 microns) band, over ~ 0.9 sq. deg. and with a 5-sigma depth of 25.3 AB in the K-band detection image. The catalog, containing approximately 188,564 (136,235) galaxies at 0.2 < z < 8.0 with stellar mass > 10$^{8}$ solar masses and K-band total magnitude K < 25.2 (24.3) AB, enables a range of extragalactic studies. We also provide photometric redshifts, corresponding redshift probability distributions, and rest-frame absolute magnitudes and colors derived using the template-fitting code eazy-py. Photometric redshift errors are less than 3 to 4 percent at z < 4 across the full brightness range in K-band and stellar mass range 10$^{8}$-10$^{12}$ solar masses. Stellar population properties (e.g., stellar mass, star-formation rate, dust extinction) are derived from the modeling of the spectral energy distributions (SEDs) using the codes FAST and Dense Basis., Comment: 26 pages, 18 figures, accepted in ApJ. Photometric Catalog, catalogs of stellar population properties derived using Eazy and Dense Basis and other high-level products available for public use at the following webpage - https://www.zaidikumail.com/feniks-uds-catalogs (DOI:10.5281/zenodo.11002298)

We present a methodology for modeling and removing light from cluster galaxies and intracluster light (ICL) from $James\ Webb\ Space\ Telescope$ ($JWST$) images of gravitational lensing clusters. We apply our method to Webb's First Deep Field the SMACS 0723 Early Release Observations and use the ICL subtracted images to select a sample of globular clusters (GCs) and dwarf galaxies within the cluster. We compare the spatial distributions of these two samples with our models of the galaxy and ICL light, finding significant similarity. Specifically we find that GCs trace the diffuse ICL, while dwarf galaxies are centrally concentrated near the cluster center We quantify the relationship between the surface density of compact sources and total cluster light, demonstrating a significant, tight correlation. We repeat our methodology and compare distributions of GCs with dark matter surface density and find a comparable result. Our findings suggest a common origin for GCs and diffuse ICL, with stripping from massive galaxies as they merge with the cluster being a plausible scenario., Comment: 16 pages, 13 figures. Submitted to ApJ. Comments welcome

We present an analysis of the galaxy stellar mass function (SMF) of 14 known protoclusters between $2.0 < z < 2.5$ in the COSMOS field, down to a mass limit of $10^{9.5}$ M$_{\odot}$. We use existing photometric redshifts with a statistical background subtraction, and consider star-forming and quiescent galaxies identified from $(NUV - r)$ and $(r - J)$ colours separately. Our fiducial sample includes galaxies within 1 Mpc of the cluster centres. The shape of the protocluster SMF of star-forming galaxies is indistinguishable from that of the general field at this redshift. Quiescent galaxies, however, show a flatter SMF than in the field, with an upturn at low mass, though this is only significant at $\sim 2\sigma$. There is no strong evidence for a dominant population of quiescent galaxies at any mass, with a fraction of $< 15\%$ at $1\sigma$ confidence for galaxies with log$M_{\ast}/M_{\odot} < 10.5$. We compare our results with a sample of galaxies groups at $1 < z < 1.5$, and demonstrate that a significant amount of environmental quenching must take place between these epochs, increasing the relative abundance of high-mass ($\rm M > 10^{10.5} M_{\odot}$) quiescent galaxies by a factor of $\gtrsim$ 2. However, we find that at lower masses ($\rm M < 10^{10.5} M_{\odot}$), no additional environmental quenching is required., Comment: 23 pages, 22 figures. Accepted for publication in MNRAS

We present a new sample of 158 galaxies at redshift $z>7.5$ selected from deep \jwst\ NIRCam imaging of five widely-separated sightlines in the CANUCS survey. Two-thirds of the pointings and 80\% of the galaxies are covered by 12 to 14 NIRCam filters, including seven to nine medium bands, providing accurate photometric redshifts and robustness against low redshift interlopers. A sample of 28 galaxies at $z>7.5$ with spectroscopic redshifts shows a low systematic offset and scatter in the difference between photometric and spectroscopic redshifts. We derive the galaxy UV luminosity function at redshifts 8 to 12, finding a slightly higher normalization than previously seen with \hst\ at redshifts 8 to 10. We observe a steeper decline in the galaxy space density from $z=8$ to $12$ than found by most \jwst\ Cycle 1 studies. In particular, we find only eight galaxies at $z>10$ and none at $z>12.5$, with no $z>10$ galaxies brighter than F277W AB=28 or $M_{\rm UV}=-20$ in our unmasked, delensed survey area of 53.4 square arcminutes. We attribute the lack of bright $z>10$ galaxies in CANUCS compared to GLASS and CEERS to intrinsic variance in the galaxy density along different sightlines. The evolution in the CANUCS luminosity function between $z=8$ and $12$ is comparable to that predicted by simulations that assume a standard star formation efficiency, without invoking any special adjustments., Comment: 30 pages, 12 figures, ApJ, in press

Early JWST observations that targeted so-called double-break sources (attributed to Lyman and Balmer breaks at $z>7$), reported a previously unknown population of very massive, evolved high-redshift galaxies. This surprising discovery led to a flurry of attempts to explain these objects' unexpected existence including invoking alternatives to the standard $\Lambda$CDM cosmological paradigm. To test these early results, we adopted the same double-break candidate galaxy selection criteria to search for such objects in the JWST images of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS), and found a sample of 19 sources over five independent CANUCS fields that cover a total effective area of $\sim60\,$arcmin$^2$ at $z\sim8$. However, (1) our SED fits do not yield exceptionally high stellar masses for our candidates, while (2) spectroscopy of five of the candidates shows that while all five are at high redshifts, their red colours are due to high-EW emission lines in star-forming galaxies rather than Balmer breaks in massive, evolved systems. Additionally, (3) field-to-field variance leads to differences of $\sim 1.5$ dex in the maximum stellar masses measured in the different fields, suggesting that the early single-field JWST observations may have suffered from cosmic variance and/or sample bias. Finally, (4) we show that the presence of even a single massive outlier can dominate conclusions from small samples such as those in early JWST observations. In conclusion, we find that the double-break sources in CANUCS are not sufficiently massive or numerous to warrant questioning the standard $\Lambda$CDM paradigm., Comment: V3 Accepted for publication in MNRAS, 13 pages (+5 in Appendix), 7 figures (+4), 1 table (+2)

We use CANUCS JWST/NIRCam imaging of galaxies behind the gravitationally-lensing cluster MACS J0417.5-1154 to investigate star formation burstiness in low-mass ($M_\star\sim10^8\ M_\odot$) galaxies at $z\sim4.7-6.5$. Our sample of 123 galaxies is selected using the Lyman break selection and photometric emission-line excess methods. Sixty per cent of the 123 galaxies in this sample have H$\alpha$-to-UV flux ratios that deviate significantly from the range of $\eta_{1500}$ values consistent with smooth and steady star formation histories. This large fraction indicates that the majority of low-mass galaxies is experiencing bursty star formation histories at high redshift. We also searched for interacting galaxies in our sample and found that they are remarkably common ($\sim40$ per cent of the sample). Compared to non-interacting galaxies, interacting galaxies are more likely to have very low H$\alpha$-to-UV ratios, suggesting that galaxy-galaxy interactions enhance star formation burstiness and enable faster quenching (with timescales of $\lesssim100$ Myr) that follows the rapid rise of star formation activity. Given the high frequency of galaxy-galaxy interactions and the rapid SFR fluctuations they appear to cause, we conclude that galaxy-galaxy interactions could be a leading cause of bursty star formation in low-mass, high-$z$ galaxies. They could thus play a significant role in the evolution of the galaxy population at early cosmological times., Comment: 20 pages, 15 figures, and 1 table. Accepted for publication in MNRAS

The recent UNCOVER survey with the James Webb Space Telescope (JWST) exploits the nearby cluster Abell 2744 to create the deepest view of our universe to date by leveraging strong gravitational lensing. In this work, we perform photometric fitting of more than 50,000 robustly detected sources out to $z \sim 15$. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of $0.2 \lesssim z \lesssim 15$. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector-$\beta$ priors on masses and star formation histories to produce joint redshift and stellar population posteriors, and additionally lensing magnification is performed on-the-fly to ensure consistency with the scale-dependent priors. We show that this approach produces excellent photometric redshifts with $\sigma_{\rm NMAD} \sim 0.03$, of a similar quality to the established photometric redshift code EAzY. In line with the open-source scientific objective of the Treasury survey, we publicly release the stellar population catalog with this paper, derived from the photometric catalog adapting aperture sizes based on source profiles. This release includes posterior moments, maximum-likelihood spectra, star-formation histories, and full posterior distributions, offering a rich data set to explore the processes governing galaxy formation and evolution over a parameter space now accessible by JWST., Comment: Corrected typos: Eq.1 should've been (1-kappa)^2, and the lens maps are normalized to D_ds/D_s=1. These errors were only in the writing; no data products or results were affected. The SPS catalogs are accessible via the UNCOVER survey webpage: https://jwst-uncover.github.io/DR2.html#SPSCatalogs, with a copy deposited to Zenodo: https://doi.org/10.5281/zenodo.8401181

We present measurements of stellar populations properties of a z = 9.1 gravitationally lensed galaxy MACS1149-JD1 using deep JWST NIRISS slitless spectroscopy as well as NIRISS and NIRCam imaging from the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). The galaxy is split into four components. Three magnified (${\mu}$ ~ 17) star-forming components are unresolved, giving intrinsic sizes of < 50pc. In addition, the underlying extended component contains the bulk of the stellar mass, formed the majority of its stars ~ 50Myr earlier than the other three components and is not the site of the most active star formation currently. The NIRISS and NIRCam resolved photometry does not confirm a strong Balmer break previously seen in Spitzer. The NIRISS grism spectrum has been extracted for the entire galaxy and shows a clear continuum and Lyman-break, with no Lyman-${\alpha}$ detected., Comment: This is the version of the article before peer review or editing, as submitted by an author to ApJLetters

The identification of sources driving cosmic reionization, a major phase transition from neutral Hydrogen to ionized plasma around 600-800 Myr after the Big Bang (Dayal et al. 2018, Mason et al. 2019, Robertson et al. 2022), has been a matter of intense debate (Robertson et al. 2022). Some models suggest that high ionizing emissivity and escape fractions ($f_{\rm esc}$) from quasars support their role in driving cosmic reionization (Madau & Haardt 2015, Mitra et al. 2018). Others propose that the high $f_{\rm esc}$ values from bright galaxies generates sufficient ionizing radiation to drive this process (Naidu et al. 2020). Finally, a few studies suggest that the number density of faint galaxies, when combined with a stellar-mass-dependent model of ionizing efficiency and $f_{\rm esc}$, can effectively dominate cosmic reionization (Finkelstein et al. 2019, Dayal et al. 2020). However, so far, low-mass galaxies have eluded comprehensive spectroscopic studies owing to their extreme faintness. Here we report an analysis of eight ultra-faint galaxies (in a very small field) during the epoch of reionization with absolute magnitudes between $M_{\rm UV}$ $\sim -17$ to $-15$ mag (down to 0.005 $L^{\star}$. We find that faint galaxies during the Universe's first billion years produce ionizing photons with log($\xi_{\mathrm{ion}}$/ Hz erg$^{-1}$) =$25.80\pm 0.14$, a factor of 4 higher than commonly assumed values (Robertson et al. 2015). If this field is representative of the large scale distribution of faint galaxies, the rate of ionizing photons exceeds that needed for reionization, even for escape fractions of order five per cent., Comment: 29 pages, 7 figures, 2 tables. Published in Nature

Motivated by spectroscopic confirmation of three overdense regions in the COSMOS field at $z\sim3.35$, we analyze the uniquely deep multi-wavelength photometry and extensive spectroscopy available in the field to identify any further related structure. We construct a three dimensional density map using the Voronoi tesselation Monte Carlo method and find additional regions of significant overdensity. Here we present and examine a set of six overdense structures at $3.205\sigma$ overdensity separated by up to 70 cMpc, suggestive of a proto-supercluster similar to the Hyperion system at $z\sim2.45$. Upcoming photometric surveys with JWST such as COSMOS-Web, and further spectroscopic follow-up will enable more extensive analysis of the evolutionary effects that such an environment may have on its component galaxies at these early times., Comment: Accepted for publication in MNRAS. 7 pages, 3 figures, 3 tables

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

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe., Comment: 28 pages, 12 figures, 4 tables, Accepted to ApJ

We use CEERS JWST/NIRCam imaging to measure rest-frame near-IR light profiles of $>$500 $M_\star>10^{10}~M_\odot$ galaxies in the redshift range $0.510^{11}~M_\odot$) quiescent galaxies between $z=0.5$ and $z=2.3$, again comparing $R_{\rm{opt}}$, \rmass~and \rmassd. We conclude that the main tenets of the size evolution narrative established over the past 20 years, based on rest-frame optical light profile analysis, still hold in the era of JWST/NIRCam observations in the rest-frame near-IR., Comment: Submitted to ApJ. Comments welcome

Many $z=1.5$ galaxies with a stellar mass ($M_{\star}$) $\geq 10^{10}\,\mathrm{M}_\odot$ are already quenched in both galaxy clusters ($>50$ per cent) and the field ($>20$ per cent), with clusters having a higher quenched fraction at all stellar masses compared to the field. A puzzling issue is that these massive quenched galaxies have stellar populations of similar age in both clusters and the field. This suggests that, despite the higher quenched fraction in clusters, the dominant quenching mechanism for massive galaxies is similar in both environments. In this work, we use data from the cosmological hydrodynamic simulations Hydrangea and EAGLE to test whether the excess quenched fraction of massive galaxies in $z = 1.5$ clusters results from fundamental differences in their halo properties compared to the field. We find that (i) at $10^{10} \leq$ $M_{\star}/\,\mathrm{M}_\odot\leq 10^{11}$, quenched fractions in the redshift range $1.5 < z < 3.5$ are consistently higher for galaxies with higher peak maximum circular velocity of the dark matter halo ($v_{\mathrm{max, peak}}$), and (ii) the distribution of $v_{\mathrm{max, peak}}$ is strongly biased towards higher values for cluster satellites compared to the field centrals. Due to this difference in the halo properties of cluster and field galaxies, secular processes alone may account for (most of) the environmental excess of massive quenched galaxies in high-redshift (proto) clusters. Taken at face value, our results challenge a fundamental assumption of popular quenching models, namely that clusters are assembled from an unbiased subset of infalling field galaxies. If confirmed, this would imply that such models must necessarily fail at high redshift, as indicated by recent observations., Comment: 11 pages, 4+2 figures (including appendices), Accepted for publication in MNRAS

We explore models of massive ($\gt 10^{10}~{\rm M}_{\odot}$) satellite quenching in massive clusters at $z\gtrsim1$ using an MCMC framework, focusing on two primary parameters: $R_{\rm quench}$ (the host-centric radius at which quenching begins) and $\tau_{\rm quench}$ (the timescale upon which a satellite quenches after crossing $R_{\rm quench}$). Our MCMC analysis shows two local maxima in the 1D posterior probability distribution of $R_{\rm quench}$ at approximately $0.25$ and $1.0~R_{\rm{200}}$. Analyzing four distinct solutions in the $\tau_{\rm quench}$-$R_{\rm quench}$ parameter space, nearly all of which yield quiescent fractions consistent with observational data from the GOGREEN survey, we investigate whether these solutions represent distinct quenching pathways and find that they can be separated between \textquote{starvation} and \textquote{core quenching} scenarios. The starvation pathway is characterized by quenching timescales that are roughly consistent with the total cold gas (H$_{2}$+H{\scriptsize I}) depletion timescale at intermediate $z$, while core quenching is characterized by satellites with relatively high line-of-sight velocities that quench on short timescales ($\sim 0.25$ Gyr) after reaching the inner region of the cluster ($\lt 0.30~R_{\rm{200}}$). Lastly, we break the degeneracy between these solutions by comparing the observed properties of transition galaxies from the GOGREEN survey. We conclude that only the \textquote{starvation} pathway is consistent with the projected phase-space distribution and relative abundance of transition galaxies at $z \sim 1$. However, we acknowledge that ram pressure might contribute as a secondary quenching mechanism., Comment: 15 pages; 8 figures; Accepted for publication in Monthly Notices of the Royal Astronomical Society

We use JWST/NIRCam medium band photometry in a single pointing of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS) to identify 118 Extreme Emission Line Galaxies (EELGs) over $1.7 \lesssim z \lesssim 6.7$, selected using a set of color cuts that target galaxies with extreme $\text{[OIII] + H}\beta$ and H$\alpha$ emission. We show that our medium band color selections are able to select galaxies based on emission line equivalent width (EW), which is advantageous to more commonly used selections since it does not require strong continuum emission, and can select galaxies with faint or red continuum fluxes. The median EWs of our sample is $EW(\text{H}\alpha) = 893 $ \AA\ and $ EW(\text{[OIII] + H}\beta) = 1255 $ \AA, and includes some objects with $EW(\text{[OIII] + H}\beta) \sim 3000$ \AA. These systems are mostly compact with low stellar mass (median $\log(M_\star/M_\odot) = 8.03$), low metallicity (median $Z = 0.14 Z_\odot$), little dust (median $A_V = 0.18$ mag) and high SSFR (median $SSFR = 1.18 \times 10^{-8}/yr$). Additionally, galaxies in our sample show increasing EW(\Ha) and EW(\OIIIHb) with redshift, an anti-correlation of EW(\Ha) with stellar mass, and no correlation between EW(\OIIIHb) and stellar mass. Finally, we present NIRSpec spectroscopy of 15 of the EELGs in our sample. These spectra confirm the redshifts and EWs of the EELGs calculated from the medium bands, which demonstrates the accuracy and efficiency of our color selections. Overall, we show that there are significant advantages to using medium band photometry to identify and study EELGs at a wide range of redshifts., Comment: Submitted to ApJL. 11 pages, 5 figures, animated Figure 1 is available to view here: https://www.youtube.com/watch?v=bTjjjyFRJys and available for download here: https://github.com/sunnawithers/EELG_colours/blob/main/Figure1_EELGcolours_Withers%2B2023.mp4

We present the first results on the spatial distribution of dust attenuation at $1.0

We report CFHT/SITELLE imaging Fourier Transform Spectrograph observations of the Brightest Cluster Galaxy (BCG) of galaxy cluster Abell 2390 at z=0.228. The BCG displays a prominent cone of emission in H$\alpha$, H$\beta$, [NII], and [OII] to the North-West with PA = 42$^o$, 4.4 arcsec in length (15.9 kpc), which is associated with elongated and asymmetric Chandra soft X-ray emission. The H$\alpha$ flux map also contains a "hook" of H$\alpha$ and [NII] emission resulting in a broadened northern edge to the cone. Using SITELLE/LUCI software we extract emission line flux, velocity, velocity dispersion, and continuum maps, and utilize them to derive flux ratio maps to determine ionization mechanisms and dynamical information in the BCG's emission line region. The Baldwin-Phillips-Terlevich diagnostics on the BCG cone indicate a composite ionization origin of photoionization due to star formation and shock. Strong LINER-like emission is seen in the nuclear region which hosts an AGN. As Abell 2390 is a cool-core cluster, we suggest that the cooling flow is falling onto the central BCG and interacting with the central AGN. The AGN produces jets that inflate "bubbles" of plasma in the ICM, as is often observed in local galaxy clusters. Furthermore, combining signs of AGN activities from radio, optical emission line and X-ray data over a large range of physical scale, we find evidence for three possible episodes of AGN activity in different epochs associated with the Abell 2390 BCG., Comment: 15 pages, 13 figures, accepted by MNRAS

We report the discovery of a low-mass $z=5.200\pm 0.002$ galaxy that is in the process of ceasing its star formation. The galaxy, MACS0417-z5PSB, is multiply imaged with magnification factors $\sim40$ by the galaxy cluster MACS J0417.5-1154, observed as part of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS). Using observations of MACS0417-z5PSB with a JWST/NIRSpec Prism spectrum and NIRCam imaging, we investigate the mechanism responsible for the cessation of star formation of the galaxy, and speculate about possibilities for its future. Using spectrophotometric fitting, we find a remarkably low stellar mass of $\rm{M_*}=4.3\pm^{0.9}_{0.8} \times 10^{7} \rm{M_{\odot}}$, less than 1% of the characteristic stellar mass at $z\sim5$. We measure a de-lensed rest-UV half-light radius in the source plane of $30\pm^{7}_{5}$ pc, and measure a star formation rate from H$\alpha$ of $0.14\pm^{0.17}_{0.12}$ $\rm{M_{\odot}/yr}$. We find that under the assumption of a double power law star formation history, MACS0417-z5PSB has seen a recent rise in star formation, peaking $\sim10-30$ Myr ago and declining precipitously since then. Together, these measurements reveal a low-mass, extremely compact galaxy which is in the process of ceasing star formation. We investigate the possibilities of mechanisms that have led to the cessation of star formation in MACS0417-z5PSB, considering stellar and AGN feedback, and environmental processes. We can likely rule out an AGN and most environmental processes, but leave open the possibility that MACS0417-z5PSB could be a star forming galaxy in the lull of a bursty star formation history., Comment: 9 pages, 4 figures, 1 table

In November 2022, the James Webb Space Telescope (JWST) returned deep near-infrared images of Abell 2744 -- a powerful lensing cluster capable of magnifying distant, incipient galaxies beyond it. Together with the existing Hubble Space Telescope (HST) imaging, this publicly available dataset opens a fundamentally new discovery space to understand the remaining mysteries of the formation and evolution of galaxies across cosmic time. In this work, we detect and measure some 60,000 objects across the 49 arcmin$^2$ JWST footprint down to a $5\,\sigma$ limiting magnitude of $\sim$30 mag in 0.32" apertures. Photometry is performed using circular apertures on images matched to the point spread function of the reddest NIRCam broad band, F444W, and cleaned of bright cluster galaxies and the related intra-cluster light. To give an impression of the photometric performance, we measure photometric redshifts and achieve a $\sigma_{\rm NMAD}\approx0.03$ based on known, but relatively small, spectroscopic samples. With this paper, we publicly release our HST and JWST PSF-matched photometric catalog with optimally assigned aperture sizes for easy use, along with single aperture catalogs, photometric redshifts, rest-frame colors, and individual magnification estimates. These catalogs will set the stage for efficient and deep spectroscopic follow-up of some of the first JWST-selected samples in Summer 2023., Comment: 28 pages, 19 figures, resubmitted to ApJS following significant data product improvements. Comments welcome. Catalogs can be accessed at https://jwst-uncover.github.io/DR2.html#PhotometricCatalogs

Using CANUCS imaging we found an apparent major merger of two $z\sim5$ ultra-low-mass galaxies ($M_\star\sim10^{7}M_\odot$ each) that are doubly imaged and magnified $\sim$12-15$\times$ by the lensing cluster MACS 0417. Both galaxies are experiencing young ($\sim$100 Myr), synchronised bursts of star formation with $\log({\rm sSFR/Gyr^{-1}} )\sim$1.3-1.4, yet SFRs of just $\sim 0.2 M_\odot\ {\rm yr}^{-1}$. They have sub-solar ($Z\sim0.2Z_\odot$) gas-phase metallicities and are connected by an even more metal-poor star-forming bridge. The galaxy that forms from the merger will have a mass of at least $M_\star\sim 2\times10^7 M_\odot$, at least half of it formed during the interaction-induced starburst. More than half of the ionizing photons produced by the system (before and during the merger) will have been produced during the burst. This system provides the first detailed look at a merger involving two high-$z$ ultra-low-mass galaxies of the type believed to be responsible for reionizing the Universe. It suggests that such galaxies can grow via a combination of mass obtained through major mergers, merger-triggered starbursts, and long-term in-situ star formation. If such high-$z$ mergers are common, then merger-triggered starbursts could be significant contributors to the ionizing photon budget of the Universe., Comment: 6 pages, 2 figures, and 1 table; accepted for publication in MNRAS Letter

We present a spectroscopic redshift catalogue of the SMACS J0723.3$-$7327 field ("Webb's First Deep Field") obtained from JWST/NIRISS grism spectroscopy and supplemented with JWST/NIRSpec and VLT/MUSE redshifts. The catalogue contains a total of 190 sources with secure spectroscopic redshifts, including 156 NIRISS grism redshifts, 123 of which are for sources whose redshifts were previously unknown. These new grism redshifts are secured with two or more spectroscopic features (64 sources), or with a single spectral feature whose identity is secured from the object's nine-band photometric redshift (59 sources). These are complemented with 17 NIRSpec and 48 MUSE redshifts, including six new NIRSpec redshifts identified in this work. In addition to the $z_{\rm cl}=0.39$ cluster galaxy redshifts (for which we provide $\sim$40 new NIRISS absorption-line redshifts), we also find three prominent galaxy overdensities at higher redshifts - at $z=1.1$, $z=1.4$, and $z=2.0$ - that were until now not seen in the JWST/NIRSpec and VLT/MUSE data. The paper describes the characteristics of our spectroscopic redshift sample and the methodology we have employed to obtain it. Our redshift catalogue is made available to the community at https://niriss.github.io/smacs0723., Comment: 19 pages, 13 figures, 3 appendices. Accepted for publication in MNRAS

In this paper we describe the survey design for the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) Cycle 1 \JWST Treasury program, which executed its early imaging component in November 2022. The UNCOVER survey includes ultradeep ($\sim29-30\mathrm{AB}$) imaging of $\sim$45 arcmin$^2$ on and around the well-studied Abell 2744 galaxy cluster at $z=0.308$ and will follow-up ${\sim}500$ galaxies with extremely deep low-resolution spectroscopy with the NIRSpec/PRISM during the summer of 2023, with repeat visits in summer 2024. We describe the science goals, survey design, target selection, and planned data releases. We also present and characterize the depths of the first NIRCam imaging mosaic, highlighting previously unparalleled resolved and ultradeep 2-4 micron imaging of known objects in the field. The UNCOVER primary NIRCam mosaic spans 28.8 arcmin$^2$ in seven filters (F115W, F150W, F200W, F277W, F356W, F410M, F444W) and 16.8 arcmin$^2$ in our NIRISS parallel (F115W, F150W, F200W, F356W, and F444W). To maximize early community use of the Treasury data set, we publicly release full reduced mosaics of public JWST imaging including 45 arcmin$^2$ NIRCam and 17 arcmin$^2$ NIRISS mosaics on and around the Abell 2744 cluster, including the Hubble Frontier Field primary and parallel footprints., Comment: 19 pages, 8 figures, 4 tables, Accepted for publication in ApJ

We present 65 extremely dust-obscured galaxies from the UltraVISTA DR3 survey of the COSMOS field at $1

Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension., Comment: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions