skip to main content


Title: An Empirical Representation of a Physical Model for the ISM [C ii], CO, and [C i] Emission at Redshift 1 ≤ z ≤ 9
Abstract Submillimeter emission lines produced by the interstellar medium (ISM) are strong tracers of star formation and are some of the main targets of line intensity mapping (LIM) surveys. In this work we present an empirical multiline emission model that simultaneously covers the mean, scatter, and correlations of [C ii ], CO J = 1–0 to J = 5–4, and [C i ] lines in the redshift range 1 ≤ z ≤ 9. We assume that the galaxy ISM line emission luminosity versus halo mass relations can be described by double power laws with redshift-dependent lognormal scatter. The model parameters are then derived by fitting to the state-of-the-art semianalytic simulation results that have successfully reproduced multiple submillimeter line observations at 0 ≤ z ≲ 6. We cross-check the line emission statistics predicted by the semianalytic simulation and our empirical model, finding that at z ≥ 1 our model reproduces the simulated line intensities with fractional error less than about 10%. The fractional difference is less than 25% for the power spectra. Grounded on physically motivated and self-consistent galaxy simulations, this computationally efficient model will be helpful in forecasting ISM emission-line statistics for upcoming LIM surveys.  more » « less
Award ID(s):
2108411
NSF-PAR ID:
10330066
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
929
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
140
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Line intensity mapping (LIM) is emerging as a powerful technique to map the cosmic large-scale structure and to probe cosmology over a wide range of redshifts and spatial scales. We perform Fisher forecasts to determine the optimal design of wide-field ground-based millimeter-wavelength LIM surveys for constraining properties of neutrinos and light relics. We consider measuring the auto-power spectra of several CO rotational lines (from J = 2–1 to J = 6–5) and the [C ii ] fine-structure line in the redshift range of 0.25 < z < 12. We study the constraints with and without interloper lines as a source of noise in our analysis, and for several one-parameter and multiparameter extensions of ΛCDM. We show that LIM surveys deployable this decade, in combination with existing cosmic microwave background (CMB; primary) data, could achieve order-of-magnitude improvements over Planck constraints on N eff and M ν . Compared to next-generation CMB and galaxy surveys, a LIM experiment of this scale could achieve bounds that are a factor of ∼3 better than those forecasted for surveys such as EUCLID (galaxy clustering), and potentially exceed the constraining power of CMB-S4 by a factor of ∼1.5 and ∼3 for N eff and M ν , respectively. We show that the forecasted constraints are not substantially affected when enlarging the parameter space, and additionally demonstrate that such a survey could also be used to measure ΛCDM parameters and the dark energy equation of state exquisitely well. 
    more » « less
  2. 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
  3. Aims . We present and study spatially resolved imaging obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of multiple 12 CO( J  = 6 − 5, 8−7, and 9−8) and two H 2 O(2 02 −1 11 and 2 11 −2 02 ) emission lines and cold dust continuum toward the gravitationally lensed dusty star-forming galaxy SPT 0346-52 at z  = 5.656. Methods . Using a visibility-domain source-plane reconstruction we probe the structure and dynamics of the different components of the interstellar medium (ISM) in this galaxy down to scales of 1 kpc in the source plane. Results . Measurements of the intrinsic sizes of the different CO emission lines indicate that the higher J transitions trace more compact regions in the galaxy. Similarly, we find smaller dust continuum intrinsic sizes with decreasing wavelength, based on observations at rest frame 130, 300, and 450 μ m. The source shows significant velocity structure, and clear asymmetry where an elongated structure is observed in the source plane with significant variations in their reconstructed sizes. This could be attributed to a compact merger or turbulent disk rotation. The differences in velocity structure through the different line tracers, however, hint at the former scenario in agreement with previous [CII] line imaging results. Measurements of the CO line ratios and magnifications yield significant variations as a function of velocity, suggesting that modeling of the ISM using integrated values could be misinterpreted. Modeling of the ISM in SPT 0346-52 based on delensed fluxes indicates a highly dense and warm medium, qualitatively similar to that observed in high-redshift quasar hosts. 
    more » « less
  4. Abstract Line intensity mapping (LIM) is a rapidly emerging technique for constraining cosmology and galaxy formation using multi-frequency, low angular resolution maps.Many LIM applications crucially rely on cross-correlations of two line intensity maps, or of intensity maps with galaxy surveys or galaxy/CMB lensing.We present a consistent halo model to predict all these cross-correlations and enable joint analyses, in 3D redshift-space and for 2D projected maps.We extend the conditional luminosity function formalism to the multi-line case, to consistently account for correlated scatter between multiple galaxy line luminosities.This allows us to model the scale-dependent decorrelation between two line intensity maps,a key input for foreground rejection and for approaches that estimate auto-spectra from cross-spectra.This also enables LIM cross-correlations to reveal astrophysical properties of the interstellar medium inacessible with LIM auto-spectra.We expose the different sources of luminosity scatter or “line noise” in LIM, and clarify their effects on the 1-halo and galaxy shot noise terms.In particular, we show that the effective number density of halos can in some cases exceed that of galaxies, counterintuitively.Using observational and simulation input, we implement this halo model for the Hα, [Oiii], Lyman-α, CO and [Cii] lines.We encourage observers and simulators to measure galaxy luminosity correlation coefficients for pairs of lines whenever possible.Our code is publicly available at https://github.com/EmmanuelSchaan/HaloGen/tree/LIM .In a companion paper, we use this halo model formalism and codeto highlight the degeneracies between cosmology and astrophysics in LIM, and to compare the LIM observables to galaxy detection for a number of surveys. 
    more » « less
  5. Abstract We present an update to the framework called Simulator of Galaxy Millimeter/submillimeter Emission ( sígame ). sígame derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a postprocessing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parameterizing the probability distribution function (PDF) of the gas density in higher-resolution simulations run with the pseudo-Lagrangian, Voronoi mesh code arepo . The parameterized PDFs are used as a look-up table, and reach higher densities than in previous work. sígame v3 is tested on redshift z = 0 galaxies drawn from the simba cosmological simulation for eight FIR emission lines tracing vastly different phases of the interstellar medium. This version of sígame includes dust radiative transfer with S kirt and high-resolution photoionization models with C loudy , the latter sampled according to the density PDF of the arepo simulations to augment the densities in the cosmological simulation. The quartile distributions of the predicted line luminosities overlap with the observed range for nearby galaxies of similar star formation rate (SFR) for all but two emission lines: [O i ]63 and CO(3–2), which are overestimated by median factors of 1.3 and 1.0 dex, respectively, compared to the observed line–SFR relation of mixed-type galaxies. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on sub-grid scales (≲200 pc) and differences in sample selection. 
    more » « less