We report the discovery of Cerberus an extremely red object detected with the MIRI Deep Imaging Survey MIDIS observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at signaltonoise ratio SN  6 with F1000W  27 mag and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey JADES fainter than the 30.030.5 mag 5 detection limits in individual bands as well as in the MIDIS F560W ultradeep data 29 mag 5. Analyzing the spectral energy distribution built with lowSN lt5 measurements in individual opticaltomidinfrared filters and higherSN 5 measurements in stacked NIRCam data we discuss the possible nature of this red NIRCamdark source using a battery of codes. We discard the possibility of Cerberus being a solar system body based on the lt0.016 proper motion in the 1 yr apart JADES and MIDIS observations. A substellar Galactic nature is deemed unlikely given that the Cerberuss relatively flat NIRCamtoNIRCam and very red NIRCamtoMIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers three possibilities 1 a z  0.4 galaxy with strong emission from polycyclic aromatic hydrocarbonsthe very low inferred stellar mass M SUBSUB  10SUP5SUP10SUP6SUP M SUBSUB makes this possibility highly improbable 2 a dusty galaxy at z  4 with an inferred stellar mass M SUBSUB  10SUP8SUP M SUBSUB and 3 a galaxy with observational properties similar to those of the reddest little red dots discovered around z  7 but Cerberus lying at z  15 with the restframe optical dominated by emission from a dusty torus or a dusty starburst.
['Galaxy formation', 'Galaxy evolution', 'High-redshift galaxies', 'Stellar populations', 'Broad band photometry', 'Galaxy ages', 'Active galactic nuclei', 'James Webb Space Telescope', '595', '594', '734', '1622', '184', '576', '16', '2291', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
A NIRCamdark Galaxy Detected with the MIRIF1000W Filter in the MIDISJADES Hubble Ultra Deep Field
{'A NIRCam-dark galaxy detected with the MIRI/F1000W filter in the MIDIS/JADES Hubble Ultra Deep Field': "Pablo G. P'erez-Gonz'alez , 1 Pierluigi Rinaldi , 2 Karina I. Caputi , 2 Javier ' Alvarez-M'arquez , 1 Marianna Annunziatella , 1 Danial Langeroodi , 3 Thibaud Moutard , 4 Leindert Boogaard , 5 Edoardo Iani , 2 Jens Melinder , 6 Luca Costantin , 1 Goran Ostlin , 6 Luis Colina , 1 Thomas R. Greve , 7, 8, 9 Gillian Wright , 10 Almudena Alonso-Herrero , 11 Arjan Bik , 6 Sarah E. I. Bosman , 12, 13 Alejandro Crespo G'omez , 1 Daniel Dicken , 10 Andreas Eckart , 14 Macarena Garc'ıa-Mar'ın , 15 Steven Gillman , 7, 8 Manuel Gudel , 16, 17 Thomas Henning , 5 Jens Hjorth , 3 Iris Jermann , 7, 8 ' Alvaro Labiano , 18 Romain A. Meyer , 19 Florian Pei β ker , 14 John P. Pye , 20 Thomas P. Ray , 21 Tuomo Tikkanen , 20 Fabian Walter , 5 and Paul P. van der Werf 22  \n1  \nCentro de Astrobiolog'ıa (CAB), CSIC-INTA, Ctra. de Ajalvir km 4, Torrej'on de Ardoz, E-28850, Madrid, Spain  \n2  \nKapteyn Astronomical Institute, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands  \n3  \nDARK, Niels Bohr Institute, University of Copenhagen, Jagtvej 155, 2200 Copenhagen, Denmark  \n- 4 European Space Agency (ESA), European Space Astronomy Centre (ESAC), Camino Bajo del Castillo s/n, 28692 Villanueva de la Ca˜nada, Madrid, Spain  \n5  \nMax Planck Institut fur Astronomie, Konigstuhl 17, D-69117, Heidelberg, Germany  \n- 6 Department of Astronomy, Stockholm University, Oscar Klein Centre, AlbaNova University Centre, 106 91 Stockholm, Sweden 7 Cosmic Dawn Center (DAWN), Denmark  \n8 DTU Space, Technical University of Denmark, Elektrovej, Building 328, 2800, Kgs. Lyngby, Denmark  \n9  \nDept. of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom  \n10  \nUK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK  \nCentro de Astrobiolog'ıa (CAB), CSIC-INTA, Camino Bajo del Castillo s/n, E-28692 Villanueva de la Ca˜nada, Madrid, Spain  \n12  \nInstitute for Theoretical Physics, Heidelberg University, Philosophenweg 12, D-69120, Heidelberg, Germany  \n13  \nMax-Planck-Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany  \nI.Physikalisches Institut der Universitat zu Koln, Zulpicher Str. 77, 50937 Koln, Germany  \n15 European Space Agency, Space Telescope Science Institute, Baltimore, Maryland, USA  \nDept. of Astrophysics, University of Vienna, Turkenschanzstr 17, A-1180 Vienna, Austria  \n14  \n16  \nETH Zurich, Institute for Particle Physics and Astrophysics, Wolfgang-Pauli-Str. 27, 8093 Zurich, Switzerland  \n17  \n- 18 Telespazio UK for the European Space Agency, ESAC, Camino Bajo del Castillo s/n, 28692 Villanueva de la Ca˜nada, Spain 19 Department of Astronomy, University of Geneva, Chemin Pegasi 51, 1290 Versoix, Switzerland\n- 20 School of Physics & Astronomy, Space Research Centre, Space Park Leicester, University of Leicester, 92 Corporation Road, Leicester LE4 5SP, UK\n- 21  \nDublin Institute for Advanced Studies, 31 Fitzwilliam Place, D02 XF86, Dublin, Ireland 22 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands", 'ABSTRACT': "We report the discovery of Cerberus , an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at S/N ∼ 6, with F1000W ∼ 27 mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag 5 σ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ( ∼ 29 mag, 5 σ ). Analyzing the spectral energy distribution built with lowS/N ( < 5) measurements in individual optical-to-midinfrared filters and higher S/N ( ≳ 5) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the < 0 . '' 016 proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-toNIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A z ∼ 0 . 4 galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, M ⋆ = 10 5 -6 M ⊙ , makes this possibility highly improbable. (2) A dusty galaxy at z ∼ 4 with an inferred stellar mass M ⋆ ∼ 10 8 M ⊙ .  \n11  \n- (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around z ∼ 7, but Cerberus lying at z ∼ 15, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst.  \nKeywords: Galaxy formation (595) - Galaxy evolution (594) - High-redshift galaxies (734) - Stellar populations (1622) - Broad band photometry (184) - Galaxy ages (576) - JWST (2291) - Active galactic nuclei (16)", '1. INTRODUCTION': "Building increasingly powerful telescopes operating at redder and redder wavelengths and (unexpectedly) discovering new types of galaxies at higher and higher redshifts is a classical industry now, established nearly 40 years ago. Indeed, deep optical surveys carried out in the late twentieth century using highly sensitive detectors were found to miss interesting galaxy populations of evolved and dust-rich galaxies at cosmological distances, both types being relatively bright at near-infrared wavelengths and easily detectable by new instruments, even though they were significantly less efficient in detecting photons and were affected by much larger backgrounds. The key for their discovery was large (red) colors (flux density ratios typically larger than an order of magnitude) between a near-infrared and an optical band, i.e., they were extremely red objects.  \nThe first Extremely Red Objects (EROs) were reported using new ground-based near-infrared telescopes built nearly 40 years ago. They were found to be members of cosmologically relevant galaxy populations, comprising a variety of galaxies ranging from high redshift quasars to very young and dust-rich galaxies, or even 'high' redshift passive galaxies (see, among many, Elston et al. 1988; McCarthy 2004; Graham & Dey 1996; Franx et al. 2003; Pozzetti & Mannucci 2000).  \nWith the development of very sensitive instruments onboard space telescopes such as Spitzer , the search was moved to longer wavelengths in the mid-infrared, at first compared with optical datasets (Wilson et al. 2004; Yan et al. 2004).  \nRemarkably, when the Hubble Space Telescope gathered extremely deep data taken in the near-infrared up to 1.6 µ m, with surveys such as CANDELS or the WFC3 Hubble Ultra Deep Field (Koekemoer et al. 2011; Grogin et al. 2011; Koekemoer et al. 2013), Spitzer /IRAC data were still shown to detect unique distinctively red objects which were only visible at wavelengths longer than 3 µ m, or at least were extremely faint in the near-infrared and optical spectral ranges. These socalled HST-dark galaxies were studied for around a decade, with significant uncertainties about their nature due to the strong sensitivity limitations of our  \nmulti-wavelength telescopes, even those operating in very different observational windows such as the (sub)millimeter range (see, e.g., Simpson et al. 2014; Franco et al. 2018; Alcalde Pampliega et al. 2019; Wang et al. 2019). JWST has started, from the earliest phases of the mission, to shed light on those HST-dark galaxies, i.e., mid-infrared bright near-infrared/optically faint sources which now constitute the bulk of the JWST galaxy exploration (P'erez-Gonz'alez et al. 2023a; Nelson et al. 2023; Barrufet et al. 2023; Williams et al. 2023a; Gillman et al. 2023).  \nThe ERO astrophysical industry is intimately tied to the discovery of submillimeter galaxies (SMGs, Hughes et al. 1998). Indeed, many EROs are dusty starbursts, which appear in mid- and far-infrared surveys in different flavors and with different names, starting from SMGs (see, e.g., Blain et al. 2002; Chapman et al. 2005) but also (hot) dust-obscured galaxies (DOGs, Dey et al. 2008; Narayanan et al. 2010; hotDOGs, Tsai et al. 2015), extremely red quasars (ERQ, Hamann et al. 2017), or, more descriptive from an observational point of view, (ultra-)luminous infrared galaxies (U/LIRGs, Sanders & Mirabel 1996; Elbaz et al. 2002; P'erez-Gonz'alez et al. 2005; Le Floc'h et al. 2005; Caputi et al. 2007; Magnelli et al. 2009; Madau & Dickinson 2014).  \nLast but not least, EROs, especially if selected at the faintest magnitudes and with appropriate bands for the color, overlap with very high-redshift galaxy samples selected using the Lyman break technique (Steidel et al. 1996) at higher and higher redshifts (e.g., Bouwens et al. 2011).  \nRemarkably, with the advent of JWST, both industries, searching for red objects and for those with higher and higher redshifts, have been merged at a profound level, with the possibility of finding double-break (Lyman and Balmer) red galaxies (Labb'e et al. 2023a), now known as little red dots (Furtak et al. 2023; Ubler et al. 2023; Matthee et al. 2023; Barro et al. 2023; Greene et al. 2023; Williams et al. 2023a; Killi et al. 2023; P'erezGonz'alez et al. 2024; Kokorev et al. 2024), or high redshift galaxies with very high equivalent width emission lines, first discovered with Spitzer (Smit et al. 2014), now appearing in many JWST works (P'erez-Gonz'alez et al.  \n2023a; Rinaldi et al. 2023). Interestingly, some of these sources have been confirmed to present broad emission line components possibly linked to an active galactic nucleus (AGN; Kocevski et al. 2023; Greene et al. 2023).  \nThe interesting subtopics of finding and characterizing dusty galaxies on the one hand, and high redshift galaxy candidates on the other, has been demonstrated to be very relevant for JWST surveys, which can easily misidentify one type of galaxy with the other (see, e.g., Zavala et al. 2023; Naidu et al. 2022).  \nThe MIRI instrument onboard the JWST (Rieke et al. 2015; Wright et al. 2023) now offers the possibility to search for red galaxies with longer color baselines, beyond the wavelengths first probed by Spitzer /IRAC and now easily accessible with impressive depths ( ∼ 30 mag and beyond) using JWST/NIRCam. Even though MIRI cannot reach NIRCam depths, very deep imaging at 5.6 - 12.8 µ m reaching magnitudes around 27 - 29 is already available and opens a whole new window for the search and characterization of high redshift galaxies (Papovich et al. 2023; Barro et al. 2023; Rinaldi et al. 2023; P'erez-Gonz'alez et al. 2024).  \nIn this paper, we report the discovery of a MIRI-bright red object appearing in the MIRI Deep Imaging Survey (MIDIS) F1000W imaging of the Hubble Ultra Deep Field ( Ostlin et al., 2024, in prep.). The source is extremely faint at NIRCam long wavelengths, i.e., it is a NIRCam-dark object lying at the edge of the detection limit ( ∼ 30 . 5 mag) of one of the deepest JWST surveys available to date, JADES (Eisenstein et al. 2023a), having not been catalogued before (Rieke et al. 2023; Eisenstein et al. 2023b). Consequently, our object presents an extremely red color between 4.4 and 10 µ m, and, in fact, the source is also red with respect to the bluest MIRI filter, as revealed by the MIDIS F560W imaging.  \nThis letter is organized as follows. Section 2 presents the deep MIRI data used in this letter, as well as ancillary observations taken with NIRCam, ALMA, and VLT. Section 3 describes our method to detect NIRCamdark galaxies, and the case for a positive detection in the Hubble Ultra Deep Field, the source we call Cerberus . Section 4 presents the spectral energy distribution of Cerberus and discusses the several possible interpretations, analyzing the implied physical properties for each one of them. Finally, Section 5 presents a summary of our findings.  \nThroughout the paper, we assume a flat cosmology with Ω M = 0 . 3 , Ω Λ = 0 . 7, and a Hubble constant H 0 = 70kms -1 Mpc -1 . We use AB magnitudes (Oke & Gunn 1983). All stellar mass and SFR estimations assume a universal Chabrier (2003) initial mass function (IMF), unless stated otherwise.", '2.1. MIRI data': "MIRI data in the F1000W filter were taken as part of the MIRI Deep Imaging Survey (MIDIS, Ostlin et al. 2024, in prep.; PID 1283) in December 2023. The observation consisted of 11 exposures, each with 100 groups, FASTR1 readout, and 10 integrations, for a total on source exposure time of 30800 s, centered on the Hubble Ultra Deep Field (HUDF). The dithering pattern was set to large-size cycling, with the 11 exposures taken in different positions on the sky separated by up to 10 '' .  \nThe reduction of the MIRI data followed the methodology described in P'erez-Gonz'alez et al. (2024), based on the official jwst pipeline (v1.13.4, jwst 1202.pmap) and offline bespoke procedures. Briefly, a super-background strategy is used to build background maps for each single image using all the other exposures (since they were taken during the same campaign), which results in a very homogeneous background in terms of level and noise. Quantitatively for the data used in this letter, the standard deviation of the background pixels for each cal.fits file obtained with the standard official pipeline is reduced by a factor of ∼ 3 . 5 by using our superbackground method, and the depth in the final mosaic improves by 0.8 mag. We refer the reader to Appendix A in P'erez-Gonz'alez et al. (2024) for more details. Known sources are masked to avoid biasing the determination of the very local background in the super-background image. Our F1000W final mosaic, reduced with a 60 milliarcsec pixels, presents an average 5 σ depth of 26.8 mag, measured in a 0 . '' 4 radius circular aperture, and taking into account noise correlation (linked to the drizzling method used in the mosaicking) and after applying the aperture correction (a 1.3 factor).  \nWe present the MIDIS F1000W data in Figure 1 with an RGB image built in combination with some NIRCam long wavelength data.", '2.2. Ancillary data': 'In this subsection, we describe other datasets in the HUDF that we have used to complement the MIRI data.', '2.2.1. Other MIRI data': "For this letter, we combined the recently acquired MIDIS F1000W data with the F560W ultra-deep observations also carried out by MIDIS in December 2022. These data will be described in Ostlin et al. (2024, in prep.) and consist of ∼ 40 hours on source taken in the HUDF, which allowed us to reach 28.7 mag 5 σ for point-like sources measured in an r=0 . '' 23 circular aper-tu  \nFigure 1. In the center of the figure, we show an RGB composition of the MIDIS F1000W field. The color images have been built with JADES data in 2 NIRCam filters, F277W and F356W, and the MIDIS MIRI F1000W filter (all images convolved to the same PSF of MIRI/F1000W). In the background, we show the HUDF JADES data in grayscale. We provide a series of zoomed-in RGB frames that lead to the NIRCam-dark MIRI source, Cerberus , whose image is located at the bottom right corner. The NIRCam-dark MIRI source Cerberus is highlighted with a square in that RGB frame. Additionally, we provide postage stamp images (2 '' × 2 '' ) of Cerberus in MIRI/F1000W, MIRI/F560W, and NIRCam/F444W (left column), showcasing that Cerberus (marked with an r=0 . '' 25 green circle) is clearly detected only at 10 µ m.  \n<!-- image -->  \nradius being chosen to ensure a ∼ 70% encircled energy).", '2.2.2. NIRCam data': "The MIRI/F1000W data were also complemented with NIRCam imaging taken by the JWST Advanced Deep Extragalactic Survey, JADES (Eisenstein et al. 2023a), Data Release 2 (Eisenstein et al. 2023b), which includes observations from the JWST Extragalactic Medium-band Survey (JEMS, Williams et al. 2023b) and the First Reionization Epoch Spectroscopically Complete Observations (FRESCO, Oesch et al. 2023). This dataset provides a total of 14 bands from 0.9 to 4.8 µ m (6 at short-wavelength, SW, and 8 at longwavelength, LW), with 5 σ depths ranging from 30.5 to 30.9 mag (measured in a 0 . '' 2 radius circular apertures). We remark that JADES is one of the deepest NIRCam surveys on the sky, only matched in depth (in some bands) by the MIDIS NIRCam-parallel project (P'erezGonz'alez et al. 2023b) and The Next Generation Deep Extragalactic Exploratory Public Near-Infrared Slitless  \nSurvey, NGDEEP (Bagley et al. 2023), which means that our search for NIRCam-dark galaxies is really probing the current depth limitations of JWST extragalactic surveys.", '2.2.3. ALMA data': "The Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) is a Cycle 4 Large Program over a 4.6 arcmin 2 scan at 1.2 mm (Decarli et al. 2020; Gonz'alezL'opez et al. 2020) and 3.0 mm (Decarli et al. 2019; Gonz'alez-L'opez et al. 2019). The ultra-deep 1.2 mm data reaches an rms sensitivity of 9.3 µ Jy/beam, with beam dimensions of 1 . '' 5 × 1 . '' 1. The 3.0 mm data reaches 1.4 µ Jy/beam, with a 1 . '' 8 × 1 . '' 5 beam.", '2.2.4. MUSE data': "The HUDF has been extensively studied in the last decade with the Multi Unit Spectroscopic Explorer (MUSE, Bacon et al. 2010) mounted on the Very Large Telescope (VLT) as part of the MOSAIC and UDF-10 fields (GTO programs 094.A-0289(B), 095.A-0010(A),  \n096.A-0045(A) and 096.A-0045(B), PI: R. Bacon) and the most recent MXDF observations (GTO Large Program 1101.A-0127, PI: R. Bacon). The MUSE data has a spectral wavelength range between 4700 - 9300 ˚ A and a spectral resolving power (R) that ranges from 1770 (4800 ˚ A) to 3590 (9300 ˚ A). The HUDF dataset was taken during 140 hours on-source with the VLT Adaptive Optics Facility (AOF) and its GALACSI adaptive optics module, as detailed in Kolb et al. (2016) and Madec et al. (2018), achieving a spatial sampling of 0 . '' 2 × 0 . '' 2, with strong variations of up to a factor of × 4. More details about these programs can be found in Bacon et al. (2017, 2023).", '3. SELECTION AND VALIDATION OF NIRCAM-DARK SOURCES': 'This letter presents a source that was discovered on the MIDIS F1000W image and found not to be included in existing NIRCam published catalogs (Rieke et al. 2023; Eisenstein et al. 2023b). In this section, we describe the selection procedure for such types of NIRCamdark sources.', '3.1. Selection of NIRCam-dark sources': "In order to select NIRCam-dark sources, we executed the sextractor software (Bertin & Arnouts 1996) directly using the MIRI/F1000W as the detection image. We followed the hot-mode extraction described in Galametz et al. (2013), which is particularly effective for detecting extremely faint sources.  \nOnce we had a hot MIRI/F1000W catalog, we anti cross-matched it (using a 0 . '' 3 search radius) with the official JADES DR2 catalog published by Eisenstein et al. (2023b, which superseded the DR1 catalog), constructed by selecting sources in a LW stack. The goal of the anti cross-match was to isolate sources uniquely detected by MIRI/F1000W. The result was a list of 17 potential NIRCam-dark sources selected at 10 µ m.  \nWe meticulously examined each of the potential NIRCam-dark sources in the final mosaic to rule out the presence of an unflagged cosmic ray or any kind of hot and/or bad pixels that could have caused a spurious detection.  \nVisual inspection of the individual exposures revealed that for 16 of the 17 candidates, there was indeed a hot pixel in one of the 11 exposures which had survived the flagging and stacking process performed to produce the final mosaic. We automated and generalized the source validation methodology by producing stacked mosaics with a limited number of exposures, excluding 1, 2, or 3 of the full set randomly. This procedure, obviously, affects the depth. Explaining the latter case in more detail as an example, we started from our collection of 11  \n* cal.fits files, and we created the 165 unique mosaics which can be obtained by excluding sets of 3 * cal.fits files for each partial mosaic. After building the 165 realizations, we inspected the 17 potential NIRCam-dark sources, concluding that only one source consistently appeared in all 165 different partial mosaics. The remaining candidates were found to vanish in some of the realizations, in fact even in the mosaics built excluding one single exposure, indicating their most-likely spurious nature. To further validate these findings, we made an in-depth analysis of each * cal.fits file, confirming the presence of unflagged hot and/or bad pixels in some of them that indeed led to false detections in a number of the realizations.  \nOne source of the 17 pre-selected candidates did not vanish in any of the partial mosaics, even when removing 3 random images. For this source, we measured photometry in a small 0 . '' 25 radius aperture (which maximized S/N on the partial mosaics), as well as the sky noise in a square 3 '' × 3 '' region around, on all the mosaics produced by excluding 1, 2, and 3 exposures. We found consistent photometry and a degradation of the S/N consistent with the differences in exposure time. For example, for the skip-3 mosaics, we found an average magnitude 27.05 ± 0.19 mag, consistent with our measurement in the full mosaic (see Section 4).  \nHowever, we found 2 contiguous exposures, numbers 7 and 8, which, if removed, produced a dimming of the source to 27.4 mag, still within 2 σ of the measurement using the full dataset. Each one of these exposures has an integration time of 2800 s, with the start of the first exposure and the end of the second separated by ∼ 5700 s. We examined these exposures in detail and found no signs of a cosmic ray hit, but a slightly enhanced signal in an extended region. This could be due to the general properties of the detector (which consistently showed a dimmer region in the super-background frames for all images). It could also be due to a cosmic ray shower, which are known to affect MIRI imaging observations (Dicken et al. 2024), with an incidence of several per exposure. Cosmic ray showers are extended persistence artifacts from cosmic ray hits and typically affect the detector for several minutes (Wright et al. 2023). However, a high flux cosmic ray hit could cause a shower lasting long enough to enhance the signal seen in these exposures. We also note that the enhancement in signal is stronger in exposure 7 compared to 8, which is also consistent with persistence (see ' Alvarez-M'arquez et al. 2023). Given that the two relatively (and slightly) brighter images were contiguous, we do not rule out a transient behavior of the source. Based on the perseverance of the signal from the source in all partial stacks  \nFigure 2 shows a stacked NIRCam image (including all filters at λ > 3 µ m) convolved to the resolution of the F1000W data (FWHM=0 . '' 328). Each calibrated individual NIRCam filter was rebinned to a 0 . '' 06 pixel scale and PSF-matched to F1000W using WebbPSF (v. 1.2.1 Perrin et al. 2014) models and the method described in Melinder et al. (2023), and then median-stacked. We measured the centroid of the source for both JWST instruments (using photutils with a circular footprint of 3 pixels radius), finding an agreement within 0.26 pixels, i.e., 0 . '' 016. Remarkably, the source is slightly extended to the NW in both images. The morphology of the source is similar for the NIRCam stack and the MIRI F1000W image, with a tadpole-like shape extending to the upper right, confirming a detection of the source with NIRCam stacked data (see next section). We also  \n<!-- image -->  \nFigure 2. Stacked NIRCam (left panel) and MIRI/F1000W (right panel) images of the region around Cerberus . The stacked NIRCam image was produced from individually PSF matched LW images. Note that the spurious pixels in the F460M image (see text for details) are removed in the median stacking procedure. The angular size of the images and zoomed inset is 3 '' × 3 '' and 0 . '' 6 × 0 . '' 6, respectively. The red (blue) circle shows the location of the detected source in F1000W and has a radius of 0 . '' 2. The green filled circle shows the resolution of the F1000W image (with a diameter of 0 . '' 328). In the zoomed insets, we also show the centroid positions with the blue and red crosses, corresponding to the F1000W image and the NIRCam stack, respectively.  \n<!-- image -->  \nof the MIRI F1000W dataset, we concluded that it was not spurious.  \nThe NIRCam-dark source we identified with this methodology was dubbed Cerberus , to acknowledge that it is a weird beast lurking in the darkest corners of the JWST surveys and with different interpretations/heads for its nature, as will be shown in Section 4. Its coordinates (J2000) are α = 03:32:38.09 (hours), δ = -27:47:11.81 (degrees). NIRCam and MIRI images of Cerberus are presented in Figure 1.", '3.2. Multi-wavelength validation and observational properties of the NIRCam-dark object Cerberus': 'In order to further validate and characterize Cerberus , we searched for detections by forcing photometric measurements in the datasets described in Section 2, more specifically in the ultra-deep observations (among or even the deepest on the sky) provided by JWST/NIRCam, ALMA, and VLT/MUSE.', '3.2.1. JWST/NIRCam measurements': 'Directly linked to our selection criteria, the NIRCamdark Cerberus object was not present in any of the catalogs of the JADES NIRCam survey (Rieke et al. 2023; Eisenstein et al. 2023b). We then proceeded to force photometric measurements in all bands (described in Section 4). Forced photometric measurements in all  \nNIRCam individual bands do not provide any measurement signal-to-noise ratio S/N > 5.  \nWe then stacked NIRCam images to search for faint emission. We constructed two different sets of stacks. One set consisted in cumulative stacks constructed by averaging NIRCam bands starting from the reddest and from the bluest. Another set consisted of disjoint stacks (to avoid noise spikes in one band affecting the measurements in cumulative stacks), combining 2, 3, 4 and 5 bands, always starting from the reddest and not repeating them in different stacks of the same set.  \nnotice a qualitative resemblance between the structure of Cerberus and the z = 10 . 6 galaxy GN-z11 (Oesch et al. 2016; Tacchella et al. 2023), with a concentrated flux distribution identified with the core of the galaxy and an extended haze component.', '3.2.2. MUSE measurements': "We extracted the MUSE 1D spectrum at the coordinates of our target within circular apertures of different sizes (0 . '' 4 and 0 . '' 8 radius) but found no significant spectral feature at the position of the NIRCam-dark candidate. This is also confirmed by the publicly available catalog of detected sources by Bacon et al. (2023).", '3.2.3. ALMA measurements': 'The source is not detected at 1.2 mm and 3.0 mm. This finding implies 5 σ upper limits on the flux density at 1.2 mm and 3.0 mm of 47 µ Jy and 7 µ Jy respectively, assuming it is effectively an unresolved point-like source.', '4.1. Multi-wavelength Photometry': "A spectral energy distribution (SED) for Cerberus was built by measuring photometry in all available NIRCam (including both individual bands and stacks) and MIRI data, forcing measurements at the position determined for Cerberus with the MIRI data. Given the faint and small nature of our source, we measured forced photometry in several circular apertures with radii ranging from 0.1 '' to 0.6 '' . We measured the background in a small square region around Cerberus (3 '' on the side) and calculated the noise using non-adjacent pixels (5 pixels apart) to take into account noise correlation, as explained in P'erez-Gonz'alez et al. (2023b), obtaining similar depth results (within 0.2 mag) to those reported in Eisenstein et al. (2023a) for the NIRCam bands. Photometric measurements are provided in Table 1.  \nNone of the measurements in individual NIRCam bands have a > 5 σ significance, which we use as our limit to define a detection in each filter.  \nGiven our magnitude measurement in the F1000W band, 27 . 13 ± 0 . 17 mag, obtained with photometry in an aperture of radius 0 . '' 3, which maximizes the S/N , and the S/N ∼ 3 measurements in individual bands and S/N > 5 in stacked data for the reddest LW NIRCam channels, we infer a very red color, F444W -F1000W = 3 . 8 ± 0.4 mag. Using the reddest disjoint stack with an actual detection (i.e., S/N > 5), we obtain a color [4 . 497] -[10 . 0] = 2 . 8 ± 0 . 3 mag.  \nThese measurements indicate that Cerberus is not robustly detected in any individual NIRCam filter, but  \nconsistent with being slightly below the 5 σ detections limits in some bands, quoted to be around 30-31 mag for the JADES DR2 data (Eisenstein et al. 2023a). We remark that the F460M filter, which provides the brightest magnitude, presents a group of 4 pixels to the North of the position of our galaxy, which we confirmed to be an artifact of the reduction, and it is not seen in any other band (even the ones with deeper limiting magnitudes). We excluded this filter from the rest of analysis.  \nWe measured photometry in the stacked data in the same apertures used for the individual bands, assigning average aperture corrections for the filters included in each stack. The 5 LW stacks with the highest number of bands all provide S/N > 5 fluxes (the average being S/N ∼ 6), ranging from 29.9 to 30.8 mag. This means that the (low-S/N) SED is characterized by a red color extending through all the S/N > 5 stacks from 2 to 4 µ m and beyond, up to 29 . 9 ± 0.2 mag for the first (reddest) stack (including F480M+F444W), which reaches ∼ 5 µ m. For the SW stacks, no significant signal is recovered in any of them. We measure a 5 σ upper limit of ∼ 31 mag for the stack adding all SW data. The comparison between SW and LW stacks implies a ≳ 0.6 mag color jump from the ∼ 1 -2 to the ∼ 3 -5 µ m ranges.  \nThe SED of Cerberus is, in fact, very similar to the little red dots (LRDs) recently presented and characterized in a number of papers (e.g., Labb'e et al. 2023a,b, Matthee et al. 2023, Barro et al. 2023, Williams et al. 2023a, P'erez-Gonz'alez et al. 2024, Kokorev et al. 2024), but shifted to longer wavelengths. Based on detections of Cerberus in stacked images, we can infer that the F277W -F444W > 1 mag general selection criterion of LRDs, and, more specifically, the 3 mag color of the example source 203749 in P'erezGonz'alez et al. (2024), compares well, if we multiply the observed wavelengths by a factor of ∼ 2, with the F444W -F1000W ∼ 3 . 5 mag, or, using the F560W 5 σ upper limit, F560W -F1000W ≳ 2 mag. The flat F150W -F200W color distinctive of LRDs (the largest value allowed by the selection is typically 0.5 mag), i.e., a change in slope at around 0.3-0.4 µ m in the rest-frame for the typical redshift of LRDs 5 < z < 9, compares well with the color between 3 and 4 µ mfor Cerberus , just multiplying again the observed wavelengths by a factor of ∼ 2. We remark that this comparison is based on S/N > 5 measurements in stacked data. That rest-frame UV color is indeed consistent with the reddest LRDs in P'erez-Gonz'alez et al. (2024) or the LRD of which a NIRSpec spectrum was presented in Killi et al. (2023). The analogy can also be established in terms of the apparent magnitude: the typical magnitude of LRDs in the blue part of their SED lies around 27-28 mag, Cerberus being  \n≳ 10 times dimmer in NIRCam LW bands. We note that the difference in luminosity distance between z ∼ 7, the median redshift of LRDs (P'erez-Gonz'alez et al. 2024) and z ∼ 14 (a factor of 2 larger, see also Section 4.2.1) will translate to a flux ratio between bands covering the same rest-frame wavelength around an order of magnitude for the same type of object.  \nFigure 3. Photometric redshift probability distribution functions of Cerberus . All distributions are normalized to unity at the peak. We show the main solutions and also secondary solutions at z < 10 and z > 10. Different codes were used in the SED analysis, described in Appendix A.  \n<!-- image -->", '4.2. Analysis of the SED of Cerberus': "In the following subsections, we discuss possible interpretations of the SED of Cerberus . The NIRCam and MIRI datasets were taken more than one year apart. However, given the spatial coincidence of the object in the NIRCam and MIRI observations (centroids separated by 0 . '' 016), we proceed by assuming that they are the same object and not a chance alignment of two different objects.", '4.2.1. Extragalactic origin': 'The SED of Cerberus was fitted to estimate a photometric redshift and obtain its physical properties in a 2-step procedure. We used several codes for this purpose, separating the two tasks to analyze facets of each method and understand the (possibly systematic) uncertainties in detail. We used the fluxes described above as well as the ALMA upper limits quoted in Section 2, shown to be important for high redshift sources (even for determining their redshift, see Meyer et al. 2024), and particularly for an ERO such as Cerberus . We provide details about the fitting methods in Appendix A.  \nThe results of our SED modeling analysis are given in Figures 3 (photometric redshift probability distributions) and 4 (SED models).  \nWe remark that the SED of Cerberus is very noisy, given that it is an ultra-faint source at the very depth edge of the JWST surveys. Therefore, the physical properties of the source are highly uncertain, starting from its redshift (see several peaks in Figure 3).  \nHowever, leaving the redshift apart, there are two physical properties for which our results are similar among all SED modelling codes, mainly because they are linked to the 2 main observational characteristics of Cerberus , i.e., its faintness (translating to relatively low luminosity) and red color. The very red color of our NIRCam-dark object can only be reproduced with significant amounts of dust, either linked to PAHs if the redshift is low, z ∼ 0 . 4, or warm/hot dust if the redshift is around z ∼ 4 or z ∼ 15. For the lowest redshift, all codes provide very low stellar masses, around ∼ 10 5 -6 M ⊙ . The very small stellar mass, relatively high dust attenuation, and strong PAH emission is an uncommon combination, at odds with the massmetallicity relationship (e.g., Tremonti et al. 2004) and the weakening of PAH features at low metallicities (e.g., Engelbracht et al. 2005). For the higher redshift solutions, z ∼ 4 and z ∼ 15, all codes provide masses around ∼ 10 7 -8 M ⊙ , and they all need large warm dust contents, which can be explained with a very compact starburst or an AGN torus (similar to those detected at high redshift by JWST, see Maiolino et al. 2023, as well as Pacucci et al. 2023; Pacucci & Loeb 2024, and also Hamann et al. 2017).', '4.2.2. Solar System origin': "Abright mid-infrared object such as Cerberus matches the expected properties of Solar System bodies, e.g., asteroids. Given the long observational MIRI/F1000W campaign, we would expect to see proper motion, which is not the case, unless the asteroid is distant and/or near aphelion. However, the fact that the NIRCam data detect Cerberus (in stacked images) and the observations from MIDIS and JADES were taken more that one year apart indicates that Cerberus cannot be a Solar System body (assuming that the NIRCam and MIRI sources are the same object and not a fortuitous superposition). In Section 3.2, we also measure the centroid shift of the object between the MIRI and NIRCam data to 0 . '' 016 which is consistent with a zero offset. This makes the source even more unlikely to be a Solar System body.  \n4.2.3. Galactic origin  \nAs shown by Langeroodi & Hjorth (2023a), Burgasser et al. (2024), and Hainline et al. (2023) for LRDs, red  \nFigure 4. Detailed analysis of the spectral energy distribution of Cerberus for the main redshift solutions. The symbols are the same as described in Figure 3. We provide best-fitting models provided by different codes that include both stellar and AGN emission. The top-left panel shows the best-fitting brown dwarf template from the ATMO2020++ models, for a Y-type dwarf with an effective temperature of ∼ 300K at a distance of ∼ 300 pc. The rest of the panels show results for z = 0 . 4 (top-right), z = 3 . 5 (bottom-left), and z = 15 (bottom-right).  \n<!-- image -->  \nsources detected by JWST in cosmological fields such as GOODS-S or A2744 can be brown dwarfs. Following the method in Langeroodi & Hjorth (2023a), we fit the photometry of Cerberus with brown dwarf atmosphere templates. We include the Sonora cloud-free models (Karalidi et al. 2021; Marley et al. 2021) as well as the ATMO2020++ cloud-free adiabat-adjusted T and Y dwarf models (Leggett et al. 2021; Meisner et al. 2023, see also Su'arez & Metchev 2022). In particular the ATMO2020++ models were included to ensure coverage of extremely cool brown dwarfs with 250 K < T eff < 500K. Compared to warmer brown dwarfs, these models more closely resemble the SED shape of Cerberus , which  \nis relatively faint at near-infrared wavelengths while getting significantly brighter at 10 µ m.  \nThe best-fit brown dwarf template is shown in the top-left panel of Figure 4. We infer a best-fit effective temperature ( T eff [K]) of 300 +100 -25 , surface gravity (log g [cm s -2 ]) of 3 . 0 +1 . 0 -0 . 5 , and distance of 332 ± 67 pc. Its effective temperature would most likely classify Cerberus as a Y-type dwarf (see Figure 12 in Langeroodi & Hjorth 2023a). Compared to the UltracoolSheet compilation of known brown dwarfs with parallax distances (Best et al. 2020), its inferred distance ( ∼ 300pc) would make Cerberus one of the farthest (top 20) brown dwarfs discovered to date (see Figure 11 in Langeroodi & Hjorth  \n2023a). This distance at the high Galactic latitude of Cerberus would place it at the edge of the thin disk (see Figure 6 in Hainline et al. 2023). Recently, three T-type dwarfs were discovered and spectroscopically confirmed with NIRSpec prism (Langeroodi & Hjorth 2023a; Burgasser et al. 2023) at estimated distances between 0.7 and 4.8 kpc, all most likely outside the thin disk; among them A2744-BD3 is a late T-type at 755 pc.  \nOur best-fit to a brown dwarf spectral template, shown in the top-left panel of Figure 4, nicely fit the F444W -F1000W color and is consistent with the nondetection in the deep MIDIS F560W data. However, that brown dwarf model is very red in both the NIRCam SW-to-LW and 3-to-4 µ m colors, which is not favored by our measurements (in stacked data), pointing to a flatter SED.  \nIn addition, both the MIRI and stacked NIRCam images (see Section 3.2) shows the source to be slightly extended and asymmetrical as compared to the PSF. The two instruments, NIRCam and MIRI, indicate a morphology for the source that makes it unlikely that it is a point source, thus disfavoring the single brown dwarf scenario.  \nThe slight extension of order 0 . '' 1 alone would however not argue against a brown dwarf binary . Even at only 100 pc distance, a semi-major axis of 0 . '' 1 for a binary of two identical brown dwarfs/planets with masses of (0 . 0003 -0 . 07) M ⊙ each would result in orbit periods of order 100-1000 years, i.e., the elongated image would not vary between the NIRCam and the MIRI observations.", '5. SUMMARY AND CONCLUSIONS': "We report the discovery of a NIRCam-dark source identified in the MIRI Deep Imaging Survey (MIDIS) of the Hubble Ultra Deep Field carried out with the F1000W filter. We call this source Cerberus . The source has a magnitude around 27 mag (6 σ ) in F1000W and is extremely faint in all the NIRCam data taken by the deepest survey on the sky, the JWST Advanced Deep Extragalactic Survey, JADES, implying a magnitude fainter than ∼ 30.5 at wavelengths ≲ 5 µ m. Cerberus is also undetected in the MIDIS F560W ultra-deep observations of the field, i.e., F560W ≳ 29 mag. Our analysis of MIRI F1000W mosaics produced with limited datasets as well as the detection in NIRCam stacks at ∼ 5 σ level at wavelengths longer than ∼ 3 µ m, both confirm Cerberus is real and qualifies as an extremely red object, with a color around F444W-F1000W ∼ 3.5 mag and F560W-F1000W ≳ 2 mag. The morphology of the source in the MIRI F1000W image matches that observed in the NIRCam stacked data (with 1 year passing  \nbetween both observations), in both cases presenting a main point-like component and similar-sized (i.e., unresolved) haze, further supporting that Cerberus is a real source not linked to a fortuitous alignment.  \nWe discuss the possible nature of Cerberus and identify 2 possible Galactic source types and 3 different kinds of galaxy that could match its spectral energy distribution.  \nGiven that the MIRI and NIRCam observations were taken with a one year epoch difference, we discard a Solar System origin. The NIRCam-to-NIRCam and NIRCam-to-MIRI colors do not match any known substellar object, even the coolest and reddest brown dwarfs present mid-infrared colors that are not compatible with the spectral energy distribution of Cerberus .  \nConcerning an extragalactic origin, we identify 3 possibilities:  \n- · A low-redshift solution, z ∼ 0 . 4, would imply that Cerberus is a dusty galaxy with strong emission from polycyclic aromatic hydrocarbons or warm dust heated by an active galactic nucleus. The low mass obtained for this solution, M ⋆ ∼ 10 5 M ⊙ , would point to an unknown type of dwarf galaxy with large dust content and no PAH depletion.\n- · The second possibility is a dusty starburst or post-starburst galaxy, or a galaxy hosting a midinfrared bright obscured active galactic nucleus, at z ∼ 4 with a stellar mass M ∼ 10 8 M ⊙ and high extinction.\n- · The third possibility is a z ∼ 15 galaxy, the red color emanating from a spectral energy distribution dominated by emission from an obscured active nucleus with a bright torus significantly contributing to the flux in the rest-frame optical range, with shorter wavelengths dominated by a M ⋆ ∼ 10 7 M ⊙ galaxy.  \nDespite vetting all the imaging data meticulously, confirming beyond doubt that the source is real, we have to accept that the available data are insufficient to establish the true nature of Cerberus with certainty. Given that the source is located in the deepest HST and JWST field, deeper imaging will likely be available in the near future and help to obtain a better understanding of this unique object. In any case, further investigation of Cerberus will challenge the capabilities of JWST, as it was the case in the past for objects discovered at the limit of photometric surveys (remarkably, with spacebased missions such as HST or Spitzer ), whose spectroscopic follow-up (typically with ground-based facilities) was tremendously exigent (but still attempted with  \nvarying success). The discovery and analysis presented in this letter demonstrates the feasibility of identifying NIRCam-dark sources when selected as extremely red objects with MIRI observations up to at least 10 µ m, maybe revealing new types of Milky Way sub-stellar bodies or of galaxy populations more easily detectable through ultra-deep photometric surveys in (some of) the reddest wavelengths probed by JWST.  \n- We thank the referee for their constructive com1\n- ments to our original manuscript. PGP-G and LC 2\n- acknowledge support from grant PID2022-139567NB3\n- I00 funded by Spanish Ministerio de Ciencia e In4\n- novaci'on MCIN/AEI/10.13039/501100011033, FEDER 5\n- Una manera de hacer Europa . KIC and EI acknowl6\n- edge funding from the Netherlands Research School 7\n- for Astronomy (NOVA). KIC acknowledges funding 8\n- from the Dutch Research Council (NWO) through 9\n- the award of the Vici Grant VI.C.212.036. MA 10\n- acknowledges financial support from Comunidad de 11\n- Madrid under Atracci'on de Talento grant 2020-T2/TIC12\n- 19971. This work was supported by research grants 13\n- (VIL16599, VIL54489) from VILLUM FONDEN. SG 14\n- acknowledges financial support from the Villum Young 15\n- Investigator grant 37440 and 13160 and the Cosmic 16\n- Dawn Center (DAWN), funded by the Danish Na17\n- tional Research Foundation (DNRF) under grant No. 18\n- 140. JPP and TT acknowledge financial support 19\n- from the UK Science and Technology Facilities Coun20\n- cil, and the UK Space Agency. TPR acknowledges 21\n- support from ERC Advanced Grant 743029 EASY. 22\n- AAH acknowledges support from grant PID202123\n- 124665NB-I00 funded by the Spanish Ministry of Sci24\n- ence and Innovation and the State Agency of Research 25\n- MCIN/AEI/10.13039/501100011033 and ERDF A way 26\n- of making Europe. JA-M, AC-G, LC acknowledge sup27\n- port by grant PIB2021-127718NB-100 from the Span28\n- ish Ministry of Science and Innovation/State Agency 29\n- of Research MCIN/AEI/10.13039/501100011033 and by 30\n- 'ERDF A way of making Europe'. SEIB is funded 31\n- by the Deutsche Forschungsgemeinschaft (DFG) under 32\n- Emmy Noether grant number BO 5771/1-1. JM, G O 33\n- and AB acknowledge support from the Swedish Na34\n- tional Space Administration (SNSA). RAM acknowl35\n- edges support from the Swiss National Science Foun36\n- dation (SNSF) through project grant 200020 207349. 37\n- Some of the data presented in this paper were ob38\n- tained from the Mikulski Archive for Space Telescopes 39\n- (MAST) at the Space Telescope Science Institute. The 40\n- specific observations analyzed can be accessed via DOI: 41\n- 10.17909/je9x-d314 DOI: 10.17909/je9x-d314. 42", 'REFERENCES': "Kolb, J., Madec, P.-Y., Arsenault, R., et al. 2016, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9909, Adaptive Optics Systems V, ed. E. Marchetti, L. M. Close, & J.-P. V'eran, 99092S, doi: 10.1117/12.2232788  \nKriek, M., & Conroy, C. 2013, ApJL, 775, L16,  \ndoi: 10.1088/2041-8205/775/1/L16  \nKroupa, P. 2001, MNRAS, 322, 231,  \ndoi: 10.1046/j.1365-8711.2001.04022.x  \n- Labb'e, I., van Dokkum, P., Nelson, E., et al. 2023a, Nature, 616, 266, doi: 10.1038/s41586-023-05786-2\n- Labb'e, I., Greene, J. E., Bezanson, R., et al. 2023b, arXiv e-prints, arXiv:2306.07320,  \ndoi: 10.48550/arXiv.2306.07320  \nLangeroodi, D., & Hjorth, J. 2023a, ApJL, 957, L27, doi: 10.3847/2041-8213/acfeec  \n-. 2023b, arXiv e-prints, arXiv:2307.06336,  \ndoi: 10.48550/arXiv.2307.06336  \n- Langeroodi, D., Hjorth, J., Chen, W., et al. 2023, ApJ, 957, 39, doi: 10.3847/1538-4357/acdbc1\n- Le Floc'h, E., Papovich, C., Dole, H., et al. 2005, ApJ, 632, 169, doi: 10.1086/432789\n- Leggett, S. K., Tremblin, P., Phillips, M. W., et al. 2021, ApJ, 918, 11, doi: 10.3847/1538-4357/ac0cfe\n- Leja, J., Carnall, A. C., Johnson, B. D., Conroy, C., & Speagle, J. S. 2019, ApJ, 876, 3, doi: 10.3847/1538-4357/ab133c\n- Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415, doi: 10.1146/annurev-astro-081811-125615  \nMadec, P. Y., Arsenault, R., Kuntschner, H., et al. 2018, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 10703, Adaptive Optics Systems VI, ed. L. M. Close, L. Schreiber, & D. Schmidt, 1070302, doi: 10.1117/12.2312428  \n- Magnelli, B., Elbaz, D., Chary, R. R., et al. 2009, A&A, 496, 57, doi: 10.1051/0004-6361:200811443\n- Maiolino, R., Scholtz, J., Curtis-Lake, E., et al. 2023, arXiv e-prints, arXiv:2308.01230,  \ndoi: 10.48550/arXiv.2308.01230  \n- Marley, M. S., Saumon, D., Visscher, C., et al. 2021, ApJ, 920, 85, doi: 10.3847/1538-4357/ac141d\n- Matthee, J., Naidu, R. P., Brammer, G., et al. 2023, arXiv e-prints, arXiv:2306.05448,  \ndoi: 10.48550/arXiv.2306.05448  \nMcCarthy, P. J. 2004, ARA&A, 42, 477,  \ndoi: 10.1146/annurev.astro.42.053102.134032  \nMeisner, A. M., Leggett, S. K., Logsdon, S. E., et al. 2023,  \nAJ, 166, 57, doi: 10.3847/1538-3881/acdb68  \n- Melinder, J., Ostlin, G., Hayes, M., et al. 2023, ApJS, 266, 15, doi: 10.3847/1538-4365/acc2b8  \nTable 1. Photometric measurements for the Cerberus source.  \nNote -The table provides photometric measurements for Cerberus in all NIRCam individual images from JADES, incremental (named red and blue stacks) and disjoint (adding 2, 3, 4, and 5 contiguous filters) stacks, and MIRI bands from MIDIS. The central wavelengths and widths are given in µ m, and the photometry in AB magnitudes, including 5 σ upper limits for data points with S/N < 3.", 'A. FITTING CODES USED TO ANALYZE CERBERUS SPECTRAL ENERGY DISTRIBUTION': "In this Appendix, we give details of the codes used to analyze the SED of Cerberus and obtain its possible redshifts if it were an extragalactic object, as well as the implications of those redshifts in its physical properties.  \nWe first used the eazy code (Brammer et al. 2008) to fit the data in 2 ways: one using actual (low S/N ) flux measurements for all filters, and another replacing all S/N < 3 measurements by 5 σ upper limits that templates were not allowed to surpass (achieved with a modified version of the code). We employed in the fits all v1.3 templates for stellar-dominated galaxies 1 , which include a dusty galaxy with a high-equivalent width emission-line spectrum. We also allowed combinations of stellar-only templates with the new models for little red dots and high-redshift AGN+torus recently added based on JWST data (Killi et al. 2023). We did not impose any prior and worked with minimum χ 2 photometric redshift estimates in the range 0 < z < 20.  \nWe also ran BAGPIPES (Carnall et al. 2018). BAGPIPES is a stellar population synthesis modeling package built on the updated Bruzual & Charlot (2003) spectral library with the 2016 version of the MILES library (Falc'onBarroso et al. 2011). It uses a Kroupa (2001) IMF. We adopted a Calzetti et al. (2000) dust attenuation allowing 0 < A V < 8 mag, and included nebular emission lines. The SFH is set to delayedτ model and the code included an AGN component as in Carnall et al. (2023).  \nWe fitted the extracted photometry of Cerberus with prospector (Johnson et al. 2021) to constrain its photometric redshift and stellar population properties. We adopted the prospector setup described in detail in Langeroodi et al. (2023) and Langeroodi & Hjorth (2023b). In brief, the SFH is modelled non-parametrically in five temporal bins with a continuity prior (see e.g., Leja et al. 2019); nebular emission is added from the cloudy (Chatzikos et al. 2023) runs compiled in Byler et al. (2017); we treat the nebular- and stellar- metallicity as independent free parameters; dust attenuation is modelled with a two-component model, one for the entire galaxy and one for the star-forming regions (Kriek & Conroy 2013). In a first prospector run, redshift is fitted as a free parameter with a flat prior in 0 < z < 20. We then re-ran the same prospector setup, fixing the redshift to the best-fit values mutually agreed using results from all the SED-fitting codes used in this work (see below).  \nComplementarily, we made use of LePHARE++ , the C++ latest evolution of the SED-fitting code LePHARE (Arnouts et al. 2002; Ilbert et al. 2006) to fit the Cerberus photometry, following the approach of Moutard et al. (2024, in prep.). In brief, we considered stellar population synthesis models from Bruzual & Charlot (2003) with a Chabrier (2003) IMF and two metallicities ( Z ⊙ , 0.5 Z ⊙ ), assuming exponential SFHs with 0 . 1 ≤ τ ≤ 30 Gyr, and delayed SFHs peaking after 1 and 3 Gyrs. To take into account the strong contribution from the nebular emission lines that can occur at very young ages, they are added following the line ratios adopted by Saito et al. (2020) and their normalization is allowed to vary by a factor of four. Aiming to take into account extremely dusty galaxies, attenuation is considered through the laws of Calzetti et al. (2000) and Arnouts et al. (2013), with 0 ≤ E ( B -V ) ≤ 1 . 5 mag, and dust IR reemission is taken into account following B'ethermin et al. (2012). We finally ran LePHARE++ on a grid between redshift z = 0 and 20.  \nApart from the previous codes, which were used to both constrain a photometric redshift and determine the physical properties of Cerberus , we also performed SED modeling with synthesizer-AGN and cigale , first fixing the redshift (to several values, see below).  \nThe synthesizer-AGN code assumes that the SED can be modeled with a composite stellar population (P'erezGonz'alez et al. 2003, 2008) and AGN emission coming from the accretion disk and the dust torus (P'erez-Gonz'alez et al. 2024). The stellar emission includes a young and a more evolved star formation event, each one described by a delayed exponential function with timescales between 1 Myr and 1 Gyr, and with ages from 1 Myr up to the age of the Universe at the redshift of the source. The attenuations of the emission from each stellar population are independent and described by a Calzetti et al. (2000) law, with A V values ranging from 0 to 10 mag for each population. The stellar emission is described by the Bruzual & Charlot (2003) models, assuming a Chabrier (2003) IMF with stellar mass limits between 0.1 and 100 M ⊙ , and the nebular emission is also considered (P'erez-Gonz'alez et al. 2003). The  \nAGN emission is modeled with a QSO average spectrum (Vanden Berk et al. 2001; Glikman et al. 2006). The dust emission from the AGN is modeled with the self-consistent templates of AGN tori presented in Siebenmorgen et al. (2015), and dust emission linked to star formation is added using the models from Draine & Li (2007).  \nCIGALE (Boquien et al. 2019) SED-fitting analysis assume an SFH modelled using a constant star-formation rate with ages ranging from 1 to 100 Myr. We adopted the stellar populations models from Bruzual & Charlot (2003) with solar metallicity, and the Chabrier (2003) IMF. We included nebular continuum and emission lines using solar metallicity, electron density of 100 cm -3 , and ionised parameter equal to log( U )= -2. The dust attenuation and far-IR emission uses the Calzetti Law (Calzetti et al. 2000) and the Draine et al. (2014) models, respectively. AGN emission is added using the Fritz et al. (2006) models following the initial parameters suggested by Ciesla et al. (2015).  \nThe obtained photometric redshifts solutions are: z ∼ 0 . 4 (provided by prospector and LePHARE++ ), z ∼ 4 (very prominent solution for by eazy , also detected weakly by prospector and LePHARE++ ), z = 14 -15 (obtained by all codes), z ∼ 17 (obtained by LePHARE++ ), and z ∼ 19 (obtained by bagpipes and LePHARE++ ). We note that the z ∼ 0 . 4 solution corresponds to a PAH line entering the F1000W passband, the z ∼ 4 solution would imply Pa α covered by the filter, for z ∼ 15, H α would contribute to the MIRI-F1000W flux, and z ∼ 19 means that the F1000W emission could be enhanced by H β and/or [O III ] λλ 4960,5008 emission.  \nWe fitted the SED shown in Figure 4, which includes all flux data points with S/N > 3 and upper limits for the rest of bands, using the codes described above. We checked that our results concerning the photometric redshift and physical properties (for the two parameters mentioned in the main text, the stellar mass and the dust attenuation) do not change significantly if we only use S/N > 5 fluxes and upper limits, or just S/N > 5 fluxes. More specifically, the photometric redshift solutions are the same, but with different statistical weights (i.e., the peaks vary their relative strength). Concerning the physical properties, variations were very similar to the systematic differences between the distinct estimations for the two properties that are consistent among all codes and techniques, given that they are constrained by the two robust observational properties of Cerberus , i.e., the low luminosity and red color. More specifically, for the stellar mass we obtained variations of 0.3 dex when using higher signal-to-noise data for z ∼ 4 and z ∼ 15, and 0.7 dex for z ∼ 0 . 4. For the dust content, we obtained differences of < 0 . 3 mag in all cases, always remaining relatively high, above 3 mag, or finding a similarly obscured torus model."}