Search

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., Comment: Submitted to ApJL

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)

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