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
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
- PAR ID:
- 10112706
- 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
-
-
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
-
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
-
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
-
ABSTRACT We present 10 main-sequence ALPINE galaxies (log (M/M⊙) = 9.2−11.1 and $${\rm SFR}=23-190\, {\rm M_{\odot }\, yr^{-1}}$$) at z ∼ 4.5 with optical [O ii] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multiwavelength sample at these redshifts, combining various measurements in the ultraviolet, optical, and far-infrared including [C ii]158 $$\mu$$m line emission and dust continuum from ALMA and H α emission from Spitzer photometry. For the first time, this unique sample allows us to analyse the relation between [O ii] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [O ii]/[C ii] and [O ii]/H α luminosity ratios at z ∼ 4.5. The [O ii]−SFR relation at z ∼ 4.5 cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $$50{{\ \rm per\ cent}}$$ solar. To explain the measured dust-corrected luminosity ratios of $$\log (L_{\rm [OII]}/L_{\rm [CII]}) \sim 0.98^{+0.21}_{-0.22}$$ and $$\log (L_{\rm [OII]}/L_{\rm H\alpha }) \sim -0.22^{+0.13}_{-0.15}$$ for our sample, ionization parameters log (U) < −2 and electron densities $$\log (\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$$ are required. The former is consistent with galaxies at z ∼ 2−3, however lower than at z > 6. The latter may be slightly higher than expected given the galaxies’ specific SFR. The analysis of this pilot sample suggests that typical log (M/M⊙) > 9 galaxies at z ∼ 4.5 to have broadly similar ISM properties as their descendants at z ∼ 2 and suggest a strong evolution of ISM properties since the epoch of reionization at z > 6.more » « less
An official website of the United States government

