We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

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 spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.

We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overline{\psi} \psi$ and pseudoscalar $\overline{\psi} i\gamma_5 \psi$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory.

Dark matter searches in gamma ray final states often make use of the fact that photons can be produced from final state muons. Modern Monte Carlo generators and DM codes include the effects of final state radiation from muons produced in the dark matter annihilation process itself, but neglect the O(1%) radiative correction that arises from the subsequent muon decay. After implementing this correction we demonstrate the effect that it can have on dark matter phenomenology by considering the case of dark matter annihilation to four muons via scalar mediator production. We first show that the AMS-02 positron excess can no longer easily be made consistent with this final state once the Fermi-LAT dwarf limits are calculated with the inclusion of radiative muon decays, and we next show that the Fermi-LAT galactic centre gamma excess can be improved with this final state after inclusion of the same effect. We provide code and tables for the implementation of this effect in the popular dark matter code micrOMEGAs, providing a solution for any model producing final state muons.

We present an improved event-level likelihood formalism for including neutrino telescope data in global fits to new physics. We derive limits on spin-dependent dark matter-proton scattering by employing the new formalism in a re-analysis of data from the 79-string IceCube search for dark matter annihilation in the Sun, including explicit energy information for each event. The new analysis excludes a number of models in the weak-scale minimal supersymmetric standard model (MSSM) for the first time. This work is accompanied by the public release of the 79-string IceCube data, as well as an associated computer code for applying the new likelihood to arbitrary dark matter models.

Weakly Interacting Massive Particles (WIMPs) may constitute a large fraction of the matter in the Universe. There are excess events in the data of DAMA/LIBRA, CoGeNT, CRESST-II, and recently CDMS-Si, which could be consistent with WIMP masses of approximately 10 GeV/c2. However, for MDM > 10 GeV/c2 null results of the CDMS-Ge, XENON, and LUX detectors may be in tension with the potential detections for certain dark matter scenarios and assuming a certain light response. We propose the use of a new class of biological dark matter (DM) detectors to further examine this light dark matter hypothesis, taking advantage of new signatures with low atomic number targets, Two types of biological DM detectors are discussed here: DNA-based detectors and enzymatic reactions (ER) based detectors. In the case of DNA-based detectors, we discuss a new implementation. In the case of ER detectors, there are four crucial phases of the detection process: a) change of state due to energy deposited by a particle; b) amplification due to the release of energy derived from the action of an enzyme on its substrate; c) sustainable but non-explosive enzymatic reaction; d) self-termination due to the denaturation of the enzyme, when the temperature is raised. This paper provides information of how to design as well as optimize these four processes.

We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling limits from the recent results of the LUX direct search dark matter experiment. After a brief review of WIMP-nucleon scattering, we derive LUX limits on the spin-dependent WIMP-nucleon couplings over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical parameters. We find that, under these and other common assumptions, LUX excludes the entire spin-dependent parameter space consistent with a dark matter interpretation of DAMA's anomalous signal, the first time a single experiment has been able to do so. We also revisit the case of spin-independent couplings, and demonstrate good agreement between our results and the published LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to place constraints on generic dark matter theories.

It is shown that a Weakly Interacting Massive dark matter Particle (WIMP) interpretation for the positron excess observed in a variety of experiments, HEAT, PAMELA, and AMS-02, is highly constrained by the Fermi/LAT observations of dwarf galaxies. In particular, this paper has focused on the annihilation channels that best fit the current AMS-02 data (Boudaud et al., 2014). The Fermi satellite has surveyed the $\gamma$-ray sky, and its observations of dwarf satellites are used to place strong bounds on the annihilation of WIMPs into a variety of channels. For the single channel case, we find that dark matter annihilation into {$b\bar{b}$, $e^+e^-$, $\mu^+\mu^-$, $\tau^+\tau^-$, 4-$e$, or 4-$\tau$} is ruled out as an explanation of the AMS positron excess (here $b$ quarks are a proxy for all quarks, gauge and Higgs bosons). In addition, we find that the Fermi/LAT 2$\sigma$ upper limits, assuming the best-fit AMS-02 branching ratios, exclude multichannel combinations into $b\bar{b}$ and leptons. The tension between the results might relax if the branching ratios are allowed to deviate from their best-fit values, though a substantial change would be required. Of all the channels we considered, the only viable channel that survives the Fermi/LAT constraint and produces a good fit to the AMS-02 data is annihilation (via a mediator) to 4-$\mu$, or mainly to 4-$\mu$ in the case of multichannel combinations.

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.

Direct detection experiments, which are designed to detect the scattering of dark matter off nuclei in detectors, are a critical component in the search for the Universe's missing matter. The count rate in these experiments should experience an annual modulation due to the relative motion of the Earth around the Sun. This modulation, not present for most known background sources, is critical for solidifying the origin of a potential signal as dark matter. In this article, we review the physics of annual modulation, discussing the practical formulae needed to interpret a modulating signal. We focus on how the modulation spectrum changes depending on the particle and astrophysics models for the dark matter. For standard assumptions, the count rate has a cosine dependence with time, with a maximum in June and a minimum in December. Well-motivated generalizations of these models, however, can affect both the phase and amplitude of the modulation. We show how a measurement of an annually modulating signal could teach us about the presence of substructure in the Galactic halo or about the interactions between dark and baryonic matter. Although primarily a theoretical review, we briefly discuss the current experimental situation for annual modulation and future experimental directions.

We analyse the sensitivity of IceCube-DeepCore to annihilation of neutralino dark matter in the solar core, generated within a 25 parameter version of the minimally supersymmetric standard model (MSSM-25). We explore the 25-dimensional parameter space using scanning methods based on importance sampling and using DarkSUSY 5.0.6 to calculate observables. Our scans produced a database of 6.02 million parameter space points with neutralino dark matter consistent with the relic density implied by WMAP 7-year data, as well as with accelerator searches. We performed a model exclusion analysis upon these points using the expected capabilities of the IceCube-DeepCore Neutrino Telescope. We show that IceCube-DeepCore will be sensitive to a number of models that are not accessible to direct detection experiments such as SIMPLE, COUPP and XENON100, indirect detection using Fermi-LAT observations of dwarf spheroidal galaxies, nor to current LHC searches.

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

We investigate the interactions of Weakly Interacting Massive Particles (WIMPs) with nuclei in the human body. We are motivated by the fact that WIMPs are excellent candidates for the dark matter in the Universe. Our estimates use a 70 kg human and a variety of WIMP masses and cross-sections. The contributions from individual elements in the body are presented and it is found that the dominant contribution is from scattering off of oxygen (hydrogen) nuclei for the spin-independent (spin-dependent) interactions. For the case of 60 GeV WIMPs, we find that, of the billions of WIMPs passing through a human body per second, roughly ~10 WIMPs hit one of the nuclei in the human body in an average year, if the scattering is at the maximum consistent with current bounds on WIMP interactions. We also study the 10-20 GeV WIMPs with much larger cross-sections that best fit the DAMA, COGENT, and CRESST data sets and find much higher rates: in this case as many as $10^5$ WIMPs hit a nucleus in the human body in an average year, corresponding to almost one a minute. Though WIMP interactions are a source of radiation in the body, the annual exposure is negligible compared to that from other natural sources (including radon and cosmic rays), and the WIMP collisions are harmless to humans.

We examine the future sensitivity of CoGeNT to the presence of dark matter streams and find that consideration of streams in the data may lead to differences in the interpretation of the results. We show the allowed particle mass and cross section for different halo parameters, assuming spin-independent elastic scattering. As an example, we choose a stream with the same velocity profile as that of the Sagittarius stream (and in the Solar neighborhood) and find that, with an exposure of $\sim$ 10 kg year, the CoGeNT results can be expected to exclude the SHM-only halo in favor of an SHM+stream halo at the 95% (99.7%) confidence level provided the stream contributes 3% (5%) of the local dark matter density. The presence of a significant stream component may result in incorrect estimates of the particle mass and cross section unless the presence of the stream is taken into account. We conclude that the CoGeNT experiment is sensitive to streams and care should be taken to include the possibility of streams when analyzing experimental results.

Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.

We extend our previous studies of the neutrino fluxes expected from neutralino LSP annihilations inside the Sun to include variants of the minimal supersymmetric extension of the Standard Model (MSSM) with squark, slepton and gaugino masses constrained to be universal at the GUT scale, but allowing one or two non-universal supersymmetry-breaking parameters contributing to the Higgs masses (NUHM1,2). As in the constrained MSSM (CMSSM) with universal Higgs masses, there are large regions of the NUHM parameter space where the LSP density inside the Sun is not in equilibrium, so that the annihilation rate may be far below the capture rate, and there are also large regions where the capture rate is not dominated by spin-dependent LSP-proton scattering. The spectra possible in the NUHM are qualitatively similar to those in the CMSSM. We calculate neutrino-induced muon fluxes above a threshold energy of 10 GeV, appropriate for the IceCube/DeepCore detector, for points where the NUHM yields the correct cosmological relic density for representative choices of the NUHM parameters. We find that the IceCube/DeepCore detector can probe regions of the NUHM parameter space in addition to analogues of the focus-point strip and the tip of the coannihilation strip familiar from the CMSSM. These include regions with enhanced Higgsino-gaugino mixing in the LSP composition, that occurs where neutralino mass eigenstates cross over. On the other hand, rapid-annihilation funnel regions in general yield neutrino fluxes that are unobservably small.

Weakly interacting massive particles (WIMPs) are amongst the most interesting dark matter (DM) candidates. Many DM candidates naturally arise in theories beyond the standard model (SM) of particle physics, like weak-scale supersymmetry (SUSY). Experiments aim to detect WIMPs by scattering, annihilation or direct production, and thereby determine the underlying theory to which they belong, along with its parameters. Here we examine the prospects for further constraining the Constrained Minimal Supersymmetric Standard Model (CMSSM) with future ton-scale direct detection experiments. We consider ton-scale extrapolations of three current experiments: CDMS, XENON and COUPP, with 1000 kg-years of raw exposure each. We assume energy resolutions, energy ranges and efficiencies similar to the current versions of the experiments, and include backgrounds at target levels. Our analysis is based on full likelihood constructions for the experiments. We also take into account present uncertainties on hadronic matrix elements for neutralino-quark couplings, and on halo model parameters. We generate synthetic data based on four benchmark points and scan over the CMSSM parameter space using nested sampling. We construct both Bayesian posterior PDFs and frequentist profile likelihoods for the model parameters, as well as the mass and various cross-sections of the lightest neutralino. Future ton-scale experiments will help substantially in constraining supersymmetry, especially when results of experiments primarily targeting spin-dependent nuclear scattering are combined with those directed more toward spin-independent interactions.

Data taken during the final shallow-site run of the first tower of the Cryogenic Dark Matter Search (CDMS II) detectors have been reanalyzed with improved sensitivity to small energy depositions. Four ~224 g germanium and two ~105 g silicon detectors were operated at the Stanford Underground Facility (SUF) between December 2001 and June 2002, yielding 118 live days of raw exposure. Three of the germanium and both silicon detectors were analyzed with a new low-threshold technique, making it possible to lower the germanium and silicon analysis thresholds down to the actual trigger thresholds of ~1 keV and ~2 keV, respectively. Limits on the spin-independent cross section for weakly interacting massive particles (WIMPs) to elastically scatter from nuclei based on these data exclude interesting parameter space for WIMPs with masses below 9 GeV/c^2. Under standard halo assumptions, these data partially exclude parameter space favored by interpretations of the DAMA/LIBRA and CoGeNT experiments' data as WIMP signals, and exclude new parameter space for WIMP masses between 3 GeV/c^2 and 4 GeV/c^2.

The determination of dark matter constraints from liquid xenon direct detection experiments depends upon the amount of scintillation light produced by nuclear recoils in the detector, a quantity that is characterized by the scintillation efficiency factor Leff. We examine how uncertainties in the measurements of Leff and the extrapolated behavior of Leff at low recoil energies (where measurements do not exist) affect the constraints from experiments such as XENON10 and XENON100, particularly in the light WIMP regions of interest for the DAMA and CoGeNT experimental results.

We consider the compatibility of DAMA/LIBRA, CoGeNT, XENON10 and XENON100 results for spin-independent (SI) dark matter Weakly Interacting Massive Particles (WIMPs), particularly at low masses (~ 10 GeV), assuming a standard dark matter halo. The XENON bounds depend on the scintillation efficiency factor Leff for which there is considerable uncertainty. Thus we consider various extrapolations for Leff at low energy. With the Leff measurements we consider, XENON100 results are found to be insensitive to the low energy extrapolation. We find the strongest bounds are from XENON10, rather than XENON100, due to the lower energy threshold. For reasonable choices of Leff and for the case of SI elastic scattering, XENON10 is incompatible with the DAMA/LIBRA 3$\sigma$ region and severely constrains the 7-12 GeV WIMP mass region of interest published by the CoGeNT collaboration.

We evaluate the neutrino fluxes to be expected from neutralino LSP annihilations inside the Sun, within the minimal supersymmetric extension of the Standard Model with supersymmetry-breaking scalar and gaugino masses constrained to be universal at the GUT scale (the CMSSM). We find that there are large regions of typical CMSSM $(m_{1/2}, m_0)$ planes where the LSP density inside the Sun is not in equilibrium, so that the annihilation rate may be far below the capture rate. We show that neutrino fluxes are dependent on the solar model at the 20% level, and adopt the AGSS09 model of Serenelli et al. for our detailed studies. We find that there are large regions of the CMSSM $(m_{1/2}, m_0)$ planes where the capture rate is not dominated by spin-dependent LSP-proton scattering, e.g., at large $m_{1/2}$ along the CMSSM coannihilation strip. We calculate neutrino fluxes above various threshold energies for points along the coannihilation/rapid-annihilation and focus-point strips where the CMSSM yields the correct cosmological relic density for tan(beta) = 10 and 55 for $\mu$ > 0, exploring their sensitivities to uncertainties in the spin-dependent and -independent scattering matrix elements. We also present detailed neutrino spectra for four benchmark models that illustrate generic possibilities within the CMSSM. Scanning the cosmologically-favored parts of the parameter space of the CMSSM, we find that the IceCube/DeepCore detector can probe at best only parts of this parameter space, notably the focus-point region and possibly also at the low-mass tip of the coannihilation strip.

The DAMA/NaI and DAMA/LIBRA annual modulation data, which may be interpreted as a signal for the existence of weakly interacting dark matter (WIMPs) in our galactic halo, are re-examined in light of new measurements of the local velocity relative to the galactic halo. In the vicinity of the Sun, the velocity of the Galactic disk has been estimated to be 250 km/s rather than 220 km/s. Our analysis is performed both with and without the channeling effect included. The best fit regions to the DAMA data are shown to move to slightly lower WIMP masses. Compatibility of DAMA data with null results from other experiments (CDMS, XENON10, and CRESST I) is investigated given these new velocities. A small region of spin-independent (elastic) scattering for 7-8 GeV WIMP masses remains at 3$\sigma$. Spin-dependent scattering off of protons is viable for 5-15 GeV WIMP masses for direct detection experiments (but has been argued by others to be further constrained by Super-Kamiokande due to annihilation in the Sun).

The DAMA/NaI and DAMA/LIBRA annual modulation data, which may be interpreted as a signal for the existence of weakly interacting dark matter (WIMPs) in our galactic halo, are examined in light of null results from other experiments. We use the energy spectrum of the combined DAMA modulation data given in 36 bins, and include the effect of channeling. Several statistical tools are implemented in our study: likelihood ratio with a global fit and with raster scans in the WIMP mass and goodness-of-fit (g.o.f.). These approaches allow us to differentiate between the preferred (global best fit) and allowed (g.o.f.) parameter regions. It is hard to find WIMP masses and couplings consistent with all existing data sets. For spin-independent (SI) interactions, the best fit DAMA regions are ruled out to the 3$\sigma$ C.L., even with channeling taken into account. However, for WIMP masses of ~8 GeV some parameters outside these regions still yield a moderately reasonable fit to the DAMA data and are compatible with all 90% C.L. upper limits from negative searches, when channeling is included. For spin-dependent (SD) interactions with proton-only couplings, a range of masses below 10 GeV is compatible with DAMA and other experiments, with and without channeling, when SuperK indirect detection constraints are included; without the SuperK constraints, masses as high as ~20 GeV are compatible. For SD neutron-only couplings we find no parameters compatible with all the experiments. Mixed SD couplings are examined: e.g. ~8 GeV mass WIMPs with a_n = +/- a_p are found to be consistent with all experiments. In short, there are surviving regions at low mass for both SI and SD interactions; if indirect detection limits are relaxed, some SD proton-only couplings at high masses also survive.

We review the uncertainties in the spin-independent and -dependent elastic scattering cross sections of supersymmetric dark matter particles on protons and neutrons. We propagate the uncertainties in quark masses and hadronic matrix elements that are related to the $\pi$-nucleon $\sigma$ term and the spin content of the nucleon. By far the largest single uncertainty is that in spin-independent scattering induced by our ignorance of the $<N | {\bar q} q | N>$ matrix elements linked to the $\pi$-nucleon $\sigma$ term, which affects the ratio of cross sections on proton and neutron targets as well as their absolute values. This uncertainty is already impacting the interpretations of experimental searches for cold dark matter. {\it We plead for an experimental campaign to determine better the $\pi$-nucleon $\sigma$ term.} Uncertainties in the spin content of the proton affect significantly, but less strongly, the calculation of rates used in indirect searches.

The model of Natural Inflation is examined in light of recent 3-year data from the Wilkinson Microwave Anisotropy Probe and shown to provide a good fit. The inflaton potential is naturally flat due to shift symmetries, and in the simplest version takes the form $V(\phi) = \Lambda^4 [1 \pm \cos(N\phi/f)]$. The model agrees with WMAP3 measurements as long as $f > 0.7 m_{Pl}$ (where $m_{Pl} = 1.22 \times 10^{19}$GeV) and $\Lambda \sim m_{GUT}$. The running of the scalar spectral index is shown to be small -- an order of magnitude below the sensitivity of WMAP3. The location of the field in the potential when perturbations on observable scales are produced is examined; for $f > 5 m_{Pl}$, the relevant part of the potential is indistinguishable from a quadratic, yet has the advantage that the required flatness is well-motivated. Depending on the value of $f$, the model falls into the large field ($f \ge 1.5 m_{Pl}$) or small field ($f < 1.5 m_{Pl}$) classification scheme that has been applied to inflation models. Natural inflation provides a good fit to WMAP3 data.

In addition to a smooth component of WIMP dark matter in galaxies, there may be streams of material; the effects of WIMP streams on direct detection experiments is examined in this paper. The contribution to the count rate due to the stream cuts off at some characteristic energy. Near this cutoff energy, the stream contribution to the annual modulation of recoils in the detector is comparable to that of the thermalized halo, even if the stream represents only a small portion (~5% or less) of the local halo density. Consequently the total modulation may be quite different than would be expected for the standard halo model alone: it may not be cosine-like and can peak at a different date than expected. The effects of speed, direction, density, and velocity dispersion of a stream on the modulation are examined. We describe how the observation of a modulation can be used to determine these stream parameters. Alternatively, the presence of a dropoff in the recoil spectrum can be used to determine the WIMP mass if the stream speed is known. The annual modulation of the cutoff energy together with the annual modulation of the overall signal provide a "smoking gun" for WIMP detection.

We obtain models for a triaxial Milky Way spheroid based on data by Newberg and Yanny. The best fits to the data occur for a spheroid center that is shifted by 3kpc from the Galactic Center. We investigate effects of the triaxiality on the microlensing optical depth to the Large Magellanic Cloud (LMC). The optical depth can be used to ascertain the number of Massive Compact Halo Objects (MACHOs); a larger spheroid contribution would imply fewer Halo MACHOs. On the one hand, the triaxiality gives rise to more spheroid mass along the line of sight between us and the LMC and thus a larger optical depth. However, shifting the spheroid center leads to an effect that goes in the other direction: the best fit to the spheroid center is_away_ from the line of sight to the LMC. As a consequence, these two effects tend to cancel so that the change in optical depth due to the Newberg/Yanny triaxial halo is at most 50%. After subtracting the spheroid contribution in the four models we consider, the MACHO contribution (central value) to the mass of the Galactic Halo varies from \~(8-20)% if all excess lensing events observed by the MACHO collaboration are assumed to be due to MACHOs. Here the maximum is due to the original MACHO collaboration results and the minimum is consistent with 0% at the 1 sigma error level in the data.

The age of the universe is obtained in a subset of Cardassian models by using WMAP data. Cardassian expansion is a modification to the Friedmann equation that allows the universe to be flat, matter dominated, and accelerating, without a vacuum component. Since this model changes the evolution of the universe, we should not a priori expect the Cardassian age to be the same as the WMAP Friedmann derived result of 13.7 +/- 0.2 Gyrs. However, in the subset of Cardassian models we consider, we discover that the age of the universe varies from 13.4 - 13.8 Gyr over the range of parameter space we explore, a result close to that of the standard Lambda Cold Dark Matter model. The Hubble constant h, which may also vary in these models, likewise varies little from the Friedmann result.

We examine whether the annual modulation found by the DAMA dark matter experiment can be explained by Weakly Interacting Massive Particles (WIMPs), in light of new null results from other experiments. CDMS II has already ruled out most WIMP-nucleus spin-independent couplings as an explanation for DAMA data. Hence we here focus on spin-dependent (axial vector; SD) couplings of WIMPs to nuclei. We expand upon previous work by (i) considering the general case of coupling to both protons and neutrons and (ii) incorporating bounds from all existing experiments. We note the surprising fact that CMDS II places one of the strongest bounds on the WIMP-neutron cross-section, and show that SD WIMP-neutron scattering alone is excluded. We also show that SD WIMP-proton scattering alone is allowed only for WIMP masses in the 5-13 GeV range. For the general case of coupling to both protons and neutrons, we find that, for WIMP masses above 13 GeV and below 5 GeV, there is no region of parameter space that is compatible with DAMA and all other experiments. In the range (5-13) GeV, we find acceptable regions of parameter space, including ones in which the WIMP-neutron coupling is comparable to the WIMP-proton coupling.