We present JWST NIRCam imaging targeting 13 zsim3 infraredluminous Lrm IRsim5times1012Lodot galaxies from the ALESS survey with uniquely deep highresolution 0.080.16 ALMA 870mum imaging. The 2.04.4mum observed frame NIRCam imaging reveals the restframe nearinfrared stellar emission in these submillimeterselected galaxies SMGs at the same subkpc resolution as the 870mum dust continuum. The newly revealed stellar morphologies show striking similarities with the dust continuum morphologies at 870mum with the centers and position angles agreeing for most sources clearly illustrating that the spatial offsets reported previously between the 870mum and HST morphologies were due to strong differential dust obscuration. The F444W sizes are 78pm21 larger than those measured at 870mum in contrast to recent results from hydrodynamical simulations that predict larger 870mum sizes. We report evidence for significant dust obscuration in F444W for the highestredshift sources emphasizing the importance of longerwavelength MIRI imaging. The majority of the sources show evidence that they are undergoing mergersinteractions including tidal tailsplumes  some of which are also detected at 870mum. We find a clear correlation between NIRCam colors and 870mum surface brightness on sim1 kpc scales indicating that the galaxies are primarily red due to dust  not stellar age  and we show that the dust structure on simkpcscales is broadly similar to that in nearby galaxies. Finally we find no strong stellar bars in the restframe nearinfrared suggesting the extended barlike features seen at 870mum are highly obscured andor gasdominated structures that are likely early precursors to significant bulge growth.
['Astrophysics - Astrophysics of Galaxies']
ALESSJWST Joint subkiloparsec JWST and ALMA imaging of zsim3 submillimeter galaxies reveals heavily obscured bulge formation events
{'ALESS-JWST: Joint (sub-)kiloparsec JWST and ALMA imaging of z ∼ 3 submillimeter galaxies reveals heavily obscured bulge formation events': "<!-- image -->  \n1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, the Netherlands ∗  \n2 International Centre for Radio Astronomy Research, University of Western Australia, 35 Stirling Hwy, Crawley 26WA 6009, Australia 3 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia 4 European Space Agency (ESA), ESA Office, Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 5 Centre for Extragalactic Astronomy, Department of Physics, Durham University, South Road, Durham DH1 3LE, UK 6 NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD, USA 7 Academia Sinica Institute of Astronomy and Astrophysics (ASIAA), No. 1, Sec. 4, Roosevelt Road, Taipei 106216, Taiwan 8 Max-Planck Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 9 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 10 Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia 11 Department of Astronomy and Astrophysics, 525 Davey Lab, The Pennsylvania State University, University Park, PA 16802, USA 12 Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 13 Department of Physics, 104 Davey Laboratory, The Pennsylvania State University, University Park, PA 16802, USA 14 German Aerospace Center (DLR), Institute of Communications and Navigation, Wessling, Germany 15 European Southern Observatory (ESO), Karl-Schwarzschild-Straße 2, 85748 Garching bei M /dieresis.ts1 unchen, Germany 16 Department of Physics and Atmospheric Science, Dalhousie University, 6310 Coburg Road, B3H 4R2, Halifax, Canada 17 National Research Council, Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, V9E 2E7, Canada 18 Department of Physics and Astronomy, University of British Columbia, 6225 Agricultural Road, Vancouver, V6T 1Z1, Canada 19 Eureka Scientific Inc, Oakland, CA 94602, USA 20 Sorbonne Universit'e, UPMC Universit'e Paris 6 and CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98bis boulevard Arago, 75014 Paris, France 21 Instituto de Astrof'ısica de Canarias, V'ıa L'actea, 39020 La Laguna (Tenerife), Spain 22 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, via Gobetti 93/3, I-40129, Bologna, Italy 23 Cosmic Dawn Center (DAWN), Denmark 24 DTU Space, Technical University of Denmark, Elektrovej, Building 328, 2800, Kgs. Lyngby, Denmark 25 Dept. of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom 26 Department of Space, Earth and Environment, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden 27 Max-Planck-Institut fur Radioastronomie, Auf dem Hugel 69, D-53121 Bonn, Germany 28  \nFaculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands  \n29 SRON - Netherlands Institute for Space Research, Niels Bohrweg 4, 2333 CA Leiden, The Netherlands 30 Department of Physics, Lancaster University, Lancaster, LA1 4YB, UK", 'ABSTRACT': "We present JWST NIRCam imaging targeting 13 z ∼ 3 infrared-luminous ( L IR ∼ 5 × 10 12 L ⊙ ) galaxies from the ALESS survey with uniquely deep, high-resolution (0.08 '' -0.16 '' ) ALMA 870 µ m imaging. The 2.0-4.4 µ m (observed frame) NIRCam imaging reveals the rest-frame near-infrared stellar emission in these submillimeter-selected galaxies (SMGs) at the same (sub-)kpc resolution as the 870 µ m dust continuum. The newly revealed stellar morphologies show striking similarities with the dust continuum morphologies at 870 µ m, with the centers and position angles agreeing for most sources, clearly illustrating that the spatial offsets reported previously between the 870 µ mand HST morphologies were due to strong differential dust obscuration. The F444W sizes are 78 ± 21% larger than those measured at 870 µ m, in contrast to recent results from hydrodynamical simulations that predict larger 870 µ m sizes. We report evidence for significant dust obscuration in F444W for the highest-redshift sources,  \nemphasizing the importance of longer-wavelength MIRI imaging. The majority of the sources show evidence that they are undergoing mergers/interactions, including tidal tails/plumes-some of which are also detected at 870 µ m. We find a clear correlation between NIRCam colors and 870 µ m surface brightness on ∼ 1 kpc scales, indicating that the galaxies are primarily red due to dust-not stellar age-and we show that the dust structure on ∼ kpc-scales is broadly similar to that in nearby galaxies. Finally, we find no strong stellar bars in the rest-frame near-infrared, suggesting the extended barlike features seen at 870 µ m are highly obscured and/or gas-dominated structures that are likely early precursors to significant bulge growth.  \nKeywords: galaxies: evolution - galaxies: star formation", '1. INTRODUCTION': "With the advent of the first infrared (IR) sky surveys, including by IRAS (Neugebauer et al. 1984) and COBE (Mather 1994), a basic fact of the Universe was revealed: Roughly half of the star formation that has occurred over cosmic time took place shrouded behind cosmic dust (e.g., Puget et al. 1996; Hauser et al. 1998). Thanks to the decades of dedicated follow-up that subsequently occurred in order to systematically resolve the IR/sub-millimeter background, we now know that this dust-obscured star formation was as important as unobscured star formation at least out to z ∼ 4 (e.g., Dudzeviˇci¯ut˙e et al. 2020; Bouwens et al. 2020; Zavala et al. 2021), and perhaps even beyond (e.g., Algera et al. 2023). The galaxies with the most significant dustobscured star formation rates are the z ∼ 2 -5 galaxies discovered in the earliest single-dish submillimeter surveys ('SMGs'; e.g., Smail et al. 1997; Hughes et al. 1998; Barger et al. 1998; Eales et al. 1999). Canonically, SMGs are thought to be tied to the evolution of QSOs (e.g., Sanders & Mirabel 1996; Hopkins et al. 2006) and ultimately to the build-up of massive elliptical galaxies (e.g., Smail et al. 2004; Swinbank et al. 2006; Cimatti et al. 2008; van Dokkum et al. 2008; Simpson et al. 2014; Toft et al. 2014, and see Casey et al. 2014 for a review). Understanding the physical processes that drove their intense dusty star-formation is thus crucial for understanding massive galaxy evolution.  \nPrior to the launch of the James Webb Space Telescope (JWST), characterizing the stellar counterparts of SMGs was notoriously difficult. Studies using, e.g., deep HST H-band imaging report conflicting results on the prevalence of merger activity (e.g., Swinbank et al. 2010; Targett et al. 2013; Chen et al. 2015). However, these dust-rich galaxies are particularly prone to structured dust extinction (e.g., Simpson et al. 2017; Cochrane et al. 2019), suggesting their rest-frame UV/optical clas-  \nsifications may be unreliable (e.g., Lang et al. 2019; Popping et al. 2022). Moreover, such studies also report high fractions of 'NIR-faint' sources ( ∼ 20% at m H ∼ 27.8 mag depth; Chen et al. 2015), presumably due to their high levels of dust attenuation (median global A V ∼ 23 mag; da Cunha et al. 2015; Dudzeviˇci¯ut˙e et al. 2020), precluding reliable morphological classification. Finally, even for the sources detected with, e.g., HST, their stellar masses typically carry large uncertainties due to the well-known degeneracy between stellar age and dust attenuation (e.g., Hainline et al. 2011; Michaglyph[suppress]lowski et al. 2014; Dudzeviˇci¯ut˙e et al. 2020). These observational challenges have contributed to the lack of theoretical consensus on the nature of SMGs and the physical mechanism(s) by which their dusty star formation was triggered (e.g., Baugh et al. 2005; Narayanan et al. 2010, 2015; Hayward et al. 2013b,a; McAlpine et al. 2019; Lovell et al. 2021).  \nIn the last decade, the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) has provided critical new insight into dust-obscured star formation at high redshift (e.g., see Hodge & da Cunha 2020 for a review). In particular, ALMA (along with the Submillimeter Array, SMA, and the Northern Extended Millimeter Array, NOEMA) have allowed large samples of the classically single-dish-selected SMGs to be located with interferometric accuracy (e.g., Younger et al. 2008, 2009; Barger et al. 2012; Smolˇci'c et al. 2012; Hodge et al. 2013; Simpson et al. 2014; Spilker et al. 2016; Hill et al. 2018; Stach et al. 2019) and resolved directly in their (sub-)millimeter continuum emission (e.g., Simpson et al. 2015; Ikarashi et al. 2015; Hodge et al. 2016; Gullberg et al. 2019). These studies have revealed that the (sub-)millimeter continuum emission in these galaxies is relatively compact and has disk-like morphology (i.e., with a S'ersic index n ∼ 1; e.g., Simpson et al. 2015; Hodge et al. 2016; Gullberg et al. 2019) and often appears spatially offset from the rest-frame UV/optical continuum (e.g., Hodge et al. 2012; Chen et al. 2015; Calistro Rivera et al. 2018), which some  \nstudies have interpreted as being due to real physical offsets between the location of the dusty star formation and the bulk of the stellar mass (Chen et al. 2015). The highest-resolution studies have further resolved the (sub-)millimeter continuum-traced dust disks in small sub-samples into distinct substructures on ∼ 0.5-1 kpc scales, uncovering features with morphologies suggestive of bars, rings, and (tidally-induced) spiral arms (Hodge et al. 2019, and see also Gullberg et al. 2019). Unfortunately, despite the depth of the available HST imaging, no meaningful comparison could be made to the stellar structure due to the faintness of the sources in the rest-frame UV/optical.  \nRecently, the launch of JWST has shed new light on the nature of submillimeter-selected sources by revealing their rest-frame near-infrared morphologies on (sub-)kpc scales (e.g., Chen et al. 2022; Cheng et al. 2023; Colina et al. 2023; Huang et al. 2023; Gillman et al. 2023, 2024; Boogaard et al. 2024; Crespo G'omez et al. 2024). Some of these early studies have reported evidence for undisturbed, disk-like morphologies that have been taken as evidence for the importance of secular growth for their stellar mass assembly (e.g., Cheng et al. 2023). Other studies have also reported evidence for disk-like morphologies but with a broad scatter in the non-parametric parameter space for morphological classification (Gillman et al. 2023), indicating a diverse population that nevertheless reflects the field population (Gillman et al. 2024). A number of studies have identified structures in the rest-frame near-infrared emission akin to bulges, spiral arms, and bars (e.g., Chen et al. 2022; Huang et al. 2023; Smail et al. 2023; Amvrosiadis et al. 2024), as well as examples of tidal features indicating dynamical interactions (e.g., Chen et al. 2022; Cathey et al. 2024). However, given the relatively narrow-field coverage of JWST in relation to the number density of SMGs, the submillimeter-selected galaxies studied so far in extragalactic deep fields have largely constituted submillimeter-fainter sources (i.e., ≲ few mJy at 850 µ m), which may undergo different triggering mechanisms than the submillimeter-brighter sources (e.g., Hayward et al. 2013b). More crucially, the current samples have lacked equally high-resolution ALMAdust continuum imaging, preventing the incorporation of the (energetically important) dust continuum emission into the physical interpretation.  \nHere, we present multi-band NIRCam imaging at 24.4 µ m (observed-frame) of 13 z ∼ 3 SMGs from the ALMA follow-up of the LABOCA Extended Chandra Deep Field South Survey (ALESS; Hodge et al. 2013; Karim et al. 2013). These submillimeter-bright (median S 870 ∼ 6 mJy) galaxies all have deep ALMA 870 µ m con-  \nimaging at 0.08 '' -0.16 '' resolution (Hodge et al. 2016, 2019), allowing us to compare their newly revealed rest-frame near-infrared emission to their dust continuum emission on the same (sub-)kpc scales. We begin in § 2 by discussing the sample selection, the JWST observations and data reduction, and the existing 870 µ m ALMA data. We present our analysis and results in § 3, including a comparison between the detected stellar and dust continuum emission, a visual classification, a curve-of-growth analysis of the effective radii, Galfit modeling, isophotal fitting to search for stellar bars, and finally a comparison of the NIRCam colors with ALMA 870 µ m. We follow this with a discussion in § 4. We end with our conclusions in § 5. Throughout this paper, we adopt a standard spatially flat ΛCDM cosmology from the Planck 2018 results with Hubble constant H 0 = 67.4 km s -1 Mpc -1 , dark energy density parameter Ω Λ = 0.685, and matter density parameter Ω M = 0.315 (Planck Collaboration et al. 2020). With these parameters, 1 '' corresponds to 7.9 kpc at z ∼ 3. We adopt a Chabrier (2003) initial mass function (IMF) and AB system magnitudes.", '2.1. Sample Selection': "The JWST observations targeted 13 SMGs 1 originally detected in the LABOCA 870 µ m survey of the ECDFS (Weiß et al. 2009) and then interferometrically identified as part of the ALESS survey (Hodge et al. 2013; Karim et al. 2013; Simpson et al. 2014; da Cunha et al. 2015; Danielson et al. 2017). These 13 galaxies were part of the sample subsequently followed up with highresolution, high-fidelity ALMA 870 µ m continuum mapping observations (Hodge et al. 2016, 2019, and see Section 2.3 for details). The sample is listed in Table 1. Eight of the targets have spectroscopic redshifts, also listed in Table 1; the remaining five have high-quality Magphys -based photometric redshifts which have been updated from those listed in da Cunha et al. (2015) to include the new JWST NIRCam+MIRI photometry (Li et al. in prep; Table 1).  \nDue to the original selection criteria for the highresolution ALMA 870 µ m follow-up, which was based on total 870 µ m flux density and the availability of (effectively randomly targeted) HST coverage (Chen et al. 2015), the sources targeted here include some of the  \nTable 1. Target Properties: Sample  \nf CO-based redshift from Westoby et al. (in prep.)  \nsubmillimeter-brightest sources from the ALESS parent sample: the median 870 µ m flux density of the sources is 6.4 mJy, with a 16th - 84th percentile range of 3.48.4 mJy. The galaxies nevertheless span a wide range of redshifts ( z ≃ 1 . 5 -4 . 5) and star formation rates ( ∼ 60-1000 M ⊙ yr -1 ). The median redshift of the targets is 3.4, with a 16th - 84th percentile range of 2.8 to 3.8. The median SFR is 350 M ⊙ yr -1 , with a 16th - 84th percentile range of 190-820 M ⊙ yr -1 . The sample is thus slightly higher-redshift and more highly starforming than the parent population of ALESS SMGs ( z phot = 2.7 ± 0.1 and SFR = 280 ± 70 M ⊙ yr -1 ; da Cunha et al. 2015). Finally, two of the sources (ALESS 17.1 and 45.1) are associated with known X-ray sources, of which one (ALESS 17.1) was identified as an X-ray AGN by Wang et al. (2013), while they found the X-  \nLESS 45.1 could be plausibly explained by star formation. The presence of a bright AGN could enhance the central brightness of the observed stellar distribution and lead to a higher S'ersic index; we note that this will not affect the majority of the sample. We also note that the X-ray source in ALESS 17.1 appears to be associated with a companion to the 870 µ m-bright SMG (see Fig. 5).", '2.2. JWST Observations & Data Reduction': "The JWST data presented here are part of a GO Cycle 1 Program (PID 2516, PIs: Hodge & da Cunha) to perform near- and mid-infrared imaging of submillimeter galaxies from the Hodge et al. (2016, 2019) samples, as described above. The 13 primary targets were selected from the Hodge et al. (2016, 2019) samples  \nsuch that they could be covered with four pointings with the NIRCam instrument. Each field is imaged in the F200W, F356W and F444W filters, using the 'INTRAMODULEBOX' 4-point dither pattern in order to cover a compact square region without gaps. Three NIRCam pointings used Module B only, the 4 th used both Modules A and B; with an on-source exposure time of 30 mins per field. The observations were executed between 2022 September and December. One observation was heavily impacted by slew artifacts on first execution, having been executed immediately following NIRCam observations of a very bright target; this was repeated in 2022 November.  \nThe data were reduced and calibrated using version 1.11.3 of the jwst calibration pipeline (Bushouse et al. 2023) and context 1119 of the JWST Calibration Reference Data System (CRDS). We used the standard JWST calibration pipeline, broadly following the recipe and modifications from the CEERS survey (Bagley et al. 2023). The first pipeline stage, performing detector calibrations, was run with snowball detections enabled and a jump step detection threshold of 6 σ . The 1 /f noise striping pattern was removed using the remstriping algorithm 2 , before proceeding with the standard settings for the calwebb image2 pipeline, which returns photometrically calibrated images for each exposure. The third pipeline stage ( calwebb image3 ) was run in stages with some modifications, following Bagley et al. (2023); this includes the background subtraction, which uses a tiered source mask approach to compute and subtract a pedestal value for the sky, and adjusts the variance arrays in the background-subtracted images accordingly. The resulting photometry for the targets is consistent with expectations based on the MUSYC catalog (Cardamone et al. 2010) and will be discussed further in Li et al. (in prep.). The PSF FWHM in the F200W, F356W, and F444W filters is 0.07 '' , 0.12 '' , and 0.14 '' , respectively. For the redshift range [median redshift] of our targets, the F200W, F356W, and F444W filters will trace restframe 0.35-0.8 µ m [0.5 µ m], 0.6-1.4 µ m [0.8 µ m], and 0.81.8 µ m [1 µ m] stellar continuum emission 3 , respectively, with the F444W resolution corresponding to physical scales of ∼ 0.9 -1.2 kpc for our range of target redshifts.  \nA common pixel sampling of 0.03 '' was used for all filters.  \nGiven the requirement for precise alignment to the ALMA data, the astrometric calibration is an important aspect of the data reduction, and this calibration step proved challenging due to the lack of point sources in our fields. The tweakreg step was used to perform both the relative and absolute alignments, the latter using the Gaia DR3 catalog (Gaia Collaboration et al. 2023). However, some modifications were used for the individual pointings to optimize the outcome of the adjustment. These include removing Gaia sources without proper motion information from the reference catalog 4 ; and removing poorly-centroided sources from the pipeline-created source catalogs and/or using PSF photometry to improve centroids, especially for partiallysaturated sources.  \nThe resulting relative and absolute astrometric accuracy was estimated from the final mosaics, by comparing measured point-source centroids between filters, and computing the average offset between the NIRCam and Gaia DR3 catalog coordinates, respectively. The resulting astrometric accuracies are shown in the Appendix. We also list the number of Gaia stars available for absolute alignment. We estimate a median absolute uncertainty in the astrometric solutions of 0.03 '' (i.e., one pixel).", '2.3. ALMA Data': "The ALMA continuum observations utilized for the primary targets in this study were taken as part of programs 2012.1.00307.S and 2016.1.00048.S and previously presented in Hodge et al. (2016, 2019), respectively. For both studies, the frequency setup was centered at 344 GHz (Band 7; 870 µ m observed frame) with 4 × 128 dual polarization channels covering the 8 GHz bandwidth in Time Division Mode (TDM). At the median redshift of our targets ( z ∼ 3), 870 µ m corresponds to 210 µ min the rest-frame. The naturally weighted continuum data were imaged with pixel scales of 0.02 '' and 0.01 '' and achieved angular resolutions of 0.17 '' × 0.15 '' and 0.10 '' × 0.07 '' for the analyses of Hodge et al. (2016) and Hodge et al. (2019), respectively (corresponding to physical scales of ∼ 1.3 × 1.2 kpc and ∼ 800 × 550 pc at a redshift of z ∼ 3). This range of spatial resolutions is similar to that achieved by NIRCam from 2.0-4.4 µ m (0.07 '' -0.14 '' ). The typical RMS noise achieved in the 870 µ m continuum images is 64 µ Jy beam -1 (for the 0.16 '' images) and 20 µ Jy beam -1 (for the 0.08 '' images).  \nFigure 1. JWST NIRCam RGB images of 8 '' × 8 '' regions ( ∼ 65 × 65 kpc at z ∼ 3) around FIR-luminous SMGs from the ALMA ALESS survey. North is up and East is to the left. The RGB images were made with the F444W (red), F356W (green), and F200W (blue) NIRCam filters. The positions of the 870 µ m sources are indicated with their ALESS ID numbers. All galaxies are clearly detected by NIRCam. ALESS 1.1, 1.2, and 1.3 are shown together in the same cutout. Note the long tidal tail revealed by the NIRCam imaging of ALESS 1.1/1.2, indicating an interaction. There is also an 870 µ m-detected companion (3.1-comp) to the South of ALESS 3.1 with the same spectroscopic redshift (see Section 3.1 for details).  \n<!-- image -->  \nGiven the angular resolution of the 870 µ m observations and the high SNR of the targets, the statistical astrometric accuracy is likely limited by the phase variations over the array to a few [10] mas for images with 0.08 '' [0.16 '' ] resolution 5 .", '3.1. Stellar emission compared to dust continuum': "Fig. 1 shows 8 '' × 8 '' ( ∼ 65 × 65 kpc at z ∼ 3) threecolor ('RGB') images of the 13 primary targets. The RGBimages were made with the NIRCam F444W (red), F356W (green), and F200W (blue) filters. The stellar counterparts to all of the SMGs are visible in the NIRCam RGB images. Interestingly, the NIRCam imaging also reveals an apparent tidal connection between ALESS 1.1 and 1.2 (which have spectroscopic redshifts of z = 4 . 674 and z = 4 . 669, respectively), providing the first direct evidence for an ongoing interaction. We note that while ALESS 1.3 also lies nearby, its photometric redshift ( z = 2.86 +0 . 46 -1 . 53 ) indicates it may be at lower redshift.  \nAlso apparent in Fig. 1 is a galaxy ∼ 5 '' south of ALESS 3.1 which was not detected at 870 µ m in the original ALESS survey (Hodge et al. 2013) but was strongly detected in the subsequent much deeper (and high-resolution) follow-up ALMA imaging utilized in this work, with an integrated flux density of S 870 =  \n1.07 ± 0.06 mJy. (This is fainter than the 3.5 σ limit of the original ALESS map which had σ 870 = 0.41 mJy beam -1 ; Hodge et al. 2013). ALESS 3.1-comp is only ∼ 40 km s -1 offset from ALESS 3.1 (based on a recent CO(5-4) detection; Westoby et al. in prep) and has a projected nuclear separation with ALESS 3.1 of ∼ 40 kpc. Given this, its strong submillimeter emission, and the tidal features apparent in the NIRCam images (Fig. 2), it is possible that this galaxy is undergoing an early-stage interaction/merger with ALESS 3.1. We therefore refer to it as '3.1-companion' (hereafter '3.1comp') and include it in the primary sample in Table 1 and for the analysis in the remainder of the paper.  \nTo investigate how the morphology of the sources changes with wavelength, Fig. 2 shows 3 '' × 3 '' cutouts ( ∼ 25 × 25 kpc at z ∼ 3) of the 13 sources in each of the NIRCam filters along with the corresponding RGB images. All of the sources are detected in at least the F356W and F444W NIRCam filters (i.e., longward of 3.5 µ mobserved-frame), including the red stellar counterpart to ALESS 17.1, which lies ∼ 0.8 '' from an optically bright spiral galaxy companion at the same redshift 6 . The majority of the sources show much weaker and/or centrally suppressed emission in NIR-  \nFig. 3 shows a comparison of the JWST NIRCam RGB imaging (zoomed to 2 '' × 2 '' ) and the high-resolution 870 µ m imaging for the 13 sources. The morphologies of the (typically) rest-frame near-infrared stellar and dust  \n<!-- image -->  \nFigure 2. JWST NIRCam cutouts (3 '' × 3 '' , or ∼ 25 × 25 kpc at z ∼ 3) of the 13 sources. For each galaxy, the columns show (from left): F200W, F356W, F444W, and the corresponding RGB image. Individual filters are scaled with a linear stretch between 0.25% and 99.5%, except for ALESS 17.1, where a maximum of 95.0% is used to highlight the faint red submillimeter galaxy to the West and the brighter companion to the East. Most of the targeted SMGs are heavily dust-attenuated even with NIRCam's 2 µ m filter but are significantly less obscured above 3.5 µ m ( ≃ 1 µ m rest-frame). Some examples of apparent tidal tails are highlighted with boxes in the filter where they are most visible, and the panels are labeled with redshift and visual classification (Disk, Indeterminate, or Merger; See Section 3.2.)  \n<!-- image -->  \nCam's F200W filter compared to the F356W/F444W filters, highlighting the red galaxy centers. Some apparent tidal features are also visible; we indicate a few prominent examples in Fig. 2 in the filter where they are most evident.  \ncontinuum images (traced by F444W and 870 µ m, respectively) show some striking similarities (as well as some notable differences), which we explore further below. Here, we note that the peak of the emission from the rest-frame near-infrared counterpart (measured from the peak pixel in the F444W filter) is consistent with the peak of the 870 µ m emission to within the F444W PSF FWHM of < 0.14 '' ( ∼ 1 kpc) for the majority (9/13) of the galaxies: ALESS 3.1, 3.1-comp, 9.1, 10.1, 15.1, 17.1,  \n<!-- image -->  \nFigure 3. JWST NIRCam versus ALMA 870 µ m imaging for the 13 SMGs. For each source, the columns show the NIRCam RGB image (constructed from F444W, F356W, and F200W; left), the ALMA 870 µ mimage (middle), and the ALMA contours in red overlaid on the RGB image (right). The ALMA images are scaled with a power-law stretch between -10 σ and 20 σ (exponent = 2.1), with ALMA contours starting at ± 2 σ and increasing in powers of √ 2. Panels are 2 '' × 2 '' , or ∼ 15 kpc at z ∼ 3. The NIRCam images reveal the rest-frame near-infrared morphologies of the galaxies on the same (sub-)kpc scales as the 870 µ m dust continuum, with the 870 µ m emission tracing the dust-reddened centers as well as some more extended features.  \n<!-- image -->  \nFigure 4. Effective radius versus wavelength for the sources (excluding 17.1 due to confusion with its optically bright companion). Individual sources are shown in different colors. The sample median is shown with the large black circles, where the error bars were calculated using bootstrapping. Including both the NIRCam and 870 µ m observations, we find that the median effective radius systematically decreases with increasing wavelength.  \n<!-- image -->  \n29.1, 45.1, and 112.1. For the remaining four sources (ALESS 1.1, 1.2, 1.3, and 76.1), the peak of the 870 µ m emission is ≥ 0.14 '' from the peak of the apparent counterpart in the F444W imaging, with the largest offset measuring 3 kpc (for ALESS 1.3, where the 870 µ memission may be more consistent with the fainter of what appear to be two marginally blended components in the NIRCam images). For the sources with significant offsets, this may indicate either that the ALMA 870 µ mand NIRCam imaging are tracing physically distinct components in these sources, or that the detected stellar emission traced by the NIRCam F444W filter is still experiencing significant dust attenuation. We explore these possibilities further below.", '3.2. Visual Classification': "We next visually classify the galaxies into three broad classes:  \nM): 'Merger': Sources with evidence for a merger or interaction;  \nI): 'Indeterminate';  \nD): 'Disk': Sources with no clear evidence for a merger or interaction, and which thus appear as undisturbed disks.  \nTo determine the visual classification for each galaxy, five coauthors independently classified the morphology of each of the 13 sources, taking into account the restframe near-infrared and 870 µ m morphologies (and further informed by the presence of a companion at the same redshift, if applicable). For each galaxy, we adopt the modal classification (preferred by at least 3/5); if no majority classification resulted from this method, then the galaxy was automatically assigned class I (indeterminate). We then calculated the fraction of sources in each of the three classes. Given the subjective nature of such classifications, we take the difference between the fraction in a given class from the above method and that from the most discrepant individual classification by any one classifier and adopt it as a rough estimate of the uncertainty.  \nWith the above method, we find the following results: M: 54 ± 23% (show evidence of mergers); I: 23 ± 31% (indeterminate); D: 23 ± 8% (appear as undisturbed disks). These results indicate that the majority of the galaxies in the current sample show signs of ongoing major/minor mergers or interactions, while only a minority show no clear evidence from the current datasets. Larger samples of galaxies spanning a wider range of 870 µ mflux density (e.g., Gillman et al. 2024) will be required to determine the importance of mergers/interactions for the SMG population as a whole.", '3.3. Flux density profiles': "To determine the effective radii of the galaxies in each filter, we create curves of growth. We first run SExtractor on the NIRCam images with a detection threshold of 3.0 times the local noise, a minimum number of pixels above the threshold of 12 at 0.03 '' sampling per pixel, and DEBLEND MINCONT set to 0.05. The latter parameter prevents deblending ALESS1.2 and ALESS10.1 into two sources each. We determine the aperture shape for each source using the results from the F444W filter (chosen to minimize extinction), with the centroid taken as the peak pixel in 870 µ m, which will be more robust against dust-obscuration effects than the NIRCam filters. We then use PhotUtils to measure the flux density in elliptical apertures as a function of distance along the major axis, with the total integrated flux density in each filter taken using the F444W maximum aperture (again to minimize extinction effects, and after ensuring the emission from all filters is covered). Finally, we use the same apertures to measure the flux density in the 870 µ m image for each source.  \nFig. 4 shows the effective radii (and median values) derived from this method versus observed wavelength, and including both the NIRCam and 870 µ m measurements. Here we exclude ALESS 17.1 due to confusion with its optically bright companion in the JWST images. The values for individual sources can be found in Table 2 and the flux density profiles themselves can be found in the Appendix.  \nWe find that the effective radius measured for the galaxies systematically decreases with increasing wavelength for the majority of the sources, with median values (and bootstrapped errors) of 0.61 '' ± 0.08 '' , 0.45 '' ± 0.06 '' , and 0.41 '' ± 0.05 '' for the F200W, F356W, and F444W filters, respectively. These values correspond to 4.5 ± 0.5 kpc, 3.4 ± 0.4 kpc, and 3.0 ± 0.3 kpc using the individual source redshifts, with 16th-84th percentile ranges of 3.0-5.3 kpc (F200W), 2.4-4.4 kpc (F356W), and 2.4-4.2 kpc (F444W). The effective radius then decreases even further for the majority of the sources at 870 µ m, with a median value of 0.23 '' ± 0.01 '' (1.8 ± 0.1 kpc, with a 16th-84th percentile range of 1.42.1 kpc). Overall, the median effective radius of the galaxies at F200W [F356W] is 48 ± 30% [10 ± 20%] larger than that measured at F444W, reflecting the relative brightness of the red galaxy centers at the reddest wavelengths. We explore whether this trend is due to dust obscuration in the following sections. Meanwhile, the median effective radius at F444W is 78 ± 21% larger than that measured with ALMA at 870 µ m. We note that if we had instead used the peak F444W pixel positions as the aperture centers, our conclusions would not change. We thus find much more compact dust continuum than rest-frame ∼ 1 µ m light, consistent with earlier indications from select sources (e.g., Chen et al. 2022; Gillman et al. 2024).", '3.4. Galfit modeling': "To compare the morphology of the NIRCam-detected stellar emission with that of the 870 µ m dust continuum in more detail, we fit both the F444W images (chosen to minimize dust extinction) and the 870 µ m images with single two-dimensional S'ersic profiles using Galfit (Peng et al. 2002, 2010). We consider only singlecomponent profiles for a more straightforward comparison between parameters. For the 870 µ m image fitting, we use the ∼ 0.16 '' 870 µ m ALMA maps available for all of our targeted galaxies and verify that we recover the results of Hodge et al. (2016), noting that Hodge et al. (2019) found S'ersic indices that agreed within the errors for sources also observed at the higher resolution (0.08 '' ) with ALMA.  \nFor the fitting of the F444W images, we obtain the most realistic PSFs by adopting the following strategy: We first use WebbPSF v1.2.1 (Perrin et al. 2014) to generate the F444W PSF oversampled by a factor of four. We then match the simulated PSF to our science data for each target using the webbpsf.setup sim to match data utility. Finally, to further improve these models, we follow Chen et al. (2022) and use Galfit to convolve the resulting PSFs with Gaussian profiles and fit them to nearby (unsaturated) stars. This results in reduced χ 2 values of ∼ 1, and we therefore adopt the best-fit (normalized) models as the final PSFs. We then run Galfit , masking the other sources in the field using the segmentation maps from SExtractor (and with the parameters presented in Section 3.3). For initial parameter guesses, we use the parameters derived from the 870 µ m image fitting, but we checked that our results were not sensitive to the exact input parameters. For ALESS 17.1, we simultaneously fit the sMG and its optically bright companion.  \nFig. 5 shows the results of running Galfit on the F444W images (where the values of the resultant S'ersic indices, axis ratios, and position angles can be found in Table 2). From the residual panels, it is evident that many sources contain structures beyond a single S'ersic profile. For example, roughly half of the sources (ALESS 3.1, 3.1-comp, 9.1, 15.1, 29.1, and 45.1) show evidence for excess emission beyond the models that is concentrated near the galaxy centers. A comparison of the Galfit residuals with the high-resolution 870 µ m imaging shows that this excess emission is aligned with the peak of the 870 µ memission in the majority of the cases. These central residuals are thus potentially evidence for bulge-like components and/or bright central AGN, as also reported in JWST imaging studies of submillimeterfainter sources (e.g., Chen et al. 2022).  \nBeyond these central stellar concentrations, there are additional residual structures evident in Fig. 5 that are not captured in the Galfit models. In particular, we see excess emission in the outskirts of several of the sources, further highlighting some of the features already apparent in Fig. 2 resembling tidal tails/plumes (e.g., see ALESS 3.1 and 3.1-comp). Interestingly, some of these potential tidal features are also partially detected in the 870 µ m continuum (see Fig. 5) for the sources with the highest-resolution 870 µ mimaging (e.g., ALESS 3.1, 9.1, and 112.1, though we note that the classification of these features can be somewhat subjective; Section 3.2).  \nFig. 6 shows a comparison of the morphological parameters derived for each galaxy from the F444W and  \nFig. 6 also shows a comparison of the best-fit effective radii ( R e ) measured in the 870 µ m and F444W images using the curve of growth analysis of Section 3.3. For the F444W images, the R e values derived from Galfit are consistent with the curve-of-growth values within the uncertainties. For the 870 µ m images, the two methods also produce consistent results for all sources except ALESS 1.1, where Galfit prefers (at > 3 σ ) a smaller radius (0.9 ± 0.1 kpc). Given that the (non-parametric) curve of growth analysis is less sensitive to morphological asymmetries, we consider these values the fiducial  \n<!-- image -->  \nFigure 5. Galfit modeling of the NIRCam F444W images for the 13 targeted SMGs. For each source, the columns show the NIRCam F444W image (1st column), the best-fit S'ersic profile model (2nd column), the residuals (3rd column), as well as a comparison of those residuals to the ALMA 870 µ m contours overlaid in red (4th column). The first and second column are scaled with a linear stretch between 0.25% and 99%, with contours showing 20%, 40%, 60%, and 80% of the peak intensity of the data (i.e., column 1). The third and fourth columns are shown with a linear stretch and the peak intensity scaled down by a factor of 5 to highlight the discrepancies. We note there is no qualitative difference to the residuals when the fits with freely varying n values are shown instead. 870 µ m contours start at ± 2 σ and increase in powers of √ 2, where σ =19 µ Jy. Panels are 2 '' × 2 '' , or ∼ 15 kpc at z ∼ 3. See Section 3.4 for details. The cyan cross overlaid on the third column for ALESS 17.1 shows the position of the known X-ray AGN (corrected for the median offset between the X-ray catalog and Gaia DR1; Luo et al. 2017), indicating it is associated with the optically bright companion rather than the SMG.  \n<!-- image -->  \n870 µ m images. (The optically bright companion to ALESS 17.1 is not included in this plot due to its lowS/N in the 870 µ m image; see Fig. 3). The S'ersic index ( n ), ellipticity ( b/a ), and position angle (PA) were taken from the best-fit Galfit models. As the uncertainties reported by Galfit are unrealistically small, we generate more realistic error bars by adding 10% fractional uncertainty in quadrature on all parameters. In addition, we fix the S'ersic indices n to 1.0 for both F444W and 870 µ m to ensure a more straightforward comparison when fitting the other relevant structural parameters ( b/a and PA). We take the difference between the parameters calculated with n held fixed and allowed to  \nvary and add this in quadrature as an additional systematic uncertainty for these parameters.  \nFigure 6. Comparison of the morphological parameters derived from the F444W images with those from the 870 µ m images. The panels showing the S'ersic index n ( upper left ), the axis ratio b/a ( bottom left ), and the major axis position angle (PA , bottom right ) were derived using Galfit with single S'ersic profile fits (Section 3.4). For the axis ratio and position angle, the S'ersic index is fixed to n = 1 to faciliate direct comparison. The effective radii R e ( upper right ) were derived non-parametrically using curves of growth (Section 3.3). All data points are color-coded by redshift. The literature SMGs in the top right panel show other sources with measurements at both wavelengths (compiled from Chen et al. 2022; Cheng et al. 2023; Kamieneski et al. 2023; Smail et al. 2023). Comparing the F444W- and 870 µ m-derived parameters, we find a strong positive correlation along the one-to-one line between the position angles, a tentative positive correlation between the effective radii, and no correlation between either the S'ersic indices or axis ratios.  \n<!-- image -->  \nvalues. However, we note that if we instead used the Galfit -derived results, it would only strengthen our finding that the dust continuum is more compact than the rest-frame ∼ 1 µ m (Section 3.3).  \nFrom Fig. 6, we see that some of the morphological parameters describing the dust continuum and detected stellar distributions are more tightly correlated than others. To quantify this effect, we calculate the Pearson correlation coefficient, ρ , between the F444W- and 870 µ m-derived parameters, which measures the strength of the linear monotonic correlation between two sets of data 7 . We find values of ρ (and corresponding proba-  \nities) of -0.14 (0.65), +0.68 (0.01), -0.02 (0.94), and +0.72 (0.005) for n , R e , b/a , and PA, respectively 8 .  \nThis exercise indicates that the position angle is the morphological parameter that shows the strongest correlation (of those tested) between the F444W and 870 µ m images. Moreover, as Fig. 6 shows, the position angle is not just strongly correlated, but it is correlated along the one-to-one line. This result indicates that the global morphologies of the rest-frame infrared and dust continuum emission are generally well-aligned. This close agreement can be attributed to the fact that the deep 870 µ m imaging - while more compact than the F444W emission for these galaxies - nevertheless still traces rest-frame near-infrared structures beyond the central  \nTable 2. Morphology with JWST and ALMA  \nbulges(/bars). In particular, the global PA agreement seen here appears to be driven in at least some cases by the (apparent) tidal features detected in the outskirts of some sources in both the 870 µ m emission as well as in the detected stellar emission (e.g., Figure 3). Meanwhile, the two clear outliers from the PA comparison in Fig. 6 are both understandable: ALESS 1.3 shows evidence for two marginally blended components in the NIRCam imaging (driving the PA), and the 870 µ m image of ALESS 3.1-comp is relatively shallow and compact compared to NIRCam.  \nIn addition to the strong correlation of the position angles, the effective radii show a tentative positive correlation. However, we also find - perhaps surprisingly that the axis ratios show no correlation. This may be due to a real physical difference in the gas and stellar distributions, dust obscuration, or the fact that the axis ratios are tracing structure at smaller radii in the 870 µ m emission.  \nWe also find no correlation between S'ersic indices in the F444W and 870 µ m emission. The Pearson's correlation coefficient suggests a mild anti-correlation, which may be caused by the influence of central dust obscuration on the F444W profiles, but it is not statistically significant. Finally, we note that the most extreme out-  \ner from this plot, ALESS 3.1-comp, shows a very high S'ersic index in the F444W image ( n = 5 . 4 ± 0 . 5), which may hint at the presence of a bright compact nucleus (either a star cluster or a type 1 AGN) in the center.", '3.5. Isophotal fitting: Searching for stellar bars': "Given the previous evidence for bar-like features in the dust continuum in SMG samples (e.g., Gullberg et al. 2019), and in particular, the elongated central features visible in the high-resolution images of some of the current targets (Hodge et al. 2019), one might expect to find evidence for stellar bars. We thus search for potential stellar bars in the sources by fitting elliptical isophotes to the F444W images using the astropy package PhotUtils . Using this method, the surface brightness, ellipticity ( e = 1 -b/a ), and position angle of the fitted ellipses are allowed to vary as a function of distance along the semi-major axis, which increases by a factor of 1.1 for each subsequent ellipse fit. The fitting continues until either the relative error in the local radial intensity gradient is > 0.5 for two consecutive ellipse fits or the outermost ellipse extends to the low signal-to-noise region (S/N ≤ 3). For comparison, we fit the highestresolution 870 µ m data available for each galaxy.  \nFollowing Guo et al. (2023), we carry out the ellipse fitting in a two-step process: (1) In the first iteration,  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 7. An example of the ellipse fitting in the F444W filter to ALESS 15.1. Panels (a) and (b) show the ellipticity and PA (respectively) in F444W of the best-fit elliptical isophotes as a function of semi-major axis. The vertical dashed lines show the scale of the respective PSFs. (Note that the errors on the fits for the second-to-last ring extend beyond the ranges plotted.) The F444W data is shown in panel (c) with log scaling, while a linear selection of best-fit isophotes are overplotted in panel (d). Panel (e) shows the best-fit elliptical isophotes for the 870 µ m data for comparison. ALESS 15.1 is the only (non-clumpy) galaxy which meets all of the criteria for hosting a rest-frame near-infrared bar in the F444W imaging (see Section 3.5), but we consider this a tentative classification due to the combination of the long bar length implied, the marginally significant rotation of the PA in the outer disk, lack of evidence in the Galfit residuals, and lack of PA correspondence between F444W and 870 µ m in the inner 'bar' region. We therefore find no evidence for strong rest-frame near-infrared bars in the targeted galaxies.  \n<!-- image -->  \nwe keep the center of each fitted ellipse as a free parameter. (2) In the second iteration, we fix the center of all ellipses to the mean center found from fitting the inner six ellipses (i.e., central 0.3 '' ) in step 1.  \nTo ensure that any bars identified are reliable, we exclude highly inclined sources (with i > 60 deg as inferred from the projected axis ratios of the outermost ellipse in F444W). This excludes only ALESS 1.1 and ALESS 1.2. We exclude also ALESS 17.1 due to confusion with its companion but include all ten remaining sources, including those that are clearly clumpy as a test of the robustness of the methodology.  \nWe then test whether the resulting radial profiles meet the following criteria, which are characteristic of stellar bars (e.g., Wozniak et al. 1995; Jogee et al. 2004; Marinova & Jogee 2007):  \n- 1. The ellipticity ( e ) rises smoothly to a global maximum e max > 0.4 (where the position of the maximum then defines the bar end) before decreasing again in the outer disk.\n- 2. The position angle remains fairly constant (∆PA < 20 deg) along the bar.\n- 3. In the transition region between the bar and outer disk, e drops by > 0.1.\n- 4. The position angle changes by > 10 deg beyond the bar end.  \nWe note that our choice of e max > 0.4 conservatively restricts the selection to strong bars, as most ( > 70%) bars in near-infrared images of local spiral galaxies have ellipticities > 0.4 (e.g., Marinova & Jogee 2007; Men'endezDelmestre et al. 2007) and Jogee et al. (2004) find that it is difficult to unambiguously identify weaker bars at  \nz > 1 (c.f., Sheth et al. 2008). See Fig. 1 in Jogee et al. (2004) for an example of radial profiles for a z ∼ 0 . 5 galaxy displaying these characteristic bar signatures.  \nWe find that no meaningful fit is possible for ALESS 10.1, which appears very clumpy/complex. The ellipse fitting also fails on the first or second iterations for ALESS 76.1 and 112.1 due to a strong dependence on the definition of the galaxy center. Given that these sources appear irregular (see the residuals from the GALFIT fitting in Fig. 5), they will not be discussed further.  \nOut of the remaining seven galaxies, only two sources meet all of the above criteria: ALESS 9.1 and ALESS 15.1. ALESS 9.1 is clearly clumpy, so we disregard it in the remaining discussion, except as a cautionary example that the criteria above are necessary, but not sufficient, to identify bars, particularly in sources with irregular morphologies. Indeed, we also find that the apparent tidal features in the outskirts of some sources (e.g., ALESS 3.1) can cause a second peak in the ellipticity at large radii, which can complicate bar identification using the above criteria, but we have verified that we do not find any additional bar candidates if we define the 'outer disk' using isophotes at smaller radii.  \nALESS 15.1 is thus the only remaining bar candidate. The ellipse fitting results for this source are shown in Fig. 7. If there is a bar in this source, the ellipse fitting suggests it has a maximum ellipticity of ∼ 0.5 and a semi-major axis length of ∼ 6.3 kpc. The latter value is within the range of bar lengths reported for z < 1 galaxies (e.g., Jogee et al. 2004) but on the large end of the range of values found for bars recently reported in z > 1 galaxies with JWST ( ∼ 2.0-7.5 kpc; Guo et al. 2023; Smail et al. 2023; Huang et al. 2023). Given 1) that this galaxy is at higher redshift ( z phot = 2 . 67) than all but one of those sources, where bar lengths should", 'Hodge et al.': "F444w data  \nbe shorter on average due to bar size evolution (e.g., Debattista & Sellwood 2000; Martinez-Valpuesta et al. 2006; Algorry et al. 2017); 2) the marginal change in position angle beyond the 'bar' end ( ∼ 10 deg; i.e., barely above the threshold for a bar classification); 3) the lack of evidence for a bar in the Galfit residuals; and 4) the lack of correspondence between the position angle of the F444W and 870 µ m isophotes in the 'bar' region (Fig. 7), where we would expect agreement between the two tracers in the case of a real bar; we consider this bar classification tentative at best. We therefore find no convincing evidence for strong stellar bars in the current rest-frame near-infrared data.", '3.6. NIRCam color-color diagrams versus ALMA 870 µ m': "Lastly, we quantitatively study the relationship between the 870 µ m dust continuum and detected stellar emission on ≲ 1 kpc scales to relate the influence of the former on the latter and the effect that dust obscuration has on our understanding of high-redshift star-forming galaxies. The methodology we utilize will be discussed further in Li et al. (in prep.). In short, we first resample the 870 µ m images to match the coarser NIRCam pixel scale. We use a 1000 × 1000 pixel grid with 30 mas pixels, fixing the coordinates of the target (Table 1) at the centre of the central pixel. We then PSF-match the pixel-aligned NIRCam F200W, F356W and ALMA 870 µ m images to the F444W filter images using PSF models from WebbPSF and 2D elliptical Gaussians for the NIRCam PSF and ALMA synthesized beam, respectively 9 . In particular, we generated a PSF-matching kernel for each image that requires convolution using the create matching kernel function in the PhotUtils package. We then convolve these images with the PSFmatching kernel using the convolve function in the astropy package after masking all other detected sources in the field of our targets.  \nTo create color-color diagrams, we adopt an aperture grid with apertures of 2 × 2 pixels (60 mas as these are Nyquist sampled). The fluxes within each aperture (i.e. resolution element) are summed and recorded in a table for all filters. The error on each measured flux density value was obtained by randomly placing 2 × 2 pixel apertures in the unmasked background region in the images. For apertures with ≥ 3 σ detections in all three NIRCam  \nfilters, we compute the F200W-F356W and F356WF444W colors and plot them in Fig. 8. The data points are color-coded by their corresponding 870 µ m flux density if the aperture has at least a 1 σ flux density in the ALMA 870 µ m map; otherwise the data points are colored black. In order to qualitatively assess the dust reddening properties in these sources, we also overplot vectors indicating the predicted impact on the colors of varying A V . We emphasize that these vectors are not fits to the data, but rather were obtained using Magphys spectral libraries (da Cunha et al. 2008, 2015) by finding the median color of stellar models within ± 0.1 dex of a given A V value, shown from A V = 0 to a maximum A V that depends on redshift.  \nFig. 8 shows that for the majority of the sources, we see an extended distribution of points in the NIRCam color-color space, where the direction of the extension roughly matches that expected from the A V vector. These distributions show that, as expected, the F200W-F356W color depends more strongly on A V in the models than the F356W-F444W color (i.e., the vectors are closer to vertical than horizontal); this is due to the fact that dust attenuation affects shorter wavelengths more strongly than longer wavelengths. By comparing the observed distributions to these vectors, we find that the maximum A V implied within each source on ∼ 1 kpc scales (for regions detected in F200W) ranges from A V < 4 in ALESS 3.1-comp (which is the only source below the detection threshold of the original ALESS catalog) to A V ∼ 8 (in ALESS 9.1 and 29.1). For comparison, the median A V of the ALESS sample from global spectral energy distribution (SED) fitting is only A V ∼ 2 mag (da Cunha et al. 2015).  \nFig. 8 also shows that for the majority of the sources, there is a correlation between redder NIRCam colors and 870 µ m surface brightness; i.e., the chance of a point being bright at 870 µ m is higher if it is red in both NIRCam colors. The clearest example of this correlation is seen in ALESS 15.1, but strong trends are also visible in ALESS 1.1, 1.2, 3.1, 9.1, 17.1, 29.1, 76.1, and 112.1 (and a weaker trend is visible in ALESS 1.3). This is a strong indicator that dust is the cause of the correlation.  \nFinally, we see that for some sources, there is a large variation in F200W-F356W color for a fixed F356WF444W color and 870 µ m surface brightness. ALESS 10.1 is the clearest example, but see also 1.2, 9.1, and 29.1. While ALESS 10.1 is visually classified as a merger/interaction, the appearance of this effect in the larger sample does not correlate with morphological classification (nor does it correlate with, e.g., F444W ellipticity). It is possible that this is due to a 'frosting' of lowA V stars in front of the bulk of the higherA V  \nFigure 8. NIRCam color-color diagrams showing F200W-F356W versus F356W-F444W colors for individual 60 × 60 mas regions within each galaxy with > 5 σ detections in all three NIRCam filters. The data points are color-coded by their corresponding 870 µ m flux surface density if the aperture has at least a 1 σ detection in the ALMA 870 µ m map; otherwise the data points are colored black. Vectors indicating the predicted impact on this color space of varying A V are shown for comparison and were obtained with Magphys (da Cunha et al. 2008, 2015, see text for details). The vectors are plotted with a shift 1 dex to the left for clarity. The majority of the sources show a correlation between redder NIRCam colors and 870 µ m surface brightness.  \n<!-- image -->  \nemission. The full implications of these diagrams for the distribution of stellar age, (dust-corrected) stellar mass, metallicity, and other derived parameters will be explored in Li et al. (in prep.).  \nTo further explore the general relation between dust in absorption and dust in emission, in Fig. 9 we plot A V versus dust mass surface density for our sources. Here the A V values were derived by taking each pixel in the color-color diagram of Fig. 8 with a > 1 σ 870 µ mflux density and projecting it onto the A V vector at the source redshift (shown in Fig. 8). To emphasize that these values are not fits to the data but purely empirical, we refer to this axis as A ' V in the figure. For the dust mass surface density, we use the calibration between 870 µ m flux density and dust mass from Dudzeviˇci¯ut˙e et al. (2020), where the dust masses were derived from Magphys . We also compare to a sample of nearby galaxies from Kreckel et al. (2013) with A V values derived from the reddening of the stellar continuum. Finally, we overplot three relations: the relation for a simple foreground dust screen between the emitter and absorber (based on the dust model from Draine & Li 2007 and using the dust-to-gas  \nratio from Draine et al. 2014); the relation for a mixed media model where stars and dust are uniformly mixed (from Calzetti et al. 2000 following Kreckel et al. 2013); and the empirical relation established in Kreckel et al. (2013) for their local galaxies (with the assumption of similar A V estimates for SMGs from stellar reddening and the Balmer decrement; Birkin 2022, Taylor et al. in prep.). The Kreckel et al. (2013) relation was determined on physical scales from 350 pc to 2 kpc, which is well-matched to the scales probed here ( ∼ 1 kpc), and it was derived by scaling the simple foreground dust screen model by a factor of ∼ 4 to account for geometric effects.  \nFig. 9 shows that the relation between dust in absorption and dust in emission for the galaxies studied in this work broadly agrees with the trend reported for local galaxies, while also extending the trend to higher dust mass surface densities and V-band extinction. The scatter observed for the current targets also appears broadly consistent with that seen in the local sample, both between and within individual galaxies; in particular, we note that NGC2146 (the only LIRG in the nearby sample) shows a similar spread to some of our sources. We  \nFigure 9. V-band extinction (inferred from the reddening of the stellar continuum) as a function of dust mass surface density for the 13 SMGs and a comparison sample of local galaxies from Kreckel et al. (2013). The black contours show the density of SMG data points, which are shown in different color circles for each galaxy. The local galaxies are shown as black squares, except for NGC2146 (the only LIRG in the sample), which is shown as the orange squares. Overplotted are three lines: a foreground dust screen model (solid red line); a mixed media model (dot-dashed black line); and the empirical relation found by Kreckel et al. (2013) for nearby galaxies (dashed red line; see their Equation 8) with the assumption of similar A V estimates for SMGs from stellar reddening and the Balmer decrement (Birkin 2022, Taylor et al. in prep.). The distribution and scatter seen for the SMGs are similar to that seen for local galaxies while extending the trend to higher dust mass surface densities. See Section 3.6 for further details.  \n<!-- image -->  \n/circledot  \ndiscuss the additional implications of these findings in the following section.", '4.1. The influence of dust on rest-frame optical/near-infrared morphologies': "The JWST NIRCam images presented here reveal the bulk of the underlying stellar emission in these highresolution ALMA-mapped submillimeter-selected galaxies. For this, observing at or beyond 3.5 µ m (i.e., rest wavelength ≃ 1 µ m at z ∼ 3) appears to be crucial. The morphologies can vary dramatically between F200W and F356W (i.e., between ∼ 0.5 and ∼ 1 µ m rest-frame), with the majority of the sources appearing either barely detected or with suppressed central emission at F200W compared to F356W (Fig. 2). However, at F356W, all sources are well detected, including several galaxies which were previously undetected in deep HST images at ≤ 1.6 µ m (ALESS 1.1, 17.1, and 76.1; Chen et al. 2015).  \nThe correspondence for the majority of sources between the morphology of the NIRCam F444W and the matched-resolution 870 µ m emission (particularly the peak positions and PAs) indicates that the spatial offsets previously reported between the 870 µ m continuum  \nand HST 1.6 µ m emission (e.g., Hodge et al. 2012; Chen et al. 2015; Calistro Rivera et al. 2018) were not due to true physical offsets, but rather that the morphologies inferred from the HST images were strongly affected by heavy structured dust obscuration, an effect also reproduced by radiative transfer modeling (Cochrane et al. 2019; Popping et al. 2022). For example, ALESS 29.1 was previously classified as a source with two separate stellar components, while four separate components were identified in ALESS 15.1 (Chen et al. 2015). By uncovering the obscured centers of the underlying disks, the NIRCam images now clearly reveal that these 'components' were simply due to differential dust obscuration in the disks. The higher A V values in the galaxy centers (Fig. 8) are thus causing the observed trend of larger effective radius at decreasing NIRCam wavelength (Fig.4). We note that other studies have also reported spatial offsets between the rest-frame UV/optical stellar and ISM emission in a variety of high-redshift populations (e.g., Cheng et al. 2020; Cochrane et al. 2021; Inami et al. 2022; Killi et al. 2024), with various theories posited beyond dust geometry, including strong stellar feedback or large-scale gas inflow/outflows. While there may certainly be other mechanisms operating in different galaxy  \npopulations, we caution that-given the strong impact of dust demonstrated here up to ≃ 1 µ m (rest-frame)structured dust obscuration should be ruled out before invoking more exotic explanations.  \nFor some of the highest-redshift sources, we see strong evidence that the galaxies are still affected by significant dust obscuration even in the longest-wavelength (F444W) NIRCam filter. In particular, examining the F444W versus 870 µ mcontinuum imaging for ALESS 1.1 ( z = 4 . 674; Fig. 3), as well as the comparison of the 870 µ m map with the F444W Galfit residuals (Fig. 5), the dust continuum appears to align with a suppression in the intensity of the detectable stellar distribution, suggesting significant dust obscuration. This suppression in the central region of the F444W image likely explains the very low F444W-derived S'ersic index ( n = 0 . 36 ± 0 . 04). We note that we find a similarly low S'ersic index at F444W for our other z > 4 . 5 source, ALESS 1.2. In these sources, F444W is probing shorter restwavelength emission (780 nm rest-frame) that is more prone to dust extinction. For these particularly highredshift/dusty sources, MIRI imaging will be important to more robustly characterize the stellar distributions.", '4.2. Implications for inclination estimates': 'While some of the most basic morphological parameters show agreement between the 870 µ m and F444W observations, others show no clear correlation (Fig. 6). The most notable parameter showing a lack of correlation is the axis ratio (i.e., ellipticity). This lack of correlation has significant implications, as the ellipticity is used to set the inclination in the most simple but commonly used dynamical mass estimates (often even without error margins), with the tacit assumption that these are perfectly flat, infinitely thin and perfectly circular disks. The lack of correlation between the ellipticites derived from different tracers shows that these inclinations are in fact very poorly determined, suggesting such dynamical mass estimates should be used with extreme caution.', '4.3. Rest-frame near-infrared vs. dust continuum sizes': 'The smaller extents measured for the ALMA 870 µ m emission compared to those inferred from the F444W emission (Fig. 4) suggests that the 870 µ m emission is more weighted to the dusty centers of the sources. This is consistent with recent results inferred for JWSTobserved SMGs compared with submillimeter continuum imaging of (largely) other samples (e.g., Chen et al. 2022), but this study extends these results to direct comparisons of these tracers within the same galaxies. For sources without evidence for significant obscuration  \nin the F444W filter, which could bias those size estimates to larger values, the smaller 870 µ m continuum sizes measured here would suggest that these galaxies are still undergoing rapid morphological transformation, with the central starbursts likely building bulges that will eventually come to dominate the stellar mass distributions (for which the F444W filter serves only as a weak proxy). This result would strengthen the challenge these sources pose to recent hydrodynamical simulations, a number of which predict the opposite trend for galaxies in this stellar mass range (i.e., submillimeter extents that are larger than the stellar mass extents; e.g., Cochrane et al. 2019; Popping et al. 2022). Given the significant obscuration still present in even the longestwavelength NIRCam images in at least some sources, this result needs to be tested via the creation of (dustcorrected) stellar mass maps (see Smail et al. 2023 and Li et al., in prep).', '4.4. Evidence for mergers and interactions': "While the NIRCam observations have now detected the dust-reddened centers of the galaxies-suggesting less irregular morphologies than implied from the heavily obscured 1.6 µ m imaging (Chen et al. 2015)-this does not imply that the galaxies all appear as undisturbed disks in F444W. On the contrary, the NIRCam images reveal evidence for mergers/interactions in the majority of the sources.  \nIn particular, three of the SMGs have direct evidence for companions: ALESS 1.1 (with 1.2), ALESS 3.1, and ALESS 17.1. The most dramatic feature is the long tidal tail/stellar bridge that is now clearly visible between ALESS 1.1 and 1.2 at z = 4 . 67 (Fig. 1). With a projected separation of ∼ 20 kpc and a velocity separation of ∼ 320 km s -1 , this suggests that ALESS 1.1 and 1.2 are undergoing an early-stage (violent) interaction. Interestingly, both sources are detected in CO(54) emission (Birkin et al. 2021), where ALESS 1.1 has a double-peaked profile with a very wide FWHM velocity width of 1300 km s -1 . This could indicate that ALESS 1.1 is itself a late-stage merger or that the gas disk is dynamically disturbed. With ALESS 1.3 lying in the same field-and showing evidence for two very close components, where the southern component is redder (Fig. 2) and may be the primary source of its submillimeter emission (Fig. 3)-this means that the original submillimeter source 'LESS 1' (i.e., the brightest source in the original LABOCA 'LESS' catalog, with S 850 = 14.5 ± 1.2 mJy; Weiß et al. 2009), was actually made up of at least two or more different and likely unrelated major merger events within its 19 '' beam (FWHM). Meanwhile, while ALESS 3.1 and 3.1-comp are also likely undergoing an interac-  \note that ALESS 3.1 also has a broad (870 km s -1 ) and asymmetric CO line profile (Birkin et al. 2021), potentially indicating that ALESS 3.1 itself is undergoing a late-stage interaction/merger. Thus the original submillimeter source 'LESS 3' (i.e., the third brightest source in the LESS catalog) also appears to consist of at least one or more mergers/interactions, implying that the brightest LESS (i.e., single-dish LABOCA) sources are mergers/interactions 10 (and see Mitsuhashi et al. 2021 for a discussion of an overdensity associated with the brightest SCUBA-2 source in COSMOS). This is a direct confirmation of trends reported previously based on the multiplicity of sources comprising bright single dish-selected sources (e.g., Stach et al. 2018; Simpson et al. 2020).  \nThe NIRCam images (and Galfit residuals) also reveal evidence for interactions via features resembling stellar tidal tails/plumes and clumps in ALESS 10.1 and 76.1, as well as potential evidence in ALESS 9.1 and 112.1. Some of the apparent tidal features were previously detected in the dust continuum in the sources with the highest resolution 870 µ m emission maps (e.g., ALESS 3.1, 112.1; Hodge et al. 2019) and interpreted as either tidally induced spiral arms or the star-forming knots in an interaction/merger. ALESS 3.1 was already discussed above and indeed shows evidence of a merger/interaction also in the NIRCam images. With the new NIRCam imaging showing possible evidence of a warped disk (Fig. 5), ALESS 112.1 may also be a latestage merger, though this case is less conclusive and is thus marked as Indeterminate. ALESS 9.1, which also has a potential tidal feature detected in both NIRCam and 870 µ m continuum (to its South), appears very clumpy in its rest-frame near-infrared emission, possibly due to an advanced interaction/merger (though again marked 'Indeterminate' due to the lack of clear evidence). In fact, of the 13 sources examined, only ALESS 15.1, 29.1 and 45.1 show no clear evidence in the current data for mergers/interactions, though we caution that these sources may still be late-stage mergers or interactions where the merger/interaction is less obvious due to projection effects or low-mass companions. Conversely, it should be noted that the departures from asymmetry used here to visually classify the galaxies as mergers/interactions in cases without a clear companion may also result from accretion from a cooling halo in gas-rich systems, which can be episodic and sporadic (e.g., Bland-Hawthorn et al. 2024).  \n10 We note that LESS 2 was not included in the high-resolution ALMA/JWST follow-up due to the sample selection.  \nKeeping these caveats in mind, the current interpretations result in a merger fraction of 54 ± 23% and a fraction of undisturbed disks of only 23 ± 8%consistent (within the errors) with the larger but typically submillimeter-fainter sample of Gillman et al. (2024), who report 20 ± 5% candidate late-stage mergers plus 40 ± 10% potential minor mergers. Confirming these interpretations will ultimately require highresolution kinematic data. However, regardless of their nature, the fact that a number of these (apparent tidal) features are detected also in the submillimeter continuum indicates that-despite the centrally concentrated 870 µ memitting regions-the high dust column densities and active star formation in these galaxies currently extends well beyond the existing central bulge, as we discuss further in Section 4.6.", '4.5. Stellar colors & the structure of the ISM': "One of the most notorious degeneracies in SED fitting is that between dust attenuation and stellar age, and this has long been a concern for SMG samples in particular (e.g., Hainline et al. 2011; Michaglyph[suppress]lowski et al. 2014; Dudzeviˇci¯ut˙e et al. 2020). Although SED fitting that requires energy balance with the FIR emission can help by providing additional constraints on the dust attenuation (e.g., da Cunha et al. 2015), the assumptions behind this energy balance have not been thoroughly tested in the high-redshift universe. In particular, the dust model used in such SED fitting codes is calibrated based on the Milky Way (e.g., Draine & Li 2007) and may not be applicable beyond the local universe and/or in extreme starbursts.  \nWith this in mind, the relations seen in Figs. 8 and 9 are notable. In particular, the correlation between NIRCam colors and 870 µ m surface brightness seen for the majority of the sources (Fig. 8) implies that the primary driver behind the red stellar colors in these SMGs is dust rather than stellar age. Fig. 9 further implies that the dust-to-stellar distributions in these SMGs are similar (on the same ∼ kpc-scales) to those in nearby star-forming galaxies, particularly the LIRG NGC2146, which exhibits similar A V for gas and stars (as also seen in high-redshift SMGs; Birkin 2022, Taylor et al. in prep.). The bulk of the SMG data points fall in a similar regime between the foreground dust screen and mixed media models as the Kreckel et al. (2013) relation, indicating a similar correction factor to the foregroundscreen model is applicable for these z ∼ 3 sources as found at z ∼ 0. This suggests that the ISM in these galaxies consists of a similar combination of mixed dust and emitting material and foreground screen components (c.f. Tomiˇci'c et al. 2017 for a higher-resolution  \nstudy in M31), with the remaining scatter likely due to geometric effects and non-uniformity in the dust distribution on scales beneath our resolution.  \nAt low dust mass surface densities ( < 10 6 M ⊙ kpc -2 ), it is noteworthy that the SMG data points remain below the limit set by the dust screen model despite the NIRCam data's sensitivity to higher A V values. This suggests that this simple model-based on Milky Waylike dust-remains valid for these z ∼ 3 sources. Meanwhile, at the highest dust mass surface densities probed by these data ( ≳ 10 7 -10 8 M ⊙ kpc -2 ), the Kreckel et al. (2013) relation would predict A V values up to ∼ 100 mag, far in excess of what we measure. However, for these particularly high-density regions, our measured A V values likely correspond to lower limits due to the maximum optical depth to the stellar continuum emission that we can detect with NIRCam. Given that the nearby galaxies used to calibrate the Kreckel et al. (2013) relation did not sample such high dust mass surface densities, this does not affect our conclusions.  \nIn summary, the ISM structure in these SMGs appears largely similar to local star-forming galaxies on ∼ kpc scales. This similarity implies that some of the most basic assumptions that go into SED-fitting codes such as, e.g., Magphys (i.e., the dust model) are still valid for these z ∼ 3 galaxies.", '4.6. Rest-frame near-infrared and gas/dust bars': "Lastly, we discuss the presence of bars in the sources. As discussed in detail in Hodge et al. (2019), there are clearly visible bar-like sub-structures in the 870 µ m images of some of the targets. While bars are often associated with secular evolution due to their tendency to arise naturally in dynamically cool disk galaxies (e.g., Athanassoula & Sellwood 1986) they can also be triggered by interactions (Barnes & Hernquist 1991; Cox et al. 2008; Hopkins et al. 2009), and thus are not mutually exclusive with the tidal signatures observed here in a large fraction of the sample. These bars are also consistent with the most recent theoretical expectations for gas-rich galaxies at high-redshift: Though previous work on bar formation found that the presence of (inert) gas reduces a stellar bar's lifetime (Villa-Vargas et al. 2010) or weakens the bars that form (Athanassoula et al. 2013), more recent work using hydrodynamic N-body simulations to predict bar formation in gas-rich disks at high redshift finds that turbulent gas accelerates bar formation, predicting bar-like phenomena even in fully gas-dominated turbulent disks (BlandHawthorn et al. 2024). Finally, these bars can also be an important mechanism to systematically drive gas inwards in sources with no current evidence for merg-  \nrs/interaction, helping to explain the intense dusty star formation (and possibly supermassive black hole growth).  \nGiven the high A V values implied in the center of these sources, it is perhaps not surprising that we see no evidence for strong bars in the rest-frame near-infared images presented here, as the stellar components could easily be hidden behind the high implied dust columns. Another possible (but more speculative) explanation comes from the work of Bland-Hawthorn et al. (2024), who predict that as the gas fraction increases, the role of any stellar bar becomes less important, with mostly gas bars emerging via radial shear flows in galaxies with the highest gas fractions. Considering the high gas fractions inferred for these sources (e.g., Swinbank et al. 2014), it is thus possible that their stellar bars are weaker than those in less gas-dominated sources. In either scenario, these dust-rich structures would necessarily be young. Notably, Bland-Hawthorn et al. (2024) find that the bars formed in the most gas-rich disks collapse to form bulges after ∼ 1.5 Gyr, which would then imply that the bar-like features observed in these sources are an immediate precursor to further significant bulge growth. Ultimately, larger samples of sources covering a range of submillimeter flux densities will be necessary to determine the importance of gas/dust and/or stellar bars in the mass assembly of the general SMG population.", '5. CONCLUSIONS': "We have presented JWST NIRCam imaging of 13 z ∼ 3 submillimeter-selected galaxies that have uniquely deep, high-resolution (0.08 '' -0.16 '' ) ALMA 870 µ mimaging which previously mapped their dust disks on ∼ 0.51 kpc scales. Our main findings are:  \n- · All of the sources are securely detected by the JWST NIRcam imaging, including several galaxies which were previously undetected in deep HST 1.6 µ m images. For this, observing at or beyond ∼ 3.5 µ m ( ≃ 1 µ m rest-frame at z ∼ 3) appears to be crucial, as the morphologies of the sources vary dramatically between F200W and F356W (i.e., between ∼ 0.5 µ m and ∼ 1 µ m rest-frame).\n- · With the dust-reddened galaxy centers now more visible, the newly revealed rest-frame near-infrared morphologies show some clear similarities to the 870 µ m dust continuum images (specifically the position angles and peaks), demonstrating that the spatial offsets previously reported between the 870 µ m and HST morphologies were due to strong differential dust obscuration.  \n- · However, we find no correlation between the axis ratio (i.e., ellipticity) derived from the rest-frame near-infrared and 870 µ m dust continuum. As the ellipticity is commonly used to set the inclination in dynamical mass estimates, this lack of correlation shows that these inclinations are in fact very poorly determined and that resulting dynamical mass estimates should be used with extreme caution.\n- · In some sources (particularly those at the highest redshifts), we see convincing evidence that the galaxies are still affected by significant dust obscuration even in the longest-wavelength (F444W) NIRCam filter. For these particularly dusty/highredshift sources, MIRI imaging will be crucial to more robustly characterize the stellar distributions.\n- · Due to the high level of central obscuration in the galaxies, we find that the median effective radius of the galaxies systematically decreases with increasing NIRCam wavelength. We nevertheless find that the F444W sizes are still 78 ± 21% larger than those measured at 870 µ m. Thanks to the unique depth and resolution of the 870 µ m images available for these targets, this strengthens the challenge posed to recent hydrodynamical simulations claiming the opposite trend (Cochrane et al. 2019; Popping et al. 2022). However, given the significant obscuration still present in even the longestwavelength NIRCam filter in at least some of the sources, this result needs to be tested via the creation of stellar mass maps (Li et al., in prep.).\n- · The NIRCam imaging reveals clear evidence for mergers/interactions (e.g., tidal tails/plumes) in the majority of the sources, with 54 ± 23% visually classified as mergers/interactions, 23 ± 31% classified as indeterminate, and only 23 ± 8% appearing as undisturbed disks (though we caution that these interpretations may miss late-stage mergers or be affected by episodic gas accretion from cooling halos). Several of the (apparently tidal) features are also detected in 870 µ memission, indicating that the active star formation in these galaxies currently extends well beyond the existing central bulges.\n- · We find a clear correlation between redder NIRCam colors and 870 µ m surface brightness on ∼ 1 kpc scales. This is a strong indicator that the SMGs have red NIRCam colors due to dust rather than stellar age. We further show that the relation  \nbetween dust in absorption and dust in emission for the galaxies studied in this work is broadly similar to the trend seen in nearby star-forming galaxies. This suggests that the ISM structure on ∼ kpc-scales in these z ∼ 3 galaxies is similar to that in z ∼ 0 sources.  \n- · We find no evidence for strong bars in our targets in the rest-frame near-infrared. This suggests that the elongated bar-like structures seen in the highresolution 870 µ m images are highly dust-obscured and/or gas-rich and thus young, implying they are the immediate precursors to further significant bulge growth.  \nTaken together, these findings suggest we are witnessing heavily obscured and largely interaction-induced bulge formation events at the centers of these massive star-forming galaxies. The present study thus demonstrates the joint power of JWST and ALMA for uncovering the morphologies and likely formation histories of high-redshift dusty galaxies. Future work will incorporate the multi-wavelength data into resolved SED fitting to explore the resolved physical properties, stellar masses, and resolved kinematics (Li et al. in prep; Westoby et al. in prep). Combined with future large samples of sources covering a range of submillimeter flux densities, such studies will be crucial to shed further light on the dominant mechanisms governing the mass assembly of massive, dusty galaxies at high-redshift and the relation of these sources to the larger high-redshift galaxy population.  \nWe thank Joss Bland-Hawthorn for his feedback on the manuscript. JH and BW acknowledge support from the ERC Consolidator Grant 101088676 ('VOYAJ'). EdC and JL acknowledge support from the Australian Research Council (projects DP240100589 and CE170100013). AMS and IRS acknowledge STFC (ST/X001075/1). CCC acknowledges support from the National Science and Technology Council of Taiwan (111-2112M-001-045-MY3), as well as Academia Sinica through the Career Development Award (AS-CDA112-M02). RD acknowledges support from the INAF GO 2022 grant 'The birth of the giants: JWST sheds light on the build-up of quasars at cosmic dawn' and by the PRIN MUR '2022935STW', RFF M4.C2.1.1, CUP J53D23001570006 and C53D23000950006. TRG acknowledge funding from the Cosmic Dawn Center (DAWN), funded by the Danish National Research Foundation (DNRF) under grant DNRF140. KK acknowledges support from the Knut and Alice Wallenberg Foundation. MR is supported by the NWO Veni project 'Under the lens' (VI.Veni.202.225). This work is based in part on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program # 2516. Support for program #2516 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.00048.S, ADS/JAO.ALMA#2012.1.00307.S, and ADS/ JAO.ALMA#2011.1.00294.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The Hubble Source Catalog is based on data from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESAC/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). SK acknowledges advice and support from Anton Koekemoer, Armin Rest, and Pablo Perez Gonzalez for the reduction of NIRCam data. This research made use of astropy , a community developed core Python package for astronomy (Astropy Collaboration et al. 2013, 2018) hosted at http://www.astropy.org/, matplotlib (Hunter 2007), numpy (van der Walt et al. 2011), scipy (Jones et al. 2001), and of topcat (Taylor 2005).  \nTable 3. NIRCam final astrometric alignment accuracy  \n- a relative accuracy is measured between F200W/F356W/F444W filters.", 'A. ADDITIONAL TABLES & FIGURES': 'Table 3 shows the astrometric accuracy we estimate for each of the NIRCam pointings. Both relative and absolute accuracy are indicated, as well as the number of Gaia stars (with catalog proper motions) available for the absolute astrometric alignment.  \nFig. 10 shows the cumulative fraction of integrated flux density versus radius in each of the three NIRCam  \nfilters for the targeted SMGs (including 3.1-comp, and excluding 17.1 due to confusion with its optically bright companion). The corresponding curve for the highestresolution 870 µ m image available for each source is also shown. These curves are not deconvolved from their (respective) PSFs, but this has a minor effect given that the sources are well-resolved. A curve corresponding to the F444W PSF, which has the lowest resolution of the NIRCam filters, is shown as an example.', 'REFERENCES': 'Algera, H. S. B., Inami, H., Oesch, P. A., et al. 2023, MNRAS, 518, 6142, doi: 10.1093/mnras/stac3195 Algorry, D. G., Navarro, J. F., Abadi, M. G., et al. 2017, MNRAS, 469, 1054, doi: 10.1093/mnras/stx1008 Amvrosiadis, A., Lange, S., Nightingale, J., et al. 2024, arXiv e-prints, arXiv:2404.01918, doi: 10.48550/arXiv.2404.01918 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, aap, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip"ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Athanassoula, E., Machado, R. E. G., & Rodionov, S. A. 2013, MNRAS, 429, 1949, doi: 10.1093/mnras/sts452 Athanassoula, E., & Sellwood, J. A. 1986, MNRAS, 221, 213, doi: 10.1093/mnras/221.2.213 Bagley, M. B., Finkelstein, S. L., Koekemoer, A. M., et al. 2023, ApJL, 946, L12, doi: 10.3847/2041-8213/acbb08 Barger, A. J., Cowie, L. L., Sanders, D. B., et al. 1998, Nature, 394, 248, doi: 10.1038/28338  \nBarger, A. J., Wang, W.-H., Cowie, L. L., et al. 2012, ApJ, 761, 89, doi: 10.1088/0004-637X/761/2/89  \nBarnes, J. E., & Hernquist, L. E. 1991, ApJL, 370, L65, doi: 10.1086/185978 Baugh, C. M., Lacey, C. G., Frenk, C. S., et al. 2005, MNRAS, 356, 1191, doi: 10.1111/j.1365-2966.2004.08553.x Birkin, J. E. 2022, PhD thesis, Durham University, UK Birkin, J. E., Weiss, A., Wardlow, J. L., et al. 2021, MNRAS, 501, 3926, doi: 10.1093/mnras/staa3862 Bland-Hawthorn, J., Tepper-Garcia, T., Agertz, O., & Federrath, C. 2024, arXiv e-prints, arXiv:2402.06060, doi: 10.48550/arXiv.2402.06060 Boogaard, L. A., Gillman, S., Melinder, J., et al. 2024, ApJ, 969, 27, doi: 10.3847/1538-4357/ad43e5 Bouwens, R., Gonz\'alez-L\'opez, J., Aravena, M., et al. 2020, ApJ, 902, 112, doi: 10.3847/1538-4357/abb830 Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2023, JWST Calibration Pipeline, 1.10.2, Zenodo, doi: 10.5281/zenodo.7829329 Calistro Rivera, G., Hodge, J. A., Smail, I., et al. 2018, ApJ, 863, 56, doi: 10.3847/1538-4357/aacffa Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682, doi: 10.1086/308692  \nFigure 10. Cumulative fraction of integrated flux density in JWST NIRCam and ALMA 870 µ m versus radius along the major axis for the targeted SMGs (including 3.1-comp, and excluding 17.1 due to confusion with its optically bright companion). The dotted line indicates the PSF in the F444W filter. The dashed lines indicate the effective radius in each filter, which we find to generally decrease with increasing NIRCam wavelength in these galaxies. Note that the curve corresponding to 870 µ m may stray above 1.0 or even decrease with radius due to the presence of noise in the ALMA imaging, indicated by the gray-toned band (representing ± 1 σ ). In the majority of the sources, the effective radius decreases with increasing NIRCam wavelength, and the effective radius of the ALMA 870 µ m emission is more compact than that from even the reddest JWST NIRCam filter.  \n<!-- image -->  \nCardamone, C. N., van Dokkum, P. G., Urry, C. M., et al. 2010, ApJS, 189, 270, doi: 10.1088/0067-0049/189/2/270  \nCasey, C. M., Narayanan, D., & Cooray, A. 2014, PhR, 541, 45, doi: 10.1016/j.physrep.2014.02.009 Cathey, J., Gonzalez, A. H., Lower, S., et al. 2024, ApJ, 967, 11, doi: 10.3847/1538-4357/ad33c9 Chabrier, G. 2003, PASP, 115, 763, doi: 10.1086/376392 Chen, C.-C., Smail, I., Swinbank, A. M., et al. 2015, ApJ, 799, 194, doi: 10.1088/0004-637X/799/2/194  \nChen, C.-C., Harrison, C. M., Smail, I., et al. 2020, A&A, 635, A119, doi: 10.1051/0004-6361/201936286  \nChen, C.-C., Gao, Z.-K., Hsu, Q.-N., et al. 2022, ApJL, 939, L7, doi: 10.3847/2041-8213/ac98c6  \nCheng, C., Ibar, E., Smail, I., et al. 2020, MNRAS, 499, 5241, doi: 10.1093/mnras/staa3036  \nCheng, C., Huang, J.-S., Smail, I., et al. 2023, ApJL, 942, L19, doi: 10.3847/2041-8213/aca9d0  \nCimatti, A., Cassata, P., Pozzetti, L., et al. 2008, A&A, 482, 21, doi: 10.1051/0004-6361:20078739  \nCochrane, R. K., Hayward, C. C., Angl´es-Alc´azar, D., et al. 2019, MNRAS, 488, 1779, doi: 10.1093/mnras/stz1736 Cochrane, R. K., Best, P. N., Smail, I., et al. 2021, MNRAS, 503, 2622, doi: 10.1093/mnras/stab467'}