We present the Model-Agnostic Dark Halo Analysis Tool (MADHAT), a numerical tool which implements a Fermi-LAT data-driven, model-independent analysis of gamma-ray emission from dwarf satellite galaxies and dwarf galaxy candidates due to dark matter annihilation, dark matter decay, or other nonstandard or unknown astrophysics. This tool efficiently provides statistical upper bounds on the number of observed photons in excess of the number expected, based on empirical determinations of foregrounds and backgrounds, using a stacked analysis of any selected set of dwarf targets. It also calculates the resulting bounds on the properties of dark matter under any assumptions the user makes regarding dark sector particle physics or astrophysics. As an application, we determine new bounds on Sommerfeld-enhanced dark matter annihilation in a set of eight dwarfs. MADHAT v1.0 includes 58 dwarfs and dwarf candidate targets, and we discuss future planned developments. MADHAT is available and will be maintained at https://github.com/MADHATdm

The DAMPE experiment has recently reported an electron spectrum that can be explained by dark matter annihilation into charged lepton pairs in a nearby dark matter clump. The accompanying bremsstrahlung may yield a gamma-ray excess with a known spectral shape that extends over an angular scale of $O(10^\circ)$. We show that such an excess is not present in Fermi-LAT data.

If dark matter interactions with Standard Model particles are $CP$-violating, then dark matter annihilation/decay can produce photons with a net circular polarization. We consider the prospects for experimentally detecting evidence for such a circular polarization. We identify optimal models for dark matter interactions with the Standard Model, from the point of view of detectability of the net polarization, for the case of either symmetric or asymmetric dark matter. We find that, for symmetric dark matter, evidence for net polarization could be found by a search of the Galactic Center by an instrument sensitive to circular polarization with an efficiency-weighted exposure of at least $50000~\text{cm}^2~\text{yr}$, provided the systematic detector uncertainties are constrained at the $1\%$ level. Better sensitivity can be obtained in the case of asymmetric dark matter. We discuss the prospects for achieving the needed level of performance using possible detector technologies.

We present a general, model-independent formalism for determining bounds on the production of photons in dwarf spheroidal galaxies via dark matter annihilation, applicable to any set of assumptions about dark matter particle physics or astrophysics. As an illustration, we analyze gamma-ray data from the Fermi Large Area Telescope to constrain a variety of nonstandard dark matter models, several of which have not previously been studied in the context of dwarf galaxy searches.

We consider the prospects for multiple dark matter direct detection experiments to determine if the interactions of a dark matter candidate are isospin-violating. We focus on theoretically well-motivated examples of isospin-violating dark matter (IVDM), including models in which dark matter interactions with nuclei are mediated by a dark photon, a Z, or a squark. We determine that the best prospects for distinguishing IVDM from the isospin-invariant scenario arise in the cases of dark photon- or Z-mediated interactions, and that the ideal experimental scenario would consist of large exposure xenon- and neon-based detectors. If such models just evade current direct detection limits, then one could distinguish such models from the standard isospin-invariant case with two detectors with of order 100 ton-year exposure.

If dark matter is composed of weakly interacting particles, Earth's orbital motion may induce a small annual variation in the rate at which these particles interact in a terrestrial detector. The DAMA collaboration has identified at a 9.3$\sigma$ confidence level such an annual modulation in their event rate over two detector iterations, DAMA/NaI and DAMA/LIBRA, each with $\sim7$ years of observations. We statistically examine the time dependence of the modulation amplitudes, which "by eye" appear to be decreasing with time in certain energy ranges. We perform a chi-squared goodness of fit test of the average modulation amplitudes measured\ by the two detector iterations which rejects the hypothesis of a consistent modulation amplitude at greater than 80\%, 96\%, and 99.6\% for the 2--4~keVee, 2--5~keVee and 2--6~keVee energy ranges, respectively. We also find that among the 14 annual cycles there are three $\gtrsim 3\sigma$ departures from the average in the 5-6~keVee energy range. In addition, we examined several phenomenological models for the time dependence of the modulation amplitude. Using a maximum likelihood test, we find that descriptions of the modulation amplitude as decreasing with time are preferred over a constant modulation amplitude at anywhere between 1$\sigma$ and 3$\sigma$, depending on the phenomenological model for the time dependence and the signal energy range considered. A time dependent modulation amplitude is not expected for a dark matter signal, at least for dark matter halo morphologies consistent with the DAMA signal. New data from DAMA/LIBRA--phase2 will certainly aid in determining whether any apparent time dependence is a real effect or a statistical fluctuation.

The density spike of dark matter (DM) in the subparsec region near the supermassive black hole at the Galactic Center can provide potentially observable gamma-ray signals coming from DM annihilations. Taking Fermi-$LAT$ data for the gamma-ray flux from the point source 3FGL J1745.6-2859c (Sgr A*), we calculate the resulting constraints on generic models of DM, allowing for the possibility of a non-negligible velocity-dependent component of the annihilation cross section. We consider a variety of selections for the astrophysical parameters that describe the spike profile and find that the gamma-ray flux is strongly dependent on these selections, particularly the modelling of spike depletion effects due to gravitational interactions with baryons, which affect the spike radius and steepness profile. We calculate constraints on the DM parameter space for both attenuated and idealized spikes, considering different choices for the steepness profiles in each case. We find that for the most conservative selection of parameters, corresponding to a depleted spike with an NFW cusp profile, the gamma-ray flux for a 100 GeV thermal relic is lower than current observational constraints by several orders of magnitude. For more optimistic choices of parameters corresponding to spikes that have not been attenuated, bounds on DM masses can be obtained for a variety of choices of steepness profiles. We then specialize to a class of simplified models of fermionic DM that annihilate dominantly through the $t-$channel exchange of two scalar mediators with arbitrary mixing angle $\alpha$, and calculate the indirect detection constraints coming from the DM spike. Along the way, we discuss constraints on the astrophysical parameters describing the DM spike, taking the Galactic Center excess as a signal of DM annihilation.

We consider a simplified model in which Majorana fermion dark matter annihilates to charged fermions through exchange of charged mediators. We consider the gamma-ray signals arising from the processes $XX \rightarrow \bar f f \gamma$, $\gamma \gamma$, and $\gamma Z$ in the most general case, including non-trivial fermion mass and non-trivial left-right mixing and $CP$-violating phase for the charged mediators. In particular, we find the most general spectrum for internal bremsstrahlung, which interpolates between the regimes dominated by virtual internal bremsstrahlung and by final state radiation. We also examine the variation in the ratio $\sigma(\gamma \gamma) / \sigma (\gamma Z)$ and the helicity asymmetry in the $XX \rightarrow \gamma \gamma$ process, each as a function of mixing angle and $CP$-violating phase. As an application, we apply these results to searches for a class of MSSM models.

We consider two potential non-accelerator signatures of generalizations of the well-studied constrained minimal supersymmetric standard model (CMSSM). In one generalization, the universality constraints on soft supersymmetry-breaking parameters are applied at some input scale $M_{in}$ below the grand unification (GUT) scale $M_{GUT}$, a scenario referred to as `sub-GUT'. The other generalization we consider is to retain GUT-scale universality for the squark and slepton masses, but to relax universality for the soft supersymmetry-breaking contributions to the masses of the Higgs doublets. As with other CMSSM-like models, the measured Higgs mass requires supersymmetric particle masses near or beyond the TeV scale. Because of these rather heavy sparticle masses, the embedding of these CMSSM-like models in a minimal SU(5) model of grand unification can yield a proton lifetime consistent with current experimental limits, and may be accessible in existing and future proton decay experiments. Another possible signature of these CMSSM-like models is direct detection of supersymmetric dark matter. The direct dark matter scattering rate is typically below the reach of the LUX-ZEPLIN (LZ) experiment if $M_{in}$ is close to $M_{GUT}$, but may lie within its reach if $M_{in} \lesssim 10^{11}$ GeV. Likewise, generalizing the CMSSM to allow non-universal supersymmetry-breaking contributions to the Higgs offers extensive possibilities for models within reach of the LZ experiment that have long proton lifetimes.

We consider searches for dark matter annihilation to monoenergetic neutrinos in the core of the Sun. We find that liquid scintillation neutrino detectors have enhanced sensitivity to this class of dark matter models, due to the energy and angular resolution possible for electron neutrinos and antineutrinos that scatter via charged-current interactions. In particular we find that KamLAND, utilizing existing data, could provide better sensitivity to such models than any current direct detection experiment for $m_X \lesssim 15$ GeV. KamLAND's sensitivity is signal-limited, and future liquid scintillation or liquid argon detectors with similar energy and angular resolution, but with larger exposure, will provide significantly better sensitivity. These detectors may be particularly powerful probes of dark matter with mass $\mathcal{O}(10)$ GeV.

Recent experimental results from the LHC have placed strong constraints on the masses of colored superpartners. The MSSM parameter space is also constrained by the measurement of the Higgs boson mass, and the requirement that the relic density of lightest neutralinos be consistent with observations. Although large regions of the MSSM parameter space can be excluded by these combined bounds, leptophilic versions of the MSSM can survive these constraints. In this paper we consider a scenario in which the requirements of minimal flavor violation, vanishing $CP$-violation, and mass universality are relaxed, specifically focusing on scenarios with light sleptons. We find a large region of parameter space, analogous to the original bulk region, for which the lightest neutralino is a thermal relic with an abundance consistent with that of dark matter. We find that these leptophilic models are constrained by measurements of the magnetic and electric dipole moments of the electron and muon, and that these models have interesting signatures at a variety of indirect detection experiments.

The DAMA experiment searches for Weakly Interacting Massive Particle (WIMP) dark matter via its expected but rare interactions within the detector, where the interaction rates will modulate throughout the year due to the orbital motion of the Earth. Over the course of more than 10 years of operation, DAMA has indeed detected a strong modulation in the event rate above the detector threshold of 2 keVee. Under standard assumptions regarding the dark matter halo and WIMP interactions, this signal is consistent with that expected of WIMPs of two different approximate masses: ~ 10 GeV and ~ 70 GeV. We examine how a lower threshold, allowed by recent upgrades to the DAMA detector, may shed light on this situation. We find that the lower threshold data should rule out one of the two mass ranges for spin-independent couplings (in the worst case, disfavoring one of the masses by still more than 2.6$\sigma$) and is likely, though not certain, to do the same for spin-dependent couplings. Furthermore, the data may indicate whether the interaction is predominantly spin-independent or spin-dependent in some cases. Our findings illustrate the importance of a low threshold in modulation searches.

In order for neutralino dark matter to avoid being overproduced in the early universe, these particles must annihilate (or coannihilate) rather efficiently. Neutralinos with sufficiently large couplings to annihilate at such high a rate (such as those resulting from gaugino-higgsino mixing, as in "well-tempered" or "focus point" scenarios), however, have become increasingly disfavored by the null results of XENON100 and other direct detection experiments. One of the few remaining ways that neutralinos could potentially evade such constraints is if they annihilate through a resonance, as can occur if 2$m_{\chi^0}$ falls within about $\sim$10% of either $m_{A/H}$, $m_h$, or $m_Z$. If no signal is observed from upcoming direct detection experiments, the degree to which such a resonance must be tuned will increase significantly. In this paper, we quantify the degree to which such a resonance must be tuned in order to evade current and projected constraints from direct detection experiments. Assuming a future rate of progress among direct detection experiments that is similar to that obtained over the past decade, we project that within 7 years the light Higgs and $Z$ pole regions will be entirely closed, while the remaining parameter space near the $A/H$ resonance will require that $2m_{\chi^0}$ be matched to the central value (near $m_A$) to within less than 4%. At this rate of progress, it will be a little over a decade before multi-ton direct detection experiments will be able to close the remaining, highly-tuned, regions of the $A/H$ resonance parameter space.

We calculate the limits on the fraction of viable dark matter minihalos in the early universe to host Population III.1 stars, surviving today as dark matter spikes in our Milky Way halo. Motivated by potential hints of light dark matter from the DAMA and CoGeNT direct dark matter searches, we consider thermal relic WIMP dark matter with masses of 5, 10, and 20 GeV, and annihilation to mu^+ mu^-, tau^+ tau^-, and q bar{q}. From this brief study we conclude that, if dark matter is light, either the typical black hole size is \lesssim 100 M_\odot (i.e. there is no significant Dark Star phase), and/or dark matter annihilates primarily to mu^+ mu^- or other final states that result in low gamma-ray luminosity, and/or that an extremely small fraction of minihalos in the early universe that seem suitable to host the formation of the first stars actually did.

We have studied Electroweak Symmetry Breaking (EWSB) fine-tuning in the context of two unified Supersymmetry scenarios: the Constrained Minimal Supersymmetric Model (CMSSM) and models with Non-Universal Higgs Masses (NUHM), in light of current and upcoming direct detection dark matter experiments. We consider both those models that satisfy a one-sided bound on the relic density of neutralinos, $\Omega_{\chi} h^2 < 0.12$, and also the subset that satisfy the two-sided bound in which the relic density is within the 2 sigma best fit of WMAP7 + BAO + H0 data. We find that current direct detection searches for dark matter probe the least fine-tuned regions of parameter-space, or equivalently those of lowest Higgs mass parameter $\mu$, and will tend to probe progressively more and more fine-tuned models, though the trend is more pronounced in the CMSSM than in the NUHM. Additionally, we examine several subsets of model points, categorized by common mass hierarchies; $M_{\chi_0} \sim M_{\chi^\pm}, M_{\chi_0} \sim M_{\stau}, M_{\chi_0} \sim M_{\stop_1}$, the light and heavy Higgs poles, and any additional models classified as "other"; the relevance of these mass hierarchies is their connection to the preferred neutralino annihilation channel that determines the relic abundance. For each of these subsets of models we investigated the degree of fine-tuning and discoverability in current and next generation direct detection experiments.

We consider both velocity-dependent and velocity-independent contributions to spin-dependent (SD) and spin-independent (SI) nuclear scattering (including one-loop corrections) of WIMPless dark matter, in the case where the dark matter candidate is a Majorana fermion. We find that spin-independent scattering arises only from the mixing of exotic squarks, or from velocity-dependent terms. Nevertheless (and contrary to the case of MSSM neutralino WIMPs), we find a class of models which cannot be detected through SI scattering, but can be detected at IceCube/DeepCore through SD scattering. We study the detection prospects for both SI and SD detection strategies for a large range of Majorana fermion WIMPless model parameters.

10 to 10^5 solar mass black holes with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today are examined in light of recent data from the Fermi Gamma-Ray Space Telescope (FGST). The dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation. We examine the signatures of annihilations into gamma-rays, electrons and positrons, and neutrinos. We find that some significant fraction of the point sources detected by FGST might be due to dark matter annihilation near black holes in our Galaxy. We obtain limits on the properties of dark matter annihilations in the spikes using the information in the FGST First Source Catalog as well as the diffuse gamma-ray flux measured by FGST. We determine the maximum fraction of high redshift minihalos that could have hosted the formation of the first generation of stars and, subsequently, their black hole remnants. The strength of the limits depends on the choice of annihilation channel and black hole mass; limits are strongest for the heaviest black holes and annhilation to $b \bar{b}$ and $W^+W^-$ final states. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted and thermonuclear burning runs its course; thus FGST observations may be used to constrain the properties of Dark Stars. Additionally, we comment on the excess positron flux found by PAMELA and its possible interpretation in terms of dark matter annihilation around these black hole spikes.

In this paper we examine whether indirect detection constraints on dark matter associated with a non-thermal history may be significantly improved when accounting for the presence of galactic substructure in the form of dark matter spikes. We find that significant constraints may be derived from the non-observation of an excess of diffuse gamma-rays and from the properties of bright gamma-ray point sources observed by the Fermi Gamma-Ray Space Telescope, but these constraints depend sensitively on the details of the formation of the first stars and their subsequent black hole remnants. However, we also find that, especially if WIMPs annihilate primarily to quarks or gauge bosons, it is possible to extract meaningful and conservative bounds on the annihilation cross section.

The very first stars likely formed from metal-free, molecular hydrogen-cooled gas at the centers of dark matter minihalos. Prior to nuclear fusion, these stars may have been supported by dark matter heating from annihilations in the star, in which case they could have grown to be quite massive before collapsing to black holes. Many remnant black holes and their surrounding dark matter density spikes may be part of our Milky Way halo today. Here we explore the gamma-ray signatures of dark matter annihilations in the dark matter spikes surrounding these black holes for a range of star formation scenarios, black hole masses, and dark matter annihilation modes. Data from the Fermi Gamma-Ray Space Telescope are used to constrain models of dark matter annihilation and the formation of the first stars.

In this paper, we compare the relative sensitivities of gamma-ray and neutrino observations to the dark matter annihilation cross section in leptophilic models such as have been designed to explain PAMELA data. We investigate whether the high energy neutrino telescope IceCube will be competitive with current and upcoming searches by gamma-ray telescopes, such as the Atmospheric Cerenkov Telescopes (ACTs) (HESS, VERITAS and MAGIC), or the Fermi Gamma Ray Space Telescope, in detecting or constraining dark matter particles annihilating in dwarf spheroidal galaxies. We find that after ten years of observation of the most promising nearby dwarfs, IceCube will have sensitivity comparable to the current sensitivity of gamma-ray telescopes only for very heavy (m_X > 7 TeV) or relatively light (m_X < 200 GeV) dark matter particles which annihilate primarily to mu+mu-. If dark matter particles annihilate primarily to tau+tau-, IceCube will have superior sensitivity only for dark matter particle masses below the 200 GeV threshold of current ACTs. If dark matter annihilations proceed directly to neutrino-antineutrino pairs a substantial fraction of the time, IceCube will be competitive with gamma-ray telescopes for a much wider range of dark matter masses.

A brief overview is given of the phenomenology of particle dark matter and the properties of some of the most widely studied dark matter candidates. Recent developments in direct and indirect dark matter searches are discussed.

We consider the potentials of the LHC and a linear e^+e^- collider (LC) for discovering supersymmetric particles in variants of the MSSM with soft supersymmetry-breaking mass parameters constrained to be universal at the GUT scale (CMSSM) or at some lower scale M_{in} (GUT-less models), as may occur in some scenarios with mirage unification. Whereas the LHC should be able to discover squarks and/or gluinos along all the CMSSM coannihilation strip where the relic neutralino LSP density lies within the range favoured for cold dark matter, many GUT-less models could escape LHC detection. In particular, if M_{in} < 10^{11} GeV, the LHC would not detect sparticles if the relic density lies within the favoured range. For any given discovery of supersymmetry at the LHC, in such GUT-less models the lightest neutralino mass and hence the threshold for sparticle pair production at a LC increases as M_{in} decreases, and the CMSSM offers the best prospects for measuring sparticles at a LC. For example, if the LHC discovers sparticles with 1 fb^{-1} of data, within the CMSSM a centre-of-mass energy of 600 GeV would suffice for a LC to to produce pairs of neutralinos, if they provide the cold dark matter, whereas over 1 TeV might be required in a general GUT-less model. These required energies increase to 800 GeV in the CMSSM and 1.4 TeV in GUT-less models if the LHC requires 10 fb^{-1} to discover supersymmetry.

We study models in which supersymmetry breaking appears at an intermediate scale, M_{in}, below the GUT scale. We assume that the soft supersymmetry-breaking parameters of the MSSM are universal at M_{in}, and analyze the morphology of the constraints from cosmology and collider experiments on the allowed regions of parameter space as M_{in} is reduced from the GUT scale. We present separate analyses of the (m_{1/2},m_0) planes for tan(beta)=10 and tan(beta)=50, as well as a discussion of non-zero trilinear couplings, A_0. Specific scenarios where the gaugino and scalar masses appear to be universal below the GUT scale have been found in mirage-mediation models, which we also address here. We demand that the lightest neutralino be the LSP, and that the relic neutralino density not conflict with measurements by WMAP and other observations. At moderate values of M_{in}, we find that the allowed regions of the (m_{1/2},m_0) plane are squeezed by the requirements of electroweak symmetry breaking and that the lightest neutralino be the LSP, whereas the constraint on the relic density is less severe. At very low M_{in}, the electroweak vacuum conditions become the dominant constraint, and a secondary source of astrophysical cold dark matter would be necessary to explain the measured relic density for nearly all values of the soft SUSY-breaking parameters and tan(beta). We calculate the neutralino-nucleon cross sections for viable scenarios and compare them with the present and projected limits from direct dark matter searches.

We consider the stochastic background of gravitational waves produced by an early generation of Population III stars coupled with a normal mode of star formation at lower redshift. The computation is performed in the framework of hierarchical structure formation and is based on cosmic star formation histories constrained to reproduce the observed star formation rate at redshift $z \la 6$, the observed chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an early reionization of the Universe at $z\sim 10-20$ as indicated by the first year results released by WMAP. We find that the normal mode of star formation produces a gravitational wave background which peaks at 300-500 Hz and is within LIGO III sensitivity. The Population III component peaks at lower frequencies (30-100 Hz depending on the model), and could be detected by LIGO III as well as the planned BBO and DECIGO interferometers.

We consider the relic neutrino background produced by Population III stars coupled with a normal mode of star formation at lower redshift. The computation is performed in the framework of hierarchical structure formation and is based on cosmic star formation histories constrained to reproduce the observed star formation rate at redshift z \la 6, the observed chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an early reionization of the Universe at z ~ 10-20. We consider both a burst and non-burst model for Population III star formation. We find that although the high redshift burst of Population III stars does lead to an appreciable flux of neutrinos at relatively low energy (E_\nu \approx 1 MeV), the observable neutrino flux is dominated by the normal mode of star formation. We also find that predicted fluxes are at the present level of the SuperK limit. As a consequence, the supernova relic neutrino background has a direct impact on models of chemical evolution and/or supernova dynamics.

Evidence from the WMAP polarization data indicates that the Universe may have been reionized at very high redshift. It is often suggested that the ionizing UV flux originates from an early population of massive or very massive stars. Depending on their mass, such stars can explode either as type II supernovae or pair-instability supernovae, or may entirely collapse into a black hole. The resulting neutrino emission can be quite different in each case. We consider here the relic neutrino background produced by an early burst of Population III stars coupled with a normal mode of star formation at lower redshift. The computation is performed in the framework of hierarchical structure formation and is based on cosmic star formation histories constrained to reproduce the observed star formation rate at redshift z \la 6, the observed chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an early reionization of the Universe at z \sim 10-20. We find that although the high redshift burst of Population III stars does lead to an appreciable flux of neutrinos at relatively low energy (E_\nu \approx 1 MeV), the observable neutrino flux is dominated by the normal mode of star formation. We also find that predicted fluxes are at the present level of the SuperK limit. As a consequence, the supernova relic neutrino background has a direct impact on models of chemical evolution and/or supernova dynamics.