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.

Scott’s formula (7.25) for χ α is<br />

where<br />

∫ t<br />

χ 2 c ln χ α = 4π<br />

0<br />

dt ′<br />

k 2 (t ′ )<br />

N e ∑<br />

i=1<br />

N i α 2 i<br />

[<br />

ln χ α1 +<br />

]<br />

ln χel i<br />

Z i<br />

. (2.306)<br />

k = p/¯h, (2.307)<br />

p = particle momentum,<br />

N i =<br />

N a ρ<br />

M p i = density of atoms of type i, (2.308)<br />

α i = αZ i /β = Z i e 2 /¯hν , (2.309)<br />

χ αi = the screening angle for atoms of type i,<br />

[<br />

= χ 2 0(1.13 + 3.76α 2 1/2<br />

i )]<br />

, (2.310)<br />

χ 0 =<br />

¯λ 0<br />

,<br />

r TF<br />

(2.311)<br />

¯λ 0 = λ 0 /2π = ¯h/p = wavelength of electron/2π, (2.312)<br />

r TF = Thomas − Fermi radius of atom<br />

= 0.885 a 0 Z −1/3<br />

i , (2.313)<br />

a 0 = Bohr radius = ¯h 2 /m e e 2 (2.314)<br />

χ el<br />

i = screening angle for the atomic<br />

atoms of type i.<br />

The next step is to let β = 1 in the α i that are in the χ αi , to delete the term with χ el<br />

i , and to let<br />

Zi 2 → Z i(Z i + ξ MS ). Recalling that p = Eβ/c, Equation 2.306 now becomes<br />

χ 2 c ln χ α = N [<br />

aρ ∑ Ne<br />

M 4πe4<br />

i=1<br />

p i Z i (Z i + ξ MS ) ln χ αi<br />

] ∫ t<br />

Since e 2 = r 0 m e c 2 , and using Equation 2.301, we obtain<br />

χ 2 c ln χ α = N [ Ne<br />

]<br />

aρ ∑<br />

M 4πr2 0 p i Z i (Z i + ξ MS ) ln χ αi<br />

i=1<br />

Dividing by Equation 2.300 and multiplying by 2, we get<br />

[ Ne<br />

ln χ 2 α = ∑<br />

i=1<br />

But using Equations 2.310 through 2.314,<br />

ln χ 2 0 = ln [<br />

ln χ 2 α i<br />

= ln χ 2 0 + ln<br />

¯h 2 m 2 ee 4<br />

p 2 (0.885) 2¯h 4 Z −2/3<br />

i<br />

p i Z i (Z i + ξ MS ) ln χ 2 α i<br />

]<br />

0<br />

dt ′<br />

E 2 β 4 . (2.315)<br />

m 2 t<br />

. (2.316)<br />

ĔMS 2 β4<br />

Z −1<br />

S<br />

. (2.317)<br />

[<br />

1.13 + 3.76(αZ i ) 2] , (2.318)<br />

] [<br />

]<br />

m 2<br />

= ln<br />

ee 4<br />

p 2¯h 2 − ln Z −2/3<br />

(0.885) 2 i , (2.319)<br />

86

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

Saved successfully!

Ooh no, something went wrong!