Most cosmic ray particles observed derive from the explosions of massive stars, which commonly produce stellar black holes in their supernova explosions. When two such black holes find themselves in a tight binary system they finally merge in a gigantic emission of gravitational waves, events that have now been detected. After an introduction (section 1) we introduce the basic concept (section 2): Cosmic rays from exploding massive stars with winds always show two cosmic ray components at the same time: (i) the weaker polar cap component only produced by Diffusive Shock Acceleration with a cut-off at the knee, and (ii) the stronger $4 \pi$ component with a down-turn to a steeper power-law spectrum at the knee, and a final cutoff at the ankle. In section 3 we use the Alpha Magnetic Spectrometer (AMS) data to differentiate these two cosmic ray spectral components. The ensuing secondary spectra can explain anti-protons, lower energy positrons, and other secondary particles. Triplet pair production may explain the higher energy positron AMS data. In section 4 we test this paradigm with a theory of injection based on a combined effect of first and second ionization potential; this reproduces the ratio of Cosmic Ray source abundances to source material abundances. In section 5 we interpret the compact radio source 41.9+58 in the starburst galaxy M82 as a recent binary black hole merger, with an accompanying gamma ray burst. This can also explain the Ultra High Energy Cosmic Ray (UHECR) data in the Northern sky. Thus, by studying the cosmic ray particles, their abundances at knee energies, and their spectra, we can learn about what drives these stars to produce the observed cosmic rays.

Up to now searches for Dark Matter (DM) detection have not been successful, either because our paradigm in how DM signals should look like are wrong or the detector sensitivity is still too low in spite of the large progress made in recent years. We discuss both possibilities starting with what we know about DM from cosmology and why Supersymmetry provides such an interesting paradigm for cosmology and particle physics in order to appreciate what it means to give up this paradigm. In addition, we compare the predicted cross sections for direct and indirect DM detection with observations with emphasis on the latest developments. Especially, we discuss the possible origins of the two hotly debated candidates for a DM annihilation signal, namely the positron excess and the Fermi GeV excess, which are unfortunately incompatible with each other and more mundane astrophysical explanations exist.

In the Next-to-Minimal-Supersymmetric-Standard-Model (NMSSM) the lightest supersymmetric particle (LSP) is a candidate for the dark matter (DM) in the universe. It is a mixture from the various gauginos and Higgsinos and can be bino-, Higgsino- or singlino-dominated. Singlino-dominated LSPs can have very low cross sections below the neutrino background from coherent neutrino scattering which is limiting the sensitivity of future direct DM search experiments. However, previous studies suggested that the combination of both, the spin-dependent (SD) and spin-independent (SI) searches are sensitive in complementary regions of parameter space, so considering both searches will allow to explore practically the whole parameter space of the NMSSM. In this letter, the different scenarios are investigated with a new scanning technique, which reveals that significant regions of the NMSSM parameter space cannot be explored, even if one considers both, SI and SD, searches.

The so-called "GeV-excess" of the diffuse Galactic gamma-ray emission is studied with a spectral template fit based on energy spectra. The spectral templates can be obtained in a data-driven way from the gamma-ray data, which avoids the use of emissivity models to subtract the standardbackground processes from the data. Instead, one can determine these backgrounds simultaneously with any "signals" in any sky direction, including the Galactic disk and the Galactic center. Using the spectral template fit two hypothesis of the "GeV-excess" were tested: the dark matter (DM) hypothesis assuming the excess is caused by DM annihilation and the molecular cloud (MC) hypothesis assuming the "GeV-excess" is related to a depletion of gamma-rays below 2 GeV, as is directly observed in the Central Molecular Zone (CMZ). Both hypotheses provide acceptable fits, if one considers a limited field-of-view centered within 20$^\circ$ around the Galactic center and applies cuts on the energy range and/or excludes low latitudes, cuts typically applied by the proponents of the DM hypothesis. However, if one considers the whole gamma-ray sky and includes gamma-ray energies up to 100 GeV we find that the MC hypothesis is preferred over the DM hypothesis for several reasons: i) The MC hypothesis provides significantly better fits; ii) The morphology of the "GeV-excess" follows the morphology of the CO-maps, a tracer of MCs, i.e. there exists a strong "GeV-excess" in the Galactic disk also at large longitudes; iii) The massive CMZ with a rectangular field-of-view of $l \times b = 3.5^{\circ} \times 0.5^{\circ}$ shows the maximum of the energy flux per log bin in the diffuse gamma-ray spectrum at 2 GeV, i.e. the "GeV-excess", already in the raw data without any analysis. The rectangular profile contradicts the spherical morphology expected for DM annihilation.

The "GeV-excess" of the diffuse gamma-rays in the halo is studied with a template fit based on energy spectra for each possible process of gamma-ray emission. Such a fit allows to determine the background and signal simultaneously, so the Galactic Disk can be included in the analysis. We find evidence that the "excess", characterised by a gamma-ray spectrum peaking at 2 GeV, is much stronger in Molecular Clouds in the disk than the "GeV-excess" observed up to now in the halo. The possible reason why the emissivity of Molecular Clouds peaks at 2 GeV are the energy losses and magnetic cutoffs inside MCs, thus depleting the low energy part of the CR spectra and shifting the maximum of the gamma-ray spectra to higher energies. This peaking of the emissivity in Molecular Clouds at 2 GeV was clearly observed from the spectrum of the Central Molecular Zone, which dominates the emission in the inner few degrees of the Galactic Centre. Although the spectrum of the Central Molecular Zone peaks at 2 GeV, it cannot be responsible for the "GeV-excess" observed in the halo, since the latitude extension of the Zone is below $\rm |b|<0.5^\circ$. However, lines-of-sight into the halo cross Molecular Clouds in the disk, so the emissivity of clouds in the disk will be observed in the halo as an apparent "GeV-excess". The fact, that this "GeV-excess" has the same morphology in the disk {\it and} in the halo as the column density of Molecular Clouds, as traced by the CO map from the Planck satellite resembling an NFW-like latitude profile, and the fact the MCs have an emissivity peaking at 2 GeV shows that the "GeV-excess" originates from Molecular Clouds in the disk, not from a process surrounding the Galactic Centre.

The Fermi-LAT data revealed giant bubbles of emission above and below the Galactic plane with an energy spectrum significantly harder than seen from other directions. How the bubbles connect to the Galactic plane is unclear. We find that bubble-like emission is not only found in the halo, but it is strongly present in the Galactic plane as well with a morphology close to the spatial distribution of the 1.8 MeV gamma-ray line from 26Al, a radioactive nucleus synthesized in SNRs. In addition, the spectral shape of this hard component coincides with the predicted spectrum from cosmic rays trapped in sources. Hence, we propose that the bubble-like emission in the plane has a hadronic origin, which arises from SCRs. The bubbles in the halo have the same energy spectrum, which suggests that they are outflows from the plane with the gamma-rays arising from hadronic interactions of protons trapped in a plasma of advected gas. Towards the Galactic center (GC) we observe the 1-3 GeV excess, but observe this excess in addition in all regions where there is strong 26Al production. This excludes the dark matter annihilation interpretation. Instead, we propose that the excess is caused by a deficit of the CRs at low rigidities in the regions of dense molecular clouds, which are characterized by higher energy losses and stronger stellar winds. The 26Al line is a tracer of such regions. If we introduce a break in the cosmic ray injection spectra in the regions of strong 26Al production to effectively take the deficit of low rigidity CRs into account, the excess in the Galactic plane and towards the GC both disappear. We thus conclude that the excess in the Galactic center is largely an artefact from the excess in the Galactic plane. The correlation between the two originates simply from the lines-of-sight towards the GC crossing the Galactic plane.

The diffuse gamma-ray sky revealed 'Bubbles' of emission above and below the Galactic Plane symmetric around the centre of the Milky Way with a height of 10 kpc in both directions. At present there is no convincing explanation for the origin. To understand the role of the Galactic Centre (GC) one has to study the Bubble spectrum inside the disc, a region which has been excluded from previous analysis because of the large foreground. From a novel template fit, which allows a simultaneous determination of the signal and foreground in any direction, we find that bubble-like emission is not only found in the halo, but in the Galactic Plane as well with a width in latitude coinciding with the molecular clouds. The longitude distribution has a width corresponding to the Galactic Bar with an additional contribution from the Scutum-Centaurus arm. The energy spectrum of the Bubbles coincides with the predicted contribution from CRs trapped in sources (SCRs). Also the energetics fits well. Hence, we conclude that the bubble-like emission has a hadronic origin, which arises from SCRs and the Bubbles in the halo arise from hadronic interactions in advected gas. Evidence for advection is provided by the ROSAT X-rays from hot gas in the Bubble region.

We highlight the differences of the dark matter sector between the constrained minimal supersymmetric SM (CMSSM) and the next-to-minimal supersymmetric SM (NMSSM) including the 126 GeV Higgs boson using GUT scale parameters. In the dark matter sector the two models are quite orthogonal: in the CMSSM the WIMP is largely a bino and requires large masses from the LHC constraints. In the NMSSM the WIMP has a large singlino component and is therefore independent of the LHC SUSY mass limits. The light NMSSM neutralino mass range is of interest for the hints concerning light WIMPs in the Fermi data. Such low mass WIMPs cannot be explained in the CMSSM. Furthermore, prospects for discovery of XENON1T and LHC at 14 TeV are given.

The recent discovery of a Higgs-like boson at the LHC with a mass of 126 GeV has revived the interest in supersymmetric models, which predicted a Higgs boson mass below 130 GeV long before its discovery. We compare systematically the allowed parameter space in the constrained Minimal Supersymmetric Standard Model (CMSSM) and the Next-to-Minimal Supersymmetric Model (NMSSM) by minimizing the chi^2 function with respect to all known constraints from accelerators and cosmology using GUT scale parameters. For the CMSSM the Higgs boson mass at tree level is below the Z^0 boson mass and large radiative corrections are needed to obtain a Higgs boson mass of 126 GeV, which requires stop squark masses in the multi-TeV range. In contrast, for the NMSSM light stop quarks are allowed, since in the NMSSM at tree level the Higgs boson mass can be above the Z^0 boson mass from mixing with the additional singlet Higgs boson. Predictions for the scalar boson masses are given in both models with emphasis on the unique signatures of the NMSSM, where the heaviest scalar Higgs boson decays in the two lighter scalar Higgs bosons with a significant branching ratio, in which case one should observe double Higgs boson production at the LHC. Such a signal is strongly suppressed in the CMSSM. In addition, since the LSP is higgsino-like, Higgs boson decays into LSPs can be appreciable, thus leading to invisible Higgs decays.

The ATLAS and CMS experiments did not find evidence for Supersymmetry using close to 5/fb of published LHC data at a center-of-mass energy of 7 TeV. We combine these LHC data with data on B_s -> mu mu (LHCb experiment), the relic density (WMAP and other cosmological data) and upper limits on the dark matter scattering cross sections on nuclei (XENON100 data). The excluded regions in the constrained Minimal Supersymmetric SM (CMSSM) lead to gluinos excluded below 1270 GeV and dark matter candidates below 220 GeV for values of the scalar masses (m_0) below 1500 GeV. For large m_0 values the limits of the gluinos and the dark matter candidate are reduced to 970 GeV and 130 GeV, respectively. If a Higgs mass of 125 GeV is imposed in the fit, the preferred SUSY region is above this excluded region, but the size of the preferred region is strongly dependent on the assumed theoretical error.

The direct searches for Superymmetry at colliders can be complemented by direct searches for dark matter (DM) in underground experiments, if one assumes the Lightest Supersymmetric Particle (LSP) provides the dark matter of the universe. It will be shown that within the Constrained minimal Supersymmetric Model (CMSSM) the direct searches for DM are complementary to direct LHC searches for SUSY and Higgs particles using analytical formulae. A combined excluded region from LHC, WMAP and XENON100 will be provided, showing that within the CMSSM gluinos below 1 TeV and LSP masses below 160 GeV are excluded (m_{1/2} > 400 GeV) independent of the squark masses.

The pure leptonic decay B_s -> mu mu is strongly suppressed in the Standard Model (SM), but can have large enhancements in Supersymmetry, especially at large values of tanbe. New limits on this decay channel from recent LHC data have been used to claim that these limits restrict the SUSY parameter space even more than the direct searches. However, direct searches are hardly dependent on tanbe, while BR(B_s -> mu mu) is proportional to tanbe^6. The relic density constraint requires large tanbe in a large region of the parameter space, which can lead to large values of B_s -> mu mu. Nevertheless, the experimental upper limit on BR(B_s -> mu mu) is not constraining the parameter space of the CMSSM more than the direct searches and the present Higgs limits, if combined with the relic density. We also observe SUSY parameter regions with negative interferences, where the B_s -> mu mu value is up to a factor three below the SM expectation, even at large values of tanbe.

Both the gas flaring and the dip in the rotation curve, which was recently reconfirmed with precise measurements using the VERA VLBI array in Japan, suggest doughnut-like substructure in the dark matter (DM) halo. A global fit to all available data shows that the data are indeed best described by an NFW DM profile complemented by two doughnut-like DM substructures with radii of 4.2 and 12.4 kpc, which coincide with the local dust ring and the Monocerus ring of stars, respectively. Both regions have been suggested as regions with tidal streams from "shredded" satellites. If real, the radial extensions of these nearby ringlike structures enhance the local dark matter density by a factor of four to about 1.3$\pm0.3$ GeV/cm$^3$. It is shown that i) this higher DM density is perfectly consistent with the local gravitational potential determining the surface density and the local matter density (Oort limit), ii) previous determinations of the surface density were biased by the assumption of a smoothly varying DM halo and iii) the s-shaped gas flaring is explained. Such a possible enhancement of the local DM density is of great interest for direct DM searches and would change the directional dependence for indirect DM searches.

The rotation curve, the total mass and the gravitational potential of the Galaxy are sensitive measurements of the dark matter halo profile. In this publication cuspy and cored DM halo profiles are analysed with respect to recent astronomical constraints in order to constrain the shape of the Galactic DM halo and the local DM density. All Galactic density components (luminous matter and DM) are parametrized. Then the total density distribution is constrained by astronomical observations: 1) the total mass of the Galaxy, 2) the total matter density at the position of the Sun, 3) the surface density of the visible matter, 4) the surface density of the total matter in the vicinity of the Sun, 5) the rotation speed of the Sun and 6) the shape of the velocity distribution within and above the Galactic disc. The mass model of the Galaxy is mainly constrained by the local matter density (Oort limit), the rotation speed of the Sun and the total mass of the Galaxy from tracer stars in the halo. It is shown from a statistical chi^2 fit to all data that the local DM density is strongly positively (negatively) correlated with the scale length of the DM halo (baryonic disc). Since these scale lengths are poorly constrained the local DM density can vary from 0.2 to 0.4 GeV/cm^3 (0.005 - 0.01 M_sun/pc^3) for a spherical DM halo profile and allowing total Galaxy masses up to 2 * 10^12 M_sun. For oblate DM halos and dark matter discs, as predicted in recent N-body simulations, the local DM density can be increased significantly.

Recently, new data on antiprotons and positrons from PAMELA, e- + e+ spectra from ATIC, FERMI and HESS up to TeV energies all indicate deviations from expectations, which has caused an interesting mix of new explanations, ranging from background, standard astrophysical sources to signals from dark matter (DM) annihilation. Unfortunately, the excess in positrons is not matched with obvious excesses in antiprotons or gamma rays, so a new class of DM scenarios with leptophilic WIMP candidates have been invoked. On the other hand, the increase in the positron fraction, which could have had any spectral shape for new physics, matches well the shape expected from proton background.

Isotropic diffusion models for Galactic cosmic ray transport put tight constraints on the maximum convection velocity in the halo. For a half halo height of 4 kpc the maximum convection speed is limited to 40 km/s in the halo, since otherwise the constraints from local secondary to primary ratios and radioactive isotopes cannot be met. The ROSAT Galactic wind observations of wind speeds up to 760 km/s therefore constitute a problem for diffusion models. It is shown that such wind speeds are possible, if the diffusion coefficient in the halo is different from the diffusion coefficient in the disk. The radial dependence of the wind velocity was taken to be proportional to the source strength, as expected from winds which are sustained by cosmic ray pressure. In this case the cosmic ray density and with it the diffuse $\gamma$-ray production from nuclear interactions are suppressed near the sources. This solves in a natural way the problem of the soft gradient in the radial dependence of the $\gamma$-ray flux. Furthermore, the large bulge over disk ratio in positron annihilation as observed by INTEGRAL, can be explained by positron escape from the disk in such a model.

Dark Matter annihilation (DMA) may yield an excess of gamma rays and antimatter particles, like antiprotons and positrons, above the background from cosmic ray interactions. The excess of diffuse Galactic Gamma Rays from EGRET shows all the features expected from DMA. The new precise measurements of the antiproton and positron fractions from PAMELA are compared with the EGRET excess. It is shown that the charged particles are strongly dependent on the propagation model used. The usual propagation models with isotropic propagation models are incompatible with the recently observed convection in our Galaxy. Convection leads to an order of magnitude uncertainty in the yield of charged particles from DMA, since even a rather small convection will let drift the charged particles in the halo to outer space. It is shown that such anisotropic propagation models including convection prefer a contribution from DMA for the antiprotons, but the rise in the positron fraction, as observed by PAMELA, is incompatible with DMA, if compared with the EGRET excess. A rise in the positron/electron ratio is expected, if the observed rise in the proton/electron ratio is carefully fitted in a propagation model, although the rise is slightly larger than expected, so additional local sources may contribute as well within the limited accuracy of the data.

Dark Matter annihilation (DMA) may yield an excess of gamma rays and antimatter particles, like antiprotons and positrons, above the background from cosmic ray interactions. Several signatures, ranging from the positron excess, as observed by HEAT, AMS-01 and PAMELA, the gamma ray excess, as observed by the EGRET spectrometer, the WMAP-haze, and constraints from antiprotons, as observed by CAPRICE, BESS and PAMELA, have been discussed in the literature. Unfortunately, the different signatures all lead to different WIMP masses, indicating that at least some of these interpretations are likely to be incorrect. Here we review them and discuss their relative merits and uncertainties. New x-ray data from ROSAT suggests non-negligible convection in our Galaxy, which leads to an order of magnitude uncertainty in the yield of charged particles from DMA, since even a rather small convection will let drift the charged particles in the halo to outer space.

The EGRET excess of diffuse Galactic gamma rays shows all the key features of dark matter annihilation (DMA) for a WIMP mass in the range 50-100 GeV, especially the distribution of the excess is compatible with a standard halo profile with some additional ringlike substructures at 4 and 13 kpc from the Galactic centre. These substructures coincide with the gravitational potential well expected from the ring of dust at 4 kpc and the tidal stream of dark matter from the Canis Major satellite galaxy at 13 kpc, as deduced from N-body simulations fitting to the Monoceros ring of stars. Strong independent support for this substructure is given by the gas flaring in our Galaxy. The gamma rays from DMA are originating predominantly from the hadronization of mono-energetic quarks, which should produce also a small, but known fraction of protons and antiprotons. Bergstrom et al. an antiproton flux far above the observed antiproton flux and they conclude that the DMA interpretation of the EGRET excess can therefore be excluded. However, they used an isotropic propagation model, i.e. the same diffusive propagation in the disk and the halo. It is shown that an anisotropic propagation model is consistent with the EGRET gamma ray excess, the antiproton flux and the ratios of secondary/primary and unstable/stable cosmic ray particles. Such an anisotropic propagation is supported by the large bulge/disk ratio of the positron annihilation line, as observed by the INTEGRAL satellite. In this case no need for new sources specific to the bulge are needed, so the claimed evidence for strong DMA in the bulge from these observations is strongly propagation model dependent.

The excess of diffuse galactic gamma rays above 1 GeV, as observed by the EGRET telescope on the NASA Compton Gamma Ray Observatory, shows all the key features from Dark Matter (DM) annihilation: (i) the energy spectrum of the excess is the same in all sky directions and is consistent with the gamma rays expected for the annihilation of WIMPs with a mass between 50-100 GeV; (ii) the intensity distribution of the excess in the sky is used to determine the halo profile, which was found to correspond to the usual profile from N-body simulations with additional substructure in the form of two doughnut-shaped structures at radii of 4 and 13 kpc; (iii) recent N-body simulations of the tidal disruption of the Canis Major dwarf galaxy show that it is a perfect progenitor of the ringlike Monoceros tidal stream of stars at 13 kpc with ring parameters in agreement with the EGRET data; (iiii) the mass of the outer ring is so large, that its gravitational effects influence both the gas flaring and the rotation curve of the Milky Way. Both effects are clearly observed in agreement with the DMA interpretation of the EGRET excess.

Some part of the relic Dark Matter is distributed in small-scale clumps which survived structure formation in inflation cosmological scenario. The annihilation of DM inside these clumps is a strong source of stable charged particles which can have a substantial density near the clump core. The streaming of the annihilation products from the clump can enhance irregularities in the galactic magnetic field. This can produce small scale variations in diffusion coefficient affecting propagation of Cosmic Rays.

The extragalactic background (EGB) of diffuse gamma rays can be determined by subtracting the Galactic contribution from the data. This requires a Galactic model (GM) and we include for the first time the contribution of dark matter annihilation (DMA), which was previously proposed as an explanation for the EGRET excess of diffuse Galactic gamma rays above 1 GeV. In this paper it is shown that the newly determined EGB shows a characteristic high energy bump on top of a steeply falling soft contribution. The bump is shown to be compatible with a contribution from an extragalactic DMA signal from weakly interacting massive particles (WIMPs) with a mass between 50 and 100 GeV in agreement with the EGRET excess of the Galactic diffuse gamma rays and in disagreement with earlier analysis. The remaining soft contribution of the EGB is shown to resemble the spectra of the observed point sources in our Galaxy.

The diffuse galactic EGRET gamma ray data show a clear excess for energies above 1 GeV in comparison with the expectations from conventional galactic models. This excess shows all the features expected from Dark Matter WIMP Annihilation: a)it is present and has the same spectrum in all sky directions, not just in the galactic plane, as expected for WIMP annihilation b) it shows an interesting substructure in the form of a doughnut shaped ring at 14 kpc from the centre of the galaxy, where a ring of stars indicated the probable infall of a dwarf galaxy. From the spectral shape of the excess the WIMP mass is estimated to be between 50 and 100 GeV, while from the intensity the halo profile is reconstructed, which is shown to explain the peculiar change of slope in the rotation curve at about 11 kpc (due to the ring of DM at 14 kpc). Recently it was claimed by Bergstrom et al. that the DM interpretation of the EGRET gamma ray excess is excluded by the antiproton fluxes, since in their propagation model with isotropic diffusion the flux of antiprotons would be far beyond the observed flux. However, the propagation can be largely anisotropic, because of the convection of particles perpendicular to the disc and inhomogeneities in the local environment. It is shown that anisotropic propagation can reduce the antiproton yield by an order of magnitude, while still being consistent with the B/C ratio. Therefore it is hard to use antiprotons to search for {\it light} DM particles, which yield a similar antiproton spectrum as the background, but the antiprotons are a perfect means to tune the many degenerate parameters in the propagation models.

Recently it was shown that the excess of diffuse Galactic gamma rays above 1 GeV traces the Dark Matter halo, as proven by reconstructing the peculiar shape of the rotation curve of our Galaxy from the gamma ray excess. This can be interpreted as a Dark Matter annihilation signal. In this paper we investigate if this interpretation is consistent with Supersymmetry. It is found that the EGRET excess combined with all electroweak constraints is fully consistent with the minimal mSUGRA model for scalars in the TeV range and gauginos below 500 GeV.

Elsaesser and Mannheim fit a contribution of Dark Matter Annihilation (DMA) to the extragalactic contribution of the galactic diffuse gamma ray flux, as deduced from the EGRET data by Strong, Moskalenko and Reimer.They find a WIMP mass of 515{+110}{-75} GeV and quote a systematic error of 30\%. However, they do not include large systematic uncertainties from the fact that the determination of the extragalactic flux (EGF) requires a model for the subtraction of the Galactic flux from the data.The data used were obtained with a model without Galactic DM, so one expects additional uncertainty in the region where DMA contributes. Including a Galactic DMA contribution reduces the significance and the WIMP mass. The latter then becomes compatible with the Galactic excess of diffuse gamma rays, which posseses all the properties of DMA with a much higher significance than the extragalactic excess.

The public data from the EGRET space telescope on diffuse Galactic gamma rays in the energy range from 0.1 to 10 GeV are reanalyzed with the purpose of searching for signals of Dark Matter annihilation (DMA). The analysis confirms the previously observed excess for energies above 1 GeV in comparison with the expectations from conventional Galactic models. In addition, the excess was found to show all the key features of a signal from Dark Matter Annihilation (DMA): a) the excess is observable in all sky directions and has the same shape everywhere, thus pointing to a common source; b) the shape corresponds to the expected spectrum of the annihilation of non-relativistic massive particles into - among others - neutral $\pi^0$ mesons, which decay into photons. From the energy spectrum of the excess we deduce a WIMP mass between 50 and 100 GeV, while from the intensity of the excess in all sky directions the shape of the halo could be reconstructed. The DM halo is consistent with an almost spherical isothermal profile with substructure in the Galactic plane in the form of toroidal rings at 4 and 14 kpc from the center. These rings lead to a peculiar shape of the rotation curve, in agreement with the data, which proves that the EGRET excess traces the Dark Matter.

The EGRET excess in the diffuse galactic gamma ray data above 1 GeV shows all features expected from Dark Matter WIMP Annihilation: a)it is present and has the same spectrum in all sky directions, not just in the galactic plane. b) The intensity of the excess shows the $1/r^2$ profile expected for a flat rotation curve outside the galactic disc with additionally an interesting substructure in the disc in the form of a doughnut shaped ring at 14 kpc from the centre of the galaxy. At this radius a ring of stars indicates the probable infall of a dwarf galaxy, which can explain the increase in DM density. From the spectral shape of the excess the WIMP mass is estimated to be between 50 and 100 GeV, while from the intensity the halo profile is reconstructed. Given the mass and intensity of the WIMPs the mass of the ring can be calculated, which is shown to explain the peculiar change of slope in the rotation curve at about 11 kpc. These results are model independent in the sense that only the known shapes of signal and background were fitted with free normalization factors, thus being independent of the model dependent flux caluclations. The statistical significance is more than 10 sigma in comparison with a fit of the conventional galactic model to the EGRET data. These signals of Dark Matter Annihilation are compatible with Supersymmetry including all electroweak constraints. The statistical significance combined with all features mentioned above provide an intriguing hint that the EGRET excess is indeed a signal from Dark Matter Annihilation.

The EGRET excess in the diffuse galactic gamma ray data above 1 GeV shows all the features expected from Dark Matter WIMP Annihilation: a)it is present and has the same spectrum in all sky directions, not just in the galactic plane. b) The intensity of the excess shows the $1/r^2$ profile expected for a flat rotation curve outside the galactic disc with additionally an interesting substructure in the disc in the form of a doughnut shaped ring at 14 kpc from the centre of the galaxy. At this radius a ring of stars indicates the probable infall of a dwarf galaxy, which can explain the increase in DM density. From the spectral shape of the excess the WIMP mass is estimated to be between 50 and 100 GeV, while from the intensity the halo profile is reconstructed. Given the mass and intensity of the WIMP's the mass of the ring can be calculated, which is shown to explain the peculiar change of slope in the rotation curve at about 11 kpc. These signals of Dark Matter Annihilation are compatible with Supersymmetry and have a statistical significance of more than 10 stnadard deviations in comparison with a fit of the conventional galactic model to the EGRET data. The statistical significance combined with all features mentioned above provide an intriguing hint that the EGRET excess is indeed a signal from Dark Matter Annihilation.

The diffuse galactic EGRET gamma ray data show a clear excess for energies above 1 GeV in comparison with the expectations from conventional galactic models. The excess is seen with the same spectrum in all sky directions, as expected for Dark Matter (DM) annihilation. This hypothesis is investigated in detail. The energy spectrum of the excess is used to limit the WIMP mass to the 50-100 GeV range, while the skymaps are used to determine the halo structure, which is consistent with a triaxial isothermal halo with additional enhancement of Dark Matter in the disc. The latter is strongly correlated with the ring of stars around our galaxy at a distance of 14 kpc, thought to originate from the tidal disruption of a dwarf galaxy. It is shown that this ring of DM with a mass of $\approx 2\cdot 10^{11} M_\odot$ causes the mysterious change of slope in the rotation curve at $R=1.1R_0$ and the large local surface density of the disc. The total mass of the halo is determined to be $3\cdot 10^{12} M_\odot$. A cuspy profile is definitely excluded to describe the gamma ray data. These signals of Dark Matter Annihilation are compatible with Supersymmetry for boost factors of 20 upwards and have a statistical significance of more than $10\sigma$ in comparison with the conventional galactic model. The latter combined with all features mentioned above provides an intriguing hint that the EGRET excess is indeed indirect evidence for Dark Matter Annihilation.

The diffuse galactic EGRET gamma ray data show a clear excess for energies above 1 GeV in comparison with the expectations from conventional galactic models. The excess is seen with the same spectrum in all sky directions, as expected for Dark Matter (DM) annihilation. This hypothesis is investigated in detail. The energy spectrum of the excess is used to limit the WIMP mass to the 50-100 GeV range, while the sky maps are used to determine the halo structure, which is consistent with a triaxial isothermal halo with additional enhancement of Dark Matter in the disc. The latter is strongly correlated with the ring of stars around our galaxy at a distance of 14 kpc, thought to originate from the tidal disruption of a dwarf galaxy. It is shown that this ring of DM with a mass of $\approx 2\cdot 10^{11} M_\odot$ causes the mysterious change of slope in the rotation curve at $R=1.1R_0$ and the large local surface density of the disc. The total mass of the halo is determined to be $3\cdot 10^{12} M_\odot$. A cuspy profile is definitely excluded to describe the gamma ray data. These signals of Dark Matter Annihilation are compatible with Supersymmetry for boost factors of 20 upwards and have a statistical significance of more than $10\sigma$ in comparison with the conventional galactic model. The latter combined with all features mentioned above provides an intriguing hint that the EGRET excess is indeed a signal from Dark Matter Annihilation.

A fit to the present data on the cosmic ray positron fraction can be considerably improved,if in addition to the positron production by nuclear interactions in the universe the possible contribution from supersymmetric dark matter annihilation is taken into account. We scan over the complete SUSY parameter space of the Constrained Minimal Supersymmetric Model (CMSSM) and find that in the acceptable regions the neutralino annihilation into b-bbar quark pairs is the dominant channel with hard positrons emerging from the semileptonic decays of the B-mesons.

We investigate, if the cosmic ray positron fraction, as reported by the HEAT and AMS collaborations, is compatible with the annihilation of neutralinos in the supergravity inspired Constrained Minimal Supersymmetric Model (CMSSM), thus complementing previous investigations, which did not consider constraints from unification, electroweak symmetry breaking and the present Higgs limits at LEP. We perform a global fit of all SUSY parameters to low energy constraints and find that in the acceptable parameter regions the neutralino annihilation into b-bbar quark pairs is the dominant channel and improves the fit to the experimental positron fraction data considerably with respect to a fit with background only. These fits are comparable to the fit for regions of parameter space, where the annihilation into W^+ W^- pairs dominates. However, these latter regions are ruled out by the present Higgs limit of 114 GeV from LEP in combination with b->s gamma data.