28.11.2014 Views

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

THE EGS5 CODE SYSTEM

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

over the range of energy transfers which still give rise to soft final state secondary electrons, as in:<br />

−<br />

( ) dE±<br />

dx Sub-Cutoff<br />

Atomic Electrons<br />

=<br />

∫ Tmax<br />

0<br />

T dΣ ±<br />

dT<br />

dT . (2.254)<br />

For values of T on the order of the atomic excitation levels, the frequencies and strengths<br />

of the atomic oscillators must be taken into account and the evaluation of this integral is quite<br />

complicated. For values of T large enough that the atomic electrons may be considered free, the<br />

Møller or Bhabha cross sections can be used to describe the scattering. If we let T med be a value<br />

of T which is sufficiently above the atomic excitation level, but is still small compared to T max , we<br />

have ( )<br />

− dE −<br />

dx<br />

=<br />

∫ Tmed<br />

0<br />

T dΣ −<br />

dT<br />

dT + ∫ Tmax<br />

T med<br />

T dΣ Møller<br />

dE<br />

dT (2.255)<br />

and (<br />

− dE +<br />

dx<br />

)<br />

=<br />

∫ Tmed<br />

0<br />

T dΣ +<br />

dT<br />

dT + ∫ Tmax<br />

T med<br />

T dΣ Bhabha<br />

dE<br />

dE (2.256)<br />

When appropriate approximations are made, Equations 2.224 and 2.225 are independent of<br />

T med . We use the formulas recommended by Berger and Seltzer [17] for restricted stopping power 7<br />

which are based on the Bethe-Bloch formula [21, 22, 36]. Note that we have corrected typographic<br />

errors in Berger and Seltzer’s Equations 22-24 for F + (τ, ∆)). Other useful references regarding this<br />

topic are Rohrlich and Carlson[140], Jauch and Rohrlich[83], Turner[173], Sternheimer[159, 160,<br />

161, 162, 163, 165], Evans[53] (for Bhabha and Møller cross sections), Armstrong and Alsmiller[10],<br />

and Messel and Crawford[103]. The formula used in PEGS for the restricted stopping power (i.e.,<br />

due to sub-cutoff electrons) is<br />

(<br />

−X 0<br />

dĔ±<br />

dx<br />

)<br />

Sub-Cutoff<br />

Atomic Electrons<br />

= X 0n2πr0 2m<br />

[<br />

β 2 ln<br />

]<br />

2(τ + 2)<br />

(Īadj/m) + F ± (τ, ∆) − δ<br />

(2.257)<br />

where γ = Ĕ 0 /m = usual relativistic factor,<br />

√<br />

(2.258)<br />

η = γ 2 − 1 = βγ = ˘p 0 c/m,<br />

√<br />

(2.259)<br />

β = 1 − γ −2 = v/c for incident particle, (2.260)<br />

T ′ E = T E /m = K.E. cutoff in electron mass units, (2.261)<br />

τ = γ − 1, (2.262)<br />

y = (γ + 1) −1 (See Bhabha formula), (2.263)<br />

T ′ max = maximum energy transfer<br />

= (τ if positron, τ/2 if electron), (2.264)<br />

∆ = restricted maximum energy transfer<br />

= min(T ′ E, T ′ max), (2.265)<br />

7 The concept of “restricted stopping power” is discussed in detail in the book by Kase and Nelson[87].<br />

74

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!