Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full
advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$alpha$ emitters). We find that obtaining $sigma_z/(1+z)lesssim0.003$ accuracy in redshifts and $gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements.
Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opp
ortunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.
Line-intensity mapping surveys probe large-scale structure through spatial variations in molecular line emission from a population of unresolved cosmological sources. Future such surveys of carbon monoxide line emission, specifically the CO(1-0) line
, face potential contamination from a disjoint population of sources emitting in a hydrogen cyanide emission line, HCN(1-0). This paper explores the potential range of the strength of HCN emission and its effect on the CO auto power spectrum, using simulations with an empirical model of the CO/HCN--halo connection. We find that effects on the observed CO power spectrum depend on modeling assumptions but are very small for our fiducial model based on our understanding of the galaxy--halo connection, with the bias in overall CO detection significance due to HCN expected to be less than 1%.
Intensity mapping provides a unique means to probe the epoch of reionization (EoR), when the neutral intergalactic medium was ionized by the energetic photons emitted from the first galaxies. The [CII] 158$mu$m fine-structure line is typically one of
the brightest emission lines of star-forming galaxies and thus a promising tracer of the global EoR star-formation activity. However, [CII] intensity maps at $6 lesssim z lesssim 8$ are contaminated by interloping CO rotational line emission ($3 leq J_{rm upp} leq 6$) from lower-redshift galaxies. Here we present a strategy to remove the foreground contamination in upcoming [CII] intensity mapping experiments, guided by a model of CO emission from foreground galaxies. The model is based on empirical measurements of the mean and scatter of the total infrared luminosities of galaxies at $z < 3$ and with stellar masses $M_{*} > 10^{8},rm M_{rm odot}$ selected in $K$-band from the COSMOS/UltraVISTA survey, which can be converted to CO line strengths. For a mock field of the Tomographic Ionized-carbon Mapping Experiment (TIME), we find that masking out the voxels (spectral-spatial elements) containing foreground galaxies identified using an optimized CO flux threshold results in a $z$-dependent criterion $m^{rm AB}_{rm K} lesssim 22$ (or $M_{*} gtrsim 10^{9} ,rm M_{rm odot}$) at $z < 1$ and makes a [CII]/CO$_{rm tot}$ power ratio of $gtrsim 10$ at $k=0.1$ $h$/Mpc achievable, at the cost of a moderate $lesssim 8%$ loss of total survey volume.
We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic (sub)millimeter background at: 250, 350, and 500 um (1200, 860, and 600 GHz) from observations made with the Balloon-borne Large Aperture Submillimeter T
elescope, BLAST; and at 1380 and 2030 um (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. The overlapping observations cover 8.6 deg^2 in an area relatively free of Galactic dust near the south ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zeldovich (SZ) effect from galaxy clusters, and to emission by radio and dusty star-forming galaxies (DSFGs), while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the BLAST analysis of clustering with an independent pipeline, and also detect correlations between the ACT and BLAST maps at over 25sigma significance, which we interpret as a detection of the DSFGs in the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a clustered signal at >4sigma, and using a model for the DSFG evolution and number counts, we successfully fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on scale. Finally, the data are compared to, and generally agree with, phenomenological models for the DSFG population. This study represents a first of its kind, and demonstrates the constraining power of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with more data will impose tight constraints on future models.
We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (Mast approx 2 times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS)
catalog on maps: at 24 {mu}m (Spitzer/MIPS); 70, 100, and 160{mu}m (Herschel/PACS); 250, 350, 500{mu}m (BLAST); and 870{mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.
We detect correlations in the cosmic far-infrared background due to the clustering of star-forming galaxies in observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST, at 250, 350, and 500 microns. We perform jackknife
and other tests to confirm the reality of the signal. The measured correlations are well fit by a power law over scales of 5-25 arcminutes, with Delta I/I = 15.1 +/- 1.7%. We adopt a specific model for submillimeter sources in which the contribution to clustering comes from sources in the redshift ranges 1.3 <= z <= 2.2, 1.5 <= z <= 2.7, and 1.7 <= z <= 3.2, at 250, 350, and 500 microns, respectively. With these distributions, our measurement of the power spectrum, P(k_theta), corresponds to linear bias parameters, b = 3.8 +/- 0.6, 3.9 +/- 0.6 and 4.4 +/- 0.7, respectively. We further interpret the results in terms of the halo model, and find that at the smaller scales, the simplest halo model fails to fit our results. One way to improve the fit is to increase the radius at which dark matter halos are artificially truncated in the model, which is equivalent to having some star-forming galaxies at z >= 1 located in the outskirts of groups and clusters. In the context of this model we find a minimum halo mass required to host a galaxy is log (M_min / M_sun) = 11.5 (+0.4/-0.1), and we derive effective biases $b_eff = 2.2 +/- 0.2, 2.4 +/- 0.2, and 2.6 +/- 0.2, and effective masses log (M_eff / M_sun) = 12.9 +/- 0.3, 12.8 +/- 0.2, and 12.7 +/- 0.2, at 250, 350, and 500 microns, corresponding to spatial correlation lengths of r_0 = 4.9, 5.0, and 5.2 +/- 0.7 h^-1 Mpc, respectively. Finally, we discuss implications for clustering measurement strategies with Herschel and Planck.
