• 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • We propose a model for non isothermal ferromagnetic phase transition based on a phase field approach, in which the phase parameter is related but not identified with the magnetization. The magnetization is split in a paramagnetic and in a ferromagnetic contribution, dependent on a scalar phase parameter and identically null above the Curie temperature. The dynamics of the magnetization below the Curie temperature is governed by the order parameter evolution equation and by a Landau-Lifshitz type equation for the magnetization vector. In the simple situation of a uniaxial magnet it is shown how the order parameter dynamics reproduces the hysteresis effect of the magnetization.
  • A Ginzburg-Landau model for the macroscopic behaviour of a shape memory alloy is proposed. The model is one-dimensional in essence, in that we consider the effect of the martensitic phase transition in terms of a uniaxial deformation along a fixed direction and we use a scalar order parameter whose equilibrium values describe the austenitic phase and the two martensitic variants. The model relies on a Ginzburg-Landau free energy defined as a function of macroscopically measurable quantities, and accounts for thermal effects; couplings between the various relevant physical aspects are established according to thermodynamic consistency. The theoretical model has been implemented within a finite-element framework and a number of numerical tests are presented which investigate the mechanical behaviour of the model under different conditions; the results obtained are analysed in relation to experimental evidences available in literature. In particular, the influence of the strain-rate and of the ambient conditions on the response of the model is highlighted.
  • 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.
  • We consider a model describing the behavior of a mixture of two incompressible fluids with the same density in isothermal conditions. The model consists of three balance equations: continuity equation, Navier-Stokes equation for the mean velocity of the mixture, and diffusion equation (Cahn-Hilliard equation). We assume that the chemical potential depends upon the velocity of the mixture in such a way that an increase of the velocity improves the miscibility of the mixture. We examine the thermodynamic consistence of the model which leads to the introduction of an additional constitutive force in the motion equation. Then, we prove existence and uniqueness of the solution of the resulting differential problem.
  • We propose a model describing the liquid-vapour phase transition according to a phase-field approach. The model takes up a setting proposed by the second author, where a phase field is introduced whose equilibrium values 0 and 1 are associated to the liquid and vapour phases. The phase field obeys Ginzburg-Landau equation and enters the constitutive relation of the density, accounting for the sudden density jump occurring at the phase transition. In this paper we concern ourselves with the extension of this model to take into account the existence of the critical point in the coexistence line.