Search

We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail., Comment: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to ApJ

We present the current state of models for the $z\sim3$ 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_{\rm clust}\lesssim70\,\mu$K$^2$ greatly improves on the indirect upper limit of $420\,\mu$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$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, $\langle{Tb}\rangle^2\lesssim50\,\mu$K$^2$, and the cosmic molecular gas density, $\rho_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise 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., Comment: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $\mu$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR., Comment: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Line Intensity Mapping (LIM) offers a novel avenue to observe and characterize our universe. LIM data of CO spectral lines are becoming available, such as those obtained by the CO Mapping Array Project (COMAP). COMAP data can be used to probe the molecular gas content of the universe from the last stages of the Epoch of Reionization (EoR) ($z < 8.0$) to $z \sim 2.5$. In this work, we examine the prospects for deriving voxel-level statistical constraints on high-redshift galaxies from COMAP data by considering the additional information available from observations of LAEs galaxies using the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope (HET). We post-process the IllustrisTNG300 galaxy-formation simulation with a set of prescriptions to consistently determine CO and Ly$\alpha$ line luminosities. The different line prescriptions span the uncertainty in the CO line luminosity according to current observations by the VLA high-z CO surveys and set the Ly$\alpha$ emission to be compatible with observational LAE luminosity functions. We produce mock observations for the two surveys over a $(300\, {\rm Mpc})^3$ volume. These are then used to formulate and test methodologies for data analysis and to predict COMAP constraints on CO emission. We use combinations of masking, stacking, voxel intensity distribution (VID), and other statistics. We find that in combination with VIRUS/HET, a voxel-level analysis of the COMAP Pathfinder survey can detect and characterize the CO signal from $z\sim3$ and improve current constraints on the $z\sim6$ signal, identify individual voxels with bright CO(1-0) emission at $z\sim3$ and probe the redshift evolution of the CO emission. This study illustrates the potential of synergies between LIM and galaxy surveys both to improve the significance of a detection and to aid the interpretation of noisy LIM data., Comment: 15 pages, 9 figures, Submitted to Astronomy & Astrophysics

Line-intensity mapping observations will find fluctuations of integrated line emission are attenuated by varying degrees at small scales due to the width of the line emission profiles. This attenuation may significantly impact estimates of astrophysical or cosmological quantities derived from measurements. We consider a theoretical treatment of the effect of line broadening on both the clustering and shot-noise components of the power spectrum of a generic line-intensity power spectrum using a halo model. We then consider possible simplifications to allow easier application in analysis, particularly in the context of inferences that require numerous, repeated, fast computations of model line-intensity signals across a large parameter space. For the CO Mapping Array Project (COMAP) and the CO(1-0) line-intensity field at $z\sim3$ serving as our primary case study, we expect a $\sim10\%$ attenuation of the spherically averaged power spectrum on average at relevant scales of $k\approx0.2$-$0.3$ Mpc$^{-1}$, compared to $\sim25\%$ for the interferometric Millimetre-wave Intensity Mapping Experiment (mmIME) targeting shot noise from CO lines at $z\sim1$-$5$ at scales of $k\gtrsim1$ Mpc$^{-1}$. We also consider the nature and amplitude of errors introduced by simplified treatments of line broadening, and find that while an approximation using a single effective velocity scale is sufficient for spherically-averaged power spectra, a more careful treatment is necessary when considering other statistics such as higher multipoles of the anisotropic power spectrum or the voxel intensity distribution., Comment: 24 pages + appendix and bibliography (33 pages total), 16 figures, 2 tables; accepted for publication in ApJ

The intergalactic medium (IGM) plays an important role in the formation and evolution of galaxies. Recent developments in upcoming radio telescopes are starting to open up the possibility of making a first direct detection of the 21 cm signal of neutral hydrogen (HI) from the warm gas of the IGM in large-scale filaments. The cosmological hydrodynamical EAGLE simulation is used to estimate the typical IGM filament signal. Assuming the same average signal for all filaments, a prediction is made for the detectability of such a signal with the upcoming mid-frequency array of the Square Kilometer Array (SKA1-mid) or the future upgrade to SKA2. The signal-to-noise (S/N) then only depends on the size and orientation of each filament. With filament spines inferred from existing galaxy surveys as a proxy for typical real filaments, we find hundreds of filaments in the region of the sky accessible to the SKA that can be detected. Once the various phases of the SKA telescope become operational, their own surveys will be able to find the galaxies required to infer the position of even more filaments within the survey area. We find that in 120 h, SKA1-mid/SKA2 will detect HI emission from the strongest filaments in the field with a S/N of the order of 10 to $\sim$150 for the most pessimistic model considered here. Some of the brighter filaments can be detected with an integration time of a few minutes with SKA1-mid and a few seconds with SKA2. Therefore, SKA2 will be capable of not only detecting but also mapping a large part of the IGM in these filaments., Comment: Accepted for publication in MNRAS, 12 pages, 6 figures, 3 tables

This paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. We argue that the proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution over a reasonable fraction of the volume of the observable Universe. In addition, the surveyed area needs to be very large to beat cosmic variance and to probe the largest scales where its easier to separate the astrophysical and cosmological contributions to the observed signal. Only, a space-borne mission can properly meet these requirements., Comment: 25 pages, 6 figures, ESA Voyage-2050 White Paper

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter., Comment: Astro2020 white paper (submitted to the NASA 2020 astronomy and astrophysics decadal survey)

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations., Comment: 62 pages

The near-infrared background between 0.5 $\mu$m to 2 $\mu$m contains a wealth of information related to radiative processes in the universe. Infrared background anisotropies encode the redshift-weighted total emission over cosmic history, including any spatially diffuse and extended contributions. The anisotropy power spectrum is dominated by undetected galaxies at small angular scales and diffuse background of Galactic emission at large angular scales. In addition to these known sources, the infrared background also arises from intra-halo light (IHL) at $z < 3$ associated with tidally-stripped stars during galaxy mergers. Moreover, it contains information on the very first galaxies from the epoch of reionization (EoR). The EoR signal has a spectral energy distribution (SED) that goes to zero near optical wavelengths due to Lyman absorption, while other signals have spectra that vary smoothly with frequency. Due to differences in SEDs and spatial clustering, these components may be separated in a multi-wavelength-fluctuation experiment. To study the extent to which EoR fluctuations can be separated in the presence of IHL, extra-galactic and Galactic foregrounds, we develop a maximum likelihood technique that incorporates a full covariance matrix among all the frequencies at different angular scales. We apply this technique to simulated deep imaging data over a 2$\times$100 deg$^2$ sky area from 0.75 $\mu$m to 5 $\mu$m in 9 bands and find that such a "frequency tomography" can successfully reconstruct both the amplitude and spectral shape for representative EoR, IHL and the foreground signals., Comment: 13 pages, 9 figures

The H-alpha line emission is an important probe for a number of fundamental quantities in galaxies, including their number density, star formation rate (SFR) and overall gas content. A new generation of low-resolution intensity mapping probes, e.g. SPHEREx and CDIM, will observe galaxies in H-alpha emission over a large fraction of the sky from the local Universe till a redshift of z ~ 6 to 10, respectively. This will also be the target line for observations by the high-resolution Euclid and WFIRST instruments in the z ~ 0.7 - 2 redshift range. In this paper, we estimate the intensity and power spectra of the H-alpha line in the z ~ 0 - 5 redshift range using observed line luminosity functions (LFs), when possible, and simulations, otherwise. We estimate the significance of our predictions by accounting for the modelling uncertainties (e.g. SFR, extinction, etc.) and observational contamination. We find that Intensity Mapping (IM) surveys can make a statistical detection of the full H-alpha emission between z ~ 0.8 - 5. Moreover, we find that the high-frequency resolution and the sensitivity of the planned CDIM surveys allow for the separation of H-alpha emission from several interloping lines. We explore ways to use the combination of these line intensities to probe galaxy properties. As expected, our study indicates that galaxy surveys will only detect bright galaxies that contribute up to a few percent of the overall H-alpha intensity. However, these surveys will provide important constraints on the high end of the H-alpha LF and put strong constraints on the AGN LF., Comment: Submitted to MNRAS

Large scale filaments, with lengths that can reach tens of Mpc, are the most prominent features in the cosmic web. These filaments have only been observed indirectly through the positions of galaxies in large galaxy surveys or through absorption features in the spectra of high redshift sources. In this study we propose to go one step further and directly detect intergalactic medium filaments through their emission in the HI 21cm line. We make use of high resolution cosmological simulations to estimate the intensity of this emission in low redshift filaments and use it to make predictions for the direct detectability of specific filaments previously inferred from galaxy surveys, in particular the Sloan Digital Sky Survey. Given the expected signal of these filaments our study shows that HI emission from large filaments can be observed by current and next generation radio telescopes. We estimate that gas in filaments of length $l \gtrsim$ 15 $h^{-1}$Mpc with relatively small inclinations to the line of sight ($\lesssim 10^\circ$) can be observed in $\sim40-100$ hours with telescopes such as GMRT or EVLA, potentially providing large improvements over our knowledge of the astrophysical properties of these filaments. Due to their large field of view and sufficiently long integration times, upcoming HI surveys with the Apertif and ASKAP instruments will be able to detect large filaments independently of their orientation and curvature. Furthermore, our estimates indicate that a more powerful future radio telescope like SKA-2 can be used to map most of these filaments, which will allow them to be used as a strong cosmological probe., Comment: 16 pages, 11 figures, Accepted for publication in MNRAS

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., Comment: 14 figures, 4 tables, re-submitted to ApJ after addressing reviewer's comments. Comments welcome

Intensity mapping is now becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the H-alpha 6563AA, [OIII]5007AA, [OII]3727AA and H-beta 4861AA lines at 0.8

We present a systematic study of the intensity mapping technique using updated models for the different emission lines from galaxies and identify which ones are more promising for cosmological studies of the post reionization epoch. We consider the emission of ${\rm Ly\alpha}$, ${\rm H\alpha}$, H$\beta$, optical and infrared oxygen lines, nitrogen lines, CII and the CO rotational lines. We then identify that ${\rm Ly\alpha}$, ${\rm H\alpha}$, OII, CII and the lowest rotational CO lines are the best candidates to be used as IM probes. These lines form a complementary set of probes of the galaxies emission spectra. We then use reasonable experimental setups from current, planned or proposed experiments to access the detectability of the power spectrum of each emission line. Intensity mapping of ${\rm Ly\alpha}$ emission from $z=2$ to 3 will be possible in the near future with HETDEX, while far-infrared lines require new dedicated experiments. We also show that the proposed SPHEREx satellite can use OII and ${\rm H\alpha}$ IM to study the large-scale distribution of matter in intermediate redshifts of 1 to 4. We found that submilimeter experiments with bolometers can have similar performances at intermediate redshifts using CII and CO(3-2)., Comment: 18 pages, 21 figures, 5 tables, published in MNRAS, typos corrected

A large fraction of the baryonic matter in the Universe is located in filaments in the intergalactic medium. However, the low surface brightness of these filaments has not yet allowed their direct detection except in very special regions in the circum-galactic medium (CGM). Here we simulate the intensity and spatial fluctuations in Lyman Alpha ${\rm (Ly\alpha)}$ emission from filaments in the intergalactic medium (IGM) and discuss the prospects for the next generation of space based instruments to detect the low surface brightness universe at UV wavelengths. Starting with a high resolution N-body simulation we obtain the dark matter density fluctuations and associate baryons with the dark matter particles assuming that they follow the same spatial distribution. The IGM thermal and ionization state is set by a model of the UV background and by the relevant cooling processes for a hydrogen and helium gas. The ${\rm Ly\alpha}$ emissivity is then estimated, taking into account recombination and collisional excitation processes. We find that the detection of these filaments through their ${\rm Ly\alpha}$ emission is well in the reach of the next generation of UV space based instruments and so it should be achieved in the next decade. The density field is populated with halos and galaxies and their ${\rm Ly\alpha}$ emission is estimated. Galaxies are treated as foregrounds and so we discuss methods to reduce their contamination from observational maps. Finally, we estimate the UV continuum background as a function of the redshift of the ${\rm Ly\alpha}$ emission line and discuss how this continuum can affect observations., Comment: 12 pages, 7 figures, submitted to MNRAS

We produce simulations of emission of the atomic CII line in large sky fields in order to determine the current prospects for mapping this line during the high redshift Epoch of Reionization. We estimate the CII line intensity, redshift evolution and spatial fluctuations using observational relations between CII emission and the SFR in a galaxy for the frequency range of 200 GHz to 300 GHz. We obtained a frequency averaged intensity of CII emission of ${\rm I_{\rm CII}=(4 \pm 2)\times10^{2}\, Jy\, \rm sr^{-1}}$ in the redshift range $z\, \sim\, 5.3\, -\, 8.5$. Observations of CII emission in this frequency range will suffer contamination from emission lines at lower redshifts, in particular from the CO rotation lines. For the relevant frequency range we estimated the CO contamination (originated in emission from galaxies at $z\, <\, 2.5$), using simulations, to be ${\rm I_{\rm CO} \approx 10^{3}\, Jy \, sr^{-1}}$ and independently confirmed the result based in observational relations. We generated maps as a function of angle and frequency using detailed simulations of the CII and CO emission across several redshifts in order to properly take into account the observational pipeline and light cone effects. In order to reduce the foreground contamination we found that we should mask galaxies below redshifts $\sim 2.5$ with a CO flux in one of the CO(J:2-1) to CO(J:6-5) lines higher than ${\rm 5\times 10^{-22}\, W\ m^{-2}}$ or a AB magnitude lower than ${\rm m_{\rm K}\, =\, 22}$. We estimate that the additional continuum contamination is of the order of ${\rm 10^{5}\, Jy\, sr^{-1}}$. It is also considered the possibility of cross correlating foreground lines with galaxies in order to probe the intensity of the foregrounds., Comment: 19 pages, 14 figures

We calculate the absolute intensity and anisotropies of the Lyman-alpha radiation field present during the epoch of reionization. We consider emission from both galaxies and the intergalactic medium (IGM) and take into account the main contributions to the production of Lyman-alpha photons: recombinations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. We find that the emission from individual galaxies dominates over the IGM with a total Lyman-alpha intensity (times frequency) of about (1.43-3.57)x10^{-8} erg s^{-1} cm^{-2} sr^{-1} at a redshift of 7. This intensity level is low so it is unlikely that the Lyman-\alpha background during reionization can be established by an experiment aiming at an absolute background light measurement. Instead we consider Lyman-\alpha intensity mapping with the aim of measuring the anisotropy power spectrum which has rms fluctuations at the level of 1 x 10^{-16} [erg s^[-1} cm^{-2} sr^{-1}]^2 at a few Mpc scales. These anisotropies could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 \mu m with fields in the order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint Lyman-\alpha emitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-\alpha intensity mapping experiment could provide a viable complimentary approach to probe reionization, when compared to 21 cm experiments, and is likely within experimental reach., Comment: 18 pages, 17 figures

We propose observations of the molecular gas distribution during the era of reionization. At $z \sim $ 6 to 8, $^{12}$CO$(J=1-0)$ line intensity results in a mean brightness temperature of about 0.5 $\mu$K with rms fluctuations of 0.1 $\mu$K at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations can be mapped with an interferometer, similar to existing and planned 21-cm background experiments, but operating at $\sim$ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics and most foregrounds in each of the 21-cm and CO$(1-0)$ line experiments. Additional instruments tuned to higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm-wavelengths targeting 158 $\mu$m CII could further improve the reionization studies with molecular gas. The combined 21-cm and CO line observations has the potential to establish the relative distribution of gas in the inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely connected to the formation of massive stars in these galaxies., Comment: To appear in ApJ Letters