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.

The CO(1-0) and [C i](1-0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium is likely to be denser, there have been limited direct comparisons of both ground-state transitions. We present a comparison of [C i](1-0) and CO(1-0) emission in 20 unlensed dusty, star-forming galaxies at z ≥ 2-5. The CO(1-0)/[C i](1-0) ratio remains constant up to z = 5, supporting the reliability of [C i](1-0) as a gas-mass tracer. We use the CO(1-0), [C i](1-0), and 3 mm dust continuum measurements to cross-calibrate their respective gas mass conversion factors, finding no dependence of these factors on either redshift or infrared luminosity. Radiative transfer modeling shows that the warmer cosmic microwave background (CMB) at high redshift can significantly affect the [C i] as well as CO emission, which can change the derived molecular gas masses by up to 70% for the coldest kinetic gas temperatures expected. Nevertheless, the magnitude of the CMB effect on the CO/[C i] ratio is within the known scatter of the L CO ′ − L [ CI ] ′ relation. Precisely determining the CMB effect on individual line intensities would require well-sampled spectral line energy distributions to robustly model the gas excitation conditions. Finally, we note that adopting a variable CO gas-mass conversion factor for different galaxy populations implies [C i](1-0) and dust conversion factors that differ from canonically assumed values. However, the revised conversion factors are consistent with expectations for (super)solar metallicities likely to be found in high-redshift dusty galaxies.

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.

The molecular gas in the interstellar medium (ISM) of star-forming galaxy populations exhibits diverse physical properties. We investigate the CO excitation of 12 dusty luminous star-forming galaxies at 2-4 by combining observations of the CO from to. The spectral line energy distribution (SLED) has a similar shape to NGC 253, M82, and local ultra-luminous infrared galaxies, with much stronger excitation than the Milky Way inner disc. By combining with resolved dust continuum sizes from high-resolution 870 m ALMA observations and dust mass measurements determined from multiwavelength spectral energy distribution fitting, we measure the relationship between the CO SLED and probable physical drivers of excitation: star-formation efficiency, the average intensity of the radiation field, and the star-formation rate surface density. The primary driver of high-CO excitation in star-forming galaxies is star-formation rate surface density. We use the ratio of the CO(3-2) and CO(6-5) line fluxes to infer the CO excitation in each source and find that the average ratios for our sample are elevated compared to observations of low-redshift, less actively star-forming galaxies and agree well with predictions from numerical models that relate the ISM excitation to the star-formation rate surface density. The significant scatter in the line ratios of a factor within our sample likely reflects intrinsic variations in the ISM properties that may be caused by other effects on the excitation of the molecular gas, such as cosmic ray ionization rates and mechanical heating through turbulence dissipation.

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.

We present a high-resolution study of the cold molecular gas as traced by CO(1-0) in the unlensed z ∼3.4 submillimeter galaxy SMM J13120+4242, using multiconfiguration observations with the Karl G. Jansky Very Large Array (JVLA). The gas reservoir, imaged on 0.″39 (∼3 kpc) scales, is resolved into two components separated by ∼11 kpc with a total extent of 16 ± 3 kpc. Despite the large spatial extent of the reservoir, the observations show a CO(1-0) FWHM linewidth of only 267 ± 64 km s-1. We derive a revised line luminosity of LCO(1-0)′ = (10 ± 3) × 1010 K km s-1 pc2 and a molecular gas mass of M gas = (13 ± 3)× 1010 (α CO/1) M ⊙. Despite the presence of a velocity gradient (consistent with previous resolved CO(6-5) imaging), the CO(1-0) imaging shows evidence for significant turbulent motions that are preventing the gas from fully settling into a disk. The system likely represents a merger in an advanced stage. Although the dynamical mass is highly uncertain, we use it to place an upper limit on the CO-to-H2 mass conversion factor α CO of 1.4. We revisit the SED fitting, finding that this galaxy lies on the very massive end of the main sequence at z = 3.4. Based on the low gas fraction, short gas depletion time, and evidence for a central AGN, we propose that SMM J13120 is in a rapid transitional phase between a merger-driven starburst and an unobscured quasar. The case of SMM J13120 highlights how mergers may drive important physical changes in galaxies without pushing them off the main sequence.