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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.

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.

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.

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

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.

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.

A measurement of primordial non-gaussianity will be of paramount importance to distinguish between different models of inflation. Cosmic microwave background (CMB) anisotropy observations have set unprecedented bounds on the non-gaussianity parameter f_NL but the interesting regime f_NL <~ 1 is beyond their reach. Brightness-temperature fluctuations in the 21-cm line during the dark ages (z ~ 30-100) are a promising successor to CMB studies, giving access to a much larger number of modes. They are, however, intrinsically non-linear, which results in secondary non-gaussianities orders of magnitude larger than the sought-after primordial signal. In this paper we carefully compute the primary and secondary bispectra of 21-cm fluctuations on small scales. We use the flat-sky formalism, which greatly simplifies the analysis, while still being very accurate on small angular scales. We show that the secondary bispectrum is highly degenerate with the primordial one, and argue that even percent-level uncertainties in the amplitude of the former lead to a bias of order Delta f_NL ~ 10. To tackle this problem we carry out a detailed Fisher analysis, marginalizing over the amplitudes of a few smooth redshift-dependent coefficients characterizing the secondary bispectrum. We find that the signal-to-noise ratio for a single redshift slice is reduced by a factor of ~5 in comparison to a case without secondary non-gaussianities. Setting aside foreground contamination, we forecast that a cosmic-variance-limited experiment observing 21-cm fluctuations over 30 < z < 100 with a 0.1-MHz bandwidth and 0.1-arcminute angular resolution could achieve a sensitivity of order f_NL[local] ~ 0.03, f_NL[equilateral] ~ 0.04, and f_NL[orthogonal] ~ 0.03.

We propose a new method to constrain elastic scattering between dark matter (DM) and standard model particles in the early Universe. Direct or indirect thermal coupling of non-relativistic DM with photons leads to a heat sink for the latter. This results in spectral distortions of the cosmic microwave background (CMB), the amplitude of which can be as large as a few times the DM-to-photon number ratio. We compute CMB spectral distortions due to DM-proton, DM-electron and DM-photon scattering for generic energy-dependent cross sections and DM mass m_DM >~ 1 keV. Using FIRAS measurements we set constraints on the cross sections for m_DM <~ 0.1 MeV. In particular, for energy-independent scattering we obtain sigma[DM-proton] <~ 10^(-24) cm^2 (keV/m_DM)^(1/2), sigma[DM-electron] <~ 10^(-27) cm^2 (keV/m_DM)^(1/2) and sigma[DM-photon] <~ 10^(-39) cm^2 (m_DM/keV). An experiment with the characteristics of PIXIE would extend the regime of sensitivity up to masses m_DM ~ 1 GeV.

A new approach to cosmological perturbation theory has been recently introduced by Bartelmann et al., relying on nonequilibrium statistical theory of classical particles, and treating the gravitational interaction perturbatively. They compute analytic expressions for the nonlinear matter power spectrum, to first order in the interaction, and at one-loop order in the linear power spectrum. The resulting power spectrum is well behaved even at large wavenumbers and seems in good agreement with results from numerical simulations. In this paper, we rederive their results concisely with a different approach, starting from the implicit integral solution to particle trajectories. We derive the matter power spectrum to first order in the interaction, but to arbitrary order in the linear power spectrum, from which the one-loop result follows. We also show that standard linear perturbation theory can only be recovered at infinite order in the gravitational interaction. At finite order in the interaction, we find that the linear power spectrum is systematically and significantly underestimated. A comprehensive study of the convergence of the theory with the order of the interaction for nonlinear scales will be the subject of future work.

Polycyclic Aromatic Hydrocarbons (PAHs) are believed to be the small-size tail of the interstellar carbonaceous dust grain population. Their vibrational emission is the most widely accepted source of the aromatic near-infrared features, and their rotational radiation is a likely explanation for the dust-correlated anomalous microwave emission (AME). Yet, no individual interstellar PAH molecule has been identified to date. It was recently recognized that quasi-symmetric planar PAHs ought to have a well identifiable comb-like rotational spectrum, and suggested to search for them in spectroscopic data with matched-filtering techniques. We report the results of the first such search, carried out with the Green Bank Telescope, and targeting the star-forming region IC348 in the Perseus molecular cloud, a known source of AME. Our observations amounted to 16.75 hours and spanned a 3 GHz-wide band extending from 23.3 to 26.3 GHz. Using frequency switching, we achieved a sensitivity of 0.4 mJy per 0.4 MHz channel (1-sigma). The non-detection of comb-like spectra allowed us to set upper bounds on the abundance of specific quasi-symmetric PAH molecules (specified uniquely by their moments of inertia) of approximately 0.1% of the total PAH abundance. This bound generically applies to PAHs with approximately 15 to 100 carbon atoms.

Fluctuations of the 21 cm brightness temperature before the formation of the first stars hold the promise of becoming a high-precision cosmological probe in the future. The growth of over densities is very well described by perturbation theory at that epoch and the signal can in principle be predicted to arbitrary accuracy for given cosmological parameters. Recently, Tseliakhovich and Hirata pointed out a previously neglected and important physical effect, due to the fact that baryons and cold dark matter (CDM) have supersonic relative velocities after recombination. This relative velocity suppresses the growth of matter fluctuations on scales $k \sim 10-10^3$ Mpc$^{-1}$. In addition, the amplitude of the small-scale power spectrum is modulated on the large scales over which the relative velocity varies, corresponding to $k \sim 0.005-1$ Mpc$^{-1}$. In this paper, the effect of the relative velocity on 21 cm brightness temperature fluctuations from redshifts $z \geq 30$ is computed. We show that the 21 cm power spectrum is affected on ${\it most}$ scales. On small scales, the signal is typically suppressed several tens of percent, except for extremely small scales ($k \gtrsim 2000$ Mpc$^{-1}$) for which the fluctuations are boosted by resonant excitation of acoustic waves. On large scales, 21 cm fluctuations are enhanced due to the non-linear dependence of the brightness temperature on the underlying gas density and temperature. The enhancement of the 21 cm power spectrum is of a few percent at $k \sim 0.1$ Mpc$^{-1}$ and up to tens of percent at $k \lesssim 0.005$ Mpc$^{-1}$, for standard $\Lambda$CDM cosmology. In principle this effect allows to probe the small-scale matter power spectrum not only through a measurement of small angular scales but also through its effect on large angular scales.

Polycyclic aromatic hydrocarbons (PAHs) are believed to be ubiquitous in the interstellar medium. Yet, to date no specific PAH molecule has been identified. In this paper, a new observational avenue is suggested to detect individual PAHs, using their rotational line emission at radio frequencies. Previous PAH searches based on rotational spectroscopy have only targeted the bowl-shaped corannulene molecule, with the underlying assumption that other polar PAHs are triaxial and have a complex and diluted spectrum unusable for identification purposes. In this paper the rotational spectrum of quasi-symmetric PAHs is computed analytically. It is shown that the asymmetry of planar, nitrogen-substituted symmetric PAHs is small enough that their rotational spectrum, when observed with a resolution of about a MHz, has the appearance of a "comb" of evenly spaced stacks of lines. The simple pattern of these "comb" spectra allows for the use of matched-filtering techniques, which can result in a significantly enhanced signal-to-noise ratio. Detection forecasts are discussed for regions harbouring "anomalous microwave emission", believed to originate from the collective PAH rotational emission. A systematic search for PAH lines in various environments is advocated. If detected, PAH "combs" would allow to the conclusive and unambiguous identification of specific, free-floating interstellar PAHs.

The annihilation or decay of Dark Matter (DM) particles could affect the thermal history of the universe and leave an observable signature in Cosmic Microwave Background (CMB) anisotropies. We update constraints on the annihilation rate of DM particles in the smooth cosmological background, using WMAP7 and recent small-scale CMB data. With a systematic analysis based on the Press-Schechter formalism, we also show that DM annihilation in halos at small redshift may explain entirely the reionization patterns observed in the CMB, under reasonable assumptions concerning the concentration and formation redshift of halos. We find that a mixed reionization model based on DM annihilation in halos as well as star formation at a redshift z~6.5 could simultaneously account for CMB observations and satisfy constraints inferred from the Gunn-Peterson effect. However, these models tend to reheat the inter-galactic medium (IGM) well above observational bounds: by including a realistic prior on the IGM temperature at low redshift, we find stronger cosmological bounds on the annihilation cross-section than with the CMB alone.

Massive neutrinos make up a fraction of the dark matter, but due to their large thermal velocities, cluster significantly less than cold dark matter (CDM) on small scales. An accurate theoretical modelling of their effect during the non-linear regime of structure formation is required in order to properly analyse current and upcoming high-precision large-scale structure data, and constrain the neutrino mass. Taking advantage of the fact that massive neutrinos remain linearly clustered up to late times, this paper treats the linear growth of neutrino overdensities in a non-linear CDM background. The evolution of the CDM component is obtained via N-body computations. The smooth neutrino component is evaluated from that background by solving the Boltzmann equation linearised with respect to the neutrino overdensity. CDM and neutrinos are simultaneously evolved in time, consistently accounting for their mutual gravitational influence. This method avoids the issue of shot-noise inherent to particle-based neutrino simulations, and, in contrast with standard Fourier-space methods, properly accounts for the non-linear potential wells in which the neutrinos evolve. Inside the most massive late-time clusters, where the escape velocity is larger than the neutrino thermal velocity, neutrinos can clump non-linearly, causing the method to formally break down. It is shown that this does not affect the total matter power spectrum, which can be very accurately computed on all relevant scales up to the present time.

As atoms formed for the first time during primordial recombination, they emitted bound-bound and free-bound radiation leading to spectral distortions to the cosmic microwave background. These distortions might become observable in the future with high-sensitivity spectrometers, and provide a new window into physical conditions in the early universe. The standard multilevel atom method habitually used to compute the recombination spectrum is computationally expensive, impeding a detailed quantitative exploration of the information contained in spectral distortions thus far. In this work it is shown that the emissivity in optically thin allowed transitions can be factored into a computationally expensive but cosmology-independent part and a computationally cheap, cosmology-dependent part. The slow part of the computation consists in pre-computing temperature-dependent effective "conductances", linearly relating line or continuum intensity to departures from Saha equilibrium of the lowest-order excited states (2s and 2p), that can be seen as "voltages". The computation of these departures from equilibrium as a function of redshift is itself very fast, thanks to the effective multilevel atom method introduced in an earlier work. With this factorization, the recurring cost of a single computation of the recombination spectrum is only a fraction of a second on a standard laptop, more than four orders of magnitude shorter than standard computations. The spectrum from helium recombination can be efficiently computed in an identical way, and a fast code computing the full primordial recombination spectrum with this method will be made publicly available soon.

This article reviews the current status of theoretical modeling of electric dipole radiation from spinning dust grains. The fundamentally simple problem of dust grain rotation appeals to a rich set of concepts of classical and quantum physics, owing the the diversity of processes involved. Rotational excitation and damping rates through various mechanisms are discussed, as well as methods of computing the grain angular momentum distribution function. Assumptions on grain properties are reviewed. The robustness of theoretical predictions now seems mostly limited by the uncertainties regarding the grains themselves, namely their abundance, dipole moments, size and shape distribution.

Chern-Simons (CS) modified gravity is an extension to general relativity (GR) in which the metric is coupled to a scalar field, resulting in modified Einstein field equations. In the dynamical theory, the scalar field is itself sourced by the Pontryagin density of the space-time. In this paper, the coupled system of equations for the metric and the scalar field is solved numerically for slowly-rotating neutron stars described with realistic equations of state and for slowly-rotating black holes. An analytic solution for a constant-density nonrelativistic object is also presented. It is shown that the black hole solution cannot be used to describe the exterior spacetime of a star as was previously assumed. In addition, whereas previous analysis were limited to the small-coupling regime, this paper considers arbitrarily large coupling strengths. It is found that the CS modification leads to two effects on the gravitomagnetic sector of the metric: (i) Near the surface of a star or the horizon of a black hole, the magnitude of the gravitomagnetic potential is decreased and frame-dragging effects are reduced in comparison to GR. (ii) In the case of a star, the angular momentum J, as measured by distant observers, is enhanced in CS gravity as compared to standard GR. For a large coupling strength, the near-zone frame-dragging effects become significantly screened, whereas the far-zone enhancement saturate at a maximum value max(Delta J) ~ (M/R) J. Using measurements of frame-dragging effects around the Earth by Gravity Probe B and the LAGEOS satellites, a weak but robust constraint is set to the characteristic CS lengthscale, xi^{1/4} <~ 10^8 km.

Anomalous microwave emission is known to exist in the Perseus cloud. One of the most promising candidates to explain this excess of emission is electric dipole radiation from rapidly rotating very small dust grains, commonly referred to as spinning dust. Photometric data obtained with the Spitzer Space Telescope have been reprocessed and used in conjunction with the dust emission model DUSTEM to characterise the properties of the dust within the cloud. This analysis has allowed us to constrain spatial variations in the strength of the interstellar radiation field ($\chi_\mathrm{ISRF}$), the mass abundances of the PAHs and VSGs relative to the BGs (Y$_\mathrm{PAH}$ and Y$_\mathrm{VSG}$), the column density of hydrogen (N$_\mathrm{H}$) and the equilibrium dust temperature (T$_\mathrm{dust}$). The parameter maps of Y$_\mathrm{PAH}$, Y$_\mathrm{VSG}$ and $\chi_\mathrm{ISRF}$ are the first of their kind to be produced for the Perseus cloud, and we used these maps to investigate the physical conditions in which anomalous emission is observed. We find that in regions of anomalous emission the strength of the ISRF, and consequently the equilibrium temperature of the dust, is enhanced while there is no significant variation in the abundances of the PAHs and the VSGs or the column density of hydrogen. We interpret these results as an indication that the enhancement in $\chi_\mathrm{ISRF}$ might be affecting the properties of the small stochastically heated dust grains resulting in an increase in the spinning dust emission observed at 33 GHz. This is the first time that such an investigation has been performed, and we believe that this type of analysis creates a new perspective in the field of anomalous emission studies, and represents a powerful new tool for constraining spinning dust models.

One of the popular modifications to the theory of general relativity is non-dynamical Chern-Simons (CS) gravity, in which the metric is coupled to an externally prescribed scalar field. Setting accurate constraints to the parameters of the theory is important owing to their implications for the scalar field and/or the underlying fundamental theory. The current best constraints rely on measurements of the periastron precession rate in the double-binary-pulsar system and place a very tight bound on the characteristic CS lengthscale k_cs^{-1} <~ 3*10^{-9} km. This paper considers several effects that were not accounted for when deriving this bound and lead to a substantial suppression of the predicted rate of periastron precession. It is shown, in particular, that the point mass approximation for extended test bodies does not apply in this case. The constraint to the characteristic CS lengthscale is revised to k_cs^{-1} <~ 0.4 km, eight orders of magnitude weaker than what was previously found.

Standard big-bang nucleosynthesis (BBN) predicts only a trace abundance of lithium and no heavier elements, but some alternatives predict a nonzero primordial metallicity. Here we explore whether CMB measurements may set useful constraints to the primordial metallicity and/or whether the standard CMB calculations are robust, within the tolerance of forthcoming CMB maps, to the possibility of primordial metals. Metals would affect the recombination history (and thus CMB power spectra) in three ways: (1) Lyman-alpha photons can be removed (and recombination thus accelerated) by photoionizing metals. (2) The Bowen resonance-fluorescence mechanism may degrade Lyman-beta photons and thus enhance the Lyman-beta escape probability and speed up recombination. (3) Metals could affect the low-redshift tail of the CMB visibility function by providing additional free electrons. The last two of these provide the strongest CMB signal. However, the effects are detectable in the Planck satellite only if the primordial metal abundance is at least a few hundredths of solar for (2) and a few tenths of solar for (3). We thus conclude that Planck will not be able to improve upon current constraints to primordial metallicity, at the level of a thousandth of solar, from the Lyman-alpha forest and ultra-metal-poor halo stars, and that the CMB power-spectrum predictions for Planck suffer no uncertainty arising from the possibility that there may be primordial metals.

We present a state-of-the-art primordial recombination code, HyRec, including all the physical effects that have been shown to significantly affect recombination. The computation of helium recombination includes simple analytic treatments of hydrogen continuum opacity in the He I 2 1P - 1 1S line, the He I] 2 3P - 1 1S line, and treats feedback between these lines within the on-the-spot approximation. Hydrogen recombination is computed using the effective multilevel atom method, virtually accounting for an infinite number of excited states. We account for two-photon transitions from 2s and higher levels as well as frequency diffusion in Lyman-alpha with a full radiative transfer calculation. We present a new method to evolve the radiation field simultaneously with the level populations and the free electron fraction. These computations are sped up by taking advantage of the particular sparseness pattern of the equations describing the radiative transfer. The computation time for a full recombination history is ~2 seconds. This makes our code well suited for inclusion in Monte Carlo Markov chains for cosmological parameter estimation from upcoming high-precision cosmic microwave background anisotropy measurements.

The calculation of a highly accurate cosmological recombination history has been the object of particular attention recently, as it constitutes the major theoretical uncertainty when predicting the angular power spectrum of Cosmic Microwave Background anisotropies. Lyman transitions, in particular the Lyman-alpha line, have long been recognized as one of the bottlenecks of recombination, due to their very low escape probabilities. The Sobolev approximation does not describe radiative transfer in the vicinity of Lyman lines to a sufficient degree of accuracy, and several corrections have already been computed in other works. In this paper, the impact of some previously ignored radiative transfer effects is calculated. First, the effect of Thomson scattering in the vicinity of the Lyman-alpha line is evaluated, using a full redistribution kernel incorporated into a radiative transfer code. The effect of feedback of distortions generated by the optically thick deuterium Lyman-alpha line blueward of the hydrogen line is investigated with an analytic approximation. It is shown that both effects are negligible during cosmological hydrogen recombination. Secondly, the importance of high-lying, non overlapping Lyman transitions is assessed. It is shown that escape from lines above Ly-gamma and frequency diffusion in Ly-beta and higher lines can be neglected without loss of accuracy. Thirdly, a formalism generalizing the Sobolev approximation is developed to account for the overlap of the high-lying Lyman lines, which is shown to lead to negligible changes to the recombination history. Finally, the possibility of a cosmological hydrogen recombination maser is investigated. It is shown that there is no such maser in the purely radiative treatment presented here.

Cosmological hydrogen recombination has recently been the subject of renewed attention because of its importance for predicting the power spectrum of cosmic microwave background anisotropies. It has become clear that it is necessary to account for a large number n >~ 100 of energy shells of the hydrogen atom, separately following the angular momentum substates in order to obtain sufficiently accurate recombination histories. However, the multi-level atom codes that follow the populations of all these levels are computationally expensive, limiting recent analyses to only a few points in parameter space. In this paper, we present a new method for solving the multi-level atom recombination problem, which splits the problem into a computationally expensive atomic physics component that is independent of the cosmology, and an ultrafast cosmological evolution component. The atomic physics component follows the network of bound-bound and bound-free transitions among excited states and computes the resulting effective transition rates for the small set of "interface" states radiatively connected to the ground state. The cosmological evolution component only follows the populations of the interface states. By pre-tabulating the effective rates, we can reduce the recurring cost of multi-level atom calculations by more than 5 orders of magnitude. The resulting code is fast enough for inclusion in Markov Chain Monte Carlo parameter estimation algorithms. It does not yet include the radiative transfer or high-n two-photon processes considered in some recent papers. Further work on analytic treatments for these effects will be required in order to produce a recombination code usable for Planck data analysis.

We present a comprehensive treatment of the spectrum of electric dipole emission from spinning dust grains, updating the commonly used model of Draine and Lazarian. Grain angular velocity distributions are computed using the Fokker-Planck equation; we revisit the drift and diffusion coefficients for the major torques on the grain, including collisions, grain-plasma interactions, and infrared emission. We use updated grain optical properties and size distributions. The theoretical formalism is implemented in the companion code, SPDUST, which is publicly available. The effect of some environmental and grain parameters on the emissivity is shown and analysed.