Relativistic Dynamics of the Compton Diffusion on a Bound Electron
Abstract
A covariant relativistic formalism for the electronphoton and nuclear dynamics is summarised making more accurate predictions in agreement with experiments for Compton scattering in shells with large electron binding energy. An exact solution for the Dirac equation for an electron in the nuclear Coulomb field is obtained, in order to write the relativistic dynamics for this QED process. This is a preparation for the calculation of the relativistic crosssection for Compton scattering on bound electrons; as a precision test for QED.
pacs:
04.20.Jb , 03.65.Pm , 12.20.m , 13.60.FzI Introduction
It is my great pleasure to write this article in honour of Monsier Directuer de recherche Frnaçois Vazeille, under the request of my student L.A Al Asfar and some of my colleagues to summarise and review the mathematical and theoretical techniques I have used in my PhD thesis for full relativistic calculation and description of Compton scattering on bound electrons my (). This description agrees completely with experimental data provided at the time I had begun writing my thesis, and up until now. It could be considered as one of the precision tests for quantum electrodynamics (QED).
It has been almost 30 years since I presented my calculations, and the literature still lack an exact full relativistic description of Compton diffusion on bound electrons, it is apparent that the work in my PhD dissertation needs to be reformulated in modern language and summarised in English so the calculations and techniques become accessible to the scientific community.
Compton diffusion is an inelastic scattering process between a free, or bound electron and a photon, and it is one of most important QED processes between electrons and photons of energy less than 1 MeV. Since materialisation dominates the processes for photons of energy higher than 1 MeV. One may consider the electron target in Compton diffusion as ’free’ if its binding energy is much smaller than the incident photon’s energy . One then uses KleinNishima equation klein1929streuung () to deduce the wellknown Compton relation for the incident photon’s wavelength and the emergent one’s :
(1) 
with ,
(2) 
the Compton wavelength of the electron. However, the formula (1) is merely an approximation when bound electrons are considered , and fails to comply with experimental data when is comparable to . In particular, taking the ratio between the differential crosssections for Compton diffusion that was experimentally measured for the Kshell electron storm1970photon (); and the one calculated from Klien Nishima formula . We find that the ratio is varying between 0.25 and 1.175, depending on the diffusion angle. Similarly for the total crosssections ratio , that varies between 0.1 and 1. Showing a fundamental difference between free and bound electron treatment in Compton diffusion.
The first exact calculation of elastic scattering on bound electrons was carried by C. Fronsdal in 1966 veigele1966compton (), As for nonelastic Compton diffusion; the relativistic cross section were first approximated by gavM72 () starting from nonrelativistic study, then carrying out dipole approximation gavrila1972compton (). Later, inclusion of electron’s spin into the calculation were made by grotch1983spin () making a quasirelativistic calculation for the photon cross sections. More modern study made a full relativistic approximation. Other studies helped deepening the understanding of this important QED process like the studies carried by Ribbers ribberfors1975relationship () on the rôle of momentum and angular distributions, along with photon’s polarisation ribberfors1975relationshipII (). A more general study were published few months before my dissertation calculated the relativistic Compton cross sections, using HartreeFock central field functions PhysRevA.37.3706 ().In this paper, I shall prepare for the calculation of the relativistic cress section from the covariant matrix elements. I shall start by analysing Compton scattering on bound electrons, derive a generalised Compton formula by taking into an account the nuclear recoil. Them derive the wavefunctions of the bound electrons, by analytically solving the Dirac equation exactly . Then the relativistic dynamics of the electron, photon and the nucleus are discussed at the end of this paper.
Ii Covariant formalism
There are two Feynman diagrams for the Compton diffusion shown in figure 1 . The first diagram shows absorption first process; where the incident photon is absorbed at vertex 1 then remitted at vertex 2. The second demonstrates emission first process, where the emergent photon is first emitted at vertex 1, then the incident photon is absorbed at vertex 2. However, these two processes are indistinguishable experimentally, hence we write the propagator of the second order process as:
(3) 
Where,  
(4a)  
(4b) 
The propagator solves the equation:
(5) 
for the free electron, and for the bound electron case :
(6) 
where , the gauge covariant derivative, and is an exterior field that satisfies Lorentz gauge condition . Note that the propagator no longer depend on the distance between and , but on the points themselves.We wish to diagonalise and write it in the form :
(7) 
The matrix , are used to calculate the cross sections, as a substitute for KlienNishima formula. This sums the covariant formalism for the free and bound electrons in Compton Diffusion, but we shall only focus on the bound electron case from now on.
Iii Kinematic States
In this section, I shall derive a generalised kinematic formula for Compton diffusion, by considering the Coulomb field and the recoil of the nucleus. From this section and on, I shall consider a units system where . Hence momenta have the units of and energies are written in the units of .
iii.1 Notation
The electron’s initial and final total energies shall be respectively denoted by :
(8a)  
(8b) 
The term , corresponds to the electron’s initial kinetic energy , and its binding energy with the nucleus . The final energy term is only given by the final kinetic energy, as the final electron is considered to be free, viz .
As for the photon’s initial and final 3momenta, they are denoted by and , where the energy is given by and . And the 4momenta satisfy the ’nulllike’ condition :
since the initial and the final photon states are real. Finally, for the nucleus states, the initial state of the nucleus of a mass is considered to be static, thus the only energy it has is the one corresponding to its mass. Whilst in the final state the nucleus will acquire a recoil kinetic energy . The 4momenta for the nucleus are written as:
(9a)  
(9b) 
iii.2 The Generalised Formula
Since the 4momenta are discussed above, we can set up the diffusion in the laboratory system of reference (LAB) as figure 2 shows. The initial scattering occurs in the plane, the incident photon approaches from the positive direction. From the figure one can write the following geometric relation:
(10) 
Then one may determine the scattering angle of the electron with respect to the nucleus by ; that takes values , and the initial angle for the electron is , that takes the values . The angels and determine the position of the nucleus in the LAB system. they are related to other angels by the relations:
(11)  
(12) 
From these geometric relations, and basic conservation laws. A generalised formula for Compton diffusion can be written in the simple form :
(13) 
The factor is written as:
(14) 
The couple determine the approximation in mind, depending on the value of as follows:

K=0 . The Classical Compton formula compton (); for free and stationary electron.
(15) 
K=1, Dumond’s formula dumond1933linear () for free but moving electron. The Doppler effect correction is included (14).
(16) 
K=2, Correction for electron’s binding with the nucleus is obtained . Matching the results of Ross and Kirkpatrick ross1934constant ().
(17) 
K=3, Corrections for moving and bound electron are obtained; the recoil of the nucleus is taken into an account. Matching the results of Veigele et al. veigele1966compton ().
(18)
As we can see, the most general case is the one obtained from taking . In some cases the nuclear recoil plays an important rôle in the scattering process that cannot be neglected. Particularly when the recoil 3momentum reaches its maximal value; this can be deduced from the geometry of the diffusion:
(19) 
The last parameters of the formula (14) can be rewritten in this case as :
(20) 
The formula (14) can be used to calculate the final energies for the photons and electron:
(21) 
However, both formulae depend on the nuclear recoil. In order to overcome this dependence, one can consider the limiting cases.Considering the cases when the photon is completely absorbed (), or when it bounces off without loosing any energy to the electron i.e , and takes its maximal value. We obtain, a formula for the maximal/minimal values of :
(22) 
Where:
(23) 
and
(24) 
The formula above is a correction to Compton’s formula for bound electrons that allows the calculation for the final 4momenta for electron in the limiting cases specified above.
The full relativistic dynamics of electron resembled by its initial 4 momentum is however more difficult to compute, it requires the calculation of the wavefunctions for the electron using Dirac equation. This shall be the task for the next section.
Iv Bound Electron’s Relativistic Wavefunction
Adopting the notation in the previous section, Dirac equation for the bound electron takes the explicit form :
(25) 
Where ^{1}^{1}1The ’particle physicist convention is used, with as the time index label, and metric signature of ’mostly minuses’., is the set of Dirac gamma matrices, with spinor indices suppressed.
The Dirac spinor can written as in van der Warden representation:
(26) 
Thus (IV) reads eliminating the spinor  :
(27) 
Where are the set of Pauli spin matrices.
Clearly (27) is separable, the solution takes the form :
(28) 
The first function is the radial function , while the second resembles the angular part . The latter function contains the coupling between the orbital and spin angular momenta for the electron, it can be written as :
(29) 
where:

are the total and orbital angular momenta quantum numbers, respectively

is the eigenvalue for the zcomponent of the angular momentum .

, is Pauli spinor corresponding the the spin eigenvalues . Such that
and . 
, the normalized spherical harmonics

, the corresponding ClebschGordan coefficients with respect to the quantum numbers mentioned above.
One can express the wavefunction in (IV) as:
(30) 
This expression is deduced from the expression from the orbital angular momentum operator representation in the sphericalpolar coordinates:
(31) 
the number resembles the angular quantum number that shall be discussed below.
Now, the spinor can be recovered, and expressed in analogous manner to , the radial function is introduced corresponding to the spinor  similar to . The radial parts of Dirac equation then becomes :
(32)  
(33) 
The wavefunction now takes the form:
(34) 
Which will be useful for later calculations of the dynamics of Compton scattering.
Now we turn to define the angular quantum number . This quantum number can be defined as the eigenvalue of the operator with the angular functions , viz:
(35) 
Hence, depends of the three quantum numbers and :
(36) 
noting it is either larger than or less than . Taking both positive and negative values, but never zero.
It is then tempting to express the angular momenta quantum numbers in terms of guided by the expression for it in (36). Moreover, writing the eigenfunctions and , more conveniently, in terms of . That is
and so on. We can also define the action of the operator on as :
(37) 
Finally we are ready to write the final form of the wavefunction, noting that the solution only depends on the quantum numbers and ; thus rewritings eq (34) using the previous result from (37) as:
(38) 
iv.1 Radial wavefunctions
iv.1.1 Continuous wavefunctions
In order to calculate the radial wavefunctions that appear in (38), one needs to find the potential. In this study, the potential is the central Coulomb potential of the nuclei , with being the finestructure constant, and is the atomic number. The argument followed here is similar to rose1961relativistic (), but with the potential taken into an account.
Taking , the spectrum of the electron’s wavefunctions should be continuous. The radial wavefunctions therefore is expressed as :
(39) 
(40) 
Where, and are the hypergeometric and gamma functions respectively. The factors appearing the expressions above are defined as the following:
(41) 
The factor is a normalisation factor and it is found to be :
(42) 
Defining the factor in similar fashion to overbo1968exact ():
(43) 
thus , depending on and .
iv.1.2 Discrete wavefunctions
For electron’s energy , the radial wavefunctions should be discrete. The take the following general form (following the lead of rose1961relativistic ()) :
(44) 
(45) 
With , being the electron’s relativistic binding energy, , with being the principal quantum number, and . The solutions and , contain the hypergeometric function with negative ’entries’, hence they can be related to Laguerre polynomials . One may express the solution in terms of Laguerre polynomials and Gamma function, in order to make numerical calculations more efficient, first defining :
(46) 
Then the radial wavefunctions are written as :
(47) 
(48) 
The ratio , satisfies the limit as , and :
(49) 
One can find a ’simple’ expression for the discrete radial functions, for taking values from to . Notice for shells with these principle quantum numbers, we have states per shell, determined fully by and . The simple form of the wavefucntions is expressed as the following:
(50a)  
(50b) 
Where and , are polynomials of of the th order :
(51) 
with
(52) 
Moreover, for , we can use the special functions representation of and to find the m th coefficients and :
(53) 
(54) 
We may also define :
(55)  
Notice that from this expansion of the radial wavefunctions, it is clear that they depend on and , but not on . The formulae above are used to calculate the radial wavefunction for the state for in the figure 3, and for state for the same element figure 4 , as an example. We observe that the radial wavefunction for the tends to be smaller and smaller as the shell gets farther and farther from the nuclear Coulomb field. On the converse, it becomes considerable for shells closer to the nucleus.
We are ready now to write the total wavefunction for the Dirac field around the nucleus. This wavefunction will be important to describe the scattering process when we shall discuss the propagators. We shall used the formula (29) to describe the angular and spin part of the total wavefunction, the formulas (44) and (45) or their ’simple’ form (50b) to describe the radial part. Noting that for Compton diffusion for bound electrons in the shells close to the nucleus, we ought to consider both radial wavefunctions and , as both of them play a rôle as we saw.
V Dynamic Electron States
Since we have the wavefunction for the bound electron , it is possible to study the dynamics of the electron in parallel to what was established in the kinematics section. The full dynamics of the electron  before the diffusion process is described by its initial momentum and its binding energy . A full relativistic calculation for these is possible since the relativistic wavefunction is found.
v.1 Relativistic binding energy
The relativistic binding energy for the electron is given by the formula rose1961relativistic () :
(56) 
keeping the same parameters from previous section.
v.2 Relativistic momentum
In order to calculate the momentum of the electron, we need fist to calculate the expectation value of its square, then take the square root of that result, viz:
The subscript in shall be dropped throughout the calculations of this section, as this argument is valid for any dynamical state for the electron. The value is calculated from:
(57) 
The kinetic energy is given by rose1961relativistic ():
(58) 
In fact, the product yields the operator . Thus substituting (58) in (57), and by the orthonormal property of the angular wavefunctions in (29), we have:
(59) 
Where :
(60) 
is the position representation for the operator obtained from (58). To calculate the integrals in (59) we use the formula for the radial wavefunctions in (50b), we obtain after lengthy but straightforward calculation :
(61) 
With:
(62) 
for , and :
(63)  
(64)  
(65) 
(66)  
(67) 
(68)  
We notice from the expression(61) that the momentum depends on the ’effective’ atomic number , with being the charge density of the other electrons. In other words, one cannot ignore the interaction between the electron of study and other electrons in the other shells gavM72 (). We can use the previous formulae to make numerical calculations for and and compare them with the nonrelativistic ones as in tables 1 and 2 for the element (lead).
Shell  K  L  L  L 

( non relativistic)  
(relativistic) 
Shell  K  L  L  L 

( non relativistic)  
(relativistic) 
From the initial binding energy and momentum, one can use the formula 14 to know the final state of the electron after the diffusion.
Vi The Nuclear Dynamics
The maximal nuclear recoil is when the final electron is state is a free electron, viz the electron is liberated from the atom due to the diffusion process. We are interested in calculating the nucleus recoil in this case. The nuclear momentum and kinetic energy are related by conservation of 4momenta by:
(69) 
Using the equation above and equation (20) we arrive to the equation that must satisfy :
(70) 
Since the initial states for the electron can determined from the previous section, and the photon’s final energy is a measurable quantity, the only quantity left is to be that can be determined a priori within the allowed ranges. The figure (5) demonstrates the nuclear recoil momentum for Compton diffusion on Kshell electron in element (Germanium) with incident photon’s energy keV
Since depends on the kinematics of the electron and the photon described by the generalised Compton formula (14), and at which degree of approximation described by .
Vii Concluding Remarks
A generalised Compton formula is derived for moving bound electrons, from a geometrical argument. This formula allows the calculation for the final electron’s energies and momenta, and predicts the final photon energy and nucleus recoil provided the incident photon’s energy is known. The covariant formalism for the diffusion process in quantum electrodynamics is discussed for bound electrons, where the binding energy from the Coulomb potential is no longer negligible.This, along with finding the electron’s relativistic wavefunction prepares for the calculation of the full relativistic cross sections, that are supposed to be more aligned with experimental results than previous calculations. The relativistic wavefunctions are obtained from analytical solution of the Dirac equation for electron in (strong) field. It is evident from this solution that for internal shells (particularly and ), the target electron possess ’large’ and ’small’ radial wavefunctions corresponding to its interaction with the virtual particles created from the (strong) electromagnetic field present in the vicinity of the nucleus. More effects are shown to take part in the diffusion process for the target electron,from the study of its dynamics. For example its interaction with other electrons in the other shells, and the nuclear recoil after the scattering process.
This paper provided a detailed discussion and calculation for this important QED process, and prepares for further exact calculation of a (generalised) relativistic cross section that are made in other paper.
Acknowledgement
This research project was supported by a grant from the ’ Research Center of the Female Scientific and Medical Colleges ’, Deanship of Scientific Research, King Saud University.
References
 (1) S. Al SalehMahrousseh. Calcul relativiste en electrodynamique quantique de la diffusion compton sur un electron lie. Theses, Université Blaise Pascal  ClermontFerrand II, 1988.
 (2) Jesse WM DuMond. The linear momenta of electrons in atoms and in solid bodies as revealed by Xray scattering. Reviews of Modern Physics, 5(1):1, 1933.
 (3) C. Fronsdal. Compton Scattering from Bound Electrons. Phys. Rev., 179:1513–1517, Mar 1969.
 (4) Mihai Gavrila. Compton Scattering by Shell Electrons. I. Nonrelativistic Theory with Retardation. Phys. Rev. A, 6:1348–1359, Oct 1972.
 (5) Mihai Gavrila. Compton scattering by Kshell electrons. II. Nonrelativistic dipole approximation. Physical Review A, 6(4):1360, 1972.
 (6) H Grotch, E Kazes, G Bhatt, and DA Owen. Spindependent Compton scattering from bound electrons: Quasirelativistic case. Physical Review A, 27(1):243, 1983.
 (7) Peter Holm. Relativistic Compton cross section for general centralfield HartreeFock wave functions. Phys. Rev. A, 37:3706–3719, May 1988.
 (8) Oskar Klein and Yoshio Nishina. Über die Streuung von Strahlung durch freie Elektronen nach der neuen relativistischen Quantendynamik von Dirac. Zeitschrift für Physik, 52(1112):853–868, 1929.
 (9) The ’particle physicist convention is used, with as the time index label, and metric signature of ’mostly minuses’.
 (10) Ingjald Øverbø, Kjell J Mork, and Haakon A Olsen. Exact calculation of pair production. Physical Review, 175(5):1978, 1968.
 (11) Roland Ribberfors. Relationship of the relativistic Compton cross section to the momentum distribution of bound electron states. II. Effects of anisotropy and polarization. Physical Review B, 12(8):3136, 1975.
 (12) Roland Ribbers. Relationship of the relativistic Compton cross section to the momentum distribution of bound electron states. Physical Review B, 12(6):2067, 1975.
 (13) Morris Edgar Rose and WH Furry. Relativistic electron theory. American Journal of Physics, 29(12):866–866, 1961.
 (14) PA Ross and Paul Kirkpatrick. The Constant in the Compton Equation. Physical Review, 45(3):223, 1934.
 (15) Lellery Storm and Harvey I Israel. Photon cross sections from 1 keV to 100 MeV for elements Z= 1 to Z= 100. Atomic Data and Nuclear Data Tables, 7(6):565–681, 1970.
 (16) William J Veigele, Philip T Tracy, and E Michael Henry. Compton effect and electron binding. American Journal of Physics, 34(12):1116–1121, 1966.