In this paper we describe the Medium Bands Mega Science JWST Cycle 2 survey JWSTGO4111 and demonstrate the power of these data to reveal both the spatially integrated and spatially resolved properties of galaxies from the local Universe to the era of cosmic dawn. Executed in 2023 November MegaScience obtained 30 arcminSUP2SUP of deep multiband NIRCam imaging centered on the z  0.3 A2744 cluster including 11 mediumband filters and the two shortestwavelength broadband filters F070W and F090W. Together MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep 2830 magSUBABSUB images in all NIRCam medium and broadband filters. This unique data set allows us to precisely constrain photometric redshifts map stellar populations and dust attenuation for large samples of distant galaxies and examine the connection between galaxy structures and formation histories. MegaScience also includes 17 arcminSUP2SUP of NIRISS parallel imaging in two broadband and four mediumband filters from 0.9 to 4.8 m expanding the footprint where robust spectral energy distribution SED fitting is possible. We provide example SEDs and multiband cutouts at a variety of redshifts and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 23. Additionally we demonstrate the spatially resolved science enabled by MegaScience by presenting maps of the O III line emission and continuum emission in three spectroscopically confirmed z gt 6 galaxies. We show that line emission in reionizationera galaxies can be clumpy extended and spatially offset from continuum emission implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields jwstuncover.github.iomegascience.
['Galaxy evolution', 'Galaxy formation', 'Galaxy structure', 'High-redshift galaxies', '594', '595', '622', '734', 'Astrophysics - Astrophysics of Galaxies']
Medium Bands Mega Science A JWSTNIRCam Mediumband Imaging Survey of A2744
{'Medium Bands, Mega Science: a JWST/NIRCam Medium-Band Imaging Survey of Abell 2744': "Katherine A. Suess, 1, ∗ John R. Weaver, 2 Sedona H. Price, 3 Richard Pan, 4 Bingjie Wang ( 王 冰 洁 ), 5, 6, 7 Rachel Bezanson, 3 Gabriel Brammer, 8, 9 Sam E. Cutler, 2 Ivo Labb'e, 10 Joel Leja, 5, 6, 7 Christina C. Williams, 11 Katherine E. Whitaker, 2, 12 Pratika Dayal, 13 Anna de Graaff, 14 Robert Feldmann, 15 Marijn Franx, 16 Yoshinobu Fudamoto, 17 Seiji Fujimoto, 18, ∗ Lukas J. Furtak, 19 Andy D. Goulding, 20 Jenny E. Greene, 20 Gourav Khullar, 3 Vasily Kokorev, 13 Mariska Kriek, 16 Brian Lorenz, 21 Danilo Marchesini, 22 Michael V. Maseda, 23 Jorryt Matthee, 24 Tim B. Miller, 25 Ikki Mitsuhashi, 26 Lamiya A. Mowla, 27 Adam Muzzin, 28 Rohan P. Naidu, 29, ∗ Themiya Nanayakkara, 30 Erica J. Nelson, 31 Pascal A. Oesch, 32, 8 David J. Setton, 20, † Heath Shipley, 33 Renske Smit, 34 Justin S. Spilker, 35 Pieter van Dokkum, 36 and Adi Zitrin 19  \n1 Kavli Institute for Particle Astrophysics and Cosmology and Department of Physics, Stanford University, Stanford, CA 94305, USA 2 Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 3 Department of Physics and Astronomy and PITT PACC, University of Pittsburgh, Pittsburgh, PA 15260, USA 4 Department of Physics and Astronomy, Tufts University, 574 Boston Ave., Medford, MA 02155, USA 5 Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 6 Institute for Computational & Data Sciences, The Pennsylvania State University, University Park, PA 16802, USA 7 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 8 Cosmic Dawn Center (DAWN), Copenhagen, Denmark 9 Niels Bohr Institute, University of Copenhagen, Jagtvej 128, Copenhagen, Denmark 10 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Melbourne, VIC 3122, Australia 11 NSF's National Optical-Infrared Astronomy Research Laboratory, 950 North Cherry Avenue, Tucson, AZ 85719, USA 12 Cosmic Dawn Center (DAWN), Denmark 13 Kapteyn Astronomical Institute, University of Groningen, 9700 AV Groningen, The Netherlands 14 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117, Heidelberg, Germany 15 Department of Astrophysics, University of Zurich, Zurich, CH-8057, Switzerland 16 Leiden Observatory, Leiden University, P.O.Box 9513, NL-2300 AA Leiden, The Netherlands 17 Center for Frontier Science, Chiba University, 1-33 Yayoi-cho, Inage-ku, Chiba 263-8522, Japan 18 Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA 19 Department of Physics, Ben-Gurion University of the Negev, P.O. Box 653, Be'er-Sheva 84105, Israel 20 Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA 21 Department of Astronomy, University of California, Berkeley, CA 94720, USA 22 Department of Physics & Astronomy, Tufts University, MA 02155, USA 23 Department of Astronomy, University of Wisconsin-Madison, 475 N. Charter St., Madison, WI 53706, USA 24 Institute of Science and Technology Austria (ISTA), Am Campus 1, 3400 Klosterneuburg, Austria 25 Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA), Northwestern University, 1800 Sherman Ave, Evanston IL 60201, USA 26 Department of Astronomy, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 27 Dunlap Institute for Astronomy & Astrophysics, Toronto, Ontario, M5S 3H4, Canada 28 Department of Physics and Astronomy, York University, 4700 Keele Street, Toronto, ON, M3J 1P3, Canada 29 MIT Kavli Institute for Astrophysics and Space Research, 77 Massachusetts Ave., Cambridge, MA 02139, USA 30 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, VIC 3122, Australia 31 Department for Astrophysical & Planetary Science, University of Colorado, Boulder, CO 80309, USA 32 Department of Astronomy, University of Geneva, Chemin Pegasi 51, 1290 Versoix, Switzerland 33 Department of Physics, Texas State University, San Marcos, TX 78666, USA 34 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK 35 Department of Physics and Astronomy and George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, Texas A&M University, 4242 TAMU, College Station, TX 77843-4242 36 Department of Astronomy, Yale University, 52 Hillhouse Ave., New Haven, CT, 06511, USA", 'ABSTRACT': "suess@stanford.edu  \nIn this paper, we describe the 'Medium Bands, Mega Science' JWST Cycle 2 survey (JWST-GO4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatiallyresolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ∼ 30arcmin 2 of deep multiband NIRCam imaging centered on the z ∼ 0 . 3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep ( ∼ 28 -30mag AB ) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ∼ 17 arcmin 2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0 . 9 -4 . 8 µ m, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multiband cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z > 6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields (jwst-uncover.github.io/megascience).  \nKeywords: Galaxy evolution (594) - Galaxy formation (595) - Galaxy structure (622) - Highredshift galaxies (608)", '1. INTRODUCTION': "The James Webb Space Telescope allows us, for the first time, to map the rest-frame optical and infrared emission from faint galaxies at z ≳ 3. In the first ∼ 1.5 years of JWST data, these new capabilities have already begun to transform our understanding of galaxy formation and evolution. Luminous z > 10 galaxies appear to be ≳ 10 times more common than expected from theoretical predictions (e.g., Naidu et al. 2022a; Adams et al. 2023; Atek et al. 2023; Austin et al. 2023; Bradley et al. 2023; Casey et al. 2023; Finkelstein et al. 2023; Robertson et al. 2023a). Galaxy formation may have begun more rapidly than expected, with massive ( ≳ 10 10 M ⊙ ) galaxies reported very early in the universe ( z ≳ 7) (e.g., Labb'e et al. 2023; Boylan-Kolchin 2023; Akins et al. 2023; Xiao et al. 2023) and massive quiescent galaxies emerging as early as z ∼ 4 . 6 (e.g., Carnall et al. 2023; Valentino et al. 2023; Setton et al. 2024). New classes of objects - particularly 'little red dots', thought to be predominantly highly dust-obscured active galactic nuclei - appear to be remarkably common in deep-field observations (e.g., Matthee et al. 2023b; Labbe et al. 2023; Furtak et al. 2023a; Greene et al. 2023; Williams  \net al. 2023a; Barro et al. 2023; Kokorev et al. 2023, 2024; Wang et al. 2024a).  \nDue to the targeted nature of spectroscopy, many of these initial discoveries rely on accurate interpretation of broad-band photometry. Spectroscopy is also significantly more expensive than imaging, meaning that photometry is also often necessary to select statistically large samples. While spectral energy distribution (SED) modeling of broad-band photometry is able to efficiently infer the physical properties of large samples of galaxies, the coarse wavelength resolution can permit contradictory physical interpretations. For example, models of photometric data alone can confuse the Lyman and Balmer breaks, obtaining dramatically incorrect photometric redshift solutions (e.g., Donnan et al. 2023; Arrabal Haro et al. 2023; Naidu et al. 2022b; Zavala et al. 2023). It can be nearly impossible to correctly disentangle the relative contribution of line and continuum emission to broad-band photometry, leading to incorrect photometric redshifts or biased stellar masses (e.g., Laporte et al. 2023; Sarrouh et al. 2024). Broad-band photometric modeling may fail to identify new classes of objects with physical properties that lie beyond the modeling assumptions. For example, standard SED fitting yielded high stellar masses and redshifts for the objects in Labb'e et al. (2023), but later spectroscopy revealed many of these objects to be a new population of reddened AGN that were not included in our pre-JWST  \nphotometric redshift template sets (e.g., Kocevski et al. 2023; Greene et al. 2023).  \nThe task of accurately inferring the physical properties of galaxies from photometry is further complicated by 'outshining': young stars so completely dominate the integrated light from galaxies that it is easy to miss up to a factor of ten more mass hidden in low-luminosity old stars (e.g., Papovich et al. 2001; for the predicted impact on z > 7 galaxies with JWST see Whitler et al. 2023; Narayanan et al. 2024; Wang et al. 2024b). Spatially-resolved SED fitting studies, which can begin to disentangle line- and continuum-dominated regions, often find up to 40% larger stellar masses and 0.3 dex older ages than spatially-integrated studies (e.g., Zibetti et al. 2009; Wuyts et al. 2012; Sorba & Sawicki 2015). These differences between integrated and resolved properties may be even more dramatic in the high-redshift universe, where high equivalent widths of 1000-3000 ˚ A are common and can bias stellar masses and ages by up to a full dex (e.g., Gim'enez-Arteaga et al. 2023; see also P'erez-Gonz'alez et al. 2023).  \nBuilding an accurate picture of a large sample of galaxies in the distant universe with JWST therefore requires spatially-resolved observations that can distinguish between nebular line and stellar continuum emission. These spatially- and spectrally-resolved observations provide a viable pathway forward towards breaking these model degeneracies and thereby producing accurate redshifts, stellar masses, star formation rates, and other stellar population properties. The most efficient way to achieve these results is to exploit one of JWST/NIRCam's most innovative features: a suite of twelve medium bands covering the full 1.5-5 µ m range (Rieke et al. 2023). These medium-band filters are effective at mapping strong emission lines, but they also have the continuum sensitivity required to probe older stellar populations in line-free regions. Together, these medium-band data allow us to efficiently map both the mass and approximate formation time of thousands of galaxies simultaneously.  \nGround-based observations have long shown the power of optical and near-infrared medium bands to finely sample the Balmer break at z ∼ 0 -4 and constrain the physical properties of distant galaxies (e.g., COMBO-17, Wolf et al. 2003; NMBS, Whitaker et al. 2011; SHARDS, P'erez-Gonz'alez et al. 2013; ZFOURGE, Straatman et al. 2016; FENIKS, Esdaile et al. 2021). While most early JWST programs used only a single medium-band filter, F410M, several more recent programs have begun to leverage the power of NIRCam's multiple medium bands: the JWST Extragalactic Medium-band Survey includes five 1-4 µ m medium bands (Williams et al.  \n2023b); parallel observations conducted as part of the CAnadian NIRISS Unbiased Cluster Survey include nine 1-4 µ m medium bands, with the three missing medium bands forthcoming in PID 3362 (PI Muzzin), and the Jades Origin Field includes six 1-3 µ m medium bands (Eisenstein et al. 2023). These datasets are already beginning to demonstrate the power of multiple medium bands to constrain the redshifts, stellar ages, masses, and emission line equivalent widths of very large samples of galaxies (e.g., Withers et al. 2023).  \nHere we present the 'Medium Bands, Mega Science' survey, designed to image the full Abell 2744 cluster field with all dozen of JWST/NIRCam's medium bands. Combined with existing broad-band imaging from the Cycle 1 UNCOVER Treasury Program (Bezanson et al. 2022), these observations make up a complete set of deep images in all of NIRCam's broad- and medium-band filters. This 20-band JWST imaging results in spatiallyresolved 'ultra low-resolution' R ∼ 15 spectrophotometry for all ∼ 70 , 000 sources in the field, with nearly continuous photometric coverage from 0 . 7 -5 µ m.  \nSection 2 describes our survey design, data reduction, and cataloguing methods. Section 3 describes our science goals, including both spatially-integrated science (including example SEDs; Section 3.1) and spatiallyresolved science (including example multi-band color images; Section 3.2). Section 4 provides a demonstration of the power of our data to localize strong line emission in distant galaxies, presenting evidence for clumpy and offset [OIII]+H β line emission in z spec ∼ 6 . 3 galaxies.  \nAlong with this publication, we publicly release full reduced mosaics and a photometric catalog. Our website, jwst-uncover.github.io/megascience, has links to all data products, and the associated Zenodo release is located at zenodo.org/doi/10.5281/zenodo.8199802.  \nThroughout this paper we assume a standard ΛCDM cosmology with Ω m = 0 . 3, Ω Λ = 0 . 7, and h = 0 . 7. Magnitudes are quoted in the AB system (Oke 1974). All ID numbers given correspond to the DR3 UNCOVER+MegaScience catalogs released along with this paper.", '2.1. Survey Design & Field': "The primary goal of the Medium Bands, Mega Science survey ('MegaScience'; PID 4111) is to obtain deep (2.3-4.6hr) imaging with the full suite of NIRCam medium-band filters over the same area as the UNCOVER survey (Bezanson et al. 2022). These mediumband data significantly increase the spectral resolution of existing broad-band observations, and allow us to perform a wealth of both spatially-integrated science (e.g.,  \nmore accurate redshifts, stellar masses, and star formation rates; see Section 3.1) and spatially-resolved science (e.g., maps of line emission and physical properties; see Section 3.2).  \nThe Abell 2744 field hosts a remarkable compilation of public data, including Hubble Space Telescope imaging through the Hubble Frontier Fields program (Lotz et al. 2017) and the BUFFALO survey (Steinhardt et al. 2020), deep MUSE spectroscopy (Mahler et al. 2018), and ALMA coverage (Fujimoto et al. 2023). Additionally, the field has been targeted by a number of early JWST programs including UNCOVER (#2561; PIs Labb'e & Bezanson), the GLASS Early Release Science program (Treu et al. 2022), DDT program #2756 (PI: Chen; Chen et al. 2022), a Cycle 2 pure parallel imaging program (#3990, PI Morishita), parallel NIRCam imaging to Cycle 2 GLASS spectroscopy (#3073; PI Castellano), as well as three Cycle 2 programs obtaining both NIRCam/grism and NIRCam/imaging (ALT, #3516, PIs Matthee & Naidu; MAGNIF, # 2883, PI Sun; and #3538, PI Iani). Footprints of pre-JWST data in Abell 2744 are shown in Bezanson et al. (2022), while footprints of JWST data including the MegaScience survey are shown in Figure 1. The vast majority of these JWST programs, including MegaScience, were taken with no proprietary time. All told, Abell 2744 has public data in both broad- and medium-band NIRCam imaging, NIRCam grism observations in multiple filters, and NIRSpec prism spectroscopy adding up to over 300hr of community investment in Cycle 1 & 2 alone. This wealth of data means that Abell 2744 is rapidly emerging as one of the preeminent public extragalactic JWST deep fields.  \nFurtak et al. (2023b) presents a state-of-the-art lensing model for Abell 2744 based off of the deep UNCOVER imaging. Out of the total ∼ 29 . 2 arcmin 2 area of our primary NIRCam pointing, 21.0 arcmin 2 (72%) has magnification µ > 1 . 5, 10.5 arcmin 2 (36%) has µ > 2, 2.3 arcmin 2 (8%) has µ > 5, and 1.1 arcmin 2 (4%) has µ > 10; nearly every source in the footprint has µ ≳ 1 . 3 1 . This lensing boosts our effective depth, allowing us to detect intrinsically fainter sources than a field survey with the same integration time. While the majority of moderately-magnified sources in the field are not significantly morphologically distorted, the small regions of high magnification near the cluster cores boast very highly magnified objects, allowing us to improve the effective spatial resolution by factors of up to ∼ 5 -10.  \nAs shown in Figure 2, the primary NIRCam pointing of MegaScience targets all medium- and broadband filters not observed as part of UNCOVER: eleven medium-band filters (F140M, F162M, F182M, F210M, F250M, F300M, F335M, F360M, F430M, F460M, and F480M; F410M was observed by UNCOVER) as well as two blue broad-band filters F070W and F090W. These blue broad-band filters can be used as dropout filters to select high-redshift candidates and are observed 'for free' during observations of NIRCam's numerous long-wavelength medium-band filters. The NIRISS parallel pointing fills in the existing broad-band UNCOVERphotometric coverage with F090W and F277W, and additionally adds four medium-band filters F140M, F158M, F430M, and F480M.  \nMegaScience was originally scheduled for four visits on November 5, 11, and 12, 2023. The first visit on November 5th failed due to guide star acquisition, and was rescheduled for November 20th with WOPR 88967. In order to obtain a suitable guide star and schedule the pointing in 2023, this rescheduled visit is slightly misaligned with the rest of the mosaic as shown in Figure 1, resulting in a slight rotation relative to the existing UNCOVER data.", '2.2. Data Reduction': "The MegaScience data are reduced following similar methods as described in Bezanson et al. (2022). All public NIRCam and NIRISS exposures in the Abell 2744 field are calibrated with the 'jwst 0995.pmap' calibration set after retrieval from MAST. This includes the MegaScience data discussed here (GO-4111, PI Suess), as well as imaging from UNCOVER (Bezanson et al. 2022), ERS-GLASS (Treu et al. 2022), DD-2756 (PI: Chen), MAGNIF (GO-2883, PI: Sun), GO-3538 (PI: Iani), and ALT (GO-3516; Naidu & Matthee in prep.). Images are aligned and co-added with the Grizli pipeline (G. Brammer in prep., version 1.9.13.dev26), which uses the updated snowblind code 2 for masking NIRCam and NIRISS snowballs and an updated NIRCam bad pixel mask. These images correspond to the version 7.2 mosaic release hosted and described in additional detail on the DAWN JWST Archive (DJA; https://dawn-cph.github.io/dja/imaging/v7/).  \nAs in UNCOVER, the brightest cluster galaxies (bCG) and associated intracluster light (ICL) are modelled and subtracted following Shipley et al. (2018) and described in detail in Weaver et al. (2024). A subsequent background subtraction is also applied, using the same methodologies as the broad-band UNCOVER images.  \nFigure 1. Schematic of the on-sky footprint of MegaScience and UNCOVER (left), other NIRCam and NIRISS imaging programs (center), and NIRCam grism programs (right); these programs are described in additional detail in the main text. The MegaScience footprint is shown in light blue in all three panels. A comparison with footprints of pre-JWST data, including HST and MUSE, is shown in Bezanson et al. (2022). For the UNCOVER parallel observations in the southeast of the image, we only show pointings which are included in the reduced mosaics released along with this paper (declination < 30 · 28).  \n<!-- image -->  \nTable 1 presents the depths of the final, bCG/ICLsubtracted mosaics.", '2.3. Source Detection and Photometry': "Multi-wavelength photometric catalogs are constructed with aperpy 3 broadly following methods used to create the UNCOVER catalog (Weaver et al. 2024). Because the depth and area of long-wavelength broadband imaging in Abell 2744 has increased significantly over the course of Cycles 1 & 2, we construct a new 56.2 arcmin 2 detection image as a noise-weighted stack of F277W, F356W, and F444W after bright cluster galaxies (bCGs) are removed. Using the same SEP (Barbary 2016) configuration as Weaver et al. (2024), we detect a total of 74,020 objects. This is a ∼ 20% increase compared to the DR2 UNCOVER catalog.  \nWe measure photometry on PSF-matched images to account for the large variation in the PSF from different filters and instruments. While F480M has the broadest PSF in our survey, we choose to match to F444W for three reasons: (1) the F444W imaging is considerably deeper (by ∼ 1 . 8mag), allowing us to construct a robust empirical PSF; (2) F480M has a significantly more complex PSF than F444W, which would introduce additional complexity in the homogenization process; and (3) matching to F444W allows us to perform 1:1 comparisons with most other early JWST surveys, which typically match to F444W. We construct matching ker-  \nls to F444W using Pypher (Boucaud et al. 2016). Figure 3 shows the performance of our matching kernels for the filters observed in MegaScience, finding good agreement with the target F444W PSF below 0.1% for most bands. The worst performance is from the four 4 µ mkernels (F410M, F430M, F460M, F480M) due to the rapid changes in the PSF structure relative to F444W, but even these long-wavelength PSFs agree well within 1%.  \nFluxes are extracted in four apertures of 0.32', 0.48', 0.70', 1.00', and 1.40' diameter and corrected to total fluxes by scaling to 'AUTO' fluxes extracted from a PSFmatched version of our LW detection image in elliptical Kron apertures. We perform an additional correction for light outside the Kron radius (Kron 1980), following the methods of Weaver et al. (2024) (see also Whitaker et al. 2011; Skelton et al. 2014; Whitaker et al. 2019). Flux errors are estimated directly from the images using 10,000 apertures placed randomly on background regions. We construct a USE PHOT flag following similar methods as Weaver et al. (2024); this flag is intended to remove both stars and imaging artifacts. Similar to the UNCOVER catalog, the vast majority (91%) of objects have USE PHOT = 1 . For ease-of-use, we provide a 'SUPER' catalog wherein photometry is extracted in the optimal aperture for each object, taking extra care to avoid failures in crowded regions of these deep images. In practice, lower redshift galaxies with extended light profiles will default to larger apertures unless they are in densely crowded regions where it is better to use a smaller aperture to avoid contamination.  \nAs a check, we compare our newly extracted photometry to UNCOVER DR2 in all common, pre-MegaScience  \nFigure 2. Existing photometric coverage in both the NIRCam primary and NIRISS parallel fields, along with new imaging from this program. MegaScience 'completes' the imaging in the primary field, with deep images in all 20 of NIRCam's broadand medium-band filters. The addition of two broad-band and four medium-band filters in the NIRISS parallel field will allow for robust SED modeling.  \n<!-- image -->  \nbands, finding zero bias at all magnitudes and a scatter commensurate with expected flux uncertainties. In other words: the photometry in our catalog agrees well with the previous public release of the UNCOVER catalog in the broad bands, despite the new segmentation map. As discussed in Weaver et al., it is difficult to assess photometric reliability (especially for the new medium bands) without comparing to an independent catalog. Weaver et al. (2024) perform a sanity check by comparing the catalog fluxes to the predicted photometry from eazy 's best-fit model. We perform a similar test using the MegaScience catalog, and, like UNCOVER, we find good agreement from all bands. We refer the reader to the documentation accompanying these catalogs for up-to-date details as to its contents, column names, and conventions.", '2.4. Photometric Redshifts': "One of the motivations behind this survey is to significantly improve photometric redshift estimates ( z phot ) and related physical parameters. As with the original UNCOVER survey, our team will provide two sets of redshift and stellar population catalogs: z phot estimates fit using eazy-py (Brammer et al. 2008), and a stellar population catalog fit using Prospector (Johnson et al. 2021), following the methods described in Wang et al. (2024c).  \nWe estimate z phot by fitting all HST and JWST photometry in the default 'SUPER' catalog with eazy-py (Brammer et al. 2008), following the same methods as Weaver et al. (2024). We provide z phot estimates run with the SFHz BLUE AGN template set, which includes  \nTable 1. Measured MegaScience imaging depths  \nNOTE: Imaging depths are calculated in the same manner as Bezanson et al. (2022), and correspond to the two-visit depths regions of the mosaic; depths are calculated using 0.08' radius apertures in the SW bands and 0.16' radius apertures in the LW bands and corrected to total assuming point sources.  \nredshift-dependent priors to ensure that galaxies cannot be fit with ages older than the age of the Universe. These templates also include an AGN torus template from (Killi et al. 2023). Priors on magnitude and β slope are switched off. We append the estimated photometric redshifts and select rest-frame fluxes to our photometric catalogs. The photometric catalog, z phot estimates, and reduced bCG-subtracted images are publicly available on our website (jwst-uncover.github.io/megascience).  \nFor the Prospector catalog, we fit the same HST and JWST photometry, but now using the methods and physical model described in Wang et al. (2024c). The catalog of best-fit redshifts is released along with this paper; a catalog of best-fit stellar population parameters (masses, star formation rates, etc) will be presented in B. Wang et al. (in prep.) along with a detailed description of the improvements in physical parameters using the medium-band photometry.  \nFigure 3. PSF growth curves for each filter added to UNCOVER before (top) and after (bottom) matching to F444W. After matching, all filters have deviations below the 1% level at the smallest aperture diameter used (0.32 '' ). Growth curves are shown relative to the F444W growth curve; a value of 1 indicates perfect matching with F444W. Dashed lines indicate the ± 1% deviations from exact matching (solid black line). Dotted lines indicate the location of 0.32 '' and 0.70 '' aperture diameters.  \n<!-- image -->  \nThe performance of both sets of photometric redshifts is discussed in more detail in Section 3.1 and Figure 5.", '3. SCIENCE OBJECTIVES': "3.1. Spatially-integrated science: improving redshifts, stellar masses, and star formation rates  \nThe MegaScience dataset samples the SEDs of all galaxies in the field at nearly the highest spectral resolution possible with NIRCam imaging 4 . This spectral resolution allows us to improve the accuracy of recovered photometric redshifts and stellar population properties while breaking degeneracies introduced by fitting broad-  \nband photometry. In this section, we demonstrate the improvement MegaScience is able to make on the spectral energy distributions of galaxies compared to UNCOVER broad-band photometry.  \nFigure 4 shows the photometry and eazy-py fits to six example galaxies, all of which have spectroscopic redshifts from either UNCOVER (IDs 61218, 34713, 49587, and 13963; S. Price et al. in prep) or ALT (IDs 30876, 15243; J. Matthee & R. Naidu et al. in prep). The left column shows the fit to UNCOVER data alone (blue); the central column shows the new MegaScience photometry and the corresponding eazy-py fit (black); the right column shows the photometric redshift posteriors of both the UNCOVER and MegaScience fits compared with the spectroscopic redshift (red). Each galaxy was chosen to represent a way that the medium-band data resolves a common 'failure mode' of broad-band SED fitting.  \nID 61218 shows a star-forming galaxy at z ∼ 1. In this redshift regime, the MegaScience data is able to both pin down the location of the Balmer break (thanks to F070W & F090W data, especially important in the portions of the field lacking HST coverage), and identify the 1.6 µ m stellar bump with the long-wavelength medium-bands. Together, these features allow for precise photometric redshift recovery and effectively rule out high-redshift solutions that are allowed by the UNCOVER fit alone. We expect the 1.6 µ m stellar bump to also help identify cluster members in Abell 2744 itself at z ∼ 0 . 33. This is especially important when pushing down the mass function towards the lowest mass cluster members, where the 1.6 micron bump is effectively the only prominent feature identifiable in the absence of adequate coverage of the Balmer break (observable at < 0 . 5 micron).  \nID 34713 shows a galaxy at cosmic noon ( z ∼ 2 . 5) with relatively low SFR ∼ 1 -10 M ⊙ yr -1 , potentially a massive galaxy caught early in the quenching process. For this class of object, the MegaScience data is able to (a) confirm low rest-UV fluxes with new 0.7-1.5 µ mdata, indicating low ongoing SFRs, (b) precisely constrain the observed wavelength of the Balmer break, and (c) constrain the shape of the spectrum redward of the break. The location and shape of the Balmer break are particularly useful for characterizing the ages of quenching and quiescent galaxies (e.g., Carnall et al. 2019; Wild et al. 2020; Suess et al. 2022).  \nID 49587 shows an example where where the eazy-py fit using UNCOVER data alone prefers a Balmer break, but the MegaScience fit reveals the (accurate) Lyman break solution. The degeneracy between Balmer and Lyman breaks is one of the most common causes of low-  \nredshift interlopers in high-redshift photometric selections. In large part, the MegaScience data is able to break this degeneracy by identifying strong line emission boosting in several of the medium-band filters; no lines are predicted in these locations for the Balmerbreak model.  \nID 13963 shows another classic photometric redshift mixup: what is fit with a rising red continuum in the UNCOVER data is revealed to be boosting by line emission in the MegaScience data. The true continuum level can be seen in the F300M filter, which is one of the few line-free filters in the entire SED < 4 µ m. JWST medium bands have been shown to be highly effective in disentangling these line vs continuum solutions for large samples of distant galaxies (e.g., Laporte et al. 2023; Desprez et al. 2023; Sarrouh et al. 2024).  \nID 30876 demonstrates the improvement in both accuracy and precision of photometric redshifts enabled by medium-band observations of sources with strong line emission. This galaxy, confirmed at z = 5 . 66 with ALT grism observations (J. Matthee & R. Naidu et al. in prep), had a relatively broad photometric redshift posterior with UNCOVER alone ( δz ∼ 0 . 43): the fit was unable to precisely constrain the location of the Lyman and Balmer breaks, and it was unclear where exactly the H β +[OIII] lines were. The MegaScience data nails down the locations of the Lyman break, the Balmer break, H β +[OIII], and H α , resulting in an accurate photometric redshift with an uncertainty of just ∼ 0.02.  \nID 15243 demonstrates a uniquely-JWST photometric redshift mixup: the H β +[OIII] and H α lines at z ∼ 5 -7 are mistaken for the 3.3 PAH feature and the Brα line at low redshift. This degeneracy appears to be relatively common in our initial eazy-py fits; however, we note that the low-redshift fits also have very low inferred stellar masses of 10 4 -6 M ⊙ - in essence, a high-redshift star-forming galaxy is mistaken for a low-redshift globular cluster. Due to the low inferred mass of this lowredshift solution, this degeneracy is unlikely to be common using codes such as Prospector which enforce a lower limit on the allowed stellar mass. We note that these very low-mass solutions present an opportunity to implemented a linked mass-dust prior to further aid disentangling degenerate solutions (e.g., Nagaraj et al. 2022; Alsing et al. 2024).  \nThe MegaScience improvements on these individual SED fits adds up to a substantial improvement in overall photometric redshift precision. In Figure 5, we show the z phot performance from both our Prospector (top) and eazy (bottom) spectrophotometric redshift catalogs relative to a sample of ∼ 300 spectroscopic redshifts from the UNCOVER survey (S. Price et al. in prep.).  \nFigure 4. Example SEDs from UNCOVER (left) vs MegaScience (center), as well as eazy-py photometric redshift distributions (right). SEDs are shown in f ν (in units of 10 30 erg s -1 cm -2 Hz -1 ). Each object is spectroscopically confirmed by UNCOVER (top four objects; Price et al. in prep.) or ALT (bottom two objects; Naidu & Matthee et al. in prep). Objects are ordered by increasing redshift; each is chosen to demonstrate one 'failure mode' of SED-fitting that is solved by the increased photometric coverage of MegaScience. Each example is discussed in additional detail in Section 3.1.  \n<!-- image -->  \nFigure 5. Comparison of photometric redshifts with a sample of ∼ 300 high-quality spectroscopic redshifts from the UNCOVER survey (S. Price et al., in prep). The top row shows our Prospector redshifts (for details, see Wang et al. 2024c), and the bottom row shows eazy-py redshifts based on photometry from the 'SUPER' catalog (for details, see Weaver et al. 2024). The left column shows redshifts from the public DR2 UNCOVER data release, while the right column shows MegaScience redshifts from the DR3 public data products accompanying this paper. MegaScience shows nearly a factor of three reduction in σ NMAD , and close to a factor of two reduction in catastrophic outliers (as in Weaver et al. 2024, defined as objects with ∆ z/ (1+ z spec ) ≥ 0 . 15; shown by dotted red lines). Prospector and eazy-py show similar performance.  \n<!-- image -->  \nWe quantify this performance using the normalized median absolute deviation, σ NMAD , as defined in Brammer et al. (2008) and Weaver et al. (2024), and the catastrophic outlier fraction, defined as the fraction of objects with | z phot -z spec | / (1 + z spec ) ≥ 0 . 15. We find that MegaScience cuts the scatter σ NMAD by a factor of three , and nearly halves the rate of catastrophic outliers. While the set of spectroscopically-confirmed galaxies from UNCOVER is a small and biased subset of the full galaxy population, this improvement in scatter and catastrophic outlier fraction is suggestive of a dramatic improvement in our catalog's overall photometric redshift precision and accuracy. Future comparisons, especially with grism redshifts, will help verify the ∼ 1 . 6% scatter quoted in this paper.  \nIn sum: the increased spectral resolution of the MegaScience observations are able to break common degeneracies in photometric redshift fitting, as well as improve the accuracy and precision of recovered physical quantities including photometric redshifts by factors of 2-3. Combining this data with the large set of ancillary data available in Abell 2744 - MUSE, HST, ALMA, and now JWST/NIRCam-grism and NIRSpec - will further improve our inferences about the physical properties of these distant objects.", '3.2. Spatially-resolved science: mapping galaxy growth across cosmic history': "The MegaScience dataset not only provides high spectral resolution photometry: the high spatial resolution provided by NIRCam allows us to map the spatiallyresolved properties of galaxies across > 10 Gyr of cosmic history. Unlike grism data, where best practices often involve filtering out continuum emission to reduce contamination (e.g., Matthee et al. 2023a), medium-band observations are sensitive to both continuum and line emission. This means that our data enable accurate 2D maps of both emission lines and, via spatially-resolved SED modeling, stellar population properties such as stellar mass and age. Figure 6 shows a range of strong spectral features that are accessible with each NIRCam medium band across redshift, from the cluster itself through cosmic dawn. Our survey strategy ('ALL the medium bands!') means that at most redshifts we cover multiple spectral features simultaneously.  \nThe remainder of this section describes the spatially resolved science achievable at various redshift epochs and provides example medium-band images of galaxies in our survey. We note that all images shown in this section have been gridded onto a common pixel scale of 0 . '' 04, but images have not been PSF-matched.", '3.2.1. The Low-redshift Universe: Star Formation, Stripping, and Dust Emission in Abell 2744': 'As shown in Figure 6, our imaging directly probes the 3.3 µ m PAH feature, Paα , Paβ , and HeI (10830 ˚ A) + Paγ in the Abell 2744 cluster itself. While Paβ and HeI+Paγ are unlikely to have high enough equivalent widths to observe via medium-band excesses, given our depth in the F250W filter (Table 1) we expect to observe Paα for any regions within cluster galaxies with SFRs ≳ 0 . 5 M ⊙ yr -1 .  \nFigure 7 shows examples of four cluster galaxies; the leftmost column shows a broad-band image constructed using the F115W, F277W, and F444W filters from UNCOVER(Bezanson et al. 2022); the remaining four color images are constructed using all dozen MegaScience medium bands as listed in the column headings. Typical cluster galaxies have minimal ongoing specific star formation rates and resemble the top object: relatively uniform in color across all medium and broad bands. However, the cluster hosts several spectacular examples of galaxies experiencing large starbursts, potentially driven by ram-pressure stripping (for examples at similar z ∼ 0 . 3, see Vulcani et al. 2024). In these three galaxies, both Paα (F250M) and the broad 3.3 µ m PAH feature (F410M & F430M) are clearly visible by eye. This emission varies across the galaxy, indicating local regions of elevated star formation and PAH emission (as discussed in detail in Spilker et al. 2023). Our observations allow us to map star formation and dust in these active cluster galaxies, as well as search for localized star-forming regions in cluster galaxies that are globally quiescent.', '3.2.2. Cosmic Noon: 1 ≲ z ≲ 3': 'At cosmic noon, our imaging is able to constrain multiple spectral features (especially Balmer lines) simultaneously; at z ∼ 1 we directly probe Pa α , Pa β , and H α +[NII], while at z ∼ 3 we directly probe Pa β , H α +[NII], H β +[OIII], and the Balmer break. Careful modeling - likely informed by the wealth of spectroscopic observations in Abell 2744 - will be required to tease out the relative contributions of H β and [OIII] (e.g., Kaasinen et al. 2017).  \nUsing these multiple line diagnostics, MegaScience is capable of simultaneously mapping star formation and dust attenuation within galaxies at cosmic noon. By combining these maps with structural fits to the galaxies, we plan to study how the relationship between galaxy star formation histories and structure builds up over time. Due to the high spatial resolution of NIRCam - and aided by strong gravitational lensing in the cluster, see Furtak et al. (2023b) - it is also possible to  \nFigure 6. Spectral features probed by our primary medium-band observations as a function of redshift.  \n<!-- image -->  \nstudy the clumpiness of both star formation and dust in massive galaxies at the height of cosmic star formation history (e.g., Shen et al. 2024; Ji et al. 2023).  \nFigure 8 shows example images of seven galaxies at 1 ≲ z ≲ 3. The leftmost column shows a broad-band image constructed using UNCOVER data; the other four columns show all NIRCam medium bands. The structures of these galaxy are extremely complex: even in seemingly-well-settled disks such as ID 56988, MegaScience observations reveal clumps of H α lighting up in the spiral arms. In many objects, star-forming clumps can be seen lighting up in multiple different strong emission lines including H β +[OIII], HeI+Pa γ , [SIII], and Pa β .', '3.2.3. Cosmic Morning: 3 ≲ z ≲ 6 . 5': "At 3 ≲ z ≲ 6 . 5, our data begin to probe the Balmer break and the [OII] and MgII lines, in addition to the H β +[OIII] and H α +NII lines. At these redshifts, our data are able to map star formation and dust attenuation, identify AGN, and potentially identify and constrain some of the first quiescent galaxies (e.g., Carnall et al. 2023; Setton et al. 2024). Importantly, including a suite of (fewer) medium bands in the restframe optical has already been shown to work well to photometrically select quiescent galaxies during this early era, even in the absence of the restframe-J band data (only probed by MIRI at these wavelengths; Alberts et al. 2023).  \nFigure 9 shows example images in this redshift range, following the same format as Figure 8. These galaxies are morphologically complex, often with multiple  \nbright clumps. Significant line boosting is present in multiple emission lines, and in some cases emission line ratios show clear spatial variation even by eye. This highlights how the power of simultaneous coverage of H β +[OIII] and H α +NII (provided by our medium bands up to z < 6 . 5) at high spatial resolution can provide key insights into the drivers of these emission lines (e.g. star formation bursts vs ionization Endsley et al. 2023) while controlling for biases present in integrated measurements (Gim'enez-Arteaga et al. 2023).  \nOur data has near-complete coverage of the H α +NII line complex in this redshift range, with important crosssection at z > 3 with Lymanα probed by other facilities. Lymanα emission (rest-frame 1216 ˚ A) experiences resonant scattering with neutral Hydrogen, making it a probe of the state of cosmic reionization (e.g. Malhotra & Rhoads 2004). However, using Lymanα as the single probe of reionization is complicated by the strong dependence of Lymanα escape on internal galaxy properties (such as dust and young stellar populations) and observational effects (such as bias in slit losses due to spatial offsets between rest-UV and Lymanα ; Matthee et al. 2016; Hoag et al. 2019). These issues motivated panchromatic surveys of Lymanα emitters at z ≲ 3 that suggest that low metallicity stars with hard ionizing UV spectra are a key source of photons that can escape to ionize the IGM (e.g. Steidel et al. 2018).  \nWith MegaScience, we can directly observe H α instead of Lymanα into the late epoch of reionization, facilitating studies of the relative ionizing photon production efficiency of galaxies over a large range in mass and lu-  \nFigure 7. Examples of four galaxies in the Abell 2744 cluster itself; each object uses the same arcsinh color stretch. While most cluster galaxies have low ongoing SFRs and relatively uniform colors (like the top object), the cluster also hosts several spectacular morphologically-disturbed galaxies with strong Pa α and 3.3 µ m PAH emission. The white scale bar shows 0 . '' 5.  \n<!-- image -->  \nnosity. This can be coupled with the deep MUSE spectroscopy available in the field (Mahler et al. 2018) covering Lymanα out to z = 6 . 7. Our full suite of filters expands the redshift range probed compared to earlier medium band studies (e.g. Simmonds et al. 2023, 2024), enabling a more complete assessment of redshift evolution. The strength of this approach is the 'untargeted' nature of both sets of observations, allowing us to study H α in NIRCam for cases where Lymanα is not detected in the MUSE data, and vice-versa. This will result in an unbiased and complete picture that is impossible with data sets that can only focus on a single observable.", '3.2.4. Reionization: 6 . 5 ≲ z < 10': "At this redshift range H α +NII has redshifted out of the NIRCam wavelength coverage; however, our obser-  \nations still probe the Balmer break, [OII], [NeIII] and H β +[OIII] (among bluer rest-UV probes). Strong emission features are common, and not all galaxies are pointlike: many still show complex, clumpy, extended morphologies (see Section 4). While the inclusion of the two bluest broad-band filters F070W and F090W allows us to select high-redshift candidates with more confidence than in the UNCOVER data alone, the suite of 3-4 µ m medium-band filters drastically expands the astrophysical constraints on the ionizing radiation from Reionization-era sources (as suggested by mock observations; Roberts-Borsani et al. 2021). In particular, the simultaneous sampling of both [OIII]+H β and the Balmer break improve constraints that discriminate between timescales of star formation history with other parameters that influence the line equivalent widths such  \nFigure 8. Examples of galaxies at cosmic noon, 1 ≲ z ≲ 3. Color images are shown on a linear stretch; the white scale bar shows 0 . '' 5, which corresponds to ∼ 4 kpc at these redshifts.  \n<!-- image -->  \nFigure 9. Example of galaxies at cosmic morning, 3 ≲ z ≲ 6. Color images are shown on a linear stretch; the white scale bar shows 0 . '' 5.  \n<!-- image -->  \nas metallicity or ionizing escape fraction (e.g. Endsley et al. 2023; Laporte et al. 2023).", '3.2.5. Cosmic Dawn: z ≳ 10': 'Beyond z > 10, NIRCam only probes the restframe UV of galaxies, yet MegaScience will reveal important diagnostics for the formation of the first stars and galaxies. These include the Lyman break, MgII, Balmer break if present, [OII] (and neighboring NeIII, which JWST has revealed appears strong in early galaxies; e.g. Castellano et al. 2024; Roberts-Borsani et al. 2024), as well as exquisite sampling of the restframe UV continuum shape, β . Together, these spectral signatures can provide improved astrophysical insight into the ages and dust content of these earliest galaxies.  \nRobust identification of candidates, however, remains a challenge with broad band photometry, the workhorse of early JWST surveys in Cycle 1. The community has now gained enormous insight into the primary interloper population in the identification of galaxies at Cosmic Dawn. Namely, galaxies at 3 < z < 6 can mimic highredshift sources through some combination of balmer break, emission line boosting and dust (e.g. Naidu et al. 2022b; Zavala et al. 2023; Arrabal Haro et al. 2023). Better sampling of the dropout features improve the discriminating power between these scenarios, with sensitivity to the difference between a sharp Lyman-alpha break and the more gradual spectral shapes, and emission line boosting, of the primary lower-redshift interloper population (see e.g. Figure 2 in Eisenstein et al. 2023). In principal, MegaScience provides medium-band based selections (or, upper limits to abundances in case of non-detection) from z ≳ 11 (F140M-dropouts).  \nAlready, medium bands are enabling improved interloper rejection, providing constraints on abundances and physical properties up to z ∼ 14 (Robertson et al. 2023b; Willott et al. 2023; Bouwens et al. 2023). A critical improvement provided by MegaScience is the larger 3 × area of coverage compared to some early medium band surveys (e.g. JEMS, JOF), since early galaxies are extremely rare and can be highly clustered. MegaScience nearly doubles the coverage from wider-area, full-filter medium band data from CANUCS, which is already providing evidence of more rapid evolution in the UV luminosity function compared to broad band surveys (Willott et al. 2023).', '4. SCIENCE DEMONSTRATION: CLUMPY, OFFSET LINE EMISSION IN REIONIZATION-ERA GALAXIES': "In this section, we demonstrate how MegaScience imaging can be used to create maps of both line and continuum emission in galaxies, focusing on H β +[OIII]  \nline emission in three spectroscopically-confirmed galaxies at z ∼ 6 . 3. Spatially-extended ionized regions have begun to be observed at z ≳ 6 even in broad-band NIRCam imaging, pointing towards the the importance of ionized galaxy outflows at high redshfit (e.g., Fujimoto et al. 2022; Zhang et al. 2023). By more narrowly sampling the line-emitting region, medium band imaging is a powerful tool to accelerate these studies of emissionline morphologies in the era of reionization and cosmic dawn. Unlike the color images shown in Section 3.2, the images we show here are PSF-matched to the F444W filter using the same methodology described in Weaver et al. (2024).  \nWe select three example galaxies from a preliminary UNCOVER spectroscopic redshift catalog (S. Price et al., in prep) which have clear H β +[OIII] emission at ∼ 3 . 6 µ m, in the center of the F360M filter. We take the neighboring medium-band filters - F335M on the blue side, and F410M on the red side - to represent the continuum emission in these galaxies. The left column of Figure 10 shows the UNCOVER spectra (black) as well as the filter curves for the F335M (blue), F360M (green), and F410M (red) filters. The spectrum confirms that the continuum is relatively flat across the 3-4 µ m range, and the continuum filters are free of any major emission lines. The central column shows color images of each of the objects using these three medium-band filters; the top source resembles a point source, while the other two have complex and clumpy morphologies.  \nWe then use these three-color images to make simple maps of the continuum and line emission in these galaxies. We take the average of the blue and red continuum bands (F335M and F410M) as the continuum map, shown in the central column. The line map, shown in the rightmost column, is constructed by subtracting the continuum map from the F360M image, which contains both line and continuum emission from the galaxy. This method of deriving the spatial distribution of the continuum may be too simple for a randomly-selected galaxy in our sample - it could easily be biased by other lines contaminating the 'continuum' bands or a strong spectral slope. However, because the spectra in Fig 10 show that the continuum is relatively flat and line-free across these filters, this simple method of spatially distinguishing line and continuum emission is likely sufficient for these three sources.  \nA grey 'x' marks the center of the color images and continuum and line maps. While the top source, ID 10522, appears to be compact in all maps (with a size consistent with being an unresolved source), IDs 12053 and 31282 have clear differences in the continuum and line maps. In both cases, the color image shows three  \nFigure 10. Left: UNCOVER spectroscopy of three z ∼ 6 . 3 objects with bright [OIII] emission; blue and red curves indicate the medium-band filters used to represent the continuum, while the green curve shows the filter containing the strong line. Center: three-color PSF-convolved images of each source. Right: continuum and line maps for each source. While ID 10522 appears point-like in both continuum and line emission, both 12053 and 31282 exhibit complex morphologies, with three separate clumps at the same redshift, and clear spatial offsets between continuum and line emission.  \n<!-- image -->  \ndistinct clumps. In ID 12053, the central clump is significantly brighter in line emission while the two southern clumps are brighter in continuum; in ID 21382, the central clump is brightest in line emission while the western clump is brightest in continuum emission. We note that ID 12053 was detected in [CII] 158 µ m by the DUALz survey (Fujimoto et al. 2023), with a ∼ 300km/s offset between the [CII] 158 µ m and [OIII] 88 µ m emission. Fujimoto et al. (2023) discusses that this kinematic offset may be due to complex gas dynamics due to an ongoing merger, an interpretation that is further supported by our observation of different morphologies between the [OIII] line emission and ∼ 3 µ m continuum emission in this object.  \nThe high sensitivity and spatial resolution ( ∼ 0 . '' 04 -0 . '' 15) achieved by the NIRCam 1-5 µ m filters in this study have revealed the detailed structures of the stellar component and ionized gas in the rest-optical wavelengths, unveiling the complex and clumpy morphologies  \nof galaxies at z ∼ 6. A growing number of studies have identified similarly resolved sources and/or sources with very nearby companions at z ≳ 6 (e.g., Matthee et al. 2023a; Hainline et al. 2023; Yan et al. 2023; Kartaltepe et al. 2023; Witten et al. 2024). These findings are consistent with pre-JWST observations finding a trend of increasing clumpiness (e.g., Livermore et al. 2012, 2015; Shibuya et al. 2015) and merging fraction (e.g., Romano et al. 2021) towards higher redshifts. Recent JWST lensing observations have also revealed high clumpiness in galaxies at z ≳ 6 (e.g., Vanzella et al. 2023; Adamo et al. 2024; Mowla et al. 2024; Fujimoto et al. 2024). While mergers play a crucial role in mass assembly, the high clumpiness within individual galaxies suggests that active star formation may be driven at least in part by disk instability (e.g., Tadaki et al. 2018; Fujimoto et al. 2024). Our observations demonstrate the complex assembly histories of galaxies even at very early cosmic times.", '5. DISCUSSION & CONCLUSIONS': "In this paper, we describe the 'Medium Bands, Mega Science' survey, a survey of Abell 2744 using all medium- and broad-band JWST/NIRCam filters, and publicly release reduced images and a photometric catalog. With 20 bands of JWST photometry spanning the full 0.7-5 µ m range, our data are able to dramatically improve photometric redshift estimates and stellar population fitting. As shown in Figure 4, the MegaScience data is able to precisely constrain break strengths and stellar ages, mitigate confusion between the Balmer and Lyman breaks, and disentangle the contribution of line and continuum emission to broad-band photometry. On a survey-wide basis, comparison to high-quality prism redshifts from UNCOVER (S. Price et al. in prep.) shows that MegaScience photometric redshifts have a factor of ∼ 3 lower scatter and a factor of ∼ 2 fewer catastrophic outliers than the UNCOVER broad-band photometric redshifts. These factors of 2-3 improvement in photometric redshifts will propagate down to improvements in other stellar population properties such as stellar masses, star formation rates, and dust contents; a full exploration of the stellar population properties of galaxies in the survey will be presented in B. Wang et al. in prep.  \nIn addition to improving our understanding of the spatially-integrated properties of galaxies in our survey, the high spatial resolution of this data combined with lensing magnifications from the cluster itself enable us to spatially-resolve galaxy growth, mapping up to ∼ 500pc physical scales across more than ten gigayears of cosmic history. These observations add to an ongoing treasury of deep data across this field: Abell 2744 has now been observed using nearly every instrument and mode on JWST, except MIRI. These data will allow us to compare and calibrate NIRCam imaging, NIRCam grism, and NIRSpec prism and medium-resolution spectroscopy, and understand the formation and evolution of the earliest galaxies.  \nWe publicly provide full reduced images as well as our photometric catalog and both eazy-py and Prospector photometric redshift catalogs on our website, jwst-uncover.github.io/megascience as well as Zenodo, zenodo.org/doi/10.5281/zenodo.8199802. Future image and photometric catalog releases will be updated on our website. Our team will also publicly provide Prospector catalogs of stellar population properties (B. Wang et al. in prep.).  \nKAS thanks Shelly Meyett, the MegaScience Program Coordinator, for invaluable assistance designing WOPR 88967 and ensuring that this program was fully observed. This proposal was conceived of and developed at the International Space Science Institute (ISSI) in Bern, through ISSI International Team project #562. This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The raw 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 2 GO program #4111, and this project has gratefully made use of a large number of public JWST programs in the Abell 2744 field including JWST-GO-2641, JWST-ERS-1324, JWST-DD2756, JWST-GO-2883, JWST-GO-3538, and JWSTGO-3516. Support for program JWST-GO-4111 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 NAS5-26555.  \nThe Cosmic Dawn Center is funded by the Danish National Research Foundation (DNRF) under grant #140. PD acknowledges 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. This work has received funding from the Swiss State Secretariat for Education, Research and Innovation (SERI) under contract number MB22.00072, as well as from the Swiss National Science Foundation (SNSF) through project grant 200020 207349. 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. TBM was supported by a CIERA Fellowship. The BGU lensing group acknowledges support by Grant No. 2020750 from the United StatesIsrael Binational Science Foundation (BSF) and Grant No. 2109066 from the United States National Science Foundation (NSF); by the Ministry of Science & Technology, Israel; and by the Israel Science Foundation Grant No. 864/23. YF acknowledges support from JSPS KAKENHI Grant Number JP22K21349 and JP23K13149. RP and DM acknowledge funding from JWST-GO-02561.013 and JWST-GO-04111.035.  \nSoftware: astropy (Astropy Collaboration et al. 2013, 2018),grizli(Brammeretal.2014),Seaborn(Waskometal. 2017), WebbPSF(Perrinetal.2014)", 'REFERENCES': "Weaver, J. R., Cutler, S. E., Pan, R., et al. 2024, ApJS,  \n- 270, 7, doi: 10.3847/1538-4365/ad07e0 Whitaker, K. E., Labb'e, I., van Dokkum, P. G., et al. 2011, ApJ, 735, 86, doi: 10.1088/0004-637X/735/2/86 Whitaker, K. E., Ashas, M., Illingworth, G., et al. 2019, ApJS, 244, 16, doi: 10.3847/1538-4365/ab3853 Whitler, L., Stark, D. P., Endsley, R., et al. 2023, MNRAS, 519, 5859, doi: 10.1093/mnras/stad004 Wild, V., Taj Aldeen, L., Carnall, A., et al. 2020, MNRAS, 494, 529, doi: 10.1093/mnras/staa674 Williams, C. C., Alberts, S., Ji, Z., et al. 2023a, arXiv e-prints, arXiv:2311.07483, doi: 10.48550/arXiv.2311.07483 Williams, C. C., Tacchella, S., Maseda, M. V., et al. 2023b, ApJS, 268, 64, doi: 10.3847/1538-4365/acf130 Willott, C. J., Desprez, G., Asada, Y., et al. 2023, arXiv e-prints, arXiv:2311.12234, doi: 10.48550/arXiv.2311.12234 Withers, S., Muzzin, A., Ravindranath, S., et al. 2023, ApJL, 958, L14, doi: 10.3847/2041-8213/ad01c0 Witten, C., Laporte, N., Martin-Alvarez, S., et al. 2024, Nature Astronomy, doi: 10.1038/s41550-023-02179-3 Wolf, C., Meisenheimer, K., Rix, H. W., et al. 2003, A&A, 401, 73, doi: 10.1051/0004-6361:20021513 Wuyts, S., Forster Schreiber, N. M., Genzel, R., et al. 2012, ApJ, 753, 114, doi: 10.1088/0004-637X/753/2/114 Xiao, M., Oesch, P., Elbaz, D., et al. 2023, arXiv e-prints, arXiv:2309.02492, doi: 10.48550/arXiv.2309.02492 Yan, H., Ma, Z., Ling, C., Cheng, C., & Huang, J.-S. 2023, ApJL, 942, L9, doi: 10.3847/2041-8213/aca80c Zavala, J. A., Buat, V., Casey, C. M., et al. 2023, ApJL, 943, L9, doi: 10.3847/2041-8213/acacfe Zhang, Y., Ouchi, M., Nakajima, K., et al. 2023, arXiv e-prints, arXiv:2306.07940, doi: 10.48550/arXiv.2306.07940  \nZibetti, S., Charlot, S., & Rix, H.-W. 2009, MNRAS, 400, 1181, doi: 10.1111/j.1365-2966.2009.15528.x"}