A measurement of the dipole anisotropy in galactic coordinates for different charged cosmic rays has been performed with the Alpha Magnetic Spectrometer (AMS) onboard the International Space Station (ISS). Results are presented for the first 7.5 years of data taking for protons, Helium, Carbon and Oxygen, and 6.5 years for positrons and electrons. All the species are found to be consistent with isotropy and upper limits to the dipole amplitude have been computed. In particular, for energies above 16 GeV a limit of delta < 1.9% and delta < 0.5% at the 95% C.I. is found for positrons and electrons respectively. For rigidities above 200 GV a limit of delta < 0.38%, delta < 0.36%, delta < 1.9% and delta < 1.7% is obtained for protons, Helium, Carbon and Oxygen.

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.

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.

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.

The ROSAT Galactic wind observations confirm that our Galaxy launches supernova driven Galactic winds with wind speeds of about 150 km/s in the Galactic plane. Galactic winds of this strength are incompatible with current isotropic models for Cosmic Ray transport. In order to reproduce our local CRs in the presence of Galactic winds, charged CRs are required to be much more localized than in the standard isotropic GALPROP models. This requires that anisotropic diffusion is the dominant diffusion mode in the interstellar medium, particularly that the diffusion in the disk and in the halo are different. In addition small scale phenomena such as trapping by molecular cloud complexes and the structure of our local environment might influence the secondary CR production rate and our local CR density gradients. We introduce an anisotropic convection driven transport model (aCDM) which is consistent with the Galactic wind observations by ROSAT. This also explains the large bulge/disk ratio as observed by INTEGRAL. Furthermore such models predict an increase in the $e^+/(e^++e^-)$-fraction as observed by PAMELA and HEAT, if the synchrotron constraints in the 408 MHz and WMAP range are taken into account. No additional contribution from Dark Matter is required. The aCDM is able to explain the absence of a positron annihilation signal from molecular clouds as observed by INTEGRAL by virtue of a mechanism which confines and isotropizes CRs between MCs. We find that the EGRET excess of diffuse $\gamma$-rays currently cannot be explained by astrophysical effects in this type of model and that the interpretation of the EGRET excess as Dark Matter annihilation is perfectly consistent with all observational constraints from local CR fluxes and synchrotron radiation.

The EGRET excess of diffuse Galactic gamma rays shows all the features expected from dark matter annihilation (DMA): a spectral shape given by the fragmentation of mono-energetic quarks, which is the same in all sky directions and an intensity distribution of the excess expected from a standard dark matter halo, predicted by the rotation curve. From the EGRET excess one can predict the flux of antiprotons from DMA. However, how many antiprotons arrive at the detector strongly depends on the propagation model. The conventional isotropic propagation models trap the antiprotons in the Galaxy leading to a local antiproton flux far above the observed flux. According to Bergstrom et. al. this excludes the DMA interpretation of the EGRET excess. Here it is shown that more realistic anisotropic propagation models, in which most antiprotons escape by fast transport in the z-direction, are consistent with the B/C ratio, the antiproton flux and the EGRET excess from DMA.

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.

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.