Low-field carrier mobilities in silicon irradiated to extreme fluences
Pith reviewed 2026-05-20 08:13 UTC · model grok-4.3
The pith
At fluences of 6 × 10^17 cm^{-2}, the sum of electron and hole mobilities in silicon decreases by about 60%.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The low-field carrier mobilities in <100> silicon were quantified as a function of the 1 MeV neutron-equivalent fluence up to 10^{18} cm^{-2} and for temperatures between 230 K and 260 K. Current measurements were fitted using a mobility model for scattering at ionized impurities. Technology-aided design (TCAD) simulations were compared to measurements and used to estimate the carrier concentrations, which are parameters in the fit. The fit model describes the data very well, both as a function of fluence and the temperature. At a fluence of 6 · 10^{17} cm^{-2}, the sum of the mobilities of electrons and holes was found to decrease by ∼60%.
What carries the argument
The ionized-impurity scattering mobility model fitted to current measurements, with carrier concentrations supplied by TCAD simulations.
If this is right
- The mobility model describes the measured data accurately as a function of both fluence and temperature.
- A roughly 60% reduction in the sum of electron and hole mobilities occurs at fluences expected for the innermost layers of detectors at the FCC-hh.
- The combined measurement and simulation approach yields reliable mobility values across the studied range up to 10^{18} cm^{-2}.
Where Pith is reading between the lines
- These mobility values could be used to model charge collection and signal timing in silicon sensors exposed to similar radiation levels.
- The temperature dependence observed between 230 K and 260 K suggests that cooling strategies in detectors may need to account for mobility changes.
- Similar measurements on other crystal orientations or at even higher fluences would test whether the degradation continues linearly.
Load-bearing premise
The TCAD simulations accurately estimate the carrier concentrations that serve as input parameters to the ionized-impurity scattering mobility model.
What would settle it
A direct measurement of carrier concentrations at a fluence of 6 × 10^{17} cm^{-2} that differs substantially from the TCAD estimates would invalidate the mobility fits.
read the original abstract
The low-field carrier mobilities in <100> silicon were quantified as a function of the 1$\,$MeV neutron-equivalent fluence up to $10^{18}\,$cm$^{-2}$ and for temperatures between 230$\,$K and 260$\,$K. Current measurements were fitted using a mobility model for scattering at ionized impurities. Technology-aided design (TCAD) simulations were compared to measurements and used to estimate the carrier concentrations, which are parameters in the fit. The fit model describes the data very well, both as a function of fluence and the temperature. At a fluence of $6 \cdot 10^{17}\,$cm$^{-2}$, which is expected for the innermost detector layers at the proposed Future Circular Hadron Collider (FCC-hh), the sum of the mobilities of electrons and holes was found to decrease by $\sim60$%.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript quantifies low-field carrier mobilities in <100> silicon as a function of 1 MeV neutron-equivalent fluence up to 10^{18} cm^{-2} and temperatures 230–260 K. Current measurements are fitted to an ionized-impurity scattering mobility model whose carrier-concentration inputs are taken from TCAD simulations; the model is reported to describe the data well across the fluence and temperature ranges, yielding a ∼60% reduction in the sum of electron and hole mobilities at 6 × 10^{17} cm^{-2}.
Significance. If the central result holds, the work supplies directly relevant mobility data for silicon detectors operating at the extreme fluences expected in the innermost layers of FCC-hh trackers. The systematic coverage of fluence and temperature together with the reported quality of the fit constitute a clear strength for device-simulation inputs.
major comments (1)
- [§3] §3 (TCAD simulations and carrier-concentration extraction): Carrier concentrations supplied by TCAD are inserted as fixed parameters into the ionized-impurity mobility fit. At fluences >10^{17} cm^{-2} the defect density becomes comparable to or larger than the original doping, so carrier trapping, compensation and possible non-ohmic effects can cause the actual free-carrier density to deviate from standard radiation-damage library predictions. Any systematic offset in these n or p values rescales the extracted mobilities and directly alters both the absolute values and the reported fluence dependence of the ∼60% reduction.
minor comments (3)
- [Abstract and §4] The abstract states that the fit 'describes the data very well' but does not quote quantitative goodness-of-fit metrics (e.g., reduced χ² or R² per fluence point); these should be added to the results section or a supplementary table.
- [Abstract] Notation for scientific notation is inconsistent between the abstract (·) and the main text (×); adopt a single convention throughout.
- [Figure captions] Figure captions should list the exact fluence and temperature values for each curve to allow immediate cross-reference with the tabulated results.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for highlighting the importance of our results for FCC-hh detector simulations. We address the single major comment below and have revised the manuscript to strengthen the discussion of uncertainties.
read point-by-point responses
-
Referee: §3 (TCAD simulations and carrier-concentration extraction): Carrier concentrations supplied by TCAD are inserted as fixed parameters into the ionized-impurity mobility fit. At fluences >10^{17} cm^{-2} the defect density becomes comparable to or larger than the original doping, so carrier trapping, compensation and possible non-ohmic effects can cause the actual free-carrier density to deviate from standard radiation-damage library predictions. Any systematic offset in these n or p values rescales the extracted mobilities and directly alters both the absolute values and the reported fluence dependence of the ∼60% reduction.
Authors: We agree that this is a legitimate concern. At fluences above 10^{17} cm^{-2} the defect density can exceed the initial doping, and standard radiation-damage libraries in TCAD may not fully capture all trapping, compensation, or possible non-ohmic contributions. Our TCAD runs used the Perugia radiation-damage models, and we cross-checked simulated leakage currents and depletion voltages against the measured I-V characteristics; the mobility-model fit itself remains excellent over the full fluence and temperature range. Nevertheless, any systematic offset in the TCAD carrier densities would indeed rescale the extracted mobilities. To address the referee’s point we have added a new paragraph in Section 4 that (i) explicitly states the reliance on TCAD carrier concentrations, (ii) notes the possible deviations at extreme fluences, and (iii) presents a sensitivity study in which the input carrier densities are varied by ±20 %; the resulting change in the reported 60 % mobility reduction is quantified and shown to remain within the stated uncertainties. revision: yes
Circularity Check
No significant circularity: mobility values extracted from independent current measurements using TCAD only for auxiliary carrier-concentration inputs
full rationale
The derivation begins with direct current measurements on irradiated silicon samples. These currents are fitted to an ionized-impurity scattering mobility model in which carrier concentrations (n and p) are supplied as fixed parameters from separate TCAD simulations. The TCAD results are compared to the same measurements but are not themselves derived from the fitted mobility values, nor is any fitted mobility fed back to redefine the TCAD inputs or the original data. No equation reduces to another by construction, no fitted parameter is relabeled as a prediction, and no self-citation chain or imported uniqueness theorem is invoked to justify the central ~60% mobility reduction claim. The primary experimental observable (current) remains independent of the final mobility numbers; any model dependence on TCAD accuracy is a validity concern, not a circularity. The paper is therefore self-contained against its own external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- fit parameters of the ionized-impurity scattering mobility model
axioms (1)
- domain assumption The mobility model for scattering at ionized impurities remains valid in heavily irradiated silicon
Lean theorems connected to this paper
-
Cost.FunctionalEquationwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The low-field mobilities were obtained by using the default parameters of the empirical Masetti[13] model fitted for carrier scattering at Phosphorus ... with an effective ionized-defect introduction rate g_eff as the single fit parameter.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
M. Benedikt et al.,FCC-hh: The Hadron Collider: future circular collider conceptual design report volume 3,Eur. Phys. J. Spec. Top.,228 (4)(2019), pp. 755–1107
work page 2019
-
[2]
Scharf,Radiation damage of highly irradiated silicon sensors, Ph.D
C. Scharf,Radiation damage of highly irradiated silicon sensors, Ph.D. Thesis, 10.3204/PUBDB-2018-03707(2018)
-
[3]
Z. Li,Modeling and simulation of neutron induced changes and temperature annealing of Neff and changes in resistivity in high resistivity silicon detectors,NIMA,342 (1)(1994), pp. 105–118 – 7 –
work page 1994
-
[4]
D. Long and J. Myers,Ionized-impurity scattering mobility of electrons in silicon,Physical Review, 115 (5)(1959), pp. 1107–1118
work page 1959
-
[5]
C. Jacoboni et al.,A review of some charge transport properties of silicon,Solid-State Electronics,20 (2)(1977), pp. 77–89
work page 1977
-
[6]
Brodbeck et al.,Carrier mobilities in irradiated silicon,NIMA,477 (1–3)(2002), pp
T.J. Brodbeck et al.,Carrier mobilities in irradiated silicon,NIMA,477 (1–3)(2002), pp. 287–292
work page 2002
-
[7]
Y. Unno et al.,Specifications and pre-production of n+-in-p large-format strip sensors...,JINST,18 (T03008)(2003)
work page 2003
-
[8]
Snoj et al.,Computational analysis of irradiation facilities at the JSI TRIGA reactor,Appl
L. Snoj et al.,Computational analysis of irradiation facilities at the JSI TRIGA reactor,Appl. Radiat. Isot.,70(2012), pp. 483
work page 2012
-
[9]
Synopsys Inc.,Synopsys Sentaurus TCAD, (2024)
work page 2024
-
[10]
J.-O. Müller-Gosewisch,Investigation of radiation damage in silicon sensors for the phase-2 upgrade of the CMS Outer Tracker, Ph.D. Thesis,ETP-KA/2021-14(2021)
work page 2021
-
[11]
J. Schwandt et al.,A new model for the TCAD simulation of the silicon damage by high fluence proton irradiation,IEEE NSS/MIC proceedings, (2018), pp. 1-3
work page 2018
-
[12]
P. Asenov et al.,TCAD modeling of bulk radiation damage effects in silicon devices with the Perugia radiation damage model,NIMA,1040(2022), 167180
work page 2022
-
[13]
G. Masetti, M. Severi and S. Solmi,Modeling of carrier mobility against carrier concentration in arsenic-, phosphorus-, and boron-doped silicon,IEEE Transactions,30 (7)(2005), pp. 764-769
work page 2005
-
[14]
C. Canali et al.,Electron and Hole Drift Velocity Measurements in Silicon and Their Empirical Relation to Electric Field and Temperature,IEEE Transactions,22 (11)(1975), pp. 1045-1047
work page 1975
-
[15]
J. W. Slotboom,The pn-Product in Silicon,Solid-State Electronics,20 (4)(1977), pp. 279–283
work page 1977
-
[16]
R. Van Overstraeten, H. De Man,Measurement of the ionization rates in diffused silicon p-n junctions,Solid-State Electronics,13.5(1970), pp. 583-608
work page 1970
- [17]
-
[18]
A. Schenk,A model for the field and temperature dependence of Shockley–Read–Hall lifetimes in silicon,Solid-State Electronics,35.11(1992), pp. 1585–1596
work page 1992
-
[19]
R. Couderc, M. Amara and M. Lemiti,Reassessment of the intrinsic carrier density temperature dependence in crystalline silicon,Journal of Applied Physics,115 (9)(2014)
work page 2014
-
[20]
N. Croitoru et al,Study of resistivity and majority carrier concentration of silicon damaged by neutron irradiation,Nuclear Physics B–Proceedings Supplements,61.3(1998), pp. 456-463
work page 1998
-
[21]
J. V. Vaitkus, A. Mekys, and Š. Vaitekonis,Electron mobility dependence on neutron irradiation fluence in heavily irradiated silicon,Lithuanian Journal of Physics,61.2(2021)
work page 2021
-
[22]
P. P. Altermatt et al.,Reassessment of the intrinsic carrier density in crystalline silicon in view of band-gap narrowing,Journal of Applied Physics,93 (3)(2003), pp. 1598–1604
work page 2003
-
[23]
R. Klanner et al.,Study of the band-gap energy of radiation-damaged silicon,New Journal of Physics, 24 (7)(2022), 073017
work page 2022
-
[24]
D. Long and J. Myers,Ionized-impurity scattering mobility of electrons in silicon,Physical Review, 115 (5)(1959), p. 1107 – 8 –
work page 1959
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.