Source author record

Manoj Kaplinghat

Manoj Kaplinghat 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

53works
9topics
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

53 published item(s)

preprint2023arXiv

Orbital Evolution of Satellite Galaxies in Self-Interacting Dark Matter Models

Dark matter self interactions can leave distinctive signatures on the properties of satellite galaxies around Milky Way--like hosts through their impact on tidal stripping, ram pressure, and gravothermal collapse. We delineate the regions of self-interacting dark matter parameter space -- specified by interaction cross section and a velocity scale -- where each of these effects dominates, and show how the relative mass loss depends on the satellite's initial mass, density profile and orbit. We obtain novel, conservative constraints in this parameter space using Milky Way satellite galaxies with notably high central densities and small pericenter distances. Our results for self-interacting dark matter models, in combination with constraints from clusters of galaxies, favor either velocity-dependent cross sections that lead to gravothermal core collapse in the densest satellites, or small cross sections which more closely resemble cold and collisionless dark matter.

preprint2022arXiv

Comparing Implementations of Self-Interacting Dark Matter in the Gizmo and Arepo Codes

Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes: Gizmo and Arepo. We begin with identical initial conditions for an isolated $10^{10}$ M$_\odot$ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in Gizmo and Arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2 g^{-1}$. Our tests are restricted to the core expansion phase where the core density decreases and core radius increases with time. We find better than 30% agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10% outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2 g^{-1}$ but finer distinctions would require further investigation.

preprint2021arXiv

Strong constraints on thermal relic dark matter from Fermi-LAT observations of the Galactic Center

The extended excess toward the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that, when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of three-dimensional inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of approximately TeV and rule out the thermal relic cross section up to approximately 300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties.

preprint2020arXiv

Accelerated core collapse in tidally stripped self-interacting dark matter halos

We use a semianalytic approach that is calibrated to N-body simulations to study the evolution of self-interacting dark matter cores in galaxies. We demarcate the regime where the temporal evolution of the core density follows a well-defined track set by the initial halo parameters and the cross section. Along this track, the central density reaches a minimum value set by the initial halo density. Further evolution leads to an outward heat transfer, inducing gravothermal core collapse such that the core shrinks as its density increases. We show that the time scale for the core collapse is highly sensitive to the outer radial density profile. Satellite galaxies with significant mass loss due to tidal stripping should have larger central densities and significantly faster core collapse compared to isolated halos. Such a scenario could explain the dense and compact cores of dwarf galaxies in the Local Group like Tucana (isolated from the Milky Way), the classical Milky Way satellite Draco, and some of the ultrafaint satellites. If the ultimate fate of core collapse is black hole formation, then the accelerated time scale provides a new mechanism for creating intermediate mass black holes.

preprint2020arXiv

Implications of a transition in the dark energy equation of state for the $H_0$ and $σ_8$ tensions

We explore the implications of a rapid appearance of dark energy between the redshifts ($z$) of one and two on the expansion rate and growth of perturbations. Using both Gaussian process regression and a parameteric model, we show that this is the preferred solution to the current set of low-redshift ($z<3$) distance measurements if $H_0=73~\rm km\,s^{-1}\,Mpc^{-1}$ to within 1\% and the high-redshift expansion history is unchanged from the $Λ$CDM inference by the Planck satellite. Dark energy was effectively non-existent around $z=2$, but its density is close to the $Λ$CDM model value today, with an equation of state greater than $-1$ at $z<0.5$. If sources of clustering other than matter are negligible, we show that this expansion history leads to slower growth of perturbations at $z<1$, compared to $Λ$CDM, that is measurable by upcoming surveys and can alleviate the $σ_8$ tension between the Planck CMB temperature and low-redshift probes of the large-scale structure.

preprint2020arXiv

Multiple Chemodynamic Stellar Populations of the Ursa Minor Dwarf Spheroidal Galaxy

We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $σ_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($σ_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($ε=0.75\pm0.03$) and the metal-poor population is closer to spherical ($ε=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

preprint2019arXiv

Dark Matter Cores and Cusps in Spiral Galaxies and their Explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

preprint2015arXiv

Dark Matter Halos as Particle Colliders: A Unified Solution to Small-Scale Structure Puzzles from Dwarfs to Clusters

Astrophysical observations spanning dwarf galaxies to galaxy clusters indicate that dark matter (DM) halos are less dense in their central regions compared to expectations from collisionless DM N-body simulations. Using detailed fits to DM halos of galaxies and clusters, we show that self-interacting DM (SIDM) may provide a consistent solution to the DM deficit problem across all scales, even though individual systems exhibit a wide diversity in halo properties. Since the characteristic velocity of DM particles varies across these systems, we are able to measure the self-interaction cross section as a function of kinetic energy and thereby deduce the SIDM particle physics model parameters. Our results prefer a mildly velocity-dependent cross section, from $σ/m \simeq 2\; {\rm cm^2/g}$ on galaxy scales to $σ/m \simeq 0.1\; {\rm cm^2/g}$ on cluster scales, consistent with the upper limits from merging clusters. Our results dramatically improve the constraints on SIDM models and may allow the masses of both DM and dark mediator particles to be measured even if the dark sector is completely hidden from the Standard Model, which we illustrate for the dark photon model.

preprint2015arXiv

Dark Matter in the Coming Decade: Complementary Paths to Discovery and Beyond

In this report we summarize the many dark matter searches currently being pursued through four complementary approaches: direct detection, indirect detection, collider experiments, and astrophysical probes. The essential features of broad classes of experiments are described, each with their own strengths and weaknesses. The complementarity of the different dark matter searches is discussed qualitatively and illustrated quantitatively in two simple theoretical frameworks. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program drawing from all four approaches.

preprint2015arXiv

Direct Detection Signatures of Self-Interacting Dark Matter with a Light Mediator

Self-interacting dark matter (SIDM) is a simple and well-motivated scenario that could explain long-standing puzzles in structure formation on small scales. If the required self-interaction arises through a light mediator (with mass $\sim 10$ MeV) in the dark sector, this new particle must be unstable to avoid overclosing the universe. The decay of the light mediator could happen due to a weak coupling of the hidden and visible sectors, providing new signatures for direct detection experiments. The SIDM nuclear recoil spectrum is more peaked towards low energies compared to the usual case of contact interactions, because the mediator mass is comparable to the momentum transfer of nuclear recoils. We show that the SIDM signal could be distinguished from that of DM particles with contact interactions by considering the time-average energy spectrum in experiments employing different target materials, or the average and modulated spectra in a single experiment. Using current limits from LUX and SuperCDMS, we also derive strong bounds on the mixing parameter between hidden and visible sector.

preprint2015arXiv

Discovery of a New Galactic Center Excess Consistent with Upscattered Starlight

We present a new extended gamma ray excess detected with the Fermi Satellite Large Area Telescope toward the Galactic Center that traces the morphology of infrared starlight emission. Combined with its measured spectrum, this new extended source is approximately consistent with inverse Compton emission from a high-energy electron-positron population with energies up to about 10 GeV. Previously detected emissions tracing the 20 cm radio, interpreted as bremsstrahlung radiation, and the Galactic Center Extended emission tracing a spherical distribution and peaking at 2 GeV, are also detected. We show that the inverse Compton and bremsstrahlung emissions are likely due to the same source of electrons and positrons. All three extended emissions may be explained within the framework of a model where the dark matter annihilates to leptons or a model with unresolved millisecond pulsars in the Galactic Center.

preprint2015arXiv

Galactic Center Excess in Gamma Rays from Annihilation of Self-Interacting Dark Matter

Observations by the Fermi-LAT telescope have uncovered a significant $γ$-ray excess toward the Milky Way Galactic Center. There has been no detection of a similar signal in the direction of the Milky Way dwarf spheroidal galaxies. Additionally, astronomical observations indicate that dwarf galaxies and other faint galaxies are less dense than predicted by the simplest cold dark matter models. We show that a self-interacting dark matter model with a particle mass of roughly 50 GeV annihilating to the mediator responsible for the strong self-interaction can simultaneously explain all three observations. The mediator is necessarily unstable and its mass must be below about 100 MeV in order to lower densities in faint galaxies. If the mediator decays to electron-positron pairs with a cross section on the order of the thermal relic value, then we find that these pairs can up-scatter the interstellar radiation field and produce the observed $γ$-ray excess. We show that this model is compatible with all current constraints and highlight detectable signatures unique to self-interacting dark matter models.

preprint2015arXiv

Inflation that runs naturally: Gravitational waves and suppression of power at large and small scales

We point out three correlated predictions of the axion monodromy inflation model: large amplitude of gravitational waves, suppression of power on horizon scales and on scales relevant for the formation of dwarf galaxies. While these predictions are likely generic to models with oscillations in the inflaton potential, the axion monodromy model naturally accommodates the required running spectral index through Planck-scale corrections to the inflaton potential. Applying this model to a combined data set of Planck, ACT, SPT, and WMAP low-$\ell$ polarization cosmic microwave background (CMB) data, we find a best-fit tensor-to-scalar ratio $r_{0.05} = 0.07^{+0.05}_{-0.04}$ due to gravitational waves, which may have been observed by the BICEP2 experiment. Despite the contribution of gravitational waves, the total power on large scales (CMB power spectrum at low multipoles) is lower than the standard $Λ$CDM cosmology with a power-law spectrum of initial perturbations and no gravitational waves, thus mitigating some of the tension on large scales. There is also a reduction in the matter power spectrum of 20-30\% at scales corresponding to $k = 10~{\rm Mpc}^{-1}$, which are relevant for dwarf galaxy formation. This will alleviate some of the unsolved small-scale structure problems in the standard $Λ$CDM cosmology. The inferred matter power spectrum is also found to be consistent with recent Lyman-$α$ forest data, which is in tension with the Planck-favored $Λ$CDM model with power-law primordial power spectrum.

preprint2015arXiv

Self-Interacting Dark Matter from a Non-Abelian Hidden Sector

There is strong evidence in favor of the idea that dark matter is self interacting, with the cross section-to-mass ratio $σ/ m \sim 1\,\mathrm{cm^2/g} \sim 1\,\mathrm{barn/GeV}$. We show that viable models of dark matter with this large cross section are straightforwardly realized with non-Abelian hidden sectors. In the simplest of such models, the hidden sector is a pure gauge theory, and the dark matter is composed of hidden glueballs with a mass around $100\,\mathrm{MeV}$. Alternatively, the hidden sector may be a supersymmetric pure gauge theory with a $\sim 10\,\mathrm{TeV}$ gluino thermal relic. In this case, the dark matter is largely composed of glueballinos that strongly self interact through the exchange of light glueballs. We present a unified framework that realizes both of these possibilities in anomaly-mediated supersymmetry breaking, where, depending on a few model parameters, the dark matter may be composed of hidden glueballinos, hidden glueballs, or a mixture of the two. These models provide simple examples of multicomponent dark matter, have interesting implications for particle physics and cosmology, and include cases where a subdominant component of dark matter may be extremely strongly self interacting, with interesting astrophysical consequences.

preprint2015arXiv

SIMPle Dark Matter: Self-Interactions and keV Lines

We consider a simple supersymmetric hidden sector: pure SU(N) gauge theory. Dark matter is made up of hidden glueballinos with mass $m_X$ and hidden glueballs with mass near the confinement scale $Λ$. For $m_X \sim 1\,\text{TeV}$ and $Λ\sim 100\,\text{MeV}$, the glueballinos freeze out with the correct relic density and self-interact through glueball exchange to resolve small-scale structure puzzles. An immediate consequence is that the glueballino spectrum has a hyperfine splitting of order $Λ^2 / m_X \sim 10\,\text{keV}$. We show that the radiative decays of the excited state can explain the observed 3.5 keV X-ray line signal from clusters of galaxies, Andromeda, and the Milky Way.

preprint2014arXiv

Astrophysical and Dark Matter Interpretations of Extended Gamma-Ray Emission from the Galactic Center

We construct empirical models of the diffuse gamma-ray background toward the Galactic Center. Including all known point sources and a template of emission associated with interactions of cosmic rays with molecular gas, we show that the extended emission observed previously in the Fermi Large Area Telescope data toward the Galactic Center is detected at high significance for all permutations of the diffuse model components. However, we find that the fluxes and spectra of the sources in our model change significantly depending on the background model. In particular, the spectrum of the central Sgr A$^\ast$ source is less steep than in previous works and the recovered spectrum of the extended emission has large systematic uncertainties, especially at lower energies. If the extended emission is interpreted to be due to dark matter annihilation, we find annihilation into pure $b$-quark and $τ$-lepton channels to be statistically equivalent goodness of fits. In the case of the pure $b$-quark channel, we find a dark matter mass of $39.4\left(^{+3.7}_{-2.9}\rm\ stat.\right)\left(\pm 7.9\rm\ sys.\right)\rm\ GeV$, while a pure $τ^{+} τ^{-}$-channel case has an estimated dark matter mass of $9.43\left(^{+0.63}_{-0.52}\rm\ stat.\right)(\pm 1.2\rm\ sys.)\ GeV$. Alternatively, if the extended emission is interpreted to be astrophysical in origin such as due to unresolved millisecond pulsars, we obtain strong bounds on dark matter annihilation, although systematic uncertainties due to the dependence on the background models are significant.

preprint2014arXiv

Running with BICEP2: Implications for Small-Scale Problems in CDM

The BICEP2 results, when interpreted as a gravitational wave signal and combined with other CMB data, suggest a roll-off in power towards small scales in the primordial matter power spectrum. Among the simplest possibilities is a running of the spectral index. Here we show that the preferred level of running alleviates small-scale issues within the $Λ$CDM model, more so even than viable WDM models. We use cosmological zoom-in simulations of a Milky Way-size halo along with full-box simulations to compare predictions among four separate cosmologies: a BICEP2-inspired running index model ($α_s$ = -0.024), two fixed-tilt $Λ$CDM models motivated by Planck, and a 2.6 keV thermal WDM model. We find that the running BICEP2 model reduces the central densities of large dwarf-size halos ($V_\mathrm{max}$ ~ 30 - 80 km s$^{-1}$) and alleviates the too-big-to-fail problem significantly compared to our adopted Planck and WDM cases. Further, the BICEP2 model suppresses the count of small subhalos by ~50% relative to Planck models, and yields a significantly lower "boost" factor for dark matter annihilation signals. Our findings highlight the need to understand the shape of the primordial power spectrum in order to correctly interpret small-scale data.

preprint2014arXiv

The dynamics of isolated Local Group galaxies

We measured velocities of 862 individual red giant stars in seven isolated dwarf galaxies in the Local Group: NGC 6822, IC 1613, VV 124 (UGC 4879), the Pegasus dwarf irregular galaxy (DDO 216), Leo A, Cetus, and Aquarius (DDO 210). We also computed velocity dispersions, taking into account the measurement uncertainties on individual stars. None of the isolated galaxies is denser than the densest Local Group satellite galaxy. Furthermore, the isolated dwarf galaxies have no obvious distinction in the velocity dispersion--half-light radius plane from the satellite galaxies of the Milky Way and M31. The similarity of the isolated and satellite galaxies' dynamics and structural parameters imposes limitations on environmental solutions to the too-big-to-fail problem, wherein there are fewer dense dwarf satellite galaxies than would be expected from cold dark matter simulations. This data set also has many other applications for dwarf galaxy evolution, including the transformation of dwarf irregular into dwarf spheroidal galaxies. We intend to explore these issues in future work.

preprint2014arXiv

Tying Dark Matter to Baryons with Self-interactions

Self-interacting dark matter (SIDM) models have been proposed to solve the small-scale issues with the collisionless cold dark matter (CDM) paradigm. We derive equilibrium solutions in these SIDM models for the dark matter halo density profile including the gravitational potential of both baryons and dark matter. Self-interactions drive dark matter to be isothermal and this ties the core sizes and shapes of dark matter halos to the spatial distribution of the stars, a radical departure from previous expectations and from CDM predictions. Compared to predictions of SIDM-only simulations, the core sizes are smaller and the core densities are higher, with the largest effects in baryon-dominated galaxies. As an example, we find a core size around 0.5 kpc for dark matter in the Milky Way, more than an order of magnitude smaller than the core size from SIDM-only simulations, which has important implications for indirect searches of SIDM candidates.

preprint2013arXiv

Are Light Sterile Neutrinos Preferred or Disfavored by Cosmology?

We find that the viability of a cosmological model that incorporates 2 sterile neutrinos with masses around 1 eV each, as favored by global neutrino oscillation analyses including short baseline results, is significantly dependent on the choice of datasets included in the analysis and the ability to control the systematic uncertainties associated with these datasets. Our analysis includes a variety of cosmological probes including the cosmic microwave background (WMAP7+SPT), Hubble constant (HST), galaxy power spectrum (SDSS-DR7), and supernova distances (SDSS and Union2 compilations). In the joint observational analysis, our sterile neutrino model is equally favored as a LCDM model when using the MLCS light curve fitter for the supernova measurements, and strongly disfavored by the data at Δχ^2 ~ 18 when using the SALT2 fitter. When excluding the supernova measurements, the sterile neutrino model is disfavored by the other datasets at Δχ^2 ~ 12, and at best becomes mildly disfavored at Δχ^2 ~ 3 when allowing for curvature, evolving dark energy, additional relativistic species, running of the spectral index, and freedom in the primordial helium abundance. No single additional parameter accounts for most of this effect. Therefore, if laboratory experiments continue to favor a scenario with roughly eV mass sterile neutrinos, and if this becomes decisively disfavored by cosmology, then a more exotic cosmological model than explored here may become necessary.

preprint2013arXiv

Dark Matter in the Coming Decade: Complementary Paths to Discovery and Beyond

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

preprint2013arXiv

Detection of a Gamma-Ray Source in the Galactic Center Consistent with Extended Emission from Dark Matter Annihilation and Concentrated Astrophysical Emission

We show the existence of a statistically significant, robust detection of a gamma-ray source in the Milky Way Galactic Center that is consistent with a spatially extended signal using about 4 years of Fermi-LAT data. The gamma-ray flux is consistent with annihilation of dark matter particles with a thermal annihilation cross-section if the spatial distribution of dark matter particles is similar to the predictions of dark matter only simulations. We find statistically significant detections of an extended source with gamma-ray spectrum that is consistent with dark matter particle masses of approximately 10 GeV to 1 TeV annihilating to b/b-bar quarks, and masses approximately 10 GeV to 30 GeV annihilating to tau+ tau- leptons. However, a part of the allowed region in this interpretation is in conflict with constraints from Fermi observations of the Milky Way satellites. The biggest improvement over the fit including just the point sources is obtained for a 30 GeV dark matter particle annihilating to b/b-bar quarks. The gamma-ray intensity and spectrum are also well fit with emission from a millisecond pulsar (MSP) population following a density profile like that of low-mass X-ray binaries observed in M31. The greatest goodness-of-fit of the extended emission is with spectra consistent with known astrophysical sources like MSPs in globular clusters or cosmic ray bremsstrahlung on molecular gas. Therefore, we conclude that the bulk of the emission is likely from an unresolved or spatially extended astrophysical source. However, the interesting possibility of all or part of the extended emission being from dark matter annihilation cannot be excluded at present.

preprint2013arXiv

Direct Detection Portals for Self-interacting Dark Matter

Dark matter self-interactions can affect the small scale structure of the Universe, reducing the central densities of dwarfs and low surface brightness galaxies in accord with observations. From a particle physics point of view, this points toward the existence of a 1-100 MeV particle in the dark sector that mediates self-interactions. Since mediator particles will generically couple to the Standard Model, direct detection experiments provide sensitive probes of self-interacting dark matter. We consider three minimal mechanisms for coupling the dark and visible sectors: photon kinetic mixing, Z boson mass mixing, and the Higgs portal. Self-interacting dark matter motivates a new benchmark paradigm for direct detection via momentum-dependent interactions, and ton-scale experiments will cover astrophysically motivated parameter regimes that are unconstrained by current limits. Direct detection is a complementary avenue to constrain velocity-dependent self-interactions that evade astrophysical bounds from larger scales, such as those from the Bullet Cluster.

preprint2013arXiv

Segue 2: The Least Massive Galaxy

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

preprint2013arXiv

Self-interacting Dark Matter Benchmarks

Dark matter self-interactions have important implications for the distributions of dark matter in the Universe, from dwarf galaxies to galaxy clusters. We present benchmark models that illustrate characteristic features of dark matter that is self-interacting through a new light mediator. These models have self-interactions large enough to change dark matter densities in the centers of galaxies in accord with observations, while remaining compatible with large-scale structure data and all astrophysical observations such as halo shapes and the Bullet Cluster. These observations favor a mediator mass in the 10 - 100 MeV range and large regions of this parameter space are accessible to direct detection experiments like LUX, SuperCDMS, and XENON1T.

preprint2013arXiv

Sterile neutrino dark matter bounds from galaxies of the Local Group

We show that the canonical oscillation-based (non-resonant) production of sterile neutrino dark matter is inconsistent at $>99$% confidence with observations of galaxies in the Local Group. We set lower limits on the non-resonant sterile neutrino mass of $2.5$ keV (equivalent to $0.7$ keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way, as well as limits of $8.8$ keV (equivalent to $1.8$ keV thermal mass) based on subhalo counts of $N$-body simulations of M 31 analogues. Combined with improved upper mass limits derived from significantly deeper X-ray data of M 31 with full consideration for background variations, we show that there remains little room for non-resonant production if sterile neutrinos are to explain $100$% of the dark matter abundance. Resonant and non-oscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos.

preprint2013arXiv

Whitepaper on Super-weakly Interacting Massive Particles for Snowmass 2013

Super-weakly interacting massive particles produced in the late decays of weakly interacting massive particles (WIMPs) are generic in large regions of supersymmetric parameter space and other frameworks for physics beyond the standard model. If their masses are similar to that of the decaying WIMP, then they could naturally account for all of the cosmological dark matter abundance. Their astrophysical consequences and collider signatures are distinct and different from WIMP candidates. In particular, they could modify Big Bang Nucleosynthesis, distort the Cosmic Microwave Background, reduce galactic substructure and lower central densities of low-mass galaxies.

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 Energy and Neutrino Masses from Future Measurements of the Expansion History and Growth of Structure

We forecast the expected cosmological constraints from a combination of probes of both the universal expansion rate and matter perturbation growth, in the form of weak lensing tomography, galaxy tomography, supernovae, and the cosmic microwave background incorporating all cross-correlations between the observables for an extensive cosmological parameter set. We allow for non-zero curvature and parameterize our ignorance of the early universe by allowing for a non-negligible fraction of dark energy (DE) at high redshifts. We find that early DE density can be constrained to 0.2% of the critical density of the universe with Planck combined with a ground-based LSST-like survey, while curvature can be constrained to 0.06%. However, these additional degrees of freedom degrade our ability to measure late-time dark energy and the sum of neutrino masses. We find that the combination of cosmological probes can break degeneracies and constrain the sum of neutrino masses to 0.04 eV, present DE density also to 0.2% of the critical density, and the equation of state to 0.01 - roughly a factor of two degradation in the constraints overall compared to the case without allowing for early DE. The constraints for a space-based mission are similar. Even a modest 1% dark energy fraction of the critical density at high redshift, if not accounted for in future analyses, biases the cosmological parameters by up to 2 sigma. Our analysis suggests that throwing out nonlinear scales (multipoles > 1000) may not result in significant degradation in future parameter measurements when multiple cosmological probes are combined. We find that including cross-correlations between the different probes can result in improved constraints by up to a factor of 2 for the sum of neutrino masses and early dark energy density.

preprint2012arXiv

Evidence for Substructure in Ursa Minor Dwarf Spheroidal Galaxy using a Bayesian Object Detection Method

We present a method for identifying localized secondary populations in stellar velocity data using Bayesian statistical techniques. We apply this method to the dwarf spheroidal galaxy Ursa Minor and find two secondary objects in this satellite of the Milky Way. One object is kinematically cold with a velocity dispersion of $4.25 \pm 0.75\ \kms$ and centered at $(9.1\arcmin \pm 1.5, 7.2\arcmin \pm 1.2)$ in relative RA and DEC with respect to the center of Ursa Minor. The second object has a large velocity offset of $-12.8^{+1.75}_{-1.5}\ \kms$ compared to Ursa Minor and centered at $(-14.0\arcmin^{+2.4}_{-5.8}, -2.5\arcmin^{+0.4}_{-1.0})$. The kinematically cold object has been found before using a smaller data set but the prediction that this cold object has a velocity dispersion larger than $2.0\ \kms$ at 95% C.L. differs from previous work. We use two and three component models along with the information criteria and Bayesian evidence model selection methods to argue that Ursa Minor has one or two localized secondary populations. The significant probability for a large velocity dispersion in each secondary object raises the intriguing possibility that each has its own dark matter halo, that is, it is a satellite of a satellite of the Milky Way.

preprint2012arXiv

The Milky Way's bright satellites as an apparent failure of LCDM

We use the Aquarius simulations to show that the most massive subhalos in galaxy-mass dark matter halos in LCDM are grossly inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies. While the best-fitting hosts of the dwarf spheroidals all have 12 < Vmax < 25 km/s, LCDM simulations predict at least ten subhalos with Vmax > 25 km/s. These subhalos are also among the most massive at earlier times, and significantly exceed the UV suppression mass back to z ~ 10. No LCDM-based model of the satellite population of the Milky Way explains this result. The problem lies in the satellites' densities: it is straightforward to match the observed Milky Way luminosity function, but doing so requires the dwarf spheroidals to have dark matter halos that are a factor of ~5 more massive than is observed. Independent of the difficulty in explaining the absence of these dense, massive subhalos, there is a basic tension between the derived properties of the bright Milky Way dwarf spheroidals and LCDM expectations. The inferred infall masses of these galaxies are all approximately equal and are much lower than standard LCDM predictions for systems with their luminosities. Consequently, their implied star formation efficiencies span over two orders of magnitude, from 0.2% to 20% of baryons converted into stars, in stark contrast with expectations gleaned from more massive galaxies. We explore possible solutions to these problems within the context of LCDM and find them to be unconvincing. In particular, we use controlled simulations to demonstrate that the small stellar masses of the bright dwarf spheroidals make supernova feedback an unlikely explanation for their low inferred densities.

preprint2011arXiv

A Complete Spectroscopic Survey of the Milky Way satellite Segue 1: Dark matter content, stellar membership and binary properties from a Bayesian analysis

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.

preprint2011arXiv

A Complete Spectroscopic Survey of the Milky Way Satellite Segue 1: The Darkest Galaxy

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

preprint2011arXiv

Primordial non-Gaussianity from the 21 cm Power Spectrum during the Epoch of Reionization

Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact of non-Gaussianity on the ionization power spectrum from 21 cm emission during the epoch of reionization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the expected constraint on the local non-Gaussianity parameter f_NL for current and next-generation 21 cm experiments. We show that experiments such as SKA and MWA could measure f_NL values of order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope like Omniscope.

preprint2011arXiv

Too big to fail? The puzzling darkness of massive Milky Way subhaloes

We show that dissipationless LCDM simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (L_V > 10^5 L_sun). These dark subhaloes have circular velocities at infall of 30-70 km/s and infall masses of [0.2-4] x 10^10 M_sun. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3 - 1 kpc.

preprint2010arXiv

Accurate masses for dispersion-supported galaxies

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 σ_los^2 r_1/2 / G ~ 4 σ_los^2 R_eff / G, where σ_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

preprint2010arXiv

Correcting velocity dispersions of dwarf spheroidal galaxies for binary orbital motion

We show that measured velocity dispersions of dwarf spheroidal galaxies from about 4 to 10 km/s are unlikely to be inflated by more than 20% due to the orbital motion of binary stars, and demonstrate that the intrinsic velocity dispersions can be determined to within a few percent accuracy using two-epoch observations with 1-2 years as the optimal time interval. The crucial observable is the threshold fraction--the fraction of stars that show velocity changes larger than a given threshold between measurements. The threshold fraction is tightly correlated with the dispersion introduced by binaries, independent of the underlying binary fraction and distribution of orbital parameters. We outline a simple procedure to correct the velocity dispersion to within a few percent accuracy by using the threshold fraction and provide fitting functions for this method. We also develop a methodology for constraining properties of binary populations from both single- and two-epoch velocity measurements by including the binary velocity distribution in a Bayesian analysis.

preprint2010arXiv

Halo Shape and Relic Density Exclusions of Sommerfeld-Enhanced Dark Matter Explanations of Cosmic Ray Excesses

Dark matter with Sommerfeld-enhanced annihilation has been proposed to explain observed cosmic ray positron excesses in the 10 GeV to TeV energy range. We show that the required enhancement implies thermal relic densities that are too small to be all of dark matter. We also show that the dark matter is sufficiently self-interacting that observations of elliptical galactic dark matter halos exclude large Sommerfeld enhancement for light force carriers. Resonant Sommerfeld enhancement does not modify these conclusions, and the astrophysical boosts required to resolve these discrepancies are disfavored, especially when significant self-interactions suppress halo substructure.

preprint2010arXiv

Sommerfeld Enhancements for Thermal Relic Dark Matter

The annihilation cross section of thermal relic dark matter determines both its relic density and indirect detection signals. We determine how large indirect signals may be in scenarios with Sommerfeld-enhanced annihilation, subject to the constraint that the dark matter has the correct relic density. This work refines our previous analysis through detailed treatments of resonant Sommerfeld enhancement and the effect of Sommerfeld enhancement on freeze out. Sommerfeld enhancements raise many interesting issues in the freeze out calculation, and we find that the cutoff of resonant enhancement, the equilibration of force carriers, the temperature of kinetic decoupling, and the efficiency of self-interactions for preserving thermal velocity distributions all play a role. These effects may have striking consequences; for example, for resonantly-enhanced Sommerfeld annihilation, dark matter freezes out but may then chemically recouple, implying highly suppressed indirect signals, in contrast to naive expectations. In the minimal scenario with standard astrophysical assumptions, and tuning all parameters to maximize the signal, we find that, for force-carrier mass m_phi = 250 MeV and dark matter masses m_X = 0.1, 0.3, and 1 TeV, the maximal Sommerfeld enhancement factors are S_eff = 7, 30, and 90, respectively. Such boosts are too small to explain both the PAMELA and Fermi excesses. Non-minimal models may require smaller boosts, but the bounds on S_eff could also be more stringent, and dedicated freeze out analyses are required. For concreteness, we focus on 4 mu final states, but we also discuss 4 e and other modes, deviations from standard astrophysical assumptions and non-minimal particle physics models, and we outline the steps required to determine if such considerations may lead to a self-consistent explanation of the PAMELA or Fermi excesses.

preprint2010arXiv

Stealth Galaxies in the Halo of the Milky Way

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects that could number in the hundreds and have typical peak circular velocities of about 8 km/s and masses within 300 pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo.

preprint2010arXiv

The Case Against Warm or Self-Interacting Dark Matter as Explanations for Cores in Low Surface Brightness Galaxies

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. We fit these dark matter models to the data and determine the halo core radii and central densities. While the minimum core size in WDM models is predicted to decrease with halo mass, we find that the inferred core radii increase with halo mass and also cannot be explained with a single value of the primordial phase space density. Moreover, if the core size is set by WDM particle properties, then even the smallest cores we infer would require primordial phase space density values that are orders of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also find that the dark matter halo core densities vary by a factor of about 30 from system to system while showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues against the core size being directly set by large self-interactions (scattering or annihilation) of dark matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or SIDM over other dark matter models.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

Hidden Charged Dark Matter

Can dark matter be stabilized by charge conservation, just as the electron is in the standard model? We examine the possibility that dark matter is hidden, that is, neutral under all standard model gauge interactions, but charged under an exact U(1) gauge symmetry of the hidden sector. Such candidates are predicted in WIMPless models, supersymmetric models in which hidden dark matter has the desired thermal relic density for a wide range of masses. Hidden charged dark matter has many novel properties not shared by neutral dark matter: (1) bound state formation and Sommerfeld-enhanced annihilation after chemical freeze out may reduce its relic density, (2) similar effects greatly enhance dark matter annihilation in protohalos at redshifts of z ~ 30, (3) Compton scattering off hidden photons delays kinetic decoupling, suppressing small scale structure, and (4) Rutherford scattering makes such dark matter self-interacting and collisional, potentially impacting properties of the Bullet Cluster and the observed morphology of galactic halos. We analyze all of these effects in a WIMPless model in which the hidden sector is a simplified version of the minimal supersymmetric standard model and the dark matter is a hidden sector stau. We find that charged hidden dark matter is viable and consistent with the correct relic density for reasonable model parameters and dark matter masses in the range 1 GeV < m_X < 10 TeV. At the same time, in the preferred range of parameters, this model predicts cores in the dark matter halos of small galaxies and other halo properties that may be within the reach of future observations. These models therefore provide a viable and well-motivated framework for collisional dark matter with Sommerfeld enhancement, with novel implications for astrophysics and dark matter searches.

preprint2009arXiv

The Dark Disk of the Milky Way

Massive satellite accretions onto early galactic disks can lead to the deposition of dark matter in disk-like configurations that co-rotate with the galaxy. This phenomenon has potentially dramatic consequences for dark matter detection experiments. We utilize focused, high-resolution simulations of accretion events onto disks designed to be Galaxy analogues, and compare the resultant disks to the morphological and kinematic properties of the Milky Way's thick disk in order to bracket the range of co-rotating accreted dark matter. We find that the Milky Way's merger history must have been unusually quiescent compared to median LCDM expectations and therefore its dark disk must be relatively small: the fraction of accreted dark disk material near the Sun is about 20% of the host halo density or smaller and the co-rotating dark matter fraction near the Sun, defined as particles moving with a rotational velocity lag less than 50 km/s, is enhanced by about 30% or less compared to a standard halo model. Such a dark disk could contribute dominantly to the low energy (of order keV for a dark matter particle with mass 100 GeV) nuclear recoil event rate of direct dectection experiments, but it will not change the likelihood of detection significantly. These dark disks provide testable predictions of weakly-interacting massive particle dark matter models and should be considered in detailed comparisons to experimental data. Our findings suggest that the dark disk of the Milky Way may provide a detectable signal for indirect detection experiments, contributing up to about 25% of the dark matter self-annihilation signal in the direction of the center of the Galaxy, lending the signal a noticeably oblate morphology.

preprint2008arXiv

CMBPol Mission Concept Study: Gravitational Lensing

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

preprint2007arXiv

Redefining the Missing Satellites Problem

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

preprint2006arXiv

A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

preprint2004arXiv

The Supernova Relic Neutrino Backgrounds at KamLAND and Super-Kamiokande

We calculate the Supernova Relic Neutrino (SRN) background flux for the KamLAND and Super-Kamiokande (SK) detectors, motivated by the reduction in background at SK and new results for the star formation history (e.g., from the Sloan Digital Sky Survey (SDSS)). Our best estimate for the flux at SK is slightly below, but very close to the current SK upper limit. The SK upper limit is already inconsistent with a range of star formation histories allowed by the SDSS data. We estimate that the SRN background should be detected (at 1-sigma) at SK with a total of about 9 years (including the existing 4 years) of data. While KamLAND is a much smaller detector compared to SK, it profits from being practically background-free and from its sensitivity to the lower energy supernova neutrinos. KamLAND could make a 1-sigma detection of the SRN with a total of about 5 years of data. Given the small expected SRN event rate, we also consider the detection of the SRN in a modified SK detector with a lower threshold and reduced background where the time to detection can be reduced by a factor of 10 relative to the existing SK estimate.

preprint2004arXiv

Weak Lensing of the CMB: Sampling Errors on B-Modes

The B modes generated by the lensing of CMB polarization are a primary target for the upcoming generation of experiments and can potentially constrain quantities such as the neutrino mass and dark energy equation of state. The net sample variance on the small scale B modes out to l=2000 exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses that generate them. It manifests itself as highly correlated band powers with correlation coefficients approaching 70% for wide bands of Delta l/l \~0.25. It will double the total variance for experiments that achieve a sensitivity of approximately 4 uK-arcmin and a beam of several arcminutes or better. This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first detection.

preprint2003arXiv

Applications of High Resolution High Sensitivity Observations of the CMB

With WMAP putting the phenomenological standard model of cosmology on a strong footing, one can look forward to mining the cosmic microwave background (CMB) for fundamental physics with higher sensitivity and on smaller scales. Future CMB observations have the potential to measure absolute neutrino masses, test for cosmic acceleration independent of supernova Ia observations, probe for the presence of dark energy at redshifts of 2 and larger, illuminate the end of the dark ages, measure the scale--dependence of the primordial power spectrum and detect gravitational waves generated by inflation.

preprint2003arXiv

Determining Neutrino Mass from the CMB Alone

Distortions of Cosmic Microwave Background (CMB) temperature and polarization maps caused by gravitational lensing, observable with high angular resolution and high sensitivity, can be used to measure the neutrino mass. Assuming two massless species and one with mass m_nu we forecast sigma(m_nu) = 0.15 eV from the Planck satellite and sigma(m_nu)=0.04 eV from observations with twice the angular resolution and about 20 times the sensitivity. A detection is likely at this higher sensitivity since the observation of atmospheric neutrino oscillations require mass-squared differences greater than about (0.04 eV)^2.

preprint2002arXiv

Testing The Friedmann Equation: The Expansion of the Universe During Big-Bang Nucleosynthesis

In conventional general relativity, the expansion rate H of a Robertson-Walker universe is related to the energy density by the Friedmann equation. Aside from the present day, the only epoch at which we can constrain the expansion history in a model-independent way is during Big-Bang Nucleosynthesis (BBN). We consider a simple two-parameter characterization of the behavior of H during BBN and derive constraints on this parameter space, finding that the allowed region of parameter space is essentially one-dimensional. We also study the effects of a large neutrino asymmetry within this framework. Our results provide a simple way to compare an alternative cosmology to the observational requirement of matching the primordial abundances of the light elements.