Cookie Policy

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

Learn more
Logo HelMod small


Screened Relativistic (SR) Treatment for NIEL Dose

Nuclear and Electronic Stopping Power Calculator

(version 10.13)

Proton High AMS02 small

Neutron Damage Function, NIEL and Hardness Parameter

Fast neutrons, as well as other types of particles, produce defects mainly by displacing silicon atoms from their lattice positions to interstitial locations, i.e., generating vacancy-interstitial pair, the so-called Frenkel pairs. Neutrons with sufficient large energies (i.e. fast neutrons) can cause the displacement of a silicon or any other atom inside the lattice, referred to as primary displacement, from its location. The displaced atom can, in turn, generate secondary displacements, following a cascading effect. In literature (see for instance[Namenson, Wolicki and Messenger (1982), Ougouag et al. (1990), ASTM E722 (2009, 2014 and 2019)]), the damage effect due to displacements induced by neutrons is expressed by the damage function D(E) in units of MeV cm2, also called the displacement kerma function, whose value, for instance, in Si at 1 MeV is the ASTM standard [ASTM E722 (2009, 2014 and 2019)] D(1 MeV)95 MeV mb. The damage function accounts for both the cross section for displacing silicon atoms and the energy released in creating displacements. It is given by:

D(E) = kσk(E) NIMB ONLY neutron0x    (1)

where E is the incoming neutron energy, σk(E) is the cross section for the kth reaction in which fk(E,ER) is the probability that a recoil energy ER is generated and, finally, P(ER) is the partition energy, i.e., the part of the recoil energy deposited in displacements calculated, for instance, in the framework of the Lindhard screened potential scattering theory, based on the Thomas-Fermi model and further developments[Lindhard, Nielsen, Scharff and Thomsen (1963), Coulter and Parkin (1977), Coulter and Parkin (1979), Coulter and Parkin (1980)]. The damage function is given in literature (e.g., see [Ougouag et al. (1990)] and references therein). The damage functions were obtained in [ASTM E722 (2009 & 2014)] for Ed = 25 eV (2009) and Ed = 20.5 eV (2014) in Si and 10 eV in GaAs, where Ed is the is the so-called displacement threshold energy

The energy density Edis is the energy per cm3 deposited through atomic displacements by neutrons, which are characterized by the spectral fluence ϕ(E) in n/cm2 MeV, and is given by:

Edis  =  N NIMB ONLY neutron1x   [MeV cm-3]    (2)

where N is the number of atoms per cm3 in the bulk silicon and Emin is the lowest neutron kinetic energy for displacement. Edis, as computed by means of Eq. (2), is the Non Ionizing Energy Loss (NIEL in units of MeV/cm) deposited by the neutrons per cm2 of the absorber (e.g., treatment in Chapter 4 of [Leroy and Rancoita (2011)], Chapter 7 of [Leroy and Rancoita (2016)]): as long as Emin is 10 keV, Edis varies slightly. In the case of the Triga reactor RC:1 neutron spectrum, the Edis value varies by no more than 0.5 % for Emin values in the energy range up to 10 keV (e.g., see [Consolandi et al. (2006)], [Leroy and Rancoita (2007)] and references therein).

In the current webcalculator the NIEL and NIEL doses are determined by means of the damage functions from [ASTM E722 (2009, 2014 and 2019)].

The displacement damages induced by neutrons can be normalized to displacement damages induced by 1 MeV neutrons by means of the hardness parameter (κ) also termed hardness factor. For fission neutrons with spectral fluence ϕ(E), the 1 MeV equivalent neutron fluence is given by:

Φeq1MeV = κ NIMB ONLY neutron2xdE = κ Φn     (3)


Φeq1MeV is the fluence of 1 MeV neutron needed for generating the same amount of displacement damage compared to the fission fluence Φn.

From Eqs. (2,3), we can derive the hardness parameter κ as:

NIMB ONLY neutron3x

with nA the number of atoms per cm3; from the above equation one finds

κ = NIMB ONLY neutron4xNIMB ONLY neutron5x-1    (4)


and, conversely,

Φeq1MeV = NIMB ONLY neutron4x    (5)



A. I. Namenson, E. A. Wolicki and G.C. Messenger (1982), IEEE Trans on Nucl. Sci. 29 1018.

A. M. Ougouag et al. (1990), IEEE Trans on Nucl. Sci. 37 2219.

ASTM E722-19 (2019). Standard Practice for Characterizing Neutron Fluence Spectra in Terms of an Equivalent Monoenergetic Neutron Fluence for Radiation-Hardness Testing of Electronics (; see also previously published standards: ASTM E722–09 (2009); ASTM E722-14 (2014).

J. Lindhard, V. Nielsen, M. Scharff and P. V. Thomsen (1963), Danske Vidensk. Selsk. Mat.-Fys. Medd. 33 vol. 10, 10

C. A. Coulter and D. M. Parkin (1980), Jou. Nucl. Mat. 88, 249.

C. A. Coulter and D. M. Parkin (1977), Trans. Am. Nucl. Soc. 27, 300.

C. A. Coulter and D. M. Parkin (1979), Jou. Nucl. Mat. 85&86, 611.

C. Leroy and P.G. Rancoita (2011),
Principles of Radiation Interaction in Matter and Detection - 3rd Edition -, World Scientific, Singapore, ISBN-978-981-4360-51-7;



 9167.coverC. Leroy and P.G. Rancoita (2016),
Principles of Radiation Interaction in Matter and Detection - 4th Edition -, World Scientific. Singapore, ISBN-978-981-4603-18-8 (printed); ISBN.978-981-4603-19-5 (ebook);; it is also partially accessible via google books.



 C. Consolandi, P.D’Angelo, G. Fallica, R. Modica, R. Mangoni, S. Pensotti and P.G. Rancoita, (2006), Systematic Investigation of Monolithic Bipolar Transistors Irradiated with Neutrons, Heavy Ions and Electrons for Space Applications, Nucl. Instr. and Meth. in Phys. Res. B 252 (2006), 276, doi:10.1016/j.nimb.2006.08.018

C. Leroy and P.G. Rancoita (2007), Particle Interaction and Displacement Damage in Silicon Devices operated in Radiation Environments Reports on Progress in Physics 70, 493-625, doi:10.1088/0034-4885/70/4/R0