We report the spectroscopic confirmation of a massive inlineformula mmlmath overflowscrollmmlmilogmmlmimmlmo stretchyfalsemmlmommlmsubmmlmrowmmlmiMmmlmimmlmrowmmlmrowmmlmommlmommlmrowmmlmsubmmlmrowmmlmo stretchytruemmlmommlmrowmmlmsubmmlmrowmmlmiMmmlmimmlmrowmmlmrowmmlmommlmommlmrowmmlmsubmmlmo stretchyfalsemmlmommlmommlmommlmn10.34mmlmnmmlmsubsupmmlmrowmmlmommlmommlmrowmmlmrowmmlmn0.07mmlmnmmlmrowmmlmrowmmlmn0.06mmlmnmmlmrowmmlmsubsupmmlmath inlineformula Hubble Space Telescopedark m SUBF150WSUB  m SUBF444WSUB  3.6 quiescent galaxy at z SUBspecSUB  3.97 in the UNCOVER survey. NIRSpecPRISM spectroscopy and a nondetection in deep Atacama Large Millimetersubmillimeter Array imaging surprisingly reveals that the galaxy is consistent with a low lt10 M SUBSUB yrSUP1SUP star formation rate SFR despite evidence for moderate dust attenuation. The F444W image is well modeled with a twocomponent Srsic fit that favors a compact r SUB e SUB  200 pc n  2.9 component and a more extended r SUB e SUB  1.6 kpc n  1.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 SFR in the center of UNCOVER 18407 require both significant A SUB v SUB  1.4 mag dust attenuation and a stellar massweighted age of 900 Myr implying 50 of the stars in the core already formed by z  7.5. Using spatially resolved annular masstolight measurements enabled by the galaxys moderate magnification inlineformula mmlmath overflowscrollmmlmimmlmimmlmommlmommlmn2.12mmlmnmmlmsubsupmmlmrowmmlmommlmommlmrowmmlmrowmmlmn0.01mmlmnmmlmrowmmlmrowmmlmn0.05mmlmnmmlmrowmmlmsubsupmmlmath inlineformula to reconstruct a radial mass profile from the bestfitting twocomponent Srsic model we infer a total massweighted inlineformula mmlmath overflowscrollmmlmsubmmlmrowmmlmirmmlmimmlmrowmmlmrowmmlmieffmmlmimmlmrowmmlmsubmmlmommlmommlmn0.74mmlmnmmlmsubsupmmlmrowmmlmommlmommlmrowmmlmrowmmlmn0.17mmlmnmmlmrowmmlmrowmmlmn0.22mmlmnmmlmrowmmlmsubsupmmlmath inlineformula kpc and log inlineformula mmlmath overflowscrollmmlmo stretchyfalsemmlmommlmsubmmlmi mathvariantnormalmmlmimmlmrowmmlmn1mmlmnmmlmspace width0.33emmmlmspacemmlmikpcmmlmimmlmrowmmlmsubmmlmspace width0.33emmmlmspacemmlmo stretchyfalsemmlmommlmsubmmlmiMmmlmimmlmommlmommlmsubmmlmspace width0.25emmmlmspacemmlmsupmmlmikpcmmlmimmlmrowmmlmommlmommlmn2mmlmnmmlmrowmmlmsupmmlmo stretchyfalsemmlmommlmo stretchyfalsemmlmommlmommlmommlmn9.65mmlmnmmlmsubsupmmlmommlmommlmn0.15mmlmnmmlmn0.12mmlmnmmlmsubsupmmlmath inlineformula. The early formation of a dense low SFR and dusty core embedded in a less attenuated stellar envelope suggests an evolutionary link between the earliestforming 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.
['High-redshift galaxies', 'Galaxy quenching', 'Galaxy evolution', 'Quenched galaxies', 'Post-starburst galaxies', 'Near infrared astronomy', 'Interstellar dust', '734', '2040', '594', '2016', '2176', '1093', '836', 'Astrophysics - Astrophysics of Galaxies']
UNCOVER NIRSpecPRISM Spectroscopy Unveils Evidence of Early Core Formation in a Massive Centrally Dusty Quiescent Galaxy at z SUBspecSUB  3.97
{'UNCOVER NIRSpec/PRISM Spectroscopy Unveils Evidence of Early Core Formation in a Massive, Centrally Dusty Quiescent Galaxy at z spec = 3 . 97': "David J. Setton, 1, ∗ Gourav Khullar, 2 Tim B. Miller, 3 Rachel Bezanson, 2 Jenny E. Greene, 1 Katherine A. Suess, 4, † Katherine E. Whitaker, 5, 6 Jacqueline Antwi-Danso, 7, 8, 5, ‡ Hakim Atek, 9 Gabriel Brammer, 10 Sam E. Cutler, 5 Pratika Dayal, 11 Robert Feldmann, 12 Seiji Fujimoto, 13, † Lukas J. Furtak, 14 Karl Glazebrook, 15 Andy D. Goulding, 1 Vasily Kokorev, 11 Ivo Labbe, 16 Joel Leja, 17, 18, 19 Yilun Ma ( 马 逸 伦 ), 1 Danilo Marchesini, 20 Themiya Nanayakkara, 15 Richard Pan, 20 Sedona H. Price, 2 Jared C. Siegel, 1, 21 Heath Shipley, 22 John R. Weaver, 5 Pieter van Dokkum, 23 Bingjie Wang ( 王 冰 洁 ), 17, 18, 19 and Christina C. Williams 24  \n1 Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA  \n2 Department of Physics and Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA  \n- 3 Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA), Northwestern University, 1800 Sherman Ave, Evanston IL 60201, USA  \n4 Kavli Institute for Particle Astrophysics and Cosmology and Department of Physics, Stanford University, Stanford, CA 94305, USA 5 Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA  \n6 Cosmic Dawn Center (DAWN), Denmark  \n- 7 David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, 50 St George Street, Toronto, ON M5S 3H4, Canada 8 Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St George Street, Toronto, ON M5S 3H4, Canada 9 Institut d'Astrophysique de Paris, CNRS, Sorbonne Universit'e, 98bis Boulevard Arago, 75014, Paris, France  \n10  \n15  \nCosmic Dawn Center (DAWN), Niels Bohr Institute, University of Copenhagen, Jagtvej 128, København N, DK-2200, Denmark  \n11  \nKapteyn Astronomical Institute, University of Groningen, 9700 AV Groningen, The Netherlands  \nDepartment of Astrophysics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland  \n13  \nDepartment of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA  \n14 Physics Department, Ben-Gurion University of the Negev, P.O. Box 653, Be'er-Sheva 84105, Israel  \nCentre for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, VIC 3122, Australia  \n16  \nCentre for Astrophysics and Supercomputing, Swinburne University of Technology, Melbourne, VIC 3122, Australia  \n- 17 Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA  \n18  \nInstitute for Computational & Data Sciences, The Pennsylvania State University, University Park, PA 16802, USA  \n19  \nInstitute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA  \n20  \nDepartment of Physics & Astronomy, Tufts University, MA 02155, USA  \n21 NSF Graduate Research Fellow  \n22 Department of Physics, Texas State University, San Marcos, TX 78666, USA  \n23 Department of Astronomy, Yale University, New Haven, CT 06511, USA  \n24 NSF's National Optical-Infrared Astronomy Research Laboratory, 950 North Cherry Avenue, Tucson, AZ 85719, USA  \n(Dated: May 14, 2024)  \nSubmitted to ApJ", 'ABSTRACT': "We report the spectroscopic confirmation of a massive (log( M ⋆ /M ⊙ ) = 10 . 34 ± 0 . 06 0 . 07 ), HST-dark ( m F150W -m 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 ⊙ yr -1 ) star formation rate despite evidence for moderate dust attenuation. The F444W image is well modeled with a two component S'ersic fit that favors a compact, r e ∼ 200 pc, n ∼ 2 . 9 component and a more extended, r e ∼ 1 . 6 kpc, n ∼ 1 . 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 ∼ 1 . 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 magnifica-  \n12  \ntion ( µ = 2 . 12 ± 0 . 05 0 . 01 ) to reconstruct a radial mass profile from the best-fitting two-component S'ersic model, we infer a total mass-weighted r eff = 0 . 72 ± 0 . 15 0 . 11 kpc and log(Σ 1kpc [M ⊙ / kpc 2 ]) = 9 . 61 ± 0 . 08 0 . 10 . 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.  \nKeywords: High-redshift galaxies (734); Galaxy quenching (2040); Galaxy evolution (594); Quenched galaxies (2016); Post-starburst galaxies (2176); Near infrared astronomy (1093); Interstellar dust (836)", '1. INTRODUCTION': "The centers of the most massive quiescent galaxies in the local universe show evidence of old ( > 10 Gyr) stellar populations (Trager et al. 2000; Thomas et al. 2005; S'anchez-Bl'azquez et al. 2006a; Graves et al. 2007; Greene et al. 2013; McDermid et al. 2015) and compact structures (e.g., Shen et al. 2003; Trujillo et al. 2004; Cappellari et al. 2013). This finding indicates that their star forming progenitors formed the bulk of their massive, dense cores in the first few Gyr of cosmic time before rapidly quenching their star formation. To link this population of red and dead galaxies to their high-z counterparts, numerous studies have attempted to connect compact star forming, sub-millimeter galaxies, post-starburst, and quiescent populations near cosmic noon to their lower-z descendants (e.g., Swinbank et al. 2006; van Dokkum et al. 2015; Ma et al. 2015; Barro et al. 2017; Suess et al. 2019, 2021) and to quantify the rate of rapid quenching across cosmic time (e.g., Wild et al. 2016; Pattarakijwanich et al. 2016; Rowlands et al. 2018; Belli et al. 2019; Park et al. 2022; Setton et al. 2023; Gould et al. 2023; Park et al. 2024). The emerging picture is that the most massive galaxies in the universe formed a significant portion of their stellar mass in the very early universe ( z ≳ 3) in dense, compact regions that evolve into the cores of present-day ellipticals (e.g., Bezanson et al. 2009), but precisely how early this formation occurs remains an open question. As such, there has been a concerted effort to find and characterize quiescent galaxies at high-redshift to understand when they form and how they quench their star formation.  \nPrior to the launch of the James Webb Space Telescope (JWST), the majority of these efforts were focused on using ground-based medium-band NIR photometry  \nto select for strong breaks at z > 3 (e.g., Whitaker et al. 2011; Straatman et al. 2014, 2016; Zaidi et al. 2024). These surveys discovered systems with a wide range of SED shapes, hinting at variety in the age and dust content of quiescent systems at this crucial epoch. However, essentially all successful spectroscopic followup programs of z > 3 quiescent galaxies have observed young, blue 'post-starburst' spectra (e.g., Marsan et al. 2015; Glazebrook et al. 2017; Schreiber et al. 2018a,b; Tanaka et al. 2019; Forrest et al. 2020; Saracco et al. 2020; Valentino et al. 2020; D'Eugenio et al. 2020; Kalita et al. 2021; Man et al. 2021; Tanaka et al. 2023; Kakimoto et al. 2024).  \nIn the absence of selection effects, the discovery of exclusively recently-quenched galaxies at ≳ 3 would indicate that the first quiescent galaxies formed the majority of their stellar mass between z =3-5 before rapidly quenching. However, continuum observations of poststarburst systems, which are brighter at fixed stellar mass due to their younger ages, at these redshifts were straining even the largest ground-based observing facilities. It has therefore been suggested that some of the apparent lack of old, quiescent systems at z ≳ 3 may be biased by the inability of ground-based facilities to robustly discover older stellar populations with higher mass-to-light ratios (Forrest et al. 2020; Antwi-Danso et al. 2023a). As such, a hidden red-and-dead population that formed in the first ∼ 1 Gyr of cosmic time (as suggested by star formation history and chemical measurements of lower-z systems that appeared maximally old, e.g., Thomas et al. 2005; Kriek et al. 2016) could still be found in deeper observations.  \nEarly JWST data has quickly confirmed this suggestion, with photometric (e.g., Carnall et al. 2023a; Valentino et al. 2023; Alberts et al. 2023) and spectroscopic (e.g., Carnall et al. 2023b; Nanayakkara et al. 2024; Glazebrook et al. 2024; de Graaff et al. 2024; Carnall et al. 2024) observations at z > 3 finding massive, quiescent galaxies with old ages that could only have resulted from extremely early-onset star formation, with  \nstar formation history fitting of deep spectroscopic samples at z ∼ 2 providing archeological evidence of similar populations (Slob et al. 2024; Park et al. 2024). Characterizing such early-forming systems has become especially important given that high-redshift studies are turning up candidate examples of rapid mass growth in the z > 7 universe (Labb'e et al. 2023; Boyett et al. 2023; Tacchella et al. 2023)-the star formation histories and structures of the earliest-forming quiescent galaxies provide an important constraint on how quickly galaxies can have assembled. As such, the detailed study of massive, old quiescent galaxies, in addition to younger co-eval post-starburst galaxies (e.g., Belli et al. 2023; D'Eugenio et al. 2023), in the first ∼ 2 Gyr of cosmic time, promises to illuminate the physical drivers of rapid mass assembly and quenching in the early universe.  \nIn this work, we present a unique example of JWSTidentified early forming galaxy-UNCOVER 18407 1 -a moderately lensed ( µ ∼ 2) massive (log( M ⋆ /M ⊙ ) ∼ 10 . 4) quiescent galaxy at z spec = 3 . 97 that was observed in seven band NIRCAM imaging and with deep PRISM spectroscopy as part of the UNCOVER Cycle 1 Treasury Program (Bezanson et al. 2022a). In contrast with the majority of spectroscopically confirmed quiescent sources at high-redshift, UNCOVER 18407 has an extremely red ( m F150W -m F444W ∼ 3 . 6) spectral energy distribution, suggesting either an old stellar population, significant dust attenuation, or some combination of the two. We leverage the high-spatial resolution of JWST complete with the added bonus from the lensing boost of Abell 2744 (Furtak et al. 2023) - to apply state-of-theart methods (e.g., Suess et al. 2019; Akhshik et al. 2020; Jafariyazani et al. 2020; Miller et al. 2023; Akhshik et al. 2023) to perform the first spatially resolved analysis of an HST-dark quiescent galaxy.  \nThis paper is laid out as follows. In Section 2, we briefly describe the UNCOVER Treasury Program (PIs: Labb'e and Bezanson, JWST-GO-2561 Bezanson et al. 2022a) and the imaging and spectroscopy we will analyze to constrain the stellar populations and structure of UNCOVER 18407. In Section 3, we describe the inference of the imaging, photometric, and spectroscopic data to infer the star formation history and structural parameters of this galaxy. In Section 4, we measure the  \nspatially resolved and global properties of UNCOVER 18407. In Section 5, we discuss the implications of our discovery of a massive, dusty, compact quiescent system at z spec = 3 . 97. Finally, in Section 6, we summarize our results and propose future avenues to build statistical samples of dusty, old quiescent galaxies at high-z to understand their cosmological significance.  \nThroughout this work, we adopt the best-fit cosmological parameters from the WMAP 9 year results (Hinshaw et al. 2013): H 0 = 69 . 32 km s -1 Mpc -1 , Ω m = 0 . 2865, and Ω Λ = 0 . 7135, utilize a Chabrier initial mass function (Chabrier 2003), and quote AB magnitudes. All effective radii ( r e ) values measured denote the semi-major axis of the best-fitting 2D S'ersic ellipse, and are corrected for the effect of lensing. Unless otherwise mentioned, all star formation rates and stellar masses are also corrected for the effect of lensing.", '2.1. Imaging and Photometry': "All JWST/NIRCam images used in this work were obtained as a part of the UNCOVER Treasury Program (PIs: Labb'e and Bezanson, JWST-GO-2561 Bezanson et al. 2022a), which took seven-band (NIRCam F115W, F150W, F200W, F277W, F356W, F410W, and F444W with AB point source depths 29-30) imaging of the Abell 2744 cluster, and Medium Bands, Mega Science (PI: Suess, JWST-GO-4111 Suess et al. 2024), which obtained NIRCam F070W and F090W broadband imaging in addition to all eleven medium bands (F140M, 162M, 182M, F210M, F250M, F300M, 335M, 360M, 430M, 460M, and F480M). The full reductions of the imaging are described in Bezanson et al. (2022a) and Suess et al. (2024), and we utilize the v7.2 imaging mosaics hosted on the DAWN JWST Archive (Valentino et al. 2023) 2 .  \nThe construction of photometric catalogs based on a combination of the full JWST/UNCOVER mosaics and archival HST imaging is described in Weaver et al. (2024). In this work, we utilize the DR3 catalog that consists of photometry extracted using circular apertures on PSF-matched images with apertures selected to be closest to the total size of the galaxy in the segmentation map. We also adopt rest frame colors measured using photometric SED fits as detailed in Wang et al. (2024a), which fits all available photometry in the catalog using the Prospectorβ model (Wang et al. 2023). To account for the lensing of the foreground Abell 2744 cluster, we correct for magnifications calculated with the  \nFigure 1. Observed m F 444 w versus m F 150 W -m F 444 W (left, Weaver et al. 2024), rest frame U -V versus V -J colors (center), and rest frame u s -g s and g s -i s colors (right, using the marginalized colors from Wang et al. 2024a) for the UNCOVER sample ( z phot > 2 and signal-to-noise in F150W and F444W > 3, black). Points with m F 150 W -m F 444 W > 3 are indicated as red, and UNCOVER 18407 is indicated as a purple star. The quiescent UVJ box and ugi s boxes are shown as in Williams et al. (2009) and Antwi-Danso et al. (2023b) respectively. UNCOVER 18407 was targeted for its unique combination of an extremely red color and quiescent photometric solution.  \n<!-- image -->  \nFigure 2. Top: NIRCAM UNCOVER (Bezanson et al. 2022a) and MegaScience (Suess et al. 2024) imaging in increasing order of wavelength for our target, with neighboring sources masked following (Weaver et al. 2024). The 0.48' (red) and 1.4' (dark blue) diameter apertures are also indicated as circles. Bottom: (Left): The extracted spectra from the three visits to UNCOVER 18407 (pre-photometric flux calibration, MSA2: green, MSA3: orange, MSA4: blue), along with the central 0.48' (red) and global 1.4' (dark blue) diameter photometry. The SED shape in the central region and in the observed spectra is considerably flatter redward of the break than in the global photometry. Key Balmer series features are labeled with dashed lines. Weak H α emission is detected in all three exposures (including MSA2, which misses the direct center of the galaxy), and the Balmer series is otherwise seen in absorption only. (Right): An RGB image of the source with the microshutter arrays labeled. Because the source is extended, we utilize a global background subtraction from empty slits.  \n<!-- image -->  \ninternal lensing model with updated spectroscopic redshifts following the methods of Furtak et al. (2023).", '2.2. Spectroscopic Targeting of UNCOVER 18407': "In addition to seven-band NIRCAM imaging, the JWST UNCOVER Treasury Program includes NIRSpec/PRISM spectroscopy (0 . 6 µ m < λ < 5 . 3 µ m) of ∼ 700 sources selected for a variety of science cases. Two stated goals of the UNCOVER program are the study of obscured/old 'HST-dark' sources (e.g., Nelson et al. 2023a; Barrufet et al. 2023; P'erez-Gonz'alez et al. 2023; Barrufet et al. 2024) and the study of the first quiescent galaxies at high-redshift; UNCOVER 18407 is a unique source within the survey that was targeted for neatly meeting both the above criteria.  \nIn Figure 1, we show the full UNCOVER dataset with z phot > 2 and signal-to-noise in F150W and F444W > 3 (black points) in observed m F150W -m F444W versus m F444W color magnitude space (left), in addition to restframe UVJ (center) and ( ugi ) s (right) plots with commonly applied quiescent selection criteria from Williams et al. (2009) and Antwi-Danso et al. (2023b), respectively. We additionally highlight the 15 galaxies with m F150W -m F444W (red), effectively selecting the reddest bright sources in the sample that would have been undetectable with the Hubble Space Telescope. The vast majority of these galaxies are best fit in Wang et al. (2024a) as dusty star forming galaxies at z ≳ 2 . 5, as suggested by their positions in the upper right in the rest frame UVJ and ugi plots. The dusty star-forming nature of these sources is supported by detections in ALMA 1.2 mm dust continuum observations in several cases (Fujimoto et al. 2023a; Mu˜noz Arancibia et al. 2023), with morphological analyses presented in Price et al. (2023).  \nHowever, one galaxy, UNCOVER 18407 (shown as a star), is unique in that it lies in the quiescent regions of UVJ and ( ugi ) s space. In the Wang et al. (2024a) catalog, it is best fit with z phot = 4 . 06 ± 0 . 13 0 . 15 , log( M ⋆ /M ⊙ )= 10 . 29 ± 0 . 06 0 . 05 , log(sSFR [yr -1 ])= -11 . 03 ± 0 . 64 1 . 12 , and a mass weighted age of 0 . 92 ± 0 . 08 0 . 12 Gyr, marking it as potentially one of the earliest-forming quiescent galaxies ever identified. The best fit model of this galaxy also suggests significant dust attenuation ( A v = 0 . 49 ± 0 . 17 0 . 14 ). The confluence of these factors with its magnification ( µ ∼ 2) made it a high priority for spectroscopic follow up to confirm its quiescent and old nature. A full presentation of the photometric identification of quiescent galaxies in the UNCOVER field will be presented in Khullar et al. in preparation.  \nUNCOVER 18407 was observed in three of the seven micro-shutter array (MSA) configurations of the UNCOVER program between July 31 and August 2, 2023.  \nMSA designs were optimized to target a wide range of target types based on the required depth of observation and source priority. Because it is extended, UNCOVER 18407 was assigned five slitlets, and observations were conducted with a 2-POINT-WITH-NIRCam-SIZE2 dither pattern. The total exposure times in each of the MSA2, MSA3, and MSA4 exposures that targeted UNCOVER 18407 were 2.6 hr, 2.3 hr, and 4.4 hr, respectively.", '2.3. Spectroscopic Reduction': "The spectroscopic reduction for each of the individual MSA exposures of UNCOVER 18407, which will further be described in Price et al. in preparation, is as follows. Stage 2 data products were downloaded from MAST reduced with msaexp (v0.6.10 Brammer et al. 2022). msaexp masks artifacts (including snowballs), and corrects for 1/f noise. Individual slits are identified and WCS solutions are applied and data are flat fielded. Because the source is extended, we do not elect to background subtract using vertically shifted spectra, but instead utilize empty slits to avoid self-subtraction. The spectra are best fit by msaexp with z spec = 3 . 97, corresponding to a magnification µ = 2 . 12 ± 0 . 05 0 . 01 using the models in Furtak et al. (2023). Throughout this work, we correct all physical measurements of mass, size, and star formation rate for the effect of this magnification.  \nIn Figure 2, we show the extracted spectra for the three MSA configurations that targeted UNCOVER 18407, in addition to photometry from two physical aperture diameters: 1.4' (purple) and 0.48' (red). Additionally, we present an RGB image of our target with the 5-slit MSA configuration for each exposure colored as in the spectrum. The spectra exhibit a range in normalization due to differences in alignment over the galaxy; for example, MSA2 is off-center enough that it misses the core of the galaxy between two shutters, and as such is significantly fainter than the other two exposures that sample that part of the galaxy. However, all three exposures, regardless of their orientation, exhibit strong Balmer breaks and weak Hα emission (see Section 3.2), lending strong credence to the quiescent photometric solution for UNCOVER 18407.  \nOur primary science goal using the spectrum is to place limits on the star formation rate of UNCOVER 18407 via the Hα line strength. We therefore elect to coadd only the red end of the spectrum to create a stacked one-dimensional spectrum of the Hα /[NII]/[SII] complex. Additionally, it is clear in Figure 2 that the MSA2 slit is missing the core of the galaxy, as evidenced by the significantly lower flux in the spectrum. As such, we elect to combine only MSA3 and MSA4 in our coadd.  \nTo calibrate the spectrum, we fit a 2nd order polynomial to synthetic photometry measured from the MSA3 and MSA4 spectra in addition to the UNCOVER DR2 0.48' diameter aperture photometry Weaver et al. (2024), which visually provides the best match to the normalization of the MSA3 and MSA4 spectra, for the F200W, F277W, F356W, F410W, and F444W imaging, in addition to the medium bands spanning that same wavelength range. We do the same to synthetic photometry in the same bands measured on each of the individual MSA exposures, and we multiply the observed spectra and error vectors by the ratio of these polynomials. We then take the inverse variance weighted average of the MSA3 and MSA4 spectra to produce the flux calibrated data product that we will utilize for this analysis.", '3. MEASURING GALAXY PROPERTIES': "The differing shapes of the spectral energy distribution within different apertures (see Figure 2, Weaver et al. 2024) hints at radially varying stellar populations or dust attenuation that we can study due to the excellent resolution and broad wavelength coverage of NIRCam in conjunction with the ∼ √ 2 boost from the gravitational lensing of Abell 2744. We therefore elect to allow radial variation in the stellar populations by adopting and modeling annular photometry following the geometry of the best-fitting S'ersic models. In this section, we demonstrate the extraction of annular photometry, our approach to radial stellar population fitting with priors informed by the weak observed Hα emission.", '3.1. Light-Weighted Structure': "To measure the structure of UNCOVER 18407, we use pysersic (Pasha & Miller 2023) to perform S'ersic fits to the 2D image of the galaxy. Because the galaxy clearly exhibits strong color gradients (see Figure 2), in addition to a traditional single component S'ersic fit, we also perform two component modeling analysis, one with two S'ersic profiles and another with a S'ersic and point source component. For these, we allow the central locations of the components to vary independently. To ensure accurate sampling, for each model we reparameterize the posterior following the procedure in Hoffman et al. (2019) using a Block Neural Auto regressive flow. (De Cao et al. 2020) We then sample the posterior using a No U-turn sampler with two chains 1000 warm-up and 1000 sampling steps (Hoffman et al. 2014; Phan et al. 2019). We ensure that the results of the sampling are robust by computing the ˆ r and effective sample size metrics which are < 1 . 01 and > 300 for all variables in all of the models (Vehtari et al. 2021). In our fitting, we adopt the empirical PSF measured in  \nWeaver et al. (2024) and smooth all intrinsic profiles to the resolution of the F444W image before fitting. The results of this fitting are presented in Table 1. Using χ 2 as a goodness-of-fit metric, it is clear that the single S'ersic model is disfavored relative to the two component models. The difference between the two S'ersic and S'ersic plus PSF models is more marginal, given how compact the central component is.  \nIn Figure 3, we demonstrate visually the preference for a two component model by showing the F444W image of the galaxy, with ellipses overlaid showing the best fitting galaxy shapes at r e for the single component fit (purple) and the two component fits (green, compact as solid and extended as dashed). We correct the axis ratio for the effect of the PSF as:  \nq corr = √ ( qr e ) 2 + r 2 PSF HWHM r 2 e + r 2 PSF HWHM (1)  \nfollowing Suess et al. (2019) and using the most compact effective radius as our benchmark r e , and r PSF HWHM = 0 . 0725'. In the middle panel of Figure 3, we show the one-dimensional surface brightness profile (extracted using the best fitting compact component in the two component model) of the galaxy in black, and of the one- and two- component S'ersic models again as purple and green respectively. The two component model, which consists of a ∼ 200 pc, S'ersic n ∼ 2.9 compact profile and a ∼ 1 . 6 kpc, n ∼ 1.7 extended profile, is better able to capture both the compact galaxy core and the extended exponential envelope, whereas the single S'ersic fit compromises by fitting a ∼ 1 kpc n ∼ 4 . 5 profile that overshoots the wings and undershoots the galaxy center, as evidenced by the residuals. We note that the two components are physically offset from one another by ∼ 1 pixel (0.04' ∼ 0.2 kpc); while this difference is well within the PSF FWHM, it is a potential hint at a disturbed morphology or the presence of strong dust lanes. In the third panel, we show the F150W - F444W surface brightness profiles (corresponding to just blueward of the break and redward of Hα at z spec = 3 . 97) measured again using the geometry of the compact S'ersic profile on images that have been matched to the F444W PSF (Weaver et al. 2024), demonstrating that the color of UNCOVER 18407 is significantly redder in the central kpc than it is in the outskirts. We note again that the two component S'ersic model is only marginally preferred over the S'ersic +point source model. However, given that the two S'ersic fit favors a resolved central component with a clear preference for slight intrinsic elongation, we carry forward with this profile as our benchmark for annular extraction, noting that none of  \nFigure 3. (Left): The F444W image of UNCOVER 18407. The plotted ellipses show the single (magenta) and two (green, component 1 solid, component 2 dashed) best fitting S'ersic models to this image at r e , accounting for the PSF smearing of the axis ratio as described in Section 3.3. (Center): The 1D surface brightness profiles of the observed galaxy (black) and the best fitting PSF-convolved one- and two- component models, along with residuals. We indicate the PSF half-width-half-maximum in the F444W image as a dashed vertical line. The two component model clearly does a better job at capturing the exponentially declining wings that are missed by the single component fit. The effective radii of the S'ersic components are labeled, with solid and dashed lines again corresponding to the compact and extended components of the two component model. (Right): The F150W-F444W color gradient, measured using the best fitting sersic ellipse on the F444W PSF-convolved images. The two component structure can be seen in a peaky central red excess that flattens outside the central kpc.  \n<!-- image -->  \nTable 1. Best fitting S'ersic parameters to the F444W image  \nNote -The best fitting structural parameters to the F444W image of UNCOVER 18407, as discussed in Section 3.1. The physical effective radii presented are corrected for the effect of gravitational lensing by dividing by √ µ .  \nour conclusions are strongly impacted by the choice to extract photometry with elliptical versus circular annuli.", '3.2. Hα and FIR Star Formation Rate Constraints': 'In order to constrain the total star formation rate, we require measurements of the Hα line luminosity. However, at the resolution of the NIRSpec/PRISM, Hα and the [NII] doublet are completely blended and appear as a single line with a joint rest-frame equiva-  \nnt width of 21 ˚ A . As such, in order to measure the line fluxes and associated errors, we jointly fit the spectrum at 0 . 63 µm < λ rest < 0 . 68 µm , modeling Hα along with the [NII] and [SII] doublets as Gaussians with tied σ (which, in this narrow region of wavelength, is a nuisance parameter that accounts for both the instrumental resolution and any intrinsic line spread) and a single free redshift. We fix the flux ratio of the [NII] doublets at  \nFigure 4. The stellar continuum-subtracted (see Section 3.2) spectrum of UNCOVER 18407 (black) in the region surrounding Hα /[NII]/[SII], along with the best fit (green). Also shown is the 1 σ confidence interval for Hα (blue) along with the best fitting Hα , [NII], and [SII] lines. At the PRISM resolution, we are not able to rule out a line that is totally dominated by Hα .  \n<!-- image -->  \nKriek et al. (2016) dust law with a dust index of 0, this corresponds to a limit of ∼ 2 . 5 -13 M ⊙ yr -1 .  \n1:3, and allow the ratio of [SII] λ 6716 /λ 6731 to vary from 0.4375 to 1.4484 following Sanders et al. (2016).  \nIn order to account for the continuum, including Hα absorption, before fitting for the lines we fit a simple stellar population with a free dust law, convolved to the wavelength dependent PRISM resolution, to the fluxcalibrated spectrum redward of 2 µ m. We mask Hα in this fit but otherwise leave the spectrum free, and we double the error vector to account for uncertainties in the small scale flux calibration. We subtract the best fitting continuum model from our spectrum, and we then utilize emcee (Foreman-Mackey et al. 2013) to fit the covariant posteriors of the line strengths.  \nWe show the results of this fitting in Figure 4, with the shaded region surrounding the best fit on Hα indicating the 1 σ constraints on that line ( f Hα = 5 . 8 ± 3 . 9 3 . 4 × 10 -19 erg / s / cm 2 after accounting for the magnification). At the NIRSpec/PRISM resolution, we are not able to deblend Hα and [NII] effectively, and the line is consistent with being essentially totally dominated by Hα . As such, we adopt the 3 σ Hα luminosity (flux), 2 . 1 × 10 41 erg / s (1 . 32 × 10 -18 erg / s / cm 2 ), corresponding to a star formation rate of 1.15 M ⊙ yr -1 with the assumption of 5 . 5 × 10 -42 M ⊙ yr -1 erg / s (Leitherer et al. 1999) after correcting for magnification, as the limiting dustunobscured star formation rate contained within the slit. Under the assumption of A v ∼ 1 -3 magnitudes and a  \nWe note that the Hα /[NII]/[SII] complex is the only region of the spectrum of UNCOVER 18407 with significantly detected emission lines. We see no evidence of strong [OIII], even in the MSA configurations that most strongly overlap with the center of the galaxy. However, given the strong measured color gradients that could arise from centrally concentrated dust, we cannot conclusively rule out a dustobscured AGN, and therefore cannot conclusively ascribe all Hα flux to star formation. The lack of ability to rule out AGN ionization is further bolstered by our constraints on line ratios of the detected features; we measure log([NII 6583 ] /H α )= -0 . 07 ± 0 . 54 0 . 54 and log([SII 6716 , 6731 ] /H α )= -0 . 21 ± 0 . 37 0 . 22 , spanning the full range of observed ratios in star forming galaxies and AGN (e.g., Kewley et al. 2006). This further motivates our choice to treat the Hα star formation rate as an upper limit rather than a firm measurement.  \nThe Hα limits on the total star formation rate even under the assumption of negligible [NII] contribution are also corroborated by the ALMA non-detection of the galaxy in the DUALZ survey (Fujimoto et al. 2023b). The 3 σ upper limit at 1.2-mm flux observed is 113 µ Jy. Assuming a FIR SED with a modified black body ( T dust = 35K and spectral index = 1.8) and accounting for the magnification due to lensing, the ALMA non-detection places a 3 σ upper limit of 8 . 5 M ⊙ yr -1 on the galaxy star formation rate. This limit is quite sensitive to the assumed dust temperature, but in conjunction with the weak observed Hα , it provides a strong constraint that the total star formation rate for UNCOVER 18407 should not exceed ∼ 10 M ⊙ yr -1 .', '3.3. Annular Stellar Population Fitting': "To further quantify the cause of these observed color gradients, we make use of the 7-band UNCOVER NIRCAM imaging to fit the radially varying spectral energy distribution of UNCOVER 18407. Using the best-fitting structural parameters from the compact component of the two-component S'ersic fit and the corrected axis ratio as in Section 3.1, we extract annular photometry in annular bins, with a minimum annulus size of 0.14' ( ∼ 0.7 kpc), the PSF FWHM of F444W PSF-matched images (Weaver et al. 2024). We extract annuli in five bins, out to ∼ 3 . 6 kpc ( > 3 r e in the single S'ersic fit). We enforce a signal-to-noise threshold of 20 on our photometry to account for systematic uncertainties in the models and the flux calibration (e.g., Wang et al. 2024a). We perform our SED fits using 19/20 available bands, omitting F140M as it exhibits discrepant colors with F150W de-  \nFigure 5. Annular prospector SED fits (left) and measured star formation histories (right) with 1 σ confidence intervals shown as a shaded region. Labels indicating the physical radius of the aperture being fit, except in the final panel, which shows the sum of the five annular fits. The solid and dashed lines on the star formation histories show the median and 1 σ constraints on the mass weighted age. The SED shape of the galaxy clearly becomes bluer at large physical radii, and the median mass weighted age correspondingly shifts to lower values, albeit with significant uncertainty. We note that while the SEDs are plotted in observed flux units, the star formation histories have been corrected for the effect of lensing.  \n<!-- image -->  \nspite the lack of any strong features in that part of the spectrum blueward of the break (see Figure 2).  \nWe utilize the non-parametric star formation histories in Prospector (Johnson & Leja 2017; Leja et al. 2017; Johnson et al. 2021) to fit the annular photometry in nine fixed bins of constant star formation. We set the widths of the nine star formation history bins in order of lookback time to 5 Myr, 25 Myr, 70 Myr, 100 Myr, 100 Myr, 150 Myr, and the final three bins are set to logarithmically fill the remaining time since the Big Bang, and, following the continuity prior in Leja et al. (2019), the logarithmic difference between the SFR in neighboring bins is sampled over with a Student-T distribution prior centered at 0 with a width of 0.3 and ν = 2. We utilize the Flexible Stellar Population Synthesis (FSPS) stellar population synthesis models (Conroy et al. 2009; Conroy & Gunn 2010), the MILES spectral library (S'anchez-Bl'azquez et al. 2006b) and MIST isochrones (Choi et al. 2016; Dotter 2016). We assume a Chabrier (2003) Initial Mass Function and fix the model redshift to the spectroscopic redshift. We assume the Kriek & Conroy (2013) dust law with a free A v and dust index (spanning [0,2.5] and [-1,0.4] respectively). Additionally, following Wild et al. (2020), we fix the optical depth of attenuation around young ( < 10 7 yr) stars to be double the optical depth for the rest of the stars. We leave stellar metallicity as a free parameter, and sample linearly in the range 0 . 01 Z ⊙ < Z < 2 Z ⊙ . We include nebular emission in our fits, with the ionization parameter log(U) varying in the range [-4,-1] and the gas phase metallicity fixed to the stellar metallicity. Because at z spec = 3 . 97 even the reddest NIRCam/NIRSpec does not constrain wavelengths ≳ 1 µ m, we fix the shape of the IR SED following the Draine et al. (2007) dust emission templates, with U min = 1.0, γ e = 0.01, and q PAH = 2.0. For similar reasons, we do not fit for any IR AGN contribution. We sample our posteriors using the dynesty nested sampling package (Speagle 2020).", '4.1. Can the core be maximally old and dust free?': "Before delving into the marginalized radial stellar population trends, we first present a simple demonstration that the core of UNCOVER 18407 must contain significant dust based on rest frame colors. Given the age of the universe at z = 3 . 97, there is only ∼ 1 . 5 Gyr over which any star formation could have occurred; using rest-frame colors, we posit the question of whether the extremely red SED of UNCOVER 18407 could result from a maximally old stellar population. In Figure 6 we show rest frame UVJ colors (Williams et al. 2009) for  \nFigure 6. Rest frame U-V and V-J colors for the central 0.7 kpc (light red), the full galaxy (green, using the 1.4' aperture photometry from Weaver et al. (2024)), and the outer 2.8-3.6 kpc (light blue) for UNCOVER 18407. Also plotted are simple stellar populations evolutionary tracks with solar metallicity generated with fsps (Conroy 2013) for a dust free (dark blue) and A v = 1 . 5 (dark red) stellar population. The model tracks are labeled with circles at 50 Myr, 100 Myr, 500 Myr, and 1 Gyr. The black vector shows the effect of A v = 1 for a 1 Gyr old population. The central regions of the galaxy cannot be produced by a dust-free, maximally old stellar population, and require significant dust to produce the observed colors, in contrast with the outer regions.  \n<!-- image -->  \nthe galaxy center (red), the full galaxy (green), and the galaxy outskirts (light blue), with solar metallicity simple stellar populations generated with fsps shown for dust-free (dark blue) and A v = 1 . 5 (dark red) spanning 0 Gyr to 1.6 Gyr, the age of the universe at z = 3 . 97.  \nIt is immediately evident that the center of these galaxy can only result from a simple stellar population with significant extinction; the dust free curve does not reach high enough values in U -V or V -J to produce the observed central colors. In fact, to produce such high U -V and U -J colors, a dust free, simple stellar population would need to be ∼ 10 Gyr old, an order of magnitude older than is physically possible at this redshift. In contrast, the integrated galaxy UVJ colors are still within the realm of possibility for a dust free SSP, albeit a quite old one, and cannot be produced in the presence of dust attenuation that is comparable to the galaxy center.  \nFigure 7. Photometric radial constraints on the specific star formation rate (top left, with the transition sSFR boundary 1 3 t H ( z =3 . 97) from Tacchella et al. (2022) and the quiescence boundary 1 5 t H ( z =3 . 97) from Pacifici et al. (2016) in green and red respectively), A v (top right), the mass-weighted age (bottom left, in units of lookback time), and the mass-to-light ratio (bottom right) for UNCOVER 18407. The error bars capture the 1 σ constraints on parameters. When relevant, the 1 σ constraints from the global star formation history are indicated as a grey shaded region around the median value. The fits indicate that observed color gradients largely result from centrally concentrated dust. We find very slight evidence of a negative age gradient, although the profile is consistent with flat at the 2 σ level. All together, these results manifest in a declining mass-to-light gradient.  \n<!-- image -->", '4.2. Radial dust and age gradients': 'The annular SED fitting allows us to fully probe the range of age and dust solutions that can produce the observed galaxy radial SEDs. The results of this fitting are shown in Figure 7 as a function of the semi-major axis (corrected for lensing) in blue. It is clear that the color gradients in the SED of the galaxy are predominantly driven by highly concentrated dust with A v ∼ 1 . 4, with the galaxy light significantly less attenuated outside the central kiloparsec. This color gradient, in light of the low-sSFR solutions at all radii, translates nearly directly to a mass-to-light gradient that indicates that the center of the galaxy hosts ∼ 1 . 5 as much mass per unit of light relative to the outer regions of the galaxy. There is significant uncertainty in this measurement, however,  \nowing to the uncertainty in the precise amount of attenuation in the central regions.  \nIn contrast with the clear gradient in A v , gradients in the mass-weighted age are more ambiguous. We measure the slope in t 50 to be -70 ± 110 100 Myr/kpc, and we are unable to rule out flat gradients due to the uncertainty in the mass weighted age measurements. We can, however, conclusively rule out any strong positive age gradients. There is some indication of radially declining specific star formation rate (averaged over the past 100 Myr in lookback time); while the galaxy is consistent with low rates of ongoing star formation at all radii, the star formation rate within the central 700 parsecs is constrained to be an average of 1 . 7 ± 2 . 4 1 . 2 M ⊙ yr -1 over the 100 Myr before observation, corresponding to log(sSFR 100Myr [yr -1 ]) = -9 . 7 ± 0 . 5 0 . 8 , slightly higher than at r=1-3 kpc but consistent at the 1 σ level.', '4.3. Global galaxy properties': "Ideally, in order to obtain global galaxy properties, one would jointly model the spectrum and photometry to increase constraining power on the total mass assembly of the galaxy. However, we note that the strong gradients measured in this section preclude the possibility of performing a combined spectrophotometric fit for this galaxy. Joint modeling of NIRSpec/PRISM spectra and photometry is already quite difficult due to the differential slit loss as a function of wavelength over a wide wavelength range (e.g., de Graaff et al. 2023), and with the added complication of color gradients, we are not confident we can robustly recover the shape of the spectrum and the corresponding photometric aperture after correcting for slit losses.  \nAs such, we obtain global constraints on the star formation history of the galaxy by combining the annular SED fits to construct a total star formation history for the galaxy by summing the posteriors of each individual fit. This method of determining the total star formation history is preferred over a global SED fit due to the strong dust gradients, which are smoothed over when all photometry is summed. From this summed posterior, we measure the total galaxy SFR averaged over the 100 Myr before observation to be 2 . 8 ± 3 . 5 1 . 7 M ⊙ yr -1 , well in line with the ALMA 3 σ limit of 8.5 M ⊙ yr -1 and the total limits from Hα (see Section 3.2). We constrain the total stellar mass within 3.5 kpc to be log( M ⋆ ( < 3 . 5kpc) /M ⊙ ) = 10 . 35 ± 0 . 06 0 . 06 . In conjunction with the measured star formation rate, this corresponds to a specific star formation rate of log(sSFR 100Myr [yr -1 ]) = -9 . 9 ± 0 . 4 0 . 5 , ∼ 1 dex below the star forming main sequence at z ∼ 4 (using the star forming main sequence in Shapley et al. 2023, scaled  \nfrom Speagle et al. 2014 under similar IMF assumptions). Definitions of 'quenched' vary in the literature, especially at z > 3 (e.g., Pacifici et al. 2016; Tacchella et al. 2022), but our constraints on the star formation history of UNCOVER 18407 show a strong preference for a declining star formation rate over the past ∼ 500 Myr (see Figure 5).  \nWe emphasize that there is a significant discrepancy between the constraints that come from our method of analyzing spatially resolved annuli and the constraints from fitting the entire galaxy SED with a single dust law. We perform an identical fit to the one performed in Section 3.3, but now on the catalog photometry within a 1.4' aperture. This global fit measures a total stellar mass log( M ⋆ /M ⊙ ) = 10 . 34 ± 0 . 06 0 . 07 , which is quite consistent with the total mass estimate from the summed annuli. However, the dust and star formation properties of the galaxy differ significantly as a result of the averaging of dust attenuation. The global fit measures A v = 0 . 66 ± 0 . 18 0 . 17 , almost a magnitude less than the attenuation we measure in the central annulus. Similarly, the star formation rate that we measure from the global fit is 1 . 2 ± 3 . 3 1 . 1 M ⊙ yr -1 , about a factor of 2 lower than the global star formation rate measured by summing the annular fits. This is a direct result of the ability for a dusty central region can hide significantly more star formation. While such a result was implied by the empirical F 150 W -F 444 W and rest-frame UVJ color gradients, only detailed modeling that marginalizes over the degeneracy between age, dust, and metallicity is capable of robustly correcting for such trends in mass and star formation rate estimates.  \nUsing the combined posterior for the total star formation history of the galaxy, we also constrain the global mass weighted age to be 880 ± 190 170 Myr. This corresponds to an age of the universe when 50% of the galaxy's stellar mass was formed of ∼ 700 Myr after the Big Bang ( z form = 7 . 5 ± 1 . 8 1 . 3 ), suggesting that UNCOVER 18407 formed the majority of its stellar mass within the first Gyr of cosmic time. We measure the global t 90 , the mass weighted age when UNCOVER 18407 formed 90% of its stars, to be 360 ± 100 120 Myr, indicating that the galaxy has not been forming significant stellar mass in the past few hundred Myr. Such ages appear to broadly line up with the formation of the first quiescent galaxies with log( M ⋆ /M ⊙ ) ∼ 10.3, which appear as quenched quite close to z ∼ 4 in simulations (e.g., Wellons et al. 2015; Hartley et al. 2023; Xie et al. 2024).", '4.4. Constraints on the Dust and Gas Mass': 'Because the spectrophotometric fit to the galaxy implies significant dust extinction, we also use the ALMA  \nBand 6 non-detection (see Section 3.2, Fujimoto et al. 2023b) to constrain the dust mass under the assumption of cold dust ( T dust = 20K), given that star formation does not seem to be extreme from SED fitting and Hα constraints. Using Equation 2 from Greve et al. (2012), we obtain a 3 σ limit of log( M dust /M ⊙ ) < 7.7. Taking the ALMA M dust limit and applying standard assumptions ( δ M H 2 /M dust =100, e.g., Scoville et al. 2016), this corresponds to a gas mass limit of log( M H 2 /M ⊙ ) < 9.7, which in turn corresponds to an upper limit on the molecular gas fraction M H 2 /M ⋆ ≲ 20%. These limits are quite similar to those found in other quiescent systems at similar redshift (Suzuki et al. 2022). However, we note that this upper limit is quite sensitive to the assumed dust temperature (e.g., Cochrane et al. 2022).  \nSimulations that include predictive dust physics forecast considerable scatter in the dust-to-gas mass ratios of quiescent systems, with the wide range of values spanning four orders of magnitude becoming evident by z = 2 and below (Whitaker et al. 2021a). While very few quiescent systems exist at z = 3 in the simulations, these early quenched galaxies instead appear to have more normal gas-to-dust mass ratios. So while dust masses can be extremely small in quiescent systems at z = 2 and below (e.g., Whitaker et al. 2021b; Caliendo et al. 2021), it may be that the dust reservoirs at z = 3 and above are instead still largely intact in the earliest quiescent galaxies and consistent with standard gas-to-dust mass ratios. With a sample of one galaxy, we cannot yet quantify how common this dusty phase of quiescence is. However, population studies using spatially resolved attenuation measurements as a proxy for dust richness provide an opportunity to constrain the evolution of this population as large-area JWST imaging surveys become public.', '5.1. A old, centrally dusty, quiescent galaxy at z spec = 3 . 97': "The question of how the most massive quiescent galaxies in the universe form is fundamentally linked to when they form. In the pre-JWST landscape, the empirical findings pointed to a when of ∼ z = 3 -5-the wide array of spectroscopic observations of quiescent galaxies at z ≳ 3 indicated that the only populations of dead galaxies that existed at those redshifts were 'post-starburst' galaxies that shut down in the few hundred million Myr before observation (e.g., Marsan et al. 2015; Glazebrook et al. 2017; Schreiber et al. 2018a,b; Tanaka et al. 2019; Forrest et al. 2020; Valentino et al. 2020; D'Eugenio et al. 2020; Saracco et al. 2020; Kalita et al. 2021; Kubo et al. 2021; Tanaka et al. 2023; Kakimoto et al. 2024). How-  \n, the presence of systems in the local Universe (e.g., Thomas et al. 2005) and at z ∼ 2 (e.g., Kriek et al. 2016; Beverage et al. 2023) with spectral features indicating earlier formation in conjunction with the ambiguous results of of following up 'old' galaxy candidates at z > 3 with ground-based facilities (e.g., Antwi-Danso et al. 2023a; Nanayakkara et al. 2024) suggested that a population of earlier forming quiescent galaxies existed just below the sensitivity of our best NIR facilities.  \nAs with many aspects of galaxy evolution, the advent of JWST has confirmed this suspicion. Early photometric samples have revealed a previously-hidden population of massive, quiescent galaxies at z > 3 with massweighted ages consistent with formation redshifts of ∼ 10 ( ∼ 500 Myr after the Big Bang, see Carnall et al. 2023a). A smaller number of high-z quiescent systems have also been the targets of follow-up spectroscopy, revealing formation times as early at z ∼ 11 (e.g., Carnall et al. 2023b; Nanayakkara et al. 2024; Glazebrook et al. 2024; de Graaff et al. 2024; Urbano Stawinski et al. 2024; Carnall et al. 2024), providing support to the argument that early universe star formation was more vigorous than previously thought (e.g., Labb'e et al. 2023; Boyett et al. 2023). UNCOVER 18407 represents an additional confirmed example of early-onset star formation; the mass weighted age of the center of the galaxy is consistent with z form = 7 . 6 ± 1 . 8 1 . 3 (700 ± 170 190 Myr after the Big Bang).  \nWe show this new post-JWST landscape of z > 3 quiescent galaxies in Figure 8, highlighting the new JWST samples (including UNCOVER 18407) as colored points and pre-JWST spectroscopically confirmed quiescent systems in grey. There is still considerable uncertainty in the formation redshift of UNCOVER 18407 owing to the coarse sensitivity of SED shape to age, especially in the presence of significant dust, but it is consistent at the 1 σ level with forming the majority of its mass in the first Gyr of cosmic time. This formation time is consistent with the emerging picture of earlierthan-expected massive galaxy formation. However, UNCOVER18407 adds another wrinkle to this literature, as while there has been clear evidence of spatially varying dust in photometric samples of HST-faint/dark galaxies (P'erez-Gonz'alez et al. 2023), no spectroscopically confirmed massive, quenching system at z > 3 to date shows such clear indications of strong central dust extinction ( A v > 1).  \nWhile it is commonly assumed that quiescent galaxies are dust poor, the presence of dust lanes in the centers of massive, local early type galaxies are quite often hosts to centrally concentrated dust (Ebneter & Balick 1985; Ebneter et al. 1988; van Dokkum & Franx 1995; Tomita et al. 2000; Tran et al. 2001; Laine et al. 2003;  \nFigure 8. The stellar mass versus formation redshift (or time of formation) for UNCOVER 18407 (purple star) and other JWST-identified (color) and pre-JWST (grey) spectroscopically confirmed quiescent systems at z > 3. We note a galaxy reported in Nanayakkara et al. (2024), Glazebrook et al. (2024), and Carnall et al. (2024) with two markers but show only the age and mass estimates from Nanayakkara et al. (2024) clarity. For objects that appear in both Schreiber et al. (2018a) and Nanayakkara et al. (2024), we show only the latter up-to-date JWST measurements. UNCOVER18407 is among the oldest quiescent systems ever identified, and as it was found in a very small field, it is likely that significant old and quiescent populations exist and have eluded prior detection.  \n<!-- image -->  \nLauer et al. 2005; Whitaker & van Dokkum 2008; Davis et al. 2013; Boizelle et al. 2017), indicating either that quiescent galaxy environments are able to retain significant leftover dust and gas for a few hundred Myr after quenching or that they are able to consistently replenish their central reservoirs. Such high attenuation is also seen in massive quiescent galaxies at z > 2 with red SED coverage. Pre-JWST work leveraging Spitzer to probe the red SED of massive, dusty galaxies has identified dusty quiescent galaxies at z > 1 (Martis et al. 2019). Similarly, Marchesini et al. (2014) finds that the typical A v in quiescent systems grows from ∼ 0 at z = 0 to ∼ 0 . 5 at z = 2 using ULTRAVISTA photometry, and many of the spectroscopically confirmed quiescent systems in Schreiber et al. (2018a) are fit with comparable A v . With the benefit of the deep NIR imaging capa-  \nty of JWST, van Dokkum et al. (2023) presents a massive lens that exhibits an integrated UVJ color consistent with a dusty, quiescent solution, and four of the z > 3 systems reported in Nanayakkara et al. (2024) have solutions with A v > 1. Such systems also appear to be common near dusty star forming galaxies; Alberts et al. (2023) reports that the quiescent companion in a massive star forming/quiescent pair located in a z ∼ 3 . 7 overdensity is likely dusty, and Kokorev et al. (2023) reports an average A v ∼ 1 in a quiescent system at z = 2 . 58 located 9 kpc from an extremely dusty star forming galaxy.  \nAdditionally, spatially resolved studies of quiescent systems have turned up dusty sub-regions similar to those at the center of UNCOVER 18407; Miller et al. (2022) identifies two UVJ -selected quiescent galaxies at z ∼ 2 that exhibit remarkably similar UVJ color gradients to UNCOVER 18407, and the GRISM spectroscopy lensed galaxies at similar redshifts in Akhshik et al. (2023) discovered a system with A v ∼ 4 in a single spatial bin (though notably, not at the galaxy center, both in the GRISM map and in direct dust detections, see Morishita et al. 2022). The fact that large samples of highA v quiescent galaxies have yet to be catalogued is likely in large part a testament to the power of spatially resolved measurements; the vast majority of catalogs of quiescent systems at high redshift (e.g., Schreiber et al. 2018a; Carnall et al. 2023a; Alberts et al. 2023) select based on star formation histories inferred from integrated spectral energy distributions. As demonstrated in Figure 6, the integrated color that averages over the obscured core and unobscured outer regions is far less extreme in UVJ space than the dusty core, and a fit to the global photometry of the galaxy that does not account for spatially varying dust/stellar populations results in constraints on the dust content of A v = 0 . 66 ± 0 . 18 0 . 17 , a typical level of attenuation seen in global SED measurements of quiescent galaxies at z ≳ 2 (Marchesini et al. 2014).  \nAs such, we propose that it is likely that many previously observed z > 3 quiescent galaxies, when imaged with JWST, may exhibit similar color gradients due to centrally concentrated dust. This inference that many high-z quiescent galaxies host central gas and dust, based on indirect attenuation measurements and an upper limit of log( M dust /M ⋆ ) < -2 . 69, squares with direct dust-continuum and CO measurements that find that a large variety in interstellar medium content in quiescent and post-starburst galaxies (e.g., Sargent et al. 2015; French et al. 2015; Suess et al. 2017; Gobat et al. 2018; Whitaker et al. 2021b; Belli et al. 2021; Williams et al. 2021; Bezanson et al. 2022b; Smercina et al. 2022;  \nWoodrum et al. 2022; Spilker et al. 2022; Otter et al. 2022; Wu et al. 2023; Donevski et al. 2023; Lee et al. 2024).  \nThe origins of this central gas and dust remain uncertain, but there are strong links to physical associations with mergers in local galaxies (Weaver et al. 2018). As mergers appear to be physically associated with quenching (e.g., Pawlik et al. 2016; Sazonova et al. 2021; Ellison et al. 2022; Verrico et al. 2023), dynamical effects may be allowing these galaxies to sustain a centrally concentrated interstellar medium in the wake of shutdown; alternatively, recent wet mergers may be responsible for briefly re-enriching a previously depleted ISM that is quickly removed (e.g., Akhshik et al. 2021). Given that there is evidence for some small residual star formation in the galaxy core, it is plausible that any remaining ISM could be removed by this, but significantly deeper sub-mm observations would be needed to quantify the remaining ISM reservoirs given that quiescent galaxies at these redshifts are consistently undetected or weakly detected, even in stacks (Schreiber et al. 2018b; Santini et al. 2019; Chworowsky et al. 2023).  \nIn the absence of deep sub-mm detections, we demonstrate that resolved stellar population studies present an opportunity to identify ISM-rich massive quiescent galaxies at cosmological distances by exploiting the sensitivity of imaging that spans the rest-optical to NIR to dust gradients. In future work, we will leverage the growing samples of spectroscopically confirmed quiescent systems to understand the diversity of ISM content in quiescent systems as a function of cosmic time and morphology.", '5.2. Implications for dense core formation in the first Gyr of cosmic time': "It is observationally well-established that quiescent galaxies out to z ∼ 3 differ significantly from their coeval star forming counterparts, whether in half light radius (e.g. van der Wel et al. 2014; Mowla et al. 2019; Kawinwanichakij et al. 2021; Hamadouche et al. 2022; Ji et al. 2024), half mass radius (Suess et al. 2019), or the central mass density (e.g. Bezanson et al. 2009; Fang et al. 2013; van Dokkum et al. 2015; Barro et al. 2017; Mosleh et al. 2017; Whitaker et al. 2017; Suess et al. 2021). The prevailing interpretation of these trends is that massive quiescent systems are the result of a quenching pathway resulting from rapid gas compaction (e.g., Dekel & Burkert 2014; Zolotov et al. 2015; Tacchella et al. 2016), leaving behind a dense core, with a more extended stellar disk component that quenches at later times (e.g., Costantin et al. 2021, 2022). Direct evidence for this two component assembly and core  \n1  \nFigure 9. The mass-weighted effective radius (left) and the central kiloparsec surface density (right) of UNCOVER 18407 (dark blue) as compared to z=2-2.5 star-forming (light blue) and quiescent (red) galaxies, measured in Suess et al. (2019). The shaded grey 'unphysical' region denotes the space where a given density is unattainable at a given stellar mass. After accounting for the centrally concentrated dust, UNCOVER 18407 clearly shows structural signatures typical of quiescent galaxies at cosmic noon, with a dense central core that indicates intense star formation on sub-kpc scales in the first Gyr of cosmic time.  \n<!-- image -->  \nformation on scales of ∼ 200 parcsecs have been seen as early as z ∼ 7 (Baker et al. 2023).  \nHowever, inferences about galaxy growth based on galaxy sructure are complicated by the presence of color gradients, which bias structural measurements toward tracing younger stellar populations due to outshining (Tortora et al. 2010; Wuyts et al. 2010; Suess et al. 2019; Miller et al. 2023). Because UNCOVER 18407 has broad rest-frame near-IR SED coverage and benefits from the mild lensing of Abell 2744 that results in a ∼ √ 2 boost in angular resolution, it represents an excellent test of whether the stellar mass distribution of a quenched system that formed in the first Gyr of cosmic time exhibits compact structure after accounting for differential dust attenuation and stellar populations.  \nIn order to measure the stellar mass distribution of UNCOVER 18407, we assume that the intrinsic M/L distribution of the galaxy follows a power law distribution:  \nlog(M / L) = α × log(a) + β (2)  \nfollowing Suess et al. (2019) where a is the semi-major axis of an ellipse following the distribution of the central component of the galaxy (see Table 1), α is a the slope of the power law, and β is the intercept. We forward model the effect of the PSF on this power law by projecting  \nTable 2. Mass-weighted structural measurements  \nNote -The total stellar mass, the mass weighted effective radius, and the central surface mass density of UNCOVER 18407 measured as described in Section 5.2. The physical radii presented here are corrected for the magnification by dividing by √ µ .  \nthe model onto a 2D grid following the aforementioned central component best fitting ellipses and convolving with the F444W PSF from Weaver et al. (2024). We then measure the mean M/L in the PSF-corrected annuli defined in Section 3.1 to compare to our measured M/L profile. We fit the intrinsic power law with α and β as free parameters using emcee (Foreman-Mackey et al. 2013) to constrain the possible range in slopes allowed by the data.  \nAs the uncertainty in the mass to light ratio is significantly larger than the uncertainty in the 2D structure of  \nthe two components of the galaxy, we then obtain a posterior in the intrinsic mass profile of the galaxy by multiplying the intrinsic 1D light profile (again measured using the best fitting compact ellipses of the central component applied to the best fitting two-component S'ersic model) with the posterior in the measured intrinsic M/L profiles, after fixing the M/L within the central pixel to the value at 1 pixel (0.04') and the M/L outside the largest radius to the value in the outer annulus following Suess et al. (2019). We measure the total mass, the half mass radius, and the central surface mass density by integrating this new intrinsic mass profile. The results of this fitting are shown in Table 2.  \nIn Figure 9, we show UNCOVER 18407 in massweighted size and Σ 1kpc (as measured from the intrinsic mass profile after accounting for lensing) versus stellar mass as compared to star forming (blue) and quiescent (red) galaxies at z=2-2.5 as measured using 3D-HST images in Suess et al. (2019). It is clear from both panels that the structure of UNCOVER 18407 is comparably dense to z ∼ 2 quiescent systems, and has a mass distribution that is more compact than typical star forming galaxies at similar mass. Since the bulk of the stars in UNCOVER 18407 formed in the first ∼ Gyr of cosmic time, the physical process that is responsible for the fairly rapid shutdown of star formation is acting in the extremely early universe; intense and compact star formation at physical scales < 1 kpc that results in a quiescent core that may later grow from the inside out via minor mergers (e.g., Bezanson et al. 2009; Belli et al. 2023).  \nThe early formation of a dense stellar core suggests that the centers of galaxies like UNCOVER 18407 could evolve directly into local massive elliptical galaxies, which have central stellar ages that are consistent with such early formation (Trager et al. 2000; Thomas et al. 2005; S'anchez-Bl'azquez et al. 2006a; Graves et al. 2007; Greene et al. 2013), although the system would need to grow by ∼ an order of magnitude in total mass over the next 10 Gyr of cosmic time to resemble massive local elliptical galaxies. This kind of mass growth appears to be quite plausible based on cosmological simulations such as IllustrisTNG, where many of the earliest-forming quiescent systems at z ∼ 4 eventually merge with an even more massive companion and evolve into brightest cluster galaxies (Hartley et al. 2023).  \nSuch early forming compact systems are suggested in cosmological simulations to be the products of secular, early dense star formation that results in a dense core with slightly older stellar populations in the center than the outskirts (e.g., Wellons et al. 2015). While the uncertainty in mass weighted age is too large to conclu-  \nsively constrain radial age gradients, the observed SEDs are not inconsistent with stellar ages that are on the order of a few 100 Myr younger in the outskirts than in the dense, central region, suggesting that the quenching of star formation could be propagating outward from the center of the galaxy as suggested in Tacchella et al. (2015).  \nThe combination of the compact mass distribution and significant dust extinction in the central regions of UNCOVER 18407 is suggestive of the galaxy's progenitor as a compact, dusty star forming galaxy in the early universe. The early formation time and high mass of UNCOVER 18407 suggests that, if it were formed in a vacuum, it must have had sustained star formation at rates comparable to the most luminous z ∼ 10 galaxies (e.g., Tacchella et al. 2023). However, due to the hierarchical assembly of galaxies, it is not necessarily true (or likely) that UNCOVER 18407 was a single galaxy at z ∼ 10 when it was rapidly forming its stars (e.g., Li et al. 2023). That said, the star formation history and structural measurements of UNCOVER 18407 are consistent with it being the descendent of individual z > 7 massive galaxy candidates that have been identified in early photometric samples; those systems exhibit extremely compact sizes at log( M ⋆ /M ⊙ ) ∼ 10 (Akins et al. 2023; Baggen et al. 2023, though their masses are quite uncertain, see Wang et al. 2024b), and the formation of ∼ 10 10 M ⊙ within a few hundred parsecs in the first Gyr of cosmic time is fully consistent with our measurements of UNCOVER 18407.  \nAdditionally, growing samples of massive, IRluminous dusty star forming galaxies that may represent a key link between rapid star formation in the early universe and quiescence at z ∼ 4 have been spectroscopically confirmed at z > 5 (Williams et al. 2023; Nelson et al. 2023b; McKinney et al. 2023; Smail et al. 2023). The detailed study of the spatial scale of that star formation and the number densities of galaxies in the first Gyr of cosmic time and beyond (e.g., Long et al. 2023) will allow for such populations to be linked to their dusty descendants such as UNCOVER 18407, as has been done with compact star forming galaxies near cosmic noon (e.g., van Dokkum et al. 2015; Barro et al. 2017; Tadaki et al. 2020).", '6. CONCLUSIONS': "In this work, we present the spectrum and resolved photometry of UNCOVER 18407, a massive, lensed, extremely red galaxy at z spec = 3 . 97. We find the following:  \n- · The rest-frame 0.9 µ m light distribution of UNCOVER 18407 is best modeled by a two-  \nponent S'ersic fit, with a compact 200 parsec, n ∼ 2 . 9 core and a more extended 1.6 kpc, n ∼ 1 . 7 component. Annular photometry indicates that the galaxy exhibits a strong F150W-F444W color gradient that drops from 4.5 in the center to 3.0 in the outskirts (see Figure 3).  \n- · The weak (equivalent width < 21 ˚ A) observed Hα in UNCOVER 18407 confirms the quiescent UVJ solutions. (see Figure 4). While deblending Hα and [NII] is difficult at the prism resolution, the 3 σ Hα limit rules out star formation rates of 1 (13) M ⊙ yr -1 under assumptions of A v =0 (3) after flux calibrating the spectrum to the photometry of the central 1.2 kpc. This upper limit is corroborated by a FIR detection (see Fujimoto et al. 2023b), which constrains the star formation rate to be less than ∼ 8 . 5 M ⊙ yr -1 at the 3 σ level.\n- · The observed central F150W-F444W and UVJ colors of the galaxy can only be modeled by an old ( ∼ 1 Gyr) stellar population with significant attenuation (see Figure 6); the universe is simply not old enough at z spec = 3 . 97 for a maximally old dust-free solution to result in the observed colors, and star forming solutions violate the constraints on the star formation rate as measured from the weak Hα emission. This implies that the core of UNCOVER 18407 is among the oldest stellar populations in the universe observed at such high redshift.\n- · Annular SED fits reveal that the color gradient in UNCOVER 18407 predominantly results from an A v gradient, with the central regions exhibiting a magnitude more extinction than the outskirts. Additionally, there are also hints of a negative age gradient, and we cannot entirely disentangle these trends due to age-dust degeneracy. The color gradients manifest in a mass-to-light gradient, suggesting that there is ∼ 1 . 5 × as much mass-perunit-light in the central ∼ kpc as there is in the less-obscured outer regions (see Figure 7).\n- · The integrated star formation history of UNCOVER 18407 is consistent with formation of ∼ 10 10 . 4 M ⊙ within the first Gyr of cosmic time. This object, in conjunction with the emerging JWST quiescent literature (e.g., Carnall et al. 2023b; Nanayakkara et al. 2024; Glazebrook et al. 2024, see Figure 8), suggests that old and dusty quiescent systems may be more common at highredshift than previously anticipated.\n- · The mass-weighted size and surface density of UNCOVER 18407 are comparable to quiescent galaxies at cosmic noon, suggesting that it represents an early example of quenching after a period of compact star formation (see Figure 9). This also suggests that the growth in size of quiescent galaxies at fixed stellar mass as a function of redshift (e.g., van der Wel et al. 2014) is dominated by growth at z < 2.  \nThe discovery of this dense, old, and dusty massive quiescent galaxy at z spec = 3 . 97 opens questions as to how common such systems are. In order to connect UNCOVER 18407 to its progenitor population, we will need to perform searches over large volumes with well understood selection functions. While the precise redshift and the constraining power of the weak Hα was useful for understanding the stellar populations in UNCOVER 18407, it is unlikely we will be able to build up well understood samples where all galaxies have spectroscopy due to the time intensive nature of performing pre-imaging followed up by targeted spectroscopy.  \nHowever, the fact that this source was identified as a potential dusty candidate in broadband imaging where the color gradient was already apparent is extremely promising for future searches. Future work can leverage the ability to select candidate quiescent galaxies in large-area imaging data to constrain how common strong color gradients are and what they can teach us about the physical processes that result in the formation of the first massive galaxies. In the coming years, a census of the quiescent population at z > 3 will allow for a detailed study of this crucial phase of evolution for some of the earliest-forming galaxies in the universe.", 'Facilities: JWST (NIRCam, NIRSpec), ALMA': 'Software: Astropy (Astropy Collaboration et al. 2013, 2018, 2022), Dynesty (Speagle 2020), grizli (Brammer 2023, github.com/gbrammer/grizli), Matplotlib (Hunter 2007), Flexible Stellar Population Synthesis (Conroy et al. 2009; Conroy & Gunn 2010), SEDPy Johnson (2019), Prospector (Johnson & Leja 2017; Leja et al. 2017; Johnson et al. 2021), Pysersic (Pasha & Miller 2023)', 'ACKNOWLEDGMENTS': "Support for this work was provided by The Brinson Foundation through a Brinson Prize Fellowship grant. This work is based in part on observations made with the NASA/ESA/CSA James Webb Space Telescope . The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with JWST Cycle 1 GO program #2561 and Cycle 3 GO program #4111. The JWST data presented in this article were obtained from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute. The specific observations analyzed can be accessed via DOI. Support for program JWSTGO-2561 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Associations of Universities for Research in Astronomy, Incorporated, under NASA contract NAS526555. DS acknowledges helpful conversations with Desika Narayanan in the early stages of assessing how interesting this object is.  \nTBM was supported by a CIERA Postdoctoral Fellowship. JAD acknowledges funding from the Dunlap Institute for Astronomy & Astrophysics. The Dunlap Institute is funded through an endowment estab-  \nlished by the David Dunlap family and the University of Toronto. RP and DM acknowledge funding from JWST-GO-2561. HA is supported by CNES, focused on the JWST mission. H.A. also acknowledges support from the Programme National Cosmology and Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, cofunded by CEA and CNES. PD acknowledge support from the NWO grant 016.VIDI.189.162 ('ODIN') and warmly thanks the European Commission's and University of Groningen's CO-FUND Rosalind Franklin program. VK acknowledges funding from the Dutch Research Council (NWO) through the award of the Vici Grant VI.C.212.036. The work of CCW is supported by NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. K.G. and T.N. acknowledge support from Australian Research Council Laureate Fellowship FL180100060. S.F. acknowledges the support from NASA through the NASA Hubble Fellowship grant HST-HF2-51505.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. Some of the data products presented herein were retrieved from the Dawn JWST Archive (DJA). DJA is an initiative of the Cosmic Dawn Center, which is funded by the Danish National Research Foundation under grant No. 140.", 'REFERENCES': 'Akhshik, M., Whitaker, K. E., Brammer, G., et al. 2020, arXiv e-prints, arXiv:2008.02276.  \nhttps://arxiv.org/abs/2008.02276  \nAkhshik, M., Whitaker, K. E., Leja, J., et al. 2021, ApJL, 907, L8, doi: 10.3847/2041-8213/abd416  \n-. 2023, ApJ, 943, 179, doi: 10.3847/1538-4357/aca677 Akins, H. B., Casey, C. M., Allen, N., et al. 2023, ApJ, 956, 61, doi: 10.3847/1538-4357/acef21 Alberts, S., Williams, C. C., Helton, J. M., et al. 2023, arXiv e-prints, arXiv:2312.12207. https://arxiv.org/abs/2312.12207 Antwi-Danso, J., Papovich, C., Esdaile, J., et al. 2023a, arXiv e-prints, arXiv:2307.09590, doi: 10.48550/arXiv.2307.09590 Antwi-Danso, J., Papovich, C., Leja, J., et al. 2023b, ApJ, 943, 166, doi: 10.3847/1538-4357/aca294 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33,  \ndoi: 10.1051/0004-6361/201322068  \nAstropy Collaboration, Price-Whelan, A. M., Sip"ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74 Baggen, J. F. W., van Dokkum, P., Labb\'e, I., et al. 2023, ApJL, 955, L12, doi: 10.3847/2041-8213/acf5ef  \nBaker, W. M., Tacchella, S., Johnson, B. D., et al. 2023,  \narXiv e-prints, arXiv:2306.02472, doi: 10.48550/arXiv.2306.02472 Barro, G., Faber, S. M., Koo, D. C., et al. 2017, ApJ, 840, 47, doi: 10.3847/1538-4357/aa6b05 Barrufet, L., Oesch, P. A., Weibel, A., et al. 2023, MNRAS, 522, 449, doi: 10.1093/mnras/stad947 Barrufet, L., Oesch, P., Marques-Chaves, R., et al. 2024, arXiv e-prints, arXiv:2404.08052, doi: 10.48550/arXiv.2404.08052 Belli, S., Newman, A. B., & Ellis, R. S. 2019, ApJ, 874, 17, doi: 10.3847/1538-4357/ab07af Belli, S., Contursi, A., Genzel, R., et al. 2021, ApJL, 909,  \nL11, doi: 10.3847/2041-8213/abe6a6  \nMosleh, M., Tacchella, S., Renzini, A., et al. 2017, ApJ, 837, 2, doi: 10.3847/1538-4357/aa5f14  \nMowla, L. A., van Dokkum, P., Brammer, G. B., et al. 2019, ApJ, 880, 57, doi: 10.3847/1538-4357/ab290a Mu˜noz Arancibia, A. M., Gonz\'alez-L\'opez, J., Ibar, E., et al. 2023, A&A, 675, A85, doi: 10.1051/0004-6361/202243528 Nanayakkara, T., Glazebrook, K., Jacobs, C., et al. 2024, Nature Scientific Reports in press, arXiv:2212.11638, doi: 10.48550/arXiv.2212.11638  \n- Nelson, E. J., Suess, K. A., Bezanson, R., et al. 2023a, ApJL, 948, L18, doi: 10.3847/2041-8213/acc1e1  \nNelson, E. J., Brammer, G., Gimenez-Arteaga, C., et al.  \n2023b, arXiv e-prints, arXiv:2310.06887,  \ndoi: 10.48550/arXiv.2310.06887  \n- Otter, J. A., Rowlands, K., Alatalo, K., et al. 2022, arXiv e-prints, arXiv:2210.12199.  \nhttps://arxiv.org/abs/2210.12199  \nPacifici, C., Kassin, S. A., Weiner, B. J., et al. 2016, ApJ, 832, 79, doi: 10.3847/0004-637X/832/1/79  \n- Park, M., Belli, S., Conroy, C., et al. 2022, arXiv e-prints, arXiv:2210.03747. https://arxiv.org/abs/2210.03747  \n-. 2024, arXiv e-prints, arXiv:2404.17945,  \ndoi: 10.48550/arXiv.2404.17945  \n- Pasha, I., & Miller, T. B. 2023, The Journal of Open Source Software, 8, 5703, doi: 10.21105/joss.05703  \nPattarakijwanich, P., Strauss, M. A., Ho, S., & Ross, N. P.  \n2016, ApJ, 833, 19, doi: 10.3847/0004-637X/833/1/19  \nPawlik, M. M., Wild, V., Walcher, C. J., et al. 2016,  \nMNRAS, 456, 3032, doi: 10.1093/mnras/stv2878  \nP\'erez-Gonz\'alez, P. G., Barro, G., Annunziatella, M., et al. 2023, ApJL, 946, L16, doi: 10.3847/2041-8213/acb3a5  \nPhan, D., Pradhan, N., & Jankowiak, M. 2019, arXiv preprint arXiv:1912.11554  \nPrice, S. H., Suess, K. A., Williams, C. C., et al. 2023, arXiv e-prints, arXiv:2310.02500,  \ndoi: 10.48550/arXiv.2310.02500  \nRowlands, K., Wild, V., Bourne, N., et al. 2018, MNRAS, 473, 1168, doi: 10.1093/mnras/stx1903  \n- S\'anchez-Bl\'azquez, P., Gorgas, J., Cardiel, N., & Gonz\'alez, J. J. 2006a, A&A, 457, 787,  \ndoi: 10.1051/0004-6361:20064842  \n- S\'anchez-Bl\'azquez, P., Peletier, R. F., Jim\'enez-Vicente, J., et al. 2006b, MNRAS, 371, 703,  \ndoi: 10.1111/j.1365-2966.2006.10699.x  \n- Sanders, R. L., Shapley, A. E., Kriek, M., et al. 2016, ApJ, 816, 23, doi: 10.3847/0004-637X/816/1/23\n- Santini, P., Merlin, E., Fontana, A., et al. 2019, MNRAS, 486, 560, doi: 10.1093/mnras/stz801\n- Saracco, P., Marchesini, D., La Barbera, F., et al. 2020, ApJ, 905, 40, doi: 10.3847/1538-4357/abc7c4'}