Cookie Policy

Welcome to SR-NIEL. This site uses cookies. Read our policy.

Learn more
Logo HelMod small

SR-NIEL – 7

Screened Relativistic (SR) Treatment for NIEL Dose

Nuclear and Electronic Stopping Power Calculator

(version 10.16)

Proton High AMS02 small

The occurance probability of events due to large energy depositions in a device - like those commonly referred to as single event effects (SEE) – is usually expressed by means of the cross section as a function of mass electronic stopping power - or mass restricted energy loss - of the device material (e.g., see Sects. 11.4-11.4.9 of [Leroy and Rancoita (2016)] and references therein).
As discussed in this webpage, one can re-express that device cross section as a function of the energy of incoming particles. In this way, the overall events occurrence can be estimated using the entire range of particle spectral fluence. In fact, nowadays the spectral fluences of low- and high-Z cosmic rays were observed up to very high energies (e.g., see here).

By exploiting experimental results, the combined effort of the physicists involved with the Galprop model for propagation in galaxy and HelMod model for the propagation in heliosphere, the local interstellar spectra (LIS) for Galactic Cosmic Rays species  up to Z=28 (Nickel) were derived. HelMod model allows one to derived the ion modulated spectra from the corresping LIS. These spectra are available as input for the calculation of the occurance probability of events due to a large energy deposition (SEE).

The following link allows one to access to the Probability of events occurrence due to a large energy deposition in devices calculator for space missions and cosmic rays fluences from HelMod model:

- Occurance probability of events due to a large energy deposition (SEE) for space missions and cosmic rays fluences at 1 AU from HelMod Model when particles impinge perpendicularly on the device front face

 

How to use this calculators for probability of events occurrence due to a large energy deposition in devices

This tool calculates the Probability of events occurrence due to a large energy deposition in devices for a spectral fluence incident on a material.

The input parameters and options for the tool are described below. When the input form has been completed, pressing the "CALCULATE" button will start the calculation and open the "Results" page (allow for pop-up in your browser settings). The result page will be also linked at the bottom of the calculator page.

 

Input Parameters:

- Device material
- SEE Cross section (Weibull function parameters)
- Energy Loss Type
- Incident particle
- Spectral fluence from HelMod Model
- Mission time period
- Shielding absorber
- Averaged traversed path
- Number of steps

Currently, as an example, the default parameters for the Weibull function are those for the device treated in this webpage.

 

Device Material

In the section "Device Material Selection" it is possible to specify an User Defined device material or a predefined Compound. User can also select the device as a gas, this is allowed only for single element and natural gas target (H, He, N, O, F, Ne, Cl, Ar, Kr, Xe, Rn).

The stopping power in target gases is usually higher than that in an equivalent solid target. The Gas/ Solid correction disappears for higher velocity ions with energies above 2 MeV/amu. But at lower velocities the effect can be quite large - almost a 2 times change in stopping bacause of the Phase effect near the Bohr velocity, 25 keV/amu.

in the User Defined section individual elements can be selected as well as the composition of the device material choosing the number of elements in the compound. The required parameters for each element are:

- Atomic number (Z)/Chemical symbol
- Stoichiometric index or element fraction

Electronic Stopping Power for User Defined Compounds

Electronic Stopping Power for User Defined Compounds can be determined by means of Bragg's additivity rule, i.e., the overall Electronic Stopping Power in units of MeV cm2/g (i.e., the mass electronic stopping power) is obtained as a weighted sum in which each material contributes proportionally to the fraction of its atomic weight. For instance, in case of a GaAs medium ones obtains (e.g., Eq. (2.20) at page 15 in [ICRUM (1993)]):

6nsnew

where 7nsnew and AGa [AAs] are the Electronic Stopping Power (in units of MeV cm2/g) and the atomic weight of Gallium [Arsenic], respectively.

For few compounds belonging to the ICRU list the parameters employed for the energy loss formula (including those for the densiity effect) are reported in Table II of Sternheimer et al. (1984).
Further discussion can be found in the help of the SR-framewrok electronic stopping (see also results page).

 

SEE Cross section

This section define parameters that are included in the Weibull function:

weibull

where x is the stopping power (or the restricted energy loss) in MeV cm2 mg-1, A is the saturation cross section, x0 is the stopping power (or LET) threshold parameter, PW is the width parameter and s is a dimensionless exponent.

 

Energy Loss type

The user can select 3 different energy loss types:

- SRIM Electronic Stopping Power
- SR-Framework Electronic Stopping
- SR-Framework Restricted Energy Loss

When restricted energy loss is selected, the Effective detectable energy W0 can be modified: the current default value is 100 keV.

For the conversion of the cross section as a function of the incoming particle energy, the input cross section is assumed to be provided as a function of the corresponding type of the selected energy loss (i.e., SRIM Electronic stopping power, SR-Framework Electronic stopping, SR-Framework Restricted energy loss). One has to remark that SRIM Electronic stopping power and SR-Framework Stopping power agree within 5% (or better) up to the maximum energy obtained from SRIM Module.exe (i.e., 5 GeV/amu), see discussion here.
As discussed at page 70 of [Leroy and Rancoita (2016)] for all incoming particle energy for which the maximum transferred energy W (see Eq. (1.29) of [Leroy and Rancoita (2016)]) is lower than or equal to W0 there is no difference in between the SR-Framework electronic stopping power and the restricted energy loss (see Figure 1).

limit plot

Figure 1. Energy limit (in MeV/amu) below which there is no difference in between the SR-Framework electronic stopping power and the restricted energy loss as a function of W0.

 

 

Incident Particle

In the calculator panel, using the pull down menu, the user can select the species of the incident particle, either a proton or one of the elemental ions.
The user can also modify the mass (in amu) of the incident particle (e.g., for all isotopes one can refer to this page): the default mass is that of the most abundant isotope (MAI), except for proton and alpha particle masses. Further information are available at the following webpage

 

Spectral Fluence from HelMod Model

This tools evaluates the modulated spectrum at Earth orbit for theorical LIS evaluated with GALPROP. The particle transport is described using the HelMod Model. For information about the computation method please refear to Monte Carlo Integration page.

The user can select the wished Cosmic Rays Species with atomic number up to Z=28.HelMod model spectra are averaged over a Carrington rotation duration (e.g., 27.2753 days).

However, incoming particle directions are assumed to be identical to those employed in measuring the SEE cross section, so that the traversed path inside the device is not modified.

perp

One has to remark (see webpage) how the physical mechanism of energy deposition per unit length in a medium is such that the resulting computed value of restricted energy loss depends on the overall traversed path to which W0 is related. The electronic stopping power, on the other hand, does not depend on the incoming particle direction with respect to the medium surface.

 

Mission time period

User can select the cosmic rays species and the time interval for the interested simulation.

 

Shielding absorber

In the section "Target Selection" it is possible to specify an User Defined target material or a predefined Compound material. User can also select the target as gas, this is allowed only for single element and natural gas target (H, He, N, O, F, Ne, Cl, Ar, Kr, Xe, Rn). 

in the User Defined section individual elements can be selected as well as the composition of the device material choosing the number of elements in the compound. The required parameters for each element are:

- Atomic number (Z)/Chemical symbol
- Stoichiometric index or element fraction

For few compounds belonging to the ICRU list the parameters employed for the energy loss formula (including those for the densiity effect) are reported in Table II of Sternheimer et al. (1984).
Further discussion can be found in the help of the SR-framewrok electronic stopping (see also results page).

 

Averaged Traversed Path

This input define the averaged traversed path by a particle fluence. The average traversed path refers to the mean path of isotropically distributed cosmic rays crossing a shielding absorber before impinging onto the device. The averaged traversed path is expressed in [g cm] and is given by:

(absorber thickness in cm) x (absorber density in g cm-3) .

User can choose to use unshielded HelMod model spectral fluence by setting the traversed path to zero (0).

 

Number of Steps

In the web Calculator, using the pull down menu, the user can select the number of steps of the calculation - i.e., the traversed path is divided by the number of the steps.
The results of each steps is used as input for the following one to obtain the final result for the total traversed path.

The minimum traversed path allowed for each step is 100 μm of Si (2.3x10-2 [g cm-2]).

 

Result

The result page contains the input spectral fluence and the SEE cross section as a function of energy. The spectral fluence is interpolated and divided in 1000 energy bins to obtain the total number of SEE. The calculation is performed with three types of stopping power (SRIM, SR electronic stopping power and Restricted energy loss). The energy, fluence, cross section and spectral number of SEE are reported in the table.