We present a comprehensive discussion about the origin of the features in the leptonic component of the cosmic-ray spectrum. Working in the framework of a up-to-date CR transport scenario tuned on the most recent AMS-02 and Voyager data, we show that the prominent features recently found in the positron and in the all-electron spectra by several experiments are explained in a scenario in which pulsar wind nebulae (PWNe) are the dominant sources of the positron flux, and nearby supernova remnants shape the high-energy peak of the electron spectrum. In particular we argue that the drop-off in positron spectrum found by AMS-02 at $\sim 300$ GeV can be explained --- under different assumptions --- in terms of a prominent PWN that provides the bulk of the observed positrons in the $\sim 100$ GeV domain, on top of the contribution from a large number of older objects. Finally, we turn our attention to the spectral softening at $\sim 1$ TeV in the all-lepton spectrum, recently reported by several experiments, showing that it requires the presence of a nearby supernova remnant at its final stage.

A dark matter overdensity around a black hole may significantly alter the dynamics of the black hole's merger with another compact object. We consider here intermediate mass-ratio inspirals of stellar-mass compact objects with intermediate-mass black holes "dressed" with dark matter. We first demonstrate that previous estimates based on a fixed dark-matter dress are unphysical for a range of binaries and dark-matter distributions by showing that the total energy dissipated by the compact object through dynamical friction, as it inspirals through the dense dark matter environment towards the black hole, is larger than the gravitational binding energy of the dark-matter dress itself. We then introduce a new formalism that allows us to self-consistently follow the evolution of the dark-matter dress due to its gravitational interaction with the binary. We show that the dephasing of the gravitational waveform induced by dark matter is smaller than previously thought, but is still potentially detectable with the LISA space interferometer. The gravitational waves from such binaries could provide powerful diagnostics of the particle nature of dark matter.

As smaller analogs of Active Galactic Nuclei (AGN), X-ray Binaries (XRBs) are also capable of launching jets that accelerate particles to high energies. In this work, we reexamine XRB jets as potential sources of high-energy cosmic rays (CRs) and explore whether they could provide a significant second Galactic component to the CR spectrum. In the most intriguing scenario, XRB-CRs could dominate the observed spectrum above the so-called "knee" feature at $\sim 3 \times 10^{15}$ eV, offering an explanation for several key issues in this transition zone from Galactic to extragalactic CRs. We discuss how such a scenario could be probed in the near future via multi-messenger observations of XRB jets, as well as diffuse Galactic neutrino flux measurements.

Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming Einstein Telescope) would set extraordinarily stringent constraints on virtually all weak-scale extensions of the Standard Model with stable relics, including those predicting a WIMP abundance much smaller than that of Dark Matter. Upcoming PBHs searches have in particular the potential to rule out almost the entire parameter space of popular theories such as the minimal supersymmetric standard model and scalar singlet Dark Matter.

A critical review of the standard paradigm for the origin of Galactic cosmic rays is presented. Recent measurements of local and far-away cosmic rays reveal unexpected behaviours, which challenge the commonly accepted scenario. These recent findings are discussed, together with long-standing open issues. Despite the progress made thanks to ever-improving observational techniques and theoretical investigations, at present our understanding of the origin and of the behaviour of cosmic rays remains incomplete. We believe it is still unclear whether a modification of the standard paradigm, or rather a radical change of the paradigm itself is needed in order to interpret all the available data on cosmic rays within a self-consistent scenario.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

Among all the high-energy environments of our Galaxy, the Galactic Center (GC) region is definitely the richest. It harbors a large amount of non-thermal emitters, including the closest supermassive black hole, dense molecular clouds, regions with strong star forming activity, multiple supernova remnants and pulsar wind nebulae, arc-like radio structures, as well as the base of what may be large-scale Galactic outflows, possibly related to the Fermi Bubbles. It also contains a strong diffuse TeV gamma-ray emission along the Galactic ridge, with a disputed origin, including the presence of a possible Pevatron, unresolved sources, and an increased relevance of the diffuse sea of cosmic rays. This very rich region will be one of the key targets for the next generation ground-based observatory for gamma-ray astronomy, the Cherenkov Telescope Array (CTA). Here we review the CTA science case for the study of the GC region, and present the planned survey strategy. These observations are simulated and we assess CTA's potential to better characterize the origin and nature of a selection of gamma-ray sources in the region.

If primordial black holes of $\mathcal{O}(1-100) M_{\odot}$ constitute a significant portion of the dark matter in the Universe, they should be very abundant in our Galaxy. We present here a detailed analysis of the radio and X-ray emission that these objects are expected to produce due to the accretion of gas from the interstellar medium. With respect to previous studies, we relax the assumption of a monochromatic mass function, and introduce an improved treatment of the physics of gas accretion onto isolated, moving compact objects, based on a set of state-of-the-art numerical simulations. By comparing our predictions with known radio and X-ray sources in the Galactic center region, we show that the maximum relic density of primordial black holes in the mass range of interest is $\sim 10^{-3}$ smaller than that of dark matter. The new upper bound is two orders of magnitude stronger with respect to previous results, based on a conservative phenomenological treatment of the accretion physics. We also provide a comprehensive critical discussion on the reliability of this bound, and on possible future developments in the field. We argue in particular that future multi-wavelength searches will soon start to probe the galactic population of astrophysical black holes.

The non-observation of Hawking radiation from primordial black holes of $10^{16}$g sets a conservative strong bound on their cosmological abundance. We revisit this bound and show how it can be improved (both in mass reach and strength) by an adequate modeling of the combined AGN and blazar emission in the MeV range. We also estimate the sensitivity to the primordial black hole abundance of a future X-ray experiment capable of identifying a significantly larger number of astrophysical sources contributing to the diffuse background in this energy range.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The formation of astrophysical and primordial black holes influences the distribution of dark matter surrounding them. Black holes are thus expected to carry a dark matter `dress' whose properties depend on their formation mechanism and on the properties of the environment. Here we carry out a numerical and analytical study of the merger of dressed black holes, and show that the distribution of dark matter around them dramatically affects the dynamical evolution of the binaries. Although the final impact on the merger rate of primordial black holes is rather small with respect to the case of `naked' black holes, we argue that our analysis places the calculation of this rate on more solid ground, with LIGO-Virgo observations ruling out a dark matter fraction of $10^{-3}$ for primordial black holes of 100 solar masses, and it paves the way to more detailed analyses of environmental effects induced by dark matter on the gravitational wave emission of binary black holes.

The existence of diffuse Galactic neutrino production is expected from cosmic ray interactions with Galactic gas and radiation fields. Thus, neutrinos are a unique messenger offering the opportunity to test the products of Galactic cosmic ray interactions up to energies of hundreds of TeV. Here we present a search for this production using ten years of ANTARES track and shower data, as well as seven years of IceCube track data. The data are combined into a joint likelihood test for neutrino emission according to the KRA$_\gamma$ model assuming a 5 PeV per nucleon Galactic cosmic ray cutoff. No significant excess is found. As a consequence, the limits presented in this work start constraining the model parameter space for Galactic cosmic ray production and transport.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Spatial variations of the average properties that characterize the hadronic component of the diffuse Galactic cosmic-ray sea, in particular the spectral slope and normalization, may unveil critical information about their confinement mechanism in the Galaxy. In the first part of this paper we perform an analysis of the Fermi-LAT gamma-ray data with the SkyFACT package, which combines image reconstruction techniques with standard template fitting, isolate the hadronic emission and decompose it into Galactocentric rings. We find a significant hardening of the hadronic spectral index towards the molecular ring. We study this hardening in different energy ranges, and assess its resilience with respect to different prescriptions in the analysis setup. In the second part we quantify the contribution to the diffuse gamma-ray flux coming from unresolved point sources with a dedicated Monte Carlo simulation, and consider whether the trend characterized in the first part can be mimicked by a progressively more relevant flux associated to this component in the inner Galaxy. We find that the observed hardening of the hadronic spectral index cannot be due to unresolved sources in the sub-TeV energy range, especially outside the molecular ring, given reasonable assumptions about the unresolved source population.

Understanding the isotopic composition of cosmic rays (CRs) observed near Earth represents a milestone towards the identification of their origin. Local fluxes contain all the known stable and long-lived isotopes, reflecting the complex history of primaries and secondaries as they traverse the interstellar medium. For that reason, a numerical code which aims at describing the CR transport in the Galaxy must unavoidably rely on accurate modelling of the production of secondary particles. In this work we provide a detailed description of the nuclear cross sections and decay network as implemented in the forthcoming release of the galactic propagation code DRAGON2. We present the secondary production models implemented in the code and we apply the different prescriptions to compute quantities of interest to interpret local CR fluxes (e.g., nuclear fragmentation timescales, secondary and tertiary source terms). In particular, we develop a nuclear secondary production model aimed at accurately computing the light secondary fluxes (namely: Li, Be, B) above 1 GeV/n. This result is achieved by fitting existing empirical or semi-empirical formalisms to a large sample of measurements in the energy range 100 MeV/n to 100 GeV/n and by considering the contribution of the most relevant decaying isotopes up to iron. Concerning secondary antiparticles (positrons and antiprotons), we describe a collection of models taken from the literature, and provide a detailed quantitative comparison.

The origins of the extraterrestrial neutrinos observed in IceCube have yet to be determined. In this study we perform a one-point fluctuation analysis of the six-year high-energy starting event (HESE) shower data, with fixed non-Poissonian contributions from atmospheric, Galactic and some extragalactic components, as well as an isotropic (and weakly non-Poissonian) template. In addition to the star-forming galaxies and blazars, our analysis suggests the presence of an additional isotropic component, not associated with any known class of sources, with best-fit intensity of $(2.8\pm0.2)\times 10^{-18}\,(E/100~{\rm TeV})^{-2.7\pm 0.5}$ ${\mathrm{cm}^{-2}\, \mathrm{s}^{-1}\,\mathrm{sr}^{-1}\,\mathrm{GeV}^{-1}}$. For the first time, we also consider high-energy extrapolations of several phenomenological models for the diffuse Galactic emission (tuned to both local cosmic-ray data and diffuse gamma-ray emission in the GeV-TeV domain). We demonstrate the potential of our framework in discriminating between different scenarios, with possible implications on the physics of cosmic ray transport in the TeV-PeV range.

The quest for the elusive dark matter (DM) that permeates the Universe (and in general the search for signatures of Physics beyond the Standard Model at astronomical scales) provides a unique opportunity and a tough challenge to the high-energy astrophysics community. In particular, the so-called DM indirect searches - mostly focused on a class of theoretically well-motivated DM candidates such as the weakly-interacting massive particles - are affected by a complex astrophysical background of cosmic radiation. The understanding and modeling of such background requires a deep comprehension of an intricate classical plasma physics problem, i.e. the interaction between high-energy charged particles, accelerated in peculiar astrophysical environments, and magneto-hydrodynamic turbulence in the interstellar medium of our Galaxy. In this review we highlight several aspects of this exciting interplay between the most recent claims of DM annihilation/decay signatures from the sky and the Galactic cosmic-ray research field. Our purpose is to further stimulate the debate about viable astrophysical explanations, discussing possible directions that would help breaking degeneracy patterns in the interpretation of current data. We eventually aim to emphasize how a deep knowledge on the physics of CR transport is therefore required to tackle the DM indirect search program at present and in the forthcoming years.

A crucial process in Galactic cosmic-ray (CR) transport is the spatial diffusion due to the interaction with the interstellar turbulent magnetic field. Usually, CR diffusion is assumed to be uniform and isotropic all across the Galaxy. However, this picture is clearly inaccurate: Several data-driven and theoretical arguments, as well as dedicated numerical simulations, show that diffusion exhibits highly anisotropic properties with respect to the direction of a background (ordered) magnetic field (i.e., parallel or perpendicular to it). In this paper we focus on a recently discovered anomaly in the hadronic CR spectrum inferred by the Fermi-LAT gamma-ray data at different positions in the Galaxy, i.e. the progressive hardening of the proton slope at low Galactocentric radii. We propose the idea that this feature can be interpreted as a signature of anisotropic diffusion in the complex Galactic magnetic field: In particular, the harder slope in the inner Galaxy is due, in our scenario, to the parallel diffusive escape along the poloidal component of the large-scale, regular, magnetic field. We implement this idea in a numerical framework, based on the DRAGON code, and perform detailed numerical tests on the accuracy of our setup. We discuss how the effect proposed depends on the relevant free parameters involved. Based on low-energy extrapolation of the few focused numerical simulations aimed at determining the scalings of the anisotropic diffusion coefficients, we finally present a set of plausible models that reproduce the behavior of the CR proton slopes inferred by gamma-ray data.

We present a systematic phenomenological study focused on the time evolution of the non-thermal radiation - from radio waves to gamma rays - emitted by typical supernova remnants via hadronic and leptonic mechanisms, for two classes of progenitors: thermonuclear and core-collapse. To this aim, we develop a numerical tool designed to model the evolution of the cosmic-ray spectrum inside a supernova remnant, and compute the associated multi-wavelength emission. We demonstrate the potential of this tool in the context of future population studies based on large collection of high-energy gamma-ray data. We discuss and explore the relevant parameter space involved in the problem, and focus in particular on their impact on the maximum energy of accelerated particles, in order to study the effectiveness and duration of the PeVatron phase. We outline the crucial role of the ambient medium through which the shock propagates during the remnant evolution. In particular, we point out the role of dense clumps in creating a significant hardening in the hadronic gamma-ray spectrum.

The flux of very high-energy neutrinos produced in our Galaxy by the interaction of accelerated cosmic rays with the interstellar medium is not yet determined. The characterization of this flux will shed light on Galactic accelerator features, gas distribution morphology and Galactic cosmic ray transport. The central Galactic plane can be the site of an enhanced neutrino production, thus leading to anisotropies in the extraterrestrial neutrino signal as measured by the IceCube Collaboration. The ANTARES neutrino telescope, located in the Mediterranean Sea, offers a favourable view on this part of the sky, thereby allowing for a contribution to the determination of this flux. The expected diffuse Galactic neutrino emission can be obtained linking a model of generation and propagation of cosmic rays with the morphology of the gas distribution in the Milky Way. In this paper, the so-called "Gamma model" introduced recently to explain the high-energy gamma ray diffuse Galactic emission, is assumed as reference. The neutrino flux predicted by the "Gamma model" depends of the assumed primary cosmic ray spectrum cut-off. Considering a radially-dependent diffusion coefficient, this proposed scenario is able to account for the local cosmic ray measurements, as well as for the Galactic gamma ray observations. Nine years of ANTARES data are used in this work to search for a possible Galactic contribution according to this scenario. All flavour neutrino interactions are considered. No excess of events is observed and an upper limit is set on the neutrino flux of $1.1$ ($1.2$) times the prediction of the "Gamma model" assuming the primary cosmic ray spectrum cut-off at 5 (50) PeV. This limit excludes the diffuse Galactic neutrino emission as the major cause of the "spectral anomaly" between the two hemispheres measured by IceCube.

Before being detected at Earth, charged cosmic rays propagate across the Solar System and undergo interactions with the turbulent solar wind and with the heliospheric magnetic field. As a result, they are subject to a series of processes that include diffusion, convection, energy losses and drifts, which significantly affect the shape and the intensity of the cosmic-ray fluxes at low energies. Here we illustrate how all these mechanisms can be realistically modelled with HelioProp, our public tool designed to treat cosmic-ray transport through the heliosphere in a charge-dependent way. We present a detailed description of the features of the code and we illustrate in a quantitative way the effects that the propagation in the heliosphere can have on the different cosmic-ray species with a particular emphasis on the antiparticle channels relevant for dark matter indirect detection.

We present a novel interpretation of the gamma-ray diffuse emission measured by H.E.S.S. in the Galactic Center (GC) region and the Galactic ridge. Our starting base is an updated analysis of PASS8 Fermi-LAT data, which allows to extend down to few GeV the spectra measured by H.E.S.S. and to infer the primary CR radial distribution above 100 GeV. We compare those results with a CR transport model assuming a harder scaling of the diffusion coefficient with rigidity in the inner Galaxy. Such a behavior reproduces the radial dependence of the CR spectral index recently inferred from Fermi-LAT measurements in the inner GP. We find that, in this scenario, the bulk of the Galactic ridge emission can be naturally explained by the interaction of the diffuse, steady-state Galactic CR sea interacting with the gas present in the Central molecular zone. The evidence of a GC PeVatron is significantly weaker than that inferred adopting a conventional (softer) CR sea.

We model the accretion of gas onto a population of massive primordial black holes in the Milky Way, and compare the predicted radio and X-ray emission with observational data. We show that under conservative assumptions on the accretion process, the possibility that ${\cal O}(10) \, M_\odot$ primordial black holes can account for all of the dark matter in the Milky Way is excluded at $5\sigma$ by a comparison with a VLA radio catalog at $1.4$ GHz, and at $\simeq 40\sigma$ by a comparison with a Chandra X-ray catalog ($0.5 - 8$ keV). We argue that this method can be used to identify such a population of primordial black holes with more sensitive future radio and X-ray surveys.

We present a novel interpretation of the $\gamma$-ray diffuse emission measured by Fermi-LAT and H.E.S.S. in the Galactic center (GC) region and the Galactic ridge (GR). In the first part we perform a data-driven analysis based on PASS8 Fermi-LAT data: we extend down to few GeV the spectra measured by H.E.S.S. and infer the primary cosmic-ray (CR) radial distribution between 0.1 and 3 TeV. In the second part we adopt a CR transport model based on a position-dependent diffusion coefficient. Such behavior reproduces the radial dependence of the CR spectral index recently inferred from the Fermi-LAT observations. We find that the bulk of the GR emission can be naturally explained by the interaction of the diffuse steady-state Galactic CR sea with the gas present in the Central Molecular Zone. Although our results leave room for a residual radial-dependent emission associated with a central source, the relevance of the large-scale background prevents from a solid evidence of a GC Pevatron.

Over the past decade, extensive studies have been undertaken to search for photon signals from dark matter annihilation or decay for dark matter particle masses above $\sim1$ GeV. However, due to the lacking sensitivity of current experiments at MeV-GeV energies, sometimes dubbed the 'MeV gap', dark matter models with MeV to sub-GeV particle masses have received little attention so far. Various proposed MeV missions (like, e.g., e-ASTROGAM or AMEGO) are aimed at closing this gap in the mid- or long-term future. This, and the absence of clear dark matter signals in the GeV-TeV range, makes it relevant to carefully reconsider the expected experimental instrumental sensitivities in this mass range. The most common two-body annihilation channels for sub-GeV dark matter are to neutrinos, electrons, pions or directly to photons. Among these, only the electron channel has been extensively studied, and almost exclusively in the context of the 511 keV line. In this work, we study the prospects for detecting MeV dark matter annihilation in general in future MeV missions, using e-ASTROGAM as reference, and focusing on dark matter masses in the range 1 MeV-3 GeV. In the case of leptonic annihilation, we emphasise the importance of the often overlooked bremsstrahlung and in-flight annihilation spectral features, which in many cases provide the dominant gamma-ray signal in this regime.

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

The measured fluxes of secondary particles produced by the interactions of Cosmic Rays (CRs) with the astronomical environment play a crucial role in understanding the physics of CR transport. In this work we present a comprehensive calculation of the secondary hadron, lepton, gamma-ray and neutrino yields produced by the inelastic interactions between several species of stable or long-lived cosmic rays projectiles (p, D, T, 3He, 4He, 6Li, 7Li, 9Be, 10Be, 10B, 11B, 12C, 13C, 14C, 14N, 15N, 16O, 17O, 18O, 20Ne, 24Mg and 28Si) and different target gas nuclei (p, 4He, 12C, 14N, 16O, 20Ne, 24Mg, 28Si and 40Ar). The yields are calculated using FLUKA, a simulation package designed to compute the energy distributions of secondary products with large accuracy in a wide energy range. The present results provide, for the first time, a complete and self-consistent set of all the relevant inclusive cross sections regarding the whole spectrum of secondary products in nuclear collisions. We cover, for the projectiles, a kinetic energy range extending from $0.1 GeV/n$ up to $100 TeV/n$ in the lab frame. In order to show the importance of our results for multi-messenger studies about the physics of CR propagation, we evaluate the propagated spectra of Galactic secondary nuclei, leptons, and gamma rays produced by the interactions of CRs with the insterstellar gas, exploiting the numerical codes DRAGON and GammaSky. We show that, adopting our cross section database, we are able to provide a good fit of a complete sample of CR observables, including: leptonic and hadronic spectra measured at Earth, the local interstellar spectra measured by Voyager, and the gamma-ray emissivities from Fermi-LAT collaboration. We also show a set of gamma-ray and neutrino full-sky maps and spectra.

The measured fluxes of secondary particles produced by the interactions of cosmic rays with the astronomical environment represent a powerful tool to infer some properties of primary cosmic rays. In this work we investigate the production of secondary particles in inelastic hadronic interactions between several cosmic rays species of projectiles and different target nuclei of the interstellar medium. The yields of secondary particles have been calculated with the FLUKA simulation package, that provides with very good accuracy the energy distributions of secondary products in a large energy range. An application to the propagation and production of secondaries in the Galaxy is presented.

Supernova remnants are believed to be the sources of galactic cosmic rays. Within this framework, diffusive shock acceleration must operate in these objects and accelerate protons all the way up to PeV energies. To do so, significant amplification of the magnetic field at the shock is required. The goal of this paper is to investigate the capability of supernova remnants to accelerate PeV protons. We present analytic estimates of the maximum energy of accelerated protons under various assumptions about the field amplification at supernova remnant shocks. We show that acceleration up to PeV energies is problematic in all the scenarios considered. This implies that either a different (more efficient) mechanism of field amplification operates at supernova remnant shocks, or that the sources of galactic cosmic rays in the PeV energy range should be searched somewhere else.

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.

The last IceCube catalog of High Energy Starting Events (HESE) obtained with a livetime of 1347 days comprises 54 neutrino events equally-distributed between the three families with energies between 25 TeV and few PeVs. Considering the homogeneous flavors distribution (1:1:1) and the spectral features of these neutrinos the IceCube collaboration claims the astrophysical origin of these events with more than $5\sigma$. The spatial distribution of cited events does not show a clear correlation with known astrophysical accelerators leaving opened both the Galactic and the extra-Galactic origin interpretations. Here, we compute the neutrino diffuse emission of our Galaxy on the basis of a recently proposed phenomenological model characterized by radially-dependent cosmic-ray (CR) transport properties. We show that the astrophysical spectrum measured by IceCube experiment can be well explained adding to the diffuse Galactic neutrino flux (obtained with this new model) a extra-Galactic component derived from the astrophysical muonic neutrinos reconstructed in the Northern hemisphere. A good agreement between the expected astrophysical neutrino flux and the IceCube data is found for the full sky as well as for the Galactic plane region.

We present a novel determination of the astrophysical uncertainties associated to the secondary antiproton flux originating from cosmic-ray spallation on the interstellar gas. We select a set of propagation models compatible with the recent B/C data from PAMELA, and find those providing minimal and maximal antiproton fluxes in different energy ranges. We use this result to determine the most conservative bounds on relevant Dark Matter (DM) annihilation channels: We find that the recent claim of a DM interpretation of a gamma-ray excess in the Galactic Center region cannot be ruled out by current antiproton data. Finally, we discuss the impact of the recently released preliminary data from AMS-02. In particular, we provide a reference model compatible with proton, helium and B/C spectra from this experiment. Remarkably, the main propagation parameters of this model are in agreement with the best fit presented in our earlier statistical analyses. We also show that the antiproton-to-proton ratio does not exhibit any significant anomaly at high energy with respect to our predictions.

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

We compute the gamma-ray and neutrino diffuse emission of the Galaxy on the basis of a recently proposed phenomenological model characterized by radially dependent cosmic-ray (CR) transport properties. We show how this model, designed to reproduce both Fermi-LAT gamma-ray data and local CR observables, naturally reproduces the anomalous TeV diffuse emission observed by Milagro in the inner Galactic plane. Above 100 TeV our picture predicts a neutrino flux that is about five (two) times larger than the neutrino flux computed with conventional models in the Galactic Center region (full-sky). Explaining in that way up to $\sim 25 \%$ of the flux measured by IceCube, we reproduce the full-sky IceCube spectrum adding an extra-Galactic component derived from the muonic neutrino flux in the northern hemisphere. We also present precise predictions for the Galactic plane region where the flux is dominated by the Galactic emission.

Cosmic-ray (CR) physics has been a prolific field of research for over a century. The open problems related to CR acceleration, transport and modulation are deeply connected with the indirect searches for particle dark matter (DM). In particular, the high-quality gamma-ray data released by Fermi-LAT are under the spotlight in the scientific community because of a recent claim about a inner Galaxy anomaly: The necessity to disentangle the astrophysical emission due to CR interactions from a possible DM signal is therefore compelling and requires a deep knowledge of several non-trivial aspects regarding CR physics. I review all these connections in this contribution. In the first part, I present a detailed overview on recent results regarding modeling of cosmic-ray (CR) production and propagation: I focus on the necessity to go beyond the standard and simplified picture of uniform and homogeneous diffusion, showing that gamma-ray data point towards different transport regimes in different regions of the Galaxy; I sketch the impact of large-scale structure on CR observables, and -- concerning the interaction with the Heliosphere -- I mention the necessity to consider a charge-dependent modulation scenario. In the second part, all these aspects are linked to the DM problem. I analyze the claim of a inner Galaxy excess and discuss the impact of the non-trivial aspects presented in the first part on our understanding of this anomaly.

Several observations from Fermi-LAT, up to few hundred GeV, and from H.E.S.S., up to $\sim$ 10 TeV, reported an intense $\gamma$-ray emission from the inner part of the Galactic plane. After the subtraction of point-like contributions, the remaining $\gamma$-ray spectrum can provide important hints about the cosmic-ray (CR) population in that region. In particular, the diffuse spectrum measured by both Fermi-LAT and H.E.S.S. in the Galactic Ridge is significantly harder with respect to the rest of the Galaxy. These results were recently interpreted in terms of a comprehensive CR transport model which, adopting a spatial dependent diffusion coefficient and convective velocity, reproduces Fermi-LAT results on the whole sky as well as local CR spectra. We showed as that model predicts a significantly harder neutrino diffuse emission compared to conventional scenarios: The predicted signal is able to account for a significant fraction of the astrophysical flux measured by IceCube. In this contribution, we use the same setup to calculate the expected neutrino flux from several windows in the inner Galactic plane and compare the results with IceCube observations and the sensitivities of Mediterranean neutrino telescopes. In particular, for the ANTARES experiment, we compare the model expectations with the upper limits obtained from a recent unblinded data-analysis focused on the galactic ridge region. Moreover, we also show the expectations from the galactic ridge for the future KM3NeT observatory, whose position is optimal to observe this portion of the sky.

As recently shown, Fermi-LAT measurements of the diffuse gamma-ray emission from the Galaxy favor the presence of a smooth softening in the primary cosmic-ray spectrum with increasing Galactocentric distance. This result can be interpreted in terms of a spatial-dependent rigidity scaling of the diffusion coefficient. The DRAGON code was used to build a model based on such feature. That scenario correctly reproduces the latest Fermi-LAT results as well as local cosmic-ray measurements from PAMELA, AMS-02 and CREAM. Here we show that the model, if extrapolated at larger energies, grasps both the gamma-ray flux measured by MILAGRO at 15 TeV and the H.E.S.S. data from the Galactic ridge, assuming that the cosmic-ray spectral hardening found by those experiments at about 250 GeV/n is present in the whole inner Galactic plane region. Moreover, we show as that model also predicts a neutrino emission which may account for a significant fraction, as well as for the correct spectral shape, of the astrophysical flux measured by IceCube above 25 TeV.

Synchrotron diffuse radiation (SDR) emission is one of the major Galactic components, in the 100 MHz up to 100 GHz frequency range. Its spectrum and sky map provide valuable measure of the galactic cosmic ray electrons (GCRE) in the relevant energy range, as well as of the strength and structure of the Galactic magnetic fields (GMF), both regular and random ones. This emission is an astrophysical sky foreground for the study of the Cosmic Microwave Background (CMB), and the extragalactic microwave measurements, and it needs to be modelled as better as possible. In this regard, in order to get an accurate description of the SDR in the Galaxy, we use - for the first time in this context - 3-dimensional GCRE models obtained by running the DRAGON code. This allows us to account for a realistic spiral arm pattern of the source distribution, demanded to get a self-consistent treatment of all relevant energy losses influencing the final synchrotron spectrum.

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

A GeV gamma-ray excess has possibly been individuated in Fermi-LAT data from the Galactic Center and interpreted in terms of Dark Matter (DM) annihilations, either in hadronic (essentially $b\bar{b}$) or leptonic channels. In order to test this tantalizing interpretation, we address two issues: (i) we improve the computation of secondary emission from DM (Inverse Compton and Bremsstrahlung) with respect to previous works, confirming it to be very relevant for determining the DM spectrum in the leptonic channels, so that any conclusion on the DM nature of the signal critically depends on this contribution; (ii) we consider the constraints from antiprotons on the DM hadronic channel, finding that the uncertainties on the propagation model, and in particular on the halo height, play a major role. Moreover, we discuss the role of solar modulation, taking into account possible charge dependent effects whose importance is estimated exploiting detailed numerical tools. The limits that we obtain severely constrain the DM interpretation of the excess in the hadronic channel, for standard assumptions on the Galactic propagation parameters and solar modulation. However, they considerably relax if more conservative choices are adopted.

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.

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.

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.

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.

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.

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.

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.

After the PAMELA finding of an increasing positron fraction above 10 GeV, the experimental evidence for the presence of a new electron and positron spectral component in the cosmic ray zoo has been recently confirmed by Fermi-LAT. We show that a simple phenomenological model which assumes the presence of a primary electron and positron extra component allows a consistent description of all available data sets. We then describe the most relevant astrophysical uncertainties which still prevent to determine the electron+positron source properties from those data and the perspectives of forthcoming experiments.

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.

In this symposium the Fermi-LAT collaboration released the results of its measurement of the $e^-$ and $e^+$ separate spectra, obtained using the Earth magnetic field. Those results confirm PAMELA finding of an increasing positron fraction above 10 GeV and allow, for the first time independently on the electron spectrum, to reject cosmic ray spallation onto the ISM as the dominant positron production mechanism above that energy. We show, for two different propagation setups, that double component models which were developed to provide a consistent interpretation of PAMELA positron fraction and Fermi-LAT $e^+ + e^-$ data, naturally reproduce also the Fermi-LAT $e^-$ and $e^+$ separate spectra.

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 Fermi Large Area Telescope (LAT) has provided the measurement of the high energy cosmic ray electrons plus positrons (CRE) spectrum with unprecedented accuracy form 20 GeV to 1 TeV. Recently this range has been extended down to ~ 7 GeV. The spectrum shows no prominent features and it is significantly harder than that inferred from several previous experiments. While the reported Fermi-LAT data alone may be interpreted in terms of a single (electron dominated) Galactic component, when combined with other complementary experimental results, specifically the CRE spectrum measured by H.E.S.S., and especially the positron fraction measured by PAMELA, an additional electron and positron component seems to be required. We show that the acceleration of electron-positron pairs in Galactic pulsars may offer a natural (though not unique) interpretation of all those results.

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.

The Fermi-LAT experiment recently reported high precision measurements of the spectrum of cosmic-ray electrons-plus-positrons (CRE) between 20 GeV and 1 TeV. The spectrum shows no prominent spectral features, and is significantly harder than that inferred from several previous experiments. Here we discuss several interpretations of the Fermi results based either on a single large scale Galactic CRE component or by invoking additional electron-positron primary sources, e.g. nearby pulsars or particle Dark Matter annihilation. We show that while the reported Fermi-LAT data alone can be interpreted in terms of a single component scenario, when combined with other complementary experimental results, specifically the CRE spectrum measured by H.E.S.S. and especially the positron fraction reported by PAMELA between 1 and 100 GeV, that class of models fails to provide a consistent interpretation. Rather, we find that several combinations of parameters, involving both the pulsar and dark matter scenarios, allow a consistent description of those results. We also briefly discuss the possibility of discriminating between the pulsar and dark matter interpretations by looking for a possible anisotropy in the CRE flux.

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.