• Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.
  • We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) power spectra, baryonic acoustic oscillation data, and local measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant $ G_{\rm eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2 T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first of these modes is consistent with the standard $\Lambda$CDM cosmology and corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} > 1.3\times10^5$. The other mode of the posterior, dubbed the "interacting neutrino mode", corresponds to neutrino decoupling occurring within a narrow redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is characterized by a high value of the effective neutrino coupling constant, together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of $H_0$, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality.
  • Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within lens galaxies is expected to be large, implying that their contribution to the lensing convergence field is approximately Gaussian and could thus be described by their power spectrum. We develop here a general formalism to compute from first principles the substructure convergence power spectrum for different populations of dark matter subhalos. As an example, we apply our framework to two distinct subhalo populations: a truncated Navarro-Frenk-White subhalo population motivated by standard CDM, and a truncated cored subhalo population motivated by self-interacting dark matter (SIDM). We study in detail how the subhalo abundance, mass function, internal density profile, and concentration affect the amplitude and shape of the substructure power spectrum. We determine that the power spectrum is mostly sensitive to a specific combination of the subhalo abundance and moments of the mass function, as well as to the average tidal truncation scale of the largest subhalos included in the analysis. Interestingly, we show that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the subhalos. In particular, the SIDM power spectrum exhibits a characteristic steepening at large wave number absent in the CDM power spectrum, opening the possibility of using this observable, if at all measurable, to discern between these two scenarios.
  • We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.
  • The amplitude of large-scale matter fluctuations inferred from the observed Sunyaev-Zeldovich (SZ) cluster mass function and from weak gravitational lensing studies, when taken at face value, is in tension with measurements of the cosmic microwave background (CMB) and baryon acoustic oscillation (BAO). In this work, we revisit whether this possible discrepancy can be attributed to new interactions in the dark matter sector. Focusing on a cosmological model where dark matter interacts with a dark radiation species until the epoch of matter-radiation equality, we find that measurements of the Lyman-alpha flux power spectrum from the Sloan Digital Sky Survey provides no support to the hypothesis that new dark matter interactions can resolve the possible tension between CMB and large-scale structure (LSS). Indeed, while the addition of dark matter-dark radiation interactions leads to an improvement of $2\Delta\ln\mathcal{L}=12$ with respect to the standard $\Lambda$ cold dark matter ($\Lambda$CDM) model when only CMB, BAO, and LSS data are considered, the inclusion of Lyman-alpha data reduces the improvement of the fit to $2\Delta\ln\mathcal{L}=6$ relative to $\Lambda$CDM. We thus conclude that the statistical evidence for new dark matter interactions (largely driven by the Planck SZ dataset) is marginal at best, and likely caused by systematics in the data. We also perform a Fisher forecast analysis for the reach of a future dataset composed of a CMB-S4 experiment combined with the Large Synoptic Survey Telescope galaxy survey. We find that the constraint on the effective number of fluid-like dark radiation species, $\Delta N_{\rm fluid}$, will be improved by an order of magnitude compared to current bounds.
  • Recently it was shown that dark matter with mass of order the weak scale can be charged under a new long-range force, decoupled from the Standard Model, with only weak constraints from early Universe cosmology. Here we consider the implications of an additional charged particle $C$ that is light enough to lead to significant dissipative dynamics on galactic times scales. We highlight several novel features of this model, which can be relevant even when the $C$ particle constitutes only a small fraction of the number density (and energy density). We assume a small asymmetric abundance of the $C$ particle whose charge is compensated by a heavy $X$ particle so that the relic abundance of dark matter consists mostly of symmetric $X$ and $\bar{X}$, with a small asymmetric component made up of $X$ and $C$. As the universe cools, it undergoes asymmetric recombination binding the free $C$s into $(XC)$ dark atoms efficiently. Even with a tiny asymmetric component, the presence of $C$ particles catalyzes tight coupling between the heavy dark matter $X$ and the dark photon plasma that can lead to a significant suppression of the matter power spectrum on small scales and lead to some of the strongest bounds on such dark matter theories. We find a viable parameter space where structure formation constraints are satisfied and significant dissipative dynamics can occur in galactic haloes but show a large region is excluded. Our model shows that subdominant components in the dark sector can dramatically affect structure formation.
  • We revisit constraints on dark matter that is charged under a $U(1)$ gauge group in the dark sector, decoupled from Standard Model forces. We find that the strongest constraints in the literature are subject to a number of mitigating factors. For instance, the naive dark matter thermalization timescale in halos is corrected by saturation effects that slow down isotropization for modest ellipticities. The weakened bounds uncover interesting parameter space, making models with weak-scale charged dark matter viable, even with electromagnetic strength interaction. This also leads to the intriguing possibility that dark matter self-interactions within small dwarf galaxies are extremely large, a relatively unexplored regime in current simulations. Such strong interactions suppress heat transfer over scales larger than the dark matter mean free path, inducing a dynamical cutoff length scale above which the system appears to have only feeble interactions. These effects must be taken into account to assess the viability of darkly-charged dark matter. Future analyses and measurements should probe a promising region of parameter space for this model.
  • The Theia Collaboration: Celine Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escude, Alexis Brandeker, Frederic Courbin, Torsten Ensslin, Antonio Falcao, Katherine Freese, Berry Holl, Lucas Labadie, Alain Leger, Fabien Malbet, Gary Mamon, Barbara McArthur, Alcione Mora, Michael Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Conrado Albertus, Stefano Bertone, Herve Bouy, Michael Boylan-Kolchin, Anthony Brown, Warren Brown, Vitor Cardoso, Laurent Chemin, Riccardo Claudi, Alexandre C. M. Correia, Mariateresa Crosta, Antoine Crouzier, Francis-Yan Cyr-Racine, Mario Damasso, Antonio da Silva, Melvyn Davies, Payel Das, Pratika Dayal, Miguel de Val-Borro, Antonaldo Diaferio, Adrienne Erickcek, Malcolm Fairbairn, Morgane Fortin, Malcolm Fridlund, Paulo Garcia, Oleg Gnedin, Ariel Goobar, Paulo Gordo, Renaud Goullioud, Nigel Hambly, Nathan Hara, David Hobbs, Erik Hog, Andrew Holland, Rodrigo Ibata, Carme Jordi, Sergei Klioner, Sergei Kopeikin, Thomas Lacroix, Jacques Laskar, Christophe Le Poncin-Lafitte, Xavier Luri, Subhabrata Majumdar, Valeri Makarov, Richard Massey, Bertrand Mennesson, Daniel Michalik, Andre Moitinho de Almeida, Ana Mourao, Leonidas Moustakas, Neil Murray, Matthew Muterspaugh, Micaela Oertel, Luisa Ostorero, Angeles Perez-Garcia, Imants Platais, Jordi Portell i de Mora, Andreas Quirrenbach, Lisa Randall, Justin Read, Eniko Regos, Barnes Rory, Krzysztof Rybicki, Pat Scott, Jean Schneider, Jakub Scholtz, Arnaud Siebert, Ismael Tereno, John Tomsick, Wesley Traub, Monica Valluri, Matt Walker, Nicholas Walton, Laura Watkins, Glenn White, Dafydd Wyn Evans, Lukasz Wyrzykowski, Rosemary Wyse
    July 2, 2017 astro-ph.IM
    In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.
  • Oct. 10, 2016 hep-th, hep-ph, gr-qc, astro-ph.CO
    This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.
  • In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.
  • We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported X-ray signals that could be due to sterile neutrino decay, we consider $7$ keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures $T \gtrsim 100$ MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark-hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at $k \simeq 10 \ h \ {\rm Mpc}^{-1}$, that are relevant to galactic structure formation.
  • In the standard model neutrinos are assumed to have streamed across the Universe since they last scattered at the weak decoupling epoch when the temperature of the standard-model plasma was ~MeV. The shear stress of free-streaming neutrinos imprints itself gravitationally on the Cosmic Microwave Background (CMB) and makes the CMB a sensitive probe of neutrino scattering. Yet, the presence of nonstandard physics in the neutrino sector may alter this standard chronology and delay neutrino free-streaming until a much later epoch. We use observations of the CMB to constrain the strength of neutrino self-interactions G_eff and put limits on new physics in the neutrino sector from the early Universe. Recent measurements of the CMB at large multipoles made by the Planck satellite and high-l experiments are critical for probing this physics. Within the context of conventional LambdaCDM parameters cosmological data are compatible with G_eff < 1/(56 MeV)^2 and neutrino free-streaming might be delayed until their temperature has cooled to as low as ~25 eV. Intriguingly, we also find an alternative cosmology compatible with cosmological data in which neutrinos scatter off each other until z~10^4 with a preferred interaction strength in a narrow region around $G_{\rm eff} \simeq 1/({\rm 10 \, MeV})^{2} \simeq 8.6\times10^8 G_{\rm F}$, where $G_{\rm F}$ is the Fermi constant. This distinct self-interacting neutrino cosmology is characterized by somewhat lower values of both the scalar spectral index and the amplitude of primordial fluctuations. While we phrase our discussion here in terms of a specific scenario in which a late onset of neutrino free-streaming could occur, our constraints on the neutrino visibility function are very general.
  • We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.
  • If all or a fraction of the dark matter (DM) were coupled to a bath of dark radiation (DR) in the early Universe we expect the combined DM-DR system to give rise to acoustic oscillations of the dark matter until it decouples from the DR. Much like the standard baryon acoustic oscillations, these dark acoustic oscillations (DAO) imprint a characteristic scale, the sound horizon of dark matter, on the matter power spectrum. We compute in detail how the microphysics of the DM-DR interaction affects the clustering of matter in the Universe and show that the DAO physics also gives rise to unique signatures in the temperature and polarization spectra of the cosmic microwave background (CMB). We use cosmological data from the CMB, baryon acoustic oscillations (BAO), and large-scale structure to constrain the possible fraction of interacting DM as well as the strength of its interaction with DR. Like nearly all knowledge we have gleaned about dark matter since inferring its existence this constraint rests on the betrayal by gravity of the location of otherwise invisible dark matter. Although our results can be straightforwardly applied to a broad class of models that couple dark matter particles to various light relativistic species, in order to make quantitative predictions, we model the interacting component as dark atoms coupled to a bath of dark photons. We find that linear cosmological data and CMB lensing put strong constraints on existence of DAO features in the CMB and the large-scale structure of the Universe. Interestingly, we find that at most ~5% of all DM can be very strongly interacting with DR. We show that our results are surprisingly constraining for the recently proposed Double-disk DM model, a novel example of how large-scale precision cosmological data can be used to constrain galactic physics and sub-galactic structure.
  • While, to ensure successful cosmology, dark matter (DM) must kinematically decouple from the standard model plasma very early in the history of the Universe, it can remain coupled to a bath of "dark radiation" until a relatively late epoch. One minimal theory that realizes such a scenario is the Atomic Dark Matter model, in which two fermions oppositely charged under a new U(1) dark force are initially coupled to a thermal bath of "dark photons" but eventually recombine into neutral atom-like bound states and begin forming gravitationally-bound structures. As dark atoms have (dark) atom-sized geometric cross sections, this model also provides an example of self-interacting DM with a velocity-dependent cross section. Delayed kinetic decoupling in this scenario predicts novel DM properties on small scales but retains the success of cold DM on larger scales. We calculate the atomic physics necessary to capture the thermal history of this dark sector and show significant improvements over the standard atomic hydrogen calculation are needed. We solve the Boltzmann equations that govern the evolution of cosmological fluctuations in this model and find in detail the impact of the atomic DM scenario on the matter power spectrum and the cosmic microwave background (CMB). This scenario imprints a new length scale, the Dark-Acoustic-Oscillation (DAO) scale, on the matter density field. This DAO scale shapes the small-scale matter power spectrum and determines the minimal DM halo mass at late times which may be many orders of magnitude larger than in a typical WIMP scenario. This model necessarily includes an extra dark radiation component, which may be favoured by current CMB experiments, and we quantify CMB signatures that distinguish an atomic DM scenario from a standard $\Lambda$CDM model containing extra free-streaming particles. [Abridged]
  • Many models of inflation predict oscillatory features in the bispectrum of primordial fluctuations. Since it has been shown that primordial non-Gaussianity can lead to a scale-dependent halo bias, we investigate the effect of oscillations in the three-point function on the clustering of dark-matter halos. Interestingly, we find that features in the inflaton potential such as oscillations or sharp steps get imprinted in the mass dependence of the non-Gaussian halo bias. In this paper, we focus on models displaying a sharp feature in the inflaton potential as well as Resonant non-Gaussianity. In both cases, we find a strong scale dependence for the non-Gaussian halo bias with a slope similar to that of the local model. In the resonant case, we find that the non-Gaussian bias oscillates with halo mass, a novel feature that is unique to this type of models. In the case of a sharp feature in the inflaton potential, we find that the clustering of halos is enhanced at the mass scale corresponding to the Fourier mode that exited the horizon when the inflaton was crossing the feature in the potential. Both of these are new effects that open the possibility of characterizing the inflationary potential with large-scale-structure surveys. We briefly discuss the prospects for detecting these non-Gaussian effects.
  • Prior to recombination photons, electrons, and atomic nuclei rapidly scattered and behaved, almost, like a single tightly-coupled photon-baryon plasma. We investigate here the accuracy of the tight-coupling approximation commonly used to numerically evolve the baryon and photon perturbation equations at early times. By solving the exact perturbations equations with a stiff solver starting deep in the radiation-dominated epoch we find the level of inaccuracy introduced by resorting to the standard first-order tight-coupling approximation. We develop a new second-order approximation in the inverse Thomson opacity expansion and show that it closely tracks the full solution, at essentially no extra numerical cost. We find the bias on estimates of cosmological parameters introduced by the first-order approximation is, for most parameters, negligible. Finally, we show that our second-order approximation can be used to reduce the time needed to compute cosmic microwave background angular spectra by as much as ~17%.
  • Reheating is an important part of inflationary cosmology. It describes the production of Standard Matter particles after the phase of accelerated expansion. We give a review of the reheating process, focusing on an in-depth discussion of the preheating stage which is characterized by exponential particle production due to a parametric resonance or tachyonic instability. We give a brief overview of the thermalization process after preheating and end with a survey of some applications to supersymmetric theories and to other issues in cosmology such as baryogenesis, dark matter and metric preheating.
  • It has recently been argued that the PAMELA, ATIC and PPB-BETS data showing an anomalous excess of high-energy cosmic ray positrons and electrons might be explained by dark matter annihilating in the Galactic halo with a cross section resonantly enhanced compared to its value in the primeval plasma. We find that with a very large annihilation cross section the flash of energetic photons and electron-positron pairs expected from dark-matter annihilation in the first protohalos that form at redshift z~40 is likely substantial and observable. As a consequence, bounds on the allowed energy injection into the primordial gas and the energy density of the diffuse gamma-ray background give rise to limits on the low-velocity dark matter cross section that can be difficult to reconcile with this interpretation of the PAMELA, ATIC and PPB-BETS results.
  • We study how the resonant decay of moduli fields arising in the Minimal Supersymmetric Standard Model (MSSM) could affect large scale curvature perturbations in the early universe. It has been known for some time that the presence of entropy perturbations in a multi-component system can act as seeds for the curvature perturbations on all scales. These entropy perturbations could be amplified exponentially if one of the moduli decays via stochastic resonance, affecting the curvature power spectrum in the process. By imposing the COBE normalization on this power spectrum, one could put constraints on the masses and couplings of the underlying particle physics model without having to rely on collider experiments. We discuss in detail the case of the MSSM but this method could be applied to other theories beyond the Standard Model.