The recent UNCOVER survey with the James Webb Space Telescope JWST exploits the nearby cluster A2744 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 50000 robustly detected sources out to z  15. We show the redshift evolution of stellar ages star formation rates and restframe colors across the full range of 0.2  z  15. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector priors on the masses and star formation histories to produce joint redshift and stellar populations posteriors. Additionally lensing magnification is performed on the fly to ensure consistency with the scaledependent priors. We show that this approach produces excellent photometric redshifts with  SUBNMADSUB  0.03 of a similar quality to the established photometric redshift code EAzY. In line with the opensource scientific objective of this Treasury survey we publicly release the stellar population catalog with this paper derived from our photometric catalog adapting aperture sizes based on source profiles. This release the catalog and all related documentation are accessible via the UNCOVER survey web page httpsjwstuncover.github.ioDR2.htmlSPSCatalogs with a copy deposited to Zenodo at doi10.5281zenodo.8401181 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.
['Abell clusters', 'Catalogs', 'Galaxy evolution', 'James Webb Space Telescope', 'Hubble Space Telescope', 'Spectral energy distribution', '9', '205', '594', '2291', '761', '2129', 'Astrophysics - Astrophysics of Galaxies']
The UNCOVER Survey A Firstlook HSTJWST Catalog of Galaxy Redshifts and Stellar Population Properties Spanning 0.2  z  15
{'The UNCOVER Survey: A First-look HST+JWST Catalog of Galaxy Redshifts and Stellar Population Properties Spanning 0 . 2 ≲ z ≲ 15': "Bingjie Wang ( 王 冰 洁 ), 1, 2, 3 Joel Leja, 1, 2, 3 Ivo Labb'e, 4 Rachel Bezanson, 5 Katherine E. Whitaker, 6, 7 Gabriel Brammer, 7 Lukas J. Furtak, 8 John R. Weaver, 6 Sedona H. Price, 5 Adi Zitrin, 8 Hakim Atek, 9 Dan Coe, 10, 11, 12 Sam E. Cutler, 6 Pratika Dayal, 13 Pieter van Dokkum, 14 Robert Feldmann, 15 Danilo Marchesini, 16 Marijn Franx, 17 Natascha Forster Schreiber, 18 Seiji Fujimoto, 19, ∗ Marla Geha, 14 Karl Glazebrook, 4 Anna de Graaff, 20 Jenny E. Greene, 21 St'ephanie Juneau, 22 Susan Kassin, 10 Mariska Kriek, 17 Gourav Khullar, 5 Michael Maseda, 23 Lamiya A. Mowla, 24 Adam Muzzin, 25 Themiya Nanayakkara, 4 Erica J. Nelson, 26 Pascal A. Oesch, 27, 7 Camilla Pacifici, 10 Richard Pan, 16 Casey Papovich, 28, 29 David J. Setton, 5 Alice E. Shapley, 30 Renske Smit, 31 Mauro Stefanon, 32, 33 Katherine A. Suess, 34, ∗ Edward N. Taylor, 4 and Christina C. Williams 35, 36  \n1 Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 2 Institute for Computational & Data Sciences, The Pennsylvania State University, University Park, PA 16802, USA 3 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 4 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Melbourne, VIC 3122, Australia 5 Department of Physics & Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA 6 Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 7 Cosmic Dawn Center (DAWN), Niels Bohr Institute, University of Copenhagen, Jagtvej 128, København N, DK-2200, Denmark 8 Department of Physics, Ben-Gurion University of the Negev, P.O. Box 653, Beer-Sheva 8410501, Israel 9 Institut d'Astrophysique de Paris, CNRS, Sorbonne Universit'e, 98bis Boulevard Arago, 75014, Paris, France 10 Space Telescope Science Institute, Baltimore, MD 21218, USA 11 Association of Universities for Research in Astronomy, Inc. for the European Space Agency 12 Center for Astrophysical Sciences, Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 13 Kapteyn Astronomical Institute, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands 14 Department of Astronomy, Yale University, New Haven, CT 06511, USA 15 Institute for Computational Science, University of Zurich, Winterhurerstrasse 190, CH-8006 Zurich, Switzerland 16 Department of Physics & Astronomy, Tufts University, Medford, MA 02155, USA 17 Leiden Observatory, Leiden University, P.O.Box 9513, NL-2300 AA Leiden, The Netherlands 18 Max-Planck-Institut fur extraterrestrische Physik, Gießenbachstraße 1, 85748 Garching, Germany 19 Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA 20 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117, Heidelberg, Germany 21 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA 22 NSF's National Optical-Infrared Astronomy Research Laboratory, 950 N. Cherry Avenue, Tucson, AZ 85719, USA 23 Department of Astronomy, University of Wisconsin-Madison, Madison, WI 53706, USA 24 Dunlap Institute for Astronomy & Astrophysics, Toronto, Ontario, M5S 3H4, Canada 25 Department of Physics & Astronomy, York University, Toronto, Ontario, ON MJ3 1P3, Canada 26 Department for Astrophysical & Planetary Science, University of Colorado, Boulder, CO 80309, USA 27 Department of Astronomy, University of Geneva, Chemin Pegasi 51, 1290 Versoix, Switzerland 28 Department of Physics & Astronomy, Texas A&M University, College Station, TX, 77843-4242 USA  \n29 George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, Texas A&M University, College Station, TX, 77843-4242 USA  \n- 30  \nDepartment of Physics & Astronomy, University of California: Los Angeles, Los Angeles, CA 90095, USA 31 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK 32 Departament d'Astronomia i Astrofisica, Universitat de Valencia, C. Dr. Moliner 50, E-46100 Burjassot, Valencia, Spain 33 Unidad Asociada CSIC 'Grupo de Astrofisica Extragalactica y Cosmologi' (Instituto de Fisica de Cantabria - Universitat de Valencia) 34 Kavli Institute for Particle Astrophysics and Cosmology and Department of Physics, Stanford University, Stanford, CA 94305, USA 35 NSF's National Optical-Infrared Astronomy Research Laboratory, Tucson, AZ 85719, USA 36 Steward Observatory, University of Arizona, Tucson, AZ 85721, USA  \nCorresponding author: Bingjie Wang  \nbwang@psu.edu", 'ABSTRACT': '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 ∼ 15. We show the redshift evolution of stellar ages, star formation rates, and rest-frame colors across the full range of 0 . 2 ≲ z ≲ 15. The galaxy properties are inferred using the Prospector Bayesian inference framework using informative Prospector -β 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 σ NMAD ∼ 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 a) includes posterior moments, maximumlikelihood 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.  \nKeywords: Abell clusters (9) - Catalogs (205) - Galaxy evolution (594) - Hubble Space Telescope (761) - James Webb Space Telescope (2291) - Spectral energy distribution (2129)', '1. INTRODUCTION': 'The Hubble Space Telescope (HST) has secured a lasting legacy in mapping the key epochs of galaxy assembly via Treasury programs (e.g., Grogin et al. 2011; Brammer et al. 2012). JWST, with its sensitive NIRCam imaging and NIRSpec spectroscopy at 1-5 µ m (Boker et al. 2023; Rieke et al. 2023; Rigby et al. 2023), promises to not only extend this legacy, but also reveal new mysteries of the distant universe out to z ∼ 15 and beyond. Thus far, data from the JWST Early Release Science programs have already shed new light on the early phases of galactic evolution (e.g., Treu et al. 2022; Finkelstein et al. 2023). The JWST extragalactic Treasury survey, Ultradeep NIRSpec and NIRCam ObserVations before the Epoch of Reionization (UNCOVER; Bezanson et al. 2022), has completed its primary NIRCam imaging observation on and around the gravitational lensing cluster Abell 2744 at z = 0 . 308 in November 2022. The images reach depths of ∼ 30 AB magnitudes. After accounting for gravitational lensing, the intrinsic depths reach ≳ 31 -32 AB magnitudes, making UNCOVER the deepest survey in Cycle 1 of JWST observations. The quest for a coherent understanding of  \nthe universe in the newly observed parameter space is just beginning.  \nFor the first time, it is possible to observe galaxies and candidates spanning 0 . 2 ≲ z ≲ 15 and infer their stellar population properties. In the context of UNCOVER, spectroscopically confirmed galaxies at z > 12 and an X-ray luminous supermassive black holes at z = 10 . 1 are studied in detail in Wang et al. (2023b) and Goulding et al. (2023), respectively, whereas a systematic search for z spec ≳ 9 sources is carried out in Fujimoto et al. (2023a). The core of this paper is to present a galaxy catalog containing key stellar population metrics over the full dynamic range probed by the survey. This work is part of the second Data Release (DR2) from UNCOVER. In accordance with the open-data intent of the Treasury survey, we have made publicly available the imaging mosaics (Bezanson et al. 2022), updated strong lensing model (Furtak et al. 2023), and first-look photometric catalogs (Weaver et al. 2024). This paper constitutes the final installation of the UNCOVER DR2: inferred galaxy properties, including redshifts, stellar masses, metallicities, ages, star formation rates (SFRs), dust attenuation values, and fractional mid-infrared active galactic nuclei (AGN) luminosities. It is accompanied by the magnification factor, and additionally the radial magnification, tangential magnification, and shear, consistent with the inferred redshifts. The full posterior distributions are also released.  \nThe main data products in this work are inferred using the Prospector Bayesian framework (Johnson et al. 2021), with two notable modifications. First, we op-  \ntimize our priors for recovering accurate photometric redshifts by including observationally-motivated, joint priors on stellar mass, stellar metallicity, and star formation history (SFH) from Prospector -β (Wang et al. 2023c). Second, we solve the magnification-redshift relationship on-the-fly within Prospector to take advantage of the mass-dependent priors, in contrast to the traditional approach where physical parameters are scaled by magnification factors post-fit.  \nInferring redshifts across a wide range of distances and galaxy properties based solely on photometric data pushes the stellar population models into new, exciting, and largely uncalibrated regimes. Comparisons between different photoz codes have clearly highlighted that accuracy in this space depends strongly on the assumptions that go into the code (e.g., Kodra et al. 2023). Given the novelty of our Prospector model, we cross-check our redshifts with those from the established template-fitting code EAzY (Brammer et al. 2008). We additionally analyze the effect of varying template sets within EAzY by comparing results from the default set fsps to sfhz . The latter incorporates redshift-dependent SFHs and a realistic emission line model at z ∼ 8 (Carnall et al. 2023).  \nThe structure of this paper is as follows. Section 2 summarizes the photometric data. Section 3 details the SED modeling. Section 4 presents the inferred parameters. Section 5 discusses known problems in the photometry and in the modeling, and how these may affect the accuracy of the inferred stellar population parameters. Section 6 concludes with a brief summary and the format of our catalogs.  \nWhere applicable, we adopt the best-fit cosmological parameters from the 9 yr results from the Wilkinson Microwave Anisotropy Probe mission: H 0 = 69 . 32 kms -1 Mpc -1 , Ω M = 0 . 2865, and Ω Λ = 0 . 7135 (Hinshaw et al. 2013), and a Chabrier initial mass function (IMF; Chabrier 2003). Unless otherwise mentioned, we report the median of the posterior, and 1 σ error bars are the 16th and 84th percentiles.', '2. DATA': "The photometry includes all public JWST/NIRCam, HST/ACS and HST/WFC3 imaging of Abell 2744 available to date. Specifically, the JWST data includes the Cycle 1 Treasury program UNCOVER covering ∼ 45 arcmin 2 (PIs Labb'e & Bezanson, JWST-GO-2561; Bezanson et al. 2022), the Early Release Science program GLASS (PI: Treu, JWST-ERS-1324; Treu et al. 2022), and a Director's Discretionary program (JWSTDD-2756, PI: Chen). These observations span ∼ 1 -5 µ m in observed frame in 8 filters: F090W, F115W, F150W,  \nF200W, F277W, F356W, F410M, and F444W. The HST data, taken from the public archive, consist of HST-GO11689 (PI: Dupke), HST-GO-13386 (PI: Rodney), HSTDD-13495 (PI: Lotz; Lotz et al. 2017), and HST-GO15117 (PI: Steinhardt; Steinhardt et al. 2020). These additional observations span ∼ 0 . 4 -1 . 6 µ m in observed frame in 7 filters: F435W, F606W, F814W, F105W, F125W, F140W, and F160W. Details of these programs including imaging depths, and ancillary information are summarized in Tables 1 and 3 in Bezanson et al. (2022). A small subsample has detections or upper limits from DUALZ (Deep UNCOVER-ALMA Legacy High-Z Survey; Fujimoto et al. 2023b). The transmission curve of this ALMA band is approximated as a top-hat function spanning 1249.9 to 1351.0 µ m.  \nAs part of UNCOVER DR2, we have released F277W+F356W+F444W-selected photometric catalogs containing total fluxes for 61,648 sources (Weaver et al. 2024). In this work we fit galaxy SEDs to the 'supercatalog,' which uses adaptive aperture selection based on their isophotal areas following Labb'e et al. (2003). We share publicly the inferred parameters for the full sample. However, it is worth noting that a subsample of 55,613 objects is deemed to have reliable photometry, i.e., use phot =1 in the photometric catalog. Furthermore, 15,861 objects in this subsample are flagged as possible blends ( flag kron = 1), and therefore their photometry is not corrected to total based on a kron ellipse but instead is corrected assuming a point-like morphology. Photometry performed on simulated galaxies with realistic sizes indicate that the total photometry of these objects may be underestimated by a factor of 2 in 0.32-0.70 '' apertures, with larger 1.00-1.40 '' apertures miss < 10% of the total flux. Consequently, physical parameters derived in this work including stellar mass and rest-frame fluxes are liable to be underestimated. Details on the selection criteria for the use phot and flag kron flags are presented in Section 4.6 in Weaver et al. (2024).  \nAtotal of ∼ 400 reliable, unblended sources in the UNCOVER photometric catalogs have spectroscopic redshifts collected from the NASA/IPAC Extragalactic Database, and from the literature. The latter consists of measurements taken with HST as part of the Grism Lens Amplified Survey from Space (GLASS; Treu et al. 2015), and with the Multi Unit Spectroscopic Explorer (MUSE; Richard et al. 2021). We utilize this sample to evaluate the accuracy of our photometric redshift recovery.", '3. SED FITTING': 'Figure 1. An overview of the key inferred parameters in our catalog. (upper panel) Joint distribution of redshifts and stellar masses. The x-axis is in linear scale at z ≤ 10, while in logarithmic scale at z > 10. The redshift at which the Abell 2744 cluster resides ( z = 0 . 308) is indicated by an orange dashed line. The secondary peak in the mass distribution at ∼ 10 6 M ⊙ is mainly due to the large population of globular clusters. (lower panel) Mass-weighted ages, star formation rates (100 Myr), and specific star formation rates (100 Myr) as functions of stellar masses in redshift bins.  \n<!-- image -->  \n<!-- image -->  \n/circledot  \n<!-- image -->  \n/circledot  \n<!-- image -->  \n/circledot  \nThis section describes the core of this work-inferring physical parameters of the galaxy populations from the photometry. An overview of the inferred parameters of our catalog is shown in Figure 1.  \nIn what follows, we first describe our modified Prospector model, which includes physically motivated priors to optimally produce joint constraints of redshifts, stellar masses, and other key stellar population metrics across any redshift range, and a consistent treatment of lensing magnification during model fitting. We then outline the adopted EAzY settings, and its two template sets, fsps and sfhz . A minimum error floor of 5% is imposed in all fits to reflect the additional systematic uncertainties-the calibration uncertainties of JWST/NIRCam at this early stage (Boyer et al. 2022), as well as uncertainties in the stellar population synthesis (SPS).', '3.1. Prospector': "The main data products in the catalog are inferred using the Prospector Bayesian inference framework (Johnson et al. 2021). The building blocks of these fits, i.e., the simple stellar populations (SSPs), come from FSPS (Conroy & Gunn 2010), where we adopt MIST isochrones (Choi et al. 2016; Dotter 2016) and MILES stellar library (S'anchez-Bl'azquez et al. 2006). The composite stellar populations (CSPs) are modeled with Prospector -β (Wang et al. 2023c), which follows Prospector -α (Leja et al. 2017) in many components. In this section, we begin with the common elements in Prospector -α , and then proceed to the new additions from Prospector -β , namely the joint priors on stellar mass, stellar metallicity, and SFH. We end with a novel modification in the likelihood calculation which allows for a self-consistent treatment for lensing magnification during model fitting.  \nThe SFH is modeled as mass formed in 7 logarithmically-spaced time bins, and assumes a continuity prior to ensure smooth transitions between bins (Leja et al. 2019a). The scheme for the age bins is refined in this work. For z < 3, we keep the conventional definition in which the first two bins are always 30 Myr and 100 Myr respectively, whereas for z ≥ 3, we only require that the first bin is always 13.47 Myr. In both cases, the last bin is 10% of the age of the universe at a given redshift, and the intervening bins are evenly spaced in logarithmic time. This new scheme ensures that no age bins are overly wide in the early universe. Nebular emission is included using a pre-computed Cloudy grid (Byler et al. 2017). Dust is described using a two-component model  \n(Charlot & Fall 2000) with a flexible dust attenuation curve (Noll et al. 2009). We also fit for the stellar metallicity, and the normalization and dust optical depth of mid-infrared AGNs. Dust emission is included in all fits (Draine & Li 2007), with the mass fraction of polycyclic aromatic hydrocarbons left free. The attenuation of the intergalactic medium (IGM) is assumed to follow Madau (1995).  \nIn contrast to previous large-scale applications of Prospector to broadband photometry, which typically fix redshift to a value determined by an external photoz code (e.g., Leja et al. 2019c), here we fit directly for redshift and use informed priors for the stellar mass, stellar metallicity, and SFH in the Prospector -β model. We opt for this approach to take the full advantage of the Prospector Bayesian inference framework, and thus to obtain consistent joint constraints on the probability distribution of the full posterior. Following Wang et al. (2023c), we first include a mass prior P(logM ⋆ | z ), constructed from the observed mass functions between 0 . 2 < z < 3 (Leja et al. 2020). For z < 0 . 2 and z > 3, we take the nearest-neighbor solution, i.e., the z = 0.2 and z = 3 mass functions. This applies the mass function as a prior where it is well-measured, and avoids relying on simulation predictions while making a reasonable null hypothesis in the absence of observational constraints. Second, we use a dynamic SFH prior, meaning that the shape of the SFH is dependent on redshift and stellar mass. The expectation value of this prior is matched to the cosmic SFR density (Behroozi et al. 2019); in other words, it encourages rising histories early in the universe, and falling histories late in the universe. This prior additionally reflects the consistent observational finding that massive galaxies form much earlier than low-mass galaxies (Cowie et al. 1996; Thomas et al. 2005) by introducing a hyper-parameter that scales the expectation value for the shape of SFH with mass, such that massive galaxies have more falling SFHs and lowmass galaxies have more rising SFHs. As noted in Wang et al. (2023c), even though the high-redshift constraints on galaxy evolution remain uncertain, these priors still represent a better expectation then the null expectation of uniform priors. Third, we place a prior based on the stellar mass-stellar metallicity relationship measured from the SDSS (Gallazzi et al. 2005). Following Leja et al. (2019a), we take the conservative approach of widening the confidence intervals from this relationship by a factor of 2 to account for potential unknown systematics or redshift evolution. This prior neglects predictions from galaxy formation models, which suggest a smaller scatter and a relatively strong redshift evolution (Ma et al. 2016; Feldmann et al. 2023).  \nTable 1. Prospector -β Parameters and Priors Adopted in This Paper  \nThe complete list of the free parameters and their associated priors is summarized in Table 1. Only the first 8 parameters in the table are used to infer physical parameters that are provided in the first public release. The other parameters are used as nuisance parameters in the great majority of cases, and will be examined further in the future.  \nIt is worth noting that we correct for lensing magnification simultaneously while fitting the full set of parameters. Gravitational lensing is in general achromatic (i.e., the deflection angle of a light ray is independent of its wavelength), meaning that colors are conserved. Therefore, while the magnification factor, µ , depends on redshift and source position, it is not conventionally accounted for in the process of SED fitting. Rather, the modeled fluxes and the scale-dependent physical parameters such as stellar mass are often divided by µ post-fit, or the observed fluxes are demagnified before fitting using external redshift information. However, scaledependent priors in Prospector -β necessitate a self-consistent treatment of µ . We devise a simple method: given a source position, we read in the convergence, κ , and shear, γ , from the 0.1 '' /pixel resolution maps provided as part of UNCOVER DR (Furtak  \net al. 2023) 1 . These maps are normalized such that D ds /D s = 1, where D ds is the angular diameter distance between the lens and the source planes, and D s is the distance between the source plane and the observer. Then in each draw of the redshift during sampling, we multiply κ and γ by the D ds /D s of the source, and magnify the model photometry by µ , which is calculated as  \nµ = 1 | (1 -κ ) 2 -γ 2 | . (1)  \nCritically, magnifying the model fluxes and demagnifying the observations are not equivalent due to the definition of likelihood in Prospector :  \nln L = N ∑ n =1 ln [ 1 √ 2 πσ 2 · exp ( xn -α ) 2 -2 σ 2 ] , (2)  \nwhere σ is the observed uncertainties, x n is the observed flux in the n th photometric band, and α is the model uncertainties. The often neglected pre-factor is properly included here because Prospector offers the functionality of estimating the observational uncertainties simultaneously with the object parameters. In the context of the magnification calculation, if we change the observational uncertainties with redshift, the pre-factor  \nFigure 2. An example from the UNCOVER field. The upper left shows the SED of the galaxy. The observed photometry is plotted as black dots, while the maximum-likelihood photometry and spectrum are shown in red. The lower left shows the inferred SFH of the galaxy. The x-axis is the lookback time in Gyr, ans the y-axis is the SFR in M ⊙ yr -1 . The RGB color composite is made from F444W, F277W, and F115W bands. The lower right is a corner plot showing the marginal and joint posterior distributions of redshift, stellar mass, stellar metallicity, star formation rate, mass-weighted age, and optical depth of the diffuse dust.  \n<!-- image -->  \nworks to minimize these uncertainties, or in other words to maximize µ . Therefore, we must magnify the model fluxes for the likelihood to behavior correctly.  \nThe posterior space is sampled with the nested sampler dynesty (Speagle 2020), and a neural net emulator, dubbed parrot , which mimics SPS models is used to decrease the runtime (Alsing et al. 2020; Mathews et al. 2023). Given that this work is the first large-scale application of the emulator beyond verification tests with the 3D-HST catalogs (Mathews et al. 2023), we test its accuracy compared to full FSPS evaluations in Appendix A. We also supply further diagnostics on photometric residuals in Appendix B. An example of an SED fit is shown in Figure 2.", '3.2. EAzY': 'EAzY is a galaxy photometric redshift code that fits the observed SEDs as a non-negative sum of templates by minimizing the χ 2 statistics (Brammer et al. 2008). EAzY offers flexibility in fitting different data sets through the modification of templates. We start by the standard publicly available fsps template set, which consists of 12 templates spanning a range of colors. As pointed out by, e.g., Steinhardt et al. (2023), the standard method of running EAzY can allow for nonphysical contributions from templates that are older than the age of the universe at a given redshift thereby biasing the best fits. The latest sfhz template set is designed  \nto mitigate this problem. Redshift dependence is introduced in the templates so that SFHs starting earlier than the age of the universe at a given epoch are disallowed. A template from the JWST/NIRSpec observation of a z ∼ 8 extreme emission-line galaxy is also added in this set (Carnall et al. 2023), in order to expand the template set to include the more exotic emission-line contributions that have been observed at the highest redshifts.  \nIn consistency with the EAzY settings in Weaver et al. (2024), we turn off the priors on the magnitude and the UV slope. We do not apply the usual methodology of an iterative application of photometric zero-point offsets either. The agreement between the spectroscopic and photometric redshifts marginally worsens after applying the correction, though we caution that the redshifts in this field have a strong selection function due to the foreground cluster and the relatively low number of spectroscopic redshifts. We defer a full exploration of the photoz accuracy in this field to future works which will exploit forthcoming medium-band and grism redshifts.', '4. INFERRED STELLAR POPULATION PARAMETERS': 'We now present the parameters contained in our stellar population catalogs, an overview of which is shown in Figure 1. Here we also compare the photometric redshifts inferred using Prospector -β , EAzY -sfhz , and EAzY -fsps as a quality check.', '4.1. Photometric Redshifts': "We show the comparison between our recovered photometric redshifts and available spectroscopic redshifts in Figure 3. The scatter in residuals is quantified using the normalized median absolute deviation (NMAD; Hoaglin et al. 1983) given its advantage of being less sensitive to outliers than standard indicators, e.g., RMS. It is defined as  \nσ NMAD = 1 . 48 × median | ∆z | , (3)  \nwhere ∆ z = ( z phot -z spec ) / (1 + z spec ). We additionally quantify an outlier fraction, f out , in which we define a catastrophic outlier as one with | ∆ z | > 0 . 15, and bias calculated using the mean bias error (MBE) as  \nδ MBE = 1 n n ∑ i =1 ∆ z. (4)  \nOverall, Prospector -β performs favorably. EAzY 's performance depends on the adopted template set: the fsps template set has a low scatter but a high outlier fraction, while the sfhz template has middle-of-the-road  \nperformance in both. This template-dependent performance from EAzY is liable to change further as the photometric extraction is improved, and more generally template-based fitting appears to be more sensitive to the choices made during photometric catalog construction. Importantly, as the spectroscopic sample unavoidably has a strong selection bias, we caution not to overinterpolate these results.  \nWeproceed to compare the photometric redshifts from Prospector -β , EAzY -sfhz , and EAzY -fsps in Figure 4. For the ease of comparing the statistics, we take results from Prospector -β as the reference point; that is, ∆ z = ( z EAzY -z Prospector -β ) / (1 + z Prospector -β ). EAzY -sfhz notably produces smaller scatter and bias than EAzY -fsps when comparing to the Prospector -β redshifts, although both agree reasonably well on a high signal-to-noise (S/N) subsample. The agreement deteriorates, however, when all the reliable photometric data is included. As seen in the lower panel of Figure 4, break confusion is not the main cause. Combining with the findings from the spectroscopic sample, model mismatch is a more likely explanation. The scheduled spectroscopic observation of UNCOVER will allow for a more definitive conclusion.  \nThe existence of z ≳ 9 sources in this catalog is apparent in both Figures 1 and 4. Possible contamination from lower-redshift objects and artifacts is discussed in Section 5.1. An analysis of the modeling uncertainties of a z ≳ 9 sample, selected based on redshift posterior distributions contained in the catalog of this paper, is presented in Wang et al. (2024). This is complementary to the study of a color-selected sample in the UNCOVER field, in which different photometric catalog and SED fitting methods are used (Atek et al. 2023).", '4.2. Stellar Masses and Other Ancillary Parameters': 'In addition to photometric redshifts, we also release stellar mass and other ancillary parameters as listed in Table 2. Stellar mass can often be robustly constrained by photometry; however, extra uncertainties are introduced by varying redshifts and magnifications. Lens models, especially for highly magnified regions, can be uncertain due to systematics (Acebron et al. 2017; Bouwens et al. 2017; Meneghetti et al. 2017; Acebron et al. 2018; Atek et al. 2018; Furtak et al. 2021). Therefore, the uncertainties in magnifications, which are not included in this catalog, can be significant (see Furtak et al. 2023 for uncertainties on the lens model used in this work). We expect to fully incorporate the lensing uncertainties in the next generations of catalogs.', '4.3. Rest-frame Colors': "z  \nEAzY-sfhz  \nEAzY-fsps  \nProspector-  \nβ  \nFigure 3. Accuracy of the photometric redshifts. (Left) The photometric redshifts inferred from the three different settings described in Section 3 are plotted against known spectroscopic redshifts. Also included are the summary statistics quantifying the scatter ( σ NMAD ; Equation 3), outlier fraction ( f out ( | ∆ z | > 0 . 15)), and bias ( δ MBE ; Equation 4). These statistics suggest that we can achieve high-quality redshifts for bright objects. (Right) The distributions of photometric redshift residuals normalized by the 1 σ width of the posteriors are illustrated by kernel density estimations. A unit Gaussian is overplotted as a gray dashed curve to guide the eye. The data are clipped to be within 3 σ of the unit Gaussian. This calibration plot suggests that all three codes tend to overestimate the uncertainties for typical objects, and underestimate the uncertainties for the outliers.  \n<!-- image -->  \n<!-- image -->  \n-  \nConsidering the common practice of using rest-frame colors to categorize star-forming and quiescent galaxies (e.g., Williams et al. 2009; Brammer et al. 2011), we provide rest-frame fluxes in UVJ filters, and also fluxes in the recently proposed synthetic ugi filters (Antwi-Danso et al. 2023), both of which are marginalized over the full Prospector -β posteriors. We show the UVJ and ugi color-color diagrams in redshift bins where the data has constraining power on the respective colors in Figure 5. Specifically, we only show the UVJ colors for galaxies in the following redshift range:  \n- · z phot > 1, if optical data (HST F606W/F814W) exists\n- · z phot > 2 . 3, otherwise\n- · z phot < 4  \nLikewise we only show the ugi colors for galaxies in the following redshift range:  \n- · z phot > 1 . 4, if optical data (HST F606W/F814W) exists\n- · z phot > 3 . 1, otherwise\n- · z phot < 6  \nThese colors are, however, provided for all objects in the catalog, including cases where extrapolations are necessary.  \nThe UVJ color-selection criteria for quiescent galaxies are  \n( U -V ) > 1 . 23 ∧ ( V -J ) < 1 . 67 ∧ ( U -V ) > ( V -J ) × 0 . 98 + 0 . 38 , (5)  \nwhereas the synthetic ugi color-selection criteria are  \n( u -g ) > 1 . 5 ∧ ( g -i ) < 1 . 8 ∧ ( u -g ) > ( g -i ) × 0 . 73 + 1 . 08 . (6)  \nBoth equations are taken from Antwi-Danso et al. 2023. The synthetic ugi colors are notably better correlated with sSFR at 1 . 5 ≲ z ≲ 3, which suggests that it may be a suitably complementary diagnostic for JWST observations. We provide further discussion on the color uncertainties and the utility of rest-frame colors calculated from the best-fit model versus posterior-averaged quantities in Appendix C.  \n(EAzY-sfhz)  \nz  \n(EAzY-sfhz)  \nz  \n20  \n10  \n8  \n6  \n4  \n2  \n0  \n20  \n10  \n8  \n6  \n4  \n2  \n0  \n0  \n0  \n2  \n4  \n6  \n8  \n10  \n20  \nz  \n(Prospector-  \nβ  \n)  \n0  \n2  \n4  \n6  \n8  \n10  \n20  \nz  \n(Prospector-  \nβ  \n)  \nFigure 4. Comparison between EAzY and Prospector redshifts. (Upper panel) Photometric redshifts inferred from Prospector -β are plotted against those fitted using the latest template set EAzY -sfhz and the standard template set EAzY -fsps , respectively. Here only the subset satisfying S/N ≥ 10 in the long-wavelength detection bands (F277W, F356W, and F444W) is shown. All the statistics are computed by taking the Prospector -β redshifts as 'truth' for the ease of comparison. The scales on both axes are switched from linear to logarithmic at z = 10. (Lower panel) Similar to the upper panel, but the photometric redshifts from the full set of sources with reliable photometry are included (i.e., use phot =1). Break confusions (a combination of Balmer, 4000 ˚ A, Lymanα , and Lyman-continuum) are indicated by the colored shades.  \n<!-- image -->  \nLWS  \n/  \nN  \nσ  \nf  \nδ  \nNMAD  \nout  \n(  \n|  \nMBE  \n∆  \nz  \n≥  \n10  \n=0.053  \n|  \n>  \n0.15) = 0.220  \n=0.062  \n2  \n4  \n6  \n8  \n10  \n20  \nz  \n(Prospector-  \nβ  \n)  \nuse phot = 1  \nσ  \nf  \nδ  \nNMAD  \nout  \n(  \n|  \nMBE  \n∆  \nz  \n=0.228  \n|  \n>  \n0.15) = 0.555  \n=0.505  \nMisidentified Lyman  \nMisidentified Balmer  \n20  \n10  \n8  \n6  \n4  \n2  \n0  \n(EAzY-fsps)  \nz  \nLWS  \n/  \nN  \nσ  \nf  \nδ  \nNMAD  \nout  \n(  \n|  \nMBE  \n∆  \nz  \n≥  \n10  \n=0.118  \n|  \n>  \n0.15) = 0.246  \n=-0.202  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 5. Here we show the rest-frame colors from marginalizing over the Prospector -β posteriors in redshift bins. The median and the 90% of the uncertainty distributions are shown as black and gray error bars, respectively, in the upper left corner. The asymmetry of the error bars are driven by the non-Gaussianity of the redshift posterior distribution. The galaxies are color-coded by sSFR. The synthetic ugi colors (lower panel) are notably better correlated with sSFR than the UVJ colors (upper panel) at 1 . 5 ≲ z ≲ 3. The best-fit versions of these plots are available in Appendix C.  \n<!-- image -->  \n-  \nAs a validation of the photometric modeling, we select a sample of z > 2 quiescent candidates and crossmatch to a sample of spectroscopically confirmed quiescent galaxies within Abell 2744 (Marchesini et al. 2023). We choose a redshift-dependent definition of quiescence, calculated using a mass doubling number as  \nD ( z ) = sSFR 100 ( z ) × t H ( z ) , (7)  \nwhich is the number of times the stellar mass doubles within the age of the universe at redshift z , t H ( z ), at the current sSFR (Tacchella et al. 2022). Focusing on the z > 2 population, we select our quiescent sample to meet all of the following criteria based on Prospector -β outputs:  \n- · 16th percentile of the redshift posteriors > 2 in the Prospector -β fits\n- · 50th percentile of the redshift posteriors < 9 in the Prospector -β fits\n- · D ( z ) ≤ 1 / 20  \nAll the candidates passing the above criteria have no data quality issues in the cutouts or SEDs upon visual inspection.  \nThe resulting sample is shown as stars on the UVJ and the ugi color plane in Figure 6. We additionally mark the two spectroscopically confirmed quiescent candidates at z > 2 on the same figure with squares (Marchesini et al. 2023), both of which are quiescent in our photometric fits.  \nThe implications from Figure 6 are consistent with previous results in the literature: not all quiescent galaxies fall in the rest-frame color-color criteria, and not all galaxies meeting rest-frame color-color criteria are quiescent. This challenge is exacerbated at higher redshifts where all stellar populations are by definition younger and bluer (Leja et al. 2019b; Carnall et al. 2023; Gould et al. 2023). This finding also suggests that full SEDfitting can bring additional value for the selection of quiescent galaxies.  \nA tailored UNCOVER catalog of quiescent candidates based on different selection criteria and model assumptions will be presented in Khullar et al in prep. More dedicated discussions on quiescence can also be found therein.", '5. KNOWN ISSUES IN THE CURRENT DR AND UPCOMING IMPROVEMENTS': 'Thus far, this paper has presented the first-generation galaxy catalog of UNCOVER, derived from photometry that reaches intrinsic depths of ≳ 31 -32 AB magnitudes  \nafter correcting for lensing magnification. We emphasize that the intention of this public release is to make the catalogs available for rapid science and is not intended to be a finalized catalog. Therefore, the user should be cautioned of the known caveats. In what follows, we discuss the known issues with the photometry and with the SPS models, as well as upcoming improvements.', '5.1. Known Challenges in the Photometric Catalog': "This work is intended to model the photometry as given; that is, it does not attempt to solve possible issues in the first-generation photometric catalogs. Here we discuss the known issues in the context of this work; further details can be found in Weaver et al. (2024).  \nInterpreting early JWST imaging data was made difficult by uncertain photometric calibration. Thanks to efforts across the community, photometric zero-points are now thought to be well understood across the detector to < 5% (Boyer et al. 2022; see also photometric calibration presented in the methods section of Labb'e et al. 2023).  \nHowever, it is still in the early stages of understanding and improving the known artifacts in JWST imaging (Rigby et al. 2023). These include 'claws' and 'snowballs,' as well as hot pixels that are particularly persistent in long-wavelength (LW) bands near the detector edges. Unidentified, these artifacts can masquerade as high-redshift galaxies, or cause spurious signals that leads to genuine high-redshift galaxies being misclassified. Such issues are difficult to identify from photometry alone, and so complicate spectral fitting. Significant effort has been expended by our team to remove these features at the reduction level; as such these issues are largely resolved in our most recent image reduction (Bezanson et al. 2022) and photometric catalog (Weaver et al. 2024).  \nThe final known issue worth discussion is that subtracting the brightest cluster galaxies leads to spurious detections of residual features. We conservatively flag objects within 3 '' of all subtracted cluster galaxy centers. However, there may remain a number of spurious sources at larger radii that are difficult to flag without simultaneously flagging robust sources. While the photometry and spectral fit of such sources may look reasonable, it is obvious when inspecting the image stamps that these sources are artifacts. While it is good practice to inspect the images of any exciting target, we especially encourage this for sources in the immediate vicinity of the three cluster cores.", '5.2. Known Unknowns in the Modeling': 'SPS models require many ingredients, including an IMF, isochrones, and stellar spectra for the construc-  \n<!-- image -->  \nV - J  \nFigure 6. All z > 2 galaxies that satisfy the quiescent criteria of the mass doubling number D ( z ) ≤ 1 / 20 are shown as stars on the two color-color planes. Additionally, all z > 2 galaxies that fall into the UVJ/ugi wedge but have D ( z ) > 1 / 20 are shown as dots. The solid squares indicate the two quiescent galaxies confirmed with NIRISS spectroscopy (Marchesini et al. 2023).  \n<!-- image -->  \ntion of SSPs, SFHs and stellar metallicity models for the construction of CSPs, a model for dust attenuation and emission, and nebular continuum and line emission (see Walcher et al. 2011; Conroy 2013 for recent reviews). Comprehensive assessments of SED fitting, based on data of exquisite quality and wavelength coverage at z ≲ 3, have already highlighted the dependence of inferred parameters on modeling assumptions (Pacifici et al. 2023). In this work, we infer a panoramic view of the galaxy populations out to z ∼ 15, reaching a parameter space where robust theoretical models have yet to undergo robust tests. Accordingly, we discuss a few salient points on the known unknowns in the stellar population modeling in this section.  \nThe stellar templates are particularly uncertain at high stellar masses and at low metallicities (e.g., Johnson et al. 2021). The evolution of massive stars is strongly affected by the choice of input physics, including the treatment of convection, close binary evolution, rotation effects, and mass-loss, each of which has its own uncertainties (Leitherer et al. 1999; Choi et al. 2016; Eldridge et al. 2017; Stanway & Eldridge 2018; Byrne et al. 2022). More importantly, current models assume a solar abundance pattern by necessity, but high-resolution spectra of quiescent galaxies have revealed variations in α -element abundances (Thomas et al. 2005; Choi et al. 2014; Conroy et al. 2014; Onodera et al. 2015; Kriek et al. 2016). These trends in abundance patterns can change the fluxes by 10-40% (Vazdekis et al. 2015; Choi et al. 2019). As for metallicity histories, we follow the established methodologies in the field that the same  \nmetallicity is assumed across different bins in the SFH (Mitchell et al. 2013; Webb et al. 2020; Tacchella et al. 2022). Although a possible time dependence has been claimed (Bellstedt et al. 2021; Thorne et al. 2022), we think that our approach is adequate for this work because the catalog is dominated by high-redshift, lowmass objects. These objects tends to have rising SFHs, whereas the metallicity histories are most impactful for star-forming objects with falling SFHs. While the exact effects of the aforementioned uncertainties on SED modeling remain unclear, they can be partly illustrated by the disparate solutions found when using different stellar templates, a variety of which exist in the literature. Those templates cover different ranges of stellar evolutionary tracks and isochrones. In some cases, it has been shown that the different choices (e.g., MIST vs. PARSEC in Whitler et al. 2023) can lead to the inferred stellar ages differing by up to an order of magnitude.  \nNebular emission can make up a significant fraction of the total flux for stars at low metallicity and at young ages (Anders & Fritze-v. Alvensleben 2003). This alone means that it plays an outsize importance at high redshift (e.g., Smit et al. 2014). Nebular emission also becomes increasingly important at high redshifts for a more technical reason: the redshifting of the spectrum causes a feature with a fixed rest-frame EW to occupy a larger fraction of the filter bandpass. It is a known challenge to model nebular emission in the early universe, due to a complex combination of non-solar abundance patterns, large uncertainties in the incident ionizing radiation, and high densities of both gas and ionizing pho-  \ntons (Schaerer & de Barros 2009; Stark et al. 2013; de Barros et al. 2014; Gutkin et al. 2016; Steidel et al. 2016; Byler et al. 2017; Strom et al. 2018; Freeman et al. 2019). Wang et al. (2024) estimate the systematic uncertainty driven by nebular physics in the inferred galaxy properties by employing a flexible nebular emission model, and find a ∼ 0 . 2 dex systematic increase in stellar mass.  \nThe initial distribution of masses for a population of stars (i.e., IMF), influences almost all the inferred properties of stellar populations-mostly notably, the total luminosity and total stellar mass. Various theories predict different shapes of the IMF (e.g., Salpeter 1955; Kroupa 2001; Chabrier 2003). While it is typically assumed to be constant, there is emerging evidence that it varies across time and environment (Conroy & van Dokkum 2012; La Barbera et al. 2013; Spiniello et al. 2014; Lyubenova et al. 2016; Lagattuta et al. 2017; van Dokkum et al. 2017). The recently proposed temperature dependence in the IMF (Steinhardt et al. 2023) and variations in a low-metallicity environment (Chon et al. 2021, 2022) may lead to a systematic departure from the universality particularly at high redshifts. Here we assume a fixed Chabrier (2003) IMF, and caution that the inferred star formation rates and stellar masses can vary by up to a factor of ten based on this assumption (Wang et al. 2024).  \nIt is worth noting that the aforementioned systematic issues apply to all SED fitting. The inflow of JWST data breathes new life into model fitting and interpretation, as is evident from an active discussion in the literature (Ferrara et al. 2023; Kannan et al. 2023; Topping et al. 2022; Adams et al. 2023; Mauerhofer & Dayal 2023; Mirocha & Furlanetto 2023; Reddy et al. 2023; Yung et al. 2024). The effects of varying SPS model assumptions including burstiness in SFH, non-universal IMF, and nebular physics, using data from this catalog, are examined in detail in Wang et al. (2024).', '5.3. Reliability of Photometric Redshifts and Future Improvements': 'Given the extensive discussion on inferring photometric redshifts in the literature (see Newman & Gruen 2022 for a recent review, and also Alsing et al. 2023; Leistedt et al. 2023 for general discussions on photometric redshift inference), we focus only on the unique issues faced in this work below.  \nConsiderable uncertainties exist in the lens models, including systematic uncertainties between different models and in the highly magnified regions within a model (e.g., Zitrin et al. 2015). Additionally, we note that the lensing maps do not cover our entire field of view, and thus all sources that fall outside of the lensing maps are  \nassigned µ = 1. In reality the peripheral areas can be magnified by µ ≲ 1 . 3. A complete set of lensing maps will be released in the future. The full effect of magnification uncertainties on redshifts and other parameters will also be studied, where we add µ as a free parameter with an informative prior from the lensing maps in the SED fitting.  \nMost photometric redshift information comes from the position of spectral breaks, e.g., the dropout technique (Steidel et al. 1996). In principle, all high-redshift objects, with specific wavelength coverage and high S/N images taken across the multiple bands, are detectable via this technique. An obvious concern is the possible contamination from lowz objects exhibiting breaks in similar locations. In particular, it can be challenging to distinguish between a Balmer break and strong emission lines (Dunlop et al. 2007; Naidu et al. 2022; Arrabal Haro et al. 2023; McKinney et al. 2023; Zavala et al. 2023). This degeneracy is mitigated, but not removed, by the adoption of the mass function prior in Prospector -β .  \nFurther uncertainties in the inferred photometric redshifts come from the modeling of Lymanα and its interaction with the IGM. The emulator used in this work (Mathews et al. 2023) is trained on a Cloudy grid (Byler et al. 2017), and hence does not accurately account for the radiative transfer process (e.g., damping wings). We defer a careful modeling of the Lymanα to future works.  \nObservationally, increasing the wavelength coverage or resolution will improve the accuracy as well as the precision of the inferred redshifts. The Abell 2744 field will be observed with all the JWST medium bands during Cycle 2 (PI Suess, JWST-GO-4111). These observations are expected to be especially helpful in determining the uncertain emission line contribution. On the modeling side, further improvements will come from an accurate estimation of the detection efficiency as a function of mass based on the flux limits in the photometry. This will inform our model about the survey volume and downweight spurious high-redshift solutions.  \nTaken together, we have put our best effort forward to create the first-look catalogs. For the photometric catalogs, Weaver et al. (2024) calibrate the photometry and remove known spurious sources as cleanly as possible; whereas for the stellar population catalogs of this paper, we adopt a full Bayesian approach, incorporating empirical priors and correcting for magnification during model fitting, to ensure the maximum science return from the early observations. However, we again caution the reader to have a healthy level of skepticism in working with first-generation JWST data, especially when using sources in more exotic parameter spaces.', '6. SUMMARY': "In this paper, we present the first-generation galaxy catalogs spanning 0 . 2 ≲ z ≲ 15, as part of the public release from the JWST extragalactic Treasury survey, UNCOVER. We adopt the Prospector Bayesian inference framework (Johnson et al. 2021), within which we use the non-parametric SFH model in Prospector -α (Leja et al. 2019a), and three observationally motivated priors on the stellar mass functions, stellar metallicities, and SFHs in Prospector -β (Wang et al. 2023c). We constrain redshifts and galaxy properties simultaneously, meaning that the commonly non-Gaussian redshift uncertainties are propagated into inferred properties of the galaxy population. We treat lensing magnification consistently within Prospector , in contrast to the conventional approach where models are corrected for magnification post-fit.  \nThis paper is accompanied by a catalog, which is derived from the photometric super-catalog assembled using adaptive aperture sizes (Weaver et al. 2024). The UNCOVER SPS catalog contains the 16th, 50th, and 84th quantiles of the posterior distributions modeled with Prospector -β for the following parameters: redshift, stellar mass, mass-weighted age, SFR, sSFR, stellar metallicity, optical depth of diffuse dust, ratio between the optical depth of birth cloud dust and diffuse dust, power law index for a Calzetti et al. (2000) attenuation curve, and ratio between the object's AGN luminosity and its bolometric luminosity. Analytic estimates of magnifications, based on the redshift posterior medians, along with radial and tangential magnifications and shears are supplied for detailed lensing analyses (Furtak et al. 2023). Complementary redshifts from EAzY are included as well. An explanation of the catalog columns can be found in Table 2. We additionally provide the maximum-likelihood spectra, SFHs, and full posterior distributions. The catalog and all related documentation are accessible via the UNCOVER survey webpage 2 with a copy deposited to Zenodo: doi:10.5281/zenodo.8401181.  \nFuture updates to the stellar population catalog are expected in accordance with releases of extended photometric and/or spectroscopic data. The main planned enhancements in parameter inference are the incorporation of a number density prior, and lensing magnification as a free parameter with informative priors from the lens model.  \nWe conclude by reiterating that the UNCOVER survey provides the deepest view into our universe to date  \nby targeting the strong lens cluster Abell 2744. The stellar population catalog presented in this paper offers a rich data set to explore the processes governing galaxy formation and evolution over a redshift range now accessible by JWST.", 'ACKNOWLEDGMENTS': "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. (AURA), under NASA contract NAS 5-03127 for JWST. These observations are associated with JWST-GO-2561, JWST-ERS-1324, and JWST-DD-2756. Support for program JWST-GO2561 was provided by NASA through a grant from the Space Telescope Science Institute under NASA contract NAS 5-26555. This research is also based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute under NASA contract NAS 5-26555. These observations are associated with programs HST-GO-11689, HST-GO-13386, HST-GO/DD-13495, HST-GO-13389, HST-GO-15117, and HST-GO/DD-17231. The specific observations included in the JWST/UNCOVER first epoch survey data can be accessed via 10.17909/zn4s0243, whereas the full list of HST/JWST observations used to produce our UNCOVER catalogs is available via 10.17909/nftp-e621.  \nBW thanks the Penn State Institute for Gravitation and the Cosmos for financial support during the course of this work. RB acknowledges support from the Research Corporation for Scientific Advancement (RCSA) Cottrell Scholar Award ID No: 27587. LF and AZ acknowledge support by Grant No. 2020750 from the United States-Israel Binational Science Foundation (BSF) and Grant No. 2109066 from the United States National Science Foundation (NSF), and by the Ministry of Science & Technology, Israel. PD acknowledges support from the NWO grant 016.VIDI.189.162 ('ODIN') and from the European Commission's and University of Groningen's CO-FUND Rosalind Franklin program. HA acknowledges support from CNES (Centre National d'Etudes Spatiales). RS acknowledges an STFC Ernest Rutherford Fellowship (ST/S004831/1). MS acknowledges support from the CIDEGENT/2021/059 grant, from project PID2019109592GB-I00/AEI/10.13039/501100011033 from the Spanish Ministerio de Ciencia e Innovaci'on - Agen-", 'Wang et al.': "Table 2. Catalog Columns  \nFigure A.1. Comparison between the emulator and the FSPS redshifts on a sample of 1,000 randomly drawn objects from the full UNCOVER catalog. These plots suggest that the emulator is well calibrated, and the difference is captured in the uncertainties reported in the catalog. (a) Shown here is the distribution of the difference between the emulator and FSPS redshifts normalized by the 1 σ width of the FSPS posterior distribution. A unit Gaussian is over-plotted as a dotted line. The fraction of > 2 σ outliers is 0.015. (b) Medians of redshift posteriors from one run using FSPS are plotted against these from a second run under identical settings. Deviation from the diagonal reflects the sampling uncertainty arising from its probabilistic nature. (c) Medians of the redshift posteriors inferred using the emulator are plotted against the medians inferred using FSPS. Only results from a subset where S/N ≥ 10 in the LW detection bands are shown. (d) Same as the third panel, but the full validation set is included. Error bars are the 16th-84th quantiles. The Lyman break confusion is indicated in blue, and the shading shows ± 0.2 in redshift.  \n<!-- image -->  \ncia Estatal de Investigaci'on, and from Proyecto ASFAE/2022/025 del Ministerio de Ciencia y Innovaci'on en el marco del Plan de Recuperac'on, Transformaci'on y Resiliencia del Gobierno de Espa˜na. The work of CCW is supported by NOIRLab, which is managed by the AURA under a cooperative agreement with the NSF.  \nComputations for this research were performed on the Pennsylvania State University's Institute for Computational and Data Sciences' Roar supercomputer. This publication made use of the NASA Astrophysical Data System for bibliographic information.  \nFacilities: HST(ACS,WFC3), JWST(NIRCam, NIRSpec)  \nSoftware: Astropy (Astropy Collaboration et al. 2013, 2018, 2022), Corner (Foreman-Mackey 2016), EAzY (Brammer et al. 2008), Dynesty (Speagle 2020), Matplotlib (Hunter 2007), NumPy (Harris et al. 2020), Prospector (Johnson et al. 2021), SciPy (Virtanen et al. 2020)", 'A. ACCURACY OF PARAMETERS INFERRED WITH THE NEURAL NET EMULATOR': 'This work is the first large-scale application of the neural net emulator for stellar population synthesis, dubbed parrot (Mathews et al. 2023), which necessitates an assessment of its accuracy in inferring parameters. Mathews et al. (2023) discusses this accuracy in both flux space and in inferred parameter space but only with fixed redshifts; a free redshift significantly complexifies the parameter space by adding multiple new posterior modes, and so in this Appendix we examine the accuracy of the emulator in this new mode. We emphasize that at peak performance, the neural net emulator will exactly replicate the FSPS results, modulo any uncertainty introduced during the sampling process by the probabilistic nature of sampling.  \nWe randomly draw 1,000 objects from the full UNCOVER sample and refit them with FSPS (Conroy & Gunn 2010). The accuracy of the most important and challenging parameter, redshift, is assessed in four aspects summarized in Figure A.1. First, we calculate the difference between the emulator and FSPS redshifts normalized by the 1 σ width in the FSPS posterior distribution. In general, we find the emulator performs well in this test. The great majority of the inferred photometric redshifts are within 1 σ of the full fits, with the fraction of > 2 σ outliers being 0.015. Although the ideal performance is a δ -function, the actual distribution of normalized residuals suggests that the emulator is well calibrated. Second, we run FSPS twice with identical settings to estimate the sampling uncertainty. The sampling  \nmethod used is nested sampling (Skilling 2004), which is better suited to sample multi-modal posteriors than other traditional techniques such as Markov chain Monte Carlo (Goodman & Weare 2010). However, it has also been shown that nested sampling does not always accurately sample the global minimum when fitting for galaxy redshifts (Wang et al. 2023a). Our finding is consistent with this earlier study. In principle, the mode-finding problem in nested sampling can be avoided by substantial increases in the accuracy settings and thus the number of models called. This comes at a cost of increased CPU-hours, which quickly becomes prohibitively expensive for large-scale applications. Here we do not attempt to completely remove the sampling error, but instead adopt realistic dynesty (Speagle 2020) settings already more strict than those typically used in the literature, which keeps the time per fit roughly under an hour. Third, we compare the emulator and FSPS redshifts on a subsample in which the LW detection bands have S/N ≥ 10. This allows for an examination of the emulator performance without the additional complication from data quality issues. The scatter and outlier fraction are comparable to the expectations from the sampling uncertainty, indicating that the emulator agrees well with FSPS on highly confident photometric data. Fourth, we compare the performance on the full validation set. The outlier fraction increases marginally, with the catastrophic outliers roughly follow the lines indicating the confusion between a Lyman and a Balmer break.  \nGiven the finding above, we infer that the difference in the redshift posterior distributions found by the emulator and FSPS is mainly caused by a combination of two systematics. First, the emulated fluxes on average have 1-4% errors compared to the exact FSPS fluxes (Mathews et al. 2023). Second, the sampler struggles to assign correct posterior mass due to the complex likelihood surface. Fitting early JWST data exacerbates this challenge; in particular in light of the newly emerged model mismatch problem suggesting photometric calibrations are yet to be complete.  \nWith that being said, the uncertainties on the redshifts in Figure A.1 accurately capture the true residuals even in the presence of emulation errors. Thus, the difference between the emulator and FSPS results are in general well described by the reported uncertainties in the catalog. An additional valuable, independent check on the fidelity of the inferred parameters can likely be obtained via an alternative inference technique, since the latter is affected by different systematics (Wang et al. 2023a), or by comparing to additional spectra. This approach will be examined in the next generations of stellar population catalogs.', 'B. PHOTOMETRIC RESIDUALS': 'An examination of the residuals between the observed photometry and its maximum-likelihood model offers important insights into the agreement between data and the model. Figure B.1 shows fractional and uncertainty-normalized residuals in the observed and the rest frames. We find general agreement between Prospector -β and EAzY fits in the observed frame. Deviations from the unit Gaussian in the uncertainty-normalized residuals suggest that the photometry and the associated uncertainty require further refinement. The rest-frame residuals from Prospector -β fits fluctuate around 0, which suggests that there is no significant bias in our model templates. The large spikes are likely due to redshift traps, since the model spectra have to be shifted by the inferred redshifts. In contrast, systematic offsets in rest-frame residuals at short wavelengths are observed in EAzY -sfhz and EAzY -fsps fits, which may indicate template mismatches. However, given the on-going efforts in the photometric calibration, we caution against an over-interpretation of the EAzY performance at the current stage.', 'C. REST-FRAME COLORS': 'The rest-frame colors released in the catalogs are marginalized over the full Prospector -β posteriors. Here we additionally show the colors estimated from maximum-likelihood spectra in Figure C.1, which have larger scatter as expected. We also include the colors from EAzY -sfhz . EAzY measures the colors directly from the best-fit template, which is simply a linear combination of templates that minimizes the χ 2 statistics. Less information is encoded in the process determining the colors this way; it is thus understandable that the EAzY colors exhibit a larger scatter.', 'REFERENCES': 'Acebron, A., Jullo, E., Limousin, M., et al. 2017, MNRAS,  \n470, 1809, doi: 10.1093/mnras/stx1330  \nAcebron, A., Cibirka, N., Zitrin, A., et al. 2018, ApJ, 858,  \n42, doi: 10.3847/1538-4357/aabe29  \nAdams, N. J., Conselice, C. J., Ferreira, L., et al. 2023, MNRAS, 518, 4755, doi: 10.1093/mnras/stac3347  \nAlsing, J., Peiris, H., Mortlock, D., Leja, J., & Leistedt, B. 2023, ApJS, 264, 29, doi: 10.3847/1538-4365/ac9583  \n<!-- image -->  \n-  \n<!-- image -->  \n-  \n<!-- image -->  \n-  \n/F  \n)  \nF  \nF  \n(  \n<!-- image -->  \n-  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n-  \n-  \n<!-- image -->  \n-  \n<!-- image -->  \nFigure B.1. The distribution of photometric residuals resulting from Prospector -β and EAzY fits. We impose a minimum 5% error floor on the photometry. The upper panel shows the following from left to right. (a) Fractional residuals in observed frame, calculated as (modeled flux-observed flux)/(modeled flux) in the F444W band; data is clipped to ± 1. (b) Uncertainty-normalized residuals in observed frame, calculated as (modeled flux-observed flux)/(uncertainty in observed flux) in the F444W band; data is clipped to ± 3 σ . The gray dashed line indicates a unit Gaussian for reference. (c) Fractional residuals in rest frame, binned logarithmically in wavelengths. (d) Error-normalized residuals in rest frame, also binned logarithmically in wavelength. Shaded regions indicate the 16-84th and the 2.5-97.5th percentiles. The middle and lower panels show the EAzY -sfhz and EAzY -fsps residuals in the same manner. The Prospector -β and EAzY fits generally agree. However, with JWST photometric calibrations still ongoing, we caution against an over-interpolation on the performance at the current stage.  \n<!-- image -->  \n<!-- image -->  \nFigure C.1. Here we compare the rest-frame colors of galaxies at 2 . 5 < z < 3 . 5 inferred by three different methods: marginalizing over the full Prospector -β posteriors, colors inferred from the maximum-likelihood (ML) spectra from Prospector -β , and from the best-fit templates of EAzY -sfhz . The medians of the uncertainty distributions are indicated by black error bars, while the 90% quantiles of the uncertainty distributions are indicated by gray error bars. The asymmetry of the error bars is driven by the non-Gaussianity of the redshift posterior distribution. Only objects of which S/N ≥ 10 in the LW detection bands are included for the ease of interpreting the difference. The marginalized colors have a smaller scatter because marginalization tamps down the photometric uncertainties, but their uncertainties are reflective of the larger scatters found by the other two methods.  \n<!-- image -->  \nAlsing, J., Peiris, H., Leja, J., et al. 2020, ApJS, 249, 5,  \ndoi: 10.3847/1538-4365/ab917f Anders, P., & Fritze-v. Alvensleben, U. 2003, A&A, 401, 1063, doi: 10.1051/0004-6361:20030151 Antwi-Danso, J., Papovich, C., Leja, J., et al. 2023, ApJ, 943, 166, doi: 10.3847/1538-4357/aca294 Arrabal Haro, P., Dickinson, M., Finkelstein, S. L., et al. 2023, ApJL, 951, L22, doi: 10.3847/2041-8213/acdd54 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy 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 Atek, H., Richard, J., Kneib, J.-P., & Schaerer, D. 2018, MNRAS, 479, 5184, doi: 10.1093/mnras/sty1820 Atek, H., Chemerynska, I., Wang, B., et al. 2023, MNRAS, 524, 5486, doi: 10.1093/mnras/stad1998 Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143, doi: 10.1093/mnras/stz1182 Bellstedt, S., Robotham, A. S. G., Driver, S. P., et al. 2021, MNRAS, 503, 3309, doi: 10.1093/mnras/stab550 Bezanson, R., Labbe, I., Whitaker, K. E., et al. 2022, arXiv e-prints, arXiv:2212.04026. https://arxiv.org/abs/2212.04026 Boker, T., Beck, T. L., Birkmann, S. M., et al. 2023, PASP,  \n135, 038001, doi: 10.1088/1538-3873/acb846  \nBouwens, R. J., Oesch, P. A., Illingworth, G. D., Ellis, R. S., & Stefanon, M. 2017, ApJ, 843, 129,  \ndoi: 10.3847/1538-4357/aa70a4  \nBoyer, M. L., Anderson, J., Gennaro, M., et al. 2022, Research Notes of the American Astronomical Society, 6, 191, doi: 10.3847/2515-5172/ac923a  \nBrammer, G. B., van Dokkum, P. G., & Coppi, P. 2008,  \nApJ, 686, 1503, doi: 10.1086/591786  \nBrammer, G. B., Whitaker, K. E., van Dokkum, P. G., et al. 2011, ApJ, 739, 24,  \ndoi: 10.1088/0004-637X/739/1/24  \nBrammer, G. B., S´anchez-Janssen, R., Labb´e, I., et al. 2012, ApJL, 758, L17, doi: 10.1088/2041-8205/758/1/L17 Byler, N., Dalcanton, J. J., Conroy, C., & Johnson, B. D. 2017, ApJ, 840, 44, doi: 10.3847/1538-4357/aa6c66  \nByrne, C. M., Stanway, E. R., Eldridge, J. J., McSwiney, L., & Townsend, O. T. 2022, MNRAS, 512, 5329,  \ndoi: 10.1093/mnras/stac807  \nCalzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682, doi: 10.1086/308692  \nCarnall, A. C., Begley, R., McLeod, D. J., et al. 2023,  \nMNRAS, 518, L45, doi: 10.1093/mnrasl/slac136  \nChabrier, G. 2003, PASP, 115, 763, doi: 10.1086/376392  \nCharlot, S., & Fall, S. M. 2000, ApJ, 539, 718,  \ndoi: 10.1086/309250  \nChoi, J., Conroy, C., & Johnson, B. D. 2019, ApJ, 872, 136, doi: 10.3847/1538-4357/aaff67  \n- Vazdekis, A., Coelho, P., Cassisi, S., et al. 2015, MNRAS, 449, 1177, doi: 10.1093/mnras/stv151\n- Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, NatMe, 17, 261, doi: 10.1038/s41592-019-0686-2\n- Walcher, J., Groves, B., Budav\'ari, T., & Dale, D. 2011, Ap&SS, 331, 1, doi: 10.1007/s10509-010-0458-z\n- Wang, B., Leja, J., Villar, V. A., & Speagle, J. S. 2023a, ApJL, 952, L10, doi: 10.3847/2041-8213/ace361\n- Wang, B., Fujimoto, S., Labb\'e, I., et al. 2023b, ApJL, 957, L34, doi: 10.3847/2041-8213/acfe07\n- Wang, B., Leja, J., Bezanson, R., et al. 2023c, ApJL, 944, L58, doi: 10.3847/2041-8213/acba99\n- Wang, B., Leja, J., Atek, H., et al. 2024, ApJ, 963, 74, doi: 10.3847/1538-4357/ad187c\n- Weaver, J. R., Cutler, S. E., Pan, R., et al. 2024, ApJS, 270, 7, doi: 10.3847/1538-4365/ad07e0\n- Webb, K., Balogh, M. L., Leja, J., et al. 2020, MNRAS, 498, 5317, doi: 10.1093/mnras/staa2752\n- Whitler, L., Stark, D. P., Endsley, R., et al. 2023, MNRAS, 519, 5859, doi: 10.1093/mnras/stad004\n- Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & Labb\'e, I. 2009, ApJ, 691, 1879, doi: 10.1088/ 0004-637X/691/2/187910.48550/arXiv.0806.0625\n- Yung, L. Y. A., Somerville, R. S., Finkelstein, S. L., Wilkins, S. M., & Gardner, J. P. 2024, MNRAS, 527, 5929, doi: 10.1093/mnras/stad3484\n- Zavala, J. A., Buat, V., Casey, C. M., et al. 2023, ApJL, 943, L9, doi: 10.3847/2041-8213/acacfe\n- Zitrin, A., Fabris, A., Merten, J., et al. 2015, ApJ, 801, 44, doi: 10.1088/0004-637X/801/1/44'}