We present a new estimator for the cross-correlation signal between line intensity maps to probe the Epoch of Reionization. The proposed estimator is the hitherto neglected antisymmetric component of the cross-correlation, under the exchange of line-of-sight positions. We consider the cross-correlation between HI and CO fluctuations, and forecast the improvement in precision on reionization parameters when the antisymmetric contribution is accounted for. As a way to break the degeneracy between astrophysics and cosmology in the intensity mapping power spectrum, we study the ratio between the antisymmetric and symmetric components. While our results depend on the highly uncertain astrophysical modelling, we show that in most standard scenarios including the antisymmetric contribution as a complementary probe can lead to a significant gain in information.

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.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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.

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.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

Large-amplitude density perturbations may have collapsed during the radiation dominated epoch of the Universe to form primordial black holes (PBHs). There are several constraints to the abundance of PBHs that stem from evaporation or gravitational effects. Due to the connection between primordial perturbations and the formation of PBHs, constraints on the present-day abundance of PBHs can be translated into limits on the primordial curvature power spectrum. We introduce several new observational and forecasted constraints to the amplitude of the primordial power spectrum and incorporate in our analysis uncertainties in the critical overdensity for collapse and considerations of ellipsoidal collapse. Our results provide the most stringent limits from PBHs on the primordial curvature power spectrum on small scales.

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.

Heat transfer between baryons and millicharged dark matter has been invoked as a possible explanation for the anomalous 21-cm absorption signal seen by EDGES. Prior work has shown that the solution requires that millicharged particles make up only a fraction $(m_\chi/{\rm MeV})\ 0.0115\% \lesssim f \lesssim 0.4\%$ of the dark matter and that their mass $m_\chi$ and charge $q_\chi$ have values $0.1 \lesssim (m_\chi/{\rm MeV})\lesssim 10$ and $10^{-6} \lesssim (q_\chi/e)\lesssim 10^{-4}$. Here we show that such particles come into chemical equilibrium before recombination, and so are subject to a constraint on the effective number $N_{\rm eff}$ of relativistic degrees of freedom, which we update using Planck 2018 data. We moreover determine the precise relic abundance $f$ that results for a given mass $m_\chi$ and charge $q_\chi$ and incorporate this abundance into the constraints on the millicharged-dark-matter solution to EDGES. With these two results, the solution is ruled out if the relic abundance is set by freeze-out.

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.

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

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.

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.

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.

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

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 LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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.

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

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.

The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$.

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.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

It has been suggested that primordial black holes (PBHs) of roughly 30 solar masses could make up the dark matter and if so, might account for the recent detections by LIGO involving binary black holes in this mass range. It has also been argued that the super-massive black holes (SMBHs) that reside at galactic centers may be surrounded by extremely-dense dark-matter (DM) spikes. Here we show that the rate for PBH mergers in these spikes may well exceed the merger rate considered before in galactic dark-matter halos, depending on the magnitudes of two competing effects on the DM spikes: depletion of PBHs due to relaxation and replenishment due to PBHs in loss cone. This may provide a plausible explanation for the current rate of detection of mergers of 30-solar-mass black holes, even if PBHs make up a subdominant contribution to the dark matter. The gravitational-wave signals from such events will always originate in galactic centers, as opposed to those from halos, which are expected to have little correlation with luminous-galaxy positions.

The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $\Lambda$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for different design sensitivities. We also include a multi-tracer analysis, distinguishing between star-forming galaxies and galaxies with an active galactic nucleus, to further improve EMU's potential. We find that EMU could measure the dark energy equation of state parameters around 35\% more precisely than existing constraints, and that the constraints on $f_{\rm NL}$ and modified gravity parameters will improve up to a factor $\sim2$ with respect to Planck and redshift space distortions measurements. With this work we demonstrate the promising potential of EMU to contribute to our understanding of the Universe.

Different inflationary models predict oscillatory features in the primordial power spectrum. These can leave an imprint on both the cosmic microwave background (CMB) and the large-scale structure (LSS) of our Universe, that can be searched for with current data. Inspired by the axion-monodromy model of inflation, we search for primordial oscillations that are logarithmic in wavenumber, using both CMB data from the Planck satellite and LSS data from the WiggleZ galaxy survey. We find that, within our search range for the new oscillation parameters (amplitude, frequency and phase), both CMB-only and CMB+LSS data yield the same best-fit oscillation frequency of $\log_{10}\omega=1.5$, improving the fit over $\Lambda$CDM by $\Delta\chi^2\! =\! -9$ and $\Delta\chi^2\! =\! -13$ (roughly corresponding to $2\sigma$ and $2.8\sigma$ significance), respectively. Bayesian evidence for the log-oscillation model versus $\Lambda$CDM indicates a very slight preference for the latter. Future CMB and LSS data will further probe this scenario.

Compact dark matter has been efficiently constrained in the M <~ 10 M_sun mass range by null searches for microlensing of stars in nearby galaxies. Here we propose to probe the mass range M >~ 10 M_sun by seeking echoes in gamma-ray-burst light curves induced by strong lensing. We show that strong gravitational lensing of gamma ray bursts (GRBs) by massive compact halo objects (MACHOs) generates superimposed GRB images with a characteristic time delay of >~ 1 ms for M >~ 10 M_sun. Using dedicated simulations to capture the relevant phenomenology of the GRB prompt emission, we calculate the signal-to-noise ratio required to detect GRB lensing events as a function of the flux ratio and time delay between the lensed images. We then analyze existing data from the Fermi/GBM and Swift/BAT instruments to assess their constraining power on the compact dark matter fraction f_DM. We find that this data is noise limited, and therefore localization-based masking of background photons is a key ingredient. Future observatories with better sensitivity will be able to probe down to the f_ DM >~ 1% level across the 10 M_sun <~ M <~ 1000 M_sun mass range.

Ultra-light hidden-photon dark matter produces an oscillating electric field in the early Universe plasma, which in turn induces an electric current in its ionized component whose dissipation results in heat transfer from the dark matter to the plasma. This will affect the global 21cm signal from the Dark Ages and Cosmic Dawn. In this work we focus on the latter, in light of the reported detection by the EDGES collaboration of an absorption signal at frequencies corresponding to redshift z~17. By measuring the 21cm global signal, a limit can be placed on the amount of gas heating, and thus the kinetic mixing strength $\varepsilon$ between the hidden and ordinary photons can be constrained. Our inferred 21cm bounds on $\varepsilon$ in the mass range $10^{-23}\,{\rm eV}\lesssim m_\chi\lesssim10^{-13}\,{\rm eV}$ are the strongest to date.

The Laser Interferometer Space Antenna (LISA) gravitational-wave (GW) observatory will be limited in its ability to detect mergers of binary black holes (BBHs) in the stellar-mass range. A future ground-based detector network, meanwhile, will achieve by the LISA launch date a sensitivity that ensures complete detection of all mergers within a volume $>\!\mathcal{O}(10)\,{\rm Gpc}^{3}$. We propose a method to use the information from the ground to revisit the LISA data in search for sub-threshold events. By discarding spurious triggers that do not overlap with the ground-based catalogue, we show that the signal-to-noise threshold $\rho_{\rm LISA}$ employed in LISA can be significantly lowered, greatly boosting the detection rate. The efficiency of this method depends predominantly on the rate of false-alarm increase when the threshold is lowered and on the uncertainty in the parameter estimation for the LISA events. As an example, we demonstrate that while all current LIGO BBH-merger detections would have evaded detection by LISA when employing a standard $\rho_{\rm LISA}=8$ threshold, this method will allow us to easily (possibly) detect an event similar to GW150914 (GW170814) in LISA. Overall, we estimate that the total rate of stellar-mass BBH mergers detected by LISA can be boosted by a factor $\sim\!4$ ($\gtrsim\!8$) under conservative (optimistic) assumptions. This will enable new tests using multi-band GW observations, significantly aided by the greatly increased lever arm in frequency.

We investigate the hypothesis that Coulomb-type interactions between dark matter (DM) and baryons explain the anomalously low 21cm brightness-temperature minimum at redshift z ~ 17 that was recently measured by the EDGES experiment. In particular, we reassess the validity of the scenario where a small fraction of the total DM is millicharged, focusing on newly derived constraints from Planck 2015 cosmic microwave background (CMB) data. Crucially, the CMB power spectrum is sensitive to DM-baryon scattering if the fraction of interacting DM is larger than (or comparable to) the fractional uncertainty in the baryon energy density. Meanwhile, there is a mass-dependent lower limit on the fraction for which the required interaction to cool the baryons sufficiently is so strong that it drives the interacting-DM temperature to the baryon temperature prior to their decoupling from the CMB. If this occurs as early as recombination, the cooling saturates. We precisely determine the viable parameter space for millicharged DM, and find that only a fraction (m_chi/MeV) 0.0115% <~ f <~ 0.4% of the entire DM content, and only for DM-particle masses between 0.5 MeV - 35 MeV, can be charged at the level needed to marginally explain the anomaly, without violating limits from SLAC, CMB, Big-Bang nucleosynthesis (BBN), or stellar and SN1987A cooling. In reality, though, we demonstrate that at least moderate fine tuning is required to both agree with the measured absorption profile and overcome various astrophysical sources of heating. Finally, we point out that a ~0.4% millicharged DM component which is tightly coupled to the baryons at recombination may resolve the current 2-sigma tension between the BBN and CMB determinations of the baryon energy density. Future CMB-S4 measurements will be able to probe this scenario directly.

We perform an improved cosmic microwave background (CMB) analysis to search for dark matter-proton scattering with a momentum-transfer cross section of the form $\sigma_0 v^n$ for $n=-2$ and $n=-4$. In particular, we present a new and robust prescription for incorporating the relative bulk velocity between the dark matter and baryon fluids into the standard linear Boltzmann calculation. Using an iterative procedure, we self-consistently include the effects of the bulk velocities in a cosmology in which dark matter interacts with baryons. With this prescription, we derive CMB bounds on the cross section, excluding $\sigma_0 > 2.3 \times 10^{-33}~\mathrm{cm}^2$ for $n=-2$ and $\sigma_0 > 1.7 \times 10^{-41}~\mathrm{cm}^2$ for $n=-4$ at $95\%$ confidence, for dark matter masses below 10 MeV. Furthermore, we investigate how these constraints change when only a subcomponent of dark matter is interacting. We show that Planck limits vanish if $\lesssim 0.4\%$ of dark matter is tightly coupled to baryons. We discuss the implications of our results for present and future cosmological observations.

There is growing evidence that intermediate-mass black holes (IMBHs), defined here as having a mass in the range M=500-10^5 Msun, are present in the dense centers of certain globular clusters (GCs). Gravitational waves (GWs) from their mergers with other IMBHs or with stellar BHs in the cluster are mostly emitted in frequencies <10 Hz, which unfortunately is out of reach for current ground-based observatories such as advanced LIGO (aLIGO). Nevertheless, we show that aLIGO measurements can be used to efficiently probe one of the possible formation mechanisms of IMBHs in GCs, namely a runaway merger process of stellar seed BHs. In this case, aLIGO will be sensitive to the lower-mass rungs of the merger ladder, ranging from the seed BH mass to masses >~50-300 Msun, where the background from standard mergers is expected to be very low. Assuming this generic IMBH formation scenario, we calculate the mass functions that correspond to the limiting cases of possible merger trees. Based on estimates for the number density of GCs and taking into account the instrumental sensitivity, we show that current observations do not effectively limit the occupation fraction f_occ of IMBHs formed by runaway mergers of stellar BHs in GCs. However, we find that a six-year run of aLIGO at design sensitivity will be able to probe down to f_occ<3% at a 99.9% confidence level, either finding evidence for this formation mechanism, or necessitating others if the fraction of GCs that harbor IMBHs is higher.

Primordial black holes (PBHs) have long been a candidate for the elusive dark matter (DM), and remain poorly constrained in the ~20-100 Msun mass range. PBH binaries were recently suggested as the possible source of LIGO's first detections. In this paper, we thoroughly revisit existing estimates of the merger rate of PBH binaries. We compute the probability distribution of orbital parameters for PBH binaries formed in the early Universe, accounting for tidal torquing by all other PBHs, as well as standard large-scale adiabatic perturbations. We then check whether the orbital parameters of PBH binaries formed in the early Universe can be significantly affected between formation and merger. Our analytic estimates indicate that the tidal field of halos and interactions with other PBHs, as well as dynamical friction by unbound standard DM particles, do not do significant work on nor torque PBH binaries. We estimate the torque due to baryon accretion to be much weaker than previous calculations, albeit possibly large enough to significantly affect the eccentricity of typical PBH binaries. We also revisit the PBH-binary merger rate resulting from gravitational capture in present-day halos, accounting for Poisson fluctuations. If binaries formed in the early Universe survive to the present time, as suggested by our analytic estimates, they dominate the total PBH merger rate. Moreover, this merger rate would be orders of magnitude larger than LIGO's current upper limits if PBHs make a significant fraction of the dark matter. As a consequence, LIGO would constrain ~10-300 Msun PBHs to constitute no more than ~1% of the dark matter. To make this conclusion fully robust, though, numerical study of several complex astrophysical processes - such as the formation of the first PBH halos and how they may affect PBH binaries, as well as the accretion of gas onto an extremely eccentric binary - is needed.

Primordial black holes (PBHs) have long been suggested as a candidate for making up some or all of the dark matter in the Universe. Most of the theoretically possible mass range for PBH dark matter has been ruled out with various null observations of expected signatures of their interaction with standard astrophysical objects. However, current constraints are significantly less robust in the 20 M_sun < M_PBH < 100 M_sun mass window, which has received much attention recently, following the detection of merging black holes with estimated masses of ~30 M_sun by LIGO and the suggestion that these could be black holes formed in the early Universe. We consider the potential of advanced LIGO (aLIGO) operating at design sensitivity to probe this mass range by looking for peaks in the mass spectrum of detected events. To quantify the background, which is due to black holes that are formed from dying stars, we model the shape of the stellar-black-hole mass function and calibrate its amplitude to match the O1 results. Adopting very conservative assumptions about the PBH and stellar-black-hole merger rates, we show that ~5 years of aLIGO data can be used to detect a contribution of >20 M_sun PBHs to dark matter down to f_PBH<0.5 at >99.9% confidence level. Combined with other probes that already suggest tension with f_PBH=1, the obtainable independent limits from aLIGO will thus enable a firm test of the scenario that PBHs make up all of dark matter.

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.

Single-field slow-roll inflation predicts a nearly scale-free power spectrum of perturbations, as observed at the scales accessible to current cosmological experiments. This spectrum is slightly red, showing a tilt $(1-n_s)\sim 0.04$. A direct consequence of this tilt are nonvanishing runnings $\alpha_s=\mathrm d n_s/\mathrm d\log k$, and $\beta_s=\mathrm d\alpha_s/\mathrm d\log k$, which in the minimal inflationary scenario should reach absolute values of $10^{-3}$ and $10^{-5}$, respectively. In this work we calculate how well future surveys can measure these two runnings. We consider a Stage-4 (S4) CMB experiment and show that it will be able to detect significant deviations from the inflationary prediction for $\alpha_s$, although not for $\beta_s$. Adding to the S4 CMB experiment the information from a WFIRST-like or a DESI-like survey improves the sensitivity to the runnings by $\sim$ 20\%, and 30\%, respectively. A spectroscopic survey with a billion objects, such as the SKA, will add enough information to the S4 measurements to allow a detection of $\alpha_s=10^{-3}$, required to probe the single-field slow-roll inflationary paradigm. We show that only a very-futuristic interferometer targeting the dark ages will be capable of measuring the minimal inflationary prediction for $\beta_s$. The results of other probes, such as a stochastic background of gravitational waves observable by LIGO, the Ly-$\alpha$ forest, and spectral distortions, are shown for comparison. Finally, we study the claims that large values of $\beta_s$, if extrapolated to the smallest scales, can produce primordial black holes of tens of solar masses, which we show to be easily testable by the S4 CMB experiment.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

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.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

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.

The possibility that part of the dark matter is made of massive compact halo objects (MACHOs) remains poorly constrained over a wide range of masses, and especially in the $20-100\, M_\odot$ window. We show that strong gravitational lensing of extragalactic fast radio bursts (FRBs) by MACHOs of masses larger than $\sim20\,M_\odot$ would result in repeated FRBs with an observable time delay. Strong lensing of an FRB by a lens of mass $M_L$ induces two images, separated by a typical time delay $\sim$ few $\times(M_L/30\, M_\odot)$ milliseconds. Considering the expected FRB detection rate by upcoming experiments, such as CHIME, of $10^4$ FRBs per year, we should observe from tens to hundreds of repeated bursts yearly, if MACHOs in this window make up all the dark matter. A null search for echoes with just $10^4$ FRBs, would constrain the fraction $f_{\rm DM}$ of dark matter in MACHOs to $f_{\rm DM}\lesssim 0.08$ for $M_L\gtrsim 20\,M_\odot$.

We demonstrate that observations lacking reliable redshift information, such as photometric and radio continuum surveys, can produce robust measurements of cosmological parameters when empowered by clustering-based redshift estimation. This method infers the redshift distribution based on the spatial clustering of sources, using cross-correlation with a reference dataset with known redshifts. Applying this method to the existing SDSS photometric galaxies, and projecting to future radio continuum surveys, we show that sources can be efficiently divided into several redshift bins, increasing their ability to constrain cosmological parameters. We forecast constraints on the dark-energy equation-of-state and on local non-gaussianity parameters. We explore several pertinent issues, including the tradeoff between including more sources versus minimizing the overlap between bins, the shot-noise limitations on binning, and the predicted performance of the method at high redshifts. Remarkably, we find that, once this technique is implemented, constraints on dynamical dark energy from the SDSS imaging catalog can be competitive with, or better than, those from the spectroscopic BOSS survey and even future planned experiments. Further, constraints on primordial non-Gaussianity from future large-sky radio-continuum surveys can outperform those from the Planck CMB experiment, and rival those from future spectroscopic galaxy surveys. The application of this method thus holds tremendous promise for cosmology.

It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.

We consider the possibility that tidal disruption events (TDEs) caused by supermassive black holes (SMBHs) in nearby galaxies can account for the ultra-high-energy cosmic-ray (UHECR) hotspot reported recently by the Telescope Array (TA) and the warm spot by Pierre Auger Observatory (PAO). We describe the expected cosmic-ray signal from a TDE and derive the constraints set by the timescale for dispersion due to intergalactic magnetic fields and the accretion time of the SMBH. We find that TDEs in M82 can explain the hotspot detected by the TA regardless of whether the UHECRs are composed of protons or heavier nuclei. We then check for consistency of the hot and warm spots from M82 and Cen A with the full-sky isotropic signal from all SMBHs within the GZK radius. This analysis applies to any scenario in which the hot/warm spots are real and due to M82 and Cen A, regardless of whether TDEs are the source of UHECRs. We find that the isotropic flux implied by the luminosity density inferred from M82 and Cen A is bigger than that observed by roughly an order of magnitude, but we provide several possible explanations, including the possibility of a local overdensity and the possibility of intermediate-mass nuclei in UHECRs, to resolve the tension.

We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.

We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources.

We explore the physics and observational consequences of tidal compression events (TCEs) of dark-matter clumps (DMCs) by supermassive black holes (SMBHs). Our analytic calculations show that a DMC approaching a SMBH much closer than the tidal radius undergoes significant compression along the axis perpendicular to the orbital plane, shortly after pericenter passage. For DMCs composed of self-annihilating dark-matter particles, we find that the boosted DMC density and velocity dispersion lead to a flaring of the annihilation rate, most pronounced for a velocity- dependent annihilation cross section. If the end products of the annihilation are photons, this results in a gamma-ray flare, detectable (and possibly already detected) by the Fermi telescope for a range of model parameters. If the end products of dark-matter annihilation are relativistic electrons and positrons and the local magnetic field is large enough, TCEs of DMCs can lead to flares of synchrotron radiation. Finally, TCEs of DMCs lead to a burst of gravitational waves, in addition to the ones radiated by the orbital motion alone, and with a different frequency spectrum. These transient phenomena provide interesting new avenues to explore the properties of dark matter.

The auto-correlation between two members of a galaxy population is symmetric under the interchange of the two galaxies being correlated. The cross-correlation between two different types of galaxies, separated by a vector $\bf{r}$, is not necessarily the same as that for a pair separated by $-\bf{r}$. Local anisotropies in the two-point cross-correlation function may thus indicate a specific direction which when mapped as a function of position trace out a vector field. This vector field can then be decomposed into longitudinal and transverse components, and those transverse components written as positive- and negative-helicity components. A locally asymmetric cross-correlation of the longitudinal type arises naturally in halo clustering, even with Gaussian initial conditions, and could be enhanced with local-type non-Gaussianity. Early-Universe scenarios that introduce a vector field may also give rise to such effects. These antisymmetric cross-correlations also provide a new possibility to seek a preferred cosmic direction correlated with the hemispherical power asymmetry in the cosmic microwave background and to seek a preferred location associated with the CMB cold spot. New ways to seek cosmic parity breaking are also possible.

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.

In the standard inflationary scenario, primordial perturbations are adiabatic. The amplitudes of most types of isocurvature perturbations are generally constrained by current data to be small. If, however, there is a baryon-density perturbation that is compensated by a dark-matter perturbation in such a way that the total matter density is unperturbed, then this compensated isocurvature perturbation (CIP) has no observable consequence in the cosmic microwave background (CMB) at linear order in the CIP amplitude. Here we search for the effects of CIPs on CMB power spectra to quadratic order in the CIP amplitude. An analysis of the Planck temperature data leads to an upper bound $\Delta_{\rm rms}^2 \leq 7.1\times 10^{-3}$, at the 68\% confidence level, to the variance $\Delta_{\rm rms}^2$ of the CIP amplitude. This is then strengthened to $\Delta_{\rm rms}^2\leq 5.0\times 10^{-3}$ if Planck small-angle polarization data are included. A cosmic-variance-limited CMB experiment could improve the $1\sigma$ sensitivity to CIPs to $\Delta^2_{\rm rms} \lesssim 9\times 10^{-4}$. It is also found that adding CIPs to the standard $\Lambda$CDM model can improve the fit of the observed smoothing of CMB acoustic peaks just as much as adding a non-standard lensing amplitude.

The search for the curl component (B mode) in the cosmic microwave background (CMB) polarization induced by inflationary gravitational waves is described. The canonical single-field slow-roll model of inflation is presented, and we explain the quantum production of primordial density perturbations and gravitational waves. It is shown how these gravitational waves then give rise to polarization in the CMB. We then describe the geometric decomposition of the CMB polarization pattern into a curl-free component (E mode) and curl component (B mode) and show explicitly that gravitational waves induce B modes. We discuss the B modes induced by gravitational lensing and by Galactic foregrounds and show how both are distinguished from those induced by inflationary gravitational waves. Issues involved in the experimental pursuit of these B modes are described, and we summarize some of the strategies being pursued. We close with a brief discussion of some other avenues toward detecting/characterizing the inflationary gravitational-wave background.

We explore the effects of elastic scattering between dark matter and baryons on the 21-cm signal during the dark ages. In particular, we consider a dark-matter---baryon interaction with a cross section of the form $\sigma = \sigma_0 v^{-4}$, in which case the effect of the drag force between the dark matter and baryon fluids grows with time. We show that, as opposed to what was previously thought, this effect heats up the baryons due to the relative velocity between dark matter and baryons. This creates an additional source of fluctuations, which can potentially make interactions easier to detect by 21-cm measurements than by using the cosmic microwave background and the Lyman-$\alpha$ forest. Our forecasts show that the magnitude of the cross section can be probed to $\sigma_0\sim 3\times 10^{-42}$ cm$^2$ for $m_{\chi}\ll 1$ GeV and $\sigma_0\sim 2 \times 10^{-41}\ (m_{\chi}/10\, \rm GeV)$ cm$^2$ for $m_{\chi}\gg 1$ GeV with next generation experiments, and improved to $\sigma_0\sim 4\times 10^{-44}$ cm$^2$ for $m_{\chi} \ll 1$ GeV and $\sigma_0\sim 4 \times 10^{-43}\ (m_{\chi}/10\, \rm GeV)$ cm$^2$ for $m_{\chi}\gg 1$ GeV with futuristic experiments.

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.

A preferred method to detect the curl-component, or B-mode, signature of inflationary gravitational waves (IGWs) in the cosmic microwave background (CMB) polarization, in the absence of foregrounds and lensing, is a prolonged integration over a single patch of sky of a few square degrees. In practice, however, foregrounds abound and the sensitivity to B modes can be improved considerably by finding the region of sky cleanest of foregrounds. The best strategy to detect B modes thus involves a tradeoff between exploration (to find lower-foreground patches) and exploitation (through prolonged integration). This problem is akin to the multi-armed bandit (MAB) problem in probability theory, wherein a gambler faces a series of slot machines with unknown winning odds and must develop a strategy to maximize his/her winnings with some finite number of pulls. While the optimal MAB strategy remains to be determined, a number of algorithms have been developed in an effort to maximize the winnings. Here, based on this resemblance, we tackle the search for IGW B modes with single frequency experiments in the presence of spatially-varying foregrounds by developing adaptive survey strategies to optimize the sensitivity to IGW B modes. We demonstrate, using realistic foreground models and taking lensing-induced B modes into account, that adaptive experiments can substantially improve the upper bound on the tensor-to-scalar ratio (by factors of 2--3 in single frequency experiments, and possibly even more). Similar techniques can be applied to other surveys, including 21-cm measurements of signatures of the epoch of reionization, searches for a stochastic primordial gravitational wave background, deep-field imaging by the James Webb Space Telescope or various radio interferometers, and transient follow-up searches.

The Planck satellite has identified several patches of sky with low polarized dust emission, obvious targets for searches for the cosmic-microwave-background (CMB) B-mode signal from inflationary gravitational waves. Still, given the Planck measurement uncertainties, the polarized dust foregrounds in these different candidate patches may differ by an order of magnitude or more. Here we show that a brief initial experiment to map these candidate patches more deeply at a single high frequency can efficiently zero in on the cleanest patch(es) and thus improve significantly the sensitivity of subsequent B-mode searches. A ground-based experiment with current detector technology operating at >~220 GHz for 3 months can efficiently identify a low-dust-amplitude patch and thus improve by up to a factor 2 or 3 on the sensitivity to cosmic B modes of the subsequent lower-frequency deep integration. A balloon experiment with current detector sensitivities covering the set of patches and operating at ~350 GHz can reach a similar result in less than two weeks. This strategy may prove crucial in accessing the smallest gravitational-wave signals possible in large-field inflation. The high-frequency data from this exploratory experiment should also provide valuable foreground templates to subsequent experiments that integrate on any of the candidate patches explored.

We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross-correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely without interacting with the intervening matter, and hence its late-time ISW signature should correlate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature brightness instead of the CMB would thus be a further check of the detection of the ISW effect, measured by different instruments at different frequencies and suffering from different systematics. We also study the ISW effect on the photons that are scattered by HI clouds. We show that a detection of the unscattered photons is achievable with planned radio arrays, while one using scattered photons will require advanced radio interferometers, either an extended version of the planned Square Kilometre Array or futuristic experiments such as a lunar radio array.

Recent developments in the search for inflationary gravitational waves in the cosmic microwave background (CMB) polarization motivate the search for new diagnostics to distinguish the Galactic foreground contribution to B modes from the cosmic signal. We show that B modes from these foregrounds should exhibit a local hexadecapolar departure in power from statistical isotropy (SI). We present a simple algorithm to search for a uniform SI violation of this sort, as may arise in a sufficiently small patch of sky. We then show how to search for these effects if the orientation of the SI violation varies across the survey region, as is more likely to occur in surveys with more sky coverage. If detected, these departures from Gaussianity would indicate some level of Galactic foreground contamination in the B-mode maps. Given uncertainties about foreground properties, though, caution should be exercised in attributing a null detection to an absence of foregrounds.

Previous claims of significant evidence for mirror-parity in the large-scale cosmic microwave background (CMB) data from the Wilkinson Microwave Anisotropy Probe (WMAP) experiment have been recently echoed in the first study of isotropy and statistics of CMB data from Planck. We revisit these claims with a careful analysis of the latest data available. We construct statistical estimators in both harmonic and pixel space, test them on simulated data with and without mirror-parity symmetry, apply different Galactic masks, and study the dependence of the results on arbitrary choices of free parameters. We confirm that the data exhibit evidence for odd mirror-parity at a significance which reaches as high as ~ 99 per cent C.L., under some circumstances. However, given the inherent biases in the pixel-based statistic and the dependence of both pixel and harmonic space statistics on the particular form of Galactic masking and other a-posteriori choices, we conclude that these results are not in significant tension with the predictions of the concordance cosmological model.

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 investigate the anisotropy in cosmic microwave background Planck maps due to the coupling between its beam asymmetry and uneven scanning strategy. Introducing a pixel space estimator based on the temperature gradients, we find a highly significant (~20 \sigma) preference for these to point along ecliptic latitudes. We examine the scale dependence, morphology and foreground sensitivity of this anisotropy, as well as the capability of detailed Planck simulations to reproduce the effect, which is crucial for its removal, as we demonstrate in a search for the weak lensing signature of cosmic defects.

We study the planar regime of curvature perturbations for single field inflationary models in an axially symmetric Bianchi I background. In a theory with standard scalar field action, the power spectrum for such modes has a pole as the planarity parameter goes to zero. We show that constraints from back reaction lead to a strong lower bound on the planarity parameter for high-momentum planar modes and use this bound to calculate the signal-to-noise ratio of the anisotropic power spectrum in the CMB, which in turn places an upper bound on the Hubble scale during inflation allowed in our model. We find that non-Gaussianities for these planar modes are enhanced for the flattened triangle and the squeezed triangle configurations, but show that the estimated values of the f_NL parameters remain well below the experimental bounds from the CMB for generic planar modes (other, more promising signatures are also discussed). For a standard action, f_NL from the squeezed configuration turns out to be larger compared to that from the flattened triangle configuration in the planar regime. However, in a theory with higher derivative operators, non-Gaussianities from the flattened triangle can become larger than the squeezed configuration in a certain limit of the planarity parameter.

An extremely large void and a cosmic texture are two possible explanations for the cold spot seen in the cosmic microwave background. We investigate how well these two hypotheses can be tested with weak lensing of 21-cm fluctuations from the epoch of reionization (EoR) measured with the Square Kilometer Array (SKA). While the void explanation for the cold spot can be tested with SKA, given enough observation time, the texture scenario requires significantly prolonged observations, at the highest frequencies that correspond to the EoR, over the field of view containing the cold spot.

We discuss the prospects to measure galaxy-cluster properties via weak lensing of 21-cm fluctuations from the dark ages and the epoch of reionization (EOR). We choose as a figure of merit the smallest cluster mass detectable through such measurements. We construct the minimum-variance quadratic estimator for the cluster mass based on lensing of 21-cm fluctuations at multiple redshifts. We discuss the tradeoff between frequency bandwidth, angular resolution, and number of redshift shells available for a fixed noise level for the radio detectors. Observations of lensing of the 21-cm background from the dark ages will be capable of detecting M>~10^12 Msun/h mass halos, but will require futuristic experiments to overcome the contaminating sources. Next-generation radio measurements of 21-cm fluctuations from the EOR will, however, have the sensitivity to detect galaxy clusters with halo masses M>~10^13 Msun/h, given enough observation time (for the relevant sky patch) and collecting area to maximize their resolution capabilities.

Studies of the peculiar velocity bulk flow based on different tools and datasets have been consistent so far in their estimation of the direction of the flow, which also happens to lie in close proximity to several features identified in the cosmic microwave background, providing motivation to use new compilations of type-Ia supernovae measurements to pinpoint it with better accuracy and up to higher redshift. Unfortunately, the peculiar velocity field estimated from the most recent Union2.1 compilation suffers from large individual errors, poor sky coverage and low redshift-volume density. We show that as a result, any naive attempt to calculate the best-fit bulk flow and its significance will be severely biased. Instead, we introduce an iterative method which calculates the amplitude and the scatter of the direction of the best-fit bulk flow as deviants are successively removed and take into account the sparsity of the data when estimating the significance of the result. Using 200 supernovae up to a redshift of z=0.2, we find that while the amplitude of the bulk flow is marginally consistent with the value expected in a LCDM universe given the large bias, the scatter of the direction is significantly low (at >= 99.5 C.L.) when compared to random simulations, supporting the quest for a cosmological origin.

We continue to investigate possible signatures of a pre-inflationary anisotropic phase in two-point and three point correlation functions of the curvature perturbation for high-momentum modes which exit the horizon well after isotropization. The late time dynamics of these modes is characterized by a non-Bunch Davies vacuum state which encodes all the information about initial anisotropy in the background space-time. We observe that, unlike the non-planar momenta, there exist regimes of planar momenta for which scale invariance of the power spectrum is strongly broken. This regime of planar momenta gives rise to enhanced non-Gaussianity in certain squeezed triangle configurations, although the enhancement of the $f_{NL}$ parameter is limited by the breakdown of linear perturbation theory at "exact planarity". Finally, we demonstrate that for the range of planar modes for which scale invariance of the power spectrum is preserved, non-Gaussianity in the curvature perturbation spectrum is naturally constrained to be extremely small.

We consider the effects of fields with suddenly changing mass on the inflationary power spectra. In this context, when a field becomes light, it will be excited. This process contributes to the tensor power spectrum. We compute these effects in a gauge-invariant manner, where we use a novel analytical method for evaluating the corrections to the tensor spectrum due to these excitations. In the case of a scalar field, we show that the net impact on the tensors is small as long as the perturbative expansion is valid. Thus, in these scenarios, measurement of tensor modes is still in one-to-one correspondence with the Hubble scale.

We search for a direction in the sky that exhibits parity symmetry under reflections through a plane. We use the natural estimator, which compares the power in even and odd $\ell+m$ multipoles, and apply minimal blind masking of outliers to the ILC map in order to avoid large errors in the reconstruction of multipoles. The multipoles of the cut sky are calculated both naively and by using the covariance inversion method and we estimate the significance of our results using $\Lambda$CDM simulations. Focusing on low multipoles, $2\leq \ell \leq \ell_{\max}$ with $\ell_{\max}=5,6$ or even 7, we find two perpendicular directions of even and odd parity in the map. While the even parity direction does not appear significant, the odd direction is quite significant -- at least a $3.6\sigma$ effect.

We find a unique direction in the CMB sky around which giant rings have an anomalous mean temperature profile. This direction is in very close alignment with the afore measured anomalously large bulk flow direction. Using Monte Carlo simulations, we estimate the significance of the giant rings at the $3\sigma$ level and the alignment with the bulk flow at $2.5\sigma$. We argue that a cosmic defect seeded by a pre-inflationary particle could explain the giant rings, the large bulk flow and their alignment.

We study some of the cosmological imprints of pre-inflationary particles. We show that each such particle provides a seed for a spherically symmetric cosmic defect. The profile of this cosmic defect is fixed and its magnitude is linear in a single parameter that is determined by the mass of the pre-inflationary particle. We study the CMB and peculiar velocity imprints of this cosmic defect and suggest that it could explain some of the large scale cosmological anomalies.

We show that models of inflection point inflation exhibit a phase transition from a region in parameter space where they are of large field type to a region where they are of small field type. The phase transition is between a universal behavior, with respect to the initial condition, at the large field region and non-universal behavior at the small field region. The order parameter is the number of e-foldings. We find integer critical exponents at the transition between the two phases.

We consider models of inflection point inflation. The main drawback of such models is that they suffer from the overshoot problem. Namely the initial condition should be fine tuned to be near the inflection point for the universe to inflate. We show that stringy realizations of inflection point inflation are common and offer a natural resolution to the overshoot problem.