Energy-loss calculator

Enter a particle, its momentum, an absorber and a thickness. This page runs the same instrument as two of our articles — the mean stopping rate from the Bethe formula, the most probable energy loss from the Landau theory (as given in the Particle Data Group review), the straggling distribution for that case, and, in silicon, the electron-hole pairs it creates per micrometre.

Reference tool

Energy-loss calculator: mean, most probable value and straggling

Your particle, your absorber. Mean rate from the Bethe formula with the density-effect correction, the most probable loss from the Landau theory as given in the PDG review, and the straggling curve for your case.

The calculator needs JavaScript. The reference values below are computed by the same code.

Method, constants and limits

Method. Mean rate: Bethe formula, ⟨−dE/dx⟩ = K z² (Z/A) β⁻² [ ½ ln(2mec²β²γ²Tmax/I²) − β² − δ/2 ], without shell or Barkas terms (they exceed a few percent only below βγ ≈ 0.1). Most probable loss: Δp = ξ [ ln(2mec²β²γ²/I) + ln(ξ/I) + 0.200 − β² − δ ], ξ = (K/2) z² (Z/A) x/β². Regime from κ = ξ/Tmax: Landau below 0.01, Vavilov to 10, Gaussian above. The Landau density is evaluated by direct quadrature of its defining integral, not an approximation. Both formulas and their notation follow the Particle Data Group review Passage of particles through matter.

Electron-hole pairs (silicon only). Pairs per micrometre = energy deposited per micrometre ÷ w, w = 3.62 eV (silicon, room temperature; PDG review / silicon-detector literature, e.g. Knoll, Radiation Detection and Measurement). Two figures are reported, not one, because they answer different questions: the mean rate uses ⟨dE/dx⟩ and does not depend on thickness (≈108/µm for a minimum-ionising particle in silicon); the most probable rate uses Δp at the thickness entered above and does depend on it, because Δp is not linear in thickness — it is what a thin detector actually reads out on a single pass, the same Landau/straggling quantity this tool otherwise reports in keV. Quoting only the mean figure for a thin-detector signal is a common error this tool deliberately avoids.

Constants. I, density, Z/A and δ(βγ) per material from NIST ESTAR (SRD 124); K = 4πNAre²mec² and mec² from CODATA; particle masses from the PDG; w = 3.62 eV for silicon as above — not an ESTAR/XCOM quantity, typed directly in the engine with its source. The Bethe engine reproduces NIST PSTAR proton stopping powers to within 1% from 50 MeV to 10 GeV in every material offered; the deviations are on record in the repository tests.

Limits. Heavy charged particles only (no electrons). βγ from 0.1 to about 2×10⁴. Thin absorbers: when the mean loss exceeds 10% of the kinetic energy the constant-velocity assumption fails and the result says so. The e-h pair figures are silicon-only: no sourced w is offered for any other material.

Reference values: a minimum-ionising muon in silicon

Muon at βγ = 3.5 (369.8 MeV/c), where the mean rate in silicon is at its minimum. Computed by the calculator above. The last two columns are the e-h pair yield per micrometre (w = 3.62 eV): the most-probable column changes with thickness because Δp is not linear in x; the mean column does not.

Thickness (µm)Most probable loss Δp (keV)Mean loss (keV)FWHM (keV)Δp / meanκe-h pairs/µm (most probable)e-h pairs/µm (mean)
5011.3119.333.8760.5858.0e-562.48106.8
10023.9638.657.7510.621.6e-466.17106.8
15037.1157.9811.630.642.4e-468.34106.8
20050.5877.315.50.6543.2e-469.87106.8
30078.2211623.250.6754.8e-472.03106.8
500135.3193.338.760.78.0e-474.75106.8
1000284386.577.510.7351.6e-378.44106.8

Method: PDG review Passage of particles through matter (Bethe formula, Landau most probable value). Constants: NIST ESTAR silicon (I = 173 eV, ρ = 2.33 g/cm³), CODATA, PDG masses. Pair-creation energy w = 3.62 eV (silicon, room temperature; PDG review / silicon-detector literature) is not an ESTAR/XCOM quantity and is typed directly in the engine.

Minimum ionisation by material

The minimum of the mean rate ⟨−dE/dx⟩ for a singly charged heavy particle, and the βγ at which it occurs.

Materialρ (g/cm³)I (eV)Z/Aβγ at minimum⟨−dE/dx⟩min (MeV cm²/g)(keV/µm)
Silicon2.331730.49853.421.6590.3865
Germanium5.3233500.44083.261.3720.7301
Carbon (amorphous)2810.49953.921.7430.3486
Aluminium2.6991660.48183.421.6110.4347
Iron7.8742860.46563.441.451.142
Copper8.963220.45643.421.4021.257
Tungsten19.37270.40253.161.1442.207
Lead11.358230.39583.021.1211.273
Argon (gas)0.0016621880.45063.161.5192.524e-4
Water1750.555141.9980.1998
Air (dry, sea level)0.00120585.70.49923.31.8152.186e-4
Plastic scintillator1.03264.70.54144.041.9560.2018
Kapton1.4279.60.51263.961.820.2584
Silicon dioxide2.32139.20.49933.641.6970.3937
Concrete (Portland)2.3135.20.50273.661.7110.3936
Bone (compact, ICRU)1.8591.90.53013.921.8490.3421
Soft tissue (ICRP)172.30.55124.021.9850.1985

Bethe formula with the NIST ESTAR density-effect correction; no shell or Barkas terms. Material constants: NIST ESTAR (SRD 124). Checked against NIST PSTAR to within 1% (50 MeV – 10 GeV protons).

Used in

This calculator is embedded, with the same engine and the same numbers, in two articles: The Landau Distribution and Energy Loss For Ionizing Particles, which explains where the most probable loss comes from, and The Straggling Function. Energy Loss Distribution of Charged Particles in Silicon Layers, which works through the silicon-detector case this page’s electron-hole pair figures are built for. Read either one for the physics; use this page for a bare instrument and a link you can share to one exact result.