02.05.2014 Views

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

C<br />

u<br />

2.3. Incoherent (Compton) scattering 47<br />

15 00<br />

m e<br />

c J( p z<br />

)<br />

1000<br />

5 00<br />

Au<br />

Al<br />

0<br />

0.001 0.010 0.100<br />

p z<br />

/ m e<br />

c<br />

Figure 2.6: Atomic Compton profiles (p z > 0) for aluminium, copper and gold. The continuous<br />

curves are numerical Hartree-Fock profiles tabulated by Biggs et al. (1975). The dashed<br />

curves represent the analytical profiles defined by eq. (2.57). (Adapted from Brusa et al.,<br />

1996.)<br />

is assumed. In the Hartree-Fock approximation for closed-shell configurations, the momentum<br />

distribution of the electrons in an atomic shell, obtained by adding the contributions<br />

of the orbitals in that shell, is isotropic. For an isotropic distribution, expression<br />

(2.39) simplifies to<br />

J i (p z ) = 2π<br />

The atomic Compton profile is given by<br />

∫ ∞<br />

|p z|<br />

p ρ i (p) dp. (2.41)<br />

J(p z ) = ∑ i<br />

f i J i (p z ), (2.42)<br />

where f i is the number of electrons in the i-th shell and J i (p z ) is the one-electron profile<br />

of this shell. The functions J(p z ) and J i (p z ) are both bell-shaped and symmetrical about<br />

p z = 0 (see fig. 2.6). Extensive tables of Hartree-Fock Compton profiles for the elements<br />

have been published by Biggs et al. (1975). These numerical profiles are adequate for<br />

bound electron shells. In the case of conductors, the one-electron Compton profile for<br />

conduction electrons may be estimated by assuming that these form a free-electron gas<br />

with ρ e electrons per unit volume. The one-electron profile for this system is (see e.g.<br />

Cooper, 1971)<br />

J feg<br />

i (p z ) = 3 ( )<br />

1 − p2 z<br />

Θ(p<br />

4p F p 2 F − |p z |),<br />

F<br />

J feg<br />

i (0) = 3<br />

4p F<br />

, (2.43)<br />

where p F ≡ ¯h(3π 2 ρ e ) 1/3 is the Fermi momentum. For scattering in a compound material,

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

Saved successfully!

Ooh no, something went wrong!