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.

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.

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.

We explore the contributions of Tidal Disruption Events (TDEs) as potential triggers for Active Galactic Nuclei (AGN). Using the latest data available from X-ray and optical observations of high-redshift galaxies, as well as the evolution of their central supermassive black holes, we calculate the extent to which TDEs trigger AGN 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 TDE triggers. However, this fraction can significantly increase at earlier cosmic times, reaching several tens of percent at $z > 2$, including about 20-50% of the observed population of Changing Look Quasars (CLQs). TDEs may comprise a significant fraction of the Compton-Thick (CT) AGN population at $z > 2$. The above findings motivate further calibration from upcoming X-ray missions and spectroscopic surveys targeting TDE-AGN.

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.

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project.

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 $\Lambda$CDM from future IM observations.

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.

Line-intensity mapping of the 21cm line is a powerful probe of large scale structure at z<6, tracing large-scale structure via neutral hydrogen content that is found within galaxies. In principle, it enables cost-efficient surveys of the matter distribution up to z~6, unlocking orders of magnitude more modes for observational cosmology. Canada has been a traditional leader in this field, having led the first detections of the cosmological 21cm signal via cross-correlations with optical galaxy surveys and having constructed the Canadian Hydrogen Intensity Mapping Experiment (CHIME). The field is now entering a new era where data is abundant, allowing studies in how to overcome systematics to be tackled in an empirical, head-on fashion. In the next few years, this will produce the first detection of the 21cm auto power spectrum, which will pave the way towards a large suite of scientific possibilities. These potentially include precision measurements on the dark energy equation of state and other LCDM parameters, constraints on how HI mass traces dark matter, a detection of neutrino effects on large-scale structure, and the use of 21cm lensing to further constrain cosmology. To turn these promising directions into reality, we recommend a sustained program of investment in 21cm cosmology, starting with funding for the Canadian Hydrogen Observatory and Radio transient Detector (CHORD), followed by small-scale development efforts targeting next-generation hardware and sustained support for theory and technical staff support. Additionally, Canada should invest in complementary line-intensity mapping efforts (such as with CO or [CII] lines) and maintain participation in next-generation international efforts such as the Packed Ultra-wideband Mapping Array (PUMA) and the Square Kilometre Array (SKA).

Mapping the unresolved intensity of 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 the present constraints on the astrophysics of HI clustering in the treatment of the brightness temperature fluctuations self-consistently into account. 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 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.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

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 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.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

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.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We use the results of previous work building a halo model formalism for the distribution of neutral hydrogen, along with experimental parameters of future radio facilities, to place forecasts on astrophysical and cosmological parameters from next generation surveys. We consider 21 cm intensity mapping surveys conducted using the BINGO, CHIME, FAST, TianLai, MeerKAT and SKA experimental configurations. We work with the 5-parameter cosmological dataset of {$\Omega_m, \sigma_8, h, n_s, \Omega_b$} assuming a flat $\Lambda$CDM model, and the astrophysical parameters {$v_{c,0}, \beta$} which represent the cutoff and slope of the HI- halo mass relation. We explore (i) quantifying the effects of the astrophysics on the recovery of the cosmological parameters, (ii) the dependence of the cosmological forecasts on the details of the astrophysical parametrization, and (iii) the improvement of the constraints on probing smaller scales in the HI power spectrum. For an SKA I MID intensity mapping survey alone, probing scales up to $\ell_{\rm max} = 1000$, we find a factor of $1.1 - 1.3$ broadening in the constraints on $\Omega_b$ and $\Omega_m$, and of $2.4 - 2.6$ on $h$, $n_s$ and $\sigma_8$, if we marginalize over astrophysical parameters without any priors. However, even the prior information coming from the present knowledge of the astrophysics largely alleviates this broadening. These findings do not change significantly on considering an extended HIHM relation, illustrating the robustness of the results to the choice of the astrophysical parametrization. Probing scales up to $\ell_{\rm max} = 2000$ improves the constraints by factors of 1.5-1.8. The forecasts improve on increasing the number of tomographic redshift bins, saturating, in many cases, with 4 - 5 redshift bins. We also forecast constraints for intensity mapping with other experiments, and draw similar conclusions.

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.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We compile available constraints on the carbon monoxide (CO) 1-0 luminosity functions and abundances at redshifts 0-3. This is used to develop a data driven halo model for the evolution of the CO galaxy abundances and clustering across intermediate redshifts. It is found that the recent constraints from the CO Power Spectrum Survey ($z \sim 3$; Keating et al. 2016), when combined with existing observations of local galaxies ($z \sim 0$; Keres et al. 2003), lead to predictions which are consistent with the results of smaller surveys at intermediate redshifts ($z \sim 1-2$). We provide convenient fitting forms for the evolution of the CO luminosity - halo mass relation, and estimates of the mean and uncertainties in the CO power spectrum in the context of future intensity mapping experiments.

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.

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<z<4$) to derive constraints on the evolution of the HI-mass halo-mass (HIHM) relation from redshift $z=4$ to $z=0$. We find that one can model the HIHM relation similar to the stellar-mass halo-mass (SHM) relation at $z \sim 0$. At $z=0$, haloes with mass $10^{11.7}$ M$_\odot$ have the highest HI mass fraction ($\sim 1\%$), which is about four times smaller than their stellar mass fraction. We model the evolution of the HIHM relation in a manner similar to that of the SHM relation. Combining this parameterisation with a redshift- and mass-dependent modified Navarro-Frenk-White (NFW) profile for the HI density within a halo, we draw constraints on the evolution of the HIHM relation from the observed HI column density, incidence rate, and clustering bias at high redshift. We compare these findings with results from hydrodynamical simulations and other approaches in the literature, and find the models to be consistent with each other at the $68\%$ confidence level.

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.

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.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We describe a combined halo model to constrain the distribution of neutral hydrogen (HI) in the post-reionization universe. We combine constraints from the various probes of HI 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 (MCMC) approach to combine the observations and place constraints on the free parameters in the model. Our best-fit model involves a relation between neutral hydrogen mass $M_{\rm HI}$ 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 \sim 2.3$. We also find indications of a possible tension between the HI column density distribution and the mass function of HI-selected galaxies at $z\sim 0$. We provide the central values of the parameters of the best-fit model so derived. We also provide a fitting form for the derived evolution of the concentration parameter of HI in dark matter haloes, and discuss the implications for the redshift evolution of the HI-halo mass relation.

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.

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.

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.

Mapping of the neutral hydrogen (HI) 21-cm intensity fluctuations across redshifts promises a novel and powerful probe of cosmology. The neutral hydrogen gas mass density, $\Omega_{\rm HI}$ and bias parameter, $b_{\rm HI}$ are key astrophysical inputs to the HI intensity fluctuation power spectrum. We compile the latest theoretical and observational constraints on $\Omega_{\rm HI}$ and $b_{\rm HI}$ at various redshifts in the post-reionization universe. Constraints are incorporated from galaxy surveys, HI intensity mapping experiments, damped Lyman-$\alpha$ system observations, theoretical prescriptions for assigning HI to dark matter halos, and the results of numerical simulations. Using a minimum variance interpolation scheme, we obtain the predicted uncertainties on the HI intensity fluctuation power spectrum across redshifts 0-3.5 for three different confidence scenarios. We provide a convenient tabular form for the interpolated values of $\Omega_{\rm HI}$, $b_{\rm HI}$ and the HI power spectrum amplitude and their uncertainties. We discuss the consequences for the measurement of the power spectrum by current and future intensity mapping experiments.

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.

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.

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.

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.

The current acceleration of the universe can be modeled in terms of a cosmological constant. We show that the extremely small value of \Lambda L_P^2 ~ 3.4 x 10^{-122}, the holy grail of theoretical physics, can be understood in terms of a new, dimensionless, conserved number CosMIn (N), which counts the number of modes crossing the Hubble radius during the three phases of evolution of the universe. Theoretical considerations suggest that N ~ 4\pi. This single postulate leads us to the correct, observed numerical value of the cosmological constant! This approach also provides a unified picture of cosmic evolution relating the early inflationary phase to the late-time accelerating phase.