• ### Deciphering the local Interstellar spectra of primary cosmic ray species with HelMod(1804.06956)

April 19, 2018 astro-ph.HE
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).
• ### HelMod in the works: from direct observations to the local interstellar spectrum of cosmic-ray electrons(1801.04059)

Jan. 12, 2018 astro-ph.HE
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.
• ### Solution of heliospheric propagation: unveiling the local interstellar spectra of cosmic ray species(1704.06337)

April 20, 2017 astro-ph.HE
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.
• ### Cosmic Rays Propagation with HelMod: Difference between forward-in-time and backward-in-time approaches(1511.00541)

Aug. 13, 2015 physics.space-ph, astro-ph.HE
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.
• ### NIEL Dose Dependence for Solar Cells Irradiated with Electrons and Protons(1312.0402)

The investigation of solar cells degradation and the prediction of its end-of-life performance is of primary importance in the preparation of a space mission. In the present work, we investigate the reduction of solar-cells' maximum power resulting from irradiations with electrons and protons. Both GaAs single junction and GaInP/GaAs/Ge triple junction solar cells were studied. The results obtained indicate how i) the dominant radiation damaging mechanism is due to atomic displacements, ii) the relative maximum power degradation is almost independent of the type of incoming particle, i.e., iii) to a first approximation, the fitted semi-empirical function expressing the decrease of maximum power depends only on the absorbed NIEL dose, and iv) the actual displacement threshold energy value (Ed=21 eV) accounts for annealing treatments, mostly due to self-annealing induced effects. Thus, for a given type of solar cell, a unique maximum power degradation curve can be determined as a function of the absorbed NIEL dose. The latter expression allows one to predict the performance of those solar cells in space radiation environment.
• ### Nuclear and Non-Ionizing Energy-Loss for Coulomb Scattered Particles from Low Energy up to Relativistic Regime in Space Radiation Environment(1011.4822)

In the space environment, instruments onboard of spacecrafts can be affected by displacement damage due to radiation. The differential scattering cross section for screened nucleus--nucleus interactions - i.e., including the effects due to screened Coulomb nuclear fields -, nuclear stopping powers and non-ionization energy losses are treated from about 50 keV/nucleon up to relativistic energies.
• ### CH+ depletion by atomic hydrogen: accuracy of new rates in photo-dominated and self-shielded environments(1306.4857)

Dec. 20, 2013 astro-ph.GA
A detailed quantum analysis of a ionic reaction with a crucial role in the ISM is carried out to generate ab initio reactive cross sections with a quantum method. From them we obtain the corresponding CH+ depletion rates over a broad range of temperatures. The new rates are further linked to a complex chemical network that shows the evolution in time of the CH+ abundance in photodissociation region (PDR) and molecular cloud (MC) environments. The evolutionary abundances of CH+ are given by numerical solutions of a large set of coupled, first-order kinetics equations by employing the new chemical package KROME. The differences found between all existing calculations from low-T experiments are explained via a simple numerical model that links the low-T cross section reductions to collinear approaches where nonadiabatic crossings dominate. The analysis of evolutionary abundance of CH+ reveals that the important region for the depletion reaction of this study is that above 100 K, hence showing that, at least for this reaction, the differences with the existing low-temperature experiments are of essentially no importance within the astrochemical environments. A detailed analysis of the chemical network involving CH+ also shows that a slight decrease in the initial oxygen abundance might lead to higher CH+ abundance since the main chemical carbon ion depletion channel is reduced in efficiency. This simplified observation might provide an alternative starting point to understand the problem of astrochemical models in matching the observed CH+ abundances.
• ### Geomagnetic Backtracing: A comparison of Tsyganenko 1996 and 2005 External Field models with AMS-02 data(1307.5192)

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.
• ### Suprathermal particle addition to solar wind pressure: possible influence on magnetospheric transmissivity of low energy cosmic rays?(1307.5195)

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.
• ### GeoMag and HelMod webmodels version for magnetosphere and heliosphere transport of cosmic rays(1307.5196)

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)
• ### Cosmic Ray Modulation studied with HelMod Monte Carlo tool and comparison with Ulysses Fast Scan Data during consecutive Solar Minima(1307.5199)

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.
• ### Latitudinal Dependence of Cosmic Rays Modulation at 1 AU and Interplanetary-Magnetic-Field Polar Correction(1212.1559)

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.
• ### Electrical Characterization of SiPM as a Function of Test Frequency and Temperature(1112.0107)

Silicon Photomultipliers (SiPM) represent a promising alternative to classical photomultipliers, for instance, for the detection of photons in high energy physics and medical physics. In the present work, electrical characterizations of test devices - manufactured by ST Microelectronics - are presented. SiPMs with an area of 3.5x3.5 micron^2 and a cell pitch of 54 micron were manufactured as arrays of 64x64 cells and exhibiting a fill factor of 31%. The capacitance of SiPMs was measured as a function of reverse bias voltage at frequencies ranging from from 20 Hz up to 1 MHz and temperatures from 300 K down to 85 K. While leakage currents were measured at temperatures from 400 K down to 85 K. Thus, the threshold voltage - i.e., voltage corresponding to that at which the multiplication regime for the leakage current begins - could be determined as a function of temperature. Finally, an electrical model suited to reproduce the dependence of the frequency dependence of capacitance is presented.
• ### Nuclear and Non-Ionizing Energy-loss of Electrons with Low and Relativistic Energies in Materials and Space Environment(1111.4042)

The treatment of the electron-nucleus interaction based on the Mott differential cross section was extended to account for effects due to screened Coulomb potentials, finite sizes and finite rest masses of nuclei for electrons above 200 keV and up to ultra high energies. This treatment allows one to determine both the total and differential cross sections, thus, subsequently to calculate the resulting nuclear and non-ionizing stopping powers. Above a few hundreds of MeV, neglecting the effect due to finite rest masses of recoil nuclei the stopping power and NIEL result to be largely underestimated; while, above a few tens of MeV the finite size of the nuclear target prevents a further large increase of stopping powers which approach almost constant values.
• ### Systematic Investigation of Solar Modulation of Galactic Protons for Solar Cycle 23 using a Monte Carlo Approach with Particle Drift Effects and Latitudinal Dependence(1110.4315)

Oct. 19, 2011 astro-ph.SR, astro-ph.EP
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.
• ### Proton Modulation in the Heliosphere for Different Solar Conditions and Prediction for AMS-02(1012.0938)

March 17, 2011 astro-ph.EP
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.