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
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 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
  2. 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
  3. 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
  4. Aims.We have implemented a novel method to create simulated [CII] emission line intensity mapping (LIM) data cubes using COSMOS 2020 galaxy catalogue data. It allows us to provide solid lower limits for previous simulation-based model predictions and the expected signal strength of upcoming surveys.

    Methods.We applied [CII]158 μm luminosity models to COSMOS 2020 to create LIM cubes covering a 1.2 × 1.2deg2sky area. These models were derived using galaxy bulk property data from the ALPINE-ALMA survey over the redshift range of 4.4 < z < 5.9, while additional models were taken from the literature. The LIM cubes cover 3.42 < z < 3.87, 4.14 < z < 4.76, 5.34 < z < 6.31, and 6.75 < z < 8.27, matched to planned observations from the EoR-Spec module of the Prime-Cam instrument in the Fred Young Submillimeter Telescope (FYST). We also created predictions including additional galaxies below current detection limits by ‘extrapolating’ from the faint end of the COSMOS 2020 luminosity function, comparing these to predictions from the literature. In addition, we computed the signal-to-noise (S/N) ratios for the power spectra, using parameters from the planned FYST survey with predicted instrumental noise levels.

    Results.We find lower limits for the expected power spectrum using the likely incomplete empirical data: when normalised by 2π2, the amplitudes atk = 1 Mpc−1are 3.06 × 107, 1.43 × 107, 9.80 × 105, 2.77 × 105 (Jy sr−1)2for the aforementioned redshift ranges. For the extrapolated sample, the power spectra are consistent with prior predictions, indicating that extrapolation is a viable method for creating mock LIM cubes. In this case, we expect a result of S/N> 1 when using FYST parameters. However, our high-redshift results remain inconclusive because of the poor completeness of COSMOS 2020 atz > 6.3. These predictions will be improved on the basis of future JWST data.

     
    more » « less
  5. Abstract

    Detecting the line-intensity mapping (LIM) signal from the galaxies of the epoch of reionization is an emerging tool to constrain their role in reionization. Ongoing and upcoming experiments target the signal fluctuations across the sky to reveal statistical and astrophysical properties of these galaxies via signal statistics, e.g. the power spectrum. Here, we revisit the [C ii]$_{158 \, \mu \text{m}}$ LIM power spectrum under non-uniform line–luminosity scatter, which has a halo-mass variation of statistical properties. Line–luminosity scatter from a cosmological hydrodynamic and radiative transfer simulation of galaxies at $z$ = 6 is considered in this study. We test the robustness of different model frameworks that interpret the impact of the line-luminosity scatter on the signal statistics. We use a simple power-law model to fit the scatter and demonstrate that the mean luminosity–halo mass correlation fit cannot preserve the mean intensity of the LIM signal (hence the clustering power spectrum) under non-uniform scatter. In our case, the mean intensity changes by ∼48 per cent compared to the mean correlation fit in contrast to the general case with semi-analytical scatter. However, we find that the prediction for the mean intensity from the most-probable fit can be modelled robustly, considering the generalized and more realistic non-uniform scatter. We also explore the possibility of diminishing luminosity bias under non-uniform scatter, affecting the clustering power spectrum, although this phenomenon might not be statistically significant. Therefore, we should adopt appropriate approaches that can consistently interpret the LIM power spectrum from observations.

     
    more » « less