Search

We forecast the ability of future-generation experiments to detect the fine-structure lines of the carbon and oxygen ions, [CII] and [OIII] in intensity mapping (IM) from the Epoch of Reionization ($z \sim 6-8$). Combining the latest empirically derived constraints relating the luminosity of the [OIII] line to the ambient star-formation rate, and using them in conjunction with previously derived estimates for the abundance of [CII] in haloes, we predict the expected auto-correlation IM signal to be observed using new experiments based on the Fred Young Submillimetre Telescope (FYST) and the balloon-borne facility, Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) over $z \sim 5.3 - 7$. We describe how improvements to both the ground-based and balloon-based surveys in the future will enable a cross-correlation signal to be detected at $\sim$ 10-30 $\sigma$ over $z \sim 5.3 - 7$. Finally, we propose a space-based mission targeting the [OIII] 88 and 52 $\mu$m lines along with the [CII] 158 $\mu$m line, configured to enhance the signal-to-noise ratio of cross-correlation measurements. We find that such a configuration can achieve a high-significance detection (hundreds of $\sigma$) in both auto- and cross-correlation modes., Comment: 11 pages, 7 figures, 2 tables; accepted for publication in MNRAS

We use the latest results from Atacama Large Millimetre/submillimetre Array (ALMA) surveys targeting the ionized carbon [CII] 158 $\mu$m and oxygen [OIII] 88 $\mu$m lines, in combination with data-driven predictions for the evolution of neutral hydrogen (HI), to illustrate the prospects for intensity mapping cross-correlations between 21 cm and submillimetre surveys over $z \sim 5-7$. We work with a dataset including the ALPINE and REBELS surveys for [CII] over $z \sim 4.5-7$, and ALMA [OIII] detections over $z \sim 6-9$. The resultant evolution of the [CII] luminosity - halo mass relation is well described by a double power law at high redshifts, with the best-fitting parameters in good agreement with the results of simulations. The data favour secure detections of the auto-power spectrum of [CII] at all redshifts with an enhanced Fred Young Submillimetre Telescope (FYST)-like configuration. Such an experiment, along with the Murchinson Widefield Array (MWA) will be able to measure the 21 cm - [CII] cross-correlation power with a signal-to-noise ratio of a few tens to a few hundreds. We find that a balloon-borne experiment improving upon the Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) should be able to detect the 21 cm - [OIII] cross-correlation with the MWA and the Square Kilometre Array (SKA)-LOW out to $z \sim 7$. Our results have implications for constraining the evolution of luminous sources during the mid-to-end stages of reionization., Comment: 13 pages, 7 figures, 3 tables; version accepted for publication in MNRAS

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}, 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

Intensity mapping (IM) of spectral lines has the potential to revolutionize cosmology by increasing the total number of observed modes by several orders of magnitude compared to the cosmic microwave background (CMB) anisotropies. In this paper, we consider IM of neutral hydrogen (HI) in the redshift range $0 \lesssim z \lesssim 3$ employing a halo model approach where HI is assumed to follow the distribution of dark matter (DM) halos. If a portion of the DM is composed of ultralight axions then the abundance of halos is changed compared to cold dark matter below the axion Jeans mass. With fixed total HI density, $\Omega_{\rm HI}$, assumed to reside entirely in halos, this effect introduces a scale-independent increase in the HI power spectrum on scales above the axion Jeans scale, which our model predicts consistent with N-body simulations. Lighter axions introduce a scale-dependent feature even on linear scales due to its suppression of the matter power spectrum near the Jeans scale. We use the Fisher matrix formalism to forecast the ability of future HI surveys to constrain the axion fraction of DM and marginalize over astrophysical and model uncertainties. We find that a HIRAX-like survey is a very reliable IM survey configuration, being affected minimally by uncertainties due to non-linear scales, while the SKA1MID configuration is the most constraining as it is sensitive to non-linear scales. Including non-linear scales and combining a SKA1MID-like IM survey with the Simons Observatory CMB, the benchmark "fuzzy DM" model with $m_a = 10^{-22}\text{ eV}$ can be constrained at the 10% level. For lighter ULAs this limit improves below 1%, and allows the possibility to test the connection between axion models and the grand unification scale across a wide range of masses., Comment: 16 pages, 14 figures, 2 tables, code to reproduce our results is publicly available at https://github.com/JurekBauer/axion21cmIM.git

We combine the latest observationally motivated constraints on stellar properties in dark matter haloes, along with data-driven predictions for the atomic (HI) and molecular (H$_2$) gas evolution in galaxies, to derive empirical relationships between the build-up of galactic components and their evolution over cosmic time. At high redshift ($z \gtrsim 4$), the frameworks imply that galaxies acquire their cold gas (both atomic and molecular) mostly by accretion, with the fraction of cold gas reaching about 20% of the cosmic baryon fraction. We infer a strong dependence of the star formation rate on the H$_2$ mass, suggesting a near-universal depletion timescale of 0.1-1 Gyr in Milky Way sized haloes (of masses $10^{12} \ M_{\odot}$ at $z = 0$). There is also evidence for a near-universality of the Kennicutt-Schmidt relation across redshifts, with very little dependence on stellar mass, if a constant conversion factor ($\alpha_{\rm CO}$) of CO luminosity to molecular gas mass is assumed. Combining the atomic and molecular gas observations with the stellar build-up illustrates that galactic mass assembly in Milky-Way sized haloes proceeds from smooth accretion at high redshifts, towards becoming merger-dominated at late times ($z \lesssim 0.6$). Our results can be used to constrain numerical simulations of the dominant growth and accretion processes of galaxies over cosmic history., Comment: 8 pages, 6 figures, accepted for publication in MNRAS

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 illustrate a novel signature of black hole binaries from their effect on the kinematics of water maser emission in their environments. With the help of simulations, we establish the condition for clumps to mase based on their coherence lengths calibrated to those of the known maser galaxy NGC 4258. This is then used to identify masing clumps around a binary black hole system, and quantify the kinematic and spectral differences relative to the single black hole case. For some generic circumstances, blue-shifted masers around a binary black hole are found to preferentially follow the Keplerian rotation curve observed in the single black hole case. The redshifted ones, however, are found to visibly deviate from this relation, and also display more scatter with a tendency towards lower absolute values of the velocity along the line-of-sight. The spectrum of the masers as a function of line-of-sight velocity also shows a double peaked structure, reminiscent of recent observations of systems such as Mrk 1. Our results motivate future prospects for identifying binary black hole candidates with the help of water maser emissions., Comment: 6 pages, 11 figures; accepted for publication in MNRAS

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, Comment: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature., Comment: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in ApJ

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., 24 pages + appendix and bibliography (33 pages total), 16 figures, 2 tables; accepted for publication in ApJ

Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop the technology for future experiments, as well as delivering early science products. With 19 receiver channels in a hexagonal focal plane arrangement on a 10.4 m antenna, and an instantaneous 26-34 GHz frequency range with 2 MHz resolution, it is ideally suited to measuring CO($J$=1-0) from $z\sim3$. In this paper we discuss strategies for designing and building the Pathfinder and the challenges that were encountered. The design of the instrument prioritized LIM requirements over those of ancillary science. After a couple of years of operation, the instrument is well understood, and the first year of data is already yielding useful science results. Experience with this Pathfinder will drive the design of the next generations of experiments., Comment: Paper 2 of 7 in series. 27 pages, 28 figures, submitted to ApJ

We illustrate the observability of the end stages of the earliest (Population III) stars at high redshifts $z \gtrsim 10$, using the recently observed transient, GN-z11-flash as an example. We find that the observed spectrum of this transient is consistent with its originating from a shock-breakout in a Population III supernova occurring in the GN-z11 galaxy at $z \sim 11$. The energetics of the explosion indicate a progenitor star of mass $\sim 300 M_{\odot}$ in that galaxy, with of order unity such events expected over an observing timescale of a few years. We forecast the expected number of such transients from $z > 10$ galaxies as a function of their host stellar mass and star formation rate. Our findings are important in the context of future searches to detect and identify the signatures of galaxies at Cosmic Dawn., 12 pages, 1 figure; to appear in Gen. Rel. Grav

We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $\chi^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers., Comment: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to ApJ

The epoch of Cosmic Dawn, when the first stars and galaxies were born, is widely considered the final frontier of observational cosmology today. Mapping the period between Cosmic Dawn and the present-day provides access to more than 90% of the baryonic (normal) matter in the Universe, and unlocks several thousand times more Fourier modes of information than available in today's cosmological surveys. We review the progress in modelling baryonic gas observations as tracers of the cosmological large-scale structure from Cosmic Dawn to the present day. We illustrate how the description of dark matter haloes can be extended to describe baryonic gas abundances and clustering. This innovative approach allows us to fully utilize our current knowledge of astrophysics to constrain cosmological parameters from future observations. Combined with the information content of multi-messenger probes, this will also elucidate the properties of the first supermassive black holes at Cosmic Dawn. We present a host of fascinating implications for constraining physics beyond the $\Lambda$CDM model, including tests of the theories of inflation and the cosmological principle, the effects of non-standard dark matter, and possible deviations from Einstein's general relativity on the largest scales., Comment: Invited review article; 64 pages, 7 figures, 3 tables. v2: typos corrected, references updated, matches published version

The technique of intensity mapping (IM) has emerged as a powerful tool to explore the universe at z < 6. IM measures the integrated emission from sources over a broad range of frequencies, unlocking significantly more information than traditional galaxy surveys. Astrophysical uncertainties, however, constitute an important systematic in our attempts to constrain cosmology with IM. I describe an innovative approach which allows us to fully utilize our current knowledge of astrophysics in order to develop cosmological forecasts from IM. This framework can be used to exploit synergies with other complementary surveys, thereby opening up the fascinating possibility of constraining physics beyond ΛCDM from future IM observations.

The coalescence of massive black hole binaries (with masses $10^4 - 10^7 M_{\odot}$) leads to gravitational wave emission that is detectable out to high redshifts ($z \sim 20$) with the forthcoming LISA observatory. We combine the theoretically derived merger rates for dark matter haloes at various redshifts, with an empirically motivated prescription that connects the mass of a dark matter halo and that of its central black hole. Using the expected constraints on the (chirp or reduced) masses of binary black holes, their mass ratios and redshift uncertainties, we forecast the measurement precision on the occupation fraction, normalization and slope of the black hole mass - halo mass relation at various redshifts, assuming a five-year LISA survey for three different confidence scenarios. We use the expected sizes of the LISA localization ellipses on the sky to estimate the number of electromagnetic counterparts to the gravitational wave sources which are detectable by future wide-field optical surveys, such as LSST., 17 pages, 3 figures; version accepted for publication in JCAP

We forecast astrophysical and cosmological parameter constraints from synergies between 21 cm intensity mapping and wide field optical galaxy surveys (both spectroscopic and photometric) over $z \sim 0-3$. We focus on the following survey combinations in this work: (i) a CHIME-like and DESI-like survey in the northern hemisphere, (ii) an LSST-like and SKA I MID-like survey and (ii) a MeerKAT-like and DES-like survey in the southern hemisphere. We work with the $\Lambda$CDM cosmological model having parameters $\{h, \Omega_m, n_s, \Omega_b, \sigma_8\}$, parameters $v_{c,0}$ and $\beta$ representing the cutoff and slope of the HI-halo mass relation in the previously developed HI halo model framework, and a parameter $Q$ that represents the scale dependence of the optical galaxy bias. Using a Fisher forecasting framework, we explore (i) the effects of the HI and galaxy astrophysical uncertainties on the cosmological parameter constraints, assuming priors from the present knowledge of the astrophysics, (ii) the improvements on astrophysical constraints over their current priors in the three configurations considered, (ii) the tightening of the constraints on the parameters relative to the corresponding HI auto-correlation surveys alone., Comment: 9 pages, 5 figures, 2 tables; accepted for publication in MNRAS

We explore the contributions of Tidal Disruption Events (TDEs) to the flares in Active Galactic Nuclei (AGN) at high redshifts. Using the latest data available from X-ray and optical observations of high-redshift galaxies, in combination with the evolution of their central supermassive black holes, we calculate the contribution of TDE to AGNs as a function of their luminosities. We find that at low redshifts ($z < 1$), a few percent of all AGN with bolometric luminosities $L_{\rm Bol} \lesssim 10^{44}$ erg s$^{-1}$ may be attributable to possible TDEs. However, this fraction can significantly increase at earlier cosmic times, including up to several tens of percent of the population of AGN at $z \gtrsim 3$. TDEs may comprise a significant fraction of the Compton-Thick (CT) AGN population at $z \gtrsim 3$. The above findings motivate further calibration with upcoming X-ray missions and spectroscopic surveys targeting TDE-AGN., 6 pages, 5 figures; accepted for publication in A&A

Square Kilometre Array: et al., arXiv:1810.02680v3, The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy, and spacetime., A. Racanelli has received funding from the People Programme (Marie Curie Actions) of the European Union H2020 Programme under REA grant agreement number 706896 (COSMOFLAGS). Funding for this work was partially provided by the Spanish MINECO under MDM-2014-0369 of ICCUB (Unidad de Excelencia ‘Maria de Maeztu’).

In the local universe, black hole masses have been inferred from the observed increase in the velocities of stars at the centres of their host galaxies. So far, masses of supermassive black holes in the early universe have only been inferred indirectly, using relationships calibrated to their locally observed counterparts. Here, we use the latest observational constraints on the evolution of stellar masses in galaxies to predict, for the first time, that the region of influence of a central supermassive black hole at the epochs where the first galaxies were formed is $\mathit{directly \ resolvable}$ by current and upcoming telescopes. We show that the existence of the black hole can be inferred from observations of the gas or stellar disc out to $> 0.5$ kpc from the host halo at redshifts $z \gtrsim 6$. Such measurements will usher in a new era of discoveries unraveling the formation of the first supermassive black holes based on subarcsecond-scale spectroscopy with the $JWST$, $ALMA$, and the $SKA$. The measured mass distribution of black holes will allow forecasting of the future detection of gravitational waves from the earliest black hole mergers., 9 pages, 3 figures; version accepted for publication in JCAP

We present a novel approach to disentangle two key contributions to the largest-scale anisotropy of the galaxy distribution: (i) the intrinsic dipole due to clustering and anisotropic geometry, and (ii) the kinematic dipole due to our peculiar velocity. Including the redshift and angular size of galaxies, in addition to their fluxes and positions allows us to measure both the direction and amplitude of our velocity independently of the intrinsic dipole of the source distribution. We find that this new approach applied to future galaxy surveys (LSST and Euclid) and a SKA radio continuum survey will allow to measure our velocity ($\beta = v/c$) with a relative error in the amplitude $\sigma(\beta)/\beta \sim (1.3 - 4.5)\%$ and in direction, $\theta_{\beta} \sim 0.9^\circ - 3.9^\circ$, well beyond what can be achieved when analysing only the number count dipole. We also find that galaxy surveys are able to measure the intrinsic large-scale anisotropy with a relative uncertainty of $\lesssim 5\%$ (measurement error, not including cosmic variance). Our method enables two simultaneous tests of the Cosmological Principle: comparing the observations of our peculiar velocity with the CMB dipole, and testing for a significant intrinsic anisotropy on large scales which would indicate effects beyond the standard cosmological model., Comment: 32 pages, 10 figures, 1 table; published in JCAP

Post-reionisation 21cm intensity mapping experiments target the spectral line of neutral hydrogen (HI) resident in dark matter haloes. According to the halo model, these discrete haloes trace the continuous dark matter density field down to a certain scale, which is dependent on the halo physical size. The halo physical size defines an exclusion region which leaves imprints on the statistical properties of HI. We show how the effect of exclusion due to the finite halo size impacts the HI power spectrum, with the physical boundary of the host halo given by the splashback radius. Most importantly, we show that the white noise-like feature that appears in the zero-momentum limit of the power spectrum is exactly cancelled when the finite halo size is taken into consideration. This cancellation in fact applies to all tracers of dark matter density field, including galaxies. Furthermore, we show that the exclusion due to finite halo size leads to a sub-Poissonian noise signature on large scales, consistent with the results from N-body simulations, Comment: 23 pages, 7 figures, one appendix. Version accepted for publication by JCAP

Mapping the unresolved intensity of the 21-cm emission of neutral hydrogen (HI) is now regarded as one the most promising tools for cosmological investigation in the coming decades. Here, we investigate, for the first time, extensions of the standard cosmological model, such as modified gravity and primordial non-Gaussianity, taking self-consistently into account the present constraints on the astrophysics of HI clustering in the treatment of the brightness temperature fluctuations. To understand the boundaries within which results thus obtained can be considered reliable, we examine the robustness of cosmological parameter estimation performed via studies of 21-cm intensity mapping, against our knowledge of the astrophysical processes leading to HI clustering. Modelling of astrophysical effects affects cosmological observables through the relation linking the overall HI mass in a bound object, to the mass of the underlying dark matter halo that hosts it. We quantify the biases in estimates of standard cosmological parameters and those describing modified gravity and primordial non-Gaussianity, that are obtained if one misconceives the slope of the relation between HI mass and halo mass, or the lower virial velocity cut-off for a dark matter halo to be able to host HI. Remarkably, we find that astrophysical uncertainties will not affect searches for primordial non-Gaussianity - one of the strongest science cases for HI intensity mapping - despite the signal being deeply linked to the HI bias., Comment: 14 pages, 8 figures, 2 tables; v4: corrected a typo in the abstract (both arXiv metadata and .pdf file)

The paper outlines research and development needed over the next decade to enable three-dimension mapping of the Universe using the 21 cm line of neutral hydrogen.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

We discuss the possible origin of the double-peaked profiles recently observed in Lyman-Alpha Emitters (LAEs) at the epoch of reionization ($z \gtrsim 6.5$) from obscured active galactic nuclei (AGN). Combining the extent of the Lyman-$\alpha$ near-zones estimated from the blue peak velocity offset in these galaxies, with the ionizing emissivity of quasars at $z \gtrsim 6$, we forecast the intrinsic UV and X-ray luminosities of the AGN needed to give rise to their double-peaked profiles. We also estimate the extent of the obscuration of the AGN by comparing their luminosities to those of similar quasar samples at these epochs. Future X-ray and radio observations, as well as those with the ${James \ Webb \ Space \ Telescope}$, will be valuable tools to test the AGN contribution to the intergalactic-scale ionization zones of high-redshift LAEs., Comment: 4 pages, 1 figure; accepted for publication in A&A Letters

Traditional excursion set based models of H II bubble growth during the epoch of reionization are known to violate photon number conservation, in the sense that the mass fraction in ionized bubbles in these models does not equal the ratio of the number of ionizing photons produced by sources and the number of hydrogen atoms in the intergalactic medium. E.g., for a Planck13 cosmology with electron scattering optical depth $\tau\simeq0.066$, the discrepancy is $\sim15$ per cent for $x_{\rm HII}=0.1$ and $\sim5$ per cent for $x_{\rm HII}=0.5$. We demonstrate that this problem arises from a fundamental conceptual shortcoming of the excursion set approach (already recognised in the literature on this formalism) which only tracks average mass fractions instead of the exact, stochastic source counts. With this insight, we build an approximately photon number conserving Monte Carlo model of bubble growth based on partitioning regions of dark matter into halos. Our model, which is formally valid for white noise initial conditions (ICs), shows dramatic improvements in photon number conservation, as well as substantial differences in the bubble size distribution, as compared to traditional models. We explore the trends obtained on applying our algorithm to more realistic ICs, finding that these improvements are robust to changes in the ICs. Since currently popular semi-numerical schemes of bubble growth also violate photon number conservation, we argue that it will be worthwhile to pursue new, explicitly photon number conserving approaches. Along the way, we clarify some misconceptions regarding this problem that have appeared in the literature., Comment: 12 pages, 7 figures; v2: revised to match version accepted in MNRAS

We combine available constraints on the local CII 158 $\mu$m line luminosity function from galaxy observations (Hemmati et al. 2017), with the evolution of the star-formation rate density and the recent CII intensity mapping measurement in Pullen et al. (2018, assuming detection), to derive the evolution of the CII luminosity - halo mass relation over $z \sim 0-6$. We develop convenient fitting forms for the evolution of the CII luminosity - halo mass relation, and forecast constraints on the CII intensity mapping power spectrum and its associated uncertainty across redshifts. We predict the sensitivities to detect the power spectrum for upcoming PIXIE-, STARFIRE-, EXCLAIM-, CONCERTO-, TIME- and CCAT-p-like surveys, as well as possible future intensity mapping observations with the ALMA facility., Comment: 10 pages, 9 figures, 2 tables; version accepted for publication in MNRAS

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 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_\text{CO}(M_\text{halo})$ model of Li et al. (2016). 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 both small and large luminosities, while the PS is more sensitive at intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between $L_\text{CO} = 10^4-10^7L_\odot$, 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., Comment: 13 pages, 5 figures. As accepted to ApJ

The evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe ($z \sim $ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts., Invited talk at the IAU Symposium 333 "Peering towards Cosmic Dawn", Dubrovnik, October 2-6, 2017; to appear in the proceedings, eds. Vibor Jelic and Thijs van der Hulst [6 pages, 2 figures]

Using hydrodynamical simulations, we explore the use of the mean and percentiles of the curvature distribution function to recover the equation of state of the high-$z$ ($2 < z < 4$) intergalactic medium (IGM). We find that the mean and percentiles of the absolute curvature distribution exhibit tight correlation with the temperatures measured at respective characteristic overdensities $\bar{\Delta}_i$'s at each redshift. Hence, they provide complementary probes of the same underlying temperature-density distribution, and can in principle be used to simultaneously recover both parameters $T_0$ and $\gamma$ of the IGM effective equation of state. We quantify the associated errors in the recovered parameters $T_0$ and $\gamma$ from the intrinsic scatter in the characteristic overdensities and the uncertainties in the curvature measurement., Comment: 5 pages, 4 figures, 1 table; version accepted for publication in MNRAS Letters

A unique feature of gravity is its ability to control the information accessible to any specific observer. We quantify the notion of cosmic information (‘CosmIn’) for an eternal observer in the universe. Demanding the finiteness of CosmIn requires the universe to have a late-time accelerated expansion. Combining the introduction of CosmIn with generic features of the quantum structure of spacetime (e.g., the holographic principle), we present a holistic model for cosmology. We show that (i) the numerical value of the cosmological constant, as well as (ii) the amplitude of the primordial, scale invariant, perturbation spectrum can be determined in terms of a single free parameter, which specifies the energy scale at which the universe makes a transition from a pre-geometric phase to the classical phase. For a specific value of the parameter, we obtain the correct results for both (i) and (ii). This formalism also shows that the quantum gravitational information content of spacetime can be tested using precision cosmology., Physics Letters B, 773, ISSN:0370-2693, ISSN:0031-9163, ISSN:1873-2445

We propose a method to extract the unlensed, intrinsic CMB temperature and polarization power spectra from the observed (i.e., lensed) spectra. Using a matrix inversion technique, we demonstrate how one can reconstruct the intrinsic CMB power spectra directly from lensed data for both flat sky and full sky analyses. The delensed spectra obtained by the technique are calibrated against the Code for Anisotropies in the Microwave Background (CAMB) using WMAP 7-year best-fit data and applied to WMAP 9-year unbinned data as well. In principle, our methodology may help in subtracting out the E-mode lensing contribution in order to obtain the intrinsic B-mode power.

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., Comment: 19 pages + appendix (31 pages total), 16 figures, 6 tables; accepted for publication in ApJ

We extend the results of previous analyses towards constraining the abundance and clustering of post-reionization ($z \sim 0-5$) neutral hydrogen (HI) systems using a halo model framework. We work with a comprehensive HI dataset including the small-scale clustering, column density and mass function of HI galaxies at low redshifts, intensity mapping measurements at intermediate redshifts and the UV/optical observations of Damped Lyman Alpha (DLA) systems at higher redshifts. We use a Markov Chain Monte Carlo (MCMC) approach to constrain the parameters of the best-fitting models, both for the HI-halo mass relation and the HI radial density profile. We find that a radial exponential profile results in a good fit to the low-redshift HI observations, including the clustering and the column density distribution. The form of the profile is also found to match the high-redshift DLA observations, when used in combination with a three-parameter HI-halo mass relation and a redshift evolution in the HI concentration. The halo model predictions are in good agreement with the observed HI surface density profiles of low-redshift galaxies, and the general trends in the the impact parameter and covering fraction observations of high-redshift DLAs. We provide convenient tables summarizing the best-fit halo model predictions., 13 pages, 12 figures, 5 tables; version accepted for publication in MNRAS

We describe a combined halo model to constrain the distribution of neutral hydrogen (H I) in the post-reionization universe. We combine constraints from the various probes of H I at different redshifts: the low-redshift 21-cm emission line surveys, intensity mapping experiments at intermediate redshifts, and the Damped Lyman-Alpha (DLA) observations at higher redshifts. We use a Markov Chain Monte Carlo approach to combine the observations and place constraints on the free parameters in the model. Our best-fitting model involves a relation between neutral hydrogen mass MHI and halo mass M with a non-unit slope, and an upper and a lower cutoff. We find that the model fits all the observables but leads to an underprediction of the bias parameter of DLAs at z ∼ 2.3. We also find indications of a possible tension between the H I column density distribution and the mass function of H I-selected galaxies at z ∼ 0. We provide the central values of the parameters of the best-fitting model so derived. We also provide a fitting form for the derived evolution of the concentration parameter of H I in dark matter haloes, and discuss the implications for the redshift evolution of the H I–halo mass relation., Monthly Notices of the Royal Astronomical Society, 464 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

The neutral hydrogen (HI) content of dark matter haloes forms an intermediate state in the baryon cycle that connects the hot shock-heated gas and cold star-forming gas in haloes. Measurement of the relationship between HI mass and halo mass therefore puts important constraints on galaxy formation models. We combine radio observations of HI in emission at low redshift ($z\sim 0$) with optical/UV observations of HI in absorption at high redshift ($1, Comment: 11 pages, 8 figures, 1 table; published in MNRAS. v3: minor typo corrected; results unchanged

The expression for the electromagnetic field of a charge moving along an arbitrary trajectory is obtained in a direct, elegant, and Lorentz invariant manner without resorting to more complicated procedures such as differentiation of the Lienard-Wiechert potentials. The derivation uses arguments based on Lorentz invariance and a physically transparent expression originally due to J.J.Thomson for the field of a charge that experiences an impulsive acceleration., The following article has been accepted by the American Journal of Physics. After it is published, it will be found at http://scitation.aip.org/ajp; 12 pages, 1 figure

We review the analytical prescriptions in the literature to model the 21-cm (emission line surveys/intensity mapping experiments) and Damped Lyman-Alpha (DLA) observations of neutral hydrogen (HI) in the post-reionization universe. While these two sets of prescriptions have typically been applied separately for the two probes, we attempt to connect these approaches to explore the consequences for the distribution and evolution of HI across redshifts. We find that a physically motivated, 21-cm based prescription, extended to account for the DLA observables provides a good fit to the majority of the available data, but cannot accommodate the recent measurement of the clustering of DLAs at $z \sim 2.3$. This highlights a tension between the DLA bias and the 21-cm measurements, unless there is a very significant change in the nature of HI-bearing systems across redshifts 0-3. We discuss the implications of our findings for the characteristic host halo masses of the DLAs and the power spectrum of 21-cm intensity fluctuations., 8 pages, 5 figures, 1 table; version accepted for publication in MNRAS

Gravity controls the amount of information that is accessible to any specific observer. We quantify the notion of cosmic information (‘CosmIn’) for the case of an eternal observer in the universe. Demanding the finiteness of CosmIn requires the universe to have a late-time accelerated expansion phase. Combined with some generic features of the quantum structure of spacetime, this leads to the determination of (i) the numerical value of the cosmological constant, as well as (ii) the amplitude of the primordial, scale invariant perturbation spectrum in terms of a single free parameter, which specifies the energy scale at which the universe makes a transition from a pre-geometric phase to the classical phase. This formalism also shows that the quantum gravitational information content of spacetime can be tested by using precision cosmology.

Monthly Notices of the Royal Astronomical Society, 447 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Observations indicate that our universe is characterized by a late-time accelerating phase, possibly driven by a cosmological constant $\Lambda$, with the dimensionless parameter $\Lambda L_P^2 \simeq 10^{-122}$, where $L_P = (G \hbar /c^3)^{1/2}$ is the Planck length. In this review, we describe how the emergent gravity paradigm provides a new insight and a possible solution to the cosmological constant problem. After reviewing the necessary background material, we identify the necessary and sufficient conditions for solving the cosmological constant problem. We show that these conditions are naturally satisfied in the emergent gravity paradigm in which (i) the field equations of gravity are invariant under the addition of a constant to the matter Lagrangian and (ii) the cosmological constant appears as an integration constant in the solution. The numerical value of this integration constant can be related to another dimensionless number (called CosMIn) that counts the number of modes inside a Hubble volume that cross the Hubble radius during the radiation and the matter dominated epochs of the universe. The emergent gravity paradigm suggests that CosMIn has the numerical value $4 \pi$, which, in turn, leads to the correct, observed value of the cosmological constant. Further, the emergent gravity paradigm provides an alternative perspective on cosmology and interprets the expansion of the universe itself as a quest towards holographic equipartition. We discuss the implications of this novel and alternate description of cosmology., Comment: Invited review to appear in Int. Jour. Mod. Phys. D; 48 pages; 5 figures

Using hydrodynamical simulations coupled to a radiative transfer code, we study the additional heating effects in the intergalactic medium (IGM) produced by $z\sim 6$ quasars in their near-zones. If helium is predominantly in HeII to begin with, both normalization ($T_0$) and slope ($\gamma$) of the IGM effective equation-of-state get modified by the excess ionization from the quasars. Using the available constraints on $T_0$ at $z\sim6$, we discuss implications for the nature and epoch of HI and HeII reionization. We study the extent of the HeIII region as a function of quasar age and show, for a typical inferred age of $z \sim 6$ quasars (i.e. $\sim 10^8$ yrs), it extends up to 80% of the HI proximity region. For these long lifetimes, the heating effects can be detected even when all the HI lines from the proximity region are used. Using the flux and curvature probability distribution functions (PDFs), we study the statistical detectability of heating effects as a function of initial physical conditions in the IGM. For the present sample size, cosmic variance dominates the flux PDF. The curvature statistics is more suited to capturing the heating effects beyond the cosmic variance, even if the sample size is half of what is presently available., Comment: 20 pages, 15 figures, 3 tables; version accepted for publication in MNRAS

We provide a detailed treatment and comparison of the weak lensing effects due to large-scale structure (LSS), or scalar density perturbations and those due to gravitational waves(GW) or tensor perturbations, on the temperature and polarization power spectra of the Cosmic Microwave Background (CMB). We carry out the analysis both in real space by using the correlation function method, as well as in the spherical harmonic space. We find an intriguing similarity between the lensing kernels associated with LSS lensing and GW lensing. It is found that the lensing kernels only differ in relative negative signs and their form is very reminiscent of even and odd parity bipolar spherical harmonic coefficients. Through a numerical study of these lensing kernels, we establish that lensing due to GW is more efficient at distorting the CMB spectra as compared to LSS lensing, particularly for the polarization power spectra. Finally we argue that the CMB B-mode power spectra measurements can be used to place interesting constraints on GW energy densities., Comment: Typos corrected. Published : Phys. Rev. D 88, 063507, 4 September 2013