The Straggling function. Energy Loss Distribution of charged particles in silicon layers

The statistical nature of the ionizing process during the passage of a fast charge particle through matter results in a large fluctuations of the energy loss (Δ) in absorber which are thin compared with the particle range.
The number of electronhole pairs (J) is related to Δ by the expression J=Δ / P, where P is a proportional factor, for the silicon equal to 3.68 eV. Both Δ and J are stochastic quantities. The probability function f(Δ) and φ(J) are usually called energy loss distribution or straggling function.
They may be characterized schematically by the position of the maximum of distribution function (Δp) and the full width at half maximum (w). The Δp is located at lower value compared to the mean energy loss obtained from BetheBloch ‹Δ› (Fig.1).
(click here to download the C++ BetheBloch calculator).
Fig.1 Stopping power for positive muons in cupper
Theoretical
calculation of this distribution
have been carried out by
Landau
and Vavilov , each of these solution
however have a different region
of applicability and the distinguishing
parameter in all these theories
is the ratio k= ‹Δ›
/ E_{max} where E_{max}
is the maximum transferable energy
in a single collision. Landau solved
this problem for all the situation
where k ≤ 0.01 deriving the expected
energy loss distribution by solving
an integral transport equation:
Here f(x,Δ) represents the distribution probability that the incident particle will lose an amount Δ of energy on traversing a layer of thickness x. W(E)dE denotes the probability per unit path length of a collision transferring energy E to an electron in the material. The function W(E)dE is not generally known but Landau was able to derive an approximate solution by using the free electron (Rutherford) cross section:
with
keV μm1 in silicon.
In the previous equation, z is the charge of incident particle, Z and A are the atomic number and weight of the material and ρ is the density. The Landau distribution is therefore given by [1]
with φ(λ) a universal function of the variable λ only
and
where CE is the Euler constant equal to 0.5772. The Landau distribution, fL(Δ), is asymmetric with a tail extending to Emax with a maximum for λ=0.229 and w=4.018ξ. This tail is mainly due to fast emitted δrays. The energy loss corresponding to the maximum of the function fL(Δ) is the most probable energy loss [2]
where I is the mean excitation potential and δ the density correction. (click here to download the Matlab Landau Fit)
Fig.2 Stopping power for positive muons in cupper Straggling functions in silicon for 500MeV pions, normalized to unity at the most probable value Δp/x . The width w is the full width at half maximum.
Subsequently Vavilov [3] derived an improved solution which takes in account the spin of the incident particle and introduced the physical limits coming from Emax. For the collision cross section Vavilov used the form
Vavilov also demonstrated that his solution tends to the Landau function for k ≤ 0.01, region where ‹Δ› is approximated to ξ. For instance for 300 um silicon detectors, the convergence is already achieved by protons with momenta larger than 550 MeV/c. For k > 10, the number of collision is very large and, for the Central Limit Theorem, the Vavilov function coincides with Gaussian distribution, where the standard deviation is
Fig. 3 Most probable energy loss in silicon, scaled to the mean loss of a minimum ionizing particle,388eV/μm (1.66 MeVcm2/g).
LandauVavilov solutions have been derived under the assumption that scatterings occur on quasifree electrons. Therefore they neglect the electron binding energies. This assumption is a valid approximation only for close collisions.
Further corrections to the theory taking into account the fact that the electrons in the material are not free have been attempted by Blunck and Leisegang [4], Shulek [5] and Bichsel [6].
In the case of solid media these calculations are aimed at the experimental observation that while the most probable energy loss agrees rather well with the prediction of the theory the width of the distribution is broader than expected and cannot be accounted for by electronic noise or imperfect resolution. The effect is particularly noticeable for very thin absorbers [7], of the order of a few hundred μm or less. The modified energy loss distribution can be improved by using a modified cross section to take into account the electron binding energy and the atomic shell structure. The modified energy straggling function can be expressed as:
.
In other words the experimentally observed energy spectrum can be calculated by convoluting the Landau distribution with a normal distribution of variance δ_{2}. The results of the convolution is a broader distribution with a peak value usually increased by a small amount compared to the Landau theory. From the experimental observations the δ2 value has been estimated:
where f_{i} is the fraction of electrons in that shell. In silicon this approaches a value
The resulting improved energy loss distribution has an overall value of w, which is roughly given by:
√[(4.018 ξ )^{^2} + ( 5.56 δ_{2} )]
As the material thickness decreases, δ2 becomes more and more the dominant term, which determines the overall w of the straggling function. Conversely, it is not expected to provide an additional broadening of the distribution at large thicknesses.
As an example, in 300 μm thick silicon detector we get ξ is ~ 5.34 keV, √(δ_{2}) is ~ 5.76 keV, and an overall w is 25.4 keV (i.e., 18% larger than Landau width) for a relativistic β~ 1 and z = 1 particle in agreement with the experimental data.
Fig. 4 Most probable value in terms of electronhole pair generated by electron in silicon, for different traversed thickness (a). The line in (b) represent the MPV of the electronhole pairs generated in the range of a 300500MeV by the electrons in function of the thickness 
References

L. Landau , On the
Energy Loss of Fast
Particles by Ionization,
J. Phys.
USSR 8 (1944)
201.

Particle Data Group
PDG,
Passage of particles
through matter,
Nuclear and Particle
Physics, vol. 33,
no. 27, pp. 258270,
July 2006.

P.V.Vavilov,
Ionization losses
of high energy heavy
particles,
Soviet Physics JETP,
5:749, 1957.

S. Meroli et al.,
Energy loss measurement
for charged particles
in very thin silicon
layers
, JINST,
6 P06013
doi:
10.1088/17480221/6/06/P06013

Blunck and S. Leisegang,
Zum Energieverlust
schneller Elektronen
in dünnen Schichten,
Z. Physik 128 (1950)
500.

P. Shulek at al.,
Fluctuations of
Ionization Loss,
Sov. J. Nucl.
Phys
4 (1967) 400.

H. Bichsel,
Straggling of Heavy
Charged Particles:
Comparison of Born
HydrogenicWaveFunction
Approximation with
FreeElectron Approximation,
Phys.
Rev. B1 (1970)
2854

H. Esbensen at al.,
Random and channeled
energy loss in thin
germanium and silicon
crystals for positive
and negative 215GeV/c
pions, kaons, and
protons,
Phys.
Rev. B18 (1978)
103939
If you like... share with your friends 