We investigate a scenario where dark matter (DM) particles can be captured and accumulate in the Sun, and subsequently annihilate into a pair of long-lived mediators. These mediators can decay further out in the Sun or outside of the Sun. Compared to the standard scenario where DM particles annihilate directly into Standard Model particles close to the solar core, here we also obtain fluxes of gamma rays and charged cosmic rays. We simulate this scenario using a full three-dimensional model of the Sun, and include interactions and neutrino oscillations. In particular, we perform a model-independent study of the complementarity between neutrino and gamma ray fluxes by comparing the recent searches from IceCube, Super-Kamiokande, Fermi-LAT, ARGO and HAWC. We find that the resulting neutrino fluxes are significantly higher at high energy when the mediators decay further out in the Sun. We also find that gamma ray searches place stronger constraints than neutrino searches on these models even in cases where the mediators decay mainly inside the Sun, except in the approximately inner 10% of the Sun where neutrino searches are more powerful. We present our results in a model-independent manner and release a new version of the WimpSim code that can be used to simulate this scenario for arbitrary mediator models.

The nature of dark matter remains one of the key science questions. Weakly Interacting Massive Particles (WIMPs) are among the best motivated particle physics candidates, allowing to explain the measured dark matter density by employing standard big-bang thermodynamics. We introduce here a radically new version of the widely used DarkSUSY package, which allows to compute the properties of such dark matter particles numerically. With DarkSUSY 6 one can accurately predict a large variety of astrophysical signals from dark matter, such as direct detection rates in low-background counting experiments and indirect detection signals through antiprotons, antideuterons, gamma rays and positrons from the Galactic halo, or high-energy neutrinos from the center of the Earth or of the Sun. For thermally produced dark matter like WIMPs, high-precision tools are provided for the computation of the relic density in the Universe today, as well as for the size of the smallest dark matter protohalos. Furthermore, the code allows to calculate dark matter self-interaction rates, which may affect the distribution of dark matter at small cosmological scales. Compared to earlier versions, DarkSUSY 6 introduces many significant physics improvements and extensions. The most fundamental new feature of this release, however, is that the code has been completely re-organized and brought into a highly modular and flexible shape. Switching between different pre-implemented dark matter candidates has thus become straight-forward, just as adding new - WIMP or non-WIMP - particle models or replacing any given functionality in a fully user-specified way. In this article, we describe the physics behind the computer package, along with the main structure and philosophy of this major revision of DarkSUSY. A detailed manual is provided together with the public release at www.darksusy.org.

Cosmic rays hitting the solar atmosphere generate neutrinos that interact and oscillate in the Sun and oscillate on the way to Earth. These neutrinos could potentially be detected with neutrino telescopes and will be a background for searches for neutrinos from dark matter annihilation in the Sun. We calculate the flux of neutrinos from these cosmic ray interactions in the Sun and also investigate the interactions near a detector on Earth that give rise to muons. We compare this background with both regular Earth-atmospheric neutrinos and signals from dark matter annihilation in the Sun. Our calculation is performed with an event-based Monte Carlo approach that should be suitable as a simulation tool for experimental collaborations. Our program package is released publicly along with this paper.

Dark matter in the form of Weakly Interacting Massive Particles (WIMPs) can be captured by the Sun and the Earth, sink to their cores, annihilate and produce neutrinos that can be searched for with neutrino telescopes. The calculation of the capture rates of WIMPs in the Sun and especially the Earth are affected by large uncertainties coming mainly from effects of the planets in the Solar System, reducing the capture rates by up to an order of magnitude (or even more in some cases). We show that the WIMPs captured by weak scatterings in the Sun also constitute an important bound WIMP population in the Solar System. Taking this population and its interplay with the population bound through gravitational diffusion into account cancel the planetary effects on the capture rates, and the capture essentially proceeds as if the Sun and the Earth were free in the galactic halo. The neutrino signals from the Sun and the Earth are thus significantly higher than claimed in the scenarios with reduced capture rates.

We show, by using an extensive sample of viable supersymmetric models as templates, that indirect detection of dark matter through gamma rays may have a large potential for identifying the nature of dark matter. This is in particular true also for models that give too weak dark matter-nucleon scattering cross sections to be probed by present and planned direct detection experiments. Also models with a mass scale too high to be accessible at CERN's LHC accelerator may show up in next-generation imaging Cherenkov telescope arrays. Based on our our findings, we therefore suggest to view indirect searches as genuine particle physics experiments, complementing other strategies to probe so far unknown regions in the parameter space of e.g. supersymmetric models, and propose a new approach that would make use of telescopes dedicated for dark matter searches. As a concrete example for the potential of such an approach, we consider an array of imaging air Cherenkov telescopes, the Dark Matter Array (DMA), and show that such an experiment could extend present-day limits by several orders of magnitude, reaching a large class of models that would remain undetected in both direct detection experiments and searches at the LHC. In addition, in a sizable part of the parameter space, signals from more than one type of dark matter detection experiment would be possible, something that may eventually be necessary in order to identify the dark matter candidate.

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

Galactic weakly interacting massive particles (WIMPs) may scatter off solar nuclei to orbits gravitationally bound to the Sun. Once bound, the WIMPs continue to lose energy by repeated scatters in the Sun, eventually leading to complete entrapment in the solar interior. While the density of the bound population is highest at the center of the Sun, the only observable signature of WIMP annihilations inside the Sun is neutrinos. It has been previously suggested that although the density of WIMPs just outside the Sun is lower than deep inside, gamma rays from WIMP annihilation just outside the surface of the Sun, in the so called WIMP halo around the Sun, may be more easily detected. We here revisit this problem using detailed Monte Carlo simulations and detailed composition and structure information about the Sun to estimate the size of the gamma-ray flux. Compared to earlier simpler estimates, we find that the gamma-ray flux from WIMP annihilations in the solar WIMP halo would be negligible; no current or planned detectors would be able to detect this flux.

With present and upcoming SUSY searches both directly, indirectly and at accelerators, the need for accurate calculations is large. We will here go through some of the tools available both from a dark matter point of view and at accelerators. For natural reasons, we will focus on public tools, even though there are some rather sophisticated private tools as well.

The possibility of explaining the positron and electron excess recently found by the PAMELA and ATIC collaborations in terms of dark matter (DM) annihilation has attracted considerable attention. Models surviving bounds from, e.g, antiproton production generally fall into two classes, where either DM annihilates directly with a large branching fraction into light leptons, or, as in the recent models of Arkani-Hamed et al., and of Nomura and Thaler, the annihilation gives low-mass (pseudo)scalars or vectors $\phi$ which then decay into $\mu^+\mu^-$ or $e^+e^-$. While the constraints on the first kind of models have recently been treated by several authors, we study here specifically models of the second type which rely on an efficient Sommerfeld enhancement in order to obtain the necessary boost in the annihilation cross section. We compute the photon flux generated by QED radiative corrections to the decay of $\phi$ and show that this indeed gives a rather spectacular broad peak in $E^2d\sigma/dE$, that for these extreme values of the cross section violate gamma-ray observations of the Galactic center for DM density profiles steeper than that of Navarro, Frenk and White. The most stringent constraint comes from the comparison of the predicted synchrotron radiation in the central part of the Galaxy with radio observations of Sgr A*. For the most commonly adopted DM profiles, the models that provide a good fit to the PAMELA and ATIC data are ruled out, unless there are physical processes that boost the local anti-matter fluxes more than one order of magnitude, while not affecting the gamma-ray or radio fluxes.

We analyze the recently released Fermi-LAT data on the sum of electrons and positrons. Compared to a conventional, pre-Fermi, background model, a surprising excess in the several hundred GeV range is found and here we analyze it in terms of dark matter models. We also compare with newly published results from PAMELA and HESS, and find models giving very good fits to these data sets as well. If this dark matter interpretation is correct, we also predict the possibility of a sharp break in the diffuse gamma ray spectrum coming from final state radiation.

The inert doublet model is a minimal extension of the standard model introducing an additional SU(2) doublet with new scalar particles that could be produced at accelerators. While there exists no LEP II analysis dedicated for these inert scalars, the absence of a signal within searches for supersymmetric neutralinos can be used to constrain the inert doublet model. This translation however requires some care because of the different properties of the inert scalars and the neutralinos. We investigate what restrictions an existing DELPHI collaboration study of neutralino pair production can put on the inert scalars and discuss the result in connection with dark matter. We find that although an important part of the inert doublet model parameter space can be excluded by the LEP II data, the lightest inert particle still constitutes a valid dark matter candidate.

The space-borne antimatter experiment PAMELA has recently reported a surprising rise in the positron to electron ratio at high energies. It has also recently been found that electromagnetic radiative corrections in some cases may boost the gamma-ray yield from supersymmetric dark matter annihilations in the galactic halo by up to three or four orders of magnitude, providing distinct spectral signatures for indirect dark matter searches to look for. Here, we investigate whether the same type of corrections can also lead to sizeable enhancements in the positron yield. We find that this is indeed the case, albeit for a smaller region of parameter space than for gamma rays; selecting models with a small mass difference between the neutralino and sleptons, like in the stau coannihilation region in mSUGRA, the effect becomes more pronounced. The resulting, rather hard positron spectrum with a relatively sharp cutoff may potentially fit the rising positron ratio measured by the PAMELA satellite. To do so, however, very large "boost factors" have to be invoked that are not expected in current models of halo structure. If the predicted cutoff would also be confirmed by later PAMELA data or upcoming experiments, one could either assume non-thermal production in the early universe or non-standard halo formation to explain such a spectral feature as an effect of dark matter annihilation. At the end of the paper, we briefly comment on the impact of radiative corrections on other annihilation channels, in particular antiprotons and neutrinos.

We compute the electromagnetic radiative corrections to all leading annihilation processes which may occur in the Galactic dark matter halo, for dark matter in the framework of supersymmetric extensions of the Standard Model (MSSM and mSUGRA), and present the results of scans over the parameter space that is consistent with present observational bounds on the dark matter density of the Universe. Although these processes have previously been considered in some special cases by various authors, our new general analysis shows novel interesting results with large corrections that may be of importance, e.g., for searches at the soon to be launched GLAST gamma-ray space telescope. In particular, it is pointed out that regions of parameter space where there is a near degeneracy between the dark matter neutralino and the tau sleptons, radiative corrections may boost the gamma-ray yield by up to three or four orders of magnitude, even for neutralino masses considerably below the TeV scale, and will enhance the very characteristic signature of dark matter annihilations, namely a sharp step at the mass of the dark matter particle. Since this is a particularly interesting region for more constrained mSUGRA models of supersymmetry, we use an extensive scan over this parameter space to verify the significance of our findings. We also re-visit the direct annihilation of neutralinos into photons and point out that, for a considerable part of the parameter space, internal bremsstrahlung is more important for indirect dark matter searches than line signals.

We investigate the sensitivity of the Gamma-ray Large Area Space Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) through the $\gamma$-ray signal that their pair annihilation produces. WIMPs are among the favorite candidates to explain the compelling evidence that about 80% of the mass in the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by various extensions of the standard model of particle physics such as Supersymmetry and Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very large energy range (20 MeV to more than 300 GeV) the main instrument on board the GLAST satellite, the Large Area Telescope (LAT), will open a new window of discovery. As our estimates show, the LAT will be able to detect an indirect DM signature for a large class of WIMP models given a cuspy profile for the DM distribution. Using the current state of the art Monte Carlo and event reconstruction software developed within the LAT collaboration, we present preliminary sensitivity studies for several possible sources inside and outside the Galaxy. We also discuss the potential of the LAT to detect UED via the electron/positron channel. Diffuse background modeling and other background issues that will be important in setting limits or seeing a signal are presented.

Weakly Interacting Massive Particles (WIMPs) are one of the main candidates for the dark matter in the Universe. If these particles make up the dark matter, then they can be captured by the Sun or the Earth, sink to the respective cores, annihilate, and produce neutrinos. Thus, these neutrinos can be a striking dark matter signature at neutrino telescopes looking towards the Sun and/or the Earth. Here, we improve previous analyses on computing the neutrino yields from WIMP annihilations in several respects. We include neutrino oscillations in a full three-flavor framework as well as all effects from neutrino interactions on the way through the Sun (absorption, energy loss, and regeneration from tau decays). In addition, we study the effects of non-zero values of the mixing angle $\theta_{13}$ as well as the normal and inverted neutrino mass hierarchies. Our study is performed in an event-based setting which makes these results very useful both for theoretical analyses and for building a neutrino telescope Monte Carlo code. All our results for the neutrino yields, as well as our Monte Carlo code, are publicly available. We find that the yield of muon-type neutrinos from WIMP annihilations in the Sun is enhanced or suppressed, depending on the dominant WIMP annihilation channel. This effect is due to an effective flavor mixing caused by neutrino oscillations. For WIMP annihilations inside the Earth, the distance from source to detector is too small to allow for any significant amount of oscillations at the neutrino energies relevant for neutrino telescopes.

We modify a stellar structure code to estimate the effect upon the main sequence of the accretion of weakly interacting dark matter onto stars and its subsequent annihilation. The effect upon the stars depends upon whether the energy generation rate from dark matter annihilation is large enough to shut off the nuclear burning in the star. Main sequence WIMP burners look much like protostars moving on the Hayashi track, although they are in principle completely stable. We make some brief comments about where such stars could be found, how they might be observed and more detailed simulations which are currently in progress. Finally we comment on whether or not it is possible to link the paradoxically young OB stars found at the galactic centre with WIMP burners.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

One way to unambiguously confirm the existence of particle dark matter and determine its mass would be to detect its annihilation into monochromatic gamma-rays in upcoming telescopes. One of the most minimal models for dark matter is the inert doublet model, obtained by adding another Higgs doublet with no direct coupling to fermions. For a mass between 40 and 80 GeV, the lightest of the new inert Higgs particles can give the correct cosmic abundance of cold dark matter in agreement with current observations. We show that for this scalar dark matter candidate, the annihilation signal of monochromatic \gamma\gamma and Z\gamma final states would be exceptionally strong. The energy range and rates for these gamma-ray line signals make them ideal to search for with the soon upcoming GLAST satellite.

Measurement of the extragalactic background (EGBR) of diffuse gamma-rays is perhaps one of the most challenging tasks for future gamma-ray observatories, such as GLAST. This is because any determination will depend on accurate subtraction of the galactic diffuse and celestial foregrounds, as well as point sources. However, the EGBR is likely to contain very rich information about the high energy-gamma ray sources of the Universe at cosmological distances. We focus on the ability of GLAST to detect a signal from dark matter in the EGBR. We present sensitivities for generic thermal WIMPs and the Inert Higgs Doublet Model. Also we discuss the various aspects of astrophysics and particle physics that determines the shape and strength of the signal, such as dark matter halo properties and different dark matter candidates. Other possible sources to the EGBR are also discussed, such as unresolved AGNs, and viewed as backgrounds.

We investigate the internal consistency of the halo dark matter model which has been proposed by de Boer et al. to explain the excess of diffuse galactic gamma rays observed by the EGRET experiment. Any model based on dark matter annihilation into quark jets, such as the supersymmetric model proposed by de Boer et al., inevitably also predicts a primary flux of antiprotons from the same jets. Since propagation of the antiprotons in the unconventional, disk-dominated type of halo model used by de Boer et al. is strongly constrained by the measured ratio of boron to carbon nuclei in cosmic rays, we investigate the viability of the model using the DarkSUSY package to compute the gamma-ray and antiproton fluxes. We are able to show that their model is excluded by a wide margin from the measured flux of antiprotons. We therefore find that a model of the type suggested by Moskalenko et al., where the intensities of protons and electrons in the cosmic rays vary with galactic position, is far more plausible to explain the gamma excess.

The question of the nature of the dark matter in the Universe remains one of the most outstanding unsolved problems in basic science. One of the best motivated particle physics candidates is the lightest supersymmetric particle, assumed to be the lightest neutralino - a linear combination of the supersymmetric partners of the photon, the Z boson and neutral scalar Higgs particles. Here we describe DarkSUSY, a publicly-available advanced numerical package for neutralino dark matter calculations. In DarkSUSY one can compute the neutralino density in the Universe today using precision methods which include resonances, pair production thresholds and coannihilations. Masses and mixings of supersymmetric particles can be computed within DarkSUSY or with the help of external programs such as FeynHiggs, ISASUGRA and SUSPECT. Accelerator bounds can be checked to identify viable dark matter candidates. DarkSUSY also computes a large variety of astrophysical signals from neutralino dark matter, such as direct detection in low-background counting experiments and indirect detection through antiprotons, antideuterons, gamma-rays and positrons from the Galactic halo or high-energy neutrinos from the center of the Earth or of the Sun. Here we describe the physics behind the package. A detailed manual will be provided with the computer package.

We perform a detailed analysis of the detection prospects of neutralino dark matter in the mSUGRA framework. We focus on models with a thermal relic density, estimated with high accuracy using the DarkSUSY package, in the range favored by current precision cosmological measurements. Direct and indirect detection rates are computed implementing two models for the dark matter halo, tracing opposite regimes for the phase of baryon infall, with fully consistent density profiles and velocity distribution functions. This has allowed, for the first time, a fully consistent comparison between direct and indirect detection prospects. We discuss all relevant regimes in the mSUGRA parameter space, underlining relevant effects, and providing the basis for extending the discussion to alternative frameworks. In general, we find that direct detection and searches for antideuterons in the cosmic rays seems to be the most promising ways to search for neutralinos in these scenarios.

Weakly Interacting Massive Particles (WIMPs) can be captured by the Earth, where they eventually sink to the core, annihilate and produce e.g. neutrinos that can be searched for with neutrino telescopes. The Earth is believed to capture WIMPs not dominantly from the Milky Way halo directly, but instead from a distribution of WIMPs that have diffused around in the solar system due to gravitational interactions with the planets in the solar system. Recently, doubts have been raised about the lifetime of these WIMP orbits due to solar capture. We here investigate this issue by detailed numerical simulations. Compared to earlier estimates, we find that the WIMP velocity distribution is significantly suppressed below about 70 km/s which results in a suppression of the capture rates mainly for heavier WIMPs (above ~100 GeV). At 1 TeV and above the reduction is almost a factor of 10. We apply these results to the case where the WIMP is a supersymmetric neutralino and find that, within the Minimal Supersymmetric Standard Model (MSSM), the annihilation rates, and thus the neutrino fluxes, are reduced even more than the capture rates. At high masses (above ~1 TeV), the suppression is almost two orders of magnitude. This suppression will make the detection of neutrinos from heavy WIMP annihilations in the Earth much harder compared to earlier estimates.

Neutralinos arise as natural dark matter candidates in many supersymmetric extensions of the Standard Model. We present a novel calculation of the neutralino relic abundance in which we include all so called coannihilation processes between neutralinos, charginos and sfermions, and, at the same time, we apply the state of the art technique to trace the freeze-out of a species in the early Universe. As a first application, we discuss here results valid in the mSUGRA framework; we enlight general trends as well as perform a detailed study of the neutralino relic densities in the mSUGRA parameter space. The emerging picture is fair agreement with previous analyses in the same framework, however we have the power to discuss it in many more details than previously done. E.g., we find that the cosmological bound on the neutralino mass is pushed up to ~565 GeV in the stau coannihilation region and to ~1500 GeV in the chargino coannihilation region.

There is mounting evidence for dark matter in the Universe and one of the favourite dark matter candidates is the neutralino, which naturally appears as the lightest supersymmetric particle (LSP) in many supersymmetric extensions of the standard model. The neutralino has the desired properties to be a good dark matter candidate and we will here review the different indirect searches for neutralino dark matter and discuss the implications on these from recent direct searches.

Using a new instrument, the HEAT collaboration has confirmed the excess of cosmic ray positrons that they first detected in 1994. We explore the possibility that this excess is due to the annihilation of neutralino dark matter in the galactic halo. We confirm that neutralino annihilation can produce enough positrons to make up the measured excess only if there is an additional enhancement to the signal. We quantify the `boost factor' that is required in the signal for various models in the Minimal Supersymmetric Standard Model parameter space, and find that a boost factor >30 provides good fits to the HEAT data. Such an enhancement in the signal could arise if we live in a clumpy halo.

The question of the nature of dark matter in the Universe remains one of the most outstanding unsolved problems in basic science. One of the best motivated particle physics candidates is the lightest supersymmetric particle, assumed to be the lightest neutralino. We here describe DarkSUSY, an advanced numerical FORTRAN package for supersymmetric dark matter calculations. With DarkSUSY one can: (i) compute masses and compositions of various supersymmetric particles; (ii) compute the relic density of the lightest neutralino, using accurate methods which include the effects of resonances, pair production thresholds and coannihilations; (iii) check accelerator bounds to identify allowed supersymmetric models; and (iv) obtain neutralino detection rates for a variety of detection methods, including direct detection and indirect detection through antiprotons, gamma-rays and positrons from the Galactic halo or neutrinos from the center of the Earth or the Sun.

We investigate the prospects of detecting weakly interacting massive particle (WIMP) dark matter by measuring the contribution to the extragalactic gamma-ray radiation induced, in any dark matter halo and at all redshifts, by WIMP pair annihilations into high-energy photons. We perform a detailed analysis of the distinctive spectral features of this signal, recently proposed in a short letter by three of the authors, with emphasis on the signature due to monochromatic gamma-ray yields: the combined effect of cosmological redshift and absorption along the line of sight produces sharp bumps, peaked at the rest frame energy of the lines and asymmetrically smeared to lower energies. The level of the flux depends both on the particle physics scenario for WIMP dark matter and on the question of how dark matter clusters. Uncertainties introduced by the latter are thoroughly discussed implementing a realistic model inspired by results of the state-of-the-art N-body simulations and semi-analytic modeling in the cold dark matter structure formation theory. We also address the question of the potential gamma-ray background originating from blazars, presenting a novel calculation. Comparing the signal with the background, we find that there are viable configurations, in the combined parameter space defined by the particle physics setup and the structure formation scenario, for which the WIMP induced extragalactic gamma-ray signal will be detectable in the new generation of gamma-ray telescopes such as GLAST.

Using a new instrument, the HEAT collaboration has confirmed the excess of cosmic ray positrons that they first detected in 1994. We explore the possibility that this excess is due to the annihilation of neutralino dark matter in the galactic halo. We confirm that neutralino annihilation can produce enough positrons to make up the measured excess only if there is an additional enhancement to the signal. We quantify the `boost factor' that is required in the signal for various models in the Minimal Supersymmetric Standard Model parameter space, and study the dependence on various parameters. We find models with a boost factor greater than 30. Such an enhancement in the signal could arise if we live in a clumpy halo. We discuss what part of supersymmetric parameter space is favored (in that it gives the largest positron signal), and the consequences for other direct and indirect searches of supersymmetric dark matter.

We propose a new signature for weakly interacting massive particle (WIMP) dark matter, a spectral feature in the diffuse extragalactic gamma-ray radiation. This feature, a sudden drop of the gamma-ray intensity at an energy corresponding to the WIMP mass, comes from the asymmetric distortion of the line due to WIMP annihilation into two gamma-rays caused by the cosmological redshift. Unlike other proposed searches for a line signal, this method is not very sensitive to the exact dark matter density distribution in halos and subhalos. The only requirement is that the mass distribution of substructure on small scales follows approximately the Press-Schechter law, and that smaller halos are on the average denser than large halos, which is a generic outcome of N-body simulations of Cold Dark Matter, and which has observational support. The upcoming Gamma-ray Large Area Space Telescope (GLAST) will be eminently suited to search for these spectral features. For numerical examples, we use rates computed for supersymmetric particle dark matter, where a detectable signal is possible.

There is mounting evidence that a self-consistent model for particle cold dark matter has to take into consideration spatial inhomogeneities on sub-galactic scales seen, for instance, in high-resolution N-body simulations of structure formation. Also in more idealized, analytic models, there appear density enhancements in certain regions of the halo. We use the results from a recent N-body simulation of the Milky Way halo and investigate the gamma-ray flux which would be produced when a specific dark matter candidate, the neutralino, annihilates in regions of enhanced density. The clumpiness found on all scales in the simulation results in very strong gamma-ray signals which seem to already rule out some regions of the supersymmetric parameter space, and would be further probed by upcoming experiments, such as the GLAST gamma-ray satellite. As an orthogonal model of structure formation, we also consider Sikivie's simple infall model of dark matter which predicts that there should exist continuous regions of enhanced density, caustic rings, in the dark matter halo of the Milky Way. We find, however, that the gamma-ray signal from caustic rings is generally too small to be detectable.

The question of the nature of the dark matter in the Universe remains one of the most outstanding unsolved problems in basic science. One of the best motivated particle physics candidates is the lightest supersymmetric particle, assumed to be the lightest neutralino. We here describe DarkSUSY, an advanced numerical FORTRAN package for supersymmetric dark matter calculations which we release for public use. With the help of this package, the masses and compositions of various supersymmetric particles can be computed, for given input parameters of the minimal supersymmetric extension of the Standard Model (MSSM). For the lightest neutralino, the relic density is computed, using accurate methods which include the effects of resonances, pair production thresholds and coannihilations. Accelerator bounds are checked to identify viable dark matter candidates. Finally, detection rates are computed for a variety of detection methods, such as direct detection and indirect detection through antiprotons, gamma-rays and positrons from the Galactic halo or neutrinos from the center of the Earth or the Sun.

A new Solar System population of Weakly Interacting Massive Particle (WIMP) dark matter has been proposed to exist. We investigate the implications of this population on indirect signals in neutrino telescopes (due to WIMP annihilations in the Earth) for the case when the WIMP is the lightest neutralino of the MSSM, the minimal supersymmetric extension of the standard model. The velocity distribution and capture rate of this new population is evaluated and the flux of neutrino-induced muons from the center of the Earth in neutrino telescopes is calculated. We show that the effects of the new population can be crucial for masses around 60-120 GeV, where enhancements of the predicted muon flux from the center of the Earth by up to a factor of 100 compared to previously published estimates occur. As a result of the new WIMP population, neutrino telescopes should be able to probe a much larger region of parameter space in this mass range.

Recently, a new Solar System population of weakly interacting massive particle (WIMP) dark matter has been proposed to exist. We investigate the implications of this population on indirect signals in neutrino telescopes (due to WIMP annihilations in the Earth) for the case when the WIMP is the lightest neutralino of the MSSM, the minimal supersymmetric extension of the standard model. The velocity distribution and capture rate of this new population is evaluated and the flux of neutrino-induced muons from the center of the Earth in neutrino telescopes is calculated. The strength of the signal is very sensitive to the velocity distribution of the new population. We analytically estimate this distribution using the approximate conservation of the component of the WIMP angular momentum orthogonal to the ecliptic plane. The non-linear problem of combining a fixed capture rate from the standard galactic WIMP population with one rising linearly with time from the new population to obtain the present-day annihilation rate in the Earth is also solved analytically. We show that the effects of the new population can be crucial for masses below around 150 GeV, where enhancements of the predicted muon flux from the center of the Earth by up to a factor of 100 compared to previously published estimates occur. As a result of the new WIMP population, the next generation of neutrino telescopes should be able to probe a much larger region of parameter space in the mass range 60-130 GeV.

The flux of cosmic ray antiprotons from neutralino annihilations in the galactic halo is computed for a large sample of models in the Minimal Supersymmetric extension of the Standard Model. We also revisit the problem of estimating the background of low-energy cosmic ray induced secondary antiprotons, taking into account their subsequent interactions (and energy loss) and the presence of nuclei in the interstellar matter. We point out that in some cases the optimal kinetic energy to search for a signal from supersymmetric dark matter is above several GeV, rather than the traditional sub-GeV region. The large astrophysical uncertainties involved do not allow the exclusion of any of the MSSM models we consider, on the basis of current data.

Neutrino telescopes of kilometer size are currently being planned. They will be two or three orders of magnitude larger than presently operating detectors, but they will have a much higher muon energy threshold. We discuss the trade-off between area and energy threshold for indirect detection of neutralino dark matter captured in the Sun and in the Earth and annihilating into high energy neutrinos. We also study the effect of a higher threshold on the complementarity of different searches for supersymmetric dark matter.

The flux of cosmic ray antiprotons from neutralino annihilations in the galactic halo is computed for a large sample of models in the MSSM (the Minimal Supersymmetric extension of the Standard Model). We also revisit the problem of estimating the background of low-energy cosmic ray induced secondary antiprotons, taking into account their subsequent interactions (and energy loss) and the presence of nuclei in the interstellar matter. We consider a two-zone diffusion model, with and without a galactic wind. We find that, given the uncertainties in the background predictions, there is no need for a primary (exotic) component to explain present data. However, allowing for a signal by playing with the uncertainties in the background estimate, we discuss the characteristic features of the supersymmetric models which give a satisfactory description of the data. We point out that in some cases the optimal kinetic energy to search for a signal from supersymmetric dark matter is above several GeV, rather than the traditional sub-GeV region. The large astrophysical uncertainties involved do not, one the other hand, allow the exclusion of any of the MSSM models we consider, on the basis of data. We present besides numerical results also convenient parameterizations of the antiproton yields of all `basic' two-body final states. We also give examples of the yield and differential energy spectrum for a set of supersymmetric models with high rates. We also remark that it is difficult to put a limit on the antiproton lifetime from present measurements, since the injection of antiprotons from neutralino annihilation can compensate the loss from decay.

Supersymmetric neutralinos are one of the most promising candidates for the dark matter in the Universe. If they exist, they should make up some fraction of the Milky Way halo. We investigate the fluxes of positrons expected at the Earth from neutralino annihilation in the halo. Positron propagation is treated in a diffusion model including energy loss. The positron source function includes contributions from both continuum and monochromatic positrons. We find that, for a "canonical" halo model and propagation parameters, the fluxes are generally too low to be visible. Given the large uncertainties in both propagation and halo structure, it is however possible to obtain observable fluxes. We also investigate the shapes of the positron spectra, including fits to a feature indicated by the results of the HEAT experiment.

Neutrino telescopes of kilometer size are currently being planned. They will be two or three orders of magnitude bigger than presently operating detectors, but they will have a much higher muon energy threshold. We discuss the trade-off between area and energy threshold for indirect detection of neutralino dark matter captured in the Sun and in the Earth and annihilating into high energy neutrinos. We also study the effect of a higher threshold on the complementarity of different searches for supersymmetric dark matter.

We investigate the possibility to detect neutralino dark matter in a scenario in which the galactic dark halo is clumpy. We find that under customary assumptions on various astrophysical parameters, the antiproton and continuum gamma-ray signals from neutralino annihilation in the halo put the strongest limits on the clumpiness of a neutralino halo. We argue that indirect detection through neutrinos from the Earth and the Sun should not be much affected by clumpiness. We identify situations in parameter space where the gamma-ray line, positron and diffuse neutrino signals from annihilations in the halo may provide interesting signals in upcoming detectors.

We investigate if the gamma ray halo, for which recent evidence has been found in EGRET data, can be explained by neutralino annihilations in a clumpy halo. We find that the measured excess gamma ray flux can be explained through a moderate amount of clumping in the halo. Moreover, the required amount of clumping implies also a measureable excess of antiprotons at low energies, for which there is support from recent measurements by the BESS collaboration. The predicted antiproton fluxes resulting from neutralino annihilations in a clumpy halo are high enough to give an excess over cosmic-ray produced antiprotons also at moderately high energies (above a few GeV). This prediction, as well as that of one or two sharp gamma lines coming from annihilations into 2 gammas or Z gamma can be tested in upcoming space-borne experiments like AMS and GLAST.

We evaluate the relic density of the lightest neutralino, the lightest supersymmetric particle, in the Minimal Supersymmetric extension of the Standard Model (MSSM). For the first time, we include all coannihilation processes between neutralinos and charginos for any neutralino mass and composition. We use the most sophisticated routines for integrating the cross sections and the Boltzmann equation. We properly treat (sub)threshold and resonant annihilations. We also include one-loop corrections to neutralino masses. We find that coannihilation processes are important not only for light higgsino-like neutralinos, as pointed out before, but also for heavy higgsinos and for mixed and gaugino-like neutralinos. Indeed, coannihilations should be included whenever $|\mu| \lsim 2 |M_1|$, independently of the neutralino composition. When $|\mu| \sim |M_1|$, coannihilations can increase or decrease the relic density in and out of the cosmologically interesting region. We find that there is still a window of light higgsino-like neutralinos that are viable dark matter candidates and that coannihilations shift the cosmological upper bound on the neutralino mass from 3 to 7 TeV.

We investigate the improvement in sensitivity to astrophysical point sources of energetic ($\gtrsim1$ GeV) neutrinos which can be achieved with angular and/or energy resolution of the neutrino-induced muon. As a specific example we consider WIMP annihilation in the Sun and in the Earth as a neutrino source. The sensitivity is improved by using the angular and energy distribution to reduce the atmospheric-neutrino background. Although the specific improvements depend on the backgrounds and assumed sources, the sensitivity to a WIMP signal may be improved, with equal exposure, by up to roughly a factor of two with good angular resolution, and by up to roughly a factor of three with good energy resolution. In case of a positive detection, energy resolution would also improve the measurement of the neutrino energy spectrum and therefore provide information on the WIMP mass and composition.