Primordial non-Gaussianity (PNG) is a key probe of the origins of primordial fluctuations in the early universe. It has been shown that multi-tracer measurements of large-scale structure can produce high-precision measurements of PNG. Future line intensity mapping surveys are well-suited to these measurements owing to their ability to rapidly survey large volumes and access the large scales at which PNG becomes important. In this paper, we explore for the first time how multi-tracer PNG measurements with intensity mapping surveys depend on the sub-galactic scale physics which drives line emission. We consider an example cross-correlation between CO maps, and find that the choice of astrophysical models has a substantial impact on $f_{\mathrm{NL}}$ measurements at sufficiently high signal-to-noise. This in effect creates a coupling between horizon-scale PNG measurements and the molecular cloud-scale interstellar medium. We discuss how these effects depend on noise level and survey design. We further find that the cross-correlation shot noise, an effect nearly unique to intensity mapping measurements, plays an important role in multi-tracer analyses and cannot be neglected.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

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

Line-intensity mapping (LIM) of emission form star-forming galaxies can be used to measure the baryon acoustic oscillation (BAO) scale as far back as the epoch of reionization. This provides a standard cosmic ruler to constrain the expansion rate of the Universe at redshifts which cannot be directly probed otherwise. In light of growing tension between measurements of the current expansion rate using the local distance ladder and those inferred from the cosmic microwave background, extending the constraints on the expansion history to bridge between the late and early Universe is of paramount importance. Using a newly derived methodology to robustly extract cosmological information from LIM, which minimizes the inherent degeneracy with unknown astrophysics, we show that present and future experiments can gradually improve the measurement precision of the expansion rate history, ultimately reaching percent-level constraints on the BAO scale. Specifically, we provide detailed forecasts for the SPHEREx satellite, which will target the H$\alpha$ and Lyman-$\alpha$ lines, and for the ground-based COMAP instrument---as well as a future stage-3 experiment---that will target the CO rotational lines. Besides weighing in on the so-called Hubble tension, reliable LIM cosmic rulers can enable wide-ranging tests of dark matter, dark energy and modified gravity.

Line-intensity mapping (LIM) provides a promising way to probe cosmology, reionization and galaxy evolution. However, its sensitivity to cosmology and astrophysics at the same time is also a nuisance. Here we develop a comprehensive framework for modelling the LIM power spectrum, which includes redshift space distortions and the Alcock-Paczynski effect. We then identify and isolate degeneracies with astrophysics so that they can be marginalized over. We study the gains of using the multipole expansion of the anisotropic power spectrum, providing an accurate analytic expression for their covariance, and find a 10%-60% increase in the precision of the baryon acoustic oscillation scale measurements when including the hexadecapole in the analysis. We discuss different observational strategies when targeting other cosmological parameters, such as the sum of neutrino masses or primordial non-Gaussianity, finding that fewer and wider bins are typically more optimal. Overall, our formalism facilitates an optimal extraction of cosmological constraints robust to astrophysics.

21 cm intensity mapping has arisen as a powerful probe of the high-redshift universe, but its potential is limited by extremely bright foregrounds and high source confusion. In this Letter, we propose a new analysis which can help solve both problems. From the combination of an intensity map with an overlapping galaxy survey we construct a new one-point statistic which is unbiased by foregrounds and contains information left out of conventional analyses. We show that our method can measure the HI mass function with unprecedented precision using observations similar to recent 21 cm detections.

Line intensity maps (LIMs) are in principle sensitive to a large amount of information about faint, distant galaxies which are invisible to conventional surveys. However, actually extracting that information from a confused, foreground-contaminated map can be challenging. In this work we present the first application of convolutional neural network (CNN) to directly determine the underlying luminosity function of a LIM, including a treatment of extragalactic foregrounds and instrumental noise. We apply the CNN to simulations of mock Carbon Monoxide (CO) line intensity maps similar to those which will be produced by the currently-active COMAP experiment. We evaluate the trained CNN on a number of noise scenarios in order to determine how robust the network predictions are for application to realistic data. We find that, in the ideal case where the mock data capture all of the features of the real data, the CNN performs comparably to or better than conventional analyses. However, the network's accuracy degrades considerably when tested on signals and systematics outside of those it was trained on. For both intensity mapping and cosmology as a whole, this motivates a broad-based study of whether simulated data can ever be generated with sufficient detail to realize the enormous potential of machine learning methods.

Current models of galaxy formation require star formation in high-mass galaxies to be limited by poorly understood mechanisms of quasar feedback. Feedback processes can be studied by examining the molecular gas content of AGN hosts through the CO rotational ladder, but the complexity of these observations means that current data are limited to only extremely CO-bright objects. Upcoming CO intensity mapping experiments offer an opportunity for a less biased probe of quasar feedback. By correlating intensity maps with spectroscopic AGN surveys, we can obtain a measurement of the mean CO luminosity of a large population of quasars simultaneously. We show that experiments like COMAP, CCAT-prime, and CONCERTO have enough sensitivity to detect this cross-correlation if existing AGN observations are representative of the whole population, and to place interesting upper limits if they are not. Future surveys will be able to increase the precision of these measurements by orders of magnitude, allowing detailed studies of quasar properties across a wide range of cosmic history.

Fast radio bursts (FRBs) are, as the name implies, short and intense pulses of radiation at wavelengths of roughly one metre. FRBs have extremely high brightness temperatures, which points to a coherent source of radiation. The energy of a single burst ranges from $10^{36}$ to $10^{39}$ erg. At the high end of the energy range, FRBs have enough energy to unbind an earth-sized planet, and even at the low end, there is enough energy to vaporise and unbind the atmosphere and the oceans. We therefore propose that FRBs are signatures of an artificial terraformer, capable of eradicating life on another planet, or even destroy the planet entirely. The necessary energy can be harvested from Wolf-Rayet stars with a Dyson sphere ($\sim 10^{38}$ erg s$^{-1}$) , and the radiation can be readily produced by astrophysical masers. We refer to this mechanism as Volatile Amplification of a Destructive Emission of Radiation (VADER). We use the observational information to constrain the properties of the apparatus. We speculate that the non-repeating FRBs are low-energy pulses used to exterminate life on a single planet, but leaving it otherwise intact, and that the stronger repeating FRB is part of an effort to destroy multiple objects in the same solar system, perhaps as a preventative measure against panspermia. In this picture, the persistent synchrotron source associated with the first repeating FRB arises from the energy harvesting process. Finally we propose that Oumuamua might have resulted from a destruction of a planet in this manner.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

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.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

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.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

The distances to fast radio bursts (FRBs) are crucial for understanding their underlying engine, and for their use as cosmological probes. In this paper, we provide three statistical estimates of the distance to ASKAP FRBs. First, we show that the number of events of similar luminosity in ASKAP does not scale as distance cubed, as one would expect, when directly using the observed dispersion measure (DM) to infer distance. Second, by comparing the average DMs of FRBs observed with different instruments, we estimated the average redshift of ASKAP FRBs to be $z\sim 0.01$ using CHIME and ASKAP, and $z\lesssim0.07$ using Parkes and ASKAP. Both values are much smaller than the upper limit $z\sim0.3$ estimated directly from the DM. Third, we cross-correlate the locations of the ASKAP FRBs with existing large-area redshift surveys, and see a 3$\sigma$ correlation with the 2MASS Redshift Survey and a 5$\sigma$ correlation with the HI Parkes All Sky Survey at $z\sim0.007$. This corresponds well with the redshift of the most likely host galaxy of ASKAP FRB 171020, which is at $z=0.00867$. These arguments combined suggest an extremely nearby origin of ASKAP FRBs and a local environment with accumulated electrons that contribute a DM of several hundred pc/cm$^3$, which should be accounted for in theoretical models.

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 the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

We present the most detailed computation to date of the 21-cm global signal and fluctuations at $z\gtrsim 500$. Our calculations include a highly precise estimate of the Wouthuysen-Field (WF) effect and the first explicit calculation of the impact of free-free processes, the two dominant components of the signal at $z\gtrsim 800$. We implement a new high-resolution Ly$\alpha$ radiative transfer calculation, coupled to a state-of-the-art primordial recombination code. Using these tools, we find a global signal from 21-cm processes alone of roughly 0.01mK at $z\sim1000$, slightly larger than it would be without the WF effect, but much weaker than previous estimates including this effect. We also find that this signal is swamped by a smooth $1-2$ mK signal due to free-free absorption at high redshift by the partially ionized gas along the line of sight. In addition, we estimate the amplitude of 21-cm fluctuations, of order $\sim 10^{-7}$ mK at $z\sim1000$. Unfortunately, we find that due to the brightness of the low-frequency sky, these fluctuations will not be observable beyond $z\sim$ a few hundred by even extremely futuristic observations. The 21 cm fluctuations are exponentially suppressed at higher redshifts by the large free-free optical depth, making this the ultimate upper redshift limit for 21-cm surveys.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

The study of molecular gas is crucial for understanding star formation, feedback, and the broader ecosystem of a galaxy as a whole. However, we have limited understanding of its physics and distribution in all but the nearest galaxies. We present a new technique for studying the composition and distribution of molecular gas in high-redshift galaxies inaccessible to existing methods. Our proposed approach is an extension of carbon monoxide intensity mapping methods, which have garnered significant experimental interest in recent years. These intensity mapping surveys target the 115 GHz $^{12}$CO (1-0) line, but also contain emission from the substantially fainter 110 GHz $^{13}$CO (1-0) transition. The method leverages the information contained in the $^{13}$CO line by cross-correlating pairs of frequency channels in an intensity mapping survey. Since $^{13}$CO is emitted from the same medium as the $^{12}$CO, but saturates at a much higher column density, this cross-correlation provides valuable information about both the gas density distribution and isotopologue ratio, inaccessible from the $^{12}$CO alone. Using a simple model of these molecular emission lines, we show that a future intensity mapping survey can constrain the abundance ratio of these two species and the fraction of emission from optically thick regions to order $\sim30\%$. These measurements cannot be made by traditional CO observations, and consequently the proposed method will provide unique insight into the physics of star formation, feedback, and galactic ecology at high redshifts.

We examine how future gravitational-wave measurements from merging black holes (BHs) can be used to infer the shape of the black-hole mass function, with important implications for the study of star formation and evolution and the properties of binary BHs. We model the mass function as a power law, inherited from the stellar initial mass function, and introduce lower and upper mass cutoff parameterizations in order to probe the minimum and maximum BH masses allowed by stellar evolution, respectively. We initially focus on the heavier BH in each binary, to minimize model dependence. Taking into account the experimental noise, the mass measurement errors and the uncertainty in the redshift-dependence of the merger rate, we show that the mass function parameters, as well as the total rate of merger events, can be measured to <10% accuracy within a few years of advanced LIGO observations at its design sensitivity. This can be used to address important open questions such as the upper limit on the stellar mass which allows for BH formation and to confirm or refute the currently observed mass gap between neutron stars and BHs. In order to glean information on the progenitors of the merging BH binaries, we then advocate the study of the two-dimensional mass distribution to constrain parameters that describe the two-body system, such as the mass ratio between the two BHs, in addition to the merger rate and mass function parameters. We argue that several years of data collection can efficiently probe models of binary formation, and show, as an example, that the hypothesis that some gravitational-wave events may involve primordial black holes can be tested. Finally, we point out that in order to maximize the constraining power of the data, it may be worthwhile to lower the signal-to-noise threshold imposed on each candidate event and amass a larger statistical ensemble of BH mergers.

In the next few years, intensity-mapping surveys that target lines such as CO, Ly$\alpha$, and CII stand to provide powerful probes of high-redshift astrophysics. However, these line emissions are highly non-Gaussian, and so the typical power-spectrum methods used to study these maps will leave out a significant amount of information. We propose a new statistic, the probability distribution of voxel intensities, which can access this extra information. Using a model of a CO intensity map at $z\sim3$ as an example, we demonstrate that this voxel intensity distribution (VID) provides substantial constraining power beyond what is obtainable from the power spectrum alone. We find that a future survey similar to the planned COMAP Full experiment could constrain the CO luminosity function to order $\sim10\%$. We also explore the effects of contamination from continuum emission, interloper lines, and gravitational lensing on our constraints and find that the VID statistic retains significant constraining power even in pessimistic scenarios.

We demonstrate how cosmic star-formation history can be measured with one-point statistics of carbon-monoxide intensity maps. Using a P(D) analysis, the luminosity function of CO-emitting sources can be inferred from the measured one-point intensity PDF. The star-formation rate density (SFRD) can then be obtained, at several redshifts, from the CO luminosity density. We study the effects of instrumental noise, line foregrounds, and target redshift, and obtain constraints on the CO luminosity density of order 10%. We show that the SFRD uncertainty is dominated by that of the model connecting CO luminosity and star formation. For pessimistic estimates of this model uncertainty, we obtain an error of order 50% on SFRD for surveys targeting redshifts between 2 and 7 with reasonable noise and foregrounds included. However, comparisons between intensity maps and galaxies could substantially reduce this model uncertainty. In this case our constraints on SFRD at these redshifts improve to roughly 5-10%, which is highly competitive with current measurements.

We address the problem of line confusion in intensity mapping surveys and explore the possibility to mitigate line foreground contamination by progressively masking the brightest pixels in the observed map. We consider experiments targeting CO(1-0) at $z=3$, Ly$\alpha$ at $z=7$, and CII at $z=7$, and use simulated intensity maps, which include both clustering and shot noise components of the signal and possible foregrounds, in order to test the efficiency of our method. We find that for CO and Ly$\alpha$ it is quite possible to remove most of the foreground contribution from the maps via only 1%-3% pixel masking. The CII maps will be more difficult to clean, however, due to instrumental constraints and the high-intensity foreground contamination involved. While the masking procedure sacrifices much of the astrophysical information present in our maps, we demonstrate that useful cosmological information in the targeted lines can be successfully retrieved.

We present a study of the feasibility of an intensity-mapping survey targeting the 115 GHz CO(1-0) rotational transition at $z\sim3$. We consider four possible models and estimate the spatial and angular power spectra of CO fluctuations predicted by each of them. The frequency bandwidths of most proposed CO intensity mapping spectrographs are too small to use the Limber approximation to calculate the angular power spectrum, so we present an alternative method for calculating the angular power spectrum. The models we consider span two orders of magnitude in signal amplitude, so there is a significant amount of uncertainty in the theoretical predictions of this signal. We then consider a parameterized set of hypothetical spectrographs designed to measure this power spectrum and predict the signal-to-noise ratios expected under these models. With the spectrographs we consider we find that three of the four models give an SNR greater than 10 within one year of observation. We also study the effects on SNR of varying the parameters of the survey in order to demonstrate the importance of carefully considering survey parameters when planning such an experiment.

We study the oscillations and stability of self-gravitating cylindrically symmetric fluid systems and collisionless systems. This is done by studying small perturbations to the equilibrium system and finding the normal modes, using methods similar to those used in astroseismology. We find that there is a single sequence of purely radial modes that become unstable if the adiabatic exponent is less than 1. Nonradial modes can be divided into p modes, which are stable and pressure-driven, and g modes, which are are gravity driven. The g modes become unstable if the adiabatic exponent is greater than the polytrope index. These modes are analogous to the modes of a spherical star, but their behavior is somewhat different because a cylindrical geometry has less symmetry than a spherical geometry. This implies that perturbations are classified by a radial quantum number, an azimuthal quantum number, and wavelength in the z direction, which can become arbitrarily large. We find that decreasing this wavelength increases the frequency of stable modes and increases the growth rate of unstable modes. We use use variational arguments to demonstrate that filaments of collisionless matter with ergodic distribution functions are stable to purely radial perturbations, and that filaments with ergodic power-law distribution functions are stable to all perturbations.