We present version X of the hammurabi package, the HEALPix-based numeric simulator for Galactic polarized emission. Improving on its earlier design, we have fully renewed the framework with modern C++ standards and features. Multi-threading support has been built in to meet the growing computational workload in future research. For the first time, we present precision profiles of hammurabi line-of-sight integral kernel with multi-layer HEALPix shells. In addition to fundamental improvements, this report focuses on simulating polarized synchrotron emission with Gaussian random magnetic fields. Two fast methods are proposed for realizing divergence-free random magnetic fields either on the Galactic scale where a field alignment and strength modulation are imposed, or on a local scale where more physically motivated models like a parameterized magneto-hydrodynamic (MHD) turbulence can be applied. As an example application, we discuss the phenomenological implications of Gaussian random magnetic fields for high Galactic latitude synchrotron foregrounds. In this, we numerically find B/E polarization mode ratios lower than unity based on Gaussian realizations of either MHD turbulent spectra or in spatially aligned magnetic fields.

Milky Way dwarf spheroidal satellites are a prime target for Dark Matter (DM) indirect searches. There have been recent reassessments of the expected DM gamma-ray signals in case of long-range interactions, commonly known as Sommerfeld enhancement. Since details of the underlying DM phase-space distribution function become critical, there are potentially large uncertainties in the final result. We provide here a first attempt towards a comprehensive investigation of these systematics, addressing the impact on the expected DM flux from Milky Way dwarfs via Bayesian inference on the available stellar kinematic datasets. After reconsidering the study case of ergodic systems, we investigate for the first time scenarios where DM particle orbits may have a radial or tangential bias. We consider both cuspy and cored parametric DM density profiles, together with the case of a non-parametric halo modelling directly connected to line-of-sight observable quantities. The main findings of our work highlight the relevance of the assumed phase-space distribution: Referring to a generalized J-factor, namely the line-of-sight convolution of the spatial part in case of velocity-dependent annihilation rate, an enhancement (suppression) with respect to the limit of isotropic phase-space distributions is obtained for the case of tangentially (radially) biased DM particle orbits. We provide new estimates for J-factors for the eight brightest Milky Way dwarfs also in the limit of velocity-independent DM annihilation, in good agreement with previous results in literature, and derive data-driven lower-bounds based on the non-parametric modelling of the halo density. The outcome of our broad study stands out as a representative of the state-of-the-art in the field, and falls within the interest of current and future experimental collaborations involved in DM indirect detection programs.

We address the problem of reconstructing the phase-space distribution function for an extended collisionless system, with known density profile and in equilibrium within an axisymmetric gravitational potential. Assuming that it depends on only two integrals of motion, namely the energy and the component of the angular momentum along the axis of symmetry $L_z$, there is a one-to-one correspondence between the density profile and the component of the distribution function that is even in $L_z$, as well as between the weighted azimuthal velocity profile and the odd component. This inversion procedure was originally proposed by Lynden-Bell and later refined in its numerical implementation by Hunter & Qian; after overcoming a technical difficulty, we apply it here for the first time in presence of a strongly flattened component, as a novel approach of extracting the phase-space distribution function for dark matter particles in the halo of spiral galaxies. We compare results obtained for realistic axisymmetric models to those in the spherical symmetric limit as assumed in previous analyses, showing the rather severe shortcomings in the latter. We then apply the scheme to the Milky Way and discuss the implications for the direct dark matter searches. In particular, we reinterpret the null results of the Xenon1T experiment for spin-(in)dependent interactions and make predictions for the annual modulation of the signal for a set of axisymmetric models, including a self-consistently defined co-rotating halo.

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.

We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions. The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios. The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well. The new version facilitates the users to include their own physical models by means of a modular C++ structure.

Dwarf satellite galaxies are ideal laboratories for identifying particle Dark Matter signals. When setting limits on particle Dark Matter properties from null searches, it becomes however crucial the level at which the Dark Matter density profile within these systems is constrained by observations. In the limit in which the spherical Jeans equation is assumed to be valid for a given tracer stellar population, we study the solution of this equation having the Dark Matter mass profile as an output rather than as a trial parametric input. Within our new formulation, we address to what level dwarf spheroidal galaxies feature a reliable mass estimator. We assess then possible extrapolation of the density profiles in the inner regions and -- keeping explicit the dependence on the orbital anisotropy profile of the tracer population -- we derive general trends on the line-of-sight integral of the density profile squared, a quantity commonly dubbed $J$-factor and crucial to estimate fluxes from prompt Dark Matter pair annihilations. Taking Ursa Minor as a study case among Milky Way satellites, we perform Bayesian inference using the available kinematical data for this galaxy. Contrary to all previous studies, we avoid marginalization over quantities poorly constrained by observations or by theoretical arguments. We find minimal $J$-factors to be about 2 to 4 times smaller than commonly quoted estimates, approximately relaxing by the same amount the limit on Dark Matter pair annihilation cross section from gamma-ray surveys of Ursa Minor. At the same time, if one goes back to a fixed trial parametric form for the density, e.g. using a NFW or Burkert profile, we show that the minimal $J$ can hardly be reduced by more than a factor of 1.5.

A spherical-symmetric gamma-ray emission from the central region of the Galaxy has been recently identified in Fermi-LAT data, and initially associated to dark matter particle annihilations. Guided by the evidence for a high gas density in the inner kpc of the Galaxy correlated with a very large Supernova rate, and hence with ongoing cosmic-ray acceleration, we investigate instead the possibility of addressing this excess in terms of ordinary cosmic-ray sources and standard steady-state diffusion. We alter the source term, and consistently the correlated gamma-ray emissions, in the context of a template-fitting analysis. We focus on a region of interest (ROI) defined as: $|l|<20^{\circ}; \hspace{5mm} 2^{\circ}<|b|<20^{\circ},$ with $l$ and $b$ the Galactic longitude and latitude coordinates. We analyze in detail the overall goodness of the fit of our framework, and perform a detailed direct comparison against data examining profiles in different directions. Remarkably, the test statistic of the fit related to our scenario turns out to be as good as the Dark Matter one in the ROI here considered.

The standard approach to cosmic-ray (CR) propagation in the Galaxy is based on the assumption that local transport properties can be extrapolated to the whole CR confining volume. Such models tend to underestimate the gamma-ray flux above few GeV measured by the Fermi Large Area Telescope towards the inner Galactic plane. We consider here for the first time a phenomenological scenario allowing for both the rigidity scaling of the diffusion coefficient and convective effects to be position-dependent. We show that within this approach we can reproduce the observed gamma-ray spectra at both low and mid Galactic latitudes - including the Galactic center - without spoiling any local CR observable.

We perform a multichannel analysis of the indirect signals for the Wino Dark Matter, including one-loop electroweak and Sommerfeld enhancement corrections. We derive limits from cosmic ray antiprotons and positrons, from continuum galactic and extragalactic diffuse gamma-ray spectra, from the absence of gamma-ray line features at the galactic center above 500 GeV in energy, from gamma-rays toward nearby dwarf spheroidal galaxies and galaxy clusters, and from CMB power-spectra. Additionally, we show the future prospects for neutrino observations toward the inner Galaxy and from antideuteron searches. For each of these indirect detection probes we include and discuss the relevance of the most important astrophysical uncertainties that can impact the strength of the derived limits. We find that the Wino as a dark matter candidate is excluded in the mass range bellow ~800 GeV from antiprotons and between 1.8 and 3.5 TeV from the absence of a gamma-ray line feature toward the galactic center. Limits from other indirect detection probes confirm the main bulk of the excluded mass ranges.

Recent advances in gamma-ray cosmic ray, infrared and radio astronomy have allowed us to develop a significantly better understanding of the galactic medium properties in the last few years. In this work using the DRAGON code, that numerically solves the CR propagation equation and calculating gamma-ray emissivities in a 2-dimensional grid enclosing the Galaxy, we study in a self consistent manner models for the galactic diffuse gamma-ray emission. Our models are cross-checked to both the available CR and gamma-ray data. We address the extend to which dark matter annihilations in the Galaxy can contribute to the diffuse gamma-ray flux towards different directions on the sky. Moreover we discuss the impact that astrophysical uncertainties of non DM nature, have on the derived gamma-ray limits. Such uncertainties are related to the diffusion properties on the Galaxy, the interstellar gas and the interstellar radiation field energy densities. Light ~10 GeV dark matter annihilating dominantly to hadrons is more strongly constrained by gamma-ray observations towards the inner parts of the Galaxy and influenced the most by assumptions of the gas distribution; while TeV scale DM annihilating dominantly to leptons has its tightest constraints from observations towards the galactic center avoiding the galactic disk plane, with the main astrophysical uncertainty being the radiation field energy density. In addition, we present a method of deriving constraints on the dark matter distribution profile from the diffuse gamma-ray spectra. These results critically depend on the assumed mass of the dark matter particles and the type of its end annihilation products.

Recently a series of indications have been put forward suggesting the presence of two gamma-ray lines at 110-130 GeV (centered at 111 and 129 GeV). Signals of these lines have been observed toward the Galactic center, at some galaxy clusters and among some of the unassociated point sources of the 2 years Fermi catalogue. Such a combination of signals could be generated by dark matter annihilations in the main dark matter halo, its substructures and nearby galaxy clusters. We discuss in this work the consistency between the number of events observed at the line energies in the sky and the predictions using results from the Via Lactea II numerical simulation and extrapolations below its mass resolution, taking into account that the annihilation cross-section to the lines can be estimated from the Galactic center signal. We find that some extrapolations to small substructures can naturally account for the point sources signal, although the hypothesis of background only cannot be rejected. We also study the morphology of the gamma-ray sky at the 2 lines energies, testing different Galactic diffuse background models to account for interstellar medium uncertainties, and different assumptions on the DM diffuse component profile. We find from template fits that within reasonable diffuse background uncertainties the presence of a spherical halo component is preferred with cuspier dark matter halo profiles being preferable even from the full sky fit. We finally check the impact of a dark disk component, suggested by cosmological simulations that include baryons and find that thin dark disks can not be disfavored, thus possibly accounting for the preferentially closer to the Galactic disk distribution of the point sources lines signal.

The complementarity between dark matter searches at colliders and in underground laboratories is an extraordinarily powerful tool in the quest for dark matter. In the vast majority of the analyses conducted so far these dark matter detection strategies have been profitably combined either to perform global fits in the context of certain particle physics models (e.g. the CMSSM) or to estimate the prospects for a direct dark matter detection given the LHC potential of discovering new physics beyond the Standard Model. In this paper we propose an alternative strategy to combine direct and collider dark matter searches: employing the potential of the upcoming generation of 1-ton direct detection experiments, we show that for certain supersymmetric configurations it is possible to translate the information encoded in an hypothetically discovered direct detection signal into classes of expected signals at the LHC. As an illustrative application of our method, we show that for a 60 GeV neutralino thermally produced via resonant annihilations and identified by a 1-ton direct detection experiment, our approach allows to forecast a clearly identifiable prediction for a LHC final state involving three leptons and missing energy. The strategy presented in this paper to systematically translate a direct detection signal into a prediction for the LHC has the potential to significantly strengthen the complementarity between these two dark matter detection strategies.

Indications for a gamma-ray line(s) signal towards the Galactic center at an energy of about 130 GeV have been recently presented. While dark matter annihilations are a viable candidate for this signal, it is generally expected that such a flux would be correlated to a gamma-ray component with continuum energy spectrum due to dark matter pair annihilating into other Standard Model particles. We use the gamma-ray data from the inner 10 x 10 degree window to derive limits for a variety of DM annihilation final states. Extending the window of observation, we discuss bounds on the morphological shape of a dark matter signal associated to the line, applying both standard templates for the dark matter profile, such as an Einasto or a NFW profile, and introducing a new more general parametrization.

We present a new determination of the local dark matter phase-space density. This result is obtained implementing, in the limit of isotropic velocity distribution and spherical symmetry, Eddington's inversion formula, which links univocally the dark matter distribution function to the density profile, and applying, within a Bayesian framework, a Markov Chain Monte Carlo algorithm to sample mass models for the Milky Way against a broad and variegated sample of dynamical constraints. We consider three possible choices for the dark matter density profile, namely the Einasto, NFW and Burkert profiles, finding that the velocity dispersion, which characterizes the width in the distribution, tends to be larger for the Burkert case, while the escape velocity depends very weakly on the profile, with the mean value we obtain being in very good agreement with estimates from stellar kinematics. The derived dark matter phase-space densities differ significantly--most dramatically in the high velocity tails--from the model usually taken as a reference in dark matter detection studies, a Maxwell-Boltzmann distribution with velocity dispersion fixed in terms of the local circular velocity and with a sharp truncation at a given value of the escape velocity. We discuss the impact of astrophysical uncertainties on dark matter scattering rates and direct detection exclusion limits, considering a few sample cases and showing that the most sensitive ones are those for light dark matter particles and for particles scattering inelastically. As a general trend, when adopting a self-consistent phase-space density, we find that rates are larger, and hence exclusion limits stronger, than with the standard Maxwell-Boltzmann approximation. Tools for applying our result on the local dark matter phase-space density to other dark matter candidates or experimental setups are provided.

The latest years have seen steady progresses in WIMP dark matter (DM) searches, with hints of possible signals suggested by both direct and indirect detection experiments. Antiprotons can play a key role validating those interpretations since they are copiously produced by WIMP annihilations in the Galactic halo, and the secondary antiproton background produced by Cosmic Ray (CR) interactions is predicted with fair accuracy and matches the observed spectrum very well. Using the publicly available numerical DRAGON code, we reconsider antiprotons as a tool to constrain DM models discussing its power and limitations. We provide updated constraints on a wide class of annihilating DM models by comparing our predictions against the most up-to-date ap measurements, taking also into account the latest spectral information on the p, He and other CR nuclei fluxes. Doing that, we probe carefully the uncertainties associated to both secondary and DM originated antiprotons, by using a variety of distinctively different assumptions for the propagation of CRs and for the DM distribution in the Galaxy. We find that the impact of the astrophysical uncertainties on constraining the DM properties can be much stronger, up to a factor of ~50, than the one due to uncertainties on the DM distribution (~2-6). Remarkably, even reducing the uncertainties on the propagation parameters derived by local observables, non-local effects can still change DM model constraints even by 50%. Nevertheless, current ap data place tight constraints on DM models, excluding some of those suggested in connection with indirect and direct searches. Finally we discuss the power of upcoming CR spectral data from the AMS-02 observatory to drastically reduce the uncertainties discussed in this paper and estimate the expected sensitivity of this instrument to some sets of DM models.

We study the high latitude (|b|>10) diffuse gamma-ray emission in the Galaxy in light of the recently published data from the Fermi collaboration at energies between 100 MeV and 100 GeV. The unprecedented accuracy in these measurements allows to probe and constrain the properties of sources and propagation of cosmic rays (CRs) in the Galaxy, as well as confirming conventional assumptions made on the interstellar medium (ISM). Using the publicly available DRAGON code, that has been shown to reproduce local measurements of CRs, we study assumptions made in the literature on HI and H2 gas distributions in the ISM, and non spatially uniform models of diffusion in the Galaxy. By performing a combined analysis of CR and gamma-ray spectra, we derive constraints on the properties of the ISM gas distribution and the vertical scale height of galactic CR diffusion, which may have implications also on indirect Dark Matter detection. We also discuss some of the possible interpretations of the break at ~230 GV in CR protons and helium spectra, recently observed by PAMELA and their impact on gamma-rays.

The spectral data on the diffuse Galactic gamma-rays, at medium and high latitudes (|b| > 10) and energies of 1-100 GeV, recently published by the Fermi Collaboration are used to produce a novel study on the gamma-ray emissivity in the Galaxy. We focus on analyzing the properties of propagation of cosmic rays (CRs), using the publicly available DRAGON code. We critically address some of the models for the interstellar HI and H2 gas distributions commonly used in the literature, as well as test a variety of propagation models. Each model assumes a distinct global profile for the diffusion and the re-acceleration of CRs. Fitting propagation parameters to well measured local CRs such as, the B/C ratio, protons, antiprotons and electron, positron fluxes, we evaluate the gamma-ray spectra at medium and high latitudes in order to place further constraints on these propagation models.

Non-thermal dark matter generation is an appealing alternative to the standard paradigm of thermal WIMP dark matter. We reconsider non-thermal production mechanisms in a systematic way, and develop a numerical code for accurate computations of the dark matter relic density. We discuss in particular scenarios with long-lived massive states decaying into dark matter particles, appearing naturally in several beyond the standard model theories, such as supergravity and superstring frameworks. Since non-thermal production favors dark matter candidates with large pair annihilation rates, we analyze the possible connection with the anomalies detected in the lepton cosmic-ray flux by Pamela and Fermi. Concentrating on supersymmetric models, we consider the effect of these non-standard cosmologies in selecting a preferred mass scale for the lightest supersymmetric particle as dark matter candidate, and the consequent impact on the interpretation of new physics discovered or excluded at the LHC. Finally, we examine a rather predictive model, the G2-MSSM, investigating some of the standard assumptions usually implemented in the solution of the Boltzmann equation for the dark matter component, including coannihilations. We question the hypothesis that kinetic equilibrium holds along the whole phase of dark matter generation, and the validity of the factorization usually implemented to rewrite the system of coupled Boltzmann equation for each coannihilating species as a single equation for the sum of all the number densities. As a byproduct we develop here a formalism to compute the kinetic decoupling temperature in case of coannihilating particles, which can be applied also to other particle physics frameworks, and also to standard thermal relics within a standard cosmology.

If dark matter particles pair annihilate into stable Standard Model particles, a population of energetic, non-thermal electrons and positrons is necessarily produced. The secondary radiation resulting from the energy losses of this non-standard population (including synchrotron and inverse Compton up-scattering of background radiation photons) in turn produces a peculiar multi-wavelength spectrum, extending from radio to gamma-ray frequencies. We give here an overview of multi-wavelength searches for dark matter, including the expected injection spectrum and production rate of electrons and positrons from dark matter annihilation, the computation of the effects of propagation and energy losses and the actual multi-wavelength emission. We then outline the application of this general framework to the case of galaxy clusters (specifically Coma, 1E 0657-56 - the so-called Bullet cluster - and Ophiuchus) and of dwarf spheroidal galaxies (including Draco, Fornax, Ursa Minor and Carina). We also review the application of multi-wavelength search strategies to our own Milky Way, and more specifically to the Galactic center environment, to dark galactic mini-halos and to the so-called WMAP haze and other radio data. We argue that multi-wavelength observations will complement gamma-ray observations as probes of particle dark matter, since the expected luminosities at different frequencies are generically comparable. The indirect search for dark matter with astronomical observations at various frequencies is therefore a crucial and indispensable element in the quest for the fundamental nature of dark matter.

We present a novel study on the problem of constructing mass models for the Milky Way, concentrating on features regarding the dark matter halo component. We have considered a variegated sample of dynamical observables for the Galaxy, including several results which have appeared recently, and studied a 7- or 8-dimensional parameter space - defining the Galaxy model - by implementing a Bayesian approach to the parameter estimation based on a Markov Chain Monte Carlo method. The main result of this analysis is a novel determination of the local dark matter halo density which, assuming spherical symmetry and either an Einasto or an NFW density profile is found to be around 0.39 GeV cm$^{-3}$ with a 1-$\sigma$ error bar of about 7%; more precisely we find a $\rho_{DM}(R_0) = 0.385 \pm 0.027 \rm GeV cm^{-3}$ for the Einasto profile and $\rho_{DM}(R_0) = 0.389 \pm 0.025 \rm GeV cm^{-3}$ for the NFW. This is in contrast to the standard assumption that $\rho_{DM}(R_0)$ is about 0.3 GeV cm$^{-3}$ with an uncertainty of a factor of 2 to 3. A very precise determination of the local halo density is very important for interpreting direct dark matter detection experiments. Indeed the results we produced, together with the recent accurate determination of the local circular velocity, should be very useful to considerably narrow astrophysical uncertainties on direct dark matter detection.

The local antimatter fluxes induced by an individual dark matter (DM) substructure can be significantly dependent on the proper motion of the source. We derive analytic solutions to the propagation equation for time-dependent positron and antiproton primary sources, finding that the static limit is a fair approximation only for very high energy particles and nearby sources. We discuss weakly interacting massive particle (WIMP) models fitting the PAMELA positron excess and the FERMI all-electron data. We show that, for a single non-static DM point-source, one cannot extract from the data, in a unique way, model independent particle physics observables, such as the WIMP mass, the pair annihilation cross section, and the annihilation yield. The gamma-ray emission associated to WIMP models inducing a significant local flux of positrons or antiprotons is found to be compatible with EGRET measurements, but it can be definitely singled out with the FERMI-LAT telescope.

We propose to test the dark matter (DM) interpretation of the positron excess observed by the PAMELA cosmic-ray (CR) detector through the identification of a Galactic diffuse gamma-ray component associated to DM-induced prompt and radiative emission. The goal is to present an analysis based on minimal sets of assumptions and extrapolations with respect to locally testable or measurable quantities. We discuss the differences between the spatial and spectral features for the DM-induced components (with an extended, possibly spherical, source function) and those for the standard CR contribution (with sources confined within the stellar disc), and propose to focus on intermediate and large latitudes. We address the dependence of the signal to background ratio on the model adopted to describe the propagation of charged CRs in the Galaxy, and find that, in general, the DM-induced signal can be detected by the Fermi Gamma-ray Space Telescope at energies above 100 GeV. An observational result in agreement with the prediction from standard CR components only, would imply very strong constraints on the DM interpretation of the PAMELA excess. On the other hand, if an excess in the diffuse emission above 100 GeV is identified, the angular profile for such emission would allow for a clean disentanglement between the DM interpretation and astrophysical explanations proposed for the PAMELA excess. We also compare to the radiative diffuse emission at lower frequencies, sketching in particular the detection prospects at infrared frequencies with the Planck satellite.

We perform a systematic study of the multi-wavelength signal induced by weakly interacting massive particle (WIMP) annihilations at the Galactic Center (GC). Referring to a generic WIMP dark matter (DM) scenario and depending on astrophysical inputs, we discuss spectral and angular features and sketch correlations among signals in the different energy bands. None of the components which have been associated to the GC source Sgr A*, nor the diffuse emission components from the GC region, have spectral or angular features typical of a DM source. Still, data-sets at all energy bands, namely, the radio, near infrared, X-ray and gamma-ray bands, contribute to place significant constraints on the WIMP parameter space. In general, the gamma-ray energy range is not the one with the largest signal to background ratio. In the case of large magnetic fields close to the GC, X-ray data give the tightest bounds. The emission in the radio-band, which is less model dependent, is very constraining as well. The recent detection by HESS of a GC gamma-ray source, and of a diffuse gamma-ray component, limits the possibility of a DM discovery with next generation of gamma-ray telescopes, like GLAST and CTA. We find that the most of the region in the parameter space accessible to these instruments is actually already excluded at other wave-lenghts. On the other hand, there may be still an open window to improve constraints with wide-field radio observations.

The generation of dark matter in late decays of quasi-stable massive particles has been proposed as a viable framework to address the excess of power found in numerical N-body simulations for cold dark matter cosmologies. We identify a convenient set of variables to illustrate which requirements need to be satisfied in any generic particle physics model to address the small scale problems and to fulfill other astrophysical constraints. As a result of this model-independent analysis, we point out that meeting these requirements in a completely natural way is inherently difficult. In particular, we re-examine the role of gravitinos and Kaluza-Klein gravitons in this context and find them disfavoured as a solution to the small-scale problems in case they are DM candidates generated in the decay of thermally produced WIMPs. We propose right-handed sneutrinos and right-handed Kaluza-Klein neutrinos as alternatives. We find that they are viable dark matter candidates, but that they can contribute to a solution of the small scale problems only in case the associated Dirac neutrino mass term appears as a subdominant contribution in the neutrino mass matrix.

We show that a recently constructed five-dimensional (5D) model with gauge-Higgs unification and explicit Lorentz symmetry breaking in the bulk, provides a natural dark matter candidate. This is the lightest Kaluza-Klein particle odd under a certain discrete Z_2 symmetry, which has been introduced to improve the naturalness of the model, and resembles KK-parity but is less constraining. The dark matter candidate is the first KK mode of a 5D gauge field and electroweak bounds force its mass above the TeV scale. Its pair annihilation rate is too small to guarantee the correct relic abundance; however coannihilations with colored particles greatly enhance the effective annihilation rate, leading to realistic relic densities.

The cluster 1ES0657-556 is an ideal astrophysical laboratory to study the distribution and the nature of Dark Matter because this last component is spatially separated from the intracluster gas. We show that microwave observations can provide crucial probes of Dark Matter in this system. We calculate the expected SZ effect from Dark Matter annihilation in the main mass concentrations of the cluster 1ES0657-556, and we estimate the sources of contamination, confusion and bias to asses its significance. We find that SZ observations at a frequency of 223 GHz can resolve both spatially and spectrally the SZ_DM signal and isolate it from the other SZ signals, and mainly from the thermal SZ effect which is null at frequencies 220-223 GHz for the case of 1ES0657-556. We conclude that SZ observations with sub-arcmin resolution and micro-K sensitivity of 1ES0657-556 are crucial, and maybe unique, to find direct astrophysical probes of the existence and of the nature of Dark Matter, or to set strong experimental limits.

We consider the split extended (N=2) supersymmetry scenario recently proposed by Antoniadis et al. [hep-ph/0507192] as a realistic low energy framework arising from intersecting brane models. While all scalar superpartners and charged gauginos are naturally at a heavy scale, the model low energy spectrum contains a Higgsino-like chargino and a neutralino sector made out of two Higgsino and two Bino states. We show that the lightest neutralino is a viable dark matter candidate, finding regions in the parameter space where its thermal relic abundance matches the latest determination of the density of matter in the Universe by WMAP. We also discuss dark matter detection strategies within this model: we point out that current data on cosmic-ray antimatter already place significant constraints on the model, while direct detection is the most promising technique for the future. Analogies and differences with respect to the standard split SUSY scenario based on the MSSM are illustrated.

We study the astrophysical implications of neutralino dark matter annihilations in galaxy clusters, with a specific application to the Coma cluster. We first address the determination of the dark halo models for Coma, starting from structure formation models and observational data, and we discuss in detail the role of sub-halos. We then perform a thorough analysis of the transport and diffusion properties of neutralino annihilation products, and investigate the resulting multi-frequency signals, from radio to gamma-ray frequencies. We also study other relevant astrophysical effects of neutralino annihilations, like the DM-induced Sunyaev-Zel'dovich effect and the intracluster gas heating. As for the particle physics setup, we adopt a two-fold approach, resorting both to model-independent bottom-up scenarios and to benchmark, GUT-motivated frameworks. We show that the Coma radio-halo data (the spectrum and the surface brightness) can be nicely fitted by the neutralino-induced signal for peculiar particle physics models and for magnetic field values, which we outline in detail. Fitting the radio data and moving to higher frequencies, we find that the multi-frequency spectral energy distributions are typically dim (with respect to the data) at EUV and X-ray frequencies, but show a non-negligible gamma-ray emission, depending on the amplitude of the Coma magnetic field. The best-fit particle physics models also produce a detectable SZ effect, but do not yield substantial heating of the intracluster gas in Coma. Due to the specific multi-frequency features of the DM-induced spectral energy distribution in Coma, we find that supersymmetric models can be significantly and optimally constrained either in the gamma-rays or at radio and microwave frequencies.

We explore the possible signatures of dark matter (DM) pair annihilations in the nearby dwarf spheroidal galaxy Draco. After investigating the mass models for Draco in the light of available observational data, we carefully model the DM density profile, taking advantage of numerical simulations of hierarchical structure formation. We then analyze the gamma-ray and electron/positron yield expected for weakly interacting DM particle (WIMP) models, including an accurate treatment of the propagation of the charged particle species. We show that unlike in larger DM structures - such as galaxy clusters - spatial diffusion plays here an important role. While Draco would appear as a point-like gamma-ray source, synchrotron emission from electrons and positrons produced by WIMP annihilations features a spatially extended structure. Depending upon the cosmic ray propagation setup and the size of the magnetic fields, the search for a diffuse radio emission from Draco can be a more sensitive indirect DM search probe than gamma rays. Finally, we show that available data are consistent with the presence of a black hole at the center of Draco: if this is indeed the case, very significant enhancements of the rates for gamma rays and other emissions related to DM annihilations are expected.

It has been recently shown that the electroweak baryogenesis mechanism is feasible in Standard Model extensions containing extra fermions with large Yukawa couplings. We show here that the lightest of these fermionic fields can naturally be a good candidate for cold dark matter. We find regions in the parameter space where the thermal relic abundance of this particle is compatible with the dark matter density of the Universe as determined by the WMAP experiment. We study direct and indirect dark matter detection for this model and compare with current experimental limits and prospects for upcoming experiments. We find, contrary to the standard lore, that indirect detection searches are more promising than direct detection ones, and they already exclude a portion of the parameter space.

We study the phenomenology of neutralino dark matter within generic supersymmetric scenarios where the Gaugino and Higgsino masses are much lighter than the scalar soft breaking masses (Split Supersymmetry). We consider a low-energy model-independent approach and show that the guidelines in the definition of this general framework come from cosmology, which forces the lightest neutralino to have a mass smaller than 2.2 TeV. The testability of the framework is addressed by discussing all viable dark matter detection techniques. Current data on cosmic rays antimatter, gamma-rays and on the abundance of primordial 6Li already set significant constraints on the parameter space. Complementarity among future direct detection experiments, indirect searches for antimatter and with neutrino telescopes, and tests of the theory at future accelerators, such as the LHC and a NLC, is highlighted. In particular, we study in detail the regimes of Wino-Higgsino mixing and Bino-Wino transition, which have been most often neglected in the past. We emphasize that our analysis may apply to more general supersymmetric models where scalar exchanges do not provide the dominant contribution to annihilation rates.

We show that charged-particles decaying in the early Universe can induce a scale-dependent or `running' spectral index in the small-scale linear and nonlinear matter power spectrum and discuss examples of this effect in minimal supersymmetric models in which the lightest neutralino is a viable cold-dark-matter candidate. We find configurations in which the neutralino relic density is set by coannihilations with a long-lived stau, and the late decay of staus partially suppresses the linear matter power spectrum. Nonlinear evolution on small scales then causes the modified linear power spectrum to evolve to a nonlinear power spectrum similar (but different in detail) to models parametrized by a constant running $\alpha_{s}=d n_{s}/d{\rm ln} k$ by redshifts of 2 to 4. Thus, Lyman-$\alpha$ forest observations, which probe the matter power spectrum at these redshifts, might not discriminate between the two effects. However, a measurement of the angular power spectrum of primordial 21-cm radiation from redshift $z \approx 30$--$200$ might distinguish between this charged-decay model and a primordial running spectral index. The direct production of a long-lived charged particle at future colliders is a dramatic prediction of this model.

We analyze the antimatter yield of supersymmetric (SUSY) models with large neutralino annihilation cross sections. We introduce three benchmark scenarios, respectively featuring bino, wino and higgsino-like lightest neutralinos, and we study in detail the resulting antimatter spectral features. We carry out a systematic and transparent comparison between current and future prospects for direct detection, neutrino telescopes and antimatter searches. We demonstrate that often, in the models we consider, antimatter searches are the only detection channel which already constrains the SUSY parameter space. Particularly large antiprotons fluxes are expected for wino-like lightest neutralinos, while significant antideuteron fluxes result from resonantly annihilating binos. We introduce a simple and general recipe which allows to assess the visibility of a given SUSY model at future antimatter search facilities. We provide evidence that upcoming space-based experiments, like PAMELA or AMS, are going to be, in many cases, the unique open road towards dark matter discovery.

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.

The thermal relic abundance of species critically depends on the assumed underlying cosmological model. In the case of neutralinos, freeze-out takes place long before Big Bang Nucleosynthesis, which provides the strongest constraint on the evolution of the Hubble parameter in the Early Universe. We show that non-standard cosmologies, such as models featuring a quintessential scalar field or primordial anisotropies, can lead to large enhancements in the neutralino relic abundance, up to six orders of magnitudes. Within these scenarios, supersymmetric models with large neutralino annihilation cross sections may account for the whole inferred amount of cold dark matter, yielding on the other hand large indirect detection rates.

The EGRET telescope has identified a gamma-ray source at the Galactic center. We point out here that the spectral features of this source are compatible with the gamma-ray flux induced by pair annihilations of dark matter weakly interacting massive particles (WIMPs). We show that the discrimination between this interpretation and other viable explanations will be possible with GLAST, the next major gamma-ray telescope in space, on the basis of both the spectral and the angular signature of the WIMP-induced component. If, on the other hand, the data will point to an alternative explanation, we prove that there will still be the possibility for GLAST to single out a weaker dark matter source at the Galactic center. The potential of GLAST has been explored both in the context of a generic simplified toy-model for WIMP dark matter, and in a more specific setup, the case of dark matter neutralinos in the minimal supergravity framework. In the latter, we find that even in the case of moderate dark matter densities in the Galactic center region, there are portions of the parameter space which will be probed by GLAST.

We investigate the enhancement of neutralino relic density in the context of a realistic cosmological scenario with Quintessence. The accurate relic density computation we perform allows us to be sensitive to both cases with shifts in the abundance at the percent level, and to enhancements as large as 10^6. We thoroughly analyze the dependence on the supersymmetric spectrum and on the mass and composition of the lightest neutralino. We point out that supersymmetric models yielding a wino or higgsino-like lightest neutralino become cosmologically appealing in the presence of Quintessence.

The detection of exotic cosmic rays due to pair annihilation of dark matter particles in the Milky Way halo is a viable techniques to search for supersymmetric dark matter candidates. The study of the spectrum of gamma-rays, antiprotons and positrons offers good possibilities to perform this search in a significant portion of the Minimal Supersymmetric Standard Model parameter space. In particular the EGRET team have seen a convincing signal for a strong excess of emission from the Galactic center that has no simple explanation with standard processes. We will review the limits achievable with the experiment GLAST taking into accounts the LEP results and we will compare this method with the antiproton and positrons experiments.

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.

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.

This paper reviews the field of gamma-ray astronomy and describes future experiments and prospects for advances in fundamental physics and high-energy astrophysics through gamma-ray measurements. We concentrate on recent progress in the understanding of active galaxies, and the use of these sources as probes of intergalactic space. We also describe prospects for future experiments in a number of areas of fundamental physics, including: searches for an annihilation line from neutralino dark matter, understanding the energetics of supermassive black holes, using AGNs as cosmological probes of the primordial radiation fields, constraints on quantum gravity, detection of a new spectral component from GRBs, and the prospects for detecting primordial black holes.

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.

We propose a novel antiparticle detector. The gaseous antiparticle spectrometer (GAPS) effects particle identification through the characteristic X-rays emitted by antiparticles when they form exotic atoms in gases. GAPS obtains particularly high grasp (effective area--solid angle product) at lower particle energies, where conventional schemes are most limited in their utility. The concept is simple and lightweight, so it can be readily employed on balloon and space-based missions. An extremely powerful potential application of GAPS is a space-based search for the neutralino through the detection of a neutralino annihilation by-product -- the antideuteron. Paradoxically, this space-based search for the neutralino is capable of achieving comparable sensitivity to as yet unrealized third generation, underground dark matter experiments. And GAPS can obtain this performance in a very modest satellite experiment. GAPS can also provide superior performance in searches for primary antiprotons produced via neutralino annihilation, black hole evaporation and in probing subdominant contributions to the antiproton flux at low energies. In a deep space mission GAPS will obtain higher sensitivity for a given weight and power than BGO calorimeters.

We consider a model for neutralino dark matter candidates arising in anomaly-mediated supersymmetry breaking schemes, and examine its testability through the search for exotic cosmic rays produced by neutralino pair annihilations in the dark halo of the Galaxy. We find that the model is already constrained by available antiprotons and positrons measurements and may be further tested in upcoming measurements of these cosmic ray species. We show also that the monochromatic gamma-ray flux from neutralino annihilations is enhanced in this model up to two orders of magnitude with respect to alternative scenarios. The gamma-ray flux detected by the Energetic Gamma Ray Experiment Telescope in the direction of the Galactic center exceeds significantly the theoretical expectation of standard emission models. We prove that if at least 10% of this excess is due to gamma-ray radiation with continuum energy spectrum from neutralino annihilations in the model under investigation, the associated gamma-ray line will be detected by upcoming gamma-ray experiments.

The past growth of the central black hole (BH) might have enhanced the density of cold dark matter halo particles at the Galactic center. We compute this effect in realistic growth models of the present (2-3)*10**6 solar mass BH from a low-mass seed BH, with special attention to dynamical modeling in a realistic galaxy environment with merger and orbital decay of a seed BH formed generally outside the exact center of the halo. An intriguing ``very-dense spike'' of dark matter has been claimed in models of Gondolo and Silk with density high enough to contradict with experimental upper bounds of neutralino annihilation radiation. This ``spike'' disappears completely or is greatly weakened when we include important dynamical processes neglected in their idealized/restrictive picture with cold particles surrounding an at-the-center zero-seed adiabaticly-growing BH. For the seed BH to spiral in and settle to the center within a Hubble time by dynamical friction, the seed mass must be at least a significant fraction of the present BH. Any subsequent at-the-center growth of the BH and steepening of the central Keplerian potential well can squeeze the halo density distribution only mildly, whether the squeezing happens adiabatically or instantaneously.

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.

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.

We investigate whether the annual modulation observed in the DAMA experiment can be due to a weakly-interacting massive particle (WIMP) with an axial-vector (spin-dependent; SD) coupling to nuclei. We evaluate the SD WIMP-proton cross section under the assumption that such scattering accounts for the DAMA modulation, and we do the same for a SD WIMP-neutron cross section. We show that SD WIMP-proton scattering is ruled out in a model-independent fashion by null searches for energetic neutrinos from WIMP annihilation in the Sun, and that SD WIMP-neutron scattering is ruled out for WIMP masses > 20 GeV by the null result with the DAMA Xe detector. A SD WIMP with mass < 20 GeV is still compatible, but only if the SD WIMP-neutron interaction is four orders of magnitude greater than the WIMP-proton interaction.

An annual modulation in the event rate of the NaI detector of the DAMA collaboration has been used to infer the existence of particle dark matter in the Galactic halo. Bounds on the WIMP mass and WIMP-nucleon cross section have been derived. These analyses have assumed that the local dark-matter velocity distribution is either isotropic or has some bulk rotation. Here we consider the effects of possible structure in the WIMP velocity distribution on the annual-modulation amplitude. We show that if we allow for a locally anisotropic velocity dispersion tensor, the interpretation of direct detection experiments could be altered significantly. We also show that uncertainties in the velocity distribution that arise from uncertainties in the radial density profile are less important if the velocity dispersion is assumed to be isotropic.

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.

More than a decade ago it was noticed that an unexpectedly large value of the measured cosmic antiproton flux at low kinetic energies could be interpreted in terms of a neutralino-induced component. An overproduction of low energy antiprotons seems however to be disfavoured by recent data from the BESS experiment. The strategy we propose here is to focus instead on the high energy antiproton spectrum. We find cases in which the signal from neutralino annihilations in a clumpy halo scenario can be unambiguously distinguished from the cosmic ray induced component. Such signatures can be tested in the near future by upcoming experiments. We investigate also the possibility that antiprotons have a finite lifetime. We find cases in which the neutralino-induced flux is consistent with present data for a lifetime as low as 0.15 Myr, in clear violation of a recent bound derived from the antiproton cosmic ray flux.

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.

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.

Recent advances in N-body simulations of cold dark matter halos point to a substantial density enhancement near the center. This means that, e.g., the $\gamma$ ray signals from neutralino dark matter annihilations would be significantly enhanced compared to old estimates based on an isothermal sphere model with large core radius. Another important development concerns new detectors, both space- and ground-based, which will cover the window between 50 and 300 GeV where presently no cosmic $\gamma$-ray data are available. Thirdly, new calculations of the $\gamma$-ray line signal (a sharp spike of $10^{-3}$ relative width) from neutralino annihilations have revealed a hitherto neglected contribution which, for heavy higgsino-like neutralinos, gives an annihilation rate an order of magnitude larger than previously predicted. We make a detailed phenomenological study of the possible detection rates given these three pieces of new information. We show that the proposed upgrade of the Whipple telescope will make it sensitive to a region of parameter space, with substantial improvements possible with the planned new generation of Air Cherenkov Telescope Arrays. We also comment on the potential of the GLAST satellite detector. An evaluation of the continuum $\gamma$-rays produced in neutralino annihilations into the main modes is also done. It is shown that a combination of high-rate models and very peaked halo models are already severely constrained by existing data.

A full one-loop calculation of neutralino S-wave annihilation into the $Z\gamma$ final state is performed in the minimal supersymmetric extension of the Standard Model. This process, like the similar one with two photons in the final state, may be of importance for the indirect detection of supersymmetric dark matter through the very narrow $\gamma$ ray line that would result from neutralino annihilations in the galactic halo. We give the complete analytical formulas for this loop-induced process and treat the case of a pure Higgsino as a first application of our expressions. Predictions for the gamma line flux are given for the halo model which is of the form suggested by Kravtsov et al. and for the profile proposed by Navarro, Frenk and White. For heavy neutralinos, the lines from $2\gamma$ and $Z\gamma$ would have indistiguishable energy in a realistic detector, making the fluxes add and facilitating discovery. For lighter neutralinos, the positions and relative strengths of the two lines would give valuable information on the nature of the supersymmetric dark matter particles.

For the first time, a full one-loop calculation of the process $\chi\chi\to 2\gamma$ is performed, where $\chi$ is the lightest neutralino in the minimal supersymmetric extension of the Standard Model. This process is of interest for dark matter detection, since it would give a sharp $\gamma$ ray line with $E_\gamma=m_\chi$. We improve upon and correct published formulas, and give cross sections for supersymmetric models with $\chi$ masses between 30 GeV and several TeV. We find a new contribution, previously neglected, which enhances the $2\gamma$ rate for TeV higgsinos by up to an order of magnitude. As a byproduct, we obtain a new expression for the related process $\chi\chi\to 2 gluons$, which on the other hand is generally smaller than previously calculated. There has been a recent claim that evidence for a 3.5 TeV higgsino annihilating into a $\gamma$ line may already exist from balloon emulsion and Air Cherenkov Telescope data. We comment on attractive features and problems with this interpretation.