Search

Primordial black holes (PBHs) would induce non-Gaussianity in the cosmic microwave background (CMB) by sourcing recombination perturbations spatially modulated by relative velocities between PBHs and the baryons they accrete. The leading non-Gaussian signatures are non-vanishing connected 4-point correlation functions, or trispectra. Earlier, we computed the CMB temperature trispectrum, and forecasted Planck to be more sensitive to it than to changes in the CMB temperature power spectrum for light enough PBHs. Excitingly, accreting PBHs would also induce non-Gaussianity in CMB polarization, and source both E and B modes, which we compute in this paper. We first calculate linear-response perturbations to the tensor-valued photon distribution function sourced by a general spatially-varying ionization history, and apply our results to accreting PBHs. We then compute linear-order perturbations to the temperature and polarization 2-point functions sourced by inhomogeneities in recombination due to accreting PBHs; we find them to be negligible relative to their counterparts sourced by homogeneous perturbations to the ionization history. Lastly, we compute all CMB trispectra including temperature, E- and B-mode polarization at linear order in the PBH abundance. We forecast that including polarization data in a 4-point-function analysis would only increase Planck's sensitivity to accreting PBHs by a factor ~2 relative to using temperature alone. As a consequence, we find that a search for PBHs using all temperature and polarization trispectra with Planck data would mostly not be competitive with current bounds from temperature and polarization power spectra. In contrast, we forecast that a CMB Stage-4 experiment would gain significant sensitivity to accreting PBHs through a 4-point-function search, in particular through the contributions of parity-odd trispectra including one B-mode field., Comment: 16 pages + appendices and refs, 3 figures

Weak magnetic fields must have existed in the early Universe, as they were sourced by the cross product of electron density and temperature gradients through the Biermann-battery mechanism. In this paper we calculate the magnetic fields generated at cosmic dawn by a variety of small-scale primordial perturbations, carefully computing the evolution of electron density and temperature fluctuations, and consistently accounting for relative velocities between baryons and dark matter. We first compute the magnetic field resulting from standard, nearly scale-invariant primordial adiabatic perturbations, making significant improvements to previous calculations. This "standard" primordial field has a root mean square (rms) of $\sim10^{-15}$ nG at $20\lesssim z \lesssim 100$, with fluctuations on $\sim$ kpc comoving scales, and could serve as the seed of present-day magnetic fields observed in galaxies and galaxy clusters. In addition, we consider early-Universe magnetic fields as a possible probe of non-standard initial conditions of the Universe on small scales $k \sim 1-10^3$ Mpc$^{-1}$. To this end, we compute the maximally-allowed magnetic fields within current upper limits on small-scale adiabatic and isocurvature perturbations. Under the current Cosmic Microwave Background spectral-distortion constraints magnetic fields could be produced with a rms of $\sim 5\times 10^{-11}$ nG at $z = 20$. Uncorrelated small-scale isocurvature perturbations within current Big-Bang Nucleosynthesis bounds could potentially enhance the magnetic field to $\sim 10^{-14}-10^{-10}$ nG at $z = 20$, depending on the specific isocurvature mode considered. While these very weak fields remain well below current observational capabilities, our work points out that magnetic fields could potentially provide an interesting window into the poorly constrained small-scale initial conditions of the Universe., Comment: 12 pages, 6 figures

Primordial magnetic fields (PMFs) may explain observations of magnetic fields on extragalactic scales. They are most cleanly constrained by measurements of cosmic microwave background radiation (CMB) anisotropies. Their effects on cosmic recombination may even be at the heart of the resolution of the Hubble tension. We present the most detailed analysis of the effects of PMFs on cosmic recombination to date. To this end we extend the public magneto-hydrodynamic code {\sl ENZO} with a new cosmic recombination routine, Monte-Carlo simulations of Lyman-$\alpha$ photon transport, and a Compton drag term in the baryon momentum equation. The resulting code allows us, for the first time, to realistically predict the impact of PMFs on the cosmic ionization history and the clumping of baryons during cosmic recombination. Our results identify the importance of mixing of Lyman-$\alpha$ photons between overdense- and underdense- regions for small PMF strength. This mixing speeds up recombination beyond the speed-up due to clumping. For non-helical PMFs with a Batchelor spectrum we find a surprising dependency of results on ultra-violet magnetic modes. The present study shall serve as a theoretical foundation for a future precise comparison of recombination with PMFs to CMB data., Comment: 22 pages, 14 figures

Elastic scattering of dark matter (DM) particles with baryons induce cosmological signals that may be detectable with modern or future telescopes. For DM-baryon scattering cross sections scaling with negative powers of relative velocity, $\sigma_{\chi b}(v) \propto v^{-2}, v^{-4}$, such interactions introduce a momentum-exchange rate that is nonlinear in DM-baryon bulk relative velocities, thus not amenable for inclusion as-is into standard linear cosmological Boltzmann codes. Linear ansatzes have been adopted in past works, but their accuracy is unknown as they do not arise from first-principles derivations. In this work, for the first time, we construct a rigorous framework for computing linear-cosmology observables as a perturbative expansion in $\sigma_{\chi b}$. We argue that this approach is accurate for Cosmic Microwave Background (CMB) angular power spectra when most or all of the DM is scattering with baryons with cross section $\sigma_{\chi b}(v) \propto v^{-2}, v^{-4}$. We derive exact formal expressions for CMB power spectra at linear order in $\sigma_{\chi b}$, and show that they only depend on a specific velocity integral of the momentum-exchange rate. Consequently, we can obtain the exact power spectra at linear order in $\sigma_{\chi b}$ by substituting the original nonlinear momentum-exchange rate with a uniquely specified linear rate. Serendipitously, we find that the exact substitution we derive from first principles precisely coincides with the most widely used linear ansatz, thus placing previous CMB-anisotropy upper bounds on a more solid footing. In addition to finally providing an exact cosmological solution to the DM-baryon scattering problem in a well-defined region of parameter space, the framework we construct opens the way to computing higher-order correlation functions, beyond power spectra, which are promising yet unexplored probes of DM-baryon scattering., Comment: 10 pages, 2 figures

We construct data-driven solutions to the Hubble tension which are perturbative modifications to the fiducial $\Lambda$CDM cosmology, using the Fisher bias formalism. Taking as proof of principle the case of a time-varying electron mass and fine structure constant, and focusing first on Planck CMB data, we demonstrate that a modified recombination can solve the Hubble tension and lower $S_8$ to match weak lensing measurements. Once baryonic acoustic oscillation and uncalibrated supernovae data are included, however, it is not possible to fully solve the tension with perturbative modifications to recombination., Comment: 5+13 pages, 2+12 figures, version accepted for publication in PRL

It is known that Primordial Black Holes (PBHs) can leave an imprint on Cosmic Microwave Background (CMB) anisotropy power spectra, due to their accretion-powered injection of energy into the recombining plasma. Here we study a qualitatively new CMB observable sourced by accreting PBHs: the temperature trispectrum or connected 4-point function. This non-Gaussian signature is due to the strong spatial modulation of the PBH accretion luminosity, thus ionization perturbations, by large-scale supersonic relative velocities between PBHs and the accreted baryons. We first derive a factorizable quadratic transfer function for free-electron fraction inhomogeneities induced by accreting PBHs. We then compute the perturbation to the CMB temperature anisotropy due to a general modification of recombination, and apply our results to accreting PBHs. We calculate a new contribution to the temperature power spectrum due to the spatial fluctuations of the ionization perturbation induced by accreting PBHs, going beyond past studies which only accounted for its homogeneous part. While these contributions are formally comparable, we find the new part to be subdominant, due to the poor correlation of the perturbed temperature field with the standard CMB anisotropy. For the first time, we compute the temperature trispectrum due to accreting PBHs. This trispectrum is weakly correlated with the local-type primordial non-Gaussianity trispectrum, hence constraints on the latter do not lead to competitive bounds on accreting PBHs. We also forecast Planck's sensitivity to the temperature trispectrum sourced by accreting PBHs. Excitingly, we find it to be more sensitive to PBHs under $\sim 10^3 M_{\odot}$ than current temperature-only power spectrum constraints. This result motivates our future work extending this study to temperature and polarization trispectra induced by inhomogeneously-accreting PBHs.

Linear cosmological observables can be used to probe elastic scattering of dark matter (DM) with baryons. Availability of high-precision data requires a critical reassessment of any assumptions that may impact the accuracy of constraints. The standard formalism for constraining DM-baryon scattering pre-recombination is based on assuming a Maxwell-Boltzmann (MB) velocity distribution for DM. This assumption is not always justified and does not allow for probing DM self-interactions in addition to its interactions with baryons. Lifting the MB assumption requires solving the full collisional Boltzmann equation (CBE), which is highly non-trivial. Earlier work proposed a more tractable Fokker-Planck (FP) approximation to the CBE, but its accuracy remained unknown. In this work, we numerically solve the exact CBE for the first time, in a homogeneous expanding background. We consider DM-baryon scattering cross-sections that are positive power-laws of relative velocity. We derive analytical expressions for the collision operator in the case of isotropic differential scattering cross-sections. We then solve the background CBE numerically and use our solution for the DM velocity distribution to compute the DM-baryon heat-exchange rate, which we compare against those obtained with the MB assumption and FP approximation. Over a broad range of DM-to-baryon mass ratios, we find that the FP approximation leads to a maximum error of 17%, significantly better than the up to 160% error introduced by the MB assumption. While our results strictly apply only to the background evolution, the accuracy of the FP approximation is likely to carry over to perturbations. This motivates its implementation into cosmological Boltzmann codes, where it can supersede the much less accurate MB assumption, and allow for a more general exploration of DM interactions with baryons and with itself., Comment: Accepted for publication in the Physical Review D; 19 pages, 4 figures

Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivity to PBHs. However, devoted analysis pipelines and theoretical modeling are required to fully leverage these novel data. The search for PBHs constitutes a low-cost, high-reward science case with significant impact on the high energy physics community., Comment: 22 pages, 4 figures, submission to the Snowmass 2021 process. Update to add references

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade., Comment: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

The Universe's initial conditions, in particular baryon and cold dark matter (CDM) isocurvature perturbations, are poorly constrained on sub-Mpc scales. In this paper, we develop a new formalism to compute the effect of small-scale baryon perturbations on the mean free-electron abundance, thus on cosmic microwave background (CMB) anisotropies. Our framework can accommodate perturbations with arbitrary time and scale dependence. We apply this formalism to four different combinations of baryon and CDM isocurvature modes, and use Planck CMB-anisotropy data to probe their initial amplitude. We find that Planck data is consistent with no small-scale isocurvature perturbations, and that this additional ingredient does not help alleviate the Hubble tension. We set upper bounds to the dimensionless initial power spectrum $\Delta_{\mathcal{I}}^2(k)$ of these isocurvature modes at comoving wavenumbers $1~\textrm{Mpc}^{-1} \le k \le 10^3$ Mpc$^{-1}$, for several parameterizations. For a scale-invariant power spectrum, our 95% confidence-level limits on $\Delta_{\mathcal{I}}^2$ are 0.023 for pure baryon isocurvature, 0.099 for pure CDM isocurvature, 0.026 for compensated baryon-CDM perturbations, and 0.009 for joint baryon-CDM isocurvature perturbations. Using a Fisher analysis generalized to non-analytic parameter dependence, we forecast that a CMB Stage-4 experiment would be able to probe small-scale isocurvature perturbations with initial power 3 to 10 times smaller than Planck limits. The formalism introduced in this work is very general and can be used more widely to probe any physical processes or initial conditions sourcing small-scale baryon perturbations., Comment: 17 pages, 7 figures. Results are available at https://github.com/nanoomlee/small-scale_baryon_CDM_isocurvature_results

Exotic electromagnetic energy injection in the early Universe may alter cosmological recombination, and ultimately cosmic microwave background (CMB) anisotropies. Moreover, if energy injection is inhomogeneous, it may induce a spatially-varying ionization fraction, and non-Gaussianity in the CMB. The observability of these signals, however, is contingent upon how far the injected particles propagate and deposit their energy into the primordial plasma, relative to the characteristic scale of energy injection fluctuations. In this study we inspect the spatial properties of energy deposition and perturbed recombination resulting from an inhomogeneous energy injection of sub-10 MeV photons, relevant to accreting primordial black holes (PBHs). We develop a novel Monte-Carlo radiation transport code accounting for all relevant photon interactions in this energy range, and including secondary electron energy deposition efficiency through a new analytic approximation. For a specified injected photon spectrum, the code outputs an injection-to-deposition Green's function depending on time and distance from the injection point. Combining this output with a linearized solution of the perturbed recombination problem, we derive time- and scale-dependent deposition-to-ionization Green's functions. We apply this general framework to accreting PBHs, whose luminosity is strongly spatially modulated by supersonic relative velocities between cold dark matter and baryons. We find that the resulting spatial fluctuations of the free-electron fraction are of the same magnitude as its mean deviation from standard recombination, from which current CMB power spectra constraints are derived. This work suggests that the sensitivity to accreting PBHs might be substantially improved by propagating these inhomogeneities to CMB anisotropy power spectra and non-Gaussian statistics, which we study in subsequent papers., Comment: Version accepted in PRD

In recent years, weak lensing of the cosmic microwave background (CMB) has emerged as a powerful tool to probe fundamental physics, such as neutrino masses, primordial non-Gaussianity, dark energy, and modified gravity. The prime target of CMB lensing surveys is the lensing potential, which is reconstructed from observed CMB temperature $T$ and polarization $E$ and $B$ fields. Until very recently, this reconstruction has been performed with quadratic estimators (QEs), which, although known to be suboptimal for high-sensitivity experiments, are numerically efficient, and useful to make forecasts and cross-check the results of more sophisticated likelihood-based methods. It is expected that ongoing and near-future CMB experiments such as AdvACT, SPT-3G and the Simons Observatory (SO), will also rely on QEs. Here, we review different QEs, and clarify their differences. In particular, we show that the Hu-Okamoto (HO02) estimator is not the absolute optimal lensing estimator that can be constructed out of quadratic combinations of $T, E$ and $B$ fields. Instead, we derive the global-minimum-variance (GMV) lensing quadratic estimator. Although this estimator can be found elsewhere in the literature, it was erroneously described as equivalent to the HO02 estimator, and has never been used in real data analyses. Here, we show explicitly that the HO02 estimator is suboptimal to the GMV estimator, with a reconstruction noise larger by up to $\sim 9\%$ for a SO-like experiment. We further show that the QE used in the Planck, and recent SPT lensing analysis are suboptimal to both the HO02 and GMV estimator, and would have a reconstruction noise up to $\sim 11\%$ larger than that of the GMV estimator for a SO-like experiment. In addition to clarifying differences between different QEs, this work should thus provide motivation to implement the GMV estimator in future lensing analyses relying on QEs., Comment: Accepted for publication in PRD

Possible interactions of dark matter (DM) with Standard Model (SM) particles can be tested with spectral distortions (SDs) of the cosmic microwave background (CMB). In particular, a non-relativistic DM particle that scatters elastically with photons, electrons or nuclei imprints a negative chemical potential $\mu$ to the CMB spectrum. This article revisits the first study of this effect, with an accurate treatment of heat exchange between DM and SM particles. We show that the instantaneous-decoupling approximation made in the original study systematically and significantly underestimates the amplitude of SDs. As a consequence, we derive tighter upper bounds to the DM-SM elastic-scattering cross section for DM masses $m_\chi \lesssim 0.1$ MeV, from the non-detection of $\mu$-distortions by FIRAS. We also show that a future instrument like PIXIE, sensitive to $|\mu| \sim 10^{-8}$, would be able to probe DM-SM cross sections much smaller than first forecasted, and orders of magnitude below current upper limits from CMB-anisotropy data, up to DM masses of $\sim 1$ GeV. Lastly, we study the sensitivity of SDs to the electric and magnetic dipole moments of the DM. Although SDs can place non-trivial constraints on these models, we find that even future SD experiments are unlikely to improve upon the best current bounds. This article is accompanied by the public code DMDIST, which allows one to compute CMB SDs for generic particle-DM models, specified by their cross sections for elastic scattering with and annihilation into SM particles., Comment: 13 pages + refs, 9 figures. DMDIST available at https://cosmo.nyu.edu/yacine/dmdist/dmdist.html

Within the next several years pulsar timing arrays (PTAs) are positioned to detect the stochastic gravitational-wave background (GWB) likely produced by the collection of inspiralling super-massive black holes binaries, and potentially constrain some exotic physics. So far most of the pulsar timing data analysis has focused on the monopole of the GWB, assuming it is perfectly isotropic. The natural next step is to search for anisotropies in the GWB. In this paper, we use the recently developed PTA Fisher matrix to gain insights into optimal search strategies for GWB anisotropies. For concreteness, we apply our results to EPTA data, using realistic noise characteristics of its pulsars. We project the detectability of a GWB whose angular dependence is assumed to be a linear combination of predetermined maps, such as spherical harmonics or coarse pixels. We find that the GWB monopole is always statistically correlated with these maps, implying a loss of sensitivity to the monopole when searching simultaneously for anisotropies. We then derive the angular distributions of the GWB intensity to which a PTA is most sensitive, and illustrate how one may use these "principal maps" to approximately reconstruct the angular dependence of the GWB. Since the principal maps are neither perfectly anisotropic nor uncorrelated with the monopole, we also develop a frequentist criterion to specifically search for anisotropies in the GWB without any prior knowledge about their angular distribution. Lastly, we show how to recover existing EPTA results with our Fisher formalism, and clarify their meaning. The tools presented here will be valuable in guiding and optimizing the computationally demanding analyses of pulsar timing data., Comment: 21 pages, 14 figures. Version accepted for publication in PRD after minor changes

The current cosmological probes have provided a fantastic confirmation of the standard $\Lambda$ Cold Dark Matter cosmological model, that has been constrained with unprecedented accuracy. However, with the increase of the experimental sensitivity a few statistically significant tensions between different independent cosmological datasets emerged. While these tensions can be in portion the result of systematic errors, the persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the need for new physics. In this Letter of Interest we will focus on the $4.4\sigma$ tension between the Planck estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we will list a few interesting new physics models that could solve this tension and discuss how the next decade experiments will be crucial., Comment: Snowmass2021 - Letter of Interest

A precise measurement of the curvature of the Universe is of primeval importance for cosmology since it could not only confirm the paradigm of primordial inflation but also help in discriminating between different early Universe scenarios. The recent observations, while broadly consistent with a spatially flat standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) model, are showing tensions that still allow (and, in some cases, even suggest) a few percent deviations from a flat universe. In particular, the Planck Cosmic Microwave Background power spectra, assuming the nominal likelihood, prefer a closed universe at more than 99\% confidence level. While new physics could be in action, this anomaly may be the result of an unresolved systematic error or just a statistical fluctuation. However, since a positive curvature allows a larger age of the Universe, an accurate determination of the age of the oldest objects provides a smoking gun in confirming or falsifying the current flat $\Lambda$CDM model., Comment: Snowmass2021 - Letter of Interest

The standard $\Lambda$ Cold Dark Matter cosmological model provides a wonderful fit to current cosmological data, but a few tensions and anomalies became statistically significant with the latest data analyses. While these anomalies could be due to the presence of systematic errors in the experiments, they could also indicate the need for new physics beyond the standard model. In this Letter of Interest we focus on the tension of the Planck data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models for solving this tension, and we discuss the importance of trying to fit with a single model a full array of data and not just one parameter at a time., Comment: Snowmass2021 - Letter of Interest

The standard $\Lambda$ Cold Dark Matter cosmological model provides an amazing description of a wide range of astrophysical and astronomical data. However, there are a few big open questions, that make the standard model look like a first-order approximation to a more realistic scenario that still needs to be fully understood. In this Letter of Interest we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances present between the different cosmological probes, as the Hubble constant $H_0$ value, the $\sigma_8 - S_8$ tension, and the anomalies present in the Planck results. Finally, we will give an overview of upgraded experiments and next-generation space-missions and facilities on Earth, that will be of crucial importance to address all these questions., Comment: Snowmass2021 - Letter of Interest

We present the new recombination code HYREC-2, holding the same accuracy as the state-of-the-art codes HYREC and COSMOREC and, at the same time, surpassing the computational speed of the code RECFAST commonly used for CMB-anisotropy data analyses. HYREC-2 is based on an effective 4-level atom model, accounting exactly for the non-equilibrium of highly excited states of hydrogen, and very accurately for radiative transfer effects with a correction to the Lyman-$\alpha$ escape rate. The latter is computed with the original HYREC, and tabulated, as a function of temperature, along with its derivatives with respect to the relevant cosmological parameters. This enables the code to keep the same accuracy as the original HYREC-2 over the full 99.7% confidence region of cosmological parameters currently allowed by Planck, while running in under one millisecond on a standard laptop. Our code leads to no noticeable bias in any cosmological parameters even in the case of an ideal cosmic-variance limited experiment up to $\ell$ = 5000. Beyond CMB anisotropy calculations, HYREC-2 will be a useful tool to compute various observables that depend on the recombination and thermal history, such as the recombination spectrum or the 21-cm signal., Comment: 14pages, 7 Figures, version accepted for publication in PRD

Pulsar timing arrays (PTAs) are currently the only experiments directly sensitive to gravitational waves with decade-long periods. Within the next five to ten years, PTAs are expected to detect the stochastic gravitational-wave background (SGWB) collectively sourced by inspiralling supermassive black hole binaries. It is expected that this background is mostly isotropic, and current searches focus on the monopole part of the SGWB. Looking ahead, anisotropies in the SGWB may provide a trove of additional information both on known and unknown astrophysical and cosmological sources. In this paper, we build a simple yet realistic Fisher formalism for anisotropic SGWB searches with PTAs. Our formalism is able to accommodate realistic properties of PTAs, and allows simple and accurate forecasts. We illustrate our approach with an idealized PTA consisting of identical, isotropically distributed pulsars. In a companion paper, we apply our formalism to current PTAs and show that it can be a powerful tool to guide and optimize real data analysis., Comment: Version accepted for publication in PRD after minor changes. Follow-up paper on applications: arXiv:2010.13958

We study if the discrepancy between the local and cosmological measurements of the Hubble constant $H_0$ can be reformulated as a tension in the cosmic microwave background (CMB) monopole temperature $T_0$. The latter is customarily fixed to the FIRAS best-fit value in CMB data analyses. Although this value was confirmed by several independent experiments, it is interesting to see how much parameter constraints depend on this prior. We first provide a detailed pedagogical description of the $T_0$ effects on cosmological observables. We show that the recombination history and transfer functions do not depend on $T_0$, provided they are parametrized by the energy scale rather than redshift, and at fixed dark matter and baryon densities per CMB photon. Thus, $T_0$ is only a property of the observer, quantifying the amount of expansion between key cosmological events and today. As a consequence, the sole effect of $T_0$ on small-scale primary CMB anisotropies is through the angular diameter distance to the epoch of last scattering, resulting in a near-perfect degeneracy between $T_0$ and $H_0$. This geometric degeneracy is partially lifted by the late-time integrated Sachs-Wolfe effect and CMB lensing. Still, Planck data alone is consistent with a broad region in the $H_0-T_0$ plane, implying that removing the FIRAS prior on $T_0$ can make Planck and SH0ES less discrepant, without introducing new physics beyond $\Lambda$CDM. One may break the degeneracy by combining Planck with SH0ES, yielding an independent measurement of $T_0$, which happens to be in $3\sigma$ tension with FIRAS. Therefore, the Hubble tension indeed can be recast into the $T_0$ tension. The agreement with FIRAS is restored when combining Planck with the baryon acoustic oscillation data instead of SH0ES. Thus, the tension between SH0ES and cosmological measurements of $H_0$ persists even if we discard the FIRAS $T_0$ prior., Comment: v2: 22 pages, 8 figures, corrected a mistake due to calibration of temperature maps, qualitative results and conclusions unchanged

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., Comment: 20 pages, white paper submitted in answer to the "Voyage 2050" call to prepare the long term plan in the ESA science programme

Cold dark matter (CDM) could be composed of primordial black holes (PBH) in addition to or instead of more orthodox weakly interacting massive particle dark matter (PDM). We study the formation of the first structures in such $\Lambda$PBH cosmologies using $N$-body simulations evolved from deep in the radiation era to redshift 99. When PBH are only a small component of the CDM, they are clothed by PDM to form isolated halos. On the other hand, when PBH make most of the CDM, halos can also grow via clustering of many PBH. We find that the halo mass function is well modelled via Poisson statistics assuming random initial conditions. We quantify the nonlinear velocities induced by structure formation and find that they are too small to significantly impact CMB constraints. A chief challenge is how best to extrapolate our results to lower redshifts relevant for some observational constraints., Comment: 19 pages, 16 figures

We examine the effects ultra-light primordial black holes (PBHs) have on the anisotropies of the cosmic microwave background (CMB). PBHs in the mass range of $10^{15}$ to $10^{17}$ g emit Hawking radiation in the early Universe, modifying the standard recombination history. This leads to a damping of small-scale temperature and polarisation anisotropies and enhances large-scale polarisation fluctuations. As some models of inflation predict PBHs with a range of masses, we investigate the impacts of extended mass distributions on PBH abundance constraints. We model PBH energy injection using a ground-up approach incorporating species-dependent deposition efficiencies. By allowing the $\Lambda$CDM parameters to vary simultaneously with the PBH fraction and mass, we show that exclusion bounds on the PBH fraction of DM $f_\text{PBH}$ are relaxed by up to an order of magnitude, compared to the case of fixed $\Lambda$CDM parameters. We also give 95% exclusion regions for $f_\text{PBH}$ for a variety of mass distributions. In particular, for a uniform mass distribution between $10^{15}$ and $10^{17}$ g, we find $f_\text{PBH} < 1.6 \times 10^{-5}$ when allowing $\Lambda$CDM parameters to vary.

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., Comment: Science white paper submitted to Decadal 2020 survey

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., Comment: White paper submitted to the Astro2020 Decadal Survey

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., Comment: 10 pages, 4 figures

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., Comment: 8 pages, 4 figures

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., Comment: 9 pages, 3 figures

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal)., Comment: Astro2020 Decadal Survey science white paper; list of endorsers available at https://github.com/veragluscevic/Astro2020-DM-Cosmology-Endorsers

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., Comment: 5 pages, 1 figure, submitted to the Astro2020 Decadal Survey Science White Paper call

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., Comment: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

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., Comment: Science White paper submitted to Astro2020 Decadal Survey

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ ., Comment: 96 pages, 22 figures, 1 table

Linear-cosmology observables, such as the Cosmic Microwave Background (CMB), or the large-scale distribution of matter, have long been used as clean probes of dark matter (DM) interactions with baryons. It is standard to model the DM as an ideal fluid with a thermal Maxwell-Boltzmann (MB) velocity distribution, in order to compute the heat and momentum-exchange rates relevant to these probes. This approximation only applies in the limit where DM self-interactions are frequent enough to efficiently redistribute DM velocities. It does not accurately describe weakly self-interacting particles, whose velocity distribution unavoidably departs from MB once they decouple from baryons. This article lays out a new formalism required to accurately model DM-baryon scattering, even when DM self-interactions are negligible. The ideal fluid equations are replaced by the collisional Boltzmann equation for the DM phase-space distribution. The collision operator is approximated by a Fokker-Planck operator, constructed to recover the exact heat and momentum exchange rates, and allowing for an efficient numerical implementation. Numerical solutions to the background evolution are presented, which show that the MB approximation can over-estimate the heat-exchange rate by factors of ~ 2-3, especially for light DM particles. A Boltzmann-Fokker-Planck hierarchy for perturbations is derived. This new formalism allows to explore a wider range of DM models, and will be especially relevant for upcoming ultra-high-sensitivity CMB probes., Comment: 18 pages, 3 figures. Version accepted for publication in PRD after minor additions

Primordial black holes (PBHs) have been brought back into the spotlight by LIGO's first direct detection of a binary-black-hole merger. One of the poorly understood properties of PBHs is how clustered they are at formation. It has important implications on the efficacy of their merging in the early Universe, as well as on observational constraints. In this work we study the initial clustering of PBHs formed from the gravitational collapse of large density fluctuations in the early Universe. We give a simple and general argument showing that, in this scenario, we do not expect clustering on very small scales beyond what is expected from a random, Poisson distribution. We illustrate this result explicitly in the case where the underlying density field is Gaussian. We moreover derive a new analytic expression for the two-point correlation function of large-threshold fluctuations, generalizing previous results to arbitrary separation, and with broader implications than the clustering of PBHs., Comment: 5 pages, 2 figures. A few clarifications, main results unchanged. Version accepted for publication in PRL

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., Comment: 14 pages, 9 figures, submitted to PRD

We present an efficient and accurate method for simulating massive neutrinos in cosmological structure formation simulations, together with an easy to use public implementation. Our method builds on our earlier implementation of the linear response approximation (LRA) for neutrinos, coupled with an N-body code for cold dark matter particles. The LRA's good behaviour at early times and in the linear regime is preserved, while better following the non-linear clustering of neutrinos on small scales. Massive neutrinos are split into initially "fast" and "slow" components. The fast component is followed analytically with the LRA all the way to redshift zero. The slow component is evolved with the LRA only down to a switch-on redshift $z_\nu = 1$, below which it is followed with the particle method, in order to fully account for its non-linear evolution. The slow neutrino particles are initialized at $z = 99$ in order to have accurate positions and velocities at the switch-on time, but are not used to compute the potential until $z \leq 1$, thus avoiding the worst effect of particle shot noise. We show that our hybrid method matches (and for small neutrino masses, exceeds) the accuracy of neutrino particle simulations with substantially lower particle load requirements., Comment: 16 pages, 13 figures. Version accepted by MNRAS

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., Comment: Version accepted for publication in PRD after minor changes

Interest in the idea that primordial black holes (PBHs) might comprise some or all of the dark matter has recently been rekindled following LIGO's first direct detection of a binary-black-hole merger. Here we revisit the effect of accreting PBHs on the cosmic microwave background (CMB) frequency spectrum and angular temperature/polarization power spectra. We compute the accretion rate and luminosity of PBHs, accounting for their suppression by Compton drag and Compton cooling by CMB photons. We estimate the gas temperature near the Schwarzschild radius, and hence the free-free luminosity, accounting for the cooling resulting from collisional ionization when the background gas is mostly neutral. We account approximately for the velocities of PBHs with respect to the background gas. We provide a simple analytic estimate of the efficiency of energy deposition in the plasma. We find that the spectral distortions generated by accreting PBHs are too small to be detected by FIRAS, as well as by future experiments now being considered. We analyze Planck CMB temperature and polarization data and find, under our most conservative hypotheses, and at the order-of-magnitude level, that they rule out PBHs with masses >~ 10^2 M_sun as the dominant component of dark matter., Comment: Version accepted in PRD. Main changes: (i) new section on local radiative feedback. (ii) replaced the "strong feedback" limit by a "photoionization" limit. Massive PBHs have a higher luminosity and tighter constraints in the latter than in the former. The most conservative bounds (derived in the "collisional ionization" limit, equivalent to the "no feedback" case in v1) remain unchanged

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., Comment: 13 pages, 8 figures, 3 tables and 1 appendix

We study the degree to which the cosmic microwave background (CMB) can be used to constrain primordial non-Gaussianity involving one tensor and two scalar fluctuations, focusing on the correlation of one polarization $B$ mode with two temperature modes. In the simplest models of inflation, the tensor-scalar-scalar primordial bispectrum is non-vanishing and is of the same order in slow-roll parameters as the scalar-scalar-scalar bispectrum. We calculate the $\langle BTT\rangle$ correlation arising from a primordial tensor-scalar-scalar bispectrum, and show that constraints from an experiment like CMB-Stage IV using this observable are more than an order of magnitude better than those on the same primordial coupling obtained from temperature measurements alone. We argue that $B$-mode non-Gaussianity opens up an as-yet-unexplored window into the early Universe, demonstrating that significant information on primordial physics remains to be harvested from CMB anisotropies., Comment: 14 pages and 8 figures. Version prepared for submission to PRD. New: paragraph pointing out analogy with gravitational lensing, paragraph on squeezed limit, paragraph on normalization of fNL, new figure (7), updated figures based on improved sampling in \ell, typos, references