Characterising the relationship between dense gas and star formation is critical for understanding the assembly of galaxies throughout cosmic history. However, due to the faintness of standard dense-gas tracers-HCN, HCO⁺, and HNC-dense gas in high-redshift galaxies remains largely unexplored. We present ALMA and NOEMA observations targeting HCN/HCO⁺/HNC (3–2) and (4–3) emission lines in 11 (mostly) gravitationally lensed dusty star-forming galaxies (DSFGs) at redshift z = 1.6 − 3.2. We detect at least one line in 10 out of 11 galaxies. Altogether, we detect 34 dense-gas transitions, more than quadrupling the number of extant high-redshift detections. Additionally, in two targets, we detect lower-abundance CO isotopologues ¹³CO and C¹⁸O, as well as CN emission. We derive excitation coefficients for HCN, HCO⁺, and HNC in DSFGs, finding them to be systematically higher than those in nearby luminous infrared galaxies. Assuming the canonical dense-mass conversion factor (α _HCN = 10), we find that DSFGs have shorter dense-gas depletion times (median 23 Myr) than nearby galaxies (≈60 Myr), with a star-forming efficiency per free-fall time of 1−2%, a factor of a few higher than in local galaxies. We find a wide range of dense-gas fractions, with HCN/CO ratios ranging between 0.01 and 0.15. Finally, we put the first constraints on the redshift evolution of the cosmic dense-gas density, which increases by a factor of 7 ± 4 between z = 0 and z = 2.5, consistent with the evolution of the cosmic molecular-gas density.

Determining the nature of emission processes at the heart of quasars is critical for understanding environments of supermassive black holes. One of the key open questions is the origin of centimetre- to millimetre-wave emission from radio-quiet quasars. The proposed mechanisms range from central star formation to dusty torus, low-power jets, or emission from the accretion-disc corona. Distinguishing between these scenarios requires probing spatial scales of ≤0.01 pc, beyond the reach of any current millimetre-wave telescope. Fortunately, in gravitationally lensed quasars, compact millimetre-wave emission might be microlensed by stars in the foreground galaxy, providing strong constraints on the source size. We report a striking change in rest-frame 1.3 mm flux ratios in RXJ1131−1231, a quadruply lensed quasar at z = 0.658 observed by the Atacama Large Millimeter/submillimeter Array (ALMA) in 2015 and 2020. Over this period, the flux ratios between the three quasar images, A, B, and C, changed by a factor of 1.6 (A/B) and 3.0 (A/C). The observed flux-ratio variability is consistent with the microlensing of a compact source with a half-light radius of ≤50 astronomical units. The compactness of the source leaves coronal emission as the most likely scenario. Furthermore, the inferred millimetre-wave and X-ray luminosities follow the Güdel-Benz relationship for stellar coronae. These observations represent the first unambiguous evidence that coronae are the dominant mechanism for centimetre- to millimetre-wave emission in radio-quiet quasars.

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.

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 report the highest-redshift detection of [O i] 63 m from a luminous quasar, J20540005, at based on the Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 9 observations. The [O i] 63 m line luminosity is, corresponding to the [O i] 63 m-to-far-infrared luminosity ratio of 6.7, which is consistent with the value obtained in the local Universe. Remarkably, [O i] 63 m is as bright as [C ii] 158 m, resulting in the [O i]-to-[C ii] line luminosity ratio of. Based on a careful comparison of the luminosity ratios of [O i] 63 m, [C ii] 158 m, and dust continuum emission to models of photodissociation regions, we find that J20540005 has a gas density and an incident far-ultraviolet radiation field of, showing that [O i] 63 m serves as an important coolant of the dense and warm gas in J20540005. A close examination of the [O i] and [C ii] line profiles suggests that the [O i] line may be partially self-absorbed; however, deeper observations are needed to verify this conclusion. Regardless, the gas density and incident radiation field are in broad agreement with the values obtained in nearby star-forming galaxies and objects with [O i] 63 m observations at -3 with the Herschel Space Observatory. These results demonstrate the power of ALMA high-frequency observations targeting [O i] 63 m to examine the properties of photodissociation regions in high-redshift galaxies.

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.

Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources-suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z = 1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (λ_rest = 1.3 mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_⊙, with star formation rates of 25-160 M_⊙ yr^-1 for the host galaxies. These QSOs have lower L′_CO than star-forming galaxies with the same LIR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90(α_CO/4) Myr. We find molecular gas masses in the range ≤2-40 × 10⁹(α_CO/4) M_⊙, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z ~ 1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population.

During the most active period of star formation in galaxies, which occurs in the redshift range 1 < z < 3, strong bursts of star formation result in significant quantities of dust, which obscures new stars being formed as their UV/optical light is absorbed and then re-emitted in the infrared, which redshifts into the mm/sub-mm bands for these early times. To get a complete picture of the high- z galaxy population, we need to survey a large patch of the sky in the sub-mm with sufficient angular resolution to resolve all galaxies, but we also need the depth to fully sample their cosmic evolution, and therefore obtain their redshifts using direct mm spectroscopy with a very wide frequency coverage. This requires a large single-dish sub-mm telescope with fast mapping speeds at high sensitivity and angular resolution, a large bandwidth with good spectral resolution and multiplex spectroscopic capabilities. The proposed 50-m Atacama Large Aperture Submillimeter Telescope (AtLAST) will deliver these specifications. We discuss how AtLAST allows us to study the whole population of high-z galaxies, including the dusty star-forming ones which can only be detected and studied in the sub-mm, and obtain a wealth of information for each of these up to z ∼ 7: gas content, cooling budget, star formation rate, dust mass, and dust temperature. We present worked examples of surveys that AtLAST can perform, both deep and wide, and also focused on galaxies in proto-clusters. In addition we show how such surveys with AtLAST can measure the growth rate f σ ₈ and the Hubble constant with high accuracy, and demonstrate the power of the line-intensity mapping method in the mm/sub-mm wavebands to constrain the cosmic expansion history at high redshifts, as good examples of what can uniquely be done by AtLAST in this research field.

We present high-resolution (arcsec = 710 pc) Atacama Large Millimetre/submillimetre Array [C ii] 158 m and dust continuum follow-up observations of REBELS-25, a [C ii]-luminous () galaxy at redshift. These high-resolution, high signal-To-noise observations allow us to study the sub-kpc morphology and kinematics of this massive () star-forming (SFR) galaxy in the Epoch of Reionization. By modelling the kinematics with BAROLO, we find it has a low-velocity dispersion (km s) and a high ratio of ordered-To-random motion (), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [C ii] distribution is well fit by a near-exponential disc model, with a Sersic index, n, of, and we see tentative evidence of more complex non-Axisymmetric structures suggestive of a bar in the [C ii] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high-redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion, and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionized gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the big bang.

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.

Our knowledge of galaxy formation and evolution has incredibly progressed through multi-wavelength observational constraints of the interstellar medium (ISM) of galaxies at all cosmic epochs. However, little is known about the physical properties of the more diffuse and lower surface brightness reservoir of gas and dust that extends beyond ISM scales and fills dark matter haloes of galaxies up to their virial radii, the circumgalactic medium (CGM). New theoretical studies increasingly stress the relevance of the latter for understanding the feedback and feeding mechanisms that shape galaxies across cosmic times, whose cumulative effects leave clear imprints into the CGM. Recent studies are showing that a – so far unconstrained – fraction of the CGM mass may reside in the cold (T < 10 ⁴ K) molecular and atomic phase, especially in high-redshift dense environments. These gas phases, together with the warmer ionised phase, can be studied in galaxies from z ∼ 0 to z ∼ 10 through bright far-infrared and sub-millimeter emission lines such as [C ii] 158 µm, [O iii] 88 µm, [C I] 609 µm, [C i] 370 µm, and the rotational transitions of CO. Imaging such hidden cold CGM can lead to a breakthrough in galaxy evolution studies but requires a new facility with the specifications of the proposed Atacama Large Aperture Submillimeter Telescope (AtLAST). In this paper, we use theoretical and empirical arguments to motivate future ambitious CGM observations with AtLAST and describe the technical requirements needed for the telescope and its instrumentation to perform such science.

Feedback and outflows in galaxies that are associated with a quasar phase are expected to be pivotal in quenching the most massive galaxies. However, observations targeting the molecular outflow phase, which dominates both the mass and momentum and removes the immediate fuel for star formation, are limited in high-z QSO hosts. Massive quiescent galaxies found at z ∼ 4 are predicted to have quenched star formation already by z ∼ 5 and undergone their most intense growth at z > 6. Here, we present two Atacama Large Millimeter/submillimeter Array (ALMA) detections of molecular outflows, traced by blueshifted absorption of the OH 119 μm doublet, from a sample of three z > 6 infrared luminous QSO hosts: J2310+1855 and P183+05. OH 119 μm is also detected in emission from P183+05, and tentatively in the third source: P036+03. Using similar assumptions as for high-z dusty star-forming galaxy outflows, we find that our QSOs drive molecular outflows with comparable mass outflow rates, which are comparably energetic except for J2310+1855's significantly lower outflow energy flux. We do not find evidence, nor require additional input from the central active galactic nucleus (AGN) to drive the molecular outflow in J2310+1855, but we cannot rule out an AGN contribution in P183+05 if a significant AGN contribution to L _FIR is assumed and/or if the outflow covering fraction is high (≥53%), which evidence from the literature suggests is unlikely in these sources. Differences observed in the blueshifted absorption spectral properties may instead be caused by the QSO hosts’ more compact dust continuums, limiting observations to lower altitude and more central regions of the outflow.

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 deep ALMA Band 3 observations of the HCN, HCO⁺, and HNC(4-3) emission in SDP.81, a well-studied z = 3.042; strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z ≥ 1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC(4-3) lines are not detected, the HCO⁺(4-3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO⁺/HCN intensity ratio of ≥2.2. Based on the modelling of the photodissociation region, the most likely explanation for the elevated HCO⁺/HCN ratio is that SDP.81 has low mechanical heating, making up less than 10% of the total energy budget, along with a sub-solar metallicity of Z ≈ 0.5 Z_⊙. While such conditions might not be representative of the general population of high-redshift dusty galaxies, a lower-than-solar metallicity might significantly impact gas masses inferred from CO observations. In addition, we report the detection of CO(0-1) absorption from the foreground lensing galaxy and CO(1-0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the south-east.

In this erratum, we correct a mistake in the derivation of OH 119 µm equivalent width in two sources: J2310+1855 and P183+05. Consequently, we also correct the molecular gas outflow mass, mass outflow rate (MOFR), momentum flux, kinetic energy flux, and depletion times as the derivation of these values involves the equivalent width. We provide an updated version of Table 3 from the published article, and of Figures 3, 5, and 6. We no longer find significantly larger OH 119 µm absorption EWs in our unobscured QSO sources with respect to the high-z DSFGs from the literature Spilker et al. (2020a, 2020b). Furthermore, the MOFR, momentum flux, and kinetic energy of the molecular outflows in J2310+1855 and P183+05, as traced by the blueshifted OH 119 µm absorption, are now all significantly offset to lower values with respect to the trends with far-infrared (FIR) luminosity seen in high-z DSFGs (Spilker et al. 2020a, 2020b). Even with an assumed 50% contribution to the FIR luminosity from the central active nucleus, both galaxies appear to have suppressed outflow properties. The star formation rate (SFR) exceeds the MOFR in both sources and is therefore the dominant mechanism responsible for depleting the molecular gas reservoir in these systems. The original conclusion of the published article is therefore unchanged and even reinforced. We would like to thank Tom Bakx for bringing this error to our attention.

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [C ii] and ∼158 continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; LIR = L⊙) at z = 7.31, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR. This is about four times the SFR estimated from an extrapolated main sequence. We also infer a [C ii]-based molecular gas mass of, implying a molecular gas depletion time of Gyr. We observe a [C ii] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [C ii] emission at large positive velocities (∼500 km s-1), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200, which is consistent with expectations for a star-formation-driven outflow. Given its large stellar mass, SFR, and molecular gas reservoir ∼700 Myr after the big bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at z ∼4.

We present a feasibility study for the high-redshift galaxy part of the Science Verification Campaign with the 220–440 GHz deshima 2.0 integrated superconducting spectrometer on the ASTE telescope. The first version of the deshima 2.0 chip has been recently manufactured and tested in the lab. Based on these realistic performance measurements, we evaluate potential target samples and prospects for detecting the [CII] and CO emission lines. The planned observations comprise two distinct, but complementary objectives: (1) acquiring spectroscopic redshifts for dusty galaxies selected in far-infrared/mm-wave surveys; (2) multi-line observations to infer physical conditions in dusty galaxies.

Integrated superconducting spectrometer (ISS) technology will enable ultra-wideband, integral-field spectroscopy for (sub)millimeter-wave astronomy, in particular, for uncovering the dust-obscured cosmic star formation and galaxy evolution over cosmic time. Here, we present the development of DESHIMA 2.0, an ISS for ultra-wideband spectroscopy toward high-redshift galaxies. DESHIMA 2.0 is designed to observe the 220–440 GHz band in a single shot, corresponding to a redshift range of z = 3.3–7.6 for the ionized carbon emission ([C II] 158 μ m). The first-light experiment of DESHIMA 1.0, using the 332–377 GHz band, has shown an excellent agreement among the on-sky measurements, the laboratory measurements, and the design. As a successor to DESHIMA 1.0, we plan the commissioning and the scientific observation campaign of DESHIMA 2.0 on the ASTE 10-m telescope in 2023. Ongoing upgrades for the full octave-bandwidth system include the wideband 347-channel chip design and the wideband quasi-optical system. For efficient measurements, we also develop the observation strategy using the mechanical fast sky-position chopper and the sky-noise removal technique based on a novel data-scientific approach. In the paper, we show the recent status of the upgrades and the plans for the scientific observation campaign.

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.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4-3) line emitting galaxies in 17 quasar fields at z ∼4 aimed at performing the first systematic search of dusty galaxies in high-z quasar environments. Our blind search of galaxies around the quasars results in five CO emitters with S/N ≥ 5.6 within a projected radius of R ∼2 1.5 h -1 cMpc and a velocity range of δv = ±1000 km s-1 around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity of 17.6-7.6+11.9 in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross-correlation length of r0,QG=8.37-2.04+2.42h-1 cMpc, assuming a fixed slope γ = 1.8. This contradicts the reported mild overdensities (x1.4) of Lyα emitters (LAEs) in the same fields at scales of R ∼2 7 h -1 cMpc, which are well described by a cross-correlation length 3.0-1.4+1.5 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at z ∼4, finding an autocorrelation length of r 0,CO = 3.14 ±1.71 h -1 cMpc (with γ = 1.8). Our work, together with the previous study of LAEs around quasars, traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multiwavelength studies, and highlighting important questions about galaxy properties in high-z dense environments.