Massive star-forming galaxies at high redshift require a supply of molecular gas from their gas reservoirs that is replenished by infall from the surrounding circumgalactic medium to sustain their immense star formation rates. Our knowledge of the extent and morphology of cold-gas reservoirs of early galaxies is still in its infancy, however. We present the results of stacking more than 80 hours of JVLA observations of CO(1–0) emission, which traces the cold molecular gas, in 19 z = 2.0−4.5 dusty star-forming galaxies from the AS2VLA survey. The visibility-plane stack reveals extended emission with a half-light radius of 3.8 ± 0.5 kpc, which is a factor of 2–3 more extended than the dust-obscured star formation and 1.4 ± 0.2× more extended than the stellar emission revealed by the JWST. Stacking the [C i](1–0) observations for 10 galaxies from our parent sample yielded a half-light radius ≤2.6 kpc, which is marginally smaller than CO(1–0). The CO(1–0) size is also comparable to that of the [C ii] haloes detected around high-redshift star-forming galaxies. This suggests that these arise from molecular gas. Photo-dissociation region modelling indicates that the extended CO(1–0) emission arises from clumpy dense clouds and not from smooth diffuse gas. Our results show that the bulk (up to 80%) of the molecular gas in these galaxies resides outside the star-forming region with only a small part directly contributing to the star formation.

We present JWST NIRCam imaging targeting 13 z ~ 3 infrared-luminous (L_IR ∼ 5 × 10¹²L_⊙) galaxies from the ALESS survey with uniquely deep, high-resolution (0 . ″ 08-0 . ″ 16) Atacama Large Millimeter/submillimeter Array 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 at the same (sub)kiloparsec 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 μm and Hubble Space Telescope 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, emphasizing 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 bar-like features seen at 870 μm are highly obscured and/or gas-dominated structures that are likely early precursors to significant bulge growth.

We present the initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO(J = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies (SFGs) at z = 2-5 with existing mid-J CO detections from the Atacama Large Millimeter/submillimeter Array and NOrthern Extended Millimeter Array, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative (∼1.5σ-2σ) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23 × 10¹⁰ (α _CO/1) M _⊙, with gas mass fractions, f _gas = M _mol/(M _*+M _mol), of 0.1-0.8 and a median depletion time of (140 ± 70) Myr. We find median CO excitation ratios of r ₃₁ = 0.75 ± 0.39 and r ₄₁ = 0.63 ± 0.44, with significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width, or Σ_SFR. We only find a tentative positive correlation between r ₄₁ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times, and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive SFGs at high redshift, and the importance of CO(1-0) observations to robustly constrain their total molecular gas content and interstellar medium properties.