skip to main content


Title: Dust Temperature Uncertainties Hamper the Inference of Dust and Molecular Gas Masses from the Dust Continuum Emission of Quiescent High-redshift Galaxies
Abstract

Single flux density measurements at observed-frame submillimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive halos from the Feedback in Realistic Environments project. Our starting point is four star-forming central galaxies at seven redshifts betweenz= 1.5 andz= 4.5. We generate modified quiescent galaxies that have been quenched for 100 Myr, 500 Myr, or 1 Gyr prior to each of the studied redshifts by reassigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature,Tdust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumedTdust∼ 25 K, we show that the ratio between dust mass and 1.3 mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non-star-forming galaxies. This can lead to an underestimation of dust mass (and, when submillimeter flux density is used as a proxy for molecular gas content and gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3 mm flux density probes rest-frame wavelengths far from the Rayleigh–Jeans regime, and hence depends superlinearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses.

 
more » « less
Award ID(s):
2009687
NSF-PAR ID:
10379708
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
939
Issue:
2
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L27
Size(s):
["Article No. L27"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Context. Recent years have seen building evidence that planet formation starts early, in the first ~0.5 Myr. Studying the dust masses available in young disks enables us to understand the origin of planetary systems given that mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. Aims. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. Methods. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.1–1.3 mm) continuum observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka -band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities measured in the image plane. Results. We find a strong linear correlation between the ALMA and VLA fluxes, demonstrating that emission at both wavelengths is dominated by dust emission. For a subsample of optically thin sources, we find a median spectral index of 2.5 from which we derive the dust opacity index β = 0.5, suggesting significant dust growth. Comparison with ALMA surveys of Orion shows that the Class I dust disk mass distribution between the two regions is similar, but that the Class 0 disks are more massive in Perseus than those in Orion. Using the DIANA opacity model including large grains, with a dust opacity value of κ 9 mm = 0.28 cm 2 g −1 , the median dust masses of the embedded disks in Perseus are 158 M ⊕ for Class 0 and 52 M ⊕ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M ⊕ and 12 M ⊕ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value κ 1.3 mm = 2.3 cm 2 g −1 . The dust masses of young Class 0 and I disks are larger by at least a factor of ten and three, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. Conclusions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of ~15%. A higher efficiency of ~30% is necessary if the planet formation is set to start in Class I disks. 
    more » « less
  2. ABSTRACT The infrared (IR) spectral energy distributions (SEDs) of main-sequence galaxies in the early Universe (z > 4) is currently unconstrained as IR continuum observations are time-consuming and not feasible for large samples. We present Atacama Large Millimetre Array Band 8 observations of four main-sequence galaxies at z ∼ 5.5 to study their IR SED shape in detail. Our continuum data (rest-frame 110 $\rm \mu m$, close to the peak of IR emission) allows us to constrain luminosity-weighted dust temperatures and total IR luminosities. With data at longer wavelengths, we measure for the first time the emissivity index at these redshifts to provide more robust estimates of molecular gas masses based on dust continuum. The Band 8 observations of three out of four galaxies can only be reconciled with optically thin emission redward of rest-frame $100\, {\rm \mu m}$. The derived dust peak temperatures at z ∼ 5.5 ($30\!-\!43\, {\rm K}$) are elevated compared to average local galaxies, however, $\sim 10\, {\rm K}$ below what would be predicted from an extrapolation of the trend at z < 4. This behaviour can be explained by decreasing dust abundance (or density) towards high redshifts, which would cause the IR SED at the peak to be more optically thin, making hot dust more visible to the external observer. From the $850{\hbox{-}}{\rm \mu m}$ dust continuum, we derive molecular gas masses between 1010 and $10^{11}\, {\rm M_{\odot }}$ and gas fractions (gas over total mass) of $30\!-\!80{{\ \rm per\ cent}}$ (gas depletion times of $100\!-\!220\, {\rm Myr}$). All in all, our results provide a first measured benchmark SED to interpret future millimetre observations of normal, main-sequence galaxies in the early Universe. 
    more » « less
  3. Abstract We present results from Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm continuum observations of a sample of 27 star-forming galaxies at z = 2.1–2.5 from the MOSFIRE Deep Evolution Field survey with metallicity and star formation rate measurements from optical emission lines. Using stacks of Spitzer, Herschel, and ALMA photometry (rest frame ∼8–400 μ m), we examine the infrared (IR) spectral energy distributions (SED) of z ∼ 2.3 subsolar-metallicity (∼0.5 Z ⊙ ) luminous infrared galaxies (LIRGs). We find that the data agree well with an average template of higher-luminosity local low-metallicity dwarf galaxies (reduced χ 2 = 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ultraluminous IR galaxies, even in the most favorable case (with reduced χ 2 = 2.8), the templates are rejected at >98% confidence. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar-metallicity galaxies may result from different grain properties or a harder/more intense ionizing radiation field that increases the dust temperature. The obscured star formation rate (SFR) indicated by the far-IR emission of the subsolar-metallicity galaxies is only ∼60% of the total SFR, considerably lower than that of the local LIRGs with ∼96%–97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data overestimate the Rayleigh–Jeans tail measurements by a factor of 2–20. These templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4 dex. At a given stellar mass or metallicity, dust masses at z ∼ 2.3 are an order of magnitude higher than z ∼ 0. Given the predicted molecular gas fractions, the observed z ∼ 2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies with similar metallicities. 
    more » « less
  4. ABSTRACT We present the first detailed study of the spatially resolved dust continuum emission of simulated galaxies at 1 < z < 5. We run the radiative transfer code skirt on a sample of submillimetre-bright galaxies drawn from the Feedback In Realistic Environments (FIRE) project. These simulated galaxies reach Milky Way masses by z = 2. Our modelling provides predictions for the full rest-frame far-ultraviolet-to-far-infrared spectral energy distributions of these simulated galaxies, as well as 25-pc resolution maps of their emission across the wavelength spectrum. The derived morphologies are notably different in different wavebands, with the same galaxy often appearing clumpy and extended in the far-ultraviolet yet an ordered spiral at far-infrared wavelengths. The observed-frame 870-$\mu$m half-light radii of our FIRE-2 galaxies are ${\sim} 0.5\rm {-}4\, \rm {kpc}$, consistent with existing ALMA observations of galaxies with similarly high redshifts and stellar masses. In both simulated and observed galaxies, the dust continuum emission is generally more compact than the cold gas and the dust mass, but more extended than the stellar component. The most extreme cases of compact dust emission seem to be driven by particularly compact recent star formation, which generates steep dust temperature gradients. Our results confirm that the spatial extent of the dust continuum emission is sensitive to both the dust mass and star formation rate distributions. 
    more » « less
  5. Abstract

    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) atz≈ 1.64–2.51. Our observations yielded one new detection of CO(3–2) emission from a galaxy atz= 2.4604 using NOEMA, associated with thez= 2.4628 DLA toward QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of5624+38% and119+26%, 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 × 1010M. This indicates that the highest-metallicity DLAs atz≈ 2 are associated with the most massive galaxies. The newly identifiedz≈ 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) × (r31/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.3Myr−1on 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.

     
    more » « less