skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

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

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Energy balance SED modelling can be effective at high redshifts regardless of UV-FIR offsets
ABSTRACT Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially ‘decoupled’). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies’ true properties. In this work, we test how well the SED fitting code magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass, and luminosity by fitting 6706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4–500 μm, and spanning 1 < z < 8) with fits based on either the starlight ($$\lambda _\mathrm{eff} \le 2.2\, \mu$$m) or dust ($$\ge 100\, \mu$$m) alone, we highlight the power of considering the full range of multiwavelength data alongside an energy balance criterion. Overall, we obtain acceptable fits for 83 per cent of the synthetic SEDs, though the success rate falls rapidly beyond z ≈ 4, in part due to the sparser sampling of the priors at earlier times since SFHs must be physically plausible (i.e. shorter than the age of the universe). We use the ground truth from the simulations to show that when the quality of fit is acceptable, the fidelity of magphys estimates is independent of the degree of UV/FIR offset, with performance very similar to that previously reported for local galaxies.  more » « less
Award ID(s):
2009687
PAR ID:
10458994
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
525
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
1535 to 1552
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    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
  2. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    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
  3. ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    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
  4. ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution thesan reionization simulations. A combination of subresolution photoionization modelling for H ii regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies (z ≳ 5.5). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity–star formation rate relation (L–SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-z ones because the ISM of reionization era galaxies is generally less metal enriched than in their low-redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross-correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign (ktransition) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that ktransition probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email