skip to main content

This content will become publicly available on April 1, 2023

Title: On Estimating the Cosmic Molecular Gas Density from CO Line Intensity Mapping Observations
Abstract The Millimeter-wave Intensity Mapping Experiment (mmIME) recently reported a detection of excess spatial fluctuations at a wavelength of 3 mm, which can be attributed to unresolved emission of several CO rotational transitions between z ∼ 1 and 5. We study the implications of these data for the high-redshift interstellar medium using a suite of state-of-the-art semianalytic simulations that have successfully reproduced many other submillimeter line observations across the relevant redshift range. We find that the semianalytic predictions are mildly in tension with the mmIME result, with a predicted CO power ∼3.5 σ below what was observed. We explore some simple modifications to the models that could resolve this tension. Increasing the molecular gas abundance at the relevant redshifts to ∼10 8 M ⊙ Mpc −3 , a value well above that obtained from directly imaged sources, would resolve the discrepancy, as would assuming a CO–H 2 conversion factor α CO of ∼1.5 M ⊙ K −1 (km s −1 ) −1 pc 2 , a value somewhat lower than is commonly assumed. We go on to demonstrate that these conclusions are quite sensitive to the detailed assumptions of our simulations, highlighting the need for more careful modeling efforts as more » more intensity mapping data become available. « less
Authors:
; ; ; ; ;
Award ID(s):
2108411
Publication Date:
NSF-PAR ID:
10330064
Journal Name:
The Astrophysical Journal
Volume:
929
Issue:
1
Page Range or eLocation-ID:
30
ISSN:
0004-637X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.
  2. Abstract We present the current state of models for the z ∼ 3 carbon monoxide (CO) line intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy–halo connection and previous CO (1–0) observations. The Pathfinder early science data spanning wavenumbers k = 0.051–0.62 Mpc −1 represent the first direct 3D constraint on the clustering component of the CO (1–0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude A clust ≲ 70 μ K 2 greatly improves on the indirect upper limit of 420 μ K 2 reported from the CO Power Spectrum Survey (COPSS) measurement at k ∼ 1 Mpc −1 . The COMAP limit excludes a subset of models from previous literature and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity–bias product, Tb 2 ≲ 50 μ K 2 , and the cosmic molecular gas density, ρ H2 < 2.5 × 10 8 Mmore »⊙ Mpc −3 (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the 5 yr Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise ratio of 9–17. Between then and now, we also expect to detect the CO–galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star formation and galaxy evolution history.« less
  3. 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 Mmore »ν , 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.« less
  4. Abstract

    We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semianalytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ∼0.1μas resolution and ∼1μJy sensitivity could access >106SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ≲2μas resolution and ≲0.5 mJy sensitivity, whereas doing the same for second-order photon rings requires ≲0.1μas resolution and ≲5μJy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as ∼5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), whereas a next-generation EHT observing at 345 GHz should have access to ∼3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experimentsmore »that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

    « less
  5. Abstract

    We present Symphony, a compilation of 262 cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from 1011–1015M. This compilation includes three existing simulation suites at the cluster and Milky Way–mass scales, and two new suites: 39 Large Magellanic Cloud-mass (1011M) and 49 strong-lens-analog (1013M) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of ≈3 × 10−7times the host virial mass and a Plummer-equivalent gravitational softening length of ≈9 × 10−4times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are ≈50% higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts’ underlying dark matter distributions. We compare our results to the semianalytic modelGalacticus, which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We useUniverseMachineto model halo and subhalo star formationmore »histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available athttp://web.stanford.edu/group/gfc/symphony.

    « less