Dusty starforming galaxies DSFGs significantly contribute to the stellar buildup in galaxies during cosmic noon the peak epoch of cosmic star formation. Major mergers and gas accretion are often invoked to explain DSFGs prodigious star formation rates SFRs and large stellar masses. We conducted a spatially resolved morphological analysis of the restframe ultravioletnearinfrared 0.251.3 m emission in three DSFGs at z  2.5. Initially discovered as carbon monoxide CO emitters by NOrthern Extended Millimeter Array NOEMA observations of a bright SSUB350 mSUB  111  10 mJy Herschel source we observed them with the James Webb Space TelescopeNIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies stellar populations and dust distributions on scales of 250 pc. Spatial variations in stellar mass SFR and dust extinction are determined in resolved maps obtained through pixelbased spectral energy distribution fitting. The CO emitters are massive MSUBstarSUB  3  3010SUP10SUP MSUBSUB dusty starburst galaxies with SFRs ranging from 340 to 2500 MSUBSUB yrSUP1SUP positioning them among the most active starforming galaxies at 2 lt z lt 3. Notably they belong to the 1.5 of the entire JWST population with extremely red colors. Their morphologies are disk like Srsic index n  1 with effective radii of 2.04.4 kpc and exhibit substructures such as clumps and spiral arms. The galaxies have dust extinctions up to ASUBVSUB  57 mag extending over several kiloparsecs with asymmetric distributions that include offcenter regions resembling bent spiral arms and clumps. The nearinfrared dustattenuation curve in these sources deviates from standard laws possibly implying different duststar geometries or dust grain properties than commonly assumed in starburst galaxies. The proximity lt 5 of galaxies with consistent redshifts strong color gradients an overall disturbed appearance asymmetric dust obscuration and widespread star formation collectively favor interactions minor mergers and flybys as the mechanism driving the CO galaxies exceptional SFRs. The galaxies large masses and rich environment hint at membership in two protostructures as initially inferred from their association with a Planckselected highz source.
['dust', 'extinction', 'galaxies: starburst', 'galaxies: star formation', 'galaxies: structure', 'submillimeter: galaxies', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
JWSTs PEARLS Resolved study of the stellar and dust components in starburst galaxies at cosmic noon
{"JWST's PEARLS: Resolved study of the stellar and dust components in starburst galaxies at cosmic noon": "M. Polletta 1 , B. L. Frye 2 , N. Garuda 2 , S. P. Willner 3 , S. Berta 4 , R. Kneissl 5 , 6 , H. Dole 7 , R. A. Jansen 8 , M. D. Lehnert 9 , S. H. Cohen 8 , J. Summers 8 , R. A. Windhorst 8 , J. C. J. D'Silva 10 , 11 , A. M. Koekemoer 12 , D. Coe 12 , 13 , 14 , C. J. Conselice 15 , S. P. Driver 10 , N. A. Grogin 12 , M. A. Marshall 11 , 16 , M. Nonino 17 , R. Ortiz III 8 , N. Pirzkal 12 , A. Robotham 10 , R. E. Ryan, Jr. 12 , C. N. A. Willmer 2 , H. Yan 18 , V. Arumugam 4 , C. Cheng 19 , 20 , H. B. Gim 21 , N. P. Hathi 12 , B. Holwerda 22 , P. Kamieneski 8 , W. C. Keel 23 , J. Li 10 , M. Pascale 24 , H. Rottgering 25 , B. M. Smith 8 , and M. S. Yun 26  \n(A ffi liations can be found after the references)  \nReceived 09 May 2024 / Accepted 29 June 2024", 'ABSTRACT': "Dusty star-forming galaxies (DSFGs) significantly contribute to the stellar buildup in galaxies during 'cosmic noon,' the peak epoch of cosmic star formation. Major mergers and gas accretion are often invoked to explain DSFGs' prodigious star formation rates (SFRs) and large stellar masses. We conducted a spatially resolved morphological analysis of the rest-frame ultraviolet / near-infrared ( ∼ 0.25-1.3 µ m) emission in three DSFGs at z ≃ 2 . 5. Initially discovered as carbon monoxide (CO) emitters by NOrthern Extended Millimeter Array (NOEMA) observations of a bright ( S 350 µ m = 111 ± 10 mJy) Herschel source, we observed them with the James Webb Space Telescope / NIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies' stellar populations and dust distributions on scales of 250 pc. Spatial variations in stellar mass, SFR, and dust extinction are determined in resolved maps obtained through pixel-based spectral energy distribution fitting. The CO emitters are massive ( M star ≃ (3 -30) × 10 10 M ⊙ ), dusty starburst galaxies with SFRs ranging from 340 to 2500 M ⊙ yr -1 , positioning them among the most active starforming galaxies at 2 < z < 3. Notably, they belong to the ∼ 1.5% of the entire JWST population with extremely red colors. Their morphologies are disk like (Sérsic index n ≃ 1), with e ff ective radii of 2.0-4.4 kpc, and exhibit substructures such as clumps and spiral arms. The galaxies have dust extinctions up to AV = 5-7 mag extending over several kiloparsecs with asymmetric distributions that include o ff -center regions resembling bent spiral arms and clumps. The near-infrared dust-attenuation curve in these sources deviates from standard laws, possibly implying di ff erent duststar geometries or dust grain properties than commonly assumed in starburst galaxies. The proximity ( < 5 '' ) of galaxies with consistent redshifts, strong color gradients, an overall disturbed appearance, asymmetric dust obscuration, and widespread star formation collectively favor interactions (minor mergers and flybys) as the mechanism driving the CO galaxies' exceptional SFRs. The galaxies' large masses and rich environment hint at membership in two proto-structures, as initially inferred from their association with a Planck -selected highz source.  \nKey words. Large scale structure - Submillimetre: galaxies - Star forming galaxies - Galaxies: clusters: general", '1. Introduction': 'The star formation rate (SFR) density of the Universe reached its maximum about 10 Gyr ago or at 2 < ∼ z < ∼ 3 (also called "Cosmic Noon;" Madau & Dickinson 2014). This epoch marks the theoretical transition between di ff erent growth regimes for galaxies, the peak of supermassive black holes\' (SMBHs\') accretion rate (Delvecchio et al. 2014), and the collapse of cosmic structures (Dekel et al. 2009b). At z ≳ 2 galaxies go through a rapid growth phase during which stars are formed \'in situ\' within the galaxy from infalling cold gas. This early growth is followed by a slower phase during which the stars that formed \'ex situ\' are accreted through merger events (e.g., Oser et al. 2010). In galaxy protoclusters, environment-related quenching mechanisms are likely e ffi cient by z ∼ 2 producing substantial quenched galaxy populations and forming the first true galaxy clusters (Alberts & Noble 2022).  \nThe main contributors ( ∼ 85%) to the cosmic SFR density at z ∼ 2-3 are dusty galaxies with extreme SFRs-hundreds to thousands of solar masses per year (Dunlop et al. 2017). These galaxies are usually identified as bright submillimeter (submm) sources and are frequently referred to as dusty star-forming  \ngalaxies (DSFGs) or submm galaxies (SMGs). SMGs were discovered in the late 1990s as bright ( > 1-10 mJy) sources at 450 and 850 µ m(Blain & Longair 1993; Clements et al. 1993; Smail et al. 1997). They are usually massive galaxies ( M star ∼ 10 11 M ⊙ ) with large molecular gas reservoirs (Frayer et al. 1998; Greve et al. 2005; Bothwell et al. 2013).  \nIn spite of the importance of DSFGs and decades of study, it is still unclear which mechanisms power their intense SFRs and how they acquire the gas to fuel these high rates. There is a substantial controversy in interpreting the observations. On the one hand, there is evidence for rapid growth through major merger-driven events, where interactions with similarly massive and gas-rich galaxies can trigger a short-lived and intense burst of star formation (Engel et al. 2010; Casey 2016; Ginolfi et al. 2020; Gómez-Guijarro et al. 2018). On the other hand, DSFGs are often large, gas-rich disk galaxies, implying secular growth through gas inflows and disk instabilities (e.g., Hodge et al. 2012; Gullberg et al. 2019; Hodge et al. 2019; Rizzo et al. 2020; Jiménez-Andrade et al. 2020; Amvrosiadis et al. 2023; Liu et al. 2024; Gillman et al. 2024). In this case, the galaxy interstellar medium (ISM) is continuously fed and star formation is sustained for a long time (1-10 Gyr) by gas accretion from the circumgalactic medium (CGM), the intergalactic medium (IGM), or other galaxies through flybys or minor merg-  \ns (Davis & Bureau 2016; Camps-Fariña et al. 2023; Narayanan et al. 2015). Current cosmological simulations (EAGLE, Crain et al. 2015; Illustris-TNG, Nelson et al. 2018; Pillepich et al. 2018; SIMBA, Davé et al. 2019; Romulus, Tremmel et al. 2017; Auriga Davé et al. 2019) are tuned to reproduce the evolutionary history of stellar masses and galaxy sizes as well as galaxy scaling-relations and their scatter (see Crain & van de Voort 2023, for a review). Even though independent simulations produce similarly realistic galaxies, they are often based on di ff erent prescriptions, in particular regarding the mechanisms regulating star formation, such as supernovae or feedback from an active galactic nucleus (AGN). These di ff erent prescriptions result in di ff erent gas accretion rates onto the galaxy and feedbackdriven outflow rates in the inner regions of massive galaxies\' haloes ( < 200 kpc) (Crain & van de Voort 2023). The study of the gas and star formation properties of massive star-forming galaxies (SFGs) is thus fundamental to providing observational constraints needed to refine theoretical models of galaxy evolution (Baugh et al. 2005; Davé et al. 2010; Hayward et al. 2013; Lagos et al. 2020).  \nDetermining the triggering and fueling mechanism of star formation in DSFGs has important implications on our understanding of this important galaxy population as well as of the growth of massive galaxies. Because environmental-related processes, such as interactions and gas feeding through cold streams (Dekel et al. 2009b), may a ff ect the star formation triggering and fueling mechanisms, it is also important to consider the environment in which a galaxy is evolving when studying its star formation. DSFGs are considered to be progenitors of the population of massive galaxies that is predominant in local galaxy clusters (Lilly et al. 1999; Chapman et al. 2005; Simpson et al. 2014) and signposts of matter density peaks and galaxy protoclusters (Oteo et al. 2018; Miller et al. 2018; Polletta et al. 2021; Calvi et al. 2023). The study of their star formation and gas properties in connection with the environment is thus necessary to assess the role environmental density plays in their growth and evolution. An overdense environment at high redshift is characterized by a higher merger probability and a larger gas reservoir than the field or an overdensity at low redshifts. These conditions might play a significant and perhaps decisive role in both driving and establishing the timing of the extreme SFRs in massive galaxies by triggering star formation and regulating the gas accretion rates onto galaxies.  \nThe various driving mechanisms are expected to leave di ff erent imprints on the molecular gas properties, the gas kinematics, the star formation e ffi ciency (SFE), the gas-depletion timescale, and in the stellar morphologies (Hopkins et al. 2013; Boogaard et al. 2024). In order to establish whether these intense star formers grow through major mergers or through gas accretion, we thus need to map their stellar, dust, and cold gas components. The maps require su ffi cient spatial resolution to reveal the components\' relative distribution. Resolved (with < 0 \'\' . 1 angular resolution) studies of these sources, at both near-infrared (NIR), where the e ff ects of dust obscuration are reduced, and at millimeter (mm) wavelengths, are necessary to measure the spatial extent of the stellar and dust components, establish the formation of disks and bulges, identify mergers, and trace gas fueling and ejection processes (Dekel et al. 2009b; Krumholz et al. 2018). At low resolution, minor mergers can mimic smooth rotating disks (e.g., Rizzo et al. 2022), and deviations from circular motion due to gas inflow or turbulent motions throughout the disk can be misinterpreted as signs of merging. Heavy obscuration in the galaxy center, if not accounted for, can create artificial peaks of stellar mass erroneously leading to a merger scenario (Sun et al.  \n2024). On the other hand, identifying a major merger depends on the spatial resolution and the source separation, and the latter depends on the merger stage.  \nAt mm wavelengths, facilities such as the NOrthern Extended Millimeter Array (NOEMA) or the Atacama Large Millimeter Array (ALMA) o ff er the sensitivity, resolution, and frequency coverage to map the molecular gas and the dust continuum, and constrain the gas mass and kinematics with su ffi cient resolution to test these di ff erent scenarios. These observations have also enabled the identification and study of their visible and NIR counterparts and the possibility to carry out statistical multiwavelength studies of large samples (e.g., Danielson et al. 2017; Aravena et al. 2019; Dudzeviˇci¯ut˙e et al. 2020; Donevski et al. 2020; Birkin et al. 2021; Berta et al. 2023; Liao et al. 2024).  \nThe advent of the James Webb Space Telescope (JWST; Gardner et al. 2006; Rieke et al. 2005; Beichman et al. 2012; Windhorst et al. 2008), with its unprecedented sensitivity and angular resolution at NIR wavelengths, o ff ers the opportunity to explore the rest-frame visible emission of even the dustiest galaxies at z ∼ 2. It is now possible to study their internal structure on scales of 10--100 pc or even less in lensed systems (e.g., Smail et al. 2023; Kamieneski et al. 2023, 2024b; Liu et al. 2024) and estimate their physical properties such as stellar and dust mass, SFR, and gas fraction. The JWST spatial resolution permits identifying signs of interaction and mapping their stellar structure and morphology and the close environment. Since its launch, several studies have been carried out with JWST on DSFGs, often combined with high resolution mm data. The subarcsec resolution JWST images of DSFGs at 1 . 4 < z < 10 show a diversity of structures and sizes and a wide range of morphologies (Gillman et al. 2023, 2024). Some DSFGs are gas-rich isolated massive disks (Huang et al. 2023; Liu et al. 2024), others show disturbances such as arms, bars, tidal features (Smail et al. 2023; Wu et al. 2023; Chen et al. 2022; Huang et al. 2023; Sugahara et al. 2024). The majority exhibit a compact stellar mass distribution (Chen et al. 2022; Smail et al. 2023) with a core mass fraction that increases at lower redshift possibly due to the growth of a bulge (Le Bail et al. 2024). Star formation activity is often spread over galaxy-scales (Kamieneski et al. 2024a; Crespo Gómez et al. 2024), but quiescent regions are sometimes present in the core or in the disk (Le Bail et al. 2024; Kamieneski et al. 2024a; Cheng et al. 2023). Many live in an overdense environments (groups or protoclusters) (Smail et al. 2023; Wu et al. 2023; Kokorev et al. 2023; Frye et al. 2024). In general they are heavily obscured (Smail et al. 2023; ÁlvarezMárquez et al. 2023), with dust obscuration more concentrated in the center (Price et al. 2023; Sun et al. 2024), but not exclusively (Kamieneski et al. 2024a). The emerging picture is consistent with inside-out galaxy evolution, with the growth with time of a centrally concentrated older stellar populations over more extended, younger and obscured star-forming regions that dominate the mass budget (Chen et al. 2022). Overall, these results suggest that star formation in DSFGs at highz is driven by gas accretion, rather than by major mergers (but see Sugahara et al. 2024). All these studies are revolutionizing our understanding of this emblematic galaxy population, but evidence in favor of the accretion scenario is still sparse and the amount of information to be delivered to modelers needs to be parsed and standardized. Last, these studies cover a broad range of redshifts and luminosities and the environment is rarely well characterized. Further progress is clearly needed on the study of DSFGs at high resolution. The premises are very promising, facilities such as ALMA, NOEMA and JWST hold the keys to unveil the mechanisms at the origin of the vigorous growth of these extreme star  \nformers and to finally untangle such a fascinating and enigmatic population.  \nThis paper presents JWST / Near-Infrared Camera (NIRCam; Rieke et al. 2003) and NOEMA observations of three DSFGs at z ∼ 2 . 5 with the goal of establishing the mechanisms that regulate their star formation activity. The target selection is summarized in Sect. 2. Sect. 3 describes NOEMA, LOFAR, and JWST / NIRCam observations and data reduction. Sect. 4 presents the integrated and the resolved spectral energy distributions (SEDs) and derived properties, including stellar masses, dust extinctions, and SFRs. The DSFGs\' morphological properties are analyzed in Sect. 5. Sect. 6 presents their star formation and molecular gas properties and compares them with samples from the literature and with scaling relations. Sect. 7 discusses dust properties and presents near-IR (NIR) attenuation curves. Sect. 8 compares our sources with other red JWST sources from the literature. Sect. 9 discusses the mechanisms that might be at the origin of the DSFGs\' prodigious star formation activity and the link with their environment, and Sect. 10 summarizes our conclusions.  \nThroughout this work we adopt a Chabrier (2003) initial mass function (IMF) and a flat Λ cold dark matter ( Λ CDM) model with cosmological parameters ΩΛ = 0 . 685; Ω M = 0 . 315; H 0 = 67 . 4 km s -1 Mpc -1 (Planck Collaboration et al. 2020). Magnitudes are given in the AB system (Oke & Gunn 1983). The prefix \'p\' in front of distance units given in kpc means physical, and \'c\' means comoving.', '2. Selected targets': "The selected sources were discovered as carbon monoxide (CO) emitters in NOEMA follow-up observations (Sect. 3.1) of the brightest Herschel 350 µ m source (Planck Collaboration et al. 2015) in a 20 ' × 20 ' field centered on the Planck highz source (known as PHz; Planck Collaboration et al. 2016) PHzG191.24 + 62.04 (G191 hereinafter). The Herschel source, known as Herschel ID 01 (H01, hereinafter), is a known COemitter at z = 2 . 55 from previous observations carried out with the Eight MIxer Receiver (EMIR) on the Institut de radioastronomie millimétrique (IRAM) 30-m telescope (Polletta et al. 2022). The total star formation rate (SFR) of H01 is ∼ 2300 M ⊙ yr -1 . This is derived from the total infrared (IR) luminosity ( L IR) and assuming the relation in Kennicutt & Evans (2012) corrected for a Chabrier IMF (i.e., SFR = 1 . 40 × 10 -10 × L IR) 1 . The factor adopted to derive the SFR from the total IR luminosity can di ff er by up to 50% even for the same IMF. The IR luminosity, log( L IR / L ⊙ ) = 13.21 ± 0.21, was estimated by fitting the Herschel / SPIRE fluxes (Table 1) with a modified black body and assuming the CO redshift.  \nThe G191 field (RA = 161.223 deg, Dec =+ 33.8482 deg) was included in the target list of the JWST Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) GTO Program (ID 1176, PI: Windhorst; Windhorst et al. 2023), as a protocluster candidate. The source was first identified as one of the Planck highz (PHz) sources with the largest SFR (i.e.,  \nTable 1: Herschel / SPIRE fluxes of the selected targets.  \nSFR ≳ 40000 M ⊙ yr -1 ; Planck Collaboration et al. 2015). The PHz sources are galaxy protocluster candidates at z ≃ 2-4. The interest in this field grew after observations with Herschel in the submm and with Spitzer in the mid-IR (MIR) revealed a significant overdensity of both red SPIRE sources (significant at 6.2 σ ; Planck Collaboration et al. 2015), and red Spitzer sources (Martinache et al. 2018). Such overdensities are expected in the case of highly star-forming protoclusters at z ≳ 2 (Chiang et al. 2013; Contini et al. 2016). Nine Herschel sources were observed with the IRAM 30-m telescope with the goal of detecting a CO line and determining their redshift, but a CO line was detected in only three sources with redshifts 1.89, 2.55 and 2.62. Because the undetected sources are all fainter at submm wavelengths than the detected ones, it is not possible to rule out that there might be several Herschel sources at the same redshift, but at the moment the protocluster nature of this PHz source is not spectroscopically confirmed. In this work, we study in detail the brightest Herschel source in the G191 field, H01, for which a redshift of 2.55 was confirmed through the detection of three CO line transitions (Polletta et al. 2022).  \nThe Herschel source H01 in G191 was observed with EMIR (Carter et al. 2012) at the IRAM / 30-m telescope with multiple tunings at 2, and 3 mm with the goal of detecting the CO emission line and thus constraining the source redshift and molecular gas properties. These observations yielded significant ( > 6 σ ) detections of the CO(3-2), and CO(4-3) transitions, and a tentative detection of the CO(5-4) transition at an average redshift of ⟨ z ⟩ = 2 . 55132 ± 0 . 0009 (see details in Polletta et al. 2022). Surprisingly, the lines become fainter at higher transitions, implying a CO spectral line energy distribution (CO SLED) with peak at J up = 3. The line ratio I CO(4 -3) / I CO(3 -2) = 0 . 41 ± 0 . 10 is much lower than observed in SMGs (i.e., 1.40 ± 0.02 in Bothwell et al. 2013, and 0.96 ± 0.21 in Birkin et al. 2021), and in the Milky Way (i.e., 1.12; Carilli & Walter 2013). Since the size of the EMIR beam decreases with frequency (the EMIR half power beam width varying from 27 '' at 91 GHz to 16 '' at 145 GHz) 2 , such a peculiar CO SLED could be artificially produced if multiple CO-emitters at the same redshift contribute to the large beam at short frequencies, but not all in the small beam at higher frequencies. It is also possible that there is only one CO-emitter but a wrong or inaccurate pointing might have caused signal loss in the higher frequency observations. Another possibility is that optically thick dust suppresses the highJ CO lines (Papadopoulos et al. 2010). Finally, gravitational lensing, a phenomenon that a ff ects some bright DSFGs and PHz (Negrello et al. 2010; Cañameras et al. 2015), can also produce a larger magnification of the di ff use gas emission seen in the low transitions than that coming from the compact gas traced by higher transitions, artificially producing a low peak in the CO SLED (see Hezaveh et al. 2012). To test all these possibilities, we carried out CO observations with NOEMA at a higher spatial resolution. In this work,  \nwe present the follow-up observations carried out in the millimeter (mm) with NOEMA and in the NIR with JWST / NIRCam.", '3.1. NOEMA observations': "We observed H01 in the G191 field with the NOEMA interferometer (project ID: W21DA, PI: Polletta) to assess whether multiplicity a ff ected the CO(3-2) emission previously detected with EMIR. The observations were carried out with 11 antennas with 24-368 m baseline range on 2022 March 26 and 27 for a total observing time of ∼ 5 hr including ∼ 2.6 hr on-source. The weather conditions were good with precipitable water vapor (PWV) ∼ 2 mm, average system temperature T sys ∼ 120 K, and wind speed ∼ 5 m / s. The quasar 3C 84 was used as receiver bandpass (RF) calibrator, the emission line star LkH α 101 and the BL Lac 2010 + 723 as flux calibrators, and the Seyfert 1 galaxy 1030 + 415 and the Blazar 1100 + 305 as phase and amplitude calibrators. The targeted CO(3-2) emission line at rest-frame frequency ν rest = 345.796 GHz is redshifted to ν obs = 97.37 GHz at z = 2.55. A single spectral tuning was used with the NOEMA PolyFiX correlator in dual polarization mode in the 3 mm band providing coverage over two frequency bandwidths, 79-86 GHz and 94.2-102 GHz with spectral resolution 2 MHz (corresponding to 6 km s -1 ). The spectral tuning covered the CO(3-2) line in the upper sideband and the continuum centered at 3.1 mm in the lower sideband. To search for multiple CO emitters and solve confusion, we chose the C configuration to achieve angular resolution ∼ 2 '' at the observed frequency.  \nCalibration and imaging of the u v data were carried out using the GILDAS 3 package clic with assistance from the IRAM sta ff . In the 3 mm band, the absolute flux calibration is accurate at the 10% level. We resampled the u v tables to 20 MHz ( ∼ 60 km s -1 ). With the calibrated visibility data, CO(3-2) and dust-continuum maps were created using the software GILDAS / MAPPING using natural weighting and the Hogbom cleaning algorithm (Högbom 1974) with a threshold at 1 σ , where σ ∼ 0 . 5 mJy beam -1 is the noise after deconvolution. The continuum u v table was created from the calibrated visibilities by averaging the line-free channels using the MAPPING task uv\\_continuum with lineemission channels filtered with the uv\\_filter task. The continuum emission was subtracted from the line u v data by fitting a linear baseline for each visibility with the uv\\_baseline task. The final synthesized beam size is 2 '' . 5 × 1 '' . 8 with position angle (PA) 75 · . The NOEMA primary beam (field of view) for NOEMA observations at the CO(3-2) frequency is 52 '' . 3. The resulting background root mean square (rms) noise level of the naturalweighted image cube is σ ∼ 0.5 mJy beam -1 per 60 km s -1 channel and ∼ 0.03 mJy beam -1 for the continuum image. A CO(3-2) intensity map was created by averaging the flux over the line velocity range ∼ 1200 km s -1 and applying a primary-beam correction.", '3.1.1. CO measurements': "The observations show three emission-line sources within 15 '' of the image center. One, here called 01a, is at the same frequency as the previous EMIR detection, and two more, designated 01b and 01c, are at a higher frequency corresponding to z ∼ 2.42, assuming the observed line is CO(3-2). Alternative transitions would imply redshifts ( z ∼ 1 . 28, or 3.56) that are incon-  \nwith the observed spectral energy distributions (SEDs) (see Sect. 4). The position of each line emitter, measured using the GILDAS / MAPPING GILDAS software, is reported in Table 2, and Fig. 1 shows their locations with respect to the Herschel 250 µ memission. The CO(3-2) line spectra, extracted with the GILDAS uv\\_fit task, are shown in Fig. 2. Table 2 gives the line intensities, central frequencies, and widths (FWHM) measured using the Cube Analysis and Rendering Tool for Astronomy (CARTA; Comrie et al. 2021). 4 After discovering the CO emission from 01b and 01c in the NOEMA image, we reanalyzed the existing EMIR spectra in the spectral regions where the COlines from 01b and 01c were expected and were able to identify a line at a frequency consistent with the CO line frequencies of 01b and 01c, and with S / N = 5.6. The parameters of the EMIR line associated with 01a, and with 01b and 01c, derived from a single Gaussian fit, are also reported in Table 2. Because the lines emitted by 01b and 01c do not overlap in frequency with the line from source 01a, we can rule out their contamination to the previous EMIR detection. On the other hand, these two lines overlap in the EMIR spectrum and even if their peaks are shifted by ∼ 490 km s -1 and the two sources are ∼ 5.7 '' apart, without sufficient spatial resolution it would be di ffi cult to identify the two separate sources. A search for other line emitters in the NOEMA primary beam yielded no additional detections.", '3.1.2. Dust continuum measurements': "Wedetected only 01a and 01b in the 3 mm continuum map. They are unresolved at a resolution of 2 '' × 3 '' . We measured their flux by fitting a point like source using the MAPPING task uv\\_fit in the u v plane. Table 3 gives the source coordinates and continuum flux densities. It is quite likely that multiple sources might be associated with a single Herschel source, especially at bright flux densities. For example, Montaña et al. (2021) and Scudder et al. (2016) found that 9-23% of Herschel sources with S 500 µ m > 3580 mJy are multiple. It is also common to find that some of these sources are at the same redshift (e.g., Ivison et al. 2013; Wang et al. 2016; Oteo et al. 2018; Coogan et al. 2018; Zhou et al. 2024). In PLCK G073.4 -57.5, the only Planck -selected field for which high spatial resolution mm observations are available, Kneissl et al. (2019) found between one and four ALMA continuum detections associated with each of the eight observed Herschel sources.  \nThe NOEMA observations demonstrate that source H01 in G191 su ff ers from multiplicity e ff ects and that its submm continuum emission is likely due to the combined emission of all three CO emitters. Because the Herschel observations have insu ffi cient spatial resolution ( FWHM = 25 '' at 350 µ m) to resolve the emission from each source, we estimate their submm flux by assuming that it is proportional to the NOEMA mm flux. This assumption is likely valid because the sources are at similar redshifts, and they presumably have similar dust temperatures and SEDs. Because 01c was not detected in the continuum by NOEMA, we assigned an upper limit to its submm flux. The measured and scaled Herschel flux densities are listed in Table 1.", '3.2. LOFAR observations': "The G191 field was observed with LOFAR at 144 MHz as part of the mosaic P160 + 35. The LOFAR data from the LOFAR  \nFig. 1: Multicolor (red: F356W + F410M + F444W, green: F200W + F277W, and blue: F090W + F115W + F150W) JWST / NIRCam images. Left panel: 1 ' × 1 ' image centered on the Herschel source H01. The Herschel contours at 250 µ m are in white, the radio LOFAR 144 MHz contours in yellow, and the IRAM / NOEMA CO-detected sources are indicated by magenta arrows and labeled. Right panel: 17 '' × 15 '' image showing the CO contours starting from 3 σ in steps of 3 σ (green) and the NOEMA beam (white ellipse in the top right corner).  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFig. 2: NOEMA CO(3-2) emission (gray histogram) from 01a at z = 2 . 5513 (left panel), 01b at z = 2 . 4278 (middle panel), and 01c at z = 2 . 4222 (right panel) and Gaussian bestfit model (red line). The source ID, line observed frequency and FWHM are annotated in each panel.  \n<!-- image -->  \nTable 2: CO(3-2) line measurements.Notes. ( a ) The gas mass is derived from the L ' CO(1 -0) luminosity as M gas = 1 . 36 α CO L ' CO(1 -0) and assuming α CO = 0.8 M ⊙ pc -2 (Kkms -1 ) -1 as typically assumed in starburst galaxies. The L ' CO(1 -0) luminosity is derived from the L ' CO(3 -2) luminosity assuming the brightness-temperature ratio measured in SMGs ( r 3 , 2 = L ' CO(3 -2) / L ' CO(1 -0) = 0.63; Birkin et al. 2021). ( b ) Line parameters obtained from the IRAM 30-m telescope / EMIR spectra (Polletta et al. 2022), and assuming the same r 3 , 2 value and the M gas -L ' CO(1 -0) relation adopted for the NOEMA data.  \nLoTSS second data release (DR2; Shimwell et al. 2017) include images at an angular resolution of 6 arcsec and rms noise 0.074 mJy beam -1 and a radio source list. Two LOFAR sources fall within the H01 beam and match the two CO emitters 01a and 01b (Fig. 1), but source 01c was not detected. The sources' LOFAR positions and flux densities are listed in Table 4.", '3.3. JWST observations': "JWST / NIRCam multiband observations of the G191 field were obtained as part of the PEARLS GTO Program (ID 1176, PI: Windhorst; Windhorst et al. 2023). Observations were carried out on UT 2023 April 24 in eight filters: F090W, F115W, F150W, F200W of the short wavelength (SW) module and  \nTable 3: NOEMA measurementsTable 4: LOFAR measurements  \nNotes. ( a ) 5 σ upper limit  \nF277W, F356W, F410M, and F444W of the long wavelength (LW) module. Exposure times ranged from 1.9 ks to 2.5 ks giving 5 σ point-source limits 28.32, 28.50, 28.52, 28.79, 28.49, 28.49, 27.99, and 28.29 in order of increasing wavelength for the eight NIRCam bands (Windhorst et al. 2023).  \nThe data were obtained from the Mikulski Archive for Space Telescopes (MAST) and processed with the STScI pipeline v1.7.2 (Bushouse et al. 2022) with reference files specified by the context file jwst\\_0995.pmap. Raw images were corrected for 1 / f noise, flat field, and cosmic-ray hits using the prescription of C. Willott. 5 We then ran the ProFound code (Robotham et al. 2018) for the sky subtraction to correct remaining residuals, detector-level o ff sets, snowballs, and wisps (Robotham et al. 2023). The frames were then corrected for image distortions, aligned to the GAIA DR3 reference catalog (Gaia Collaboration et al. 2023) 6 , and drizzled into final science mosaics with resampled pixel size of 30 mas, following similar methodology to that first described by Koekemoer et al. (2011) updated to use the JWST pipeline 7 . Windhorst et al. (2023) gave more details on the observations, data reduction, source detection, and photometry.  \nThe single and multiband images of the three CO emitters are shown in Figs. 3 and 4. To measure the flux densities of the three galaxies, we PSF-matched all the images to that of F444W, which has the coarsest resolution. The PSFs came from W ebb PSF models publicly provided by STScI, 8 rotated to match the orientations of the observed images. Convolution kernels were created for each F444W-other band PSF pair following the convolution theorem by computing the ratio of the Fourier transforms and inverse transforming (e.g., Aniano et al. 2011; Boucaud et al. 2016). This approach su ff ers at high spatial frequencies, where the Fourier transform of the PSF is small and dominated by noise (Aniano et al. 2011). To eliminate spurious high frequency modes induced by the noise, we applied a simple cosine-bell window function where 80% of the array values are tapered in the kernel creation. Each image was then convolved using the appropriate kernel following the procedure described  \nTable 5: JWST / NIRCam measurements.  \nNotes. α JWST and δ JWST are J2000 equatorial coordinates in h:m:s, and · : ' : '' respectively, and S filter are flux densities in the NIRCam filters in µ Jy.  \nby Pascale et al. (2022). The PSF-matched images were then stacked, and a segmentation map was created from the stacked image by selecting all pixels within a specific radius per source, from 22 to 30 pixels (equivalent to 0 '' . 66-0 '' . 9), with a value above a threshold. The threshold was derived from a wide region close to the selected sources where we computed the 3 σ -clipped mean value ( µ bck) and the standard deviation ( σ bck). The threshold value was defined as µ bck + 10 × σ bck. Nearby sources were flagged and removed from the segmentation mask. For 01a, we identified two clumps (called A1 and A2) whose association with the main galaxy is dubious in spite of their proximity because of the clumps' much bluer colors (Fig. 3). Their spectral energy distributions and photometric redshifts (Appendix A) indicate that they are distinct sources, and therefore they were masked.  \nThe total flux density of each source per band was obtained as the sum of the signal in each image of all associated pixels in the segmentation map. The associated uncertainty per source and band was computed from the variance defined as  \nσ 2 src = S src -¯ B +       NS + N 2 S NB       , (1)  \nwhere S src is the total source flux density, ¯ B the average background value obtained from seven clean circular regions with 0 '' . 72 radius around each source, NS is the number of source pixels, and NB the number of background pixels. 9 The JWST flux densities are listed in Table 5.", '4.1. Whole galaxy SED fits': "To determine the physical properties of the three CO emitters, we modeled their SEDs comprising the JWST / NIRCam data, the submm-mm data from Herschel and NOEMA, and the radio LOFAR data. The SED modeling was carried out with the Code Investigating Galaxy Emission ( CIGALE ; Boquien et al. 2019). The code provides estimates of the stellar mass, the dust mass, the dust extinction, and the SFR and constraints on the stellar-population age. We fixed the redshifts at CO(3-2) values and the metallicities to the solar value Z = 0.02. Following Osborne & Salim (2024), we used the BC03 (Bruzual & Charlot 2003) stellar population models and chose a star formation history (SFH) modeled with a double exponential with a recent  \n<!-- image -->  \n<!-- image -->  \nFig. 3: Multicolor (red: F356W + F410M + F444W, green: F200W + F277W, and blue: F090W + F115W + F150W) JWST / NIRCam images of individual CO counterparts. Panels left to right show 01a, 01b, and 01c. The size of each panel is noted at top, and north is up, east to the left. The white outlined regions in 01a indicate likely distinct sources. The center of each galaxy is marked with a black cross. Rectangular boxes mark regions used to analyze how the key physical parameters vary across each galaxy (Sect. 4.2).  \n<!-- image -->  \n<!-- image -->  \nFig. 4: Negative 3 '' × 3 '' thumbnail JWST / NIRCam images. Panels top to bottom show 01a, 01b, and 01c and left to right show increasing wavelengths as labeled. The images grayscale stretches from 0 nJy to 3.2 nJy in the first three bandpasses, and to 6.3 nJy in F200W for all three sources. The maximum values adopted in the LW filters increase with wavelength (from F277W to F444W) and di ff er for each source (from 10.6 to 42 nJy in 01a, from 10.6 to 21 nJy in 01b, and from 6.3 to 11.4 nJy in 01c). North is up and east to the left.  \nburst with constant SFR ( τ burst = 20 Gyr). These choices produce the most accurate stellar masses and SFRs among all those that are possible within CIGALE (Michałowski et al. 2014; Osborne & Salim 2024). Dust emission was modeled by Draine et al. (2014) models and dust attenuation by the two-component model of Charlot & Fall (2000, CF00 hereinafter). The CF00 model assumes a di ff erent attenuation law and strength for the birth clouds and the interstellar medium (ISM) as it might be the case when the visible-NIR and FIR emissions are not cospatial (e.g., Chen et al. 2015; Hodge et al. 2016; Smail et al. 2023). Finally, we included a radio component to model the synchrotron emission produced by the interaction of relativistic electrons from supernovae with the local magnetic field. This radio component is related to the galaxy star formation activity, and modeled as a power law with spectral index α radio = 0 . 8 ( F ν ∝ ν -α radio ; Brienza et al. 2017). Its strength is bound to the estimated SFR through the radio-FIR correlation (Helou et al. 1985). We did not include an AGN component in the model be-  \nere is no bright central source as expected in case of an AGN (Ortiz et al. 2024) and the CO spectral-line energy distribution (SLED) in 01a peaks at a low transition, contrary to what is observed in AGN (Carilli & Walter 2013). While an unobscured AGN is ruled out by the shape of the rest-frame visibleNIR SED, an obscured AGN might instead be present and contribute to the mid-IR and radio emission. No nebular component was included because all the major emission lines are outside the NIRCam filter bandpasses. The best fit was determined as the template with the lowest χ 2 , and the bestfit parameters and associated uncertainties are the likelihood-weighted means and standard deviations, respectively. To assess the parameters' reliability and choose the best model, we compared the results obtained with di ff erent SFHs (i.e., delayed and double exponential, with and without a burst), extinction curves (i.e., CF00 and the Calzetti et al. (2000) law), and dust models. A burst greatly improves the fits (it halves the reducedχ 2 ), while the choice of underlying SFH, of extinction curves, and of dust model pro-duce onl  \nFig. 5: SEDs of three CO emitters (filled magenta circles): 01a (top panel), 01b (middle panel), and 01c (bottom panel). Downward arrows represent 5 σ upper limits. The CIGALE bestfit model is shown with a solid black line. Dotted lines show the stellar light (red for the older population and blue for the younger) before dust attenuation, and the solid red and blue lines show the attenuated stellar light. The orange line shows the dust component, and the green line shows the radio synchrotron emission.  \n<!-- image -->  \nminor variations. The parameters derived assuming a double exponential or delayed SFHs are all consistent. The main parameters derived from the best fit model obtained with a double exponential SFH are listed in Table 6, and the SEDs and bestfit models are shown in Fig. 5. Based on the SED modeling, the three galaxies are massive, highly star-forming, and heavily obscured. Interestingly, the SED of source 01b exhibits a small excess of radio emission with respect to that predicted by the best fit model and based on the estimated SFR. This excess is a factor of ∼ 4 and might be an indication of synchrotron emission from a radio jet powered by an active galactic nucleus (AGN). This possibility will be discussed in Sect. 6.1.", '4.2. Resolved SED modeling': "The selected targets have NIR sizes (based on the segmentation map described in Sect. 3.3) ranging from 0 '' . 8 to 3 '' corresponding to ∼ 7-25 pkpc. Their NIRCam images show asymmetric color gradients across them (Fig. 3). These color gradients can be driven by non uniform dust distributions or variations in  \nrelative stellar ages. In order to analyze the origin of these color gradients, we built, for each source, a pixel-by-pixel SED using the eight NIRCam images and including only the pixels in the segmentation maps with 3 σ detections in at least four bands, but all pixels ended up having at least five band detections. The pixel flux density uncertainty was computed using eq. 1 and assuming N S = 1. These pixel-based 0.9-4.4 µ m SEDs were then modeled with CIGALE using the same model components as for the integrated SEDs but without the radio component. We thus obtained for each pixel with a su ffi ciently sampled SED a bestfit model and associated physical parameters. Fig. 6 shows resolved maps of the main physical parameters. Only pixels with AV , and M star measured at > 2 σ ( M star > 1 . 3 σ for 01c), and with SFR at > 1 σ are shown. This is because the SFRs are poorly constrained when the SEDs do not extend to long wavelengths (50 µ m-3mm). However, the smoothness of these maps indicate that the large uncertainties do not significantly a ff ect the bestfit values, and that they can be used to qualitatively analyze the parameters' spatial variations.  \nThe resolved maps in Fig. 6 show that the derived physical parameters have an asymmetric distribution across each galaxy, and there are several substructures that are not symmetrically distributed with respect to the galaxy centers. To quantify these variations, we computed physical parameters in rectangular boxes, shown in Fig. 3. For 01a and 01b, two boxes are needed to capture variations at di ff erent azimuths. Fig. 7 shows how parameters vary along the box. There are significant and asymmetric variations of dust extinction across all three galaxies, with the largest values measured away from the center in all of them. This is well illustrated in the top panels of Figs. 6 and 7. Extinction values AV ≳ 7 mag are measured in the central and western regions of 01a, AV ∼ 4-5 mag to the southeast of the center in 01b, and AV gradually increases to AV ≳ 5 mag from the south to the north in 01c. The M star and SFR profiles peak instead at the center despite the asymmetric dust extinction distribution. The substructures visible in the NIRCam images show up as relatively small variations on top of quite smooth profiles.  \nResolved color maps can identify regions within a galaxy experiencing di ff erent activity levels (e.g., quiescent, star-forming, or extremely dusty and star-forming). Two diagnostic diagrams, rest-frame NUV -r vs. r -K and U -V vs. V -J (Williams et al. 2009; Arnouts et al. 2013), are widely used to separate quiescent galaxies from SFGs up to z ≲ 5 (Shahidi et al. 2020) and help distinguish dust reddening from age (Miller et al. 2022, 2023). To derive the rest-frame magnitudes necessary to compute these colors, we convolved the filter transmission curves with the pixel-based bestfit models obtained using CIGALE . All pixels in the three sources fall in the star-forming locus. A subset of pixels-35% in 01a, 5% in 01b, and 38% in 01chave red colors, that is, ( NUV -r ) > 9 . 2 -2 . 1( r -K ) and ( U -V ) > 3 . 8 -( V -J ), consistent with being highly extincted. The position of each pixel in the three galaxies is shown in Fig. 8. Fig. 9 shows the resolved maps with the color-based classification. Based on this analysis, we do not identify any quiescent region within the three galaxies. As previously pointed out, the most obscured regions in our sources are not always in the center but are preferentially on one side of each galaxy (see similar cases in Smail et al. 2023; Kamieneski et al. 2024a). Furthermore, their distribution covers wide, contiguous regions extending over several kiloparsecs.  \nPrevious studies have reported that color gradients in starforming galaxies at 0 . 5 < z < 2 are mostly due to higher dust extinction, which is usually higher in their center, and in highly inclined galaxies (Liu et al. 2016, 2017; Miller et al. 2022, 2023).  \nTable 6: CIGALE physical properties.  \nNotes. CIGALE fits are for a double-exponential SFH with an additional burst. ( a ) SFR averaged over the past 100 Myrs. ( b ) Instantaneous SFR. ( c ) Ratio between the estimated instantaneous SFR and the SFR expected for a galaxy with the same stellar mass and redshift on the main sequence (Popesso et al. 2023).  \nFig. 6: Pixel-based maps of CIGALE bestfit parameters. Panels from top to bottom show AV , stellar mass surface density, and instantaneous SFR surface density. From left to right, columns show source 01a, 01b, and 01c. In all maps north is up and east to the left, and a scale bar shows 1 '' or ∼ 8 pkpc. The green cross marks the galaxy center. Only pixels with parameters measured at > 2 σ are shown in the AV and Σ M star maps, and at > 1 σ significance in the Σ SFR maps. This is because the SFRs are poorly constrained when the SEDs do not extend to long wavelengths (50 µ m-3mm), where dust emits. The white circles and ellipses in the left panel represent regions masking distinct galaxies from 01a.  \n<!-- image -->  \nThis interpretation is supported by the color-color diagrams, such as the UVJ diagram (Wang et al. 2017), by the Balmer decrement (Nelson et al. 2016), and by di ff erent dust and stellar mass distributions (Tadaki et al. 2020; Gómez-Guijarro et al. 2022). The color gradients in our galaxies can be also attributed  \nto dust reddening, but the highest extinction values are not always in the center. Furthermore, based on the inclination of our sources, we can claim that the highest extinction values do not necessarily require a high inclination. For example, the extinction is higher in the face-on system 01a than in the edge-on  \nFig. 7: Physical parameters in the boxes shown in Fig. 3. Columns show the three galaxies as labeled, and rows from top to bottom show mean AV , M star, and instantaneous SFR. The parameters were derived from the resolved maps in bins of 2-3 pixels within 0 '' . 21 wide boxes spanning 2 '' , crossing the center and sampling the most dramatic color gradients across the whole galaxy. The abscissa shows the distance of each bin from the galaxy center, indicated by the dashed vertical line, from southeast to northwest. For 01a and 01b, the two boxes with di ff erent orientations are indicated by blue diamonds and a blue line for box1 and red triangles and a red line for box2. (One box su ffi ced for 01c.)  \n<!-- image -->  \ngalaxy 01c. Similar conclusions were drawn from other studies of optically faint sources (Gómez-Guijarro et al. 2023; Lorenz et al. 2023). This implies that the extinction in some highz SFGs may be due to dust in dense star-forming regions rather than to dust in the di ff use interstellar medium (ISM).", '5. Morphological analysis': "To further investigate the origin of the color gradients in our galaxies and compare their properties with other SFGs and JWST red sources, we performed a detailed structural study of their stellar component (e.g., Gillman et al. 2023; Kokorev et al. 2023; Smail et al. 2023; Price et al. 2023; Nelson et al. 2023). Fig. 4 shows the NIRCam single-band images of the three selected galaxies. The galaxies' appearances change significantly as a function of the wavelength. In the rest-frame UV (F090W, F115W), all three galaxies are extremely faint. In the 01a restUV images, the light from three distinct sources (App. A) in the vicinity dominate. The galaxy 01b is visible but much fainter than at longer wavelengths, and 01c is invisible in F090W and barely hinted at in F115W. In the rest-frame visible (F150WF277W), all galaxies are much brighter and are clumpy. Source 01a shows multiple clumps, extended tails to the south and to the north (the latter is more visible in Fig. 3), and two substructures that mimic spiral arms appear at λ ≥ 2 . 7 µ m ( λ rest ≥ 0 . 78 µ m). Overall, 01a resembles a disturbed dusty spiral. The appearance of 01b is consistent with a face-on spiral galaxy, but at short wavelengths only the center and one spiral arm are visible. Source 01c is an edge-on disk with several clumps along the disk major axis.  \nBecause of the dust extinction, a visible-light or single-band study of these galaxies would be a ff ected by serious observational biases. The best tracer of the stellar component and the least biased toward more massive stars is the F444W band (restframe NIR ∼ 1.3 µ m). Nebular emission does not a ff ect the analysis because the strongest emission lines at the sources' redshifts  \nfall between the selected NIRCam filters. Therefore the observed images are tracers of stellar radiation.  \nTo characterize the galaxies' morphology in a quantitative way and compare them with other galaxies, we employed two codes that model surface brightness and return morphological parameters commonly used in the literature. The two codes are statmorph 10 (Rodriguez-Gomez et al. 2019) and Galfit (Peng et al. 2002, 2010). We used the WebbPSF software (Perrin et al. 2014) to generate F444W PSFs on the same 30 mas pixel scale as our image mosaics. Image cutouts with sizes 200 × 200 pixel (6 '' × 6 '' ) are large enough to contain all of the light from the galaxy of interest and to characterize the noise. We ran statmorph only on the F444W cutouts using the segmentation map described in Sect. 3.3. The code returns the half-light radius, the Sérsic index ( n ), and the quantitative morphological indicators concentration ( C ), asymmetry ( A ), and clumpiness ( S ) (CAS; Abraham et al. 2003; Lotz et al. 2008b,a; Conselice 2014) plus Gini ( G ) and M 20. Higher CAS values indicate more concentrated, asymmetric, and clumpier light profiles, respectively. The Gini parameter defines the light distribution across the pixels. G = 1 means that the light is concentrated in a single pixel, while a uniform surface brightness gives G = 0. The M 20 parameter measures the moment of a galaxy brightest regions containing 20% of the total flux, normalised by the total light moment for all pixels. Negative M 20 values indicate a high concentration of light, although not necessarily at the center. More detailed descriptions of each morphological parameter were given by Lotz et al. (2004) and Snyder et al. (2015).  \nWe ran galfit on all seven NIRCam wide-band images, fixing the center and the position angle (PA) of all objects to the F444W values. In the case of 01a, we fit simultaneously 01a and the two bright nearby objects A1 and A3 to model their relative contributions in the overlapping regions. Object A2 is not included in the fit, because the code was not able to model it.  \nFig. 10 shows the e ff ective radii for the three galaxies. The COemitters have e ff ective radii ranging from 2.0 to 4.4 kpc, consistent with those of other sources at similar redshift and with the mass-size relation of SFGs at z = 2 . 5 (Ward et al. 2024). Fig. 11 shows the comparison with Sérsic index. All sources share similar stellar profiles with Sérsic indices n ∼ 1, as expected for a disk galaxy, and relatively compact sizes, R e ff ∼ 1-10 kpc. In summary, the stellar light of the CO emitters is well described by an exponential disk with size consistent with other SFGs at similar redshifts, and no signs of on-going major merging activity. Nonetheless, substructures such as clumps and arcs appear in the residual images. While these indicate disturbances, the substructures are not captured by the global morphological parameters, implying that their contribution to the light profile is small. The low values of the clumpiness parameter S also imply that these substructures contain a small fraction of the total light (Conselice et al. 2011). On the other hand, some studies find that clumpiness is low (consistent with local ellipticals) when derived from restframe NIR images (Baes et al. 2020) or in high redshift galaxies (Gillman et al. 2023) because of resolution e ff ects.  \n<!-- image -->  \nFig. 8: Pixel-based rest-frame NUVrK, and UVJ diagrams (Williams et al. 2009; Arnouts et al. 2013) of 01a (open circles), 01b (open triangles), and 01c (open squares). In each diagram the regions where quiescent and star-forming systems lie are outlined by a solid black line and labeled. The diagonal dashed line represents the criterion to identify the highly extincted pixels (red-shaded symbols), that is, NUV -r > 9.2 -2.1(r -K) or U -V > 3.8 -(V -J). The solid black arrow shows the change in colors in the case of AV = 1 extinction.  \n<!-- image -->  \nSince A2 is fainter than all the other objects and located in the periphery of 01a, its presence should not a ff ect the fit of 01a, and indeed it is present in the residual map. The remaining galaxies, as well as those around 01b and 01c, were masked. Both codes ran successfully and returned robust measurements (based on the quality flag for statmorph and on the χ 2 for galfit ). The uncertainty on the statmorph values were estimated by running statmorph 100 times, each time adding random noise to the cutout of magnitude equal to the original rms noise measured in the cutout. The final properties (and uncertainty) for each galaxy are defined as the median (and standard deviation) of these 100 iterations. The measured morphological parameters derived with statmorph are listed in Table 7, and the reference images are shown in Fig. C.1. The galfit input images, bestfit model and residuals are shown in Fig. C.3.  \nThe statmorph parameters can be used to classify galaxies as mergers, ellipticals, or disk galaxies (Lotz et al. 2008a; Conselice 2003; Bershady et al. 2000). The comparison between the statmorph parameters and the classification regions defined by  \nLotz et al. (2008a), Conselice (2003), and Bershady et al. (2000) is shown in Sect. C. All in all, the morphological parameters suggest the three sources are late-type galaxies with no obvious signs of mergers ( A < 0 . 35). The parameter values are consistent with those obtained for a large SMG sample using JWST images in the rest-frame NIR (Gillman et al. 2023). A recently introduced morphological indicator, sensitive to faint features in a galaxy's outer region, is the shape asymmetry As (Pawlik et al. 2016). Values of As > 0 . 2 imply strongly asymmetric and disturbed systems, 0 . 1 < As < 0 . 2 are expected in mildly asymmetric and disturbed objects, and As < 0 . 1 indicates symmetric and undisturbed systems. Sources 01b and 01c have As values well below the 0.1 threshold, implying that they are undisturbed. Source 01a has As = 0 . 19 ± 0 . 01 and would be considered mildly disturbed.", '6. Star formation properties': "The SFRs of SFGs correlate with stellar mass according to the so-called main sequence (MS) relation. We refer to Popesso et al. (2023) for a recent compilation. This relation varies with redshift: for a given galaxy mass, the SFR is higher at higher redshift. Galaxies with SFRs that are > 4 times higher than the MS relation are classified as starbursts (Rodighiero et al. 2011). Fig. 12 shows the SFRs of the CO sources with respect to the MS relation.  \nAll three CO sources have SFR > 6 times the MS value (Table 6), implying that they are undergoing a starburst phase. For two galaxies, the ratio is > 4 even considering only SFR100Myr, and it is ∼ 3 in 01c. This suggests that the starburst phase might have lasted at least 100 Myr, consistent with the model predicted duration (e.g., Hopkins et al. 2013). The rates occupy the upper envelope of the SFR distribution at all masses and are among the highest known for DSFGs at 2 < z < 3 (Fig. 12). The MS o ff sets are similar to those of the AS2COSPEC sources (median SFR ∼ 1400 M ⊙ yr -1 ; Liao et al. 2024) which by selection are the brightest ( S 870 µ m > 12.4 mJy) submm sources in the COSMOS field. Interestingly, these sources are also characterized by large extinction values ( AV = 1.9-4.4) (Liao et al. 2024). Finding four AS2COSPEC sources at 2 < z < 3 with SFR ∼ 1000 M ⊙ yr -1 in the ∼ 1 . 4 deg 2 COSMOS field implies a surface density of ∼ 3 deg -2 . The AS2UDS survey (fluxlimited survey at S 850 -870 µ m > 3.6 mJy; Dudzeviˇci¯ut˙e et al. 2020) found ∼ 10 DSFGs per square degree at z = 2-3 (Dudzeviˇci¯ut˙e  \nFig. 9: Resolved maps of the three CO sources, 01a ( left panel ), 01b ( middle panel ), and 01c ( right panel ). Colors indicate the classification of each pixel based on the NUV rK and UVJ diagrams (Williams et al. 2009; Arnouts et al. 2013). There are no pixels classified as quiescent. Those classified as star-forming are shown in blue, cyan, or purple, and those classified as star-forming with high extinction-NUV -r > 9 . 2 -2 . 1( r -K ) or U -V > 3 . 8 -( V -J )-are in red, magenta, or orange.  \n<!-- image -->  \nTable 7: F444W statmorph morphological parameters.  \nNotes. Morphological parameters derived from the F444W image using statmorph , R half circ : half-light radius that contains half of the light emitted by the galaxy for a circular aperture, n : Sérsic index, C : concentration, A : asymmetry, S : clumpiness, Gini: Gini coe ffi cient, M 20: spatial moment of the brightest quintile of pixel flux values, relative to the total moment, and As : shape asymmetry.  \nFig. 11: Stellar e ff ective (half-light) radius in F444W as a function of the Sérsic index. Filled stars show our three CO emitters identified by color as indicated in the legend. Filled diamonds show sources at 2 < z < 3 from the literature (lilac: ALMA-selected galaxies from Boogaard et al. 2024; pink: SMGs in the SMACS J0723 -7327 field from Cheng et al. 2022; sea green: some of the SMGs in PEARLS from Cheng et al. 2023; royal blue: median values for the SMG sample from Gillman et al. 2023 and brown: SMGs in the Abell 2744 field from Price et al. 2023).  \n<!-- image -->  \nFig. 10: Stellar e ff ective (half-light) radius in F444W as a function of the stellar mass. Filled stars show CO emitters identified by color as indicated in the legend. Filled diamonds show a variety of sources at 2 < z < 3 (lilac: ALMA-selected galaxies from Boogaard et al. 2024; pink: SMGs in the SMACS J0723 -7327 field from Cheng et al. 2022; sea green: some of the SMGs in PEARLS from Cheng et al. 2023; orange: HST-dark galaxies from Nelson et al. 2023; and brown: SMGs in the Abell 2744 field from Price et al. 2023). The magenta line shows the average relation derived for a sample of SFGs at z = 2 . 5 by Ward et al. (2024).  \n<!-- image -->  \net al. 2020) at all SFRs, making the AS2COSPEC sources only ∼ 30% of the total. Finding three DSFGs, of which two have SFR > 1000 M ⊙ yr -1 , within 8 '' of each other suggests that our sources might reside in an overdensity (see Sect. 9.2), or that  \nsuch sources might be more common than previously thought but need deep NIR or mm observations to identify them.  \nFig. 12: SFR as a function of the stellar mass. Filled green stars show the three selected targets. Comparison DSFGs from the literature include: the ASPECS sample of CO emitters from the Hubble Ultra Deep Field (red triangles; Aravena et al. 2019), the AS2COSPEC sample of the brightest ( S 870 µ m > 12.4 mJy) submm sources in the COSMOS field (blue triangles; Liao et al. 2024); and the AS2UDS (small orange diamonds; Dudzeviˇci¯ut˙e et al. 2020) and the ALESS (small green circles; Danielson et al. 2017) flux-limited SMG samples ( S 850 -870 µ m > 3 . 6 and 2 mJy, respectively) with follow-up ALMA observations. CO observations are available for a subset (small purple circles; Birkin et al. 2021) of these samples. Labeled magenta triangles indicate three DSFGs studied in detail by JWST: ALMA J010748.3 -173028 (aka J0107a) at z = 2 . 467 (Huang et al. 2023), Arc 1 in the PLCKG165.7 + 67.0 field at z = 2 . 24 (Frye et al. 2024) and El Anzuelo at z = 2 . 3 (Kamieneski et al. 2023). The solid and dashed black lines represent the main sequence of star formation at z = 2 . 5 and its scatter (Popesso et al. 2023).  \n<!-- image -->", '6.1. SFR assessment': "Because the Herschel beam size ( FWHM = 25 '' at 350 µ m; Gri ffi n et al. 2010) is larger than the source separations, the thermal dust emission from the CO emitters is poorly constrained. In addition to the CO sources, there might be other sources in the Herschel beam that contribute to the total submm emission. The LOFAR 144MHz and NOEMA 3mm observations provide additional tracers of star formation activity, and hence of submm emission, in their ∼ 50 '' fields of view. No sources were detected in the NOEMA continuum map and in the LOFAR image except 01a and 01b, not even 01c, for which the bestfit model estimates S 250 µ m ≃ 3 mJy and S 144 MHz ≃ 0 . 1 mJy. While other sources might be present and not detected, their contribution is likely ≲ 3% of the total submm and radio fluxes. Therefore most of the Herschel flux is due to 01a and 01b, and their summed SFRs might be overestimated by a few percent at most.  \nAn obscuration-free SFR estimate independent of Herschel comes from the radio flux density, assuming it is due only to star formation. The radio-SFR relations from the literature (Best et al. 2023; Gürkan et al. 2018; Smith et al. 2021) yield SFRradio estimates with up to 50% di ff erences (e.g., SFRradio = 1507 + 264 -224 M ⊙ yr -1 for 01a assuming the relation in Best et al. (2023), and SFRradio = 2327 + 232 -207 M ⊙ yr -1 assuming that in Gürkan et al. (2018)). Furthermore, the SFR derived from the bestfit model depends on the timescale where it is computed. To  \nassess the agreement between the SFR derived from the bestfit model and the measured radio emission, we thus compare the observed (see Table 4) and the predicted radio flux densities. The bestfit models predict a 144 MHz flux density of 1.0, 0.8, and 0.12 mJy for 01a, 01b, and 01c, respectively. As illustrated in Fig. 5, the predicted radio flux density agrees with the observed one for 01a, and 01c, while it is a factor of 2.5 smaller (8 σ ) than the observed one for 01b.  \nThe excess of radio emission in 01b could be explained if the bestfit SED underestimated the SFR or if an AGN contributes to the radio flux. The bright, red nucleus seen by NIRCam also suggests an AGN, although there is no evidence of unresolved PSF like features in its core as observed in other AGN candidates in JWST / NIRCam images (Ortiz et al. 2024). The possibility that 01b hosts an AGN is of particular interest because of the frequent connection between AGN-driven radio activity and an overdense environment (Miley & De Breuck 2008) and several observational facts suggesting that 01b might reside in a proto-structure (see Sect. 9.2). Despite the uncertainties in the radio-FIR relation, the comparison with the radio emission tends to confirm the large estimated SFRs and the reliability of the SED modeling.  \nFig. 13: Molecular gas properties (gas depletion time in the top panel and gas-to-stellar mass ratio in the bottom panel) as a function of the redshift of the CO sources (filled green stars). Comparison samples are ASPECS (red triangles: Aravena et al. 2019), AS2COSPEC (blue triangles: Liao et al. 2024), PEARLS z = 4 . 6 SMGs in Abell 1489 (magenta triangles: Smail et al. 2023), and other CO-detected SMGs (light purple circles: Birkin et al. 2021). Solid lines show the relations derived by Liu et al. (2019) for SFGs with di ff erent stellar masses ( M star = 10 9 . 6 , 10 10 . 6 , and 10 11 . 6 M ⊙ in yellow, orange, and red, respectively).  \n<!-- image -->", '6.2. Depletion times and gas-to-stellar mass ratios': "The NOEMA CO observations enable computations of the gas depletion time for the three sources. The depletion time is defined as the timescale required to exhaust the molecular gas mass with the current SFR (i.e., τ depl ≡ M gas / SFR) and is equivalent to the inverse of the SFE. The gas mass was estimated from the CO luminosity L ' CO(1 -0) as M gas = 1 . 36 α CO L ' CO(1 -0) , where α CO is the CO-to-H2 factor in units of M ⊙ pc -2 (Kkms -1 ) -1 , and the factor 1.36 takes into account the contribution of helium to the molecular gas mass (Bolatto et al. 2013). The CO-H2 conversion factor in galaxies above the MS is smaller than in normal SFGs (e.g., eq. 2 of Castignani et al. 2020). For starburst galaxies, a typically adopted value is α CO ∼ 0 . 8 (Casey et al. 2014). Adopting this value allows comparison with other SMGs in the literature. However, the correct α CO might di ff er from object to object because it depends on gas metallicity and molecular gas surface brightness (Narayanan et al. 2012). A di ff erent choice of α CO value could be that assumed for normal SFGs, that is, α CO = 3.5 M ⊙ pc -2 (Kkms -1 ) -1 (Magdis et al. 2017). This choice would result in four times higher gas masses. The L ' CO(1 -0) luminosity was derived from L ' CO(3 -2) assuming the brightness temperature ratio measured in SMGs: r 3 , 2 ≡ L ' CO(3 -2) / L ' CO(1 -0) = 0 . 63 (Birkin et al. 2021). The estimated gas masses are listed in Table 2. Fig. 13 shows the gas-depletion times as a function of the redshift. Not surprisingly, the depletion times, such as those of the AS2COSPEC sample, are much shorter (implying higher SFEs) than those typically measured in SFGs of similar mass (Liu et al. 2019). Based on these values, our sources would exhaust all their gas and stop forming stars within ∆ z ≈ -0 . 1 unless the gas is replenished.  \nAnother important parameter for assessing whether an excess of cold gas leads to the large SFRs and SFEs of our sources is the molecular gas-to-stellar mass ratio ( M gas / M star). This parameter is the fraction by which the stellar mass would grow if all the existing gas forms new stars. The estimated gas-to-stellar mass ratios are shown in the bottom panel of Fig. 13. The three sources have gas-to-stellar mass ratios consistent with the SFG scaling relation (Liu et al. 2019). Despite having gas-to-stellar mass ratios, stellar masses, and sizes (Fig. 10) similar to those of normal SFGs, our sources are more e ffi cient and powerful in forming stars. Consistent with other SMG studies (Liao et al. 2024), this suggests that what drives a galaxy to move above the MS must be greater SFE, not the amount of cold gas or stellar mass (see also Thorp et al. 2022, for a similar study in galaxy mergers). In Sect. 9, we discuss the possible mechanisms and conditions that might enhance the SFEs and SFRs in our sources.", '7.1. Dust extinction and morphology': "Our three CO sources are remarkable for their large extinction values-up to AV = 7 across kiloparsec-wide regions and AV = 2 . 9-3.8 on average for the whole systems (Table 6). An extinction AV = 7 corresponds to a molecular hydrogen column density N H2 ≃ 3 × 10 21 cm -2 (Predehl & Schmitt 1995; Nowak et al. 2012; Skalidis et al. 2024). Such large gas column densities are not typical of the di ff use ISM but of giant molecular clouds (GMCs) or of the molecular torus at the center of an AGN, both of which have sizes of hundreds of parsecs at most. Similarly high extinction values are usually confined to the central 1-2 kpc of local, luminous IR galaxies (Mayya et al. 2004; Scoville et al. 2015). Typical dusty SFGs at 2 < z < 3 have average extinc-  \nn AV ∼ 2.5 (Knudsen et al. 2005; da Cunha et al. 2015; Dudzeviˇci¯ut˙e et al. 2020), although recently several systems with exceptionally high values have been found, e.g., the SMG 850.1 at z = 4 . 26 with AV = 5.1 ± 0.2 (Smail et al. 2023), AS2COS0001.1 at z = 4 . 62 with AV = 5.83 ± 0.14 (Liao et al. 2024), about 10% of the SMGs studied by Dudzeviˇci¯ut˙e et al. (2020) have AV > 5 (and up to AV = 7.7), and 5% of Herschel -selected starburst galaxies at 0 . 5 < z < 0 . 9 have AV = 17-24 (Calabrò et al. 2019). When the extinction in SMGs is derived from the dust mass and a fixed gas-to-dust ratio, even higher values are found (median AV ≃ 540 + 80 -40 mag; Simpson et al. 2017), implying the presence of optically thick dust at submm wavelengths and no light escaping at NIR wavelengths (Papadopoulos et al. 2010). This is an intriguing possibility as it might explain the faint highJ CO emission in 01a (Polletta et al. 2022). Optically thick dust would cause more absorption at higher frequencies (higher CO transitions) than at lower frequencies (lower CO transitions). The difference in absorption of the CO emission at di ff erent transitions might be further enhanced if the denser and more excited gas is mostly located in the optically thick regions, and the less excited gas is more extended (see Fig. 9 in Papadopoulos et al. 2010). The presence of dust with such optical depths could have important repercussions on estimates of key properties such as stellar mass and SFR.  \nAssuming that dust extinction is a good tracer of dust distribution (this requires the presence of an UV-visible light source behind or mixed in with the dust), we can infer the dust morphology from the extinction map. The CO sources' extinction maps (top panels of Fig. 6) show the presence of substructures. These substructures are also perceptible in the SFR maps (bottom panels of Fig. 6), in spite of those maps' irregular coverage. These similarities support the dust related origin of these structures and rule out the possibility that they might be caused by older stellar ages. As further validation of this interpretation, the stellar age maps of the three sources show only small variations that do not match those of the extinction maps. The negligible variations in the stellar ages across each galaxy imply that the stellar populations have uniform intrinsic colors. Given constant intrinsic colors, a more detailed image of dust attenuation can be obtained simply by dividing any two images. The most leverage is to use the longest available wavelength, F444W, and a shorter wavelengths with su ffi cient signal across most of the galaxy. Fig. 14 shows the results using F277W. (We do not show the results for 01c because, with the exception of some clumps along the major axis, its edge-on orientation makes the substructure details indistinguishable.) The substructures revealed in the extinction maps are clearly visible. They resemble spiral arms with significant bends, depicted by changes in their pitch angle, and with widths of 0 '' . 12-0 '' . 17 (equivalent to 1.0-1.4 pkpc), consistent with the width of the spiral arms in nearby grand-design spiral galaxies such as the Whirlpool Galaxy (Marchuk et al. 2024). The arms contain clumps (with radii ∼ 0 '' . 05-0 '' . 09, equivalent to 0.40-0.75 pkpc) and also straight sections (up to 0 '' . 30 '' . 4 or 2.5-3.3 pkpc in length) with notable deflections ( ∼ 50 · -70 · ). The sizes and shapes of these substructures are consistent with spiral arms, but their bends are indicative of disturbances, likely due to interactions.", '7.2. Dust attenuation curve': "Attenuation corrections are critical for inferring key physical properties such as SFR, stellar mass, and dust content (Kriek & Conroy 2013; Salim & Narayanan 2020; Robertson et al. 2024; Lo Faro et al. 2017; Hamed et al. 2021; Boquien et al. 2022;  \n<!-- image -->  \nFig. 14: 3 '' × 3 '' pseudo-extinction images of 01a ( left ) and 01b ( right ). Lighter areas correspond to higher extinction, which was derived as -2 . 5 × log( S F277W / S F444W) and only for pixels with > 3 σ detections in both bands. The green crosses mark the galaxy centers.  \n<!-- image -->  \n<!-- image -->  \nFig. 15: Multicolor 3 '' × 3 '' (corresponding to 25 pkpc × 25 pkpc) JWST / NIRCam images of 01a ( left ) and 01b ( right ). Colors are as in Fig. 3. The areas used to compute the dust-attenuation curve are outlined in blue for the less attenuated regions and in red for the more attenuated. In 01a, the ring has minimum radius of 0 '' . 18 and maximum of 0 '' . 42, and in 01b, it has minimum radius of 0 '' . 18 and maximum of 0 '' . 38. The center of each galaxy is indicated by a black cross. North is up and east to the left.  \n<!-- image -->  \nMarkov et al. 2023). Meaningful corrections require distinguishing between a galaxy's attenuation curve and the dust's extinction curve (Gordon et al. 1997). The latter is directly related to the physical properties of the dust grains, such as their composition and size distribution (Fitzpatrick 1999; Nozawa 2016; Hsu et al. 2023). The attenuation curve, defined as the change in the galaxy's flux due to the presence of dust, depends on the extinction curve but also on the spatial distributions of the stars and dust (Witt et al. 1992). At high redshifts, we can expect di ff erent extinction curves than observed in the local Universe because the processes a ff ecting the formation, growth, and disruption of dust grains might di ff er. Likewise, extremely dusty galaxies might be characterized by ISM conditions resulting in atypical attenuation laws. JWST is already providing new insights on the variety of attenuation curves by investigating systems at high redshift (Witstok et al. 2023; Markov et al. 2023) and along di ff erent sightlines through resolved studies of nearby galaxies (Keel et al. 2023; Robertson et al. 2024). Here, we take advantage of the asymmetric dust extinction distribution of 01a and 01b to measure their dust attenuation curves.  \nTo measure the attenuation curve, we adopted a method similar to the pair method (Stecher 1965) but using two regions of the same galaxy instead of two stars of the same spectral type. The attenuation curve can be directly derived from the ratio between more and less attenuated SEDs of two regions with the same  \nFig. 16: Normalized dust attenuation values A λ/ AV for 01a (filled red stars) and for 01b (filled blue stars). Solid lines show various attenuation laws from the literature: the CF00 attenuation law in brown (Charlot & Fall 2000); the law measured in the inner regions of starburst galaxies in magenta (Calzetti et al. 2000, Calzetti00); a Milky Way extinction law in orange (Cardelli et al. 1989, CCM89); an extinction law based on red clump stars in green (Wang & Chen 2019); an attenuation law from a sample of galaxies at 0 . 75 < z < 1 . 5 in cyan (Battisti et al. 2022) and the extinction law measured in the Small Magellanic Cloud bar in purple (Gordon et al. 2003). The dashed red and blue lines represent the bestfit cubic polynomial to the 01a and 01b values, respectively.  \n<!-- image -->  \nintrinsic spectrum, provided the attenuation law is the same for the chosen pair of regions. If the underlying stellar populations, the dust mixtures, and the radiation fields are the same in the two paired regions, any SED di ff erence between the two regions must come from di ff erent dust extinctions. In both 01a and 01b, the western and eastern sides exhibit exhibit di ff erent apparent colors (Fig. 3) that resulted in di ff erent extinction values from the resolved SED fitting procedure (see top-left and top-middle panels of Fig. 7 and top panels in Fig. 6). In each galaxy, we selected two regions, one per side, as shown in Fig. 15. The two regions were selected at the same distance from the center of each galaxy and of approximately the same size. The attenuation values were derived as  \nA λ = -2 . 5 × log( S red / S blue ) , (2)  \nwhere S red is the flux density of the more reddened SED, and S blue that of the less reddened SED (following Battisti et al. 2022) normalized by the region size. Fig. 16 shows the normalized attenuation curves of 01a and 01b derived for all NIRCam bands in all locations where there is > 5 σ signal. The 01a curve is not constrained at rest-frame λ < 0 . 5 µ m because the western region is undetected in the F090W, F115W, and F150W bands. In the visible-NIR ( λ rest > 0 . 3 µ m), the dust attenuation curve can be expressed by a cubic polynomial (Battisti et al. 2022). The bestfit models for 01a and 01b are also shown in the figure.  \nThe attenuation curves measured in 01a and in 01b di ff er significantly from each other. The 01b curve is flatter than most curves, while that of 01a is as steep as the extinction law measured in the bar of the Small Magellanic Cloud (Gordon et al. 2003). The observed diversity could be due to di ff erences in the dust-to-star geometry, a varying grain size distribution, or a combination of both (Cardelli et al. 1989; Fitzpatrick 1999; Kriek & Conroy 2013; Nozawa 2016; Hsu et al. 2023). Shallow attenuation curves, as derived for 01b, have been also measured in other highly extincted SFGs at z ∼ 1 . 5-3.0 (Salmon et al. 2016; Lo Faro et al. 2017; Hamed et al. 2023). The slope of the attenuation  \ncurve slope in these SFGs correlates with total extinction (with shallower curves in more-extincted SFGs), a result attributed to the e ff ects of scattering, star-dust geometry, and dust grain size (Salmon et al. 2016). A shallow attenuation curve can be the result of a patchy dust distribution because the emerging visibleNIR light of a galaxy is dominated by the least attenuated stars (Witt et al. 1992; Witt & Gordon 2000). With 01a more extincted than 01b, it is surprising that 01a exhibits a steeper attenuation curve. This behavior might indicate that the assumptions made to derive the attenuation law are not valid or that the correlation between extinction and the slope of the attenuation law breaks in some specific cases.  \nAs stated above, the attenuation curve depends on both the physical properties of the dust grains (i.e., composition and grain size) and the grains' spatial distribution with respect to the stars and the clumpiness of the ISM. Our sources exhibit extended regions with extremely high extinction values. It is therefore plausible that these regions might be characterized by dust properties and geometries that are di ff erent from those measured along the sightlines used for the 'standard' laws. Notably, the semianalytic galaxy-formation model SHARK (Lagos et al. 2018) is capable of reproducing several observed properties of SMGs (e.g., SEDs, number counts, and redshift distributions) thanks to the adoption of attenuation curves that scale with the dust surface density (Lagos et al. 2019). Furthermore, Donnan et al. (2024) showed that di ff erential extinction is necessary to reproduce the spectra of highly obscured galaxies, where neither simple screens nor uniformly mixed dust distributions fit the data. To interpret our attenuation curves, it would be useful to have predictions from simulations for a variety of sightlines, dust grain distributions, and dust formation and evolution models (see e.g., Cochrane et al. 2024).  \nAdopting a uniform 'standard' law in retrieving a galaxy's stellar properties might bias the estimates of stellar mass and SFR. A steeper attenuation curve at longer wavelengths translates into smaller stellar masses compared to standard laws and vice versa. The attenuation curve might vary within the same galaxy if it depends on the extinction, for example, in the di ff use ISM versus molecular clouds. High redshift DSFGs, with their heavily obscured sightlines, o ff er the opportunity to explore such a relation. DSFGs are important contributors to the cosmic SFR and are thought to be the precursors of the most massive local galaxies, and therefore properly quantifying their dust attenuation and overall properties is fundamental to understand galaxy evolution.", '7.3. Dust-to-stellar mass ratio': "The dust-to-stellar mass ratio (DSR) might provide indications on a galaxy's evolutionary stage and on dust-destruction processes (Calura et al. 2017; Popping et al. 2017; Donevski et al. 2020). To compute the DSRs of our sources, we considered the dust and stellar masses derived from fitting the whole SED with CIGALE (Table 6 and Sect. 4). SMGs show a wide range of DSRs, as shown in Fig. 17. The two dex range is much larger than the range in dusty MS and SB galaxies (Donevski et al. 2020), and the three CO emitters cover more than half of the large SMG range. The vast majority of SMGs exhibit higher DSRs than those predicted for galaxies with log( M star / M ⊙ ) = 10 . 5 -11 . 5 by the fiducial model of Popping et al. (2017). High DSRs are usually associated with e ffi cient and rapid dust grain formation (Donevski et al. 2020). The wide range for SMGs might be due to di ff erences in baryonic and stellar masses, dust extent, or evolutionary phase. DSRs depend on local conditions  \nFig. 17: M dust / M star versus redshift of our sources (green stars), and of various SMG samples from the literature (AS2COSPEC (blue triangles; Liao et al. 2024), PEARLS SMGs in Abell 1489 (magenta triangles; Smail et al. 2023), AS2UDS (orange diamonds; Dudzeviˇci¯ut˙e et al. 2020), ALESS (teal circles; Danielson et al. 2017), and the ALMA SMGs from Birkin et al. (2021) (light purple circles). The average values and 68% probability ranges for main sequence (MS) and starburst (SB) dusty galaxies from the sample of DSFGs in Donevski et al. (2020) are shown as shaded cyan and salmon areas, respectively. The gray region represents the predicted DSRs for galaxies with log( M star / M ⊙ ) = 10 . 5 -11 . 5 from the fiducial model of Popping et al. (2017).  \n<!-- image -->  \nand can vary by large factors within a galaxy, but these variations are averaged out when considering integrated quantities. The large spread of DSRs in SMGs likely reflects 'local conditions' occupying a greater fraction of the galaxies. Additional disparities might be introduced by the methods employed to estimate the dust masses (e.g., from visible-NIR SED fitting, from the mm continuum flux density, from fitting the far-IR SED with a gray body model; Berta et al. 2016), and by uncertainties associated with the stellar masses. Our dust-mass estimates are poorly constrained because of the large uncertainties associated with the submm flux densities (Table 3). New ALMA observations of 01a will enable a more accurate estimate of the dust extent and mass and a study of the resolved DSR. Such a study might reveal variations associated with the level of star formation activity and dust grain formation processes.", '8.1. The nature of the JWST red sources': "JWST has revealed an abundance of red galaxies, some invisible at shorter wavelengths even in HST images. In the literature, they are often referred to as little red dots (LRD; Kokorev et al. 2024; Pérez-González et al. 2024; Matthee et al. 2024; Wang et al. 2024a; Kocevski et al. 2024), HST-dark or HST-faint galaxies (Pérez-González et al. 2023; Barrufet et al. 2023; Frye et al. 2024; Kamieneski et al. 2024b), optically dark or optically faint galaxies (OFGs; Gómez-Guijarro et al. 2023), or extremely red objects (EROs; Barro et al. 2024). These objects include DSFGs, obscured AGN, and z > ∼ 5 galaxies. According to Rodighiero et al. (2023), these red NIRCam sources are either massive DSFGs at 3 < z < 7, quiescent galaxies at 3 < z < 5, or galaxies  \nwith dwarf masses and high extinction (e.g., AV ∼ 5 . 5 mag, Xiao et al. 2023; Bisigello et al. 2023). Kocevski et al. (2024) found instead that they are heavily reddened AGN, whereas PérezGonzález et al. (2024) claimed that they are mostly extremely intense and compact starburst galaxies with some contribution from an obscured AGN. Followup studies of individual sources (Fujimoto et al. 2023; Killi et al. 2023; Furtak et al. 2023; Kokorev et al. 2023; Akins et al. 2023; Killi et al. 2023; Kokorev et al. 2024; Zavala et al. 2023; Meyer et al. 2024) have found that the objects' redness is mostly due to heavy obscuration and to high redshifts. This claim is supported by simulations that reproduce these objects' properties by assuming that they are at z = 4-7 and heavily obscured ( AV = 2-4; Cochrane et al. 2024).  \nThe three CO emitters are bright and red at near- and mid-IR wavelengths and extremely faint at visible wavelengths, putting them in the category of EROs. Fig. 18 compares our targets' colors and magnitudes with those of various HST-dark galaxies, submm / mm-selected SFGs, and the full galaxy population from the JWST UNCOVER program (PIs: Labbé and Bezanson; Wang et al. 2024b). The CO emitters have colors similar to those of the HST-dark galaxies and DSFGs, which span about 3 mag. These galaxies' colors are distinct from the bulk of JWST sources, having F200W -F444W > 0 . 7, F150W -F356W > 1, and F200W -F356W > 0 . 5. Only 1.5% of the full UNCOVER sample (with > 5 σ detections) have such red colors. The CO emitters are among the sources with the brightest F444W magnitudes, ∼ 20.6-22.2 versus 22-26 for the DSFGs with missing or > 2 . 5 redshift and the majority of HST-dark galaxies. The DSFG El Anzuelo is intrinsically ∼ 2.5 mag fainter than shown in Fig. 18 once the magnification from strong gravitational lensing is taken into account (Kamieneski et al. 2023).  \nGottumukkala et al. (2024) and Cochrane et al. (2024) suggested that the JWST red sources constitute a new population. A few sources with similar red colors had already been discovered prior to JWST but at lower redshifts and higher luminosities than the bulk of the JWST red population. Some examples are the BAL QSO Mrk 231 at z = 0 . 042 (Berta 2005; Polletta et al. 2007) and the heavily obscured QSOs SWIRE J104409.95 + 585224.8 at z = 2 . 54 (Polletta et al. 2006). The predicted colors of these sources for 0 . 1 < z < 5 overlap with those of the red JWST sources, as shown in Fig. 18, which also shows the expected colors derived from the average SED of the ALESS SMGsample (Danielson et al. 2017). The ALESS colors overlap with the red JWST sources only for z > 2 . 5. The visible-NIR colors of the red JWST sources are thus consistent with those of SMGs at z > 2 . 5 or with heavily obscured AGNs that were known before JWST. In the color-magnitude diagram (Fig. 18 bottom), the three CO sources fall among the mm-selected DSFGs. Most of these have F444W magnitudes < 22, in contrast to the optically selected, HST-dark red sources, most of which have F444W magnitudes between 22 and 26. Both groups have similar F200W -F444W colors, much redder than the UNCOVER galaxies. One possibility is that the HST-dark galaxies and highz DSFGs represent the same population with the latter being typically more luminous.  \nIn spite of the success of obscured-AGN templates in reproducing the 2-5 µ m colors of the red JWST sources and recent spectroscopic evidence of AGN activity in a subset characterized by blue UV spectra (Greene et al. 2024), AGN emission might not be the dominant luminosity source. Williams et al. (2024) found that the SEDs of a sample of LRDs di ff ered from the Mrk 231 template at observed wavelengths λ > 7 µ m. These authors claimed that the LRDs are moderately dusty and old galaxies rather than obscured AGN (see also Pérez-González et al.  \nFig. 18: JWST color-color and color-magnitude diagrams for a variety of galaxies. Filled green stars represent the three CO targets. Small dots represent galaxies from the UNCOVER program (Wang et al. 2024b) with their colors indicating the photometric redshift: purple for z < 0 . 5, blue for 0 . 5 < z < 1 . 0, yellow for 1 . 0 < z < 3 . 0, cyan for 3 . 0 < z < 5 . 0, and red for z > 5. Other symbols represent HST-dark galaxies (purple upside-down triangles; Barro et al. 2024; Nelson et al. 2023; Barrufet et al. 2023; Pérez-González et al. 2023; Frye et al. 2023; Wang et al. 2024a) and DSFGs (magenta triangles; Peng et al. 2023; Zavala et al. 2023; Huang et al. 2023; Cheng et al. 2023; McKay et al. 2024; Smail et al. 2023; Kamieneski et al. 2023; Sun et al. 2024; Meyer et al. 2024). Filled symbols indicate sources with spectroscopic redshifts. Dashed lines mark the typical colors of JWST red sources. DSFGs with spectroscopic redshifts ∼ 2.5 are labeled in the bottom panel (i.e., El Anzuelo at z = 2 . 3 (Kamieneski et al. 2023), J0107a at z = 2 . 467 (Huang et al. 2023), and four galaxies at 2 . 1 < z < 2 . 7 from Price et al. (2023). The lilac, salmon, and teal curves show tracks for the average SED from the ALESS SMG sample (Danielson et al. 2017), the Torus template (Polletta et al. 2007), and Mrk 231, respectively, with the solid curve showing 0 . 1 < z < 2 . 5 and the dashed curve 2 . 5 < z < 5.  \n<!-- image -->  \n2024). This claim is supported by our sources because they do not show any obvious signs of AGN activity.", '8.2. The CO sources and the SMG population': 'Based on the estimated 850 µ m fluxes from our bestfit SED model (i.e., S 850 µ m = 11 and 7 mJy for 01a and 01b, respectively), as well as on the typical 870 µ m / 3 mm flux ratio for SMGs (80 ± 40; Birkin et al. 2021), sources 01a and 01b can be considered SMGs. The comparison between their main properties and those of various SMG samples from the literature (Sect. 6) indicates that 01a and 01b are among the SMGs with the highest SFRs and shortest depletion times at 2 < z < 3. However, they do not stand out in other properties such as stellar mass and gas-to-stellar mass ratio. Their morphologies (Sect. 5) are also consistent with those of other SMG samples from the literature (Gillman et al. 2023; Cheng et al. 2023; Price et al. 2023; Boogaard et al. 2024). Interestingly, there is some evidence of an overdense environment around source 01b because of its proximity to 01c, in line with the SMG expected environmental properties. Sect. 4.2 showed that the extinction varies significantly across the galaxies. This result emphasizes the importance of resolved studies of the stellar, dust, and gaseous components of this kind of sources in order to understand their formation and evolution (e.g., Cheng et al. 2022; Smail et al. 2023; Kamieneski et al. 2023; Zavala et al. 2023; Kokorev et al. 2023) and of establishing a common framework for meaningful comparisons. This study adds one additional piece of information to the SMG picture and lays the basis for followup mm observations.', '9.1. The powering mechanism of the estimated high SFRs': "The sources examined in this work are starburst galaxies with some of the highest estimated SFRs in the galaxy population at cosmic noon. They are also quite peculiar because of their heavy o ff -center obscuration and large extinction gradients (Sect. 4.2). At the same time, they have relatively smooth and symmetric stellar-mass distributions, well described by a disk like profile with similar sizes and Sérsic indices as SMGs and normal SFGs at z ∼ 2 . 5 (Sect. 5). The residual maps obtained after subtracting an exponential disk show substructures, including arcs mimicking spiral arms and clumps, although the substructures are not dominant. Indeed modeling the observed light with a Sérsic n ∼ 1 profile yields a good fit, and the galaxies have a low clumpiness parameter. Notably, many highly SFGs at 1 < z < 3 and SMGs display disks with substructures similar to our sources (e.g., Glazebrook et al. 1995; Abraham et al. 1996; van den Bergh et al. 1996; Cowie et al. 1995; Elmegreen et al. 2007, 2009; Förster Schreiber et al. 2009; Cheng et al. 2022; Wu et al. 2023; Liu et al. 2024). In our sources, these substructures seem to be associated with variations in the obscuration level. Asymmetric dust distributions have also been reported in other SMGs such as COSBO7 (Ling et al. 2024), COSMOS 1648673 (also known as PACS-819; Liu et al. 2024), El Anzuelo (Kamieneski et al. 2023), and Arc 1 in the PLCK G165.7 + 67.0 field (Frye et al. 2024). In summary, the properties of our sources, their large SFRs, stellar disk profiles, substructures, and asymmetric obscuration are common in the 1 < z < 3 SMG population.  \nStarburst galaxies in the local Universe are the result of a major merger (mass ratio ≥ 1:4; Armus et al. 1987; Sanders & Mirabel 1996; Kartaltepe et al. 2010; Ellison et al. 2013; Engel et al. 2010; Casey 2016). However, this does not seem to be always the case at high redshifts, where di ff erent conditions, such as the availability of more gas might result in higher gas accretion rates and a di ff erent mode of star formation (Scoville et al.  \n2016; Liu et al. 2019). Starburst activity may thus be driven secularly by compressive, gravitationally unstable, gas-rich disks (e.g., Ceverino et al. 2010; Tadaki et al. 2018). The transition between di ff erent growth regimes is expected to occur at 2 ≲ z ≲ 3 (Oser et al. 2010).  \nOur CO emitters, with their exceptional SFRs, are ideal laboratories to test the predictions from the various proposed star formation driving mechanisms. If a major merger were responsible for their vigorous star formation, we would expect to see double nuclei or a galaxy pair. The lack of those features as well the CAS-based morphological diagnostics (Sect. 5) do not support the major merger scenario. However, identifying a major merger depends on the source separation and the merger stage. Morphological indicators can be quite incomplete in this regard, especially in the post-coalescence phase (Bignone et al. 2017). Similar to our sources, Engel et al. (2010) found that most SMGs are compact, single sources, but those authors argued that this is consistent with end-stage, coalesced major merger products. This interpretation is supported by merger simulations that predict that the SMG phase occurs nearly a Gyr after the last major merger (Mihos & Hernquist 1996; Springel & Hernquist 2005), and by observations of late stage mergers with low clumpiness (Calabrò et al. 2019). According to Calabrò et al. (2019), low clumpiness can be observed in late stage mergers after coalescence because of the clumps' rapid destruction by strong stellar radiation or AGN feedback, or of their migration toward the center.  \nDuring a major interaction, orbits of stars are violently perturbed, gas is torqued to the center, shocked, and compressed by tidal e ff ects boosting the SFR (Sparre et al. 2022). The highly dissipative gas collapse (due to tidal torques) causes a contraction of the stellar component and a severe reduction of the disk mass fraction (Barnes & Hernquist 1996; Bournaud et al. 2011). A major merger would then destroy the disk in favor of the buildup of a spheroid. The remnant galaxy would be an early type galaxy with a small radius and a high Sérsic index (Bournaud et al. 2011; Wuyts et al. 2010). This scenario is not consistent with the properties of our sources. However, Moreno et al. (2019) have shown, through simulations at high resolution, that the turbulence created by a merger makes the ISM stable against collapse and leads to enhanced fragmentation into cold clouds and to SFR enhancement that is radially extended. Also Springel & Hernquist (2005) demonstrated that the disk is not destroyed in a major merger between disk galaxies with large gas fractions. Cooling can quickly reform a disk, yielding a remnant that, structurally and kinematically, more closely resembles a spiral galaxy than an elliptical, but tidal torques form a central, rotating bulge. Hence the remnant is a galaxy with a bulge and an extended star-forming disk. The low Sérsic index, the centrally peaked SFR profile, and the negligible variations in stellar ages across our galaxies disfavor this possibility. Another prediction from simulations is the onset of AGN activity following a major merger (Hopkins et al. 2008). Because major mergers favor the transport of gas into the galaxy central region, the central super massive black hole (SMBH) can grow fed by this gas and eventually reach a su ffi ciently high luminosity to a ff ect the galaxy gas and star formation activity through feedback. The lack of AGN activity in our sources might thus be explained by the limited gas supply to the galaxy center predicted by this picture. Although AGN activity can be present in massive disk galaxies with no sign of interaction or mergers (e.g., Rubin's galaxy; Holwerda et al. 2021). Based on these considerations the major merger scenario as driving mechanism in our sources seems disfavored.  \nAt the same time, we cannot rule out that a major merger might have occurred hundreds of Myrs in the past. To support this claim, we refer to a beautiful example of a late-stage merger ( ∼ 500 Myrs in the past), the IR-luminous galaxy NGC 3256 at z = 0 . 01. The galaxy hosts two nuclei separated by ∼ 0 . 8 kpc (Sakamoto et al. 2014). The secondary nucleus is visible only in high resolution radio images because it is heavily obscured ( AV = 16; Lira et al. 2008). The galaxy exhibits asymmetric dust extinction that gives it an appearance at visible wavelengths that is strikingly similar to 01a and 01b. 11 It was only through the extensive multiwavelength spectroscopic and photometric analysis carried out at high resolution that its merging nature was uncovered.  \nAn alternative mechanism that can drive star formation and produce high SFRs is gas accretion. Ho et al. (2019) used EAGLE simulations to examine the e ff ects of gas accretion on star formation on galaxies' disks. Because the accreted gas lacks angular momentum to maintain circular orbits, it flows radially inward and ends up with an anisotropic distribution and forms streams and disk like structures. Most of the gas that will be converted into stars is expected to reside first in the spiral arms and then concentrate near the galaxy center. In the case of gas accretion, we thus expect to see star formation throughout the disk, detecting it to large radii (e.g., Dannerbauer et al. 2017) and along spiral arms, and also a dusty star-forming core. This picture is qualitatively consistent with the appearance of our sources. Inflowing gas might stem from the circumgalactic medium (CGM), the intergalactic medium (IGM), galactic halos and groups of galaxies (e.g., Westmeier et al. 2005), may be stripped from infalling satellites, or may infall directly via a smooth cold, flow along filamentary structures. Interactions (flybys and minor mergers) can also supply gas and induce extended star formation (Pan et al. 2019; Thorp et al. 2019). Source 01a might be currently interacting with the neighbor (at 1 '' . 5) galaxy A3 and exhibits a tail of stars to the north consistent with being stripped from the disk by the forces at play during a strong interaction. Sources 01b and 01c are close in projection (5 '' . 7 apart, or ≲ 50 pkpc) and at similar redshifts and they have three nearby galaxies ( < 5 '' ) with consistent photometric redshifts, also supporting on-going interactions (Sect. 9.2). Also the width of the observed CO(3-2) lines suggests a complex kinematics consistent with disturbances from an interaction or minor mergers.  \nObservational evidence of inflowing gas is often elusive as the gas is expected to be relatively cold (at T ∼ 10 4 K) and faint. On large scales, such gas can be traced by Ly α filaments (Daddi et al. 2022), and it is often associated with massive halos and overdensities at the core of galaxy protoclusters (Valentino et al. 2016; Daddi et al. 2022; Pensabene et al. 2024). Accreting gas streams might be also traced by inflowing and extended cold molecular gas (e.g., Herrera-Camus et al. 2020; Ginolfi et al. 2017; Berta et al. 2021; Vidal-García et al. 2021). Indeed, cooling and gravitational collapse in gas streams may lead to clump condensation out of the CGM and to star formation (Dekel et al. 2009a; Bouché et al. 2013; Pallottini et al. 2014; Ceverino et al. 2016; Nelson et al. 2016). However, signatures of infalling molecular gas are notoriously di ffi cult to observe (Berta et al. 2021). Moreover, even when observed, infalling molecular gas can have many interpretations such as unresolved satellites and merging companions or gas that was previously ex-  \nTable 8: Selected galaxies around 01b and 01c.  \nd and is falling back into the CGM, where it is able to cool (Vidal-García et al. 2021). To probe the gas-accretion scenario and definitively untangle the mechanisms that triggered and sustain star formation in our sources, we would need to measure the resolved star formation and the molecular gas kinematics (e.g., Liu et al. 2024).", '9.2. The environment of the CO emitters': "The large stellar masses of the three CO emitters and the stellarto-halo mass relation (Girelli et al. 2020) suggest that they live in massive halos ( M halo = (2 . 9 ± 0 . 6) × 10 12 M ⊙ for 01b and (2 . 4 ± 0 . 4) × 10 13 M ⊙ for 01a). Other models in the literature (e.g., Behroozi et al. 2019) would imply halo masses up to four times larger. According to Chiang et al. (2013), a halo with mass above 10 13 M ⊙ at z ≃ 2 . 5 will evolve into a cluster with mass above 10 14 M ⊙ by z = 0. Assuming the halo masses estimated above, 01a might reside in a galaxy protocluster, while 01b might be in a proto-group.  \nQuite a bit is already known about the environment of sources 01b and 01c. They are at similar redshifts ( δv < 500 km s -1 ) and at a projected distance of only ∼ 5 '' . 7, equivalent to ≲ 50 pkpc. The apparent proximity supports the proto-group or protocluster hypothesis. As stated earlier, the G191 field was selected as a protocluster candidate because of its bright Planck signal and the presence of an overdensity of Herschel sources and of red Spitzer sources. Our sources have the red mid-IR colors of members of such overdensities and are associated with a Herschel source. Furthermore, SMGs are considered tracers of overdensities (Calvi et al. 2023). Several other examples of hyperluminous SMGs with no obvious signs of a major merger but located in large galaxy overdensities are known. Examples include SMM J084933 at z = 2 . 410 (Ivison et al. 2013) and GN20 at z = 4 . 055 (Hodge et al. 2012), although recent JWST observations at mid-IR wavelengths of the latter reveal a double nucleus consistent with a late-stage merger (Colina et al. 2023). Our sources and other luminous, non binary SMGs such as those mentioned may live in overdense environments, where related conditions might provide support to their star-forming disks and large SFRs. Environment might be the common trait of all these sources that explains their extreme star formation activity.  \nTo further investigate the environment of 01b and 01c and assess the possibility that they might be members of a largescale structure, we searched for nearby galaxies that might be at z ≃ 2 . 42. There are five galaxies that are either < 2 '' away or up to 5 '' away but with similarly red colors: three around 01b and two near 01c (Fig. 19). Table 8 gives the galaxies' designated source names, coordinates, and photometric redshifts, the latter measured with with CIGALE based on NIRCam flux densities measured using the method described in Sect. 3.3.  \nThe three galaxies close to 01b have similar appearance and colors. Their photometric redshifts are all consistent with that of 01b, as expected if they were all members of a proto-group  \n<!-- image -->  \nFig. 19: Multicolor JWST / NIRCam images (with filter assignments as in Fig. 3) centered on 01b ( left ) and on 01c ( right ). Images sizes are 10 '' × 10 '' and 4 '' . 5 × 4 '' . 5, respectively. Photometric redshifts were derived for the nearby sources marked with green circles and labeled to search for potential members of an overdensity at z ≃ 2 . 42. The galaxies designated as 01b.A, 01b.B and 01b.C in the left panel have all photometric redshifts consistent with the redshifts of 01b and 01c, while 01c.A and 01c.B in the right panel have lower redshifts.  \n<!-- image -->  \nor protocluster. The 01c neighbors have bluer SEDs and lower photometric redshifts than 01b and 01c, implying that they are distinct sources. This analysis suggests that 01b and 01c might be members of a proto-structure with five members in a 10 '' × 10 '' ( ∼ 200 ckpc) region. Only the densest protocluster cores contain a similar number of members in such a small region: four members in the Distant Red Core at z = 4 . 0 (Oteo et al. 2018), five in SPT-CL J2349 -5638 at z = 4 . 3 (Miller et al. 2018), and seven in ClG J1001 + 0220 at z = 2 . 51 (Xu et al. 2023). The densest regions of most z ∼ 2 . 5 protoclusters contain from one to three members across 10 '' (Darvish et al. 2020; Polletta et al. 2021; Pérez-Martínez et al. 2023, 2024). If confirmed, the protostructure with 01b and 01c may have an exceptionally dense core.  \nThe G191 field also contains one of the most significant overdensities of red Herschel sources in the PHz sample. With 01b being associated with Herschel 01, it is possible that an overdense structure might extend to distances ∼ 5 ' , where the other Herschel sources are located (Planck Collaboration et al. 2015; Polletta et al. 2022). Spectroscopic observations of the potential members around 01b and 01c and of the galaxies in the two overdensities (the red MIR sources and the red Herschel sources) in the G191 field are necessary to confirm and characterize such a structure.", '10. Summary and conclusions': 'This work has analyzed the stellar, dust, and cold gas components of three starburst galaxies at z ∼ 2 . 5 using JWST / NIRCam imaging from the PEARLS GTO program and NOEMA observations. Their visible-radio SEDs imply SFRs of hundreds to thousands of M ⊙ yr -1 , stellar masses M star ≃ 10 11 M ⊙ , and average extinction AV ∼ 3-4 mag. Compared to the star-forming main sequence (MS) at their redshifts, their SFRs are from 6 to 17 times higher, placing these galaxies among the starburst population and among the most powerful SMGs at 2 < z < 3. Our main results are summarized below.  \n- · Morphology and structural components: The selected sources are extremely red with faint rest-frame UV fluxes and bright rest-frame visible emission. Their red colors are  \nattributed to dust obscuration. Their stellar emission has surface-brightness profiles consistent with an exponential disk (Sérsic index n ∼ 1) with e ff ective radii of a few kiloparsecs. These stellar sizes are in good agreement with the e ff ective sizes of mass-matched SFGs at z = 2 . 5. Two galaxies are almost face-on disks, while one is an edge-on disk. All three systems have however a disturbed appearance with several substructures such as clumps and spiral arms across them. These substructures exhibit extinction values up to AV ≃ 57, extend over several kiloparsecs and have shapes consistent with spiral arms that have undergone some major disturbances causing them to bend. The morphology and the physical conditions derived for these substructures suggest that interactions have formed them and that they contribute to the enhanced star formation activity of these galaxies.  \n- · Star formation properties: All three galaxies have short depletion timescales ( ≤ 100 Myr), consistent with the most powerful SFGs at their redshift. Despite their large gas reservoirs ( M gas ≳ 10 11 M ⊙ ), their gas-to-stellar mass ratios ( M gas / M star ≃ 0 . 4-1.3) are consistent with those found in the SFG population with similar stellar masses and redshifts. This result suggests that what has driven these galaxies above the MS must be linked to higher e ffi ciency and not to the amount of cold gas or stellar mass, in agreement with other SMGstudies (Liao et al. 2024). Because of their large SFEs, the galaxies are expected to deplete their large gas reservoir and quench in ≲ 100 Myr unless the gas is replenished.\n- · Dust properties: The analyzed sources are characterized by elevated extinction levels that can reach AV = 7 even in regions several kiloparsecs in size. Interestingly, these high extinction values are observed away from the galaxy centers and only on one side of each galaxy, suggesting that extinction in some DSFGs may be due to dust in dense starforming regions rather than to dust in the di ff use ISM. We computed dust attenuation curves in the rest-frame wavelength range of 0.4-1.2 µ mfor two of the systems. They differ from most of the standard attenuation curves (Cardelli et al. 1989; Calzetti et al. 2000; Wang & Chen 2019; Battisti et al. 2022), possibly implying that the dust-star geometry or dust properties (e.g., size distribution, composition) in these galaxies or along heavily extincted sightlines might be different than what commonly observed. It is beyond the scope of this work to determine new attenuation laws, but we have shown that JWST observations may be used for this purpose and for exploring the dust properties. Our analysis implies that in galaxies with extreme dust extinction a di ff erential attenuation curve with an extinction-dependent shape might be necessary to correctly estimate unattenuated fluxes and ensuing physical properties, such as the stellar mass.\n- · Comparison with other red JWST sources: Our three sources belong to the population of JWST red galaxies that usually includes DSFGs, and HST-dark objects. Their red colors are attributed to their high redshift, z ≳ 2 . 5, and heavy extinction. Sources with similarly red colors represent ∼ 1.5% of the whole JWST population. Based on the analysis of galaxy templates of various types, we find that similarly red colors are reproduced by heavily obscured AGN or by DSFGs at z > 2 . 5, and that there is no need for a new population.\n- · Star formation powering mechanism: The galaxy-scale star formation activity, the presence of a single nucleus, the lack of a prominent bulge, and the classification based on the non parametric morphological parameters disfavor the major merger scenario as powering mechanism in our sources.  \nTheir disk appearance, extended SF, signs of disturbance such as tidal tails, spiral arms with significant changes of pitch angle, clumps and an asymmetric dust extinction distribution favor instead gas inflows at the origin of their powerful SFRs and enhanced SFEs. The presence of galaxies near in projection to sources 01a and 01b and with consistent photometric redshifts (see App. A and Table 8) indicates that minor mergers and flybys might provide the necessary gas to sustain their star formation activity. This picture might also explain the lack of AGN activity because of the limited supply of gas that reaches the galaxy center, contrary to what would be expected in case of a major merger.  \n- · Environment: The CO sources are associated with two overdensities (of red Herschel sources and of red Spitzer sources) spread over several arcmin. In addition, their stellar masses imply that they are hosted by massive dark matter haloes. All sources have neighbors with similar (spectroscopic or photometric) redshifts. In summary, there are multiple indications that these CO sources might be members of a galaxy protocluster or proto-group, as first suggested by their association with a Planck -selected high redshift source.  \nIn conclusion, these three powerful SFGs at z ∼ 2 . 5 stand out with respect to the average SFG population at similar stellar masses and redshifts because of their highly e ffi cient and powerful star formation activity. Their exceptional properties are attributed to accretion of gas most likely through minor mergers and flybys that sustain their star formation and supports the disk. Instabilities of the gas can explain the spatially widespread star formation and the formation of substructures with diverse dust obscuration. These properties, also observed in other luminous non binary SMGs, might be associated with an overdense environment in a massive halo. For the future, the positions of these sources relative to any underlying protocluster will be investigated through the analysis of the galaxy population over a wide area ( ∼ 5 \' ) where overdensities of Herschel and red Spitzer sources are located. In addition, sources at similar redshifts will be searched through ground-based NIR spectroscopic observations, and CO observations at mm wavelengths. Finally, recently obtained ALMA observations of source 01a at high spatial resolution will be able to test the proposed fueling scenario. This can be achieved by measuring the gas turbulence in the disk and through the search of nonrotating gas components with a resolved study of the gas kinematics.  \nAcknowledgements. It is a pleasure to thank the anonymous referee for a timely and positive response that helped to improve the presentation of this work. We kindly thank Professors Ian Smail and Adolf N. Witt for valuable comments and discussions. We thank Zhaoxuan Liu for sharing his work and useful comments. M.P. acknowledges financial support from INAF mini-grant 2023 "Galaxy growth and fueling in highz structures." B.L.F. obtained student support through a Faculty Challenge Grant for Increasing Access to Undergraduate Research, and the Arthur L. and Lee G. Herbst Endowment for Innovation and the Science Dean\'s Innovation and Education Fund, both obtained at the University of Arizona. R.A.W., S.H.C. and R.A.J. acknowledge support from NASA JWST Interdisciplinary Scientist grants NAG5-12460, NNX14AN10G and 80NSSC18K0200 from GSFC. We thank the JWST Project at NASA GSFC and JWST Program at NASA HQ for their many-decades long dedication to make the JWST mission a success. We especially thank Tony Roman, the JWST scheduling group and Mission Operations Center sta ff at STScI for their continued dedicated support to get the JWST observations scheduled. This work is based on observations made with the NASA / ESA / CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-03127 for JWST. These observations are associated with JWST programs 1176. We are very grateful to the IRAM sta ff for preparing and carrying out the observations, and for their support and help in reducing the data. Based on observations carried out with the IRAM Interferometer NOEMA under project  \nW21DA. IRAM is supported by INSU / CNRS (France), MPG (Germany) and IGN (Spain). The Herschel spacecraft was designed, built, tested, and launched under a contract to ESA managed by the Herschel / Planck Project team by an industrial consortium under the overall responsibility of the prime contractor Thales Alenia Space (Cannes), and including Astrium (Friedrichshafen) responsible for the payload module and for system testing at spacecraft level, Thales Alenia Space (Turin) responsible for the service module, and Astrium (Toulouse) responsible for the telescope, with in excess of a hundred subcontractors. SPIRE has been developed by a consortium of institutes led by Cardi ff University (UK) and including Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); and Caltech, JPL, NHSC, Univ. Colorado (USA). This development has been supported by national funding agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC, UKSA (UK); and NASA (USA). LOFAR data products were provided by the LOFAR Surveys Key Science project (LSKSP; https://lofar-surveys.org/ ) and were derived from observations with the International LOFAR Telescope (ILT). LOFAR (van Haarlem et al. 2013) is the Low Frequency Array designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, which are owned by various parties (each with their own funding sources), and which are collectively operated by the ILT foundation under a joint scientific policy. The e ff orts of the LSKSP have benefited from funding from the European Research Council, NOVA, NWO, CNRS-INSU, the SURF Co-operative, the UK Science and Technology Funding Council and the Jülich Supercomputing Centre. This work has made use of data from the European Space Agency (ESA) mission Gaia ( https://www. cosmos.esa.int/gaia ), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/ consortium ). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Software: This research made use of astropy, a community developed core Python package for astronomy (Astropy Collaboration et al. 2018), APLpy, an open-source plotting package for Python (Robitaille & Bressert 2012), the IDL Astronomy Library (Landsman 1993), CIGALE (Boquien et al. 2019), starmorph (Rodriguez-Gomez et al. 2019), Galfit (Peng et al. 2002), GILDAS (https: // www.iram.fr / IRAMFR / GILDAS), CARTA (https: // cartavis.org / ; Comrie et al. 2021), and TOPCAT (http: // www.starlink.ac.uk / topcat / ). This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France (Ochsenbein 1996). The original description of the VizieR service was published in Ochsenbein et al. (2000).', 'References': "- Abraham, R. G., van den Bergh, S., Glazebrook, K., et al. 1996, ApJS, 107, 1 Abraham, R. G., van den Bergh, S., & Nair, P. 2003, ApJ, 588, 218 Akins, H. B., Casey, C. M., Allen, N., et al. 2023, ApJ, 956, 61 Alberts, S. & Noble, A. 2022, Universe, 8, 554 Álvarez-Márquez, J., Crespo Gómez, A., Colina, L., et al. 2023, A&A, 671, A105\n- Amvrosiadis, A., Wardlow, J. L., Birkin, J. E., et al. 2023, MNRAS subm., arXiv:2312.08959\n- Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. 2011, PASP, 123, 1218\n- Aravena, M., Decarli, R., Gónzalez-López, J., et al. 2019, ApJ, 882, 136 Armus, L., Heckman, T., & Miley, G. 1987, AJ, 94, 831 Arnouts, S., Le Floc'h, E., Chevallard, J., et al. 2013, A&A, 558, A67 Astropy Collaboration, Price-Whelan, A. M., Sipocz, B. M., et al. 2018, AJ, 156, 123\n- Baes, M., Nersesian, A., Casasola, V., et al. 2020, A&A, 641, A119 Barnes, J. E. & Hernquist, L. 1996, ApJ, 471, 115\n- Barro, G., Pérez-González, P. G., Kocevski, D. D., et al. 2024, ApJ, 963, 128\n- Barrufet, L., Oesch, P. A., Weibel, A., et al. 2023, MNRAS, 522, 449\n- Battisti, A. J., Bagley, M. B., Baronchelli, I., et al. 2022, MNRAS, 513, 4431\n- Baugh, C. M., Lacey, C. G., Frenk, C. S., et al. 2005, MNRAS, 356, 1191\n- Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143  \nBeichman, C. A., Rieke, M., Eisenstein, D., et al. 2012, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 8442, Space Telescopes and Instrumentation 2012: Optical, Infrared, and Millimeter Wave, ed. M. C. Clampin, G. G. Fazio, H. A. MacEwen, & J. Oschmann, Jacobus M., 84422N  \nBershady, M. A., Jangren, A., & Conselice, C. J. 2000, AJ, 119, 2645  \n- Berta, S. 2005, PhD thesis, University of Padua, Department of Physics and Astronomy\n- Berta, S., Lutz, D., Genzel, R., Förster-Schreiber, N. M., & Tacconi, L. J. 2016, A&A, 587, A73  \nBerta, S., Stanley, F., Ismail, D., et al. 2023, A&A, 678, A28 Berta, S., Young, A. J., Cox, P., et al. 2021, A&A, 646, A122 Best, P. N., Kondapally, R., Williams, W. L., et al. 2023, MNRAS, 523, 1729 Bignone, L. A., Tissera, P. B., Sillero, E., et al. 2017, MNRAS, 465, 1106 Birkin, J. E., Weiss, A., Wardlow, J. L., et al. 2021, MNRAS, 501, 3926 Bisigello, L., Gandolfi, G., Grazian, A., et al. 2023, A&A, 676, A76 Blain, A. W. & Longair, M. S. 1993, MNRAS, 264, 509 Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 Boogaard, L. A., Gillman, S., Melinder, J., et al. 2024, ApJ, 969, 27 Boquien, M., Buat, V., Burgarella, D., et al. 2022, A&A, 663, A50 Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, A&A, 622, A103 Bothwell, M. S., Smail, I., Chapman, S. C., et al. 2013, MNRAS, 429, 3047 Boucaud, A., Bocchio, M., Abergel, A., et al. 2016, A&A, 596, A63 Bouché, N., Murphy, M. T., Kacprzak, G. G., et al. 2013, Science, 341, 50 Bournaud, F., Chapon, D., Teyssier, R., et al. 2011, ApJ, 730, 4 Brienza, M., Godfrey, L., Morganti, R., et al. 2017, A&A, 606, A98 Bruzual, A. G. & Charlot, S. 2003, MNRAS, 344, 1000  \n- Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2022, JWST Calibration Pipeline  \nCañameras, R., Nesvadba, N. P. H., Guery, D., et al. 2015, A&A, 581, A105 Calabrò, A., Daddi, E., Fensch, J., et al. 2019, A&A, 632, A98 Calura, F., Pozzi, F., Cresci, G., et al. 2017, MNRAS, 465, 54 Calvi, R., Castignani, G., & Dannerbauer, H. 2023, A&A, 678, A15 Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 Camps-Fariña, A., Sánchez-Blázquez, P., Roca-Fàbrega, S., & Sánchez, S. F. 2023, A&A, 678, A65  \nCardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245  \nCarilli, C. L. & Walter, F. 2013, ARA&A, 51, 105  \nCarter, M., Lazare  \nff  \n, B., Maier, D., et al. 2012, A&A, 538, A89  \nCasey, C. M. 2016, ApJ, 824, 36  \nCasey, C. M., Narayanan, D., & Cooray, A. 2014, Phys. Rep., 541, 45 Castignani, G., Jablonka, P., Combes, F., et al. 2020, A&A, 640, A64 Ceverino, D., Dekel, A., & Bournaud, F. 2010, MNRAS, 404, 2151  \n- Ceverino, D., Sánchez Almeida, J., Muñoz Tuñón, C., et al. 2016, MNRAS, 457, 2605\n- Chabrier, G. 2003, PASP, 115, 763\n- Chapman, S. C., Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772 Charlot, S. & Fall, S. M. 2000, ApJ, 539, 718\n- Chen, C.-C., Gao, Z.-K., Hsu, Q.-N., et al. 2022, ApJ, 939, L7\n- Chen, C.-C., Smail, I., Swinbank, A. M., et al. 2015, ApJ, 799, 194\n- Cheng, C., Huang, J.-S., Smail, I., et al. 2023, ApJ, 942, L19\n- Cheng, C., Yan, H., Huang, J.-S., et al. 2022, ApJ, 936, L19\n- Chiang, Y.-K., Overzier, R., & Gebhardt, K. 2013, ApJ, 779, 127\n- Clements, D. L., Andreani, P., & Chase, S. T. 1993, MNRAS, 261, 299\n- Cochrane, R. K., Anglés-Alcázar, D., Cullen, F., & Hayward, C. C. 2024, ApJ, 961, 37\n- Colina, L., Crespo Gómez, A., Álvarez-Márquez, J., et al. 2023, A&A, 673, L6 Comrie, A., Wang, K.-S., Hsu, S.-C., et al. 2021, CARTA: The Cube Analysis and Rendering Tool for Astronomy\n- Conroy, C. & Gunn, J. E. 2010, ApJ, 712, 833\n- Conselice, C. J. 2003, ApJS, 147, 1\n- Conselice, C. J. 2014, ARA&A, 52, 291\n- Conselice, C. J., Bluck, A. F. L., Ravindranath, S., et al. 2011, MNRAS, 417, 2770\n- Contini, E., De Lucia, G., Hatch, N., Borgani, S., & Kang, X. 2016, MNRAS, 456, 1924  \nCoogan, R. T., Daddi, E., Sargent, M. T., et al. 2018, MNRAS, 479, 703  \nCowie, L. L., Hu, E. M., & Songaila, A. 1995, AJ, 110, 1576  \nCrain, R. A., Schaye, J., Bower, R. G., et al. 2015, MNRAS, 450, 1937  \nCrain, R. A. & van de Voort, F. 2023, ARA&A, 61, 473  \n- Crespo Gómez, A., Colina, L., Álvarez-Márquez, J., et al. 2024, A&A subm., arXiv:2402.18672  \nda Cunha, E., Walter, F., Smail, I. R., et al. 2015, ApJ, 806, 110  \n- Daddi, E., Rich, R. M., Valentino, F., et al. 2022, ApJ, 926, L21\n- Danielson, A. L. R., Swinbank, A. M., Smail, I., et al. 2017, ApJ, 840, 78  \nDannerbauer, H., Lehnert, M. D., Emonts, B., et al. 2017, A&A, 608, A48  \nDarvish, B., Scoville, N. Z., Martin, C., et al. 2020, ApJ, 892, 8  \nDavé, R., Anglés-Alcázar, D., Narayanan, D., et al. 2019, MNRAS, 486, 2827  \nDavé, R., Finlator, K., Oppenheimer, B. D., et al. 2010, MNRAS, 404, 1355  \nDavis, T. A. & Bureau, M. 2016, MNRAS, 457, 272  \nDekel, A., Birnboim, Y., Engel, G., et al. 2009a, Nature, 457, 451  \nDekel, A., Sari, R., & Ceverino, D. 2009b, ApJ, 703, 785  \nDelvecchio, I., Gruppioni, C., Pozzi, F., et al. 2014, MNRAS, 439, 2736  \n- Donevski, D., Lapi, A., Małek, K., et al. 2020, A&A, 644, A144\n- Donnan, F. R., García-Bernete, I., Rigopoulou, D., et al. 2024, MNRAS, 529, 1386  \nDraine, B. T., Aniano, G., Krause, O., et al. 2014, ApJ, 780, 172 Dudzeviˇci¯ut˙e, U., Smail, I., Swinbank, A. M., et al. 2020, MNRAS, 494, 3828 Dunlop, J. S., McLure, R. J., Biggs, A. D., et al. 2017, MNRAS, 466, 861  \n- Ellison, S. L., Mendel, J. T., Scudder, J. M., Patton, D. R., & Palmer, M. J. D. 2013, MNRAS, 430, 3128 Elmegreen, D. M., Elmegreen, B. G., Ferguson, T., & Mullan, B. 2007, ApJ, 663, 734 Elmegreen, D. M., Elmegreen, B. G., Marcus, M. T., et al. 2009, ApJ, 701, 306 Engel, H., Tacconi, L. J., Davies, R. I., et al. 2010, ApJ, 724, 233 Fitzpatrick, E. L. 1999, PASP, 111, 63 Förster Schreiber, N. M., Genzel, R., Bouché, N., et al. 2009, ApJ, 706, 1364 Frayer, D. T., Ivison, R. J., Scoville, N. Z., et al. 1998, ApJ, 506, L7 Frye, B. L., Pascale, M., Foo, N., et al. 2023, ApJ, 952, 81 Frye, B. L., Pascale, M., Pierel, J., et al. 2024, ApJ, 961, 171 Fujimoto, S., Finkelstein, S. L., Burgarella, D., et al. 2023, ApJ, 955, 130 Furtak, L. J., Zitrin, A., Plat, A., et al. 2023, ApJ, 952, 142 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1 Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 Gillman, S., Gullberg, B., Brammer, G., et al. 2023, A&A, 676, A26 Gillman, S., Smail, I., Gullberg, B., et al. 2024, A&A subm., arXiv:2406.03544 Ginolfi, M., Jones, G. C., Béthermin, M., et al. 2020, A&A, 643, A7 Ginolfi, M., Maiolino, R., Nagao, T., et al. 2017, MNRAS, 468, 3468 Girelli, G., Pozzetti, L., Bolzonella, M., et al. 2020, A&A, 634, A135 Glazebrook, K., Ellis, R., Santiago, B., & Gri ffi ths, R. 1995, MNRAS, 275, L19 Gómez-Guijarro, C., Elbaz, D., Xiao, M., et al. 2022, A&A, 658, A43 Gómez-Guijarro, C., Magnelli, B., Elbaz, D., et al. 2023, A&A, 677, A34 Gómez-Guijarro, C., Toft, S., Karim, A., et al. 2018, ApJ, 856, 121 Gordon, K. D., Calzetti, D., & Witt, A. N. 1997, ApJ, 487, 625 Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wol ff , M. J. 2003, ApJ, 594, 279 Gottumukkala, R., Barrufet, L., Oesch, P. A., et al. 2024, MNRAS, 530, 966 Greene, J. E., Labbe, I., Goulding, A. D., et al. 2024, ApJ, 964, 39 Greve, T. R., Bertoldi, F., Smail, I., et al. 2005, MNRAS, 359, 1165 Gri ffi n, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 Gullberg, B., Smail, I., Swinbank, A. M., et al. 2019, MNRAS, 490, 4956 Gürkan, G., Hardcastle, M. J., Smith, D. J. B., et al. 2018, MNRAS, 475, 3010 Hamed, M., Ciesla, L., Béthermin, M., et al. 2021, A&A, 646, A127 Hamed, M., Małek, K., Buat, V., et al. 2023, A&A, 674, A99 Hayward, C. C., Narayanan, D., Kereš, D., et al. 2013, MNRAS, 428, 2529 Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, L7 Herrera-Camus, R., Sturm, E., Graciá-Carpio, J., et al. 2020, A&A, 633, L4 Hezaveh, Y. D., Marrone, D. P., & Holder, G. P. 2012, ApJ, 761, 20 Ho, S. H., Martin, C. L., & Turner, M. L. 2019, ApJ, 875, 54 Hodge, J. A., Carilli, C. L., Walter, F., et al. 2012, ApJ, 760, 11 Hodge, J. A., Smail, I., Walter, F., et al. 2019, ApJ, 876, 130 Hodge, J. A., Swinbank, A. M., Simpson, J. M., et al. 2016, ApJ, 833, 103 Högbom, J. A. 1974, A&AS, 15, 417 Holwerda, B. W., Wu, J. F., Keel, W. C., et al. 2021, ApJ, 914, 142 Hopkins, P. F., Cox, T. J., Hernquist, L., et al. 2013, MNRAS, 430, 1901 Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356 Hsu, Y.-M., Hirashita, H., Lin, Y.-H., Camps, P., & Baes, M. 2023, MNRAS, 519, 2475 Huang, S., Kawabe, R., Kohno, K., et al. 2023, ApJ, 958, L26 Ivison, R. J., Swinbank, A. M., Smail, I., et al. 2013, ApJ, 772, 137 Jiménez-Andrade, E. F., Zavala, J. A., Magnelli, B., et al. 2020, ApJ, 890, 171 Kamieneski, P. S., Frye, B. L., Pascale, M., et al. 2023, ApJ, 955, 91 Kamieneski, P. S., Frye, B. L., Windhorst, R. A., et al. 2024a, ApJ accepted, arXiv:2404.08058 Kamieneski, P. S., Yun, M. S., Harrington, K. C., et al. 2024b, ApJ, 961, 2 Kartaltepe, J. S., Sanders, D. B., Le Floc'h, E., et al. 2010, ApJ, 721, 98 Keel, W. C., Windhorst, R. A., Jansen, R. A., et al. 2023, AJ, 165, 166 Kennicutt, R. C. & Evans, N. J. 2012, ARA&A, 50, 531 Killi, M., Watson, D., Brammer, G., et al. 2023, A&A subm., arXiv:2312.03065 Kneissl, R., Polletta, M. d. C., Martinache, C., et al. 2019, A&A, 625, A96 Knudsen, K. K., van der Werf, P., Franx, M., et al. 2005, ApJ, 632, L9 Kocevski, D. D., Finkelstein, S. L., Barro, G., et al. 2024, ApJ submitted, arXiv:2404.03576 Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, ApJS, 197, 36 Kokorev, V., Caputi, K. I., Greene, J. E., et al. 2024, ApJ, 968, 38 Kokorev, V., Fujimoto, S., Labbe, I., et al. 2023, ApJ, 957, L7 Kriek, M. & Conroy, C. 2013, ApJ, 775, L16 Krumholz, M. R., Burkhart, B., Forbes, J. C., & Crocker, R. M. 2018, MNRAS, 477, 2716 Lagos, C. d. P., da Cunha, E., Robotham, A. S. G., et al. 2020, MNRAS, 499, 1948 Lagos, C. d. P., Robotham, A. S. G., Trayford, J. W., et al. 2019, MNRAS, 489, 4196 Lagos, C. d. P., Tobar, R. J., Robotham, A. S. G., et al. 2018, MNRAS, 481, 3573 Landsman, W. B. 1993, in Astronomical Society of the Pacific Conference Series, Vol. 52, Astronomical Data Analysis Software and Systems II, ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 246\n- Le Bail, A., Daddi, E., Elbaz, D., et al. 2024, A&A accepted, arXiv:2307.07599 Liao, C.-L., Chen, C.-C., Wang, W.-H., et al. 2024, ApJ, 961, 226  \n- Lilly, S. J., Eales, S. A., Gear, W. K. P., et al. 1999, ApJ, 518, 641\n- Ling, C., Sun, B., Cheng, C., et al. 2024, ApJ, 969, L28\n- Lira, P., Gonzalez-Corvalan, V., Ward, M., & Hoyer, S. 2008, MNRAS, 384, 316\n- Liu, D., Schinnerer, E., Groves, B., et al. 2019, ApJ, 887, 235\n- Liu, F. S., Jiang, D., Faber, S. M., et al. 2017, ApJ, 844, L2\n- Liu, F. S., Jiang, D., Guo, Y., et al. 2016, ApJ, 822, L25\n- Liu, Z., Silverman, J. D., Daddi, E., et al. 2024, ApJ, 968, 15\n- Lo Faro, B., Buat, V., Roehlly, Y., et al. 2017, MNRAS, 472, 1372  \nLorenz, B., Kriek, M., Shapley, A. E., et al. 2023, ApJ, 951, 29  \n- Lotz, J. M., Davis, M., Faber, S. M., et al. 2008a, ApJ, 672, 177  \nLotz, J. M., Jonsson, P., Cox, T. J., & Primack, J. R. 2008b, MNRAS, 391, 1137  \nLotz, J. M., Primack, J., & Madau, P. 2004, AJ, 128, 163  \n- Madau, P. & Dickinson, M. 2014, ARA&A, 52, 415  \nMagdis, G. E., Rigopoulou, D., Daddi, E., et al. 2017, A&A, 603, A93  \n- Marchuk, A. A., Chugunov, I. V., Gontcharov, G. A., et al. 2024, MNRAS, 528, 1276  \nMarkov, V., Gallerani, S., Pallottini, A., et al. 2023, A&A, 679, A12  \nMartinache, C., Rettura, A., Dole, H., et al. 2018, A&A, 620, A198  \n- Matthee, J., Naidu, R. P., Brammer, G., et al. 2024, ApJ, 963, 129\n- Mayya, Y. D., Bressan, A., Rodríguez, M., Valdes, J. R., & Chavez, M. 2004, ApJ, 600, 188\n- McKay, S. J., Barger, A. J., & Cowie, L. L. 2024, ApJ, 962, 128\n- Meyer, R. A., Barrufet, L., Boogaard, L. A., et al. 2024, A&A, 681, L3  \nMichałowski, M. J., Hayward, C. C., Dunlop, J. S., et al. 2014, A&A, 571, A75  \n- Mihos, J. C. & Hernquist, L. 1996, ApJ, 464, 641\n- Miley, G. & De Breuck, C. 2008, A&A Rev., 15, 67\n- Miller, T. B., Chapman, S. C., Aravena, M., et al. 2018, Nature, 556, 469\n- Miller, T. B., van Dokkum, P., & Mowla, L. 2023, ApJ, 945, 155\n- Miller, T. B., Whitaker, K. E., Nelson, E. J., et al. 2022, ApJ, 941, L37\n- Mok, A., Chandar, R., & Fall, S. M. 2020, ApJ, 893, 135\n- Montaña, A., Zavala, J. A., Aretxaga, I., et al. 2021, MNRAS, 505, 5260  \nMoreno, J., Torrey, P., Ellison, S. L., et al. 2019, MNRAS, 485, 1320  \nNarayanan, D., Krumholz, M. R., Ostriker, E. C., & Hernquist, L. 2012, MN-  \nRAS, 421, 3127  \nNarayanan, D., Turk, M., Feldmann, R., et al. 2015, Nature, 525, 496 Negrello, M., Hopwood, R., De Zotti, G., et al. 2010, Science, 330, 800 Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624 Nelson, E. J., Suess, K. A., Bezanson, R., et al. 2023, ApJ, 948, L18  \n- Nelson, E. J., van Dokkum, P. G., Förster Schreiber, N. M., et al. 2016, ApJ, 828, 27\n- Nowak, M. A., Neilsen, J., Marko ff , S. B., et al. 2012, ApJ, 759, 95\n- Nozawa, T. 2016, Planet. Space Sci., 133, 36\n- Ochsenbein, F. 1996, The VizieR database of astronomical catalogues  \nOchsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23  \n- Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713\n- Ortiz, Rafael, I., Windhorst, R. A., Cohen, S. H., et al. 2024, ApJ subm., arXiv:2404.10709\n- Osborne, C. & Salim, S. 2024, ApJ, 962, 59\n- Oser, L., Ostriker, J. P., Naab, T., Johansson, P. H., & Burkert, A. 2010, ApJ, 725, 2312\n- Oteo, I., Ivison, R. J., Dunne, L., et al. 2018, ApJ, 856, 72\n- Pallottini, A., Ferrara, A., Gallerani, S., Salvadori, S., & D'Odorico, V. 2014, MNRAS, 440, 2498\n- Pan, H.-A., Lin, L., Hsieh, B.-C., et al. 2019, ApJ, 881, 119\n- Papadopoulos, P. P., van der Werf, P., Isaak, K., & Xilouris, E. M. 2010, ApJ, 715, 775\n- Pascale, M., Frye, B. L., Dai, L., et al. 2022, ApJ, 932, 85\n- Pawlik, M. M., Wild, V., Walcher, C. J., et al. 2016, MNRAS, 456, 3032\n- Peng, B., Vishwas, A., Stacey, G., et al. 2023, ApJ, 944, L36\n- Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266  \nPeng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097  \nPensabene, A., Cantalupo, S., Cicone, C., et al. 2024, A&A, 684, A119  \n- Pérez-González, P. G., Barro, G., Annunziatella, M., et al. 2023, ApJ, 946, L16\n- Pérez-González, P. G., Barro, G., Rieke, G. H., et al. 2024, ApJ, 968, 4\n- Pérez-Martínez, J. M., Dannerbauer, H., Kodama, T., et al. 2023, MNRAS, 518, 1707\n- Pérez-Martínez, J. M., Kodama, T., Koyama, Y., et al. 2024, MNRAS, 527, 10221\n- Perrin, M. D., Sivaramakrishnan, A., Lajoie, C.-P., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, ed. J. Oschmann, Jacobus M., M. Clampin, G. G. Fazio, & H. A. MacEwen, 91433X  \nPillepich, A., Springel, V., Nelson, D., et al. 2018, MNRAS, 473, 4077 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 596, A100 Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&A, 641, A6 Planck Collaboration, Aghanim, N., Altieri, B., et al. 2015, A&A, 582, A30 Polletta, M., Dole, H., Martinache, C., et al. 2022, A&A, 662, A85 Polletta, M., Soucail, G., Dole, H., et al. 2021, A&A, 654, A121 Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81  \n- Polletta, M., Wilkes, B. J., Siana, B., et al. 2006, ApJ, 642, 673 Popesso, P., Concas, A., Cresci, G., et al. 2023, MNRAS, 519, 1526 Popping, G., Somerville, R. S., & Galametz, M. 2017, MNRAS, 471, 3152 Predehl, P. & Schmitt, J. H. M. M. 1995, A&A, 293, 889 Price, S. H., Suess, K. A., Williams, C. C., et al. 2023, ApJ subm., arXiv:2310.02500\n- Ren, J., Liu, F. S., Li, N., et al. 2024, ApJ, 969, 4\n- Rieke, M. J., Baum, S. A., Beichman, C. A., et al. 2003, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 4850, IR Space Telescopes and Instruments, ed. J. C. Mather, 478-485\n- Rieke, M. J., Kelly, D., & Horner, S. 2005, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5904, Cryogenic Optical Systems and Instruments XI, ed. J. B. Heaney & L. G. Burriesci, 1-8  \nRizzo, F., Kohandel, M., Pallottini, A., et al. 2022, A&A, 667, A5  \nRizzo, F., Vegetti, S., Powell, D., et al. 2020, Nature, 584, 201  \n- Robertson, C., Holwerda, B. W., Young, J., et al. 2024, AJ, 167, 263\n- Robitaille, T. & Bressert, E. 2012, APLpy: Astronomical Plotting Library in Python\n- Robotham, A., Davies, L., Driver, S., et al. 2018, ProFound: Source Extraction and Application to Modern Survey Data\n- Robotham, A. S. G., D'Silva, J. C. J., Windhorst, R. A., et al. 2023, PASP, 135, 085003\n- Rodighiero, G., Bisigello, L., Iani, E., et al. 2023, MNRAS, 518, L19\n- Rodighiero, G., Daddi, E., Baronchelli, I., et al. 2011, ApJ, 739, L40  \nRodriguez-Gomez, V., Snyder, G. F., Lotz, J. M., et al. 2019, MNRAS, 483, 4140  \nSakamoto, K., Aalto, S., Combes, F., Evans, A., & Peck, A. 2014, ApJ, 797, 90  \nSalim, S. & Narayanan, D. 2020, ARA&A, 58, 529  \n- Salmon, B., Papovich, C., Long, J., et al. 2016, ApJ, 827, 20\n- Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749\n- Scoville, N., Sheth, K., Aussel, H., et al. 2016, ApJ, 820, 83  \nScoville, N., Sheth, K., Walter, F., et al. 2015, ApJ, 800, 70  \n- Scudder, J. M., Oliver, S., Hurley, P. D., et al. 2016, MNRAS, 460, 1119\n- Shahidi, A., Mobasher, B., Nayyeri, H., et al. 2020, ApJ, 897, 44\n- Shimwell, T. W., Röttgering, H. J. A., Best, P. N., et al. 2017, A&A, 598, A104\n- Simpson, J. M., Smail, I., Swinbank, A. M., et al. 2017, ApJ, 839, 58\n- Simpson, J. M., Swinbank, A. M., Smail, I., et al. 2014, ApJ, 788, 125\n- Skalidis, R., Goldsmith, P. F., Hopkins, P. F., & Ponnada, S. B. 2024, A&A, 682, A161\n- Smail, I., Dudzeviˇci¯ut˙e, U., Gurwell, M., et al. 2023, ApJ, 958, 36\n- Smail, I., Ivison, R. J., & Blain, A. W. 1997, ApJ, 490, L5\n- Smith, D. J. B., Haskell, P., Gürkan, G., et al. 2021, A&A, 648, A6\n- Snyder, G. F., Torrey, P., Lotz, J. M., et al. 2015, MNRAS, 454, 1886  \nSorba, R. & Sawicki, M. 2018, MNRAS, 476, 1532  \n- Sparre, M., Whittingham, J., Damle, M., et al. 2022, MNRAS, 509, 2720\n- Springel, V. & Hernquist, L. 2005, ApJ, 622, L9  \nStecher, T. P. 1965, ApJ, 142, 1683  \n- Sugahara, Y., Álvarez-Márquez, J., Hashimoto, T., et al. 2024, ApJ subm., arXiv:2403.17133\n- Sun, F., Helton, J. M., Egami, E., et al. 2024, ApJ, 961, 69\n- Tadaki, K., Iono, D., Yun, M. S., et al. 2018, Nature, 560, 613\n- Tadaki, K.-i., Belli, S., Burkert, A., et al. 2020, ApJ, 901, 74\n- Thorp, M. D., Ellison, S. L., Pan, H.-A., et al. 2022, MNRAS, 516, 1462\n- Thorp, M. D., Ellison, S. L., Simard, L., Sánchez, S. F., & Antonio, B. 2019, MNRAS, 482, L55\n- Tohill, C., Ferreira, L., Conselice, C. J., Bamford, S. P., & Ferrari, F. 2021, ApJ, 916, 4  \nTremmel, M., Karcher, M., Governato, F., et al. 2017, MNRAS, 470, 1121 Valentino, F., Daddi, E., Finoguenov, A., et al. 2016, ApJ, 829, 53 van den Bergh, S., Abraham, R. G., Ellis, R. S., et al. 1996, AJ, 112, 359 van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2 Vidal-García, A., Falgarone, E., Arrigoni Battaia, F., et al. 2021, MNRAS, 506, 2551  \n- Wang, B., de Graa ff , A., Davies, R. L., et al. 2024a, ApJ subm.,  \narXiv:2403.02304 Wang, B., Leja, J., Labbé, I., et al. 2024b, ApJS, 270, 12 Wang, S. & Chen, X. 2019, ApJ, 877, 116 Wang, S., Liu, J., Qiu, Y., et al. 2016, ApJS, 224, 40 Wang, W., Faber, S. M., Liu, F. S., et al. 2017, MNRAS, 469, 4063 Ward, E., de la Vega, A., Mobasher, B., et al. 2024, ApJ, 962, 176 Westmeier, T., Braun, R., & Thilker, D. 2005, A&A, 436, 101 Whitney, A., Ferreira, L., Conselice, C. J., & Duncan, K. 2021, ApJ, 919, 139 Williams, C. C., Alberts, S., Ji, Z., et al. 2024, ApJ, 968, 34 Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & Labbé, I. 2009, ApJ, 691, 1879 Windhorst, R. A., Cohen, S. H., Jansen, R. A., et al. 2023, AJ, 165, 13 Windhorst, R. A., Hathi, N. P., Cohen, S. H., et al. 2008, Advances in Space Research, 41, 1965 Witstok, J., Shivaei, I., Smit, R., et al. 2023, Nature, 621, 267 Witt, A. N. & Gordon, K. D. 2000, ApJ, 528, 799  \nWitt, A. N., Thronson, Harley A., J., & Capuano, John M., J. 1992, ApJ, 393,  \n611 Wu, Y., Cai, Z., Sun, F., et al. 2023, ApJ, 942, L1 Wuyts, S., Cox, T. J., Hayward, C. C., et al. 2010, ApJ, 722, 1666 Xiao, M., Oesch, P., Elbaz, D., et al. 2023, Nature subm., arXiv:2309.02492 Xu, C., Wang, T., Gu, Q., et al. 2023, ApJ, 951, L21 Zavala, J. A., Buat, V., Casey, C. M., et al. 2023, ApJ, 943, L9 Zhou, L., Wang, T., Daddi, E., et al. 2024, A&A, 684, A196  \n- 1 INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, Via A. Corti 12, I-20133 Milano, Italy\n- 2 Steward Observatory, University of Arizona, 933 N Cherry Ave, Tucson, AZ, 85721-0009, USA\n- 3 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA\n- 4 Institut de Radioastronomie Millimétrique, 300 rue de la piscine, F38406 Saint-Martin-d'Hères, France\n- 5 European Southern Observatory, ESO Vitacura, Alonso de Córdova 3107, Vitacura, 19001, Casilla, Santiago, Chile\n- 6 Joint ALMA Observatory, Alonso de Córdova 3107, Vitacura 7630355, Santiago, Chile\n- 7 Université Paris-Saclay, CNRS, Institut d'Astrophysique Spatiale, 91405, Orsay, France\n- 8 School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287-1404, USA\n- 9 Univ Lyon, Univ Lyon1, Ens de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F-69230 Saint-Genis-Laval, France\n- 10 International Centre for Radio Astronomy Research (ICRAR) and the International Space Centre (ISC), The University of Western Australia, M468, 35 Stirling Highway, Crawley, WA 6009, Australia\n- 11 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia\n- 12 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA\n- 13 Association of Universities for Research in Astronomy (AURA) for the European Space Agency (ESA), STScI, Baltimore, MD 21218, USA\n- 14 Center for Astrophysical Sciences, Department of Physics and Astronomy, The Johns Hopkins University, 3400 N Charles St. Baltimore, MD 21218, USA\n- 15 Jodrell Bank Centre for Astrophysics, Alan Turing Building, University of Manchester, Oxford Road, Manchester M13 9PL, UK\n- 16 National Research Council of Canada, Herzberg Astronomy & Astrophysics Research Centre, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada\n- 17 INAF-Osservatorio Astronomico di Trieste, Via Bazzoni 2, 34124 Trieste, Italy\n- 18 Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA\n- 19 Chinese Academy of Sciences South America Center for Astronomy, National Astronomical Observatories, CAS, Beijing 100101, China\n- 20 CAS Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China\n- 21 Department of Physics, Montana State University, Bozeman, MT 59717, USA\n- 22 Department of Physics and Astronomy, University of Louisville, 102 Natural Science Building, Louisville KY 40292, USA\n- 23 Dept. of Physics and Astronomy, University of Alabama, Box 870324, Tuscaloosa, AL 35404, USA\n- 24 Department of Astronomy, University of California, 501 Campbell Hall #3411, Berkeley, CA 94720, USA\n- 25 Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands\n- 26 Department of Astronomy, University of Massachusetts at Amherst, Amherst, MA 01003, USA  \nTable A.1: Masked sources in 01a.Notes. ( a ) Radius of circular aperture or major and minor semi-axis of ellipse in arcsec centered on coordinates α and δ .", "Appendix A: Neighbors' contamination in 01a": "Three sources, designated A1, A2, and A3, lie within 0 '' . 5-1 '' . 5 of the center of 01a. These sources seem to be located (Fig. 3) along an arc or spiral arm extending from the major axis of A3 to the center of 01a, suggesting that the A sources might be the result of an on-going interaction with 01a and among themselves. If so, they would be at the same redshift, and at least A1 and A2 would be part of 01a, but their bluer colors suggest that they are distinct sources. To assess this hypothesis, we built their SEDs by measuring all the flux within specific apertures and used CIGALE to estimate their photometric redshifts. The centers and sizes of the three photometric apertures are listed in Table A.1 and shown in the left panel of Fig. 3.  \nWe fit the galaxies' SEDs with the redshifts unconstrained (Table A.1) and also fixed to the spectroscopic redshift of 01a. The two bestfit models are shown in Fig. A.1. The reduced χ 2 of A1 and A2 are 2.3-2.5 times smaller at the best photometric redshifts than at z = z 01a, but z = z 01a yields also a good fit for A2 and A3 ( χ 2 ν ≤ 0 . 23). The χ 2 ν for A1 at z = z 01a makes it a foreground source. Source A3 is 1 '' . 5 from 01a, and if it were at the same redshift of 01a, it would be a companion galaxy with a mass ratio of ∼ 1:22. The association between A1 and 01a is the most complex because A1 could be a super star cluster (SSC) or a star-forming region within 01a. The F444W flux density of A1 is ∼ 11% of the whole 01a galaxy. At z = z 01a, A1's stellar mass would be ≲ 3 × 10 10 M ⊙ , and its SFR would be ∼ 30 M ⊙ yr -1 . A1's size, brightness, SFR, and M star are extremely large for a SSC or a star-forming region (see e.g., Mok et al. 2020), implying that it is more likely a distinct galaxy rather than part of 01a. We have assumed that A1 and A2 are in the foreground of 01a, and we have masked them in the analysis of 01a, but we cannot confidently rule out that these sources might be at 01a's redshift and undergoing an interaction. Source A3 was also masked because it is a distinct galaxy, although likely interacting with 01a. To fully capture the properties of 01a, it is important to establish whether A1 and A2 are part of it or distinct galaxies along the line of sight or interacting with it, although this information would not change the results reported in this work. It is nonetheless important to estimate a potential contribution to the submm and CO emission that we might have wrongly assigned to 01a from these distinct sources. The bestfit SED models and the recently obtained ALMA data rule out any significant contribution to 01a' submm and CO emission from these two sources.", 'Appendix B: Comparison between total and pixel-based parameters': "Physical parameters derived from resolved SED studies do not necessarily recover the values obtained from integrated SED fitting (e.g., Sorba & Sawicki 2018). This is a concern here because we could use only the eight NIRCam data points in the  \nFig. A.1: SED of the sources within 1.5 '' of source 01a and considered distinct sources (filled red circles; A1 in the top panel, A2 in the middle panel, and A3 in the bottom panel). In each panel, the best fit model, associated photometric redshift, and reducedχ 2 are shown in blue, and the model at the same redshift of 01a and corresponding reducedχ 2 are shown in magenta.  \n<!-- image -->  \nresolved SED, while we included constraints from the submm, mm, and radio data when fitting the integrated SED. To limit this potential problem, we adopted the same components and model parameters, with the exception of the radio component, in fitting both the resolved and the integrated SEDs. Fig. B.1 compares the main parameters analyzed through resolved SED fitting obtained from the two approaches. The resolved SED fitting yields stellar masses that are, on average, a factor of 2.7 (0.4 dex) higher than those derived from the integrated SED. Such a bias had been previously reported and estimated to increase with the specific SFRs (sSFRs; Sorba & Sawicki 2018). Because of the large sSFRs of our sources, they might thus su ff er by a large bias in the global stellar mass estimates. Interestingly, the largest discrep-  \nancy is observed in 01a, the source with the highest sSFR of the three, whereas the agreement between the two mass estimates is excellent in 01c that has a much lower sSFR. This comparison suggests that the global mass estimates of 01a and 01b might need to be multiplied by a factor of 4.9 and 2.3, respectively. These corrections would not change the overall results obtained for these sources.  \nContrary to the stellar mass, the sum of the SFRs from the resolved SED fitting is on average 0.6 times (0.2 dex) lower than the SFR from the unresolved SED. The two estimates agree for 01c, and similar o ff sets are measured for 01a and 01b. As previously stated, SFRs derived from the resolved SED fitting are highly uncertain due to the lack of constraints at submm / mm / radio wavelengths that are instead available in the unresolved SED. The resolved maps of SFR are thus to be considered with caution, but it should be safe to analyze the SFR distribution across each galaxy assuming that the o ff set between resolved and unresolved SFRs in each galaxy is systematic.  \nAs a final check we also compared the extinctions and the mmflux densities from the two approaches. These quantities are both consistent implying that the distribution of extinction in the resolved maps can be safely analyzed.", 'Appendix C: Morphological analysis': 'Fig. C.1 shows the input masked F444W images and the output images obtained with statmorph along with key morphological parameters. Table 7 gives the full list of derived morphological parameters. Fig. C.2 compares the statmorph parameters with the classification regions defined by Lotz et al. (2008a); Conselice (2003); Bershady et al. (2000). This morphological classification is based on nearby galaxies and on images at 5500Å in the rest frame, and therefore it might not work as well for highz sources because of the di ff erent rest-frame wavelength (i.e., ∼ 1.3 µ m), spatial sampling and surface brightness limits. Some works find that galaxies become more compact at increasing redshifts, and accordingly the average concentration C increases from 3.0 to 3.7 and the asymmetry A from 0.22 to 0.30 from z ≃ 1 . 25 to z ≃ 2 . 75 (Whitney et al. 2021; Tohill et al. 2021), regardless of stellar mass. The separation between disks and spheroids in the A -C plane becomes less evident with increasing redshift and disappears for z ≳ 2, which corresponds to our sources. Other studies (Ren et al. 2024), using a large sample of simulated galaxies at 0 . 5 < z < 3, find no significant variations in the morphological indicators for SFGs when considering rest-frame wavelengths > 1 . 1 µ m. Since we carry out this analysis at rest-frame λ ∼ 1 . 3 µ m, the standard criteria can be safely applied.  \nGalfit was run on all seven wide NIRCam bands, but only the LW images, F277W, F356W, and F444W, gave valid results for all three sources. For 01a, we modeled simultaneously the bright nearby clump A1, and the nearby galaxy A3 (App. A). Fig. C.3, shows the results. For 01a, the residuals show several clumps that seem to be aligned along three spiral arms, and some di ff use emission to the north. The residual map of 01b shows emission in the center, an arc or spiral arm to the west, a clump to the southeast, and a straight feature, probably another galaxy, to the east. Source 01c contains four clumps along the galaxy major axis. A faint feature, similar to a tidal tail, seems to connect the southern edge of the galaxy with a nearby, faint source to the southeast.  \nFig. B.1: Comparison between the sum of the resolved parameters and the parameters derived from fitting the integrated SEDs. Filled red, green, and blue points represent, respectively, sources 01a, 01b, and 01c. The panels show, from top to bottom, the stellar mass, the instantaneous SFR, the extinction, and the 3 mm flux density. The parameter mean( AV ) is the weighted mean of the pixel-derived values. The leftpointing arrows represent 5 σ upper limits (Table 3). The dashed black line shows the 1:1 relation, and the dotted magenta line the average ratio for the three galaxies (2.7, 0.6, 1.0, and 0.8 from top to bottom). The resolved sums include only pixels with a value above 1-2 σ as was done in building the resolved maps shown in Fig. 6.  \n<!-- image -->', 'Polletta et al.: JWST starburst galaxies': "Fig. C.1: Input and output F444W images from statmorph . Pairs of rows from top to bottom show the three targets 01a, 01b, and 01c. Images within each set of six are labeled, and the 01a images show the areas masked because of sources A1 and A2. All images are 6 '' on a side.  \n<!-- image -->  \nFig. C.2: Morphological parameters of the three CO emitters (filled red circles). Dashed lines and shaded regions show the division between di ff erent morphological types as defined by Lotz et al. (2008a) in the GiniM 20 plane (left panel), by Conselice (2003) in the S -C plane (middle panel), and by Bershady et al. (2000) in the A -C plane (right panel). The blue diamond represents the median value measured on JWST NIR images for a sample of SMGs (Gillman et al. 2023) with error bars indicating the sample standard deviation of each parameter.  \n<!-- image -->  \nFig. C.3: Galfit input image (6 '' × 6 '' ), Sérsic model, and residuals in three bands, F277W (top), F356W (central), and F444W (bottom) and for sources 01a (left), 01b (middle), and 01c (right). The source center is marked by a red cross.  \n<!-- image -->"}