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.

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 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 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 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.

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 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.

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.