skip to main content


Title: On the dust temperatures of high-redshift galaxies
Abstract Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimetre broad-band flux to total infrared luminosity (LIR), and hence star formation rate, in high-redshift galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g. redshift and star formation rate. In this paper, we analyse the dust temperatures of massive ($M_{\rm star} \gt 10^{10}\, \mathrm{M}_{\odot }$) $z$ = 2–6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (fire) project. At $z$ ∼ 2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature, which is required to accurately estimate the total dust mass, does not strongly evolve with redshift over $z$ = 2–6 at fixed IR luminosity but is tightly correlated with LIR at fixed $z$. We also analyse an ‘equivalent’ dust temperature for converting (sub)millimetre flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature (‘warmer dust’) do not necessarily have higher mass-weighted temperatures. A ‘two-phase’ picture for interstellar dust can explain the different scaling relations of the various dust temperatures.  more » « less
Award ID(s):
1715847 1713353 1715216 1715101
NSF-PAR ID:
10112706
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
489
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
1397 to 1422
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT

    Recent years have seen growing interest in post-processing cosmological simulations with radiative transfer codes to predict observable fluxes for simulated galaxies. However, this can be slow, and requires a number of assumptions in cases where simulations do not resolve the interstellar medium (ISM). Zoom-in simulations better resolve the detailed structure of the ISM and the geometry of stars and gas; however, statistics are limited due to the computational cost of simulating even a single halo. In this paper, we make use of a set of high-resolution, cosmological zoom-in simulations of massive ($M_{\star }\gtrsim 10^{10.5}\, \rm {M_{\odot }}$ at z = 2), star-forming galaxies from the FIRE suite. We run the skirt radiative transfer code on hundreds of snapshots in the redshift range 1.5 < z < 5 and calibrate a power-law scaling relation between dust mass, star formation rate, and $870\, \mu \rm {m}$ flux density. The derived scaling relation shows encouraging consistency with observational results from the sub-millimetre-selected AS2UDS sample. We extend this to other wavelengths, deriving scaling relations between dust mass, stellar mass, star formation rate, and redshift and sub-millimetre flux density at observed-frame wavelengths between $\sim \! 340$ and $\sim \! 870\, \mu \rm {m}$. We then apply the scaling relations to galaxies drawn from EAGLE, a large box cosmological simulation. We show that the scaling relations predict EAGLE sub-millimetre number counts that agree well with previous results that were derived using far more computationally expensive radiative transfer techniques. Our scaling relations can be applied to other simulations and semi-analytical or semi-empirical models to generate robust and fast predictions for sub-millimetre number counts.

     
    more » « less
  3. ABSTRACT We report the detection of the far-infrared (FIR) fine-structure line of singly ionized nitrogen, [N ii] 205 $\mu$m , within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ‘The Red Radio Ring’; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (Jup = 1, 5, 8), and the FIR spectral energy distribution (SED), to explore the multiphase interstellar medium (ISM) properties of the RRR. All line profiles suggest that the H ii regions, traced by [N ii] 205 $\mu$m , and the (diffuse and dense) molecular gas, traced by CO, are cospatial when averaged over kpc-sized regions. Using its mid-IR-to-millimetre (mm) SED, we derive a non-negligible dust attenuation of the [N ii] 205 $\mu$m line emission. Assuming a uniform dust screen approximation results a mean molecular gas column density >1024 cm−2, with a molecular gas-to-dust mass ratio of 100. It is clear that dust attenuation corrections should be accounted for when studying FIR fine-structure lines in such systems. The attenuation corrected ratio of $L_{\rm N\,{\small II}205} / L_{\rm IR(8\!-\!1000\, \mu m)} = 2.7 \times 10^{-4}$ is consistent with the dispersion of local and z > 4 SFGs. We find that the lower limit, [N ii] 205 $\mu$m -based star formation rate (SFR) is less than the IR-derived SFR by a factor of 4. Finally, the dust SED, CO line SED, and $L_{\rm N\,{\small II}205}$ line-to-IR luminosity ratio of the RRR is consistent with a starburst-powered ISM. 
    more » « less
  4. 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
  5. ABSTRACT

    We post-process galaxies in the IllustrisTNG simulations with skirt radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at z ≥ 4. The rest-frame K- and z-band galaxy luminosity functions from TNG are overall consistent with observations, despite ${\sim}0.5\, \mathrm{dex}$ underprediction at z = 4 for MK ≲ −25 and Mz ≲ −24. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. Based on theoretical estimations, we show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\, {\rm arcmin}^{2}$ at z = 6 (z = 8). As opposed to the consistency in the UV, optical, and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of z ≥ 6 galaxies are overestimated by about $20\, {\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.

     
    more » « less