skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Unveiling the warm and dense ISM in z  > 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-zremain largely unexplored. In this work, we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward threez > 6 IR-bright quasarsJ2310+1855,J1148+5251, andJ0439+1634targeted 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 ofTkin = 50 K,nH2 ∼ 104.5 − 105 cm−3) molecular medium with a water vapor column density ofNH2O ∼ 2 × 1017 − 3 × 1018 cm−3. High-JH2O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component at temperatures ofTdust ∼ 80 − 190 K that account for < 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-zultra-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 observedLH2O/LTIRratios in ourz > 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.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Astronomy and Astrophysics
Date Published:
Journal Name:
Astronomy & Astrophysics
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    We present a multiline survey of the interstellar medium (ISM) in two z  > 6 quasar host galaxies, PJ231−20 ( z  = 6.59) and PJ308−21 ( z  = 6.23), and their two companion galaxies. Observations were carried out using the Atacama Large (sub-)Millimeter Array (ALMA). We targeted 11 transitions including atomic fine-structure lines (FSLs) and molecular lines: [NII] 205 μm , [CI] 369 μm , CO ( J up  = 7, 10, 15, 16), H 2 O 3 12  − 2 21 , 3 21  − 3 12 , 3 03  − 2 12 , and the OH 163 μm doublet. The underlying far-infrared (FIR) continuum samples the Rayleigh-Jeans tail of the respective dust emission. By combining this information with our earlier ALMA [CII] 158 μm observations, we explored the effects of star formation and black hole feedback on the ISM of the galaxies using the CLOUDY radiative transfer models. We estimated dust masses, spectral indexes, IR luminosities, and star-formation rates from the FIR continuum. The analysis of the FSLs indicates that the [CII] 158 μm and [CI] 369 μm emission arises predominantly from the neutral medium in photodissociation regions (PDRs). We find that line deficits agree with those of local luminous IR galaxies. The CO spectral line energy distributions (SLEDs) reveal significant high- J CO excitation in both quasar hosts. Our CO SLED modeling of the quasar PJ231−20 shows that PDRs dominate the molecular mass and CO luminosities for J up  ≤ 7, while the J up  ≥ 10 CO emission is likely driven by X-ray dissociation regions produced by the active galactic nucleus (AGN) at the very center of the quasar host. The J up  > 10 lines are undetected in the other galaxies in our study. The H 2 O 3 21  − 3 12 line detection in the same quasar places this object on the L H 2 O  −  L TIR relation found for low- z sources, thus suggesting that this water vapor transition is predominantly excited by IR pumping. Models of the H 2 O SLED and of the H 2 O-to-OH 163 μm ratio point to PDR contributions with high volume and column density ( n H  ∼ 0.8 × 10 5 cm −3 , N H  = 10 24 cm −2 ) in an intense radiation field. Our analysis suggests a less highly excited medium in the companion galaxies. However, the current data do not allow us to definitively rule out an AGN in these sources, as suggested by previous studies of the same objects. This work demonstrates the power of multiline studies of FIR diagnostics in order to dissect the physical conditions in the first massive galaxies emerging from cosmic dawn. 
    more » « less
  2. Abstract

    We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope's SPT-SZ survey and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array. This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans 1.9 <z< 6.9 and covers rest-frame frequencies of 240–800 GHz. Combining this data with low-JCO observations from the Australia Telescope Compact Array, we detect multiple bright line features from12CO, [Ci], and H2O, as well as fainter molecular transitions from13CO, HCN, HCO+, HNC, CN, H2O+, and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high-redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright12CO emission lines, and contain either warmer and more excited dense gas tracers or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (LIR> 1012L), and will inform studies of chemical evolution before the peak epoch of star formation atz= 2–3.

    more » « less
  3. Abstract

    The [Cii] 158μm emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [Cii] emission lines and FIR continua of 31 luminous quasars atz> 6.5 using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array at sub-arcsec resolution. This survey more than doubles the number of quasars with [Cii] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [Cii] emission in 27 quasar hosts with a luminosity range ofL[CII]= (0.3–5.5) × 109Land detect the FIR continuum of 28 quasar hosts with a luminosity range ofLFIR= (0.5–13.0) × 1012L. BothL[CII]andLFIRare correlated (ρ≃ 0.4) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of ∼5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [Cii] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yieldMdyn= (0.1–7.5) × 1011M. Considering our findings alongside those of literature studies, we found that the ratio betweenMBHandMdynis about 10 times higher than that of localMBHMdynrelation on average but with substantial scatter (the ratio difference ranging from ∼0.6 to 60) and large uncertainties.

    more » « less
  4. Exploiting the sensitivity of the IRAM NOrthern Extended Millimeter Array (NOEMA) and its ability to process large instantaneous bandwidths, we have studied the morphology and other properties of the molecular gas and dust in the star forming galaxy, H-ATLAS J131611.5+281219 (HerBS-89a), at z = 2.95. High angular resolution (0 . ″3) images reveal a partial 1 . ″0 diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12 CO(9−8) and H 2 O(2 02  − 1 11 ). Together with lower angular resolution (0 . ″6) images, we report the detection of a series of molecular lines including the three fundamental transitions of the molecular ion OH + , namely (1 1  − 0 1 ), (1 2  − 0 1 ), and (1 0  − 0 1 ), seen in absorption; the molecular ion CH + (1 − 0) seen in absorption, and tentatively in emission; two transitions of amidogen (NH 2 ), namely (2 02  − 1 11 ) and (2 20  − 2 11 ) seen in emission; and HCN(11 − 10) and/or NH(1 2  − 0 1 ) seen in absorption. The NOEMA data are complemented with Very Large Array data tracing the 12 CO(1 − 0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present Hubble Space Telescope imaging that reveals the foreground lensing galaxy in the near-infrared (1.15  μ m). Together with photometric data from the Gran Telescopio Canarias, we derive a photometric redshift of z phot = 0.9 −0.5 +0.3 for the foreground lensing galaxy. Modeling the lensing of HerBS-89a, we reconstruct the dust continuum (magnified by a factor μ  ≃ 5.0) and molecular emission lines (magnified by μ  ∼ 4 − 5) in the source plane, which probe scales of ∼0 . ″1 (or 800 pc). The 12 CO(9 − 8) and H 2 O(2 02  − 1 11 ) emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter, which reveals two compact rotating components with sizes of ≈1 kpc that are likely merging, more naturally accounts for the broad line widths observed in HerBS-89a. In the core of HerBS-89a, very dense gas with n H 2  ∼ 10 7 − 9 cm −3 is revealed by the NH 2 emission lines and the possible HCN(11 − 10) absorption line. HerBS-89a is a powerful star forming galaxy with a molecular gas mass of M mol  = (2.1 ± 0.4) × 10 11   M ⊙ , an infrared luminosity of L IR  = (4.6 ± 0.4) × 10 12   L ⊙ , and a dust mass of M dust  = (2.6 ± 0.2) × 10 9   M ⊙ , yielding a dust-to-gas ratio δ GDR  ≈ 80. We derive a star formation rate SFR = 614 ± 59  M ⊙ yr −1 and a depletion timescale τ depl  = (3.4 ± 1.0) × 10 8 years. The OH + and CH + absorption lines, which trace low (∼100 cm −3 ) density molecular gas, all have their main velocity component red-shifted by Δ V  ∼ 100 km s −1 relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow toward the center of a galaxy, indicating that HerBS-89a is accreting gas from its surroundings. 
    more » « less
  5. Abstract

    We present a search for extremely red, dust-obscured,z> 7 galaxies with JWST/NIRCam+MIRI imaging over the first 20 arcmin2of publicly available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W−F444W (∼2.5 mag) and detection in MIRI/F770W (∼25 mag), we identify two galaxies, COS-z8M1 and CEERS-z7M1, that have best-fit photometric redshifts ofz=8.40.4+0.3and7.60.1+0.1, respectively. We perform spectral energy distribution fitting with a variety of codes (includingbagpipes,prospector,beagle, andcigale) and find a >95% probability that these indeed lie atz> 7. Both sources are compact (Reff≲ 200 pc) and highly obscured (AV∼ 1.5–2.5) and, at our best-fit redshift estimates, likely have strong [Oiii]+Hβemission contributing to their 4.4μm photometry. We estimate stellar masses of ∼1010Mfor both sources; by virtue of detection in MIRI at 7.7μm, these measurements are robust to the inclusion of bright emission lines, for example, from an active galactic nucleus. We identify a marginal (2.9σ) Atacama Large Millimeter/submillimeter Array detection at 2 mm within 0.″5 of COS-z8M1, which, if real, would suggest a remarkably high IR luminosity of ∼1012L. These two galaxies, if confirmed atz∼ 8, would be extreme in their stellar and dust masses and may be representative of a substantial population of highly dust-obscured galaxies at cosmic dawn.

    more » « less