(Abridged) The direct precise measurements of spectra of cosmic rays (CR) species in the wide energy range form a basis for propagation models, for interpretation of gamma-ray and microwave observations, and for disentangling possible signatures of new phenomena. Composition and spectra of CR species are vital for studies of galactic nucleosynthesis and high-energy processes in variety of environments and on different scales. Since the beginning of its operation, AMS-02 has delivered outstanding quality measurements of the spectra of CR antiprotons, electrons, positrons, and nuclei: H-O, Ne, Mg, Si. However, the spectra of heavier nuclei are expected to be available only later in the mission, which is currently extended through the end of the International Space Station (ISS) in 2030. Meanwhile, a comparison of published AMS-02 results with data from earlier experiments, such as HEAO-3-C2, indicate that earlier data may be affected by significant undocumented systematic errors. This creates an awkward situation: the new data are still incomplete, while using the old data to make up for the lack of AMS-02 measurements could result in significant errors. In this paper we are using the existing AMS-02 data to test HEAO-3-C2 results. We show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by the AMS-02 collaboration. We are also updating our already published LIS to provide a complete set of LIS from H-Ni in the energy range from 1 MeV/nucleon to ~100-500 TeV/nucleon thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework that is proved to be a reliable tool in deriving the LIS of CR antiprotons, electrons, and nuclei from H to O.

Local interstellar spectra (LIS) of secondary cosmic ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species that may imply a primary lithium component. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for the model-data comparison. The proposed LIS accommodate the low-energy interstellar spectra measured by Voyager 1, HEAO-3, and ACE/CRIS as well as the high-energy observations by PAMELA, AMS-02, and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.

HelMod is a Monte Carlo code developed to describe the transport of Galactic Cosmic Rays (GCRs) through the heliosphere from the interstellar space to the Earth. In the current HelMod version 4 the modulation process, based on Parker's equation, is applied to the propagation of GCRs in the inner and outer heliosphere, i.e., including the heliosheath. HelMod was proved to reproduce protons, nuclei and electrons cosmic rays spectra observed during solar cycles 23-24 by several detectors, for instance, PAMELA, BESS and AMS-02. In particular, the unprecedented accuracy of AMS-02 observations allowed one a better tuning of the description regarding the solar modulation mechanisms implemented in HelMod. In addition, HelMod demonstrated to be capable of reproducing the uxes observed by the Voyager probes in the inner and outer regions of heliosphere up to its border.

Local interstellar spectra (LIS) of primary cosmic ray (CR) nuclei, such as helium, oxygen, and mostly primary carbon are derived for the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined into a single framework that is used to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. The developed iterative maximum-likelihood method uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The interstellar and heliospheric propagation parameters derived in this study are consistent with our prior analyses using the same methodology for propagation of CR protons, helium, antiprotons, and electrons. The resulting LIS accommodate a variety of measurements made in the local interstellar space (Voyager 1) and deep inside the heliosphere at low (ACE/CRIS, HEAO-3) and high energies (PAMELA, AMS-02).

The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997-2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.

The heliospheric modulation model HelMod solve the transport-equation for the Galactic Cosmic Rays propagation through the heliosphere down to Earth. It is based on a 2-D Monte Carlo approach, that includes a general description of the symmetric and antisymmetric parts of the diffusion tensor, thus, properly treating the particle drift effects as well as convection within the solar wind and adiabatic energy loss. The model was tuned in order to fit 1) the data observed outside the ecliptic plane at several distances from the Earth and 2) the spectra observed near the Earth for both, high and low solar activity periods. Great importance was given to description of polar regions of the heliosphere. We present the flux for protons, antiprotons and helium nuclei computed along the solar cycle 23-24 in comparison with experimental observations and prediction for the full solar cycle 24.

Local interstellar spectra (LIS) for protons, helium and antiprotons are built using the most recent experimental results combined with the state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod that provides modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameters optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997-2015. The proposed LIS accommodate both the low energy interstellar CR spectra measured by Voyager 1 and the high energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

The heliospheric modulation model \helmod{} is a two dimensional treatment dealing with the helio-colatitude and radial distance from Sun and is employed to solve the transport-equation for the GCR propagation through the heliosphere down to Earth. This work presents the current version 3 of the \helmod{} model and reviews how main processes involved in GCR propagation were implemented.~The treatment includes the so-called particle drift effects --~e.g., those resulting, for instance, from the extension of the neutral current sheet inside the heliosphere and from the curvature and gradient of the IMF --, which affect the transport of particles entering the solar cavity as a function of their charge sign.~The \helmod{} model is capable to provide modulated spectra which well agree within the experimental errors with those measured by AMS-01, BESS, PAMELA and AMS-02 during the solar cycles 23 and 24.~Furthermore, the counting rate measured by Ulysses at +/- 80degree of solar latitude and 1 to 5\,AU was also found in agreement with that expected by \helmod{} code version 3.

Propagation of cosmic rays (CRs) from their sources to the observer is described mainly as plain diffusion at high energies, while at lower energies there are other physical processes involved, both in the interstellar space and in the heliosphere. The latter was a subject of considerable uncertainty until recently. New data obtained by several CR missions can be used to and the local interstellar spectra (LIS) of CR species that would significantly reduce the uncertainties associated with the heliospheric propagation. In this paper we present the LIS of CR protons and helium outside the heliospheric boundary. The proposed LIS are tuned to accommodate both, the low energy CR spectra measured by Voyager 1, and the high energy observations publicly released by BESS, Pamela, AMS-01 and AMS-02. The proton and helium LIS are derived by combining CR propagation in the Galaxy, as described by GALPROP, with the heliospheric modulation computed using the HelMod Monte Carlo Tool. The proposed LIS are tuned to reproduce the modulated spectra for both, high and low, levels of solar activity.

AMS-02 is a wide acceptance high-energy physics experiment installed on the International Space Station in May 2011 and operating continuously since then. Using the largest number of detected particles in space of any space-borne experiment, it performs precision measurements of galactic cosmic rays fluxes. Detailed time variation studies of Protons, Heliums, Electron and Positron fluxes were presented. The low-rigidity range exhibits a decreasing general trend strongly related to the increase of solar activity, as well local decreases associated with strong solar events.

HelMod is a code evaluating the transport of Galactic cosmic rays through the inner heliosphere down to Earth. It is based on a 2-D Monte Carlo approach and includes a general description of the symmetric and antisymmetric parts of the diffusion tensor, thus, properly treating the particle drift effects. The model has been tuned in order to fit the data observed outside the ecliptic plane at several distances from the Earth and the spectra observed near the Earth for both, high and low solar activity levels. A stand-alone python module, fully compatible with GalProp, was developed for a comprehensive calculation of solar modulation effects, resulting in a newly suggested set of local interstellar spectra.

The cosmic rays modulation inside the heliosphere is well described by a transport equation introduced by Parker in 1965. To solve this equation several approaches were followed in the past. Recently the Monte Carlo approach becomes widely used in force of his advantages with respect to other numerical methods. In the Monte Carlo approach, the transport equation is associated to a fully equivalent set of Stochastic Differential Equations. This set is used to describe the stochastic path of a quasi-particle from a source, e.g., the interstellar medium, to a specific target, e.g., a detector at Earth. In this work, we present both the Forward-in-Time and Backward-in-Time Monte Carlo solutions. We present an implementation of both algorithms in the framework of HelMod Code showing that the difference between the two approach is below 5\% that can be quoted as the systematic uncertain of the Method itself.

We investigate, in terms of production from pulsars and their nebulae, the cosmic ray positron and electron fluxes above $\sim10$ GeV, observed by the AMS-02 experiment up to 1 TeV. We concentrate on the Vela-X case. Starting from the gamma-ray photon spectrum of the source, generated via synchrotron and inverse Compton processes, we estimated the electron and positron injection spectra. Several features are fixed from observations of Vela-X and unknown parameters are borrowed from the Crab nebula. The particle spectra produced in the pulsar wind nebula are then propagated up to the Solar System, using a diffusion model. Differently from previous works, the omnidirectional intensity excess for electrons and positrons is obtained as a difference between the AMS-02 data and the corresponding local interstellar spectrum. An equal amount of electron and positron excess is observed and we interpreted this excess (above $\sim$100 GeV in the AMS-02 data) as a supply coming from Vela-X. The particle contribution is consistent with models predicting the gamma-ray emission at the source. The input of a few more young pulsars is also allowed, while below $\sim$100 GeV more aged pulsars could be the main contributors.

The AMS-02 experiment confirms the excess of positrons in cosmic rays (CRs) for energy above 10 GeV with respect to the secondary production of positrons in the interstellar medium. This is interpreted as evidence of the existence of a primary source of these particles. Possible candidates are dark matter or astrophysical sources. In this work we discuss the possible contribution due to pulsars and their nebulae. Our key assumption is that the primary spectrum of electrons and positrons at the source is the same of the well known photon spectrum observed from gamma-rays telescopes. Using a diffusion model in the Galaxy we propagate the source spectra up to the Solar System. We compare our results with the recent experiments and with the LIS model

We implemented a website to deal with main effects on Cosmic Ray access to the Earth, i.e. the Solar Modulation and the Geomagnetic Field effect. In helmod.org the end user can easily access a web interface to results catalog of the HelMod Monte Carlo Code. This Model uses a Monte Carlo Approach to solves the Parker Transport Equation, obtaining a modulated proton flux for a period (monthly average) between January 1990 and december 2007. geomagsphere.org is instead based on GeoMag Backtracing Code, that solves the Lorentz equation with a Runge-Kutta method of 6th order, and, reversing charge sign and velocity, reconstruct particle trajectories in the Earth Magnetosphere back in time. We use last models of internal (IGRF-11) and external (Tsyganenko 1996 -T96- and 2005 -T05-) field components valid up to 2015. Particles are backtraced to the outer (magnetopause) or inner boundary to separate Primary (allowed trajectory) from Secondary (forbidden) Cosmic Rays. This code has been used both for reproducing known effects as East-West effect and rigidity cutoff calculations. In geomagsphere.org the user can choose the external field model from Tsyganenko (T96 or T05) and obtain for a fixed position and date from 1st Jan. 1968 (T96) and 1st Jan. 1995 (T05) respectively till 31$^{st}$ Dec 2012, the vertical rigidity cutoff estimation obtained with the backtracing technique with a rigidity step of 0.1 GV. For a more precise calculation (0.01 GV), requiring more CPU time, results are sent to the user by email (mail model)

The Cosmic Rays propagation was studied in details using the HelMod-2D Monte Carlo code, that includes a general description of the diffusion tensor, and polar magnetic-field. The Numerical Approach used in this work is based on a set of Stochastic Differential Equations fully equivalent to the well know Parker Equation for the transport of Cosmic Rays. In our approach the Diffusion tensor in the frame of the magnetic field turbolence does not depends explicitly by Solar Latitude but varies with time using a diffusion parameter obtained by Neutron Monitors. The parameters of the Model were tuned using data during the solar Cycle 23 and Ulysses latitudinal Fast Scan in 1995. The actual parametrization is able to well reproduce the observed latitudinal gradient of protons and the southward shift of the minimum of latitudinal intensity. The description of the model is also available online at website www.helmod.org. The model was then applied on Pamela/Ulysses proton intensity from 2006 up to 2009. The model during this 4-year continous period agree well with both PAMELA (at 1 AU) and Ulysses data (at various solar distance and solar latitude). The agreement improves when considering the ratio between this data. Studies done also with particles with different charge (e.g. electrons) allow us to explain the presence (or not) of protons and electrons latitudinal gradients observed by Ulysses during the Latitudinal Fast Scan in 1995 and 2007.

We used a backtracing code to reconstruct particle trajectory inside the Earth Magnetosphere during the last solar active period (2011 and 2012) when very high Solar Wind pressure values were measured. We compared our results on AMS-02 proton and electron data with 2 different External Field models, namely Tsyganenko 1996 (T96) and 2005 (T05), both for quiet (defined as the periods when the solar wind pressure is below the average value, set at 2nPa) and active periods. Although T05 has been specifically designed for storm events, at high energy the particle trajectory is similar for the two models. For instance at rigidities larger than 50 GV, the RMS of angular difference between reconstructed asymptotic direction outside the Magnetosphere is of the order of the millirad, while it increases at intermediate energies. We also confirmed, as a function of the pointing direction, the well known East-West effect on the trajectory of primary particles and on the access solid angle on the AMS detector.

Several cosmic ray experiments have measured the positron fraction up to few hundred GeV. Their data have revealed an excess of positrons above 10 GeV that is not consistent with the secondary production of these particles in the interstellar medium. A primary source like dark matter or astrophysical sources (e.g pulsars and their nebulae) were considered to account for such an excess. In this paper we analyse the possibility of a primary positron production due to pulsars. Under the assumption of equal initial spectra at the source for positrons, electrons, and gamma-rays we study the propagation of particle spectra using a diffusion model in the Galaxy. We focused our analysis on the Vela and Crab pulsars and their associated nebulae, which are well observed in gamma-rays. Comparison with experimental data is reported. The propagated positron and electron spectra generated from these sources result in a positron ratio, which is largely inconsistent with the excess observed by PAMELA and AMS

Energetic (suprathermal) solar particles, accelerated in the interplanetary medium, contribute to the solar wind pressure, in particular during high solar activity periods. We estimated the effect of the increase of solar wind pressure due to suprathermal particles on magnetospheric transmissivity of galactic cosmic rays in the case of one recent solar event.

The cosmic rays differential intensity inside the heliosphere, for energy below 30 GeV/nuc, depends on solar activity and interplanetary magnetic field polarity. This variation, termed solar modulation, is described using a 2-D (radius and colatitude) Monte Carlo approach for solving the Parker transport equation that includes diffusion, convection, magnetic drift and adiabatic energy loss. Since the whole transport is strongly related to the interplanetary magnetic field (IMF) structure, a better understanding of his description is needed in order to reproduce the cosmic rays intensity at the Earth, as well as outside the ecliptic plane. In this work an interplanetary magnetic field model including the standard description on ecliptic region and a polar correction is presented. This treatment of the IMF, implemented in the HelMod Monte Carlo code (version 2.0), was used to determine the effects on the differential intensity of Proton at 1\,AU and allowed one to investigate how latitudinal gradients of proton intensities, observed in the inner heliosphere with the Ulysses spacecraft during 1995, can be affected by the modification of the IMF in the polar regions.

The differential intensities of Cosmic Rays at Earth were calculated using a 2D stochastic Montecarlo diffusion code and compared with observation data. We evaluated the effect of stretched and compressed heliospheres on the Cosmic Ray intensities at the Earth. This was studied introducing a dependence of the diffusion parameter on the heliospherical size. Then, we found that the optimum value of the heliospherical radius better accounting for experimental data. We also found that the obtained values depends on solar activity. Our results are compatible with Voyager observations and with models of heliospherical size modulation.

A propagation model of galactic cosmic protons through the Heliosphere was implemented using a 2-D Monte Carlo approach to determine the differential intensities of protons during the solar cycle 23. The model includes the effects due to the variation of solar activity during the propagation of cosmic rays from the boundary of the heliopause down to Earth's position. Drift effects are also accounted for. The simulated spectra were found in agreement with those obtained with experimental observations carried out by BESS, AMS and PAMELA collaborations. In addition, the modulated spectrum determined with the present code for the year 1995 exhibits the latitudinal gradient and equatorial southward offset minimum found by Ulysses fast scan in 1995.

Spectra of Galactic Cosmic Rays (GCRs) measured at the Earth are the combination of several processes: sources production and acceleration, propagation in the interstellar medium and propagation in the heliosphere. Inside the solar cavity the flux of GCRs is reduced due to the solar modulation, the interaction which they have with the interplanetary medium. We realized a 2D stochastic simulation of solar modulation to reproduce CR spectra at the Earth, and evaluated the importance in our results of the Local Interstellar Spectrum (LIS) model and its agreement with data at high energy. We show a good agreement between our model and the data taken by AMS-01 and BESS experiments during periods with different solar activity conditions. Furthermore we made a prediction for the flux which will be measured by AMS-02 experiment.

We implemented a quasi time-dependent 2D stochastic model of solar modulation describing the transport of cosmic rays (CR) in the heliosphere. Our code can modulate the Local Interstellar Spectrum (LIS) of a generic charged particle (light cosmic ions and electrons), calculating the spectrum at 1AU. Several measurements of CR antiparticles have been performed. Here we focused our attention on the CR antiproton component and the antiproton over proton ratio. We show that our model, using the same heliospheric parameters for both particles, fit the observed anti-p/p ratio. We show a good agreement with BESS-97 and PAMELA data and make a prediction for the AMS-02 experiment.

Galactic Cosmic Rays (GCRs) are mainly protons confined in the galactic magnetic field to form an isotropic flux inside the galaxy. Before reaching the Earth orbit they enter the Heliosphere and undergo diffusion, convection, magnetic drift and adiabatic energy loss. The result is a reduction of particles flux at low energy (below 10 GeV), called solar modulation. We realized a quasi time-dependent 2D Stochastic Simulation of Solar Modulation that is able to reproduce CR spectra once known the Local Interstellar Spectrum (LIS). We were able to estimate the different behaviors associated to the polarity dependence of the Heliospheric modulation for particles as well as for antiparticles. We show a good agreement with the antiproton/proton ratio measured by AMS-01, Pamela, BESS, Heat and Caprice and we performed a prediction for the AMS-02 Experiment.

The solar modulation, a combination of diffusion, convection, magnetic drift and energy loss inside the heliosphere is usually seen as a depletion in the Galactic cosmic ray (CR) flux at low energy (less than 10 GeV/nuc). Antiparticles such as antiprotons or positrons undergo the same processes of respective particles but with a different magnitude depending on the Solar magnetic field polarity. For electrons and positrons, due to the small mass, energy loss mechanisms as inverse compton, synchrotron, bremsstrahlung and ionization have to be taken into account, together with the typical adiabatic losses considered in the heliosphere. We developed a Monte Carlo stochastic simulation with the aim to compare the solar modulation of particles and antiparticles in the same observation period. We are able to estimate the different behaviours associated to the charge sign dependent processes of the heliospheric modulation. We compared the simulated positron fraction with measurements performed by AMS-01 and PAMELA. We also present the prediction for the AMS-02 experiment.