Joint Power Spectrum and Voxel Intensity Distribution Forecast on the CO Luminosity Function with COMAP

H. T. Ihle, D. Chung, G. Stein, M. Alvarez, J. R. Bond, P. C. Breysse, K. A. Cleary, H. K. Eriksen, M. K. Foss, J. O. Gundersen, S. Harper, N. Murray, H. Padmanabhan, M. P. Viero, I. K. Wehus

Research output: Contribution to journalArticlepeer-review

13 Scopus citations


We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID) and apply this to simulations of CO Mapping Array Pathfinder (COMAP), a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the L CO (M halo ) model of Li et al. We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at large luminosities, while the PS and the VID are both competitive at small and intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between , and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58% and 30% for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available.

Original languageEnglish (US)
Article number75
JournalAstrophysical Journal
Issue number1
StatePublished - Jan 20 2019


  • diffuse radiation
  • galaxies: high-redshift
  • large-scale structure of universe

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science


Dive into the research topics of 'Joint Power Spectrum and Voxel Intensity Distribution Forecast on the CO Luminosity Function with COMAP'. Together they form a unique fingerprint.

Cite this