The search for antideuterons in cosmic rays has been proposed as a promising channel for dark matter indirect detection, especially for dark matter particles with a low or intermediate mass. With the current operational phase of the AMS-02 experiment and the ongoing development of a future dedicated experiment, the General Antiparticle Spectrometer (GAPS), there are exciting prospects for a dark matter detection in the near future. In this paper we develop a detailed and complete re-analysis of the cosmic-ray antideuteron signal, by discussing the main relevant issues related to antideuteron production and propagation through the interstellar medium and the heliosphere. In particular, we first critically revisit the coalescence mechanism for antideuteron production in dark matter annihilation processes. Then, since antideuteron searches have their best prospects of detection at low kinetic energies where the effect of the solar wind and magnetic field are most relevant, we address the impact of solar modulation modeling on the antideuteron flux at the Earth by developing a full numerical 4D solution of cosmic rays transport in the heliosphere. We finally use these improved predictions to provide updated estimates of the reaching capabilities for AMS-02 and GAPS, compatible with the current constraints imposed by the antiprotons measurements of PAMELA. After the antiproton bound is applied, prospects of detection of up to about 15 events in GAPS LDB+ and AMS-02 missions are found, depending on the dark matter mass, annihilation rate and production channel from one side, and on the coalescence process, galactic and solar transport parameters on the other.

We model the transport of cosmic ray nuclei in the Galaxy by means of a new numerical code. Differently from previous numerical models we account for a generic spatial distribution of the diffusion coefficient. We found that in the case of radially uniform diffusion, the main secondary/primary ratios (B/C, N/O and sub-Fe/Fe) and the modulated antiproton spectrum match consistently the available observations. Convection and re-acceleration do not seem to be required in the energy range we consider: $1 < E < 10^3$ GeV/nucleon. We generalize these results accounting for radial dependence of the diffusion coefficient, which is assumed to trace that of the cosmic ray sources. While this does not affect the prediction of secondary/primary ratios, the simulated longitude profile of the diffuse $\gamma$-ray emission is significantly different from the uniform case and may agree with EGRET measurements without invoking ad hoc assumptions on the galactic gas density distribution.

Cosmic-ray antiprotons represent an important channel for dark matter indirect-detection studies. Current measurements of the antiproton flux at the top of the atmosphere and theoretical determinations of the secondary antiproton production in the Galaxy are in good agreement, with no manifest deviation which could point to an exotic contribution in this channel. Therefore, antiprotons can be used as a powerful tool for constraining particle dark matter properties. By using the spectrum of PAMELA data from 50 MV to 180 GV in rigidity, we derive bounds on the dark matter annihilation cross section (or decay rate, for decaying dark matter) for the whole spectrum of dark matter annihilation (decay) channels and under different hypotheses of cosmic-rays transport in the Galaxy and in the heliosphere. For typical models of galactic propagation, the constraints are significantly strong, setting a lower bound on the dark matter mass of a "thermal" relic at about 50-90 GeV for hadronic annihilation channels. These bounds are enhanced to about 150 GeV on the dark matter mass, when large cosmic-rays confinement volumes in the Galaxy are considered, and are reduced to 4-5 GeV for annihilation to light quarks (no bound for heavy-quark production) when the confinement volume is small. Bounds for dark matter lighter than few tens of GeV are due to the low energy part of the PAMELA spectrum, an energy region where solar modulation is relevant: to this aim, we have implemented a detailed solution of the transport equation in the heliosphere, which allowed us not only to extend bounds to light dark matter, but also to determine the uncertainty on the constraints arising from solar modulation modeling. Finally, we estimate the impact of soon-to-come AMS-02 data on the antiproton constraints.

The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations.

The PAMELA collaboration recently released the $e^+$ absolute spectrum between 1 and 300 GeV in addition to the positron fraction and the $e^-$ spectrum previously measured in the same period. We use the newly developed three-dimensional upgrade of the DRAGON package to model those data. This code allows us to consider a realistic spiral arm source distribution in the Galaxy, which impacts the high-energy shape of the propagated spectra. At low energy we treat solar modulation with the HelioProp code and compare its results with those obtained using the usual force-field approximation. We show that all PAMELA data sets can be consistently, and accurately, described in terms of a standard background on top of which a charge symmetric $e^+$ + $e^-$ extra component with harder injection spectrum is added; this extra contribution is peaked at about 1 - 10 TeV and may originate from a diffuse population of sources located in the Galactic arms. For the first time, we compute the energy required to sustain such a relevant positron flux in the Galaxy, finding that it is naturally compatible with an astrophysical origin. We considered several reference propagation setups; we find that models with a low (or null) reacceleration - tuned against light nuclei data - nicely describe both PAMELA leptonic and hadronic data with no need to introduce a low-energy break in the proton and helium spectra, as it would be required for high reacceleration models. We also compare our models with the preliminary $e^-$ and $e^+$ absolute spectra recently measured by AMS-02. We find that those data, differently from what is inferred from the positron fraction alone, favor a high energy cutoff (about 10 TeV) of the extra component if this is uniquely generated in the Galactic arms. A lower cutoff may be allowed if a relevant contribution from powerful $e^-$ + $e^+$ nearby accelerators is invoked.

The emergence of a classical spacetime from any quantum gravity model is still a subtle and only partially understood issue. If indeed spacetime is arising as some sort of large scale condensate of more fundamental objects then it is natural to expect that matter, being a collective excitations of the spacetime constituents, will present modified kinematics at sufficiently high energies. We consider here the phenomenology of the dissipative effects necessarily arising in such a picture. Adopting dissipative hydrodynamics as a general framework for the description of the energy exchange between collective excitations and the spacetime fundamental degrees of freedom, we discuss how rates of energy loss for elementary particles can be derived from dispersion relations and used to provide strong constraints on the base of current astrophysical observations of high energy particles.

We model the CR leptonic fluxes above 20 GeV in terms of a superposition of a standard and a charge symmetric extra component, which we generically describe as power-laws in momentum. We investigate under these hypotheses the compatibility between AMS-02, Fermi-LAT and PAMELA datasets on positron fraction, electron+positron spectrum and electron spectrum respectively. We find that it is possible to reconcile AMS and Fermi-LAT data within uncertainties, if energy-dependent effects are present in Fermi-LAT systematics. We also study possible deviations from charge symmetry in the extra component and find no compelling evidence for them.

The search for an excess of antideuterons in the cosmic rays flux has been proposed as a very promising channel for dark matter indirect detection, especially for WIMPs with a low or intermediate mass. With the development of the AMS experiment and the proposal of a future dedicated experiment, i.e. the General Antiparticle Spectrometer (GAPS), there are exciting possibilities for a dark matter detection in the near future. We give an overview on the principal issues related both to the antideuterons production in dark matter annihilation reactions and to their propagation through the interstellar medium and the heliosphere, with a particular focus on the impact of various solar modulation models on the flux at Earth. Lastly, we provide an updated calculation of the reaching capabilities for current and future experiments compatible with the constraints on the dark matter annihilation cross section imposed by the antiproton measurements of PAMELA.

We study the compatibility of Alpha Magnetic Spectrometer (AMS-02) data on the cosmic-ray (CR) positron fraction with data on the CR electron and positron spectra provided by PAMELA and Fermi LAT. We do that in terms of a novel propagation model in which sources are distributed in spiral arm patterns in agreement with astrophysical observations. While former interpretations assumed an unrealistically steep injection spectrum for astrophysical background electrons, the enhanced energy losses experienced by CR leptons due to the larger average source distance from Earth allow us to reproduce the data with harder injection spectra as expected in a shock acceleration scenario. Moreover, we show that in this approach, and accounting for AMS-02 results, the contribution of nearby accelerators to the fluxes at very high energy can be significantly reduced, thus avoiding any tension with anisotropy upper limits.

The Fermi-LAT collaboration recently confirmed a discrepancy between the observed longitudinal profile of gamma-ray diffuse emission from the Galaxy and that computed with numerical codes assuming that Cosmic Rays (CRs) are produced by Galactic supernova remnants; the accurate Fermi-LAT measurements make this anomaly hardly explainable in terms of conventional diffusion schemes. Moreover, experimental data from both Muon detector and Extensive Air Shower experiments about the large scale dipole anisotropy of CRs can hardly be compatible with model predictions within the framework of conventional isotropic and homogeneous propagation. We argue that, accounting for a well physically motivated correlation between the CR escape time and the spatially dependent magnetic turbulence power, it is possible to solve both problems at the same time in a very natural way. Indeed, by exploiting this correlation, we find propagation models that fit a wide set of CR primary and secondary spectra, and consistently reproduce the CR anisotropy in the energy range 100 GeV - 10 TeV and the gamma-ray longitude distribution recently measured by Fermi-LAT.

We present here a major upgrade of DRAGON, a numerical package that computes the propagation of a wide set of CR species from both astrophysical and exotic origin in the Galaxy in a wide energy range from tens of MeV to tens of TeV. DRAGON takes into account all relevant processes in particular diffusion, convection, reacceleration, fragmentation and energy losses. For the first time, we present a full 3D version of DRAGON with anisotropic position-dependent diffusion. In this version, the propagation is calculated within a 3D cartesian grid and the user is able to implement realistic and structured three dimensional source, gas and regular magnetic field distributions. Moreover, it is possible to specify an arbitrary function of position and rigidity for the diffusion coefficients in the parallel and perpendicular direction to the regular magnetic field of the Galaxy. The code opens many new possibilities in the study of CR physics. In particular, we can study for the first time the impact of the spiral arm structure on the leptonic spectra: taking into account the fact that we live in an interarm region, far from most sources, we obtain - due to increased energy losses - a steeper electron spectrum compared to the assumption of a smooth source term. We discuss the implications of these results on our understanding on leptonic spectra and we briefly mention future studies that can be performed with our new 3D code.

We perform a consistent modeling of cosmic ray electrons, positrons and of the radio emission of the Galaxy. For the time we reproduce all relevant data sets between 1 GeV and 1 TeV including the recent AMS-02 positron fraction results. We show that below few GeV cosmic ray and radio data require that electron primary spectrum to be drastically suppressed and the propagated spectrum be dominated by secondary particles. Above 10 GeV an electron + positron extra-component with a hard spectrum is required. The positron spectrum measured below few GeV is consistently reproduced only within low reacceleration models. We also constrain the scale-height of the cosmic-ray distribution showing that a thin halo ($z_t \lsim 2 \kpc$) is excluded.

We take advantage of the cosmic-ray (CR) proton data collected by the PAMELA detector in several campaigns covering the period 2006-2009 to derive the local interstellar CR proton spectrum (LIS), outside the solar system. We describe the propagation of CR protons in the solar system by means of a detailed model of the charge-sign dependent drifts occurring while CRs diffuse in the irregular solar magnetic field. We fit PAMELA data for each year of data taking and find that they are well described by a unique time-independent LIS. We then discuss the consequences of this LIS on galactic propagation models. We find that diffusive reacceleration is strongly constrained, thus confirming previous results derived from different channels.

A multichannel analysis of cosmic ray electron and positron spectra and of the diffuse synchrotron emission of the Galaxy is performed by using the DRAGON code. This study is aimed at probing the interstellar electron source spectrum down to E ~ 1 GeV and at constraining several propagation parameters. We find that above 4 GeV the electron source spectrum is compatible with a power-law of index -2.5. Below 4 GeV instead it must be significantly suppressed and the total lepton spectrum is dominated by secondary particles. The positron spectrum and fraction measured below a few GeV are consistently reproduced only within low reacceleration models. We also constrain the scale-height zt of the cosmic-ray distribution using three independent (and, in two cases, original) arguments, showing that values of z_t < 2 kpc are excluded. This result may have strong implications for particle dark matter searches.

We compute cosmic ray (CR) nuclei, proton, antiproton, electron and positron spectra below 1 TeV at Earth by means of a detailed transport description in the galaxy and in the solar system. CR spectra below 10 GeV are strongly modified by charge-sign dependent propagation effects. These depend on the polarity of the solar magnetic field and therefore vary with the solar cycle. The puzzling discrepancy between the low-energy positron fraction measured by PAMELA and AMS-01 is then easily explained by their different data-taking epochs. We reproduce the observed spectra of CR light nuclei within the same galactic and solar-system propagation model.

Horava-Lifshitz gravity models contain higher order operators suppressed by a characteristic scale, which is required to be parametrically smaller than the Planck scale. We show that recomputed synchrotron radiation constraints from the Crab nebula suffice to exclude the possibility that this scale is of the same order of magnitude as the Lorentz breaking scale in the matter sector. This highlights the need for a mechanism that suppresses the percolation of Lorentz violation in the matter sector and is effective for higher order operators as well.

Very high energy (VHE, energy >~ 100 GeV) {\gamma}-rays originating from extragalactic sources interact with low energy photons of background radiation fields and produce electron-positron pairs. Alternatively, in the presence of ambient magnetic fields, they can convert into hypothetical axion-like particles (ALPs), pseudo-scalar spin-0 bosons, predicted by extensions of the standard model. These particles propagate unimpeded over cosmological distances. Here, the effect of photon-ALP oscillations in magnetic fields of galaxy clusters and the Milky Way on VHE {\gamma}-ray spectra is studied. It is shown that this mechanism can lead to a substantial enhancement of the VHE flux and a spectral hardening, thus effectively reducing the opacity of the Universe to VHE {\gamma}-rays.

A fraction of AGN producing VHE gamma-rays are located in galaxy clusters. The magnetic field present in the intra-cluster medium would lead to conversions of VHE photons into axion-like particles (ALPs), which are a generic prediction of several extensions of the Standard Model. ALPs produced in this way would traverse cosmological distances unaffected by the extragalactic background light at variance with VHE photons which undergo a substantial absorption. Eventually, a nontrivial fraction of ALPs would re-convert into VHE photons in the magnetic field of the Milky Way. This mechanism produces a significant hardening of the VHE spectrum of AGN in galaxy clusters. As a specific example we consider the energy spectra of two observed VHE gamma-ray sources located in galaxy clusters, namely 1ES 0414+009 at redshift z=0.287 and Mkn 501 at z=0.034. We find that the hardening in the observed spectra becomes relevant at E > 1 TeV. The detection of this signature would allow to indirectly probe the existence of ultra-light ALPs with mass m_a < 10^{-8} eV and photon-ALP coupling g_{a gamma} < 10^{-10} GeV^{-1} with the presently operating Imaging Atmospheric Cherenkov Telescopes like H.E.S.S., MAGIC, VERITAS and CANGAROO-III and even more likely with the planned detectors like CTA, HAWC and HiSCORE. An independent laboratory check of ultra-light ALPs invoked in this mechanism can be performed with the planned upgrade of the photon regeneration experiment ALPS at DESY and with the next generation solar axion detector IAXO.

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.

Version 2.0 of CRPropa is public software to model the extra-galactic propagation of ultra-high energy nuclei of atomic number Z<26 through structured magnetic fields and ambient photon backgrounds taking into account all relevant particle interactions. CRPropa covers the energy range 6*10^16 < E/eV < A*10^22 where A is the nuclear mass number. CRPropa can also be used to track secondary \gamma-rays and neutrinos which allows the study of their link with the charged primary nuclei -- the so called multi-messenger connection. After a general introduction we present several sample applications of current interest concerning the physics of extragalactic ultra-high energy radiation.

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.

Multichannel Cosmic Ray (CR) spectra and the large scale CR anisotropy can hardly be made compatible in the framework of conventional isotropic and homogeneous propagation models. These models also have problems explaining the longitude distribution and the radial emissivity gradient of the $\gamma$-ray galactic interstellar emission. We argue here that accounting for a well physically motivated correlation between the CR escape time and the spatially dependent magnetic turbulence power can naturally solve both problems. Indeed, by exploiting this correlation we find propagation models that fit a wide set of CR primary and secondary spectra, and consistently reproduce the CR anisotropy in the energy range $10^2 - 10^4 \GeV$ and the $\gamma$-ray longitude distribution recently measured by Fermi-LAT.

The mixing of photons with axion-like particles (ALPs) in the large-scale magnetic field $B$ changes the polarization angle of a linearly polarized photon beam from active galactic nuclei in radio galaxies as it propagates over cosmological distances. Using available ultraviolet polarization data concerning these sources we derive a new bound on the product of the photon-ALP coupling $g_{a\gamma}$ times $B$. We find $g_{a\gamma} B \lesssim 10^{-11}$ GeV$^{-1}$ nG for ultralight ALPs with $m_a \lesssim 10^{-15}$ eV. We compare our new bound with the ones present in the literature and we comment about possible improvements with observations of more sources.

We study the dipole anisotropy in the arrival directions of high energy CR electrons and positrons (CRE) of Dark Matter (DM) origin. We show that this quantity is very weakly model dependent and offers a viable criterion to discriminate among CRE from DM or from local discrete sources, like e.g. pulsars. In particular, we find that the maximum anisotropy which DM can provide is to a very good approximation a universal quantity and, as a consequence, if a larger anisotropy is detected, this would constitute a strong evidence for the presence of astrophysical local discrete CRE sources, whose anisotropy, instead, can be naturally larger than the DM upper limit. We further find that the main source of anisotropy from DM is given by the fluctuation in the number density of DM sub-structures in the vicinity of the observer and we thus devote special attention to the study of the variance in the sub-structures realization implementing a dedicated Montecarlo simulation. Such scenarios will be probed in the next years by Fermi-LAT, providing new hints, or constraints, about the nature of DM.

The Fermi-LAT collaboration recently confirmed EGRET finding of a discrepancy between the observed longitudinal profile of gamma-ray diffuse emission from the Galaxy and that computed with GALPROP assuming that cosmic rays are produced by Galactic supernova remnants. The accurate Fermi-LAT measurements make this anomaly hardly explainable in terms of conventional diffusion schemes. Here we use DRAGON numerical diffusion code to implement a physically motivated scenario in which the diffusion coefficient is spatially correlated to the source density. We show that under those conditions we are able to reproduce the observed flat emissivity profile in the outer Galaxy with no need to change the source term, the diffusion halo height, or the CO-H2 conversion factor (XCO) with respect to their preferred values/distributions. We also show that our models are compatible with gamma-ray longitudinal profiles measured by Fermi-LAT, and still provide a satisfactory fit of all observed secondary-to-primary ratios, such as B/C and antiprotons/protons.

Motivated by scenarios of quantum gravity, Planck-suppressed deviations from Lorentz invariance are expected at observable energies. Ultra-High-Energy Cosmic Rays, the most energetic particles ever observed in nature, yielded in the last two years strong constraints on deviations suppressed by $O(E^{2}/\Mpl^{2})$ and also, for the first time, on space-time foam, stringy inspired models of quantum gravity. We review the most important achievements and discuss future outlooks.

Motivated by experimental indications of a significant presence of heavy nuclei in the cosmic ray flux at ultra high energies ($\gtrsim 10^{19} \eV$), we consider the effects of Planck scale suppressed Lorentz Invariance Violation (LIV) on the propagation of cosmic ray nuclei. In particular we focus on LIV effects on the photodisintegration of nuclei onto the background radiation fields. After a general discussion of the behavior of the relevant quantities, we apply our formalism to a simplified model where the LIV parameters of the various nuclei are assumed to kinematically result from a single LIV parameter for the constituent nucleons, $\eta$, and we derive constraints on $\eta$. Assuming a nucleus of a particular species to be actually present at $10^{20}$ eV the following constraints can be placed: $-3\times10^{-2} \lesssim \eta \lesssim 4$ for $^{56}$Fe, $-2\times10^{-3} \lesssim \eta \lesssim 3\times10^{-2}$ for $^{16}$O and $-7\times10^{-5} \lesssim \eta \lesssim 1\times10^{-4}$ for $^{4}$He, respectively.

We try to interpret the recently updated measurement of the cosmic ray electron (CRE) spectrum observed by Fermi-LAT, together with PAMELA data on positron fraction, in a single-component scenario adopting different propagation setups; we find that the model is not adequate to reproduce the two datasets, so the evidence of an extra primary component of electrons and positrons is strengthened. Instead, a double component scenario computed in a Kraichnan-like diffusion setup (which is suggested by B/C and $\bar{p}$ data) gives a satisfactory fit of all exisiting measurements. We confirm that nearby pulsars are good source candidates for the required $e^\pm$ extra-component and we show that the predicted CRE anisotropy in our scenario is compatible with Fermi-LAT recently published constraints.

The Fermi Large Area Telescope (LAT) collaboration recently released the updated results of the measurement of the cosmic ray electron (CRE) spectrum and published its first constraints on the CRE anisotropy. With respect to the previous Fermi-LAT results, the CRE spectrum measurement was extended down from 20 to 7 GeV, thus providing a better lever arm to discriminate theoretical models. Here we show that the new data strengthen the evidence for the presence of two distinct electron and positron spectral components. Furthermore, we show that under such hypothesis most relevant CRE and positron data sets are remarkably well reproduced. Consistent fits of cosmic-ray nuclei and antiproton data, which are crucial to validate the adopted propagation setup(s) and to fix the solar modulation potential, are obtained for the Kraichnan and plain-diffusion propagation setups, while the Kolmogorov one is disfavored. We then confirm that nearby pulsars are viable source candidates of the required $e^\pm$ extra-component. In that case, we show that the predicted CRE anisotropy is compatible with Fermi-LAT constraints and that a positive detection should be at hand of that observatory. Models assuming that only nearby supernova remnants contribute to the high energy tail of the observed CRE spectrum are in contrast with anisotropy limits.

We use our numerical code, DRAGON, to study the implications of recent data on our knowledge of the propagation properties of cosmic ray nuclei in the Galaxy. We show that B/C (as well as N/O and C/O) data, including those recently taken by CREAM, and ${\bar p}/p$ data, especially including recent PAMELA results, can consistently be fitted within a unique diffusion-reacceleration model. The requirement that light nuclei and $\bar p$ data are consistently reproduced within experimental uncertainties places significant limits on the main propagation parameters. In particular, we find the allowed range of the diffusion coefficient spectral index to be $0.3 < \delta < 0.6$ at 95% confidence level and that Kraichnan type diffusion is favored with respect to Kolmogorov. While some amount of reacceleration is required, only a limited range of the Alfv\'en velocity value ($10 \simleq v_A \simleq 20 \km \s^{-1}$) is allowed by a combined analysis of nuclear and antiproton data, which we perform here for the first time. If antiproton data are not used to constrain the propagation parameters, a larger set of models is allowed. In this case, we determine which combinations of the relevant parameters maximize and minimize the antiproton flux under the condition of still fitting light nuclei data at 95% C.L. These models may then be used to constrain a possible extra antiproton component arising from astrophysical or exotic sources (e.g. dark matter annihilation or decay).

The time delays between gamma-rays of different energies from extragalactic sources have often been used to probe quantum gravity models in which Lorentz symmetry is violated. It has been claimed that these time delays can be explained by or at least put the strongest available constraints on quantum gravity scenarios that cannot be cast within an effective field theory framework, such as the space-time foam, D-brane model. Here we show that this model would predict too many photons in the ultra-high energy cosmic ray flux to be consistent with observations. The resulting constraints on the space-time foam model are much stronger than limits from time delays and allow for Lorentz violations effects way too small for explaining the observed time delays.

High energy electrons and positrons from decaying dark matter can produce a significant flux of gamma rays by inverse Compton off low energy photons in the interstellar radiation field. This possibility is inevitably related with the dark matter interpretation of the observed PAMELA and FERMI excesses. The aim of this paper is providing a simple and universal method to constraint dark matter models which produce electrons and positrons in their decay by using the Fermi LAT gamma-ray observations in the energy range between 0.5 GeV and 300 GeV. We provide a set of universal response functions that, once convolved with a specific dark matter model produce the desired constraint. Our response functions contain all the astrophysical inputs such as the electron propagation in the galaxy, the dark matter profile, the gamma-ray fluxes of known origin, and the Fermi LAT data. We study the uncertainties in the determination of the response functions and apply them to place constraints on some specific dark matter decay models that can well fit the positron and electron fluxes observed by PAMELA and Fermi LAT. To this end we also take into account prompt radiation from the dark matter decay. We find that with the available data decaying dark matter cannot be excluded as source of the PAMELA positron excess.

We study, within an effective field theory framework, $O(E^{2}/\Mpl^{2})$ Planck-scale suppressed Lorentz invariance violation (LV) effects in the neutrino sector, whose size we parameterize by a dimensionless parameter $\eta_{\nu}$. We find deviations from predictions of Lorentz invariant physics in the cosmogenic neutrino spectrum. For positive O(1) coefficients no neutrino will survive above $10^{19} \eV$. The existence of this cutoff generates a bump in the neutrino spectrum at energies of $10^{17} \eV$. Although at present no constraint can be cast, as current experiments do not have enough sensitivity to detect ultra-high-energy neutrinos, we show that experiments in construction or being planned have the potential to cast limits as strong as $\eta_{\nu} \lesssim 10^{-4}$ on the neutrino LV parameter, depending on how LV is distributed among neutrino mass states. Constraints on $\eta_{\nu} < 0$ can in principle be obtained with this strategy, but they require a more detailed modeling of how LV affects the neutrino sector.

We investigate the consequences of higher dimension Lorentz violating, CPT even kinetic operators that couple standard model fields to a non-zero vector field in an Effective Field Theory framework. Comparing the ultra-high energy cosmic ray spectrum reconstructed in the presence of such terms with data from the Pierre Auger observatory allows us to establish two sided bounds on the coefficients of the mass dimension five and six operators for the proton and pion. Our bounds imply that for both protons and pions, the energy scale of Lorentz symmetry breaking must be well above the Planck scale. In particular, the dimension five operators are constrained at the level of $10^{-3} M_{\rm Planck}^{-1}$. The magnitude of the dimension six proton coefficient is bounded at the level of $10^{-6} M_{Planck}^{-2}$ except in a narrow range where the pion and proton coefficients are both negative and nearly equal. In this small area, the magnitude of the dimension six proton coefficient must only be below $10^{-3} M_{\rm Planck}^{-2}$. Constraints on the dimension six pion coefficient are found to be much weaker, but still below $M_{\rm Planck}^{-2}$.

In the last months several ballon and satellite experiments improved significantly our knowledge of cosmic rays (CR) spectra at high energy. In particular CREAM allowed to measure B/C, C/O and N/O up to 1 TeV and PAMELA the anti-p/p ratio up to 100 GeV with unprecedented accuracy. These measurements offer a valuable probe of CR propagation properties. We performed a statistical analysis to test the compatibility of these results, as well as other most significant experimental data, with the predictions of a new numerical CR diffusion package (DRAGON). We found that above 1 GeV all data are consistent with a diffusion scenario in a well defined range of values of the diffusion coefficient energy power index and normalization.

We review the main theoretical motivations and observational constraints on Planck scale suppressed violations of Lorentz invariance. After introducing the problems related to the phenomenological study of quantum gravitational effects, we discuss the main theoretical frameworks within which possible departures from Lorentz invariance can be described. In particular, we focus on the framework of Effective Field Theory, describing several possible ways of including Lorentz violation therein and discussing their theoretical viability. We review the main low energy effects that are expected in this framework. We discuss the current observational constraints on such a framework, focusing on those achievable through high-energy astrophysics observations. In this context we present a summary of the most recent and strongest constraints on QED with Lorentz violating non-renormalizable operators. Finally, we discuss the present status of the field and its future perspectives.

Using recent polarimetric observations of the Crab Nebula in the hard X-ray band by INTEGRAL, we show that the absence of vacuum birefringence effects constrains O(E/M) Lorentz violation in QED to the level |\xi| < 9x10^{-10} at three sigma CL, tightening by more than three orders of magnitude previous constraints. We show that planned X-ray polarimeters have the potential the potential to probe |\xi|~ 10^{-16} by detecting polarization in active galaxies at red-shift ~1.

We show how the argument exploited by Galaverni & Sigl in Phys. Rev. Lett., 100, 021102 (2008) (see also arXiv:0708.1737) to constrain Lorentz invariance violation (LV) using Ultra-High-Energy photon non observation by the AUGER experiment, can be extended to QED with Planck-suppressed LV (at order $O(E/M)$ and $O(E^2/M^2)$). While the original constraints given by Galaverni & Sigl happen to be weakened, we show that, when used together with other EFT reactions and the expected detection of photons at $E > 10^{19}$ eV by AUGER, this method has the potentiality not only to basically rule out order $O(E/M)$ corrections but also to strongly constrain, for the first time, the CPT-even $O(E^2/M^2)$ LV QED.

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

We set constraints on O(E/M) Lorentz Violation in QED in an effective field theory framework. A major consequence of such assumptions is the modification of the dispersion relations for electrons/positrons and photons, which in turn can affect the electromagnetic output of astrophysical objects. We compare the information provided by multiwavelength observations with a full and self-consistent computation of the broad-band spectrum of the Crab Nebula. We cast constraints of order 10^{-5} at 95% confidence level on the lepton Lorentz Violation parameters.

In this contribution we will discuss recent results concerning the intensity and the angular distribution of the gamma-ray and neutrino emissions as should be originated from the hadronic scattering of cosmic rays (CR) with the interstellar medium (ISM). We assumed that CR sources are supernova remnants (SNR) and estimated the spatial distribution of primary nuclei by solving numerically the diffusion equation. For the ISM, we considered recent models for the 3D spatial distributions of molecular hydrogen. Respect to previous results, we find the secondary gamma-ray and neutrino emissions to be more peaked along the galactic equator and in the galactic centre which improves significantly the perspectives of a positive detection. We compare our predictions with the experimental limits/observations by MILAGRO and TIBET (for the gamma-rays) and by AMANDA-II (for the neutrinos) and discuss the detection perspectives for a km3 neutrino telescope to be built in the North hemisphere.

We simulate the neutrino and $\gamma$-ray emissions of the Galaxy which are originated from the hadronic scattering of cosmic rays (CR) with the interstellar medium (ISM). Rather than assuming a uniform CR density, we estimate the spatial distribution of CR nuclei by means of numerical simulations. We consider several models of the galactic magnetic field and of the ISM distribution finding only a weak dependence of our results on their choice. We found that by extrapolating the predicted $\gamma$-ray spectra down to few GeV we get a good agreement with EGRET measurements. Then, we can reliably compare our predictions with available observations above the TeV both for the $\gamma$-rays and the neutrinos. We confirm that the excesses observed by MILAGRO in the Cygnus region and by HESS in the Galactic Centre Ridge cannot be explained without invoking significant CR over-densities in those regions. Finally, we discuss the perspectives that a km$^3$ neutrino telescope based in the North hemisphere has to measure the diffuse emission from the inner Galaxy.

The recent detection of $\gamma$-rays with energy up to 10 TeV from dense regions surrounding some Supernova Remnants (SNR) provides strong, though still not conclusive, evidence that the nucleonic component of galactic Cosmic Rays is accelerated in the supernova outflows. Neutrino telescopes could further support the validity of such scenario by detecting neutrinos coming from the same regions. We re-evaluate the TeV range neutrino-photon flux ratio to be expected from pion decay, finding small differences respect to previous derivations. We apply our results and the recent HESS measurements of the very high energy $\gamma$-ray flux from the molecular cloud complex in the Galactic Centre, to estimate the neutrino flux from that region discussing the detectability perspectives for Mediterranean Neutrino Telescopes. We also discuss under which conditions neutron decay may give rise to a significant TeV antineutrino flux from a SNR embedded in Molecular Cloud complexes.

Sgr A East is a supernova remnant located within few parsecs from the Galactic Centre (GC). There are good reasons to believe that this object is the source of the gamma-ray excess detected by HESS in the direction of the GC meaning that Sgr A East is likely to be an efficient Cosmic Ray accelerator. Some observations suggest that strong magnetic fields may be present in that region allowing the acceleration of composite nuclei in Sgr A East beyond the EeV. We show that, if this is case, EeV neutrons should be effectively produced by the photo-disintegration of Ultra High Energy nuclei onto the IR photon background (with temperature $\sim 40$ K) in which Sgr A East is embedded. Neutrons with such an energy can reach the Earth before decaying and may be detectable under the form of a CR point-like excess in the direction of the GC. We determine the expected energy spectrum and the amplitude of this signal showing that it may be measurable by the AUGER observatory.