NOEMA Research Articles

JWST’s PEARLS: Resolved study of the stellar and dust components in starburst galaxies at cosmic noon

Dusty star-forming galaxies (DSFGs) significantly contribute to the stellar buildup in galaxies during “cosmic noon,” the peak epoch of cosmic star formation. Major mergers and gas accretion are often invoked to explain DSFGs’ prodigious star formation rates (SFRs) and large stellar masses. We conducted a spatially resolved morphological analysis of the rest-frame ultraviolet/near-infrared (∼0.25–1.3 μm) emission in three DSFGs at z ≃ 2.5. Initially discovered as carbon monoxide (CO) emitters by NOrthern Extended Millimeter Array (NOEMA) observations of a bright (S350 μm = 111 ± 10 mJy) Herschel source, we observed them with the James Webb Space Telescope/NIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies’ stellar populations and dust distributions on scales of 250 pc. Spatial variations in stellar mass, SFR, and dust extinction are determined in resolved maps obtained through pixel-based spectral energy distribution fitting. The CO emitters are massive (Mstar ≃ (3 − 30)×1010 M⊙), dusty starburst galaxies with SFRs ranging from 340 to 2500 M⊙ yr−1, positioning them among the most active star-forming galaxies at 2 < z < 3. Notably, they belong to the ∼1.5% of the entire JWST population with extremely red colors. Their morphologies are disk like (Sérsic index n ≃ 1), with effective radii of 2.0–4.4 kpc, and exhibit substructures such as clumps and spiral arms. The galaxies have dust extinctions up to AV = 5–7 mag extending over several kiloparsecs with asymmetric distributions that include off-center regions resembling bent spiral arms and clumps. The near-infrared dust-attenuation curve in these sources deviates from standard laws, possibly implying different dust–star geometries or dust grain properties than commonly assumed in starburst galaxies. The proximity (< 5″) of galaxies with consistent redshifts, strong color gradients, an overall disturbed appearance, asymmetric dust obscuration, and widespread star formation collectively favor interactions (minor mergers and flybys) as the mechanism driving the CO galaxies’ exceptional SFRs. The galaxies’ large masses and rich environment hint at membership in two proto-structures, as initially inferred from their association with a Planck-selected high-z source.

NOEMA formIng Cluster survEy (NICE): Characterizing eight massive galaxy groups at 1.5 < z < 4 in the COSMOS field

The NOrthern Extended Millimeter Array (NOEMA) formIng Cluster survEy (NICE) is a NOEMA large programme targeting 69 massive galaxy group candidates at z > 2 over six deep fields with a total area of 46 deg2. Here we report the spectroscopic confirmation of eight massive galaxy groups at redshifts 1.65 ≤ z ≤ 3.61 in the Cosmic Evolution Survey (COSMOS) field. Homogeneously selected as significant overdensities of red IRAC sources that have red Herschel colours, four groups in this sample are confirmed by CO and [CI] line detections of multiple sources with NOEMA 3 mm observations, three are confirmed with Atacama Large Millimeter Array (ALMA) observations, and one is confirmed by Hα emission from Subaru/FMOS spectroscopy. Using rich ancillary data in the far-infrared and sub-millimetre, we constructed the integrated far-infrared spectral energy distributions for the eight groups, obtaining a total infrared star formation rate (SFR) of 260–1300 M⊙ yr−1. We adopted six methods for estimating the dark matter masses of the eight groups, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass densities. We find that the radial stellar mass densities of the eight groups are consistent with a NFW profile, supporting the idea that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are log(Mh/M⊙) = 12.8 − 13.7 with a general uncertainty of 0.3 dex. Based on the halo mass estimates, we derived baryonic accretion rates (BARs) of (1 − 8)×103 M⊙/yr for this sample. Together with massive groups in the literature, we find a quasi-linear correlation between the integrated SFR/BAR ratio and the theoretical halo mass limit for cold streams, Mstream/Mh, with SFR/BAR = 10−0.46 ± 0.22(Mstream/Mh)0.71 ± 0.16 with a scatter of 0.40 dex. Furthermore, we compared the halo masses and the stellar masses with simulations, and find that the halo masses of all structures are consistent with those of progenitors of Mh(z = 0) > 1014 M⊙ galaxy clusters, and that the most massive central galaxies have stellar masses consistent with those of the brightest cluster galaxy progenitors in the TNG300 simulation. Above all, the results strongly suggest that these massive structures are in the process of forming massive galaxy clusters via baryonic and dark matter accretion.

PRODIGE - Planet-forming disks in Taurus with NOEMA. II. Modeling the CO (2-1) isotopologue emission of the Class II T Tauri disks in Taurus

To understand how planets form in protoplanetary disks, it is necessary to characterize their gas and dust distribution and masses. This requires a combination of high-resolution dust continuum and molecular line interferometric observations, coupled with advanced theoretical models of protoplanetary disk physics, chemical composition, and radiative transfer. We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at $0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution \, PIs: P. Caselli Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use these data to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C$ $O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. The best-fit disk gas masses range between 0.005 and $0.04 M_ sun $. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements via, for example, modeling of the wavelength dependence of the dust continuum emission radii, and HD and CO isotopologue emission.

PRODIGE – planet-forming disks in Taurus with NOEMA

Context. The physics and chemistry of planet-forming disks are far from being fully understood. To make further progress, both broad line surveys and observations of individual tracers in a statistically significant number of disks are required. Aims. Our aim is to perform a line survey of eight planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Methods. Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from > 20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210–280 GHz with a 1σ rms sensitivity of ~8–12 mJy beam−1 at a 0.9″ and 0.3 km s−1 resolution will be completed in 2024. The uv visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. Results. In this first paper, the combined 12CO, 13CO, and C18O J = 2–1 data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrained dynamical stellar masses, geometries, temperatures, the CO column densities, and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17–37 K, and decrease radially with the power-law exponent q ~ 0.05–0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2–3.1. The gas masses estimated from 13CO (2–1) are ~0.001–0.2 M⊙. Conclusions. Using NOEMA, we confirm the presence of temperature gradients in our disk sample. The best-fit CO column densities point to severe CO freeze-out in these disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions. The high 13CO/C18O ratio of ~23 in DM Tau could be indicative of strong selective photodissociation of C18O in this disk.

Noema formIng Cluster survEy (NICE): Discovery of a starbursting galaxy group with a radio-luminous core at z = 3.95

The study of distant galaxy groups and clusters at the peak epoch of star formation is limited by the lack of a statistically and homogeneously selected and spectroscopically confirmed sample. Recent discoveries of concentrated starburst activities in cluster cores have opened a new window to hunt for these structures based on their integrated IR luminosities. Here, we carry out a large NOEMA (NOrthern Extended Millimeter Array) program targeting a statistical sample of infrared-luminous sources associated with overdensities of massive galaxies at z > 2, the Noema formIng Cluster survEy (NICE). We present the first result from the ongoing NICE survey, a compact group at z = 3.95 in the Lockman Hole field (LH-SBC3), confirmed via four massive (M⋆ ≳ 1010.5 M⊙) galaxies detected in the CO(4–3) and [CI](1–0) lines. The four CO-detected members of LH-SBC3 are distributed over a 180 kpc physical scale and the entire structure has an estimated halo mass of ∼1013 M⊙ and total star formation rate of ∼4000 M⊙ yr−1. In addition, the most massive galaxy hosts a radio-loud active galactic nucleus with L1.4 GHz, rest = 3.0 × 1025 W Hz−1. The discovery of LH-SBC3 demonstrates the feasibility of our method to efficiently identify high-z compact groups or cluster cores undergoing formation. The existence of these starbursting cluster cores up to z ∼ 4 provides critical insights into the mass assembly history of the central massive galaxies in clusters.

The NOEMA observations of GN-z11: constraining the neutral interstellar medium and dust formation in the heart of cosmic reionization at z = 10.6

ABSTRACT We present results of dust continuum and [C ii]$\, 158\, {\rm \mu m}$ emission line observations of a remarkably UV luminous (MUV = −21.6) galaxy at z = 10.603: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $\lambda _{\rm rest}=160\, {\rm \mu m}$ continuum $1\, \sigma$ sensitivity of $13.0\, \rm {\mu Jy \, beam}^{ -1}$ and a [C ii] emission line $1\, \sigma$ sensitivity of $31\, \rm {mJy\, beam^{ -1}\, km \, s}^{ -1}$ using $50\, \rm {km \, s}^{ -1}$ binning with a $\sim 2\, {\rm arcsec}$ synthesized beam. Neither dust continuum nor [C ii]$\, 158\, {\rm \mu m}$ line emission are detected at the expected frequency of ν[C ii]$= 163.791\, \rm {GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in LIR nor L[C ii], with a dust mass $3\, \sigma$ limit of ${\rm log}\, (M_{\rm dust}/{\rm {\rm M}_{\odot }}) \, \lt\, 6.5-6.9$ and with a [C ii] based molecular gas mass $3\, \sigma$ limit of log (Mmol, [C ii]$/{\rm {\rm M}_{\odot }}) \, \lt \, 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue colour in the UV continuum of GN-z11 (βUV = −2.4), and found that ≳3 × deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high-redshift galaxies to constrain multiple phases in the interstellar medium.

The z-GAL survey

The z-GAL survey observed 137 bright Herschel-selected galaxies with the IRAM Northern Extended Millimeter Array (NOEMA) to measure their spectroscopic redshift. It detected several multiple sources and secured the spectroscopic redshift of a total of 165 individual galaxies over the range 0.8 < z < 6.5. The wealth of information provided by the NOEMA spectra allowed us to study the physical properties of these galaxies. By modelling their spectral energy distributions, we determined their dust mass, Mdust, temperature, Tdust, and emissivity index, βdust and found a significant anti-correlation between Tdust and βdust. Using the detected 12CO lines, we derived the molecular gas mass of the z-GAL sources. Combining the molecular lines and the dust continuum emission, we computed their gas depletion timescales, τdep: the z-GAL sample consists mainly of powerful starbursts with τdep in the range between 0.1 and 1.0 Gyr, while only ~ 25% belongs to the ‘main sequence’ of star forming galaxies.

SUNRISE: The rich molecular inventory of high-redshift dusty galaxies revealed by broadband spectral line surveys

Understanding the nature of high-redshift dusty galaxies requires a comprehensive view of their interstellar medium (ISM) and molecular complexity. However, the molecular ISM at high redshifts is commonly studied using only a few species beyond 12C16O, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOrthern Extended Millimeter Array (NOEMA) targeting two strongly lensed dusty galaxies observed when the Universe was less than 1.8 Gyr old: APM 08279+5255, a quasar at redshift z = 3.911, and NCv1.143 (H-ATLAS J125632.7+233625), a z = 3.565 starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330 to 550 GHz for both galaxies. We report the detection of 38 and 25 emission lines in APM 08279+5255 and NCv1.143, respectively. These lines originate from 17 species, namely CO, 13CO, C18O, CN, CCH, HCN, HCO+, HNC, CS, C34S, H2O, H3O+, NO, N2H+, CH, c-C3H2, and the vibrationally excited HCN and neutral carbon. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the analysis of the molecular gas excitation, we find that the physical properties and the chemical imprints of the ISM are different: the molecular gas is more excited in APM 08279+5255, which exhibits higher molecular gas temperatures and densities compared to NCv1.143; the molecular abundances in APM 08279+5255 are akin to the values of local active galactic nuclei (AGN), showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or the X-ray-dominated regions, while NCv1.143 more closely resembles local starburst galaxies. The most significant differences between the two sources are found in H2O: the 448 GHz ortho-H2O(423 − 330) line is significantly brighter in APM 08279+5255, which is likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that, at such high column densities, far-ultraviolet radiation is less important in regulating the ISM, while cosmic rays (and/or X-rays and shocks) are the key players in shaping the molecular abundances and the initial conditions of star formation. Both our observed CO isotopologs line ratios and the derived extreme ISM conditions (high gas temperatures, densities, and cosmic-ray ionization rates) suggest the presence of a top-heavy stellar initial mass function. From the ∼330–550 GHz continuum, we also find evidence of nonthermal millimeter flux excess in APM 08279+5255 that might be related to the central supermassive black hole. Such deep spectral line surveys open a new window into the physics and chemistry of the ISM and the radiation field of galaxies in the early Universe.

A high HDO/H2O ratio in the Class I protostar L1551 IRS5

Context. Water is a very abundant molecule in star-forming regions. Its deuterium fractionation provides an important tool for understanding its formation and evolution during the star and planet formation processes. While the HDO/H2O abundance ratio has been determined toward several young Class 0 protostars and comets, the number of studies toward Class I protostars is limited. Aims. Our aim is to study the water deuteration toward the Class I protostellar binary L1551 IRS5 and to investigate the effect of evolutionary stage and environment on variations in the water D/H ratio. Methods. Observations were carried out toward L1551 IRS5 using the NOrthern Extended Millimeter Array (NOEMA) interferometer. The HDO 31, 2–22, 1 transition at 225.9 GHz and the H218O 31, 3–22, 0 transition at 203.4 GHz were covered with a spatial resolution of 0.5″× 0.8″, while the HDO 42, 2–42, 3 transition at 143.7 GHz was observed with a resolution of 2.0″ × 2.5″. We constrained the water D/H ratio using both local thermodynamic equilibrium (LTE) and non-LTE models. Results. The three transitions are detected. The line profiles display two peaks, one at ∼6 km s−1 and one at ∼9 km s−1. We derive an HDO/H2O ratio of (2.1 ± 0.8) × 10−3 for the redshifted component and a lower limit of > 0.3 × 10−3 for the blueshifted component. This lower limit is due to the blending of the blueshifted H218O component with redshifted CH3OCH3 emission. Conclusions. The HDO/H2O ratio in L1551 IRS5 is similar to the values in Class 0 isolated sources and in the disk of the Class I protostar V883 Ori, while it is significantly higher than in the previously studied clustered Class 0 sources and the comets. This result suggests that the chemistry of protostars that belong to molecular clouds with relatively low source densities, such as L1551, share more similarities with the isolated sources than the protostars of very dense clusters. If the HDO/H2O ratios in Class 0 protostars with few sources around are comparable to those found to date in isolated Class 0 objects, it would mean that there is little water reprocessing from the Class 0 to Class I protostellar stage.

Unveiling the gravitationally unstable disc of a massive star-forming galaxy using NOEMA and MUSE

ABSTRACT Using new high-resolution data of CO (2–1), H$\alpha$ and H$\beta$ obtained with the Northern Extended Millimeter Array (NOEMA) and the Multi-Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope, we have performed a Toomre Q disc stability analysis and studied star formation, gas depletion times and other environmental parameters on sub-kpc scales within the z ∼ 0 galaxy SDSS J125013.84+073444.5 (LARS 8). The galaxy hosts a massive, clumpy disc and is a proto-typical analogue of main-sequence galaxies at z ∼ 1 − 2. We show that the massive (molecular) clumps in LARS 8 are the result of an extremely gravitationally unstable gas disc, with large scale instabilities found across the whole extent of the rotating disc, with only the innermost 500 pc being stabilized by its bulge-like structure. The radial profiles further reveal that – contrary to typical disc galaxies – the molecular gas depletion time decreases from more than 1 Gyr in the centre to less than ∼100 Myr in the outskirts of the disc, supporting the findings of a Toomre-unstable disc. We further identified and analysed 12 individual massive molecular clumps. They are virialized and follow the mass–size relation, indicating that on local (cloud/clump) scales the stars form with efficiencies comparable to those in Milky Way clouds. The observed high star formation rate must thus be the result of triggering of cloud/clump formation over large scales due to disc instability. Our study provides evidence that ‘in-situ’ massive clump formation (as also observed at high redshifts) is very efficiently induced by large-scale instabilities.

New constraints on the presence of debris disks around G 196-3 B and VHS J125601.92–125723.9 b

Context. The existence of warm (protoplanetary) disks around very young isolated planetary and brown dwarf mass objects is known based on near- and mid-infrared flux excesses and millimeter observations. These disks may later evolve into debris disks or rings, although none have been observed or confirmed so far. Little is known about circum(sub)stellar and debris disks around substellar objects. Aims. We aim to investigate the presence of debris disks around two of the closest (~20 pc), young substellar companions, namely G196-3 B and VHS J125601.92–125723.9 b (VHS J1256–1257 b), whose masses straddle the borderline between planets and brown dwarfs. Both are companions at wide orbits (≥100 au) of M-type dwarfs and their ages (50–100 Myr and 150–300 Myr, respectively) are thought to be adequate for the detection of second-generation disks. Methods. We obtained deep images of G196-3 B and VHS J1256–1257 b with the NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm. These data were combined with recently published Atacama Large Millimeter Array (ALMA) and Very Large Array (VLA) data of VHS J1256–1257 b at 0.87 mm and 0.9 cm, respectively. Results. Neither G196-3 B nor VHS J1256–1257 b were detected in the NOEMA, ALMA, and VLA data. At 1.3 mm, we imposed flux upper limits of 0.108 mJy (G196-3 B) and 0.153 mJy (VHS J1256–1257 b) with a 3-σ confidence. Using the flux upper limits at the millimeter and radio wavelength regimes, we derived maximum values of 1.38×10−2 MEarth and 5.46 × 10−3 MEarth for the mass of any cold dust that might be surrounding G196-3 B and VHS J1256–1257 b, respectively. Conclusions. We put our results in the context of other deep millimeter observations of free-floating and companion objects with substellar masses smaller than 20 MJup and ages between approximately one and a few hundred million years. Only two very young (2–5.4 Myr) objects are detected out of a few tens of them. This implies that the disks around these very low-mass objects must have small masses, and possibly reduced sizes, in agreement with findings by other groups. If debris disks around substellar objects scale down (in mass and size) in a similar manner as protoplanetary disks do, millimeter observations of moderately young brown dwarfs and planets must be at least two orders of magnitude deeper to be able to detect and characterize their surrounding debris disks.

Witnessing the fragmentation of a filament into prestellar cores in Orion B/NGC 2024

Context. Recent Herschel observations of nearby molecular clouds have shown that filamentary structures are ubiquitous and that most prestellar cores form in dense filaments. Probing the detailed density and velocity structure of molecular filaments is therefore crucial for improving our observational understanding of the star formation process. Aims. We aim to characterize both the density and the velocity field of a typical molecular filament in the process of fragmenting into cores. Methods. We mapped a portion of the NGC 2024 region in the Orion B molecular cloud with the Nobeyama 45m telescope, in the 12CO (J = 1–0), 13CO (J = 1–0), C18O (J = 1–0), and H13CO+ (J = 1–0) lines, and the southwestern part of NGC 2024, corresponding to the NGC 2024S filament, with the NOrthern Extended Millimeter Array (NOEMA) interferometer in H13CO+ (J = 1–0). Results. The maps of 13CO, C18O, and H13CO+ emission trace at least part of the filamentary structure seen in the 8″ resolution ArTéMiS+Herschel data. The median radial column density profile of the NGC 2024S filament as derived from ArTéMiS+Herschel dust emission data is well fitted by a Plummer profile with a half-power diameter DHPPlummer=0.081±0.014 pc, which is similar to the findings of previous studies of nearby molecular filaments with Herschel. On the other hand, the half-power diameters of NGC 2024S as measured from the Nobeyama 13CO and C18O data are broader than, and the half-power diameter derived from the H13CO+ data narrower than the filament diameter measured with Herschel. These results suggest that the 13CO and C18O data trace only the (low-density) outer part of the Herschel filament and the H13CO+ data only the (dense) inner part. We identify four cores in the portion of the Herschel map covered by NOEMA and find that each Herschel core corresponds to a single core detected in the combined NOEMA+45m H13CO+ data cube. The Nobeyama H13CO+ centroid velocity map reveals velocity gradients along both the major and the minor axis of the NGC 2024S filament, as well as velocity oscillations with a period λ ~0.2 pc along the major axis. Comparison between the centroid velocity and the column density distribution shows a tentative λ/4 phase shift in H13CO+ or C18O. This λ/4 shift is not simultaneously observed for all cores in any single tracer but is tentatively seen for each core in either H13CO+ or C18O. The difference between the H13CO+ and C18O velocity patterns may arise from differences in the range of densities probed by H13CO+ and C18O. We produced a toy model that takes into account the three velocity-field components: a transverse velocity gradient, a longitudinal velocity gradient, and a longitudinal oscillation mode caused by fragmentation. Examination of synthetic data shows that the longitudinal oscillation component produces an oscillation pattern in the velocity structure function of the model. Since the velocity structure function of the Nobeyama H13CO+ centroid velocity data does show an oscillation pattern, we suggest that our observations are partly tracing core-forming motions and the fragmentation of the NGC 2024S filament into cores. We also find that the mean core mass in NGC 2024S corresponds to the effective Bonnor-Ebert mass in the filament. This is consistent with a scenario in which higher-mass cores form in higher-line-mass filaments.

A NOEMA Molecular Line Scan of the Hubble Deep Field North: Improved Constraints on the CO Luminosity Functions and Cosmic Density of Molecular Gas

We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z ∼ 6 based on an 8.5 arcmin2 spectral scan survey at 3 mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distributions exploiting the extensive multiwavelength data for the HDF-N. We identify the seven highest-fidelity sources as CO emitters at 1 < z < 6, including the well-known submillimeter galaxy HDF 850.1 at z = 5.18. Four high-fidelity 3 mm continuum sources are found to be radio galaxies at z ≤ 1, plus HDF 850.1. We constrain the CO LFs in the HDF-N out to z ∼ 6, including a first measurement of the CO(5–4) LF at 〈z〉 = 5.0. The relatively large area and depth of the NOEMA HDF-N survey extends the existing LFs at 1 < z < 4 above the knee, yielding a somewhat lower density by 0.15–0.4 dex at the overlap region for the CO(2–1) and CO(3–2) transitions, attributed to cosmic variance. We perform a joint analysis of the CO LFs in the HDF-N and Hubble Ultra Deep Field from ASPECS, finding that they can be well described by a single Schechter function. The evolution of the cosmic molecular gas density from a joint analysis is in good agreement with earlier determinations. This implies that the impact of cosmic field-to-field variance on the measurements is consistent with previous estimates, adding to the challenges for simulations that model galaxies from first principles.

Young Radio Sources Expanding in Gas-Rich ISM: Using Cold Molecular Gas to Trace Their Impact

We present an overview of the results obtained from the study of the resolved distribution of molecular gas around eight young (≲106yr), peaked-spectrum radio galaxies. Tracing the distribution and kinematics of the gas around these radio sources allows us to trace the interplay between the jets and the surrounding medium. For three of these sources, we present new CO(1-0) observations, obtained with the Northern Extended Millimeter Array (NOEMA) with arcsecond resolution. In two of these targets, we also detected CN lines, both in emission and absorption. Combining the new observations with already published data, we discuss the main results obtained. Although we found that a large fraction of the cold molecular gas was distributed in disc-like rotating structures, in the vast majority of the sources, high turbulence and deviations from purely quiescent gas (including outflows) were observed in the region co-spatial with the radio continuum emission. This suggests the presence of an interaction between radio plasma and cold molecular gas. In particular, we found that newly born and young radio jets, even those with low power i.e., Pjet &lt; 1045 erg s−1), are able to drive massive outflows of cold, molecular gas. The outflows are, however, limited to the sub-kpc regions and likely short lived. On larger scales (a few kpc), we observed cases where the molecular gas appears to avoid the radio lobes and, instead, wraps around them. The results suggest the presence of an evolutionary sequence, which is consistent with previous simulations, where the type of impact of the radio plasma changes as the jet expands, going from a direct jet-cloud interaction able to drive gas outflows on sub-kpc scales to a more gentle pushing aside of the gas, increasing its turbulence and likely limiting its cooling on kpc scales. This effect can be mediated by the cocoon of shocked gas inflated by the jet–cloud interactions. Building larger samples of young and evolved radio sources for observation at a similar depth and spatial resolution to test this scenario is now needed and may be possible thanks to more data becoming available in the growing public archives.

Observations of comet C/2020 F3 (NEOWISE) with IRAM telescopes

We present the results of millimetre-wave spectroscopic and continuum observations of the comet C/2020 F3 (NEOWISE) undertaken with the Institut de RadioAstronomie Millimétrique (IRAM) 30-m and the NOrthern Extended Millimeter Array (NOEMA) telescopes on 22, 25–27 July, and 7 August 2020. Production rates of HCN, HNC, CH3OH CS, H2CO, CH3CN, H2S, and CO were determined with upper limits on six other species. The comet shows abundances within the range observed for other comets. The CO abundance is low (3.2% relative to water), while H2S is relatively abundant (1.1% relative to water). The H2CO abundance shows a steep variation with heliocentric distance, possibly related to a distributed production from the dust or macro-molecular source. The CH3OH and H2S production rates show a slower decrease post-perihelion than water. There was no detection of the nucleus point source contribution based on the interferometric map of the continuum (implying a size of r &lt; 4.7 km), but this yielded an estimate of the dust production rate, leading to a relatively low dust-to-gas ratio of 0.7 ± 0.3 on 22.4 July 2020.

Unveiling the warm and dense ISM in z &gt; 6 quasar host galaxies via water vapor emission

Water vapor (H2O) is one of the brightest molecular emitters after carbon monoxide (CO) in galaxies with high infrared (IR) luminosity, allowing us to investigate the warm and dense phase of the interstellar medium (ISM) where star formation occurs. However, due to the complexity of its radiative spectrum, H2O is not frequently exploited as an ISM tracer in distant galaxies. Therefore, H2O studies of the warm and dense gas at high-z remain largely unexplored. In this work, we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward three z &gt; 6 IR-bright quasars J2310+1855, J1148+5251, and J0439+1634 targeted in their multiple para- and ortho-H2O transitions (312 − 303, 111 − 000, 220 − 211, and 422 − 413), as well as their far-IR (FIR) dust continuum. By combining our data with previous measurements from the literature, we estimated the dust masses and temperatures, continuum optical depths, IR luminosities, and star formation rates (SFR) from the FIR continuum. We modeled the H2O lines using the MOLPOP-CEP radiative transfer code, finding that water vapor lines in our quasar host galaxies are primarily excited in the warm, dense (with a gas kinetic temperature and density of Tkin = 50 K, nH2 ∼ 104.5 − 105 cm−3) molecular medium with a water vapor column density of NH2O ∼ 2 × 1017 − 3 × 1018 cm−3. High-J H2O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component at temperatures of Tdust ∼ 80 − 190 K that account for &lt; 5 − 10% of the total dust mass. In the case of J2310+1855, our analysis points to a relatively high value of the continuum optical depth at 100 μm (τ100 ∼ 1). Our results are in agreement with expectations based on the H2O spectral line energy distribution of local and high-z ultra-luminous IR galaxies and active galactic nuclei (AGN). The analysis of the Boltzmann diagrams highlights the interplay between collisions and IR pumping in populating the high H2O energy levels and it allows us to directly compare the excitation conditions in the targeted quasar host galaxies. In addition, the observations enable us to sample the high-luminosity part of the H2O–total-IR (TIR) luminosity relations (LH2O − LTIR). Overall, our results point to supralinear trends that suggest H2O–TIR relations are likely driven by IR pumping, rather than the mere co-spatiality between the FIR continuum- and line-emitting regions. The observed LH2O/LTIR ratios in our z &gt; 6 quasars do not show any strong deviations with respect to those measured in star-forming galaxies and AGN at lower redshifts. This supports the notion that H2O can be likely used to trace the star formation activity buried deep within the dense molecular clouds.

A tentative ∼1000 km s−1 offset between the [CII] 158 m and Lyα line emission in a star-forming galaxy at z = 7.2

GN-108036 is a star-forming galaxy at z = 7.21, and one of the most distant known sources in the Northern hemisphere. Based on observations from the NOrthern Extended Millimeter Array (NOEMA), here we report the tentative detection of the [C II] line at ≈4σ significance. The integrated [C II] line emission is spatially offset about ∼4 kpc from the rest-frame ultraviolet (UV) emission. The total [C II] luminosity (L[CII] = 2.7 × 108 L⊙) is consistent with the relation between [C II] luminosity and star formation rate (SFR) observed in nearby and high-z star forming galaxies. More interestingly, the [C II] line is blueshifted with respect to the Lyα line by 980 ± 10 km s−1. If confirmed, this corresponds to the largest velocity offset reported to date between the Lyα line and a non-resonant line at z ≳ 6. According to trends observed in other high redshift galaxies, the large Lyα velocity offset in GN-108036 is consistent with its low Lyα equivalent width and high UV absolute magnitude. Based on Lyα radiative transfer models of expanding shells, the large Lyα velocity offset in GN-108036 could be interpreted as the presence of a large column density of hydrogen gas, and/or an outflow with a velocity of vout ∼ ΔvLyα/2 ∼ 500 km s−1. We also report the 3σ detection of a potential galaxy companion located ∼30 kpc east of GN-108036, at a similar systemic velocity, and with no counterpart rest-frame UV emission.

The [CII] and FIR properties ofz&gt; 6 radio-loud quasars

There are only five radio-loud quasars currently known within 1 Gyr from the Big Bang (z &gt; 6) and the properties of their host galaxies have not been explored in detail. We present a NOrthern Extended Millimeter Array (NOEMA) survey of [CII] (158 μm) and underlying continuum emission of fourz &gt; 6 radio-loud quasars, revealing their diverse properties. J0309+2717 (z = 6.10) has a bright [CII] line and underlying continuum, implying that the host galaxy is a starburst with a star-formation rate SFR = 340–1200M⊙yr−1. J1429+5447 (z = 6.18) has a SFR = 520 − 870M⊙yr−1and its [CII] profile is consistent with two Gaussians, which could be interpreted as a galaxy merger. J1427+3312 (z = 6.12) has a moderate SFR = 30–90M⊙yr−1. Notably, this is a broad absorption line quasar and we searched for the presence of high-velocity outflows in the host galaxy. Although the NOEMA data tentatively reveal a broad component of the [CII] line as wide as ∼1400 km s−1, the sensitivity of our current data are not sufficient to confirm it. Finally, P172+18 (z = 6.82) is undetected in both [CII] and the continuum, implying a SFR &lt; 22–40M⊙yr−1. The broad range of SFRs is similar to what is observed in radio-quiet quasars at similar redshifts. If radio jets do not significantly contribute to both [CII] and IR luminosities, this suggests there is no feedback from the jet on the star formation in the host galaxy.

A Massive, Dusty, Hi Absorption–Selected Galaxy at z ≈ 2.46 Identified in a CO Emission Survey

We report a NOrthern Extended Millimeter Array (NOEMA) and Atacama Large Millimeter/submillimeter Array search for redshifted CO emission from the galaxies associated with seven high-metallicity ([M/H] ≥ −1.03) damped Lyα absorbers (DLAs) at z ≈ 1.64–2.51. Our observations yielded one new detection of CO(3–2) emission from a galaxy at z = 2.4604 using NOEMA, associated with the z = 2.4628 DLA toward QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of % and %, respectively, in the fields of DLAs with [M/H] > −0.3 and [M/H] < −0.3. Further, the Hi–selected galaxies associated with five DLAs with [M/H] > −0.3 all have high molecular gas masses, ≳5 × 1010 M ⊙. This indicates that the highest-metallicity DLAs at z ≈ 2 are associated with the most massive galaxies. The newly identified z ≈ 2.4604 Hi–selected galaxy, DLA0201+365g, has an impact parameter of ≈7 kpc to the QSO sightline, and an implied molecular gas mass of (5.04 ± 0.78) × 1010 × (α CO/4.36) × (r 31/0.55) M ⊙. Archival Hubble Space Telescope Wide Field and Planetary Camera 2 imaging covering the rest-frame near-ultraviolet (NUV) and far-ultraviolet (FUV) emission from this galaxy yield nondetections of rest-frame NUV and FUV emission, and a 5σ upper limit of 2.3 M ⊙ yr−1 on the unobscured star formation rate (SFR). The low NUV-based SFR estimate, despite the very high molecular gas mass, indicates that DLA0201+365g either is a very dusty galaxy, or has a molecular gas depletion time that is around 2 orders of magnitude larger than that of star-forming galaxies at similar redshifts.

SOLIS

Context.Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass ejection provides constraints on the mass accretion history and on the nature of the driving source.Aims.We characterize the time-variability of the mass-ejection phenomena at work in the class 0 protostellar phase in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history.Methods.Using the NOrthern Extended Millimeter Array (NOEMA) interferometer, we have observed the emission of the CO 2–1 and SONJ= 54–43rotational transitions at an angular resolution of 1.0″ (820 au) and 0.4″ (330 au), respectively, toward the intermediate-mass class 0 protostellar system Cep E.Results.The CO high-velocity jet emission reveals a central component of ≤400 au diameter associated with high-velocity molecular knots that is also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to be accelerated along the main axis over a length scaleδ0~ 700 au, while its diameter gradually increases up to several 1000 au at 2000 au from the protostar. The jet is fragmented into 18 knots of mass ~10−3M⊙, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km s−1close to the protostar. This is well below the jet terminal velocities in the northern (+ 65 km s−1) and southern (−125 km s−1) lobes. The knot interval distribution is approximately bimodal on a timescale of ~50–80 yr, which is close to the jet-driving protostar Cep E-A and ~150–20 yr at larger distances &gt;12″. The mass-loss rates derived from knot masses are steady overall, with values of 2.7 × 10−5M⊙yr−1and 8.9 × 10−6M⊙yr−1in the northern and southern lobe, respectively.Conclusions.The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet. This accounts for the higher mass-loss rate in the northern lobe. The jet dynamics are well accounted for by a simple precession model with a period of 2000 yr and a mass-ejection period of 55 yr.