Source author record

Annika H. G. Peter

Annika H. G. Peter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

(Not as) Big as a Barn: Upper Bounds on Dark Matter-Nucleus Cross Sections

Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $σ_{χA} \propto A^2 μ_A^2 σ_{χN}\propto A^4 σ_{χN}$, where the $A^2$ comes from coherence and the $μ_A^2\simeq A^2 m_N^2$ from kinematics for $m_χ\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $σ_{χA} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $σ_{χN}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $σ_{χN}>10^{-32}\;\text{cm}^2$ with $m_χ\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.

preprint2022arXiv

Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo

Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $\lesssim 10^8 M_\odot$ subhalos in lens-mass systems ($\sim 10^{13}M_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhalos with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhalos ($\lesssim 1/1000$ host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhalos to core-collapse, finding that it is nearly impossible to drive core-collapse in subhalos in SIDM models with constant cross sections. Any discovery of ultra-compact dark substructures with future substructure lensing observations favors additional degrees of freedom, such as velocity-dependence, in the cross section.

preprint2022arXiv

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

preprint2022arXiv

Snowmass2021 Theory Frontier White Paper: Data-Driven Cosmology

Over the past few decades, astronomical and cosmological data sets firmly established the existence of physics beyond the Standard Model of particle physics by providing strong evidence for the existence of dark matter, dark energy, and non-zero neutrino mass. In addition, the generation of primordial perturbations most likely also relies on physics beyond the Standard Model of particle physics. Theory work, ranging from models of the early universe in string theory that that led to novel phenomenological predictions to the development of effective field theories to large-scale cosmological simulations that include the physics of galaxy evolution, has played a key role in analyzing and interpreting these data sets and in suggesting novel paths forward to isolate the origins of this new physics. Over the next decade, even more sensitive surveys are beginning to take data and are being planned. In this white paper, we describe key areas of the theory program that will be needed to optimize the physics return on investment from these new observational opportunities.

preprint2022arXiv

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities.

preprint2022arXiv

The Milky Way satellite velocity function is a sharp probe of small-scale structure problems

Twenty years ago, the mismatch between the observed number of Milky Way satellite galaxies and the predicted number of cold dark matter (CDM) subhalos was dubbed the ``missing satellites problem". Although mostly framed since in terms of satellite counts in luminosity space, the missing satellites problem was originally posed in velocity space. Importantly, the stellar velocity dispersion function encodes information about the density profile of satellites as well as their abundance. In this work, we completeness correct the MW satellite stellar velocity dispersion function down to its ultrafaint dwarfs ($L \gtrsim 340$ L$_\odot$). Our most conservative completeness correction is in good agreement with a simple CDM model in which massive, classical satellites (M$_{\rm 200} \gtrsim 5 \times 10^8~$M$_\odot$) have baryon-driven cores, while lower mass, ultrafaint satellites inhabit cuspy halos that are not strongly tidally stripped. Tidal destruction of satellites by the MW's disk must be minimal, otherwise the completeness-corrected velocity function exceeds any plausible CDM prediction -- a ``too many satellites" problem. We rule out non-core-collapsing self-interacting dark matter models with a constant cross section $\gtrsim$ 0.1 cm$^2$/g. Constraints on warm dark matter are stronger than those based on the luminosity function due to its additional sensitivity to subhalo central densities, which suppresses number counts by up to an additional 30%. A thermal relic mass $\gtrsim$ 6 keV is preferred. Reducing uncertainties on stellar velocity dispersion measurements and the amount of tidal stripping experienced by the faintest dwarfs is key to determining the severity of the too many satellites problem.

preprint2020arXiv

Luminosity Functions and Host-to-Host Scatter of Dwarf Satellite Systems in the Local Volume

Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation (SBF) distances to recently cataloged candidate dwarf satellites around 10 massive hosts within $D<12$ Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to $M_g{\sim}-9$ to $-10$, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting a SHMR that matches hydrodynamic simulations, the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function amongst LV hosts. Contrary to recent results, we find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, necessitating a higher normalization of the SHMR around halo masses of $10^{11}$ \msun\ than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy-subhalo connection in the low-mass regime.

preprint2019arXiv

The Case for Strangulation in Low-Mass Hosts: DDO 113

We investigate the case for environmental quenching of the Fornax-mass satellite DDO 113, which lies only 9 kpc in projection from its host, the Large-Magellanic-Cloud-mass galaxy NGC 4214. DDO 113 was quenched about 1 Gyr ago and is virtually gas-free, while analogs in the field are predominantly star-forming and gas-rich. We use deep imaging obtained with the Large Binocular Telescope to show that DDO 113 exhibits no evidence of tidal disruption to a surface brightness of $μ_V\sim29$ mag $\text{arcsec}^{-2}$, based on both unresolved emission and resolved stars. Mass-analogs of DDO 113 in Illustris-1 with similar hosts, small projected separations, and no significant tidal stripping first fell into their host halo 2--6 Gyr ago, showing that tidal features (or lack thereof) can be used to constrain infall times in systems where there are few other constraints on the orbit of the satellite. With the infall time setting the clock for environmental quenching mechanisms, we investigate the plausibility of several such mechanisms. We find that strangulation, the cessation of cold gas inflows, is likely the dominant quenching mechanism for DDO 113, requiring a time-averaged mass-loading factor of $η=6-11$ for star-formation-driven outflows that is consistent with theoretical and observational constraints. Motivated by recent numerical work, we connect DDO 113's strangulation to the presence of a cool circumgalactic medium (CGM) around NGC 4214. This discovery shows that the CGM of low-mass galaxies can affect their satellites significantly and motivates further work on understanding the baryon cycle in low-mass galaxies.

preprint2016arXiv

Enhanced tidal stripping of satellites in the galactic halo from dark matter self-interactions

We investigate the effects of self-interacting dark matter (SIDM) on the tidal stripping and evaporation of satellite galaxies in a Milky Way-like host. We use a suite of five zoom-in, dark-matter-only simulations, two with velocity-independent SIDM cross sections, two with velocity-dependent SIDM cross sections, and one cold dark matter simulation for comparison. After carefully assigning stellar mass to satellites at infall, we find that stars are stripped at a higher rate in SIDM than in CDM. In contrast, the total bound dark matter mass loss rate is minimally affected, with subhalo evaporation having negligible effects on satellites for viable SIDM models. Centrally located stars in SIDM haloes disperse out to larger radii as cores grow. Consequently, the half-light radius of satellites increases, stars become more vulnerable to tidal stripping, and the stellar mass function is suppressed. We find that the ratio of core radius to tidal radius accurately predicts the relative strength of enhanced SIDM stellar stripping. Velocity-independent SIDM models show a modest increase in the stellar stripping effect with satellite mass, whereas velocity-dependent SIDM models show a large increase in this effect towards lower masses, making observations of ultra-faint dwarfs prime targets for distinguishing between and constraining SIDM models. Due to small cores in the largest satellites of velocity-dependent SIDM, no identifiable imprint is left on the all-sky properties of the stellar halo. While our results focus on SIDM, the main physical mechanism of enhanced tidal stripping of stars apply similarly to satellites with cores formed via other means.

preprint2016arXiv

First Observation of Time Variation in the Solar-Disk Gamma-Ray Flux with Fermi

The solar disk is a bright gamma-ray source. Surprisingly, its flux is about one order of magnitude higher than predicted. As a first step toward understanding the physical origin of this discrepancy, we perform a new analysis in 1-100 GeV using 6 years of public Fermi-LAT data. Compared to the previous analysis by the Fermi Collaboration, who analyzed 1.5 years of data and detected the solar disk in 0.1-10 GeV, we find two new and significant results: 1. In the 1-10 GeV flux (detected at $>5σ$), we discover a significant time variation that anticorrelates with solar activity. 2. We detect gamma rays in 10-30 GeV at $>5σ$, and in 30-100 GeV at $> 2σ$. The time variation strongly indicates that solar-disk gamma rays are induced by cosmic rays and that solar atmospheric magnetic fields play an important role. Our results provide essential clues for understanding the underlying gamma-ray production processes, which may allow new probes of solar atmospheric magnetic fields, cosmic rays in the solar system, and possible new physics. Finally, we show that the Sun is a promising new target for ground-based TeV gamma-ray telescopes such as HAWC and LHAASO.

preprint2016arXiv

First results from the MADCASH Survey: A Faint Dwarf Galaxy Companion to the Low Mass Spiral Galaxy NGC 2403 at 3.2 Mpc

We report the discovery of the faintest known dwarf galaxy satellite of an LMC stellar-mass host beyond the Local Group, based on deep imaging with Subaru/Hyper Suprime-Cam. MADCASH J074238+652501-dw lies $\sim$35 kpc in projection from NGC 2403, a dwarf spiral galaxy at $D$$\approx$3.2 Mpc. This new dwarf has $M_{g} = -7.4\pm0.4$ and a half-light radius of $168\pm70$ pc, at the calculated distance of $3.39\pm0.41$ Mpc. The color-magnitude diagram reveals no evidence of young stellar populations, suggesting that MADCASH J074238+652501-dw is an old, metal-poor dwarf similar to low luminosity dwarfs in the Local Group. The lack of either detected HI gas ($M_{\rm HI}/L_{V} < 0.69 M_\odot/L_\odot$, based on Green Bank Telescope observations) or $GALEX$ NUV/FUV flux enhancement is consistent with a lack of young stars. This is the first result from the MADCASH (Magellanic Analog Dwarf Companions And Stellar Halos) survey, which is conducting a census of the stellar substructure and faint satellites in the halos of Local Volume LMC analogs via resolved stellar populations. Models predict a total of $\sim$4-10 satellites at least as massive as MADCASH J074238+652501-dw around a host with the mass of NGC 2403, with 2-3 within our field of view, slightly more than the one such satellite observed in our footprint.

preprint2016arXiv

The connection between the host halo and the satellite galaxies of the Milky Way

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

preprint2015arXiv

Identifying the Theory of Dark Matter with Direct Detection

Identifying the true theory of dark matter depends crucially on accurately characterizing interactions of dark matter (DM) with other species. In the context of DM direct detection, we present a study of the prospects for correctly identifying the low-energy effective DM-nucleus scattering operators connected to UV-complete models of DM-quark interactions. We take a census of plausible UV-complete interaction models with different low-energy leading-order DM-nuclear responses. For each model (corresponding to different spin-, momentum-, and velocity-dependent responses), we create a large number of realizations of recoil-energy spectra, and use Bayesian methods to investigate the probability that experiments will be able to select the correct scattering model within a broad set of competing scattering hypotheses. We conclude that agnostic analysis of a strong signal (such as Generation-2 would see if cross sections are just below the current limits) seen on xenon and germanium experiments is likely to correctly identify momentum dependence of the dominant response, ruling out models with either "heavy" or "light" mediators, and enabling downselection of allowed models. However, a unique determination of the correct UV completion will critically depend on the availability of measurements from a wider variety of nuclear targets, including iodine or fluorine. We investigate how model-selection prospects depend on the energy window available for the analysis. In addition, we discuss accuracy of the DM particle mass determination under a wide variety of scattering models, and investigate impact of the specific types of particle-physics uncertainties on prospects for model selection.

preprint2014arXiv

Core Formation in Dwarf Halos with Self Interacting Dark Matter: No Fine-Tuning Necessary

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range σ/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, σ/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (σ/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

preprint2014arXiv

Cosmological Simulations of Decaying Dark Matter: Implications for Small-scale Structure of Dark Matter Halos

We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Milky Way-mass halos, we use high-resolution, zoom-in simulations to explore the effects of decays on Galactic substructure. In general, halos with circular velocities comparable to the magnitude of kick velocity are most strongly affected by decays. We show that decaying dark matter models with lifetimes comparable to Hubble time and recoil speeds about 20-40 km/s can significantly reduce both the abundance of Galactic subhalos and the internal densities of the subhalos. We also compare subhalo circular velocity profiles with observational constraints on the Milky Way dwarf satellite galaxies. Interestingly, we find that decaying dark matter models that do not violate current astrophysical constraints, can significantly mitigate both the well-documented "missing satellites problem" and the more recent "too big to fail problem" associated with the abundances and densities of Local Group dwarf satellite galaxies. A relatively unique feature of late decaying dark matter models is that they predict significant evolution of halos as a function of time. This is an important consideration because at high redshifts, prior to decays, decaying models exhibit the same sequence of structure formation as cold dark matter. We conclude that models of decaying dark matter make predictions that are relevant for the interpretation of observations of small galaxies in the Local Group and can be tested or constrained by the kinematics of Local Group dwarf galaxies as well as by forthcoming large-scale surveys.

preprint2014arXiv

Effect of Gravitational Focusing on Annual Modulation in Dark-Matter Direct-Detection Experiments

The scattering rate at dark-matter direct-detection experiments should modulate annually due to the Earth's orbit around the Sun. The rate is typically thought to be extremized around June 1, when the relative velocity of the Earth with respect to the dark-matter wind is maximal. We point out that gravitational focusing can alter this modulation phase. Unbound dark-matter particles are focused by the Sun's gravitational potential, affecting their phase-space density in the lab frame. Gravitational focusing can result in a significant overall shift in the annual-modulation phase, which is most relevant for dark matter with low scattering speeds. The induced phase shift for light O(10) GeV dark matter may also be significant, depending on the threshold energy of the experiment.

preprint2014arXiv

Too Many, Too Few, or Just Right? The Predicted Number and Distribution of Milky Way Dwarf Galaxies

We predict the spatial distribution and number of Milky Way dwarf galaxies to be discovered in the DES and LSST surveys, by completeness correcting the observed SDSS dwarf population. We apply most massive in the past, earliest forming, and earliest infall toy models to a set of dark matter-only simulated Milky Way/M31 halo pairs from Exploring the Local Volume In Simulations (ELVIS). The observed spatial distribution of Milky Way dwarfs in the LSST-era will discriminate between the earliest infall and other simplified models for how dwarf galaxies populate dark matter subhalos. Inclusive of all toy models and simulations, at 90% confidence we predict a total of 37-114 L $\gtrsim 10^3$L$_{\odot}$ dwarfs and 131-782 L $\lesssim 10^3$L$_{\odot}$ dwarfs within 300 kpc. These numbers of L $\gtrsim 10^3$L$_{\odot}$ dwarfs are dramatically lower than previous predictions, owing primarily to our use of updated detection limits and the decreasing number of SDSS dwarfs discovered per sky area. For an effective $r_{\rm limit}$ of 25.8 mag, we predict: 3-13 L $\gtrsim 10^3$L$_{\odot}$ and 9-99 L $\lesssim 10^3$L$_{\odot}$ dwarfs for DES, and 18-53 L $\gtrsim 10^3$L$_{\odot}$ and 53-307 L $\lesssim 10^3$L$_{\odot}$ dwarfs for LSST. These enormous predicted ranges ensure a coming decade of near-field excitement with these next generation surveys.

preprint2014arXiv

Understanding WIMP-baryon interactions with direct detection: A Roadmap

We study prospects of dark-matter direct-detection searches for probing non-relativistic effective theory for WIMP-baryon scattering. We simulate a large set of noisy recoil-energy spectra for different scattering scenarios (beyond the standard momentum-independent contact interaction), for Generation 2 and futuristic experiments. We analyze these simulations and quantify the probability of successfully identifying the operator governing the scattering, if a WIMP signal is observed. We find that the success rate depends on a combination of factors: the WIMP mass, the mediator mass, the type of interaction, and the experimental energy window. For example, for a 20 GeV WIMP, Generation 2 is only likely to identify the right operator if the interaction is Coulomb-like, and is unlikely to do so in any other case. For a WIMP with a mass of 200 GeV or higher, success is almost guaranteed. We also find that, regardless of the scattering model and the WIMP parameters, a single Generation 2 experiment is unlikely to successfully discern the momentum dependence of the underlying operator on its own, but prospects improve drastically when experiments with different target materials and energy windows are analyzed jointly. Furthermore, we examine the quality of parameter estimation and degeneracies in the multi-dimensional parameter space of the effective theory. We find in particular that the resulting WIMP mass estimates can be severely biased if data are analyzed assuming the standard (momentum-independent) operator while the actual operator has momentum-dependence. Finally, we evaluate the ultimate reach of direct detection, finding that the prospects for successful operator selection prior to reaching the irreducible backgrounds are excellent, if the signal is just below the current limits, but slim if Generation 2 does not report WIMP detection.

preprint2013arXiv

Cold dark matter: controversies on small scales

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

preprint2013arXiv

Lyman-alpha Forest Constraints on Decaying Dark Matter

We present an analysis of high-resolution N-body simulations of decaying dark matter cosmologies focusing on the statistical properties of the transmitted Lyman-alpha forest flux in the high-redshift intergalactic medium. In this type of model a dark matter particle decays into a slightly less massive stable dark matter daughter particle and a comparably light particle. The small mass splitting provides a non-relativistic kick velocity V_k to the daughter particle resulting in free-streaming and subsequent damping of small-scale density fluctuations. Current Lyman-alpha forest power spectrum measurements probe comoving scales up to ~ 2-3 h^-1 Mpc at redshifts z ~ 2-4, providing one of the most robust ways to probe cosmological density fluctuations on relatively small scales. The suppression of structure growth due to the free-streaming of dark matter daughter particles also has a significant impact on the neutral hydrogen cloud distribution, which traces the underlying dark matter distribution well at high redshift. We exploit Lyman-alpha forest power spectrum measurements to constrain the amount of free-streaming of dark matter in such models and thereby place limits on decaying dark matter based only on the dynamics of cosmological perturbations without any assumptions about the interactions of the decay products. We find that SDSS 1D Lyman-alpha forest power spectrum data place a lifetime-dependent upper limit V_k < 30-70 km/s for decay lifetimes < 10 Gyr. This is the most stringent model-independent bound on invisible dark matter decays with small mass splittings. For large mass splittings (large V_k), Lyman-alpha forest data restrict the dark matter lifetime to Gamma^-1 > 40 Gyr. Forthcoming BOSS data should be able to provide more stringent constraints on exotic dark matter, mainly because the larger BOSS quasar spectrum sample will significantly reduce statistical errors.

preprint2013arXiv

WIMP physics with ensembles of direct-detection experiments

The search for weakly-interacting massive particle (WIMP) dark matter is multi-pronged. Ultimately, the WIMP-dark-matter picture will only be confirmed if different classes of experiments see consistent signals and infer the same WIMP properties. In this work, we review the ideas, methods, and status of direct-detection searches. We focus in particular on extracting WIMP physics (WIMP interactions and phase-space distribution) from direct-detection data in the early discovery days when multiple experiments see of order dozens to hundreds of events. To demonstrate the essential complementarity of different direct-detection experiments in this context, we create mock data intended to represent the data from the near-future Generation 2 experiments. We consider both conventional supersymmetry-inspired benchmark points (with spin-independent and -dependent elastic cross sections just below current limits), as well as benchmark points for other classes of models (inelastic and effective-operator paradigms). We also investigate the effect on parameter estimation of loosening or dropping the assumptions about the local WIMP phase-space distribution. We arrive at two main conclusions. Firstly, teasing out WIMP physics with experiments depends critically on having a wide set of detector target materials, spanning a large range of target nuclear masses and spin-dependent sensitivity. It is also highly desirable to obtain data from low-threshold experiments. Secondly, a general reconstruction of the local WIMP velocity distribution, which will only be achieved if there are multiple experiments using different target materials, is critical to obtaining a robust and unbiased estimate of the WIMP mass.

preprint2012arXiv

Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass σ/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in σ/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with σ/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

preprint2012arXiv

Cosmological Simulations with Self-Interacting Dark Matter II: Halo Shapes vs. Observations

If dark matter has a large self-interaction scattering cross section, then interactions among dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. Work in the past has used this effect to rule out velocity-independent, elastic cross sections larger than sigma/m ~ 0.02 cm^2/g based on comparisons to the shapes of galaxy cluster lensing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show that these constraints were off by more than an order of magnitude because (a) they did not properly account for the fact that the observed ellipticity gets contributions from the triaxial mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and (c) the core region retains more of its triaxial nature than estimated before. Including these effects properly shows that the same observations now allow dark matter self-interaction cross sections at least as large as sigma/m = 0.1 cm^2/g. We show that constraints on self-interacting dark matter from strong-lensing clusters are likely to improve significantly in the near future, but possibly more via central densities and core sizes than halo shapes.

preprint2012arXiv

Dark Matter: A Brief Review

From astronomical observations, we know that dark matter exists, makes up 23% of the mass budget of the Universe, clusters strongly to form the load-bearing frame of structure for galaxy formation, and hardly interacts with ordinary matter except gravitationally. However, this information is not enough to identify the particle specie(s) that make up dark matter. As such, the problem of determining the identity of dark matter has largely shifted to the fields of astroparticle and particle physics. In this talk, I will review the current status of the search for the nature of dark matter. I will provide an introduction to possible particle candidates for dark matter and highlight recent experimental astroparticle- and particle-physics results that constrain the properties of those candidates. Given the absence of detections in those experiments, I will advocate a return of the problem of dark-matter identification to astronomy, and show what kinds of theoretical and observational work might be used to pin down the nature of dark matter once and for all. This talk is intended for a broad astronomy audience.

preprint2012arXiv

Fundamental statistical limitations of future dark matter direct detection experiments

We discuss irreducible statistical limitations of future ton-scale dark matter direct detection experiments. We focus in particular on the coverage of confidence intervals, which quantifies the reliability of the statistical method used to reconstruct the dark matter parameters, and the bias of the reconstructed parameters. We study 36 benchmark dark matter models within the reach of upcoming ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover the true values of the WIMP parameters, and are hence conservative. We evaluate the probability that unavoidable statistical fluctuations in the data might lead to a biased reconstruction of the dark matter parameters, or large uncertainties on the reconstructed parameter values. We show that this probability can be surprisingly large, even for benchmark models leading to a large event rate of order a hundred counts. We find that combining data sets from two different targets leads to improved coverage properties, as well as a substantial reduction of statistical bias and uncertainty on the dark matter parameters.

preprint2012arXiv

Probing the Local Velocity Distribution of WIMP Dark Matter with Directional Detectors

We explore the ability of directional nuclear-recoil detectors to constrain the local velocity distribution of weakly interacting massive particle (WIMP) dark matter by performing Bayesian parameter estimation on simulated recoil-event data sets. We discuss in detail how directional information, when combined with measurements of the recoil-energy spectrum, helps break degeneracies in the velocity-distribution parameters. We also consider the possibility that velocity structures such as cold tidal streams or a dark disk may also be present in addition to the Galactic halo. Assuming a carbon-tetrafluoride detector with a 30-kg-yr exposure, a 50-GeV WIMP mass, and a WIMP-nucleon spin-dependent cross-section of 0.001 pb, we show that the properties of a cold tidal stream may be well constrained. However, measurement of the parameters of a dark-disk component with a low lag speed of ~50 km/s may be challenging unless energy thresholds are improved.

preprint2012arXiv

WIMP astronomy and particle physics with liquid-noble and cryogenic direct-detection experiments

Once weakly-interacting massive particles (WIMPs) are unambiguously detected in direct-detection experiments, the challenge will be to determine what one may infer from the data. Here, I examine the prospects for reconstructing the local speed distribution of WIMPs in addition to WIMP particle-physics properties (mass, cross sections) from next-generation cryogenic and liquid-noble direct-detection experiments. I find that the common method of fixing the form of the velocity distribution when estimating constraints on WIMP mass and cross sections means losing out on the information on the speed distribution contained in the data and may lead to biases in the inferred values of the particle-physics parameters. I show that using a more general, empirical form of the speed distribution can lead to good constraints on the speed distribution. Moreover, one can use Bayesian model-selection criteria to determine if a theoretically-inspired functional form for the speed distribution (such as a Maxwell-Boltzmann distribution) fits better than an empirical model. The shape of the degeneracy between WIMP mass and cross sections and their offset from the true values of those parameters depends on the hypothesis for the speed distribution, which has significant implications for consistency checks between direct-detection and collider data. In addition, I find that the uncertainties on theoretical parameters depends sensitively on the upper end of the energy range used for WIMP searches. Better constraints on the WIMP particle-physics parameters and speed distribution are obtained if the WIMP search is extended to higher energy (~ 1 MeV).

preprint2011arXiv

Infall Times for Milky Way Satellites From Their Present-Day Kinematics

We analyze subhalos in the Via Lactea II (VL2) cosmological simulation to look for correlations among their infall times and z = 0 dynamical properties. We find that the present day orbital energy is tightly correlated with the time at which subhalos last crossed into the virial radius. This energy-infall correlation provides a means to infer infall times for Milky Way satellite galaxies. Assuming that the Milky Way's assembly can be modeled by VL2, we show that the infall times of some satellites are well constrained given only their Galactocentric positions and line-of-sight velocities. The constraints sharpen for satellites with proper motion measurements. We find that Carina, Ursa Minor, and Sculptor were all accreted early, more than 8 Gyr ago. Five other dwarfs, including Sextans and Segue 1, are also probable early accreters, though with larger uncertainties. On the other extreme, Leo T is just falling into the Milky Way for the first time while Leo I fell in \sim 2 Gyr ago and is now climbing out of the Milky Way's potential after its first perigalacticon. The energies of several other dwarfs, including Fornax and Hercules, point to intermediate infall times, 2 - 8 Gyr ago. We compare our infall time estimates to published star formation histories and find hints of a dichotomy between ultrafaint and classical dwarfs. The classical dwarfs appear to have quenched star formation after infall but the ultrafaint dwarfs tend to be quenched long before infall, at least for the cases in which our uncertainties allow us to discern differences. Our analysis suggests that the Large Magellanic Cloud crossed inside the Milky Way virial radius recently, within the last \sim 4 billion years.

preprint2010arXiv

Angular Momentum Evolution in Dark Matter Halos

We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.

preprint2010arXiv

Comments on recent work on dark-matter capture in the Solar System

Recently, several papers have appeared that examine the process of capturing dark-matter particles from the Galactic halo to orbits bound to the Solar System. The authors of these papers predict large enhancements to the local dark-matter density via gravitational three-body interactions with planets. However, these conclusions are wrong; these papers do not include the inverse process to capture, namely the ejection of dark-matter particles by three-body gravitational encounters. We emphasize previous work that shows that by including both capture and ejection of dark matter from the Solar System, the density of dark matter bound to the Solar System is small compared to the local Galactic dark-matter density.

preprint2010arXiv

Constraints on Decaying Dark Matter

We explore a dark-matter model in which there are two dark-matter species nearly degenerate in mass, with epsilon = Delta M/M << 1. The heavier particle undergoes two-body decay with a half-life tau, to the lighter dark-matter particle and a noninteracting massless particle. Unlike previous work on decaying dark matter, we explore the regime tau > 100 Myr and non-relativistic kick speeds vk / c = epsilon. Using a set of N-body simulations of isolated dark-matter halos, we show how halos change as a function of tau and vk. We find that tau < 40 Gyr is ruled out for vk > 20 km s^{-1} (epsilon > 10^{-4}) when we compare the simulations to observations of dwarf-galaxy- to cluster-mass dark matter halos. We highlight which set of observations should provide better future constraints for decays and other types of dark-matter physics.

preprint2010arXiv

Dark-matter decays and Milky Way satellite galaxies

We consider constraints on a phenomenological dark-matter model consisting of two nearly degenerate particle species using observed properties of the Milky Way satellite galaxy population. The two parameters of this model, assuming the particle masses are >~ GeV, are v_k, the recoil speed of the daughter particle, and tau, the lifetime of the parent particle. The satellite constraint that spans the widest range of v_k is the number of satellites that have a mass within 300 pc M300 > 5 x 10^6 solar masses, although constraints based on M300 in the classical dwarfs and the overall velocity function are competitive for v_k >~ 50 km/s. In general, we find that tau <~ 30 Gyr is ruled out for 20 km/s <~ v_k <~ 200 km/s, although we find that the limits on tau for fixed v_k can change constraints by a factor of ~3 depending on the star-formation histories of the satellites. We advocate using the distribution of M300 in Milky Way satellites determined by next-generation all-sky surveys and follow-up spectroscopy as a probe of dark-matter properties.

preprint2010arXiv

Dark-Matter Decays and Self-Gravitating Halos

We consider models in which a dark-matter particle decays to a slightly less massive daughter particle and a noninteracting massless particle. The decay gives the daughter particle a small velocity kick. Self-gravitating dark-matter halos that have a virial velocity smaller than this velocity kick may be disrupted by these particle decays, while those with larger virial velocities will be heated. We use numerical simulations to follow the detailed evolution of the total mass and density profile of self-gravitating systems composed of particles that undergo such velocity kicks as a function of the kick speed (relative to the virial velocity) and the decay time (relative to the dynamical time). We show how these decays will affect the halo mass-concentration relation and mass function. Using measurements of the halo mass-concentration relation and galaxy-cluster mass function to constrain the lifetime--kick-velocity parameter space for decaying dark matter, we find roughly that the observations rule out the combination of kick velocities greater than 100 km/s and decay times less than a few times the age of the Universe.

preprint2010arXiv

Getting the astrophysics and particle physics of dark matter out of next-generation direct detection experiments

The next decade will bring massive new data sets from experiments of the direct detection of weakly interacting massive particle (WIMP) dark matter. The primary goal of these experiments is to identify and characterize the dark-matter particle species. However, mapping the data sets to the particle-physics properties of dark matter is complicated not only by the considerable uncertainties in the dark-matter model, but by its poorly constrained local distribution function (the "astrophysics" of dark matter). In this Letter, I propose a shift in how to do direct-detection data analysis. I show that by treating the astrophysical and particle physics uncertainties of dark matter on equal footing, and by incorporating a combination of data sets into the analysis, one may recover both the particle physics and astrophysics of dark matter. Not only does such an approach yield more accurate estimates of dark-matter properties, but may illuminate how dark matter coevolves with galaxies.

preprint2010arXiv

Mapping the allowed parameter space for decaying dark matter models

I consider constraints on a phenomenological decaying-dark-matter model, in which two weakly-interacting massive particle (WIMP) species have a small mass splitting, and in which the heavier particle decays to the lighter particle and a massless particle on cosmological timescales. The decay parameter space is parameterized by $v_k$, the speed of the lighter particle in the center-of-mass frame of the heavier particle prior to decay, and the decay time $τ$. Since I consider the case in which dark-matter halos have formed before there has been significant decay, I focus on the effects of decay in already-formed halos. I show that the $v_k-τ$ parameter space may be constrained by observed properties of dark-matter halos. I highlight which set of observations is likely to yield the cleanest constraints on $v_k-τ$ parameter space, and calculate the constraints in those cases in which the effect of decay on the observables can be calculated without N-body simulations of decaying dark matter. I show that for $v_k \gtrsim 5\times 10^3$ km s$^{-1}$, the z=0 galaxy cluster mass function and halo mass-concentration relation constrain $τ\gtrsim$ 40 Gyr, and that precise constraints on $τ$ for smaller $v_k$ will require N-body simulations.

preprint2009arXiv

Dark Matter Disc Enhanced Neutrino Fluxes from the Sun and Earth

As disc galaxies form in a hierarchical cosmology, massive merging satellites are preferentially dragged towards the disc plane. The material accreted from these satellites forms a dark matter disc that contributes 0.25 - 1.5 times the non-rotating halo density at the solar position. Here, we show the importance of the dark disc for indirect dark matter detection in neutrino telescopes. Previous predictions of the neutrino flux from WIMP annihilation in the Earth and the Sun have assumed that Galactic dark matter is spherically distributed with a Gaussian velocity distribution, the standard halo model. Although the dark disc has a local density comparable to the dark halo, its higher phase space density at low velocities greatly enhances capture rates in the Sun and Earth. For typical dark disc properties, the resulting muon flux from the Earth is increased by three orders of magnitude over the SHM, while for the Sun the increase is an order of magnitude. This significantly increases the sensitivity of neutrino telescopes to fix or constrain parameters in WIMP models. The flux from the Earth is extremely sensitive to the detailed properties of the dark disc, while the flux from the Sun is more robust. The enhancement of the muon flux from the dark disc puts the search for WIMP annihilation in the Earth on the same level as the Sun for WIMP masses < 100 GeV.

preprint2008arXiv

Dynamics of WIMPs in the solar system and implications for direct and indirect detection

Semi-analytic treatments of the evolution of orbits of weakly interacting massive particles (WIMPs) in the solar system suggest that the WIMPs bound to the solar system may enhance the direct detection rate relative to that of the unbound population by up to a factor of order unity, and boost the flux of neutrinos from WIMP annihilation in the Earth by up to two orders of magnitude. To test these important but uncertain results, we perform a suite of numerical orbit integrations to explore the properties of the bound WIMP population as a function of the WIMP mass and the scattering cross section with baryonic matter. For regions of WIMP parameter space presently allowed by experiments, we find that (i) the bound WIMP population enhances the direct detection rate by at most ~1% relative to the rate from unbound halo WIMPs; (ii) it is unlikely that planned km^3-scale neutrino telescopes will detect neutrinos from WIMP annihilation in the Earth; (iii) the event rate from neutrinos produced by WIMP annihilation in the Sun may be much smaller than implied by the usual calculations, which assume that WIMPs scattered onto bound orbits are rapidly thermalized in the Sun.