Search

We explore the physical properties of five massive quiescent galaxies at $z\sim2.5$, revealing the presence of non-negligible dust reservoirs. JWST NIRSpec observations were obtained for each target, finding no significant line emission; multiple star formation tracers independently place upper limits between $0.1-10~M_\odot / \mathrm{yr}$. Spectral energy distribution modeling with Prospector infers stellar masses between $\log_{10}[M / M_\odot] \sim 10-11$ and stellar mass-weighted ages between $1-2$ Gyr. The inferred mass-weighted effective radii ($r_{eff}\sim 0.4-1.4$ kpc) and inner $1$ kpc stellar surface densities ($\log_{10}[\Sigma / M_\odot \mathrm{kpc}^2 ]\gtrsim 9$) are typical of quiescent galaxies at $z \gtrsim 2$. The galaxies display negative color gradients (redder core and bluer outskirts); for one galaxy, this effect results from a dusty core, while for the others it may be evidence of an "inside-out" growth process. Unlike local quiescent galaxies, we identify significant reddening in these typical cosmic noon passive galaxies; all but one require $A_V \gtrsim 0.4$. This finding is in qualitative agreement with previous studies but our deep 20-band NIRCam imaging is able to significantly suppress the dust-age degeneracy and confidently determine that these galaxies are reddened. We speculate about the physical effects that may drive the decline in dust content in quiescent galaxies over cosmic time., Comment: 18 pages, 8 figures, submitted to ApJ

Hubble Space Telescope imaging shows that most star-forming galaxies at cosmic noon -- the peak of cosmic star formation history -- appear disk-dominated, leaving the origin of the dense cores in their quiescent descendants unclear. With the James Webb Space Telescope's (JWST) high-resolution imaging to 5 {\mu}m, we can now map the rest-frame near-infrared emission, a much closer proxy for stellar mass distribution, in these massive galaxies. We selected 70 star-forming galaxies with 10$<$log(M)$<$12 and 1.5$<$z$<$3 in the CEERS survey and compare their morphologies in the rest-frame optical to those in the rest-frame near-IR. While the bulk of these galaxies are disk-dominated in 1.5 {\mu}m (rest-frame optical) imaging, they appear more bulge-dominated at 4.4 {\mu}m (rest-frame near-infrared). Our analysis reveals that in massive star-forming galaxies at z$\sim$2, the radial surface brightness profiles steepen significantly, from a slope of $\sim$0.3/dex at 1.5 {\mu}m to $\sim$1.4/dex at 4.4 {\mu}m within radii $<$ 1 kpc. Additionally, we find their total flux contained within the central 1 kpc is approximately 7 times higher in F444W than in F150W. In rest-optical emission, a galaxy's central surface density appears to be the strongest indicator of whether it is quenched or star-forming. Our most significant finding is that at redder wavelengths, the central surface density ratio between quiescent and star-forming galaxies dramatically decreases from $\sim$10 to $\sim$1. This suggests the high central densities associated with galaxy quenching are already in place during the star-forming phase, imposing new constraints on the transition from star formation to quiescence., Comment: 9 pages, 5 figures, resubmitted to ApJ for publication

We present \textsc{Mathpop}, a novel method to infer the globular cluster (GC) counts in ultra-diffuse galaxies (UDGs) and low-surface brightness galaxies (LSBGs). Many known UDGs have a surprisingly high ratio of GC number to surface brightness. However, standard methods to infer GC counts in UDGs face various challenges, such as photometric measurement uncertainties, GC membership uncertainties, and assumptions about the GC luminosity functions (GCLFs). \textsc{Mathpop} tackles these challenges using the mark-dependent thinned point process, enabling joint inference of the spatial and magnitude distributions of GCs. In doing so, \textsc{Mathpop} allows us to infer and quantify the uncertainties in both GC counts and GCLFs with minimal assumptions. As a precursor to \textsc{Mathpop}, we also address the data uncertainties coming from the selection process of GC candidates: we obtain probabilistic GC candidates instead of the traditional binary classification based on the color--magnitude diagram. We apply \textsc{Mathpop} to 40 LSBGs in the Perseus cluster using GC catalogs from a \textit{Hubble Space Telescope} imaging program. We then compare our results to those from an independent study using the standard method. We further calibrate and validate our approach through extensive simulations. Our approach reveals two LSBGs having GCLF turnover points much brighter than the canonical value with Bayes' factor being $\sim4.5$ and $\sim2.5$, respectively. An additional crude maximum-likelihood estimation shows that their GCLF TO points are approximately $0.9$~mag and $1.1$~mag brighter than the canonical value, with $p$-value $\sim 10^{-8}$ and $\sim 10^{-5}$, respectively., Comment: 8 figures, 5 tables; submitted to ApJ, comments are welcomed

We report the spectroscopic discovery of a massive quiescent galaxy at $z_{\rm spec}=7.29\pm0.01$, just $\sim700\,$Myr after the Big Bang. RUBIES-UDS-QG-z7 was selected from public JWST/NIRCam and MIRI imaging from the PRIMER survey and observed with JWST/NIRSpec as part of RUBIES. The NIRSpec/PRISM spectrum reveals one of the strongest Balmer breaks observed thus far at $z>6$, no emission lines, but tentative Balmer and Ca absorption features, as well as a Lyman break. Simultaneous modeling of the NIRSpec/PRISM spectrum and NIRCam and MIRI photometry (spanning $0.9-18\,\mu$m) shows that the galaxy formed a stellar mass of log$(M_*/M_\odot)=10.23^{+0.04}_{-0.04}$ in a rapid $\sim 100-200\,$Myr burst of star formation at $z\sim8-9$, and ceased forming stars by $z\sim8$ resulting in $\log \rm{sSFR/yr}^{-1}<-10$. We measure a small physical size of $209_{-24}^{+33}\,{\rm pc}$, which implies a high stellar mass surface density within the effective radius of $\log(\Sigma_{*,\rm e}/{\rm M_\odot\,kpc}^{-2})=10.85_{-0.12}^{+0.11}$ comparable to the densities measured in quiescent galaxies at $z\sim2-5$. The 3D stellar mass density profile of RUBIES-UDS-QG-z7 is remarkably similar to the central densities of local massive ellipticals, suggesting that at least some of their cores may have already been in place at $z>7$. The discovery of RUBIES-UDS-QG-z7 has strong implications for galaxy formation models: the estimated number density of quiescent galaxies at $z\sim7$ is $>100\times$ larger than predicted from any model to date, indicating that quiescent galaxies have formed earlier than previously expected., Comment: 21 pages, 7 figures, submitted to ApJ

Early JWST studies found an apparent population of massive, compact galaxies at redshifts $z\gtrsim7$. Recently three of these galaxies were shown to have prominent Balmer breaks, demonstrating that their light at $\lambda_{\rm rest} \sim 3500$ $\r{A}$ is dominated by a stellar population that is relatively old ($\sim$200 Myr). All three also have broad H$\beta$ emission with $\sigma > 1000 \,\rm km s^{-1}$, a common feature of such `little red dots'. From S\'ersic profile fits to the NIRCam images in F200W we find that the stellar light of galaxies is extremely compact: the galaxies have half-light radii of $r_{\rm e}\sim$ 100 pc, in the regime of ultra compact dwarfs in the nearby Universe. Their masses are uncertain, as they depend on the contribution of possible light from an AGN to the flux at $\lambda_{\rm rest}>5000$ $\r{A}$. If the AGN contribution is low beyond the Balmer break region, the masses are $M_* \sim 10^{10}-10^{11}\mathrm{M}_{\odot}$, and the central densities are higher than those of any other known galaxy population by an order of magnitude. Interestingly, the implied velocity dispersions of $\sim$1500 kms$^{-1}$ are in very good agreement with the measured H$\beta$ line widths. We suggest that some of the broad lines in `little red dots' are not due to AGNs but simply reflect the kinematics of the galaxies, and speculate that the galaxies are observed in a short-lived phase where the central densities are much higher than at later times. We stress, however, that the canonical interpretation of AGNs causing the broad H$\beta$ lines also remains viable.

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 spatially resolved, simple stellar population equivalent ages, stellar metallicities, and abundance ratios for 456 massive ($10.3\lesssim\log(\mathrm{M}_*/\mathrm{M}_\odot)\lesssim11.8$) quiescent galaxies at $0.6\lesssim z\lesssim1.0$ from the Large Early Galaxy Astrophysics Census, derived using full-spectrum models. Typically, we find flat age and [Mg/Fe] gradients, and negative [Fe/H] gradients, implying iron-rich cores. We also estimate intrinsic [Fe/H] gradients via forward modelling. We examine the observed gradients in three age bins. Younger quiescent galaxies typically have negative [Fe/H] gradients and positive age gradients, possibly indicating a recent central starburst. Additionally, this finding suggests that photometrically measured flat colour gradients in young quiescent galaxies are the result of the positive age and negative metallicity gradients cancelling each other. For older quiescent galaxies, the age gradients become flat and [Fe/H] gradients weaken, though remain negative. Thus, negative colour gradients at older ages are likely driven by metallicity gradients. The diminishing age gradient may result from the starburst fading. Furthermore, the persistence of the [Fe/H] gradients may suggest that the outskirts are simultaneously built up by mergers with lower metallicity satellites. On the other hand, the gradients could be inherited from the star-forming phase, in which case mergers may not be needed to explain our findings. This work illustrates the need for resolved spectroscopy, instead of just photometry, to measure stellar population gradients. Extending these measurements to higher redshift is imperative for understanding how stellar populations in quiescent galaxies are assembled over cosmic time., Comment: Accepted for publication in MNRAS; minor typesetting corrections after copyediting

The stellar initial mass function (IMF) is thought to be bottom-heavy in the cores of the most massive galaxies, with an excess of low mass stars compared to the Milky Way. However, studies of the kinematics of quiescent galaxies at 27, to reduce tensions with galaxy formation models. Here we explore 'ski slope' IMFs that are simultaneously bottom-heavy, with a steep slope at low stellar masses, and top-heavy, with a shallow slope at high masses. We derive a form of the IMF for massive galaxies that is consistent with measurements in the local universe and yet produces relatively low M/L ratios at high redshift. This concordance IMF has slopes $\gamma_1=2.40\pm0.09$, $\gamma_2=2.00\pm0.14$, and $\gamma_3=1.85\pm0.11$ in the regimes 0.08-0.5 Msun, 0.5-1 Msun, and >1 Msun respectively. The IMF parameter $\alpha$, the mass excess compared to a Milky Way IMF, ranges from $\log(\alpha)\approx+0.3$ for present-day galaxies to $\log(\alpha)\approx-0.1$ for their star forming progenitors. The concordance IMF applies only to the central regions of the most massive galaxies, with velocity dispersions ~300 km/s, and their progenitors. However, it can be generalized using a previously-measured relation between $\alpha$ and $\sigma$. We arrive at the following modification to the Kroupa (2001) IMF for galaxies with $\sigma\gtrsim 160$ km/s: $\gamma_1\approx1.3+4.3\log\sigma_{160}$; $\gamma_2\approx2.3-1.2\log\sigma_{160}$; and $\gamma_3\approx2.3-1.7\log\sigma_{160}$, with $\sigma_{160}=\sigma/160$ km/s. If galaxies grow primarily inside-out, so that velocity dispersions are relatively stable, these relations should also hold at high redshift., Comment: Accepted for publication in ApJ Letters. Figure 3 shows the proposed IMF for massive galaxies. For those not liking 3D figures: the information in Fig. 2 is mostly illustrative

The Dragonfly Telephoto Array employs a unique design to detect very large and diffuse galaxies, which might be missed with conventional telescopes. The Dragonfly Ultrawide Survey (DFUWS) is a new wide-field survey which will cover 10,000 deg$^2$ of the northern sky, and it provides an ideal dataset to find these large diffuse galaxies. From 3100 deg$^2$ of DFUWS data, we identified eleven large, low surface brightness galaxies as a pilot sample for spectroscopic follow-up. These are the largest galaxies in the examined area that appear smooth and isolated, with effective radii of 12"-27". Eight are below 24 $\mathrm{mag\,arcsec^{-2}}$ in central $g$-band surface brightness. Keck Cosmic Web Imager (KCWI) spectra of the diffuse light show that all eleven galaxies in this sample are quiescent, and seven qualify as ultra-diffuse galaxies (UDGs). Eight galaxies have distances between 15 and 30 Mpc, while the other three are in the Pegasus cluster at 50 Mpc. Their spectra show evidence of a $\sim 1$Gyr old stellar population in addition to an even older stellar population. The intermediate-age component is present in group and satellite galaxies but not in the Pegasus cluster UDGs. All galaxies in this sample are detected in both Dragonfly and Legacy imaging, and the sample partially overlaps with existing UDG catalogs. This pilot sample provides an excellent training set for our analysis of the upcoming full 10,000 deg$^2$ DFUWS data, from which we may expect to discover even larger, previously-unknown galaxies., Comment: Submitted to ApJ

With upcoming wide field surveys from the ground and space the number of known dwarf galaxies at $\lesssim 25$ Mpc is expected to dramatically increase. Insight into their nature and analyses of these systems' intrinsic properties will rely on reliable distance estimates. Traditional techniques such as tip of the red giant branch (TRGB) or surface brightness fluctuations (SBF) are limited in their widespread applicability, especially in the semi-resolved regime. In this work we turn to the rapidly growing field of simulation based inference to infer distances, and other physical properties, of dwarf galaxies directly from multi-band images. We introduce Silkscreen: a code leveraging neural posterior estimation to infer the posterior distribution of parameters while simultaneously training a convolutional neural network to extract summary statistics from the images. Utilizing this combination of machine learning and Bayesian inference, we demonstrate the method's ability to recover accurate distances from ground-based survey images for a set of nearby galaxies ($2 < D ({\rm Mpc)} < 12$) with measured SBF or TRGB distances. We discuss caveats of the current implementation along with future prospects, focusing on the goal of amortized inference. For surveys like LSST, Silkscreen, once amortized, will facilitate inference for new dwarf galaxies in a matter of seconds using only broadband cutouts. While we focus here on dwarf galaxies, we note that this method can be generalized to more luminous systems as well., Comment: Submitted to AAS Journals, Comments Welcome! Code available here: https://github.com/tbmiller-astro/silkscreen

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.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 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 that form at higher redshift have higher [Mg/Fe] and lower [Fe/H] and [Mg/H], irrespective of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily explained by lower redshift samples naturally including galaxies with longer star-formation timescales. On the other hand, the lower [Mg/H] can be explained by galaxies forming at earlier epochs expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is also lower at z=1-3 compared to the lower redshift relation, though the slopes are similar. 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: Submitted to ApJ; 18 pages, 6 figures, 1 table

The Dragonfly Spectral Line Mapper is an innovative all-refracting telescope designed to carry out ultra-low surface brightness wide-field mapping of visible wavelength line emission. Equipped with ultranarrowband (0.8 nm bandwidth) filters mounted in Dragonfly Filter-Tilter instrumentation, the Dragonfly Spectral Line Mapper maps H$\alpha$, [NII]$\lambda$6583, and [OIII]$\lambda$5007 line emission produced by structures with sizes ranging from $\sim$1 to 1000 kpc in the local Universe. These spatial scales encompass that of the exceedingly diffuse and faintly radiating circumgalactic medium, which is singularly difficult to detect with conventional mirror-based telescope instrumentation. Extremely careful control of systematics is required to directly image these large scale structures, necessitating high fidelity sky background subtraction, wavelength calibration, and specialized flat-fielding methods. In this paper, we discuss the on-sky performance of the Dragonfly Spectral Line Mapper with these methods in place., Comment: 10 pages, 5 figures, SPIE Astronomical Telescopes and Instrumentation 2024 Proceedings

The Dragonfly Spectral Line Mapper (DSLM) is a semi-autonomous, distributed-aperture based telescope design, featuring a modular setup of 120 Canon telephoto lenses, and equal numbers of ultra-narrowband filters, detectors, and other peripherals. Here we introduce the observatory software stack for this highly-distributed system. Its core is the Dragonfly Communication Protocol (DCP), a pure-Python hardware communication framework for standardized hardware interaction. On top of this are 120 REST-ful FastAPI web servers, hosted on Raspberry Pis attached to each unit, orchestrating command translation to the hardware and providing diagnostic feedback to a central control system running the global instrument control software. We discuss key features of this software suite, including docker containerization for environment management, class composition as a flexible framework for array commands, and a state machine algorithm which controls the telescope during autonomous observations., Comment: 18 pages, presented at the SPIE Astronomical Telescopes and Instrumentation conference in Yokohama, Japan

The Dragonfly Spectral Line Mapper is a mosaic telescope comprising 120 Canon telephoto lenses, based on the design of the Dragonfly Telephoto Array. With a wide field of view, and the addition of the "Dragonfly Filter-Tilter" instrumentation holding ultra narrow bandpass filters in front of each lens, the Dragonfly Spectral Line mapper is optimized for ultra low surface brightness imaging of visible wavelength line emission. The Dragonfly Spectral Line Mapper was constructed and commissioned in four phases from March 2022 to November 2023. During this time, four individual mounts of 30 lenses each were constructed and commissioned. The commissioning of the telescope included the deployment of the "Dragonfly StarChaser" which carries out image stabilization corrections in the telephoto lens, to enable hour-long exposures to be taken. In addition, we introduced new instrumentation such as a film to cover the optics to keep the filters clean. Here we describe the updated design of the complete 120-lens array, and the implementation of the instrumentation described above. Additionally, we present updated characterization of the cameras and filter transmission for the full array. Finally, we reflect on the construction and commissioning process of the complete 120-lens array Dragonfly Spectral Line Mapper, and remark on the feasibility of a larger 1000-lens array., Comment: 12 pages, 9 figures, SPIE conference proceedings

CDG-1 is a tight grouping of four likely globular clusters in the Perseus cluster, and a candidate dark galaxy with little or no diffuse light. Here we provide new constraints on the luminosity of any underlying stellar emission, using HST/UVIS F200LP imaging. No diffuse emission is detected, with a 2$\sigma$ upper limit of F200LP>28.1 mag/arcsec$^2$ on the 5'' scale of CDG-1. This surface brightness limit corresponds to a 2$\sigma$ lower limit of >0.5 for the fraction of the total luminosity that is in the form of globular clusters. The most likely alternative, although improbable, is that CDG-1 is a chance grouping of four globular clusters in the halo of the Perseus galaxy IC312., Comment: To appear in Research Notes of the American Astronomical Society

We present spatially resolved Keck Cosmic Web Imager stellar spectroscopy of the Virgo cluster dwarf galaxies VCC 9 and VCC 1448. These galaxies have similar stellar masses and large half-light radii but very different globular cluster (GC) system richness ($\sim$25 vs. $\sim$99 GCs). Using the KCWI data, we spectroscopically confirm 10 GCs associated with VCC 1448 and one GC associated with VCC 9. We make two measurements of dynamical mass for VCC 1448 based on the stellar and GC velocities respectively. VCC 1448's mass measurements suggest that it resides in a halo in better agreement with the expectation of the stellar mass -- halo mass relationship than the expectation from its large GC counts. For VCC 9, the dynamical mass we measure agrees with the expected halo mass from both relationships. We compare VCC 1448 and VCC 9 to the GC-rich galaxy Dragonfly 44 ($\sim74$ GCs), which is similar in size but has $\sim 1$ dex less stellar mass than either Virgo galaxy. In dynamical mass -- GC number space, Dragonfly 44 and VCC 1448 exhibit richer GC systems given their dynamical mass than that of VCC 9 and other `normal' galaxies. We also place the galaxies in kinematics -- ellipticity space finding evidence of an anticorrelation between rotational support and the fraction of a galaxy's stellar mass in its GC system. i.e., VCC 9 is more rotationally supported than VCC 1448, which is more rotationally supported than Dragonfly 44. This trend may be expected if a galaxy's GC content depends on its natal gas properties at formation., Comment: 17 pages, 12 figures + 1 in the Appendix, Accepted for publication in MNRAS

We present a study of the growth of the quiescent galaxy population between 0.5 < z < 3 by tracing the number density and structural evolution of a sample of 4518 old and 583 young quiescent galaxies with log($M_*$/$M_{\odot}$)>10.4, selected from the COSMOS2020 catalog with complementary HST/F160W imaging from the 3D-DASH survey. Among the quiescent population at z$\sim$2, roughly 50% are recently quenched galaxies; these young quiescent galaxies become increasingly rare towards lower redshift, supporting the idea that the peak epoch of massive galaxy quenching occurred at z>2. Our data show that while the effective half-light radii of quiescent galaxies generally increases with time, young quiescent galaxies are significantly smaller than their older counterparts at the same redshift. In this work we investigate the connection between this size difference and other structural properties, including axis ratios, color gradients, stellar mass, and the intrinsic scatter in effective radii. We demonstrate that the size difference is driven by the most massive sub-population (log($M_*$/$M_{\odot}$)>11) and does not persist when restricting the sample to intermediate mass galaxies (10.4

The identification of red, apparently massive galaxies at $z>7$ in early JWST photometry suggests a strongly accelerated timeline compared to standard models of galaxy growth. A major uncertainty in the interpretation is whether the red colors are caused by evolved stellar populations, dust, or other effects such as emission lines or AGN. Here we show that three of the massive galaxy candidates at $z=6.7-8.4$ have prominent Balmer breaks in JWST/NIRSpec spectroscopy from the RUBIES program. The Balmer breaks demonstrate unambiguously that stellar emission dominates at $\lambda_{\rm rest} = 0.4\,\mu$m, and require formation histories extending hundreds of Myr into the past in galaxies only 600--800 Myr after the Big Bang. Two of the three galaxies also show broad Balmer lines, with H$\beta$ FWHM $>2500~{\rm km\,s^{-1}}$, suggesting that dust-reddened AGN contribute to, or even dominate, the SEDs of these galaxies at $\lambda_{\rm rest}\gtrsim 0.6\,\mu$m. All three galaxies have relatively narrow [O III] lines, seemingly ruling out a high-mass interpretation if the lines arise in dynamically-relaxed, inclined disks. Yet, the inferred masses also remain highly uncertain. We model the high-quality spectra using Prospector to decompose the continuum into stellar and AGN components, and explore limiting cases in stellar/AGN contribution. This produces a wide range of possible stellar masses, spanning $M_\star \sim 10^9 - 10^{11}\,{\rm M_{\odot}}$. Nevertheless, all fits suggest a very early and rapid formation, most of which follow with a truncation in star formation. Potential origins and evolutionary tracks for these objects are discussed, from the cores of massive galaxies to low-mass galaxies with over-massive black holes. Intriguingly, we find all of these explanations to be incomplete; deeper and redder data are needed to understand the physics of these systems., Comment: Accepted for publication in ApJL; 20 pages, 10 figures, 1 table

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 present results from the JWST First Reionization Epoch Spectroscopically Complete Observations survey (FRESCO) on the star forming sequence of galaxies at $1.09.5 M_\odot$ that are lower than found in many earlier studies by up to 0.6 dex, but in good agreement with recent results obtained with the Prospector fitting framework. The difference log(SFR(Pa$\alpha$)-SFR(Prospector)) is -0.09 $\pm$ 0.04 dex at $10^{10-11} M_\odot$. We also measure the empirical relation between Paschen-$\alpha$ luminosity and rest-frame H band magnitude and find that the scatter is only 0.04 dex lower than that of the SFR-M* relation and is much lower than the systematic differences among relations in the literature due to various methods of converting observed measurements to physical properties. We additionally identify examples of sources -- that, with standard cutoffs via the UVJ diagram, would be deemed quiescent -- with significant, typically extended, Paschen-$\alpha$ emission. Our results may be indicative of the potential unification of methods used to derive the star forming sequence with careful selection of star forming galaxies and independent star formation rate and stellar mass indicators., Comment: 17 pages, 9 figures. Accepted for publication in ApJ

We report the spectroscopic confirmation of a massive ($\log(M_\star/M_\odot)=10.34 \pm_{0.07}^{0.06}$), HST-dark ($m_\mathrm{F150W} - m_\mathrm{F444W} = 3.6$) quiescent galaxy at $z_{spec}=3.97$ in the UNCOVER survey. NIRSpec/PRISM spectroscopy and a non-detection in deep ALMA imaging surprisingly reveals that the galaxy is consistent with a low ($<$10 $M_\odot \ \mathrm{yr^{-1}}$) star formation rate despite evidence for moderate dust attenuation. The F444W image is well modeled with a two component \sersic fit that favors a compact, $r_e\sim200$ pc, $n\sim2.9$ component and a more extended, $r_e\sim1.6$ kpc, $n\sim1.7$ component. The galaxy exhibits strong color gradients: the inner regions are significantly redder than the outskirts. Spectral energy distribution models that reproduce both the red colors and low star formation rate in the center of UNCOVER 18407 require both significant ($A_v\sim1.4$ mag) dust attenuation and a stellar mass-weighted age of 900 Myr, implying 50\% of the stars in the core already formed by $z=7.5$. Using spatially resolved annular mass-to-light measurements enabled by the galaxy's moderate magnification ($\mu=2.12\pm_{0.01}^{0.05}$) to reconstruct a radial mass profile from the best-fitting two-component \sersic model, we infer a total mass-weighted $r_\mathrm{eff} = 0.72 \pm_{0.11}^{0.15}$ kpc and log$(\Sigma_\mathrm{1 kpc} \ [\mathrm{M_\odot/kpc^2}]) = 9.61 \pm_{0.10}^{0.08}$. The early formation of a dense, low star formation rate, and dusty core embedded in a less attenuated stellar envelope suggests an evolutionary link between the earliest-forming massive galaxies and their elliptical descendants. Furthermore, the disparity between the global, integrated dust properties and the spatially resolved gradients highlights the importance of accounting for radially varying stellar populations when characterizing the early growth of galaxy structure., Comment: 17 pages, 9 figures, 2 tables. Resubmitted to ApJ after response to referee and update to include new medium band imaging from the JWST MEGASCIENCE program. Comments welcome!

We present the elemental abundances and ages of 19 massive quiescent galaxies at $z\sim1.4$ and $z\sim2.1$ from the Keck Heavy Metal Survey. The ultra-deep LRIS and MOSFIRE spectra were modeled using a full-spectrum stellar population fitting code with variable abundance patterns. The galaxies have iron abundances between [Fe/H] = -0.5 and -0.1 dex, with typical values of $-0.2$ [$-0.3$] at $z\sim1.4$ [$z\sim2.1$]. We also find a tentative $\log\sigma_v$-[Fe/H] relation at $z\sim1.4$. The magnesium-to-iron ratios span [Mg/Fe] = 0.1--0.6 dex, with typical values of $0.3$ [$0.5$] dex at $z\sim1.4$ [$z\sim2.1$]. The ages imply formation redshifts of $z_{\rm form}=2-8$. Compared to quiescent galaxies at lower redshifts, we find [Fe/H] was $\sim0.2$ dex lower at $z=1.4-2.1$. We find no evolution in [Mg/Fe] out to $z\sim1.4$, though the $z\sim2.1$ galaxies are $0.2$ dex enhanced compared to $z=0-0.7$. A comparison of these results to a chemical evolution model indicates that galaxies at higher redshift form at progressively earlier epochs and over shorter star-formation timescales, with the $z\sim2.1$ galaxies forming the bulk of their stars over 150 Myr at $z_{\rm form}\sim4$. This evolution cannot be solely attributed to an increased number of quiescent galaxies at later times; several Heavy Metal galaxies have extreme chemical properties not found in massive galaxies at $z\sim0.0-0.7$. Thus, the chemical properties of individual galaxies must evolve over time. Minor mergers also cannot fully account for this evolution as they cannot increase [Fe/H], particularly in galaxy centers. Consequently, the build-up of massive quiescent galaxies since $z\sim2.1$ may require further mechanisms such as major mergers and/or central star formation., Comment: Submitted to ApJ (21 pages, 11 figures)

Over the past year, JWST has uncovered galaxies at record-breaking distances up to $z \sim 13$. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster Abell 2744 over $\sim 45$ arcmin$^{2}$ down to $\sim 29.5$ AB mag. Here, we present a robust ultraviolet (UV) luminosity function derived through lensing clusters at $99$ in recent JWST studies. In particular, a double power-law (DPL) describes better the bright-end of the luminosity function compared to the classical Schechter form. The number density of these bright galaxies is 10-100 times larger than theoretical predictions and previous findings based on Hubble Space Telescope (HST) observations. Additionally, we measure a star formation rate density of $\rho_{\rm SFR} = 10^{-2.64}$ M$_{\odot}$~yr$^{-1}$~Mpc$^{-3}$ at these redshifts, which is 4 to 10 times higher than galaxy formation models that assume a constant star formation efficiency. Future wide-area surveys and accurate modeling of lensing-assisted observations will reliably constrain both the bright and the dim end of the UV luminosity function at $z>9$, which will provide key benchmarks for galaxy formation models., Comment: Published in MNRAS

Astronomical images often have regions with missing or unwanted information, such as bad pixels, bad columns, cosmic rays, masked objects, or residuals from imperfect model subtractions. In certain situations it can be essential, or preferable, to fill in these regions. Most existing methods use low order interpolations for this task. In this paper a method is described that uses the full information that is contained in the pixels just outside masked regions. These edge pixels are extrapolated inwards, using iterative median filtering. This leads to a smoothly varying spatial resolution within the filled-in regions, and ensures seamless transitions between masked pixels and good pixels. Gaps in continuous, narrow features can be reconstructed with high fidelity, even if they are large. The method is implemented in maskfill, an open-source MIT licensed Python package. Its performance is illustrated with several examples, and compared to several alternative interpolation schemes., Comment: Accepted for publication in PASP. The code is available at https://github.com/dokkum/maskfill . More information, including example runs, at https://maskfill.readthedocs.io

In this paper, we present the Heavy Metal Survey, which obtained ultradeep medium-resolution spectra of 21 massive quiescent galaxies at $1.3

We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emission from NIRCam and Dragonfly imaging at smaller distances, we trace the surface brightness and color profiles of this galaxy over the entire extent of its predicted dark matter halo, from 0.1-100 kpc. Controlled N-body simulations have predicted that minor mergers create "accretion shelves" in the surface brightness profile at large radius. We observe such a feature in Ark 227 at 10-20 kpc, which, according to models, could be caused by a merger with total mass ratio 1:10. The metallicity declines over this radial range, further supporting the minor merger scenario. There is tentative evidence of a second shelf at $\mu_V\approx 35$ mag arcsec$^{-2}$ extending from 50-100 kpc, along with a corresponding drop in metallicity. The stellar mass in this outermost envelope is $\approx10^7M_\odot$. These results suggest that Ark 227 experienced multiple mergers with a spectrum of lower-mass galaxies -- a scenario that is broadly consistent with the hierarchical growth of structure in a cold dark matter-dominated universe. Finally, we identify an ultra-faint dwarf associated with Ark 227 with $M_\ast\approx10^5 M_\odot$ and $\mu_{V,e}=28.1$ mag arcsec$^{-2}$, demonstrating that JWST is capable of detecting very low-mass dwarfs to distances of at least ~30 Mpc., Comment: 10 pages, 9 figures, submitted to ApJ

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

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 a uniform analysis of the stellar initial mass function (IMF) from integrated light spectroscopy of 15 compact stellar systems (11 globular clusters in M31 and 4 ultra compact dwarfs in the Virgo cluster, UCDs) and two brightest Coma cluster galaxies (BCGs), covering a wide range of metallicities ($-$1.7 $<$ [Fe/H] $<$ 0.01) and velocity dispersions (7.4 km~s$^{-1}$ $< \sigma <$ 275 km~s$^{-1}$). The S/N $\sim 100$ \AA$^{-1}$ Keck LRIS spectra are fitted over the range $4000<\lambda/\mbox{\AA}<10,000$ with flexible, full-spectrum stellar population synthesis models. We use the models to fit simultaneously for ages, metallicities, and individual elemental abundances of the population, allowing us to decouple abundance variations from variations in IMF slope. We show that compact stellar systems do not follow the same trends with physical parameters that have been found for early-type galaxies. Most globular clusters in our sample have an IMF consistent with that of the Milky Way, over a wide range of [Fe/H] and [Mg/Fe]. There is more diversity among the UCDs, with some showing evidence for a bottom-heavy IMF, but with no clear correlation with metallicity, abundance, or velocity dispersion. The two Coma BCGs have similar velocity dispersion and metallicity, but we find the IMF of NGC~4874 is consistent with that of the Milky Way while NGC~4889 presents evidence for a significantly bottom-heavy IMF. For this sample, the IMF appears to vary between objects in a way that is not explained by a single metallicity-dependent prescription., Comment: Accepted for publication in MNRAS

NGC1052-DF4 was found to be the second "galaxy lacking dark matter" in the NGC1052 group, based on its velocity dispersion of $\sigma_{\rm gc}=4.2^{+4.4}_{-2.2}$ km/s as measured from the radial velocities of seven of its globular clusters. Here we verify this result by measuring the stellar velocity dispersion of the galaxy. We observed the diffuse stellar light in NGC1052-DF4 with the Keck Cosmic Web Imager (KCWI) in its highest resolution mode, with $\sigma_{\mathrm{instr}}\approx 7$ km/s. With a total science + sky exposure time of 34hrs, the resulting spectrum is exceptional both in its spectral resolution and its S/N ratio of 23\r{A}$^{-1}$. We find a stellar velocity dispersion of $\sigma_{\rm stars} = 8.0^{+2.3}_{-1.9}$ km/s, consistent with the previous measurement from the globular clusters. Combining both measurements gives a fiducial dispersion of $\sigma_{\rm f} = 6.3_{-1.6}^{+2.5}$ km/s. The implied dynamical mass within the half-light radius is $8_{-4}^{+6} \times 10^7 M_{\odot}$. The expected velocity dispersion of NGC1052-DF4 from the stellar mass alone is $7 \pm 1$ km/s, and for an NFW halo that follows the stellar mass -- halo mass relation and the halo mass -- concentration relation, the expectation is $\sim 30$ km/s. The low velocity dispersion rules out a normal NFW dark matter halo, and we confirm that NGC1052-DF4 is one of at least two galaxies in the NGC1052 group that have an anomalously low dark matter content. While any viable model for their formation should explain the properties of both galaxies, we note that NGC1052-DF4 now poses the largest challenge as it has the most stringent constraints on its dynamical mass., Comment: Accepted to ApJ, refereed version

One of the surprising results from HST was the discovery that many of the most massive galaxies at z~2 are very compact, having half-light radii of only 1-2 kpc. The interpretation is that massive galaxies formed inside-out, with their cores largely in place by z~2 and approximately half of their present-day mass added later through minor mergers. Here we present a compact, massive, quiescent galaxy at $z_{\rm phot}=1.94^{+0.13}_{-0.17}$ with a complete Einstein ring. The ring was found in the JWST COSMOS-Web survey and is produced by a background galaxy at $z_{\rm phot}=2.98^{+0.42}_{-0.47}$. Its 1.54" diameter provides a direct measurement of the mass of the "pristine" core of a massive galaxy, observed before mixing and dilution of its stellar population during the 10 Gyr of galaxy evolution between z=2 and z=0. We find a mass of $M_{\rm lens}=6.5^{+3.7}_{-1.5} \times 10^{11}$ Msun within a radius of 6.6 kpc. The stellar mass within the same radius is $M_{\rm stars}= 1.1^{+0.2}_{-0.3} \times 10^{11}$ Msun for a Chabrier initial mass function (IMF), and the fiducial dark matter mass is $M_{\rm dm} = 2.6^{+1.6}_{-0.7} \times 10^{11}$ Msun. Additional mass is needed to explain the lensing results, either in the form of a higher-than-expected dark matter density or a bottom-heavy IMF., Comment: Nature Astronomy, in press. As it happens this is the highest redshift lens currently known. Found by eyeballing the wonderful COSMOS-Web data

We present the survey design and initial results of the ALMA Cycle 9 program of DUALZ, which aims to establish a joint ALMA and JWST public legacy field targeting the massive galaxy cluster Abell 2744. DUALZ features a contiguous $4'\times6'$ ALMA 30-GHz-wide mosaic in Band 6, covering areas of $\mu>2$ down to a sensitivity of $\sigma=32.7~\mu$Jy. Through a blind search, we identified 69 dust continuum sources at S/N $\gtrsim5.0$ with median redshift and intrinsic 1.2-mm flux of $z=2.30$ and $S_{\rm 1.2mm}^{\rm int}=0.24$~mJy. Of these, 27 have been spectroscopically confirmed, leveraged by the latest NIRSpec observations, while photometric redshift estimates are constrained by the comprehensive HST, NIRCam, and ALMA data for the remaining sources. With priors, we further identify a [CII]158 $\mu$m line emitter at $z=6.3254\pm0.0004$, confirmed by the latest NIRSpec spectroscopy. The NIRCam counterparts of the 1.2-mm continuum exhibit undisturbed morphologies, denoted either by disk or spheroid, implying the triggers for the faint mm emission are less catastrophic than mergers. We have identified 8 HST-dark galaxies (F150W$>$27mag, F150W$-$F444W$>$2.3) and 2 JWST-dark (F444W$>$30mag) galaxy candidates among the ALMA continuum sources. The former includes face-on disk galaxies, hinting that substantial dust obscuration does not always result from inclination. We also detect a marginal dust emission from an X-ray-detected galaxy at $z_{\rm spec}=10.07$, suggesting an active co-evolution of the central black hole and its host. We assess the infrared luminosity function up to $z\sim10$ and find it consistent with predictions from galaxy formation models. To foster diverse scientific outcomes from the community, we publicly release reduced ALMA mosaic maps, cubes, and the source catalog., Comment: 33 pages, 16 figures, and 5 tables. Submitted to ApJS. The ALMA products are fully available from here: https://jwst-uncover.github.io/DR2.html#DUALZ

JWST is revealing a new population of dust-reddened broad-line active galactic nuclei (AGN) at redshifts $z\gtrsim5$. Here we present deep NIRSpec/Prism spectroscopy from the Cycle 1 Treasury program UNCOVER of 15 AGN candidates selected to be compact, with red continua in the rest-frame optical but with blue slopes in the UV. From NIRCam photometry alone, they could have been dominated by dusty star formation or AGN. Here we show that the majority of the compact red sources in UNCOVER are dust-reddened AGN: $60\%$ show definitive evidence for broad-line H$\alpha$ with FWHM$\, >2000$ km/s, for $20\%$ current data are inconclusive, and $20\%$ are brown dwarf stars. We propose an updated photometric criterion to select red $z>5$ AGN that excludes brown dwarfs and is expected to yield $>80\%$ AGN. Remarkably, among all $z_{\rm phot}>5$ galaxies with F277W$-$F444W$>1$ in UNCOVER at least $33\%$ are AGN regardless of compactness, climbing to at least $80\%$ AGN for sources with F277W$-$F444W$>1.6$. The confirmed AGN have black hole masses of $10^7-10^9$ M$_{\odot}$. While their UV-luminosities ($-16>M_{\rm UV}>-20$ AB mag) are low compared to UV-selected AGN at these epochs, consistent with percent-level scattered AGN light or low levels of unobscured star formation, the inferred bolometric luminosities are typical of $10^7-10^9$ M$_{\odot}$ black holes radiating at $\sim 10-40\%$ of Eddington. The number densities are surprisingly high at $\sim10^{-5}$ Mpc$^{-3}$ mag$^{-1}$, 100 times more common than the faintest UV-selected quasars, while accounting for $\sim1\%$ of the UV-selected galaxies. While their UV-faintness suggest they may not contribute strongly to reionization, their ubiquity poses challenges to models of black hole growth., Comment: 23 pages, 9 figures, 5 tables, submitted to ApJ

We report the detection of a 13$\sigma$ H$\alpha$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$\alpha$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$\alpha$ emission appears to arise entirely from the less-obscured southern component and shows a similar $\Delta$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$\alpha$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$\alpha$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures., Comment: Submitted to MNRAS. 20 pages, 10 figures and 8 tables (including appendices)

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

We present JWST NIRSpec prism spectroscopy of gravitationally lensed galaxies at $z\gtrsim9$ found behind the massive galaxy cluster Abell 2744 in the UNCOVER Cycle 1 Treasury Program. We confirm the source redshift via emission lines and/or the Ly$\alpha$ break feature for ten galaxies at z=8.50-13.08 down to $M_{\rm UV}=-17.3$. We achieve a high confirmation rate of 100\% for $z>9$ candidates reported in Atek et al. (2023). Using six sources with multiple emission line detections, we find that the offset of the redshift estimates between the lines and the Ly$\alpha$ break alone with prism can be as large as $\pm0.2$, raising caution in designing future follow-up spectroscopy for the break-only sources. With spec-$z$ confirmed sources in UNCOVER and the literature, we derive lower limits on the rest-frame ultraviolet (UV) luminosity function (LF) at $z\simeq9$-12 and find these lower limits to be consistent with recent photometric measurements. We identify at least two unambiguous and several possible active galactic nucleus (AGN) systems based on X-ray emission, broad line (BL) H$\beta$, high ionization line (e.g., NIV]1487, CIV1549) detections, and excess in UVLF. This requires the AGN LFs at $z\simeq$ 9-10 to be comparable or even higher than the X-ray AGN LF estimated at $z\sim6$ and indicates a plausible cause of the high abundance of $z>9$ galaxies claimed in recent photometric studies may be AGNs. One UV-luminous source is confirmed at the same redshift as a dusty BL AGN at $z=8.50$ with a physical separation of 380 kpc in the source plane. These two sources show blueward Ly$\alpha$ line or continuum emission, suggesting that they reside in the same ionized bubble with a radius of $7.69\pm0.18$ pMpc. Our results imply that AGNs have a non-negligible contribution to cosmic reionization., Comment: 27 pages, 11 figures, 4 tables, submitted to ApJ (See also arXiv:2308.11610)

We report JWST/NIRSpec spectra of three distant T-type brown dwarfs identified in the Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER) survey of the Abell 2744 lensing field. One source was previously reported as a candidate T dwarf on the basis of NIRCam photometry, while two sources were initially identified as candidate active galactic nuclei. Low-resolution 1--5 $\mu$m spectra confirm the presence of molecular features consistent with T dwarf atmospheres, and comparison to spectral standards infers classifications of sdT1, T6, and T8--T9. The warmest source, UNCOVER-BD-1, shows evidence of subsolar metallicity, and atmosphere model fits indicates T$_{eff}$ = 1300 K and [M/H] $\sim$ $-$1.0, making this one of the few spectroscopically-confirmed T subdwarfs known. The coldest source, UNCOVER-BD-3, is near the T/Y dwarf boundary with T$_{eff}$ = 550 K, and our analysis indicates the presence of PH$_3$ in the 3--5~$\mu$m region, favored over CO$_2$ and a possible indicator of subsolar metallicity. We estimate distances of 0.9--4.5 kpc from the Galactic midplane, making these the most distant brown dwarfs with spectroscopic confirmation. Population simulations indicate high probabilities of membership in the Galactic thick disk for two of these brown dwarfs, and potential halo membership for UNCOVER-BD-1. Our simulations indicate that there are approximately 5 T dwarfs and 1--2 L dwarfs in the Abell 2744 field down to F444W = 30 AB mag, roughly one-third of which are thick disk members. These results highlight the utility of deep JWST/NIRSpec spectroscopy for identifying and characterizing the oldest metal-poor brown dwarfs in the Milky Way., Comment: revised, accepted by ApJ 22 Nov 2023

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

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

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies., Comment: ApJL in press; 16 pages, 6 figures, 2 tables

The James Webb Space Telescope is now detecting early black holes (BHs) as they transition from "seeds" to supermassive BHs. Recently Bogdan et al. (2023) reported the detection of an X-ray luminous supermassive BH, UHZ-1, with a photometric redshift at $z > 10$. Such an extreme source at this very high redshift provides new insights on seeding and growth models for BHs given the short time available for formation and growth. Harnessing the exquisite sensitivity of JWST/NIRSpec, here we report the spectroscopic confirmation of UHZ-1 at $z = 10.073 \pm 0.002$. We find that the NIRSpec/Prism spectrum is typical of recently discovered z~10 galaxies, characterized primarily by star-formation features. We see no clear evidence of the powerful X-ray source in the rest-frame UV/optical spectrum, which may suggest heavy obscuration of the central BH, in line with the Compton-thick column density measured in the X-rays. We perform a stellar population fit simultaneously to the new NIRSpec spectroscopy and previously available photometry. The fit yields a stellar mass estimate for the host galaxy that is significantly better constrained than prior photometric estimates ($M_*\sim 1.4^{+0.3}_{-0.4} \times 10^8 M_\odot$). Given the predicted BH mass ($M_{\rm BH}\sim10^7-10^8 M_\odot$), the resulting ratio of $M_{\rm BH}/M_*$ remains two to three orders of magnitude higher than local values, thus lending support to the heavy seeding channel for the formation of supermassive BHs within the first billion years of cosmic evolution., Comment: 10 pages, 4 figures, ApJL in-press. Updated to accepted version incl. X-ray image

We analyse the joint distribution of dust attenuation and projected axis ratios, together with galaxy size and surface brightness profile information, to infer lessons on the dust content and star/dust geometry within star-forming galaxies at 0 < z <2.5. To do so, we make use of large observational datasets from KiDS+VIKING+HSC-SSP and extend the analysis out to redshift z = 2.5 using the HST surveys CANDELS and 3D-DASH. We construct suites of SKIRT radiative transfer models for idealized galaxies observed under random viewing angles with the aim of reproducing the aforementioned distributions, including the level and inclination dependence of dust attenuation. We find that attenuation-based dust mass estimates are at odds with constraints from far-infrared observations, especially at higher redshifts, when assuming smooth star and dust geometries of equal extent. We demonstrate that UV-to-near-IR and far-infrared constraints can be reconciled by invoking clumpier dust geometries for galaxies at higher redshifts and/or very compact dust cores. We discuss implications for the significant wavelength- and redshift-dependent differences between half-light and half-mass radii that result from spatially varying dust columns within -- especially massive -- star-forming galaxies., Comment: Accepted for publication in MNRAS

Unbiased sky background modeling is crucial for the analysis of deep wide-field images, but it remains a major challenge in low surface brightness astronomy. Traditional image processing algorithms are often designed to produce artificially flat backgrounds, erasing astrophysically meaningful structures. In this paper, we present three ideas that can be combined to produce wide-field astronomical data that preserve accurate representations of the background sky: (1) Use of all-sky infrared/sub-mm data to remove the large-scale time-varying components while leaving the scattered light from Galactic cirrus intact, with the assumptions of (a) the underlying background has little power on small scales, and (b) the Galactic cirrus in the field is optically thin on large scales; (2) Censoring of frames contaminated by anomalously prominent wings in the wide-angle point-spread function; and (3) Incorporation of spatial covariance in image stacking that controls the local background consistency. We demonstrate these methods using example datasets obtained with the Dragonfly Telephoto Array, but these general techniques are prospective to be applied to improve sky models in data obtained from other wide-field imaging surveys, including those from the upcoming Vera Rubin Telescope., Comment: 28 pages, 17 figures, accepted for publication in ApJ

Characterising the prevalence and properties of faint active galactic nuclei (AGN) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H$\alpha$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 H$\alpha$ lines at $z=4.2-5.5$ that have broad components with line widths from $\sim1200-3700$ km s$^{-1}$, contributing $\sim30-90$ % of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses $\sim10^{7-8}$ M$_{\odot}$. In the UV luminosity range M$_{\rm UV}=-21$ to $-18$, we measure number densities of $\approx10^{-5}$ cMpc$^{-3}$. This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions. Yet, such AGN are found in only $<1$ % of star-forming galaxies at $z\sim5$. The SMBH mass function agrees with large cosmological simulations. In two objects we detect narrow red- and blue-shifted H$\alpha$ absorption indicative, respectively, of dense gas fueling SMBH growth and outflows. We may be witnessing early AGN feedback that will clear dust-free pathways through which more massive blue quasars are seen. We uncover a strong correlation between reddening and the fraction of total galaxy luminosity arising from faint AGN. This implies that early SMBH growth is highly obscured and that faint AGN are only minor contributors to cosmic reionization., Comment: 26 pages, 20 figures. Accepted for publication in ApJ. Various improvements and updated citations. Main additions: colours (Fig 11), broad Halpha LF (Fig 18, Sect 5.2), comparisons to simulations (Sect 5.3)

The first JWST data revealed an unexpected population of red galaxies that appear to have redshifts of $z\sim 7-9$ and high masses of $M_*$ $\sim$ 10$^{10}$ M$_{\odot}$ (Labb\'e et al. 2023). Here we fit S\'ersic profiles to the F200W NIRCam images of the 13 massive galaxy candidates of Labb\'e et al., to determine their structural parameters. Satisfactory fits were obtained for nine galaxies. We find that their effective radii are extremely small, ranging from $r_{\rm e}\sim 80$ pc to $r_{\rm e} \sim 300$ pc, with a mean of $\langle r_{\rm e}\rangle \approx 150$ pc. For their apparent stellar masses, the galaxies are smaller than any other galaxy population that has been observed at any other redshift. We use the fits to derive circularized three-dimensional stellar mass profiles of the galaxies, and compare these to the mass profiles of massive quiescent galaxies at $z\sim$2.3 and nearby elliptical galaxies. We find that, despite the high redshift galaxies having $10-20$ times smaller half-light radii, the central stellar densities are comparable to those of their putative descendants at later times. The most straightforward interpretation is that the dense compact inner regions of the most massive ellipticals today were already in place $\sim 600$ Myr after the Big Bang. We caution that the redshifts and masses of the galaxies remain to be confirmed, and that the complex NIRCam point spread function is not yet fully characterized., Comment: Accepted for publication in ApJ Letters

We present the results of a search for high-redshift ($z>9$) galaxy candidates in the JWST UNCOVER survey, using deep NIRCam and NIRISS imaging in 7 bands over $\sim45$ arcmin$^2$ and ancillary HST observations. The NIRCam observations reach a $5-\sigma$ limiting magnitude of $\sim 29.2$ AB. The identification of high$-z$ candidates relies on a combination of a dropout selection and photometric redshifts. We find 16 candidates at $99$. A few galaxies at $z\sim9-10$ appear to show a clear Balmer break between the F356W and F444W/F410M bands, which helps constrain their stellar mass. We estimate blue UV continuum slopes between $\beta=-1.8$ and $-2.3$, typical for early galaxies at $z>9$ but not as extreme as the bluest recently discovered sources. We also find evidence for a rapid redshift-evolution of the mass-luminosity relation and a redshift-evolution of the UV continuum slope for a given range of intrinsic magnitude, in line with theoretical predictions. These findings suggest that deeper JWST observations are needed to reach the fainter galaxy population at those early epochs, and follow-up spectroscopy will help better constrain the physical properties and star formation histories of a larger sample of galaxies., Comment: Submitted to MNRAS

Observations across the electromagnetic spectrum of radiative processes involving interstellar dust -- emission, extinction, and scattering -- are used to constrain the parameters of dust models and more directly to aid in foreground removal of dust for extragalactic and cosmology observations. The more complementary observations, the better. Here, we quantify the relationship between scattered light and thermal emission from dust in a diffuse (cirrus) intermediate latitude cloud, Spider, using data from the Dragonfly Telephoto Array and the Herschel Space Observatory. A challenge for optical observations of faint cirrus is accurate removal of a contaminating spatially varying sky background. We present a technique to analyse two images of the same cirrus field concurrently, correlating pixel values to capture the relationship and simultaneously fitting the sky background as a complex non-correlating additive component. For the Spider, we measure a $g-r$ color of 0.644$\pm 0.024$ and a visible wavelength to 250 $\mu$m intensity ratio of $10^{-3} \times (0.855 \pm0.025)$ and $10^{-3} \times (1.55\pm0.08)$ for $g$ and $r$-band respectively. We show how to use any dust model that matches the thermal dust emission to predict an upper limit to the amount of scattered light. The actual brightness of the cirrus will be fainter than this limit because of anisotropic scattering by the dust combined with anisotropy of the incident interstellar radiation field (ISRF). Using models of dust and the ISRF in the literature we illustrate that the predicted brightness is indeed lower, though not as faint as the observations indicate., Comment: 12 pages, 6 figures, The Astrophysical Journal

This Research Note presents VLT B-band imaging of a candidate runaway supermassive black hole that was recently discovered in HST/ACS F606W+F814W imaging. The ACS data show an extremely thin, linear feature at z=0.964 that points toward a compact galaxy at the same redshift. There is a gap between the feature and the compact galaxy, which means that the proposed causal connection between the two objects is not definitive. We show here that the linear feature extends all the way to the compact galaxy in the B-band, with no gap. The B-band morphology is difficult to reconcile with models where the compact galaxy and the linear feature are independent objects, and in particular with the proposal of Sanchez Almeida et al. (2023) that the linear feature is an edge-on disk galaxy., Comment: To appear in RNAAS. A counterpoint to arXiv:2304.12344, although the debate will undoubtably continue. Upcoming data from the Cycle 30 program HST-GO-17301 should provide a definitive answer to this particular question