* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Relativistic Description of Two-body Scattering
Scalar field theory wikipedia , lookup
Bohr–Einstein debates wikipedia , lookup
Schrödinger equation wikipedia , lookup
Renormalization group wikipedia , lookup
Wave function wikipedia , lookup
Path integral formulation wikipedia , lookup
Two-body Dirac equations wikipedia , lookup
History of quantum field theory wikipedia , lookup
Canonical quantization wikipedia , lookup
Geiger–Marsden experiment wikipedia , lookup
Molecular Hamiltonian wikipedia , lookup
Double-slit experiment wikipedia , lookup
Feynman diagram wikipedia , lookup
Renormalization wikipedia , lookup
Particle in a box wikipedia , lookup
Identical particles wikipedia , lookup
Wave–particle duality wikipedia , lookup
Hydrogen atom wikipedia , lookup
Matter wave wikipedia , lookup
Dirac equation wikipedia , lookup
Atomic theory wikipedia , lookup
Symmetry in quantum mechanics wikipedia , lookup
Quantum electrodynamics wikipedia , lookup
Elementary particle wikipedia , lookup
Theoretical and experimental justification for the Schrödinger equation wikipedia , lookup
Relativistic Description of Two-body Scattering Reactions Gashaw Adera ([email protected]) African Institute for Mathematical Sciences (AIMS) Supervised by Dr. Brandon Van Der Ventel University of Stellenbosch June 4, 2007 Abstract In this paper, two-body elastic scattering is treated in the relativistic quantum mechanics framework. By using lowest-order Feynman diagrams, detailed derivations of the invariant matrix element and hence the differential cross-section for unpolarized electron-proton scattering are made. First, the proton is approximated as a spin- 12 point particle which allows a consistent quantum electrodynamical description of the scattering process. The Feynman rules and trace algebra have been employed in constructing the relativistic quantum mechanical expression of invariant amplitude. Moreover, by using a suitable basis for the second rank tensor, the hadronic tensor for point proton is generalised to include the electromagnetic form factors which lead us to treat the proton as an extended object within the finite volume. The calculated differential cross section for a point proton is compared to the Rutherford and Mott predictions at laboratory angles between 0o and 180o and initial electron energy between 1 MeV and 1 GeV by using numerical simulations. The results are plotted against scattering angle in the laboratory frame. In the simulation it is shown that the calculated differential cross-section agrees with the Mott prediction. It also shows the expected deviation from the Rutherford prediction. i Contents Abstract i List of Figures iv 1 Introduction 1 2 Kinematics of Scattering Processes 2 2.1 2.2 2.3 Kinematics of Non-relativistic Scattering . . . . . . . . . . . . . . . . . . . . . 2 2.1.1 Galilean Transformation . . . . . . . . . . . . . . . . . . . . . . . . . . 2 The Kinetic Energies in Laboratory and Center-of-mass Frames . . . . . . . . . 4 2.2.1 Scattering Cross-section in Laboratory and Center-of-mass Frames . . . . 6 Kinematics of a Relativistic Particle . . . . . . . . . . . . . . . . . . . . . . . . 7 2.3.1 Energy-momentum Four-vector . . . . . . . . . . . . . . . . . . . . . . 7 2.3.2 Lorentz Transformation (Boost) . . . . . . . . . . . . . . . . . . . . . . 8 2.3.3 Mandelstam Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 3 Relativistic Dirac Equation 12 3.1 The Klein-Gordon (KG) Equation . . . . . . . . . . . . . . . . . . . . . . . . . 12 3.2 Dirac Free Particle Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 3.3 Free Particle Solution of Dirac Equation . . . . . . . . . . . . . . . . . . . . . 14 Normalisation of Dirac Spinor . . . . . . . . . . . . . . . . . . . . . . . 15 3.3.1 4 Quantum Electrodynamics of Electron-proton Elastic Scattering 4.1 16 Feynman Diagrams . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 4.1.1 Components of Feynman Diagram . . . . . . . . . . . . . . . . . . . . . 16 4.2 Feynman Rules for Tree Diagram . . . . . . . . . . . . . . . . . . . . . . . . . 17 4.3 Fermi’s Golden Rule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 4.3.1 Transition Rate . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 4.3.2 Cross-Section for Two-body Scattering . . . . . . . . . . . . . . . . . . 19 Elastic Scattering of Electron from a Point Proton . . . . . . . . . . . . . . . . 19 4.4 ii 4.4.1 4.5 Evaluation of Invariant Amplitude of Unpolarised Electron-proton Scattering 20 Proton Form Factors and Electron-proton Elastic Scattering . . . . . . . . . . . 5 Simulation, Discussion and Conclusion 24 28 5.1 Simulation and Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 5.2 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 A Background Mathematics 31 A.1 The Dirac Delta Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 A.2 Relativistic Notation and Four-vector Formalism . . . . . . . . . . . . . . . . . 32 A.3 Dirac Algebra of Gamma Matrix and Trace Theorems . . . . . . . . . . . . . . 33 A.3.1 Pauli Matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 A.3.2 Dirac Matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 A.3.3 Trace Theorems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 B The Python Code 36 Bibliography 40 iii List of Figures 2.1 The kinematics of two-body elastic collision A + B −→ C + D process in (b) as viewed from the lab system; (b) as viewed from the C.M. frame. . . . . . . . . . 2 Illustration of the geometrical representation vL and v with directions (θL , φL) and (θ, φ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 A two-body elastic collision A + B → C + D (a) in the lab frame (b) in the centre-of-mass frame . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.4 Illustration of scattering from a target into a solid angle element dΩL = sin θL dθL dφL 6 4.1 (a) Wavy line representing the outgoing photons and (b) straight lines representing the incoming fermion. If we change the position of the dots to the other ends, we get lines that represent (a) the incoming photon and (b) the outgoing fermion. We use dots to represent lepton-photon vertices . . . . . . . . . . . . . . . . . 16 4.2 Primitive diagram that illustrates the emission of photon from electron . . . . . 17 4.3 The lowest-order tree diagram . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 4.4 The lowest-order Feynman diagram for electron-proton elastic scattering. Each factor of the invariant amplitude are associated to the corresponding element of the Feynman diagram. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 Comparison of the calculated cross-section with (a) Rutherford differential crosssection and (b) Mott differential cross-section for 10 MeV, 50 MeV, and 90 MeV incident energies of electron. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 Comparison of the calculated cross-section with (a) Rutherford differential crosssection and (b) Mott differential cross-section for 100 MeV, 150 MeV, 350 MeV and 700 MeV incident energies of electron. . . . . . . . . . . . . . . . . . . . . 29 (a) Angular distribution of differential cross-section versus energy of the initial electron. (b) Simulation for energy the scattered electron versus the scattering angle at incident energies 100 MeV, 150 MeV, 200 MeV, and 250 MeV . . . . . 29 2.2 2.3 5.1 5.2 5.3 iv 1. Introduction The study of elementary particles and their interaction has brought about a better understanding of our universe from sub-nuclear scales to cosmological scales. Scattering experiments has made possible the discovery of fundamental constituents of matter and the kind of interactions they exhibit. In other words, scattering experiments form the basis of the study of the structure of particles. There are four fundamental interactions observed in nature: gravitational, electromagnetic, weak and strong nuclear forces. The electromagnetic interaction is mediated by the exchange of virtual photons. The weak nuclear interaction is responsible for beta decay and the exchanged particles are W ± and Z 0 . For the strong nuclear interaction between quarks, the exchanged particles are gluons. The theories that describe electromagnetic, weak nuclear and strong nuclear interactions are Quantum Electrodynamics (QED), the theory of Glashow, Weinberg and Salam (GWS theory), and Quantum Chromodynamics (QCD), respectively. It can be seen that fundamental particles have been classified based on their response to these interactions as well as their structure. The two main groups are hadrons and leptons. Hadrons experience all four the fundamental forces and they are composite particles. Leptons do not experience the strong force, and they are elementary particles, for example the electrons. The particles which mediate the different forces are generally called gauge bosons [Ryd96]. Scattering is a process in which two or more particles from a distant past momentarily come together and interact by exchanging energy and momentum via the field quanta and then separate from one another to a distant future. However, this essay is aimed at describing the scattering process of two fermions, namely, the electron and the proton, which interact electromagnetically by exchanging a single virtual photon, which has a non-zero q 2 [Ron94]. The particles scattered off the target emerge with a new state which contains important information regarding the target particle. We in turn use the information to investigate the structure and properties of the target particle from which the incident particle is scattered. The quantity that is used to characterise the scattering process is called the scattering cross-section, which is the measure of the tendency for the incoming particle to scatter in certain direction after the interaction. The main focus of this study is the derivation of the invariant amplitude and the differential cross-section of unpolarized electron-proton elastic scattering by using Feynman diagrams. A study of this kind is of paramount importance in comparing the contemporary theories of Nuclear and Particle Physics with scattering experiments. In the second chapter non-relativistic and relativistic kinematics of elastic scattering process are discussed by using the collinear reference frames. Since e-p scattering is due to electromagnetic interaction, our study is made using the framework of Quantum Electrodynamics which describes elementary particles as the quantum states of electromagnetic field. The Klein-Gordon and Dirac equations for relativistic particles are also briefly presented in Chapter 3. In the fourth chapter more attention is given to the derivations of the scattering cross-section and the invariant amplitude. Firstly we consider the proton to be a structureless fermion. As the next step we take into account the proton structure by means of the form factor description. Simulations and interpretation of our results as well as inference of the study are contained in the fifth chapter. 1 2. Kinematics of Scattering Processes 2.1 Kinematics of Non-relativistic Scattering The kinematics of a system of particles are derived based on energy and momentum conservation laws [BJ95]. The derivation of the differential cross-section requires kinematic descriptions. Thus we set up the kinematic framework by first defining the frames of reference that are essential in the study of scattering in general, in particular elastic scattering in which the total energy and momentum are conserved. To begin with, we consider the elastic collision of two particles A and B having masses mA and mB . The laboratory (lab) frame of reference is the coordinate system in which the target particle B is initially at rest. The centre-of-mass (C.M.) frame of reference is a system in which the centre of mass is always at rest before and after collision [Joa84]. In the non-relativistic regime, we relate energies, scattering angles and cross-sections measured in lab frame to their C.M. counterparts using a Galilean transformation, GT . z C (pA )L A pC z (pC )L C A B pA D D pB B pD (pD )L (a) (b) Figure 2.1: The kinematics of two-body elastic collision A + B −→ C + D process in (b) as viewed from the lab system; (b) as viewed from the C.M. frame. 2.1.1 Galilean Transformation For a system of two non-relativistic particles A and B, the space-time four-vector coordinates with respect to the lab and the C.M. frames can be given as (tL , xL , yL , zL ) and (t, x, y, z), respectively. For convenience, we choose the z-axis as the incident direction of particle A (see figure 2.1). A centre of mass of the two-body system always retains its uniform motion relative to any inertial frame. If we assume that VL = VL ẑ be the velocity of the centre of mass with respect to the lab frame, then the explicit equations of the Galilean transformation from the lab frame to the C.M. frame can be written as t = tL ; x = xL ; y = yL ; 2 z = zL − VL tL (2.1) Section 2.1. Kinematics of Non-relativistic Scattering Page 3 z θL θ VL vL v y φL = φ x Figure 2.2: Illustration of the geometrical representation vL and v with directions (θL , φL) and (θ, φ) Similarly, we can go from the C.M. frame to the lab frame by only changing −VL to +VL tL = t; xL = x; yL = y; zL = z + VL t (2.2) Here, we can also relate the velocity in the lab frame to its counterpart in the C.M. frame. If we denote the velocities of the incident particle in the collision process by vL and v with the corresponding directions specified by (θL , φL ) and (θ, φ) in the lab and C.M. frames, respectively, they are related as vL = v + VL (2.3) So by simply referring to figure 2.2 and applying vectorial analysis, one can get the relation between the lab frame angle θL and the C.M. frame angle θ. That is, tan θL = or cos θL = sin θ (cos θ + τ ) sin θ 2 (sin θ + (cos θ + τ )2 ) where τ= 1 2 VL |v| = (2.4) sin θ 1 (1 + 2τ cos θ + τ 2 ) 2 (2.5) (2.6) Section 2.2. The Kinetic Energies in Laboratory and Center-of-mass Frames 2.2 Page 4 The Kinetic Energies in Laboratory and Center-of-mass Frames The total momentum of the system of, say, two particles A and B in the lab frame is equal to the velocity VL of the centre of mass of the system times the total mass M = mA + mB , that is PL = MVL (2.7) The initial kinetic energy in lab system is given by TiL = P2 (pA )2L = L 2mA 2mA (2.8) where (pA )L is the momentum of the incident particle A in the lab frame. We can establish the relation between the kinetic energies in the lab and the C.M. systems by using the relative momentum, which is invariant under the Galilean transformation, and the reduced mass of the system of colliding particles [Joa84]. If we consider a system of two particles 1 and 2, the relative momentum p in any inertial frame can be defined as m2 p1 − m1 p2 (2.9) p= m1 + m2 and the reduced mass µ is defined as µ= m1 m2 m1 + m2 (2.10) where m1 and m2 are the respective masses and p1 and p2 the corresponding momenta of the two particles. Thus in the C.M. frame pi is given by pi = pA = −pB (2.11) Then the initial kinetic energy Ti in the C.M. frame is given by Ti = p2 p2 p2A + B = i 2mA 2mB 2µi (2.12) where pA and pB , respectively, are the three-momenta of particles A and B, which are situated in the initial channel, with µi as their reduced mass. The relative momentum pi is pi = mB (pA )L mA + mB (2.13) Therefore, by plugging Eq.(2.13) into Eq.(2.12) and recalling Eq.(2.8), the relation between kinetic energy Ti in the C.M. system and its counterpart (Ti )L in the lab system of the initial channel for the collision process of particles A and B can be set up as Ti = mB (Ti )L mA + mB (2.14) Section 2.2. The Kinetic Energies in Laboratory and Center-of-mass Frames Page 5 We can get the same expression for Ti by taking out the kinetic energy of the centre of mass from the total kinetic energy available in lab system for its motion is irrelevant in the description of scattering process [BJ95]. That is, P2L 2M where M is the total mass mA + mB of the two-body system. (2.15) Ti = (Ti )L − Based on figure 2.3 which is the illustration of the binary rearrangement collision, A+B → C +D, we can also set up the relation between the final kinetic energies (Tf )L and Tf in the two frames for particles C and D of the final channel. That is (Tf )L = (p2C )L (p2D )L + 2mC 2mD (2.16) and the total momentum in the lab system is PL = (pC )L + (pD )L (2.17) where (pC )L and (pD )L , respectively, are momenta of particle C and D in the lab frame. Here again, the GT invariant relative momentum pf of the final channel is defined in the same way as Eq.(2.11), that is pf = pC = −pD (2.18) and pf = ΩL ≡ (θL , φL ) (pC )L C (pA )L A mD (pC )L − mC (pD )L mC + mD D C A pA φL D (pD )L (a) Ω ≡ (θ, φ) pC z B (2.19) pB B z φ pD (b) Figure 2.3: A two-body elastic collision A + B → C + D (a) in the lab frame (b) in the centre-of-mass frame Thus the final kinetic energy in the C.M. system is p2f p2D p2C + = Tf = 2mC 2mD 2µf (2.20) Section 2.2. The Kinetic Energies in Laboratory and Center-of-mass Frames Page 6 where µf the reduced mass of particles C and D. By applying the same definition as that of Eq.(2.15), we can relate the final channel kinetic energies available in the lab and the C.M. systems as Tf = (Tf )L − 2.2.1 (P2L) 2M (2.21) Scattering Cross-section in Laboratory and Center-of-mass Frames Before we proceed to the kinematics of cross-section, its better to have an idea of what crosssection refers to. Cross-section is a quantity that provides the information regarding the interaction in the scattering process. That is, based on the available initial and final states of colliding particles, it tells the probability for a particular scattering to take place. Note that this does not refer to the geometry of the target particle, rather it refers to the effective area over which two or more particles interact such that they make a transition from an initial state to a final state. So it represents the intrinsic scattering probability of the particular process [HM84]. Here we define a differential cross-section dσ in laboratory frame in terms of number of particles appearing in the initial and final states of the scattering process. If we consider a particle scattered off a target, which is placed at the origin of the lab frame, and then passes through small surface element subtending a solid angle element dΩL = sin θL dθL dφL at a polar angle θL and azimuthal angle φL with respect to the incident direction (see figure (2.4), then the differential cross-section is Ns (θL , φL )dΩ dσ = (2.22) Nb · Nt where Nb is incident flux, Nt is number of particles placed at the target, and Ns (θL , φL) is the flux of the scattered particle through a unit solid angle dΩ about the direction (θL , φL). xL dΩL φL incident beam dφL dθL θL zL target yL Figure 2.4: Illustration of scattering from a target into a solid angle element dΩL = sin θL dθL dφL Section 2.3. Kinematics of a Relativistic Particle Page 7 In Eq.(2.22) we can see that the same number of particles scattered into solid angle dΩ about the direction (θ, φ) in the C.M. frame as are scattered into dΩL about (θL , φL ) in the lab frame [Joa84]. So for the scattered particle C in the reaction A + B → C + D we have that, or dσC dσC (θL , φL )dΩL = (θ, φ)dΩ dΩL dΩ (2.23) dσC dσC (θL , φL ) sin θL dθL dφL = (θ, φ) sin θdθdφ dΩL dΩ (2.24) Scattering cross-section has dimensions of area-cm2 . More conveniently we use “barns”. 1 barn = 10−24 cm2 ; differential cross-sections dσ/dΩ are given in barns per steradian [Gri87]. If we refer back to figure (2.2) we notice that dφL = dφ, then for the scattered particle C Eq.(2.5) and (2.6) can be rewritten as cos θ + τ (2.25) cos θL = 1 (1 + 2τC cos θ + τC2 ) 2 with τC = VL | vC | (2.26) Then by substituting Eq.(2.25) into Eq.(2.24) we get 3 dσC (1 + τC2 + 2τC cos θ) 2 dσC (θL , φL) = (θ, φ) dΩL |1 + τC cos θ| dΩ (2.27) The same is true for particle D which is also scattered in the direction of (π − θ, π + φ) in the C.M. frame. 2.3 2.3.1 Kinematics of a Relativistic Particle Energy-momentum Four-vector From now on we use relativistic notations and four-vector formalism briefly discussed in appendix A.2. For any relativistic massive particle its energy E and three-momentum p are related as E 2 = m2 c4 + p 2 c2 (2.28) where m is rest-mass of the particle and c is the speed of light. The relativistic three-momentum and energy for a particle with 3-vector velocity are given by p = mγv (2.29) and E = mγc2 or E = mγc c (2.30) Section 2.3. Kinematics of a Relativistic Particle Page 8 where γ is Lorentz factor which we define later. Thus the energy-momentum four-vector pµ of the particle is defined as E pµ = ( , p) c (2.31) which gives E2 − p · p = m2 c2 (2.32) c2 where g µν is the metric tensor with zero off-diagonal and (+1, −1, −1, −1) diagonal elements (see appendix A.2). In natural units p2 = p · p = pµ gµν pν = pµ pµ = p2 = E 2 − p 2 = m2 (2.33) Besides, the relativistic energy of the particle is related to its kinetic energy T by E =T +m= p2 +m 2m (2.34) We often use the notation either p · x or pµ xµ to represent the scalar product or the contraction, p · x = pµ xµ = Et − p · x (2.35) Note also that we use either p or pµ interchangeably to represent the four-vector momentum. 2.3.2 Lorentz Transformation (Boost) In Minkowski space, in which special relativity is well defined, we have to take into account the distortions of the longitudinal space and the time components of the four-vectors, such as space-time coordinates of a point or an event xµ = (t, x) and energy-momentum four-vector pµ = (E, p), in the transformation from one inertial frame to another. For the study of scattering processes, the laboratory frame and the centre-of-mass frame are the two collinear frames we choose to transform kinematic and dynamic quantities. If we consider the C.M. frame which is moving with uniform velocity VL = VL ẑ along z-axis with respect to the lab frame, then we use the Lorentz transformation (LT ) to relate the coordinates of four-vectors such as, xµL = (t, xL ) and pµL = (EL , pL ) in the lab frame to their corresponding four-vectors xµ = (t, x) and pµ = (E, p) in the C.M. frame. In matrix form, the Lorentz transformation of space-time coordinate of a given event is written as tL γ 0 0 −γβ t x 0 1 0 0 xL = (2.36) y 0 0 1 0 yL zL −γβ 0 0 γ z or xµ = Λµν xνL (2.37) Section 2.3. Kinematics of a Relativistic Particle Page 9 where Λµν is the 4 × 4 matrix in Eq.(2.36). The corresponding Lorentz transformation of the space-time coordinates from the C.M. frame to the lab frame is xLν = Λµν xµ VL 2 where γ = 1 − 2 c in natural units. − 21 is called Lorentz factor and β = (2.38) VL . Notice that we are still working c The energy-momentum four-vector also transforms in the same manner. That is, pµ = Λµν pνL and pL ν = Λµν pµ (2.39) If we consider a particle that moves along z-axis, the second expression of Eq.(2.39) can explicitly be written as E γ 0 0 γβ EL 0 0 1 0 0 0 (2.40) 0 = 0 0 1 0 0 pz γβ 0 0 γ pz L But the total three-vector momentum p in the C.M. system is always zero; and this implies that pz = 0. From this argument we get that EL = γE and pL z = γβE from which we obtain pL VL = z (2.41) EL For the two-particle scattering case mentioned earlier, the VL becomes pL z (2.42) VL = EL + m where m is the mass of the target particle initially at rest in the lab frame. In what follows, we see how the three-vector velocity transforms. For convenience, we constrain the particle to xz-plane. To find the velocity that are related in terms of Lorentz transformation we need invariant time with respect to which we differentiate the coordinates of the space-time of the particle. We use the proper time recorded in the rest-mass frame of the particle to differentiate the coordinates of space-time. Then we get vz + VL vx (vz )L = and (vx )L = (2.43) (1 + VL vz ) γ(1 + VL vz ) By using Eq.(2.43) and referring back to figure(2.2), we get sin θ tan θL = (2.44) γ(cos θ + τ ) where VL (2.45) τ= |v| Relativistic quantum mechanics is in favour of quantities that are invariant under Lorentz transformation. Scalar product of a contravariant four-vector with its covariant counterpart is Lorentz invariant. Crucial examples are, xµ xµ = (t2 − x · x) pµ pµ = (E 2 − p · p) (2.46) (2.47) Section 2.3. Kinematics of a Relativistic Particle 2.3.3 Page 10 Mandelstam Variables We introduce three Lorentz invariant variables s, t, and u which are used to describe the reaction kinematics of two-body scattering; they are defined as, s = (pA + pB )2 = (pC + pD )2 t = (pA − pC )2 = (pB − pD )2 u = (pA − pD )2 = (pB − pC )2 (2.48) (2.49) (2.50) Their dependence on one another is constrained by the relation s + t + u = m2A + m2B + m2C + m2D (2.51) which can be proved by using the mass-shell condition p2A = m2A , p2B = m2B , p2C = m2C , and p2D = m2D . The variable s is equivalent to the square of energy available in the C.M. system and t is the square of the four-momentum transfer. In the C.M. frame we can redefine the initial and the final relative three-momenta pi and pf in terms of s which is expressed in the C.M. frame as 2 q q 2 2 2 2 2 (mA + |pi | ) + (mB + |pi | ) , (2.52) s = (EA + EB ) = that is q 1 |pi | = √ (s − m2A − m2B )2 − 4m2A m2B = 2 s q 1 (s − m2C − m2D )2 − 4m2C m2D = |pf | = √ 2 s q 1 √ (λ(s, m2A , m2B )) 2 s q 1 √ (λ(s, m2C , m2D )) 2 s (2.53) (2.54) where Eq.(2.54) is obtained from s = (EC +ED )2 and we have introduced the “triangle function”, λ(x, y, z) = (x − y − z)2 − 4yz. (2.55) Thus the energies can be expressed as 1 EA = √ (s + m2A − m2B ), 2 s 1 EB = √ (s − m2A + m2B ) 2 s (2.56) 1 EC = √ (s + m2C − m2D ), 2 s 1 ED = √ (s − m2C + m2D ) 2 s (2.57) and We can also get expressions for the other two Mandelstam variables t and u t = (pA − pC )2 = p2A + p2C − 2pA · pC = p2A + p2C − 2(E EC − |pA ||pC | cos θ) A θ = −4|p|2 sin 2 2 (2.58) Section 2.3. Kinematics of a Relativistic Particle and in a similar manner Page 11 θ u = −4|p| cos 2 2 2 (2.59) where we have set mA = mC , mB = mD , and |pf | = |pi | = |p|. Finally, the differential cross-sections for relativistic collision transform between the collinear frames in the same way as we have seen for the non-relativistic case in section 2.2.1. But we have to take into account the relativistic correction. Then for particle C in the final channel, the transformation relation is given by 2 3 γ (cos θ + τC )2 + sin2 θC 2 dσC dσC (θL , φL ) = (θ, φ) dΩL γ |1 + τC cos θ| dΩ (2.60) − 12 VL2 . The same holds for particle D that appears in the final channel. where γ = 1 − 2 c 3. Relativistic Dirac Equation In non-relativistic quantum mechanics, the Schrödinger equation is the fundamental differential equation that governs the behavior of quantum particles. For relativistic quantum mechanics the Klein-Gordon and Dirac equations describe non-interacting spin-0 and spin- 21 particles, respectively. However, when pair-production becomes kinematically allowed, then a simple one-particle differential equation is no longer valid, and we have to use quantum field theory as the theoretical framework [Gri87]. 3.1 The Klein-Gordon (KG) Equation Beginning with the relativistic energy and three-momentum relation E 2 = m2 c4 + p 2 c2 we can ∂ obtain the KG equation by using the corresponding quantum mechanical hermitian operators i ∂t ~ of E and p, respectively. By introducing the wavefunction ψ the Klein-Gordon equation and −i∇ is written as ∂2 ~ 2 − m2 )ψ ψ = (∇ (3.1) ∂t2 where m is the rest-mass energy of the particle in natural units. In relativistic notation Eq.(3.1) is rewritten as (pµ pµ − m2 )ψ = 0 (3.2) ∂ ~ ∂ ~ [Ryd96], hence where pµ = (E, −p) = i∂µ = i( , ∇) and pµ = (E, p) = i∂ µ = i( , −∇) ∂t ∂t pµ pµ = −( ∂2 ~ 2) −∇ ∂t2 (3.3) The continuity equation for the KG equation can be obtained from Eq.(3.1) and written as ∂ρ ~ ~ +∇·J =0 ∂t where the probability density ρ and the probability current density J~ are: ∂ψ ∗ ∗ ∂ψ ρ=i ψ − ψ ∂t ∂t and ∗~ ∗ ~ ~ J = −i ψ ∇ψ − ∇ψ φ (3.4) (3.5) (3.6) where ψ ∗ is the complex conjugate of ψ. Since the KG equation is second order in t, it has got positive and negative energy solutions and corresponding positive and negative probability densities. However, the negative probability density does not make sense for it violates the laws of statistical mechanics. This means that it is always expected to be scalar quantity. Dirac successfully derived the relativistic wave equation which is first order in both t and x and with positive definite probability density [BJ95]. 12 Section 3.2. Dirac Free Particle Equation 3.2 Page 13 Dirac Free Particle Equation The Dirac equation is written as i ∂ ∂ ∂ ∂ψ = [−i(α1 1 + α2 2 + α3 3 ) + βm]ψ ∂t ∂x ∂x ∂x ~ + βm)ψ = (−i~ α·∇ (3.7) In order to determine what the coefficients αi and β stand for, we can take the square of the operators in Dirac’s equation and compare it with KG equation. The Dirac equation with squared operator should be equivalent to the KG equation. So we can find the conditions set on αi and β. Squaring the operators on both sides gives 2 ∂ ~ + βm)(−i~ ~ + βm)ψ i ψ = (−i~ α·∇ α·∇ ∂t P P3 ∂2 ∂2 (3.8) = − 3i=1 αi2 (α α + α α ) ψ − i j j i i,j=1,i>j (∂xi )2 (∂xi ∂xj ) P ∂ + β 2 m2 ψ −im 3i=1 (αi β + βαi ) ∂xi So we assume that ψ in Eq.(3.7) must satisfy the KG equation, − ∂2 ~ 2 + m2 )ψ ψ = (−∇ ∂t2 (3.9) Then the conditions imposed on αi and β due to the above requirement are αi αj + αj αi = 0 i, j = 1, 2, 3 i 6= j αi β + βαi = 0 i, j = 1, 2, 3 αi2 = β 2 = I i = 1, 2, 3 (3.10) Using this set of equations we can construct an anticommutation relation [BJ95] {αµ , αν } = αµ αν + αν αµ = 2δµν I, ν, µ = 0, 1, 2, 3 (3.11) where α0 = β. For the coefficients αi and β to satisfy the anticommutation relation given in Eq.(3.11), they must possibly be square matrices with the minimum even dimension of 4. Thus the familiar choices for the αi0 s and β are the block matrices given by I 0 0 σi (3.12) ; β= αi = 0 −I σi 0 where I is the 2 × 2 identity matrix, 0 is also 2 × 2 matrix of zero elements, and σi are 2 × 2 Pauli matrices (see appendix A.3). As an immediate consequence, the wavefunction ψ must be a four-component column matrix, ψ1 ψ2 (3.13) ψ= ψ3 ψ4 Section 3.3. Free Particle Solution of Dirac Equation Page 14 By using Einstein summation notation, the Lorentz covariance form of the Dirac equation is written as (iγ µ ∂µ − m)ψ = 0 (3.14) where γ µ is a Dirac gamma matrix which replace β and αi , and is defined as γ 0 = β; γ k = βαk , k = 1, 2, 3 (3.15) or equivalently γ µ = (γ 0 , ~γ ) = (β, β~ α) (3.16) Similarly, the gamma matrices satisfy the anticommutation relation [BJ95] γ µ γ ν + γ ν γ µ = 2g µν I (3.17) and γ 0 γ 0 = −γ i γ i = I. The continuity equation for the Dirac equation becomes ∂µ j µ = 0 (3.18) where j µ = (ρ, J~ ) = ψ̄γ µ ψ is a four-vector current, and we have used ψ̄ = ψ † γ 0 which is Lorentz invariant. Thus the probability density becomes † ρ=ψ ψ= 4 X n=1 | ψn |2 (3.19) which is a positive (scalar) density as it is expected to be. If we multiply j µ by the electron charge −e it can be regarded as the electric current density or current-charge four-vector of the electron. That is, jeµ = −eψ̄γ µ ψ (3.20) 3.3 Free Particle Solution of Dirac Equation We begin by choosing a solution to be a plane wavefunction which has a form ψ(x, p, s) = Ne−ip·x u(p, s) (3.21) where p · x = pµ xµ = Et − p · x is a contraction of four-vectors p and x, s is the spin state of the particle and u(p, s) is four-component column matrix and is called the Dirac spinor; it can be written in terms of two two-component spinors, χ u(p, s) = C (3.22) φ where N and C are normalisation factors. It is worth noting the context in which we have used the Mandelstam variable s and the spin state s in order to avoid confusion. By plugging Eq.(3.22) Section 3.3. Free Particle Solution of Dirac Equation Page 15 into Eq.(3.21) and then Eq.(3.21) into Eq.(3.7) as well as recalling the 2 × 2 block expressions of β and αi0 s we get χ mI σ · p χ E = (3.23) φ σ · p −mI φ From which we obtain positive energy and negative energy solutions u(p, s) and v(p, s) given by ! σ·p s ! χs φ σ·p s (3.24) and v(p, s) = C u(p, s) = C E+m χ s φ E+m 1 where in the two expressions χs andφs each represents either of the two possible spin states ± 2 1 0 of a particle, that is either spin up or spin down . 0 1 In terms of gamma matrices and free particle spinors, the Dirac equation can be rewritten as (γ µ pµ − m)u(p, s) = 0 (3.25) Here we introduce the slash notation. If A is a four-vector then γ · A = γ µ Aµ ≡ A 6 (3.26) (6p − m)u(p, s) = 0 (3.27) Thus Eq.(3.25) becomes Moreover, the electron current can also be rewritten as jeµ = −eūγ µ u, which also satisfies the current conservation law ∂µ jeµ = 0. The Dirac equation for a spin- 12 particle, say, an electron interacting with external electromagnetic field having a four-vector potential Aµ = (φ, A), conσ of the electron that gives rise to tains information regarding the intrinsic magnetic moment −e 2m magnetic interaction. Here φ and A are the scalar and vector potentials of the field, respectively. 3.3.1 Normalisation of Dirac Spinor For the Dirac equation of the free particle, we can choose the plane wave solution of the form ψ = Nu(p, s)e−ip·x , s = 1, 2 (3.28) In order to determine the normalization factors, we use the covariant normalization [DD64] ū(p, s)u(p, s) = 1 or equivalently u† (p, s)u(p, s) = So by using this convention we get r m N= EV and (3.29) E m C= (3.30) r E m (3.31) where V is the normalization volume. This convention is especially useful when we deal with a system of massive particles. Moreover, it helps to get Lorentz invariant expression of the differential cross-section. 4. Quantum Electrodynamics of Electron-proton Elastic Scattering In this chapter, more emphasis is given to the dynamic description of the scattering process. Electromagnetic Interaction between spin- 12 Dirac particles is well described by the QED. In the study of two particle scattering, we need to know the rate of transition from the initial state to the final state. One way of addressing this problem is evaluating the Feynman diagram in order to find the invariant amplitude, which contains all the dynamics of the scattering under question. In other words, the physics of any scattering processes can be extracted from the corresponding Feynman diagram. 4.1 Feynman Diagrams In Quantum Electrodynamics, we can make the calculations of the relativistic scattering process easier by employing the Feynman diagrams (F D). A Feynman diagram is a graphical representation of a scattering process [Ron94]. It gives schematic description of the scattering processes that arise from the interaction of particles. In this study, we focus only on the lowest-order Feynman diagram which is often called as the tree diagram to represent the electron-proton elastic scattering due to the electromagnetic interaction. The interaction takes place via force carrier called virtual or “off-mass shell” photon which is specified by four-momentum transfer q 2 6= 0. The photon is spin-1 particle. 4.1.1 Components of Feynman Diagram Feynman diagrams are treated in Minkowski space-time. Basically, we need three elements to construct F D for any lowest-order scattering in EM fields: Straight lines for representing the free incoming and outgoing particles, wavy lines for virtual photon, and vertices where two incoming and outgoing solid lines and a single wavy line meet at the moment of interaction (see figure 4.1). (a) outgoing photon (b) incoming fermion Figure 4.1: (a) Wavy line representing the outgoing photons and (b) straight lines representing the incoming fermion. If we change the position of the dots to the other ends, we get lines that represent (a) the incoming photon and (b) the outgoing fermion. We use dots to represent lepton-photon vertices 16 Section 4.2. Feynman Rules for Tree Diagram Page 17 The building block for constructing F D of scattering, such as tree diagrams and other higher order Feynman diagrams, contains two straight (external) lines, a single wavy (internal) line and a vertex, which is a point where these three lines meet and such a graph is often called a primitive diagram (see figure 4.2). In EM two-body interactions, however, we need at least two particlephoton vertices, four incoming and outgoing solid lines and a single wavy line; such a diagram is called the lowest-order F D or tree diagram (see figure 4.3). e− γ e− Figure 4.2: Primitive diagram that illustrates the emission of photon from electron In evaluating the cross-section, we take our first step based on the assumption that both the electron and proton are spin- 12 elementary particles. p e− γ p e− Figure 4.3: The lowest-order tree diagram 4.2 Feynman Rules for Tree Diagram By using Feynman rules, we can extract multiplicative factors from a F D representing a particular scattering process in order to construct expressions for the invariant amplitude and the scattering cross-section which describe the process. The F D is drawn such that each of its elements corresponds to a factor in the express of invariant scattering amplitude. If we consider F D for fermion-fermion scattering due to EM interaction, the external (incoming and outgoing) lines are associated with Dirac spinors u(p, s) (for incoming particle) ū(p, s) (for outgoing particle) (4.1) (4.2) where p and s, respectively, are four-momenta and spin states (± 12 ) of a particle under study. The internal (wavy) line representing a virtual photon has got a propagator, −ig µν q2 (4.3) Section 4.3. Fermi’s Golden Rule Page 18 associated with it, where g µν is the metric tensor and q 2 is the square of four-momentum transfer. With each point vertex of F D we associate a factor −igγ µ (4.4) √ where g is a coupling constant. For EM interaction it is the charge of the positron e = 2πα where α = 1/137 is the fine-structure constant. Besides that, the four-momentum must be conserved at the vertices [Gri87]. In Quantum Electrodynamics, we represent a vertex by a dot if the incoming and outgoing particles are point-like, otherwise, we use a “blob”. 4.3 Fermi’s Golden Rule Expressions of cross-section have two basic ingredients: the invariant amplitude and phase space factor, which contains all kinematic aspects of scattering process [Gri87]. Here we develop the important expression for transition rate and cross-section using Fermi’s golden rule. 4.3.1 Transition Rate To begin with, we consider the scattering process of two particles 1 + 2 → 10 + 20 . The transition rate per unit volume, from the initial state to the final state [AH03], is given by Wf i = |Tf i |2 VT (4.5) where T is time interval of observation, V is the spatial volume of the interaction region and |Tf i |2 is a transition probability which is defined in terms of the transition matrix element Mf i as (4.6) Tf i = −iN1 N2 N10 N20 (2π)4 δ 4 (p1 + p2 − p01 − p02 )Mf i Thus (N1 N2 N10 N20 )2 [(2π)4 δ 4 (p1 + p2 − p01 − p02 )]2 |Mf i|2 Wf i = VT (4.7) Now we use the delta function property in the well defined spacial volume V and time interval T in space-time [GR03] [(2π)4 δ 4 (p1 + p2 − p01 − p02 )]2 = (2π)4 δ(0)(2π)4δ 4 (p1 + p2 − p01 − p02 ) = V T (2π)4 δ 4 (p1 + p2 − p01 − p02 ) (4.8) Wf i = (N1 N2 N10 N20 )2 (2π)4 δ 4 (p1 + p2 − p01 − p02 )|Mf i |2 (4.9) then we get It is worth noting that in our derivations thus far, we have indirectly used Dirac’s plane wavefunctions which represent particles in Quantum Field Theory. Section 4.4. Elastic Scattering of Electron from a Point Proton 4.3.2 Page 19 Cross-Section for Two-body Scattering We now define the differential cross-section dσ in terms of the transition rate per unit volume [HM84], that is Wf i dσ = (number of final states) (4.10) initial flux Experiments, however, measure the differential cross-section for scattering into a particular solid angle dΩ in the momentum space interval d3 p01 d3 p02 , where p01 and p02 are the final three-vector momenta of the electron and the proton, respectively [Ryd96]. Initial flux is the product of incident flux and the density of the target particle and hence given by 1 1 ((p1 · p2 )2 − m21 m22 ) 2 |v1 − v2 | = 2 (4.11) V2 V E1 E2 where v1 and v2 are the velocities, E1 and E2 are energies and m1 and m2 are masses of the incident and the target particles, respectively. The number of final states in momentum space interval d3 p01 d3 p02 is given by V V d3 p01 d3 p02 (4.12) 3 3 (2π) (2π) and thus the differential cross-section for a transition from initial state to final state is expressed as E1 E2 V 2 V V d3 p01 d3 p02 dσ = 1 Wf i 3 3 2 2 2 (2π) (2π) 2 ((p1 · p2 ) − m1 m2 ) m1 m2 2 4 4 0 0 = (4.13) 1 |Mf i | (2π) δ (p1 + p2 − p1 − p2 ) 2 2 2 2 ((p1 · p2 ) − m1 m2 ) m1 d3 p01 m2 d3 p02 × (2π)3 E10 (2π)3 E20 In a more convenient and compact form dσ = where 4m1 m2 ((p1 · p2 )2 1 − m21 m22 ) 2 |Mf i |2 dLips(s; p01 , p02 ) dLips(s; p01 , p02 ) = (2π)4 δ 4 (p1 + p2 − p01 − p02 ) m1 d3 p01 m2 d3 p02 , (2π)3 2E10 (2π)3 2E20 (4.14) (4.15) is the Lorentz invariant phase space factor for a system of two particles and s = (p1 + p2 )2 = (p01 + p02 )2 is a Mandelstam variable. We notice that other factors in Eq.(4.14) are also Lorentz invariant. 4.4 Elastic Scattering of Electron from a Point Proton In this section we derive the invariant amplitude Mf i by putting together the factors that we associated with F D earlier in section 4.2; we also take the initial and final spin states into consideration. Section 4.4. Elastic Scattering of Electron from a Point Proton 4.4.1 Page 20 Evaluation of Invariant Amplitude of Unpolarised Electron-proton Scattering p e− ū20 ū10 −igγ µ u1 −ig µν q2 −igγ µ u2 p e− Figure 4.4: The lowest-order Feynman diagram for electron-proton elastic scattering. Each factor of the invariant amplitude are associated to the corresponding element of the Feynman diagram. The process can be labelled as e+p→e+p 1 + 2 → 1 0 + 20 (4.16) (4.17) The initial and final free electron and free proton are described by wave functions of the form r m ψ= u(p, s)e−ip·x (4.18) EV r m ψ̄ = ū(p, s)e+ip·x (4.19) EV We start with the general way of expressing the invariant amplitude from the F D based on the Feynman rules. By referring to the Feynman diagram in figure 4.4, we begin by writing down the quantum mechanical expression of Mf i . That is Mf i = hFinal States|Interaction Operator|Inital Statesi = hPf , Sf |Ô|Pi, Si i (4.20) where Pi and Pf are the total 4-momenta and Si and Sf are the total spins states of the electron and the proton before and after scattering, respectively, in the scattering process. By using the Kronecker (tensor) product we can rewrite Eq.(4.20) as Mf i = g2 [ū(p01 , s01 ) ⊗ ū(p02 , s02 )] [γµ ⊗ γ µ ] [u(p1 , s1 ) ⊗ u(p2 , s2 )] q2 (4.21) Now we use the property of the Kronecker product for matrices, say, A, B, C, and D (A ⊗ B)(C ⊗ D) = AC ⊗ BD Thus we get Mf i = g2 [ū10 γµ u1 ][ū20 γ µ u2 ] q2 (4.22) (4.23) Section 4.4. Elastic Scattering of Electron from a Point Proton Page 21 where we have used a shorthand notation ui = u(pi , si ). This expression reminds us that Mf i is the product of the electron and the proton currents and also note that ū10 γµ u1 is simply a complex number. The complex conjugate of Mf i then becomes g2 [ū10 γµ u1 ]∗ [ū20 γ µ u2 ]∗ q2 g2 = 2 [ū1 γν u10 ][ū2 γ ν u20 ] q g2 = 2 [ū1 γν u10 ][ū2 γ ν u20 ] q M∗f i = (4.24) where we have used the fact that the complex conjugate of any scalar number is equal to its hermitian conjugate and with the help of the gamma matrices obey such properties as γ µ = γ µ , † γ 0 = γ 0 , and (γ 0 )2 = I. Hence the square of the invariant amplitude becomes 2 2 g 2 ∗ |Mf i| = Mf i Mf i = (4.25) [(ū10 γµ u1 )(ū20 γ µ u2 )][(ū1 γν u10 )(ū2 γ ν u20 )] q2 For unpolarised incoming and outgoing particles, however, we have to average over the initial spin states and sum over the final spin states in order to get the spin-averaged invariant amplitude. Then it becomes independent of spin states [Bra05, Gin04, Gri87]. That is, 1 X |Mf i |2 (4.26) |Mf i|2 = 4 0 0 s1 ,s2 ,s1 ,s2 After plugging Eq.(4.25) into Eq.(4.26) and collecting together factors that belong to the same particle, we may get a compact expression 2 2 g 2 |Mf i| = Lµν W µν (4.27) q2 where Lµν = 1X [(ū10 γµ u1 )(ū1 γν u10 )], 2 0 (4.28) s1 s1 W µν = 1X [(ū20 γ µ u2 )(ū2 γ ν u20 )] 2 0 (4.29) s2 s2 are lepton and hadronic tensors, respectively. For convenience, let’s first evaluate W µν . Here we rather choose to write it by using matrix indices as 1X (4.30) [(ū20 )α (γ µ )αρ (u2 )ρ (ū2 )ε (γ ν )εη (u20 )η ] W µν = 2 0 s2 s2 So now we are dealing with elements of matrices, thus we can change their position without affecting the expression. That is 1X W µν = (4.31) [(u20 )η (ū20 )α (γ µ )αρ (u2 )ρ (ū2 )ε (γ ν )εη ] 2 η Section 4.4. Elastic Scattering of Electron from a Point Proton Page 22 As we can see in Eq.(4.31), the spin sum is reduced to matrix trace. We also introduce the positive energy projection operator Λ(p, s)αβ = X u(p, s)α ū(p, s)β = s (6p + m) 2m (4.32) Thus P 1P 0 0 µ ν Λ(p , s ) γ ηα 2 2 αρ s02 Λ(p2 , s2 )ρε γεη 2 s2 0 6p2 + m2 1 P 6p 2 + m2 µ ν γαρ γεη = η 2 2m2 2m 2 ηα ρε 1 P 0 µ ν = ((6p 2 + m2 )γ (6p2 + m2 )γ )ηη 8m22 η 1 Tr[(6p 20 + m2 )γ µ (6p2 + m2 )γ ν ] = 8m22 W µν = (4.33) The method of reducing the spin sum to matrix trace is commonly called as Casimir’s trick [Gri87]. Now we are left only with the evaluation of the trace without using the explicit forms of Dirac spinors and gamma matrices. By applying trace theorems of gamma matrices (see appendix A.3), Tr[γ µ γ ν ] = 4g µν Trace of any product of an odd number of γµ 0 s is zero Tr[γ µ γ ν γ β γ α ] = 4[g µν g βα − g µβ g να + g µα g νβ ] then the hadronic tensor becomes 4 m22 0 µ ν µ ν µν W = Tr[6p 2 γ 6p2 γ ] + Tr[γ γ ] 8m22 4 1 0µ ν µν 0 2 µν = ([pµ2 p0ν 2 + p2 p2 − g (p2 · p2 )] + m2 g ) 2 2m2 1 µ 0ν µν 2 0 = [p p + pν2 p0µ 2 + g (m2 − p2 · p2 )] 2 2 2 2m2 (4.34) (4.35) q 2 µν g . 2 Similarly the electron tensor can be written as It is simple to show that g µν (m22 − p2 · p02 ) = Lµν = 1 [(p1 )µ (p01 )ν + (p1 )ν (p01 )µ − gµν (p1 · p01 − m21 )] 2 2m1 (4.36) After performing the contraction of lepton and hadronic tensors, we get 4 1 e [(p1 ·p2 )(p01 ·p02 )+(p1 ·p02 )(p2 ·p01 )−m21 p2 ·p02 −m22 p1 ·p01 +2m22 m21 ] (4.37) |Mf i |2 = 2 2 4 q 2m1 m2 This expression reminds us that the spin-averaged amplitude depends only on the initial and the final momenta of the electron and the proton. Section 4.4. Elastic Scattering of Electron from a Point Proton Page 23 Having calculated the square of the invariant amplitude, we proceed to evaluate the differential cross-section. The e-p scattering experiments are often carried out in the setting where the proton is initially at rest. As a consequence, we choose to work in the lab frame. Thus the four-momenta are given by p2 = (E2 , 0) = (m2 , 0), p1 = (E1 , p1 ), p01 = (E10 , p01 ) (4.38) We notice that p01 and p02 are not independent variables; they appeared in Eq.(4.13) as integration variables. So in order to restrict their values we have to integrate over d3 p01 and d3 p02 . However, we do not integrate over dΩ for we want to obtain the angular distribution of the differential cross-section. Moreover, we have to express dp01 in terms of spherical coordinates in momentum space as d3 p01 = |p01 |2 d|p01 |dΩe = |p01 |E10 dE10 dΩe d3 p01 = |p01 |dE10 dΩe ⇒ E10 (4.39) (4.40) where we have used the relativistic energy and momentum relation |p01 |2 + m21 = E102 . The subscript e is used because we prefer to look only at the scattered electron. Based on Eq.(4.38) in the lab frame the incoming flux factor is simplified as m1 m2 m1 m2 m1 1 = 1 = |p01 | ((p1 · p2 )2 − m21 m22 ) 2 (m22 (E12 − m21 )) 2 (4.41) Thus the invariant differential cross-section becomes Z m21 m2 1 d3 p02 2 δ 4 (p + p − p0 − p0 )|p0 |dE 0 dΩe dσ̄e = |M | fi 1 2 1 2 1 1 |p1 | E20 2π 2 (4.42) Note that we have used the same symbol dσ¯e before and after integration [Gri87]. The three dimensional delta function enables us to integrate over d3 p02 which is evaluated at p02 = p1 + p2 − p01 that gives rise to the three-vector momentum conservation and after integration we are left with 0 Z q 1 m21 m2 |p1 |dE10 dΩe 0 0 2 2 2 dσ̄e = |M | δ E + E − E − (p + p − p ) + m fi 1 2 1 2 1 1 2 |p1 | E20 (2π)2 (4.43) Here by reminding ourselves that we are working in the lab frame where the proton was initially at rest (i.e., p2 = 0), hence we write E20 as 1 p p 0 0 2 02 2 2 2 2 2 02 2 (4.44) E2 (E1 ) = E1 + E 1 − 2m1 − 2 E1 − m1 E1 − m1 cos θ + m2 In Eq.(4.44) E20 (E10 ) is just to indicate the dependence of E20 on E10 . Then Z 1 m21 m2 0 0 0 2 dσ̄e = 2 |Mf i | δ[f (E1 )]h(E1 )dE1 dΩe |p1 | (2π) (4.45) where f (E10 ) = E1 + E2 − E10 21 q q 2 2 2 02 2 2 02 2 − E1 − 2m1 + m2 + E1 − 2 E1 − m1 E1 − m1 cos θ (4.46) Section 4.5. Proton Form Factors and Electron-proton Elastic Scattering Page 24 p E102 − m21 . Now if we introduce the property of Dirac delta function (see appendix and h(E10 ) = E20 (E10 ) A.1): Z df (x) −1 df (x) −1 δ[f (x − xo )] = δ(x − xo ) and h(x)δ[f (x)] = h(xo ) (4.47) dx x=xo dx x=xo where xo is a root of f (x) = 0. Eq.(4.45) demands the first derivative as well as the zero of f (E10 ). Obviously, f (E10 ) = 0 does not only enforce the conservation of energy but also it does fix p01 and hence p02 . If we denote E1∗ as the positive energy solution for f (E10 ) = 0, after some steps we can write 1 2AB − 2A3 + C [B 2 − 2A2 B + A4 + m21 (C 2 − 4A2 )] 2 = (4.48) (C 2 − 4A2 ) p where we have denoted A = E1 + E2 , B = E12 − 2m21 + m22 and C = 2 E12 − m21 cos θ; and i h 0 02 2 − 21 E 2 − C(E − m ) 0 1 1 1 df (E1 ) (4.49) = −1 − h i1 0 p dE1 02 02 2 2 2 E1 + B − C E1 − m1 E1∗ As we can see Eq.(4.48) relates E10 to the initial electron energy and scattering angle. Note also that Eq.(4.49) is evaluated at E1∗ . Thus the differential cross-section given in Eq.(4.43) for the electron scattered through the solid angle element dΩe becomes dσ̄e 1 h(E1∗ ) m2 m2 2 |M | = 1 2p 2 fi df (E10 ) dΩe (2π) E1 − m21 dE 0 0 ∗ 1 E =E 1 (4.50) 1 Thus Eq.(4.50) is the one we compare to the Mott cross-section for the electron scattered from the point proton. Apparently, since most scattering experiments use unpolarized incident and target particles, we can also compare the outcomes of our calculation to the results of those appropriate experiments. But in this project the comparison is made only with the Rutherford and the Mott predictions for the scattering from a stationary point charge. The Mott cross-section is given by α2 dσe θ 2 2 (4.51) = 1 − β sin 4 θ 2 2 dΩe 2 4|p1 | β sin ( 2 ) where β = 4.5 |p1 | in natural units. It reduces to the Rutherford formula in the limit β → 0. E1 Proton Form Factors and Electron-proton Elastic Scattering The real proton is not a point particle and that is why Mott cross-section does not fit to experimental data for high energy probe by electrons [CR56, Jaf98]. So this suggests that we need Section 4.5. Proton Form Factors and Electron-proton Elastic Scattering Page 25 to add some other information regarding the proton structure to our calculations. In the previous section, we constructed the spin-averaged lepton and hadronic tensors Lµν and W µν for the unpolarized e-p scattering in the framework of QED. Unlike the lepton tensor for the electronphoton vertex, the hadronic tensor for proton-photon vertex is not completely understood and specified by QED. In this section, therefore, we generalize from point proton to composite proton by writing down the most general form of the second-rank tensor. So it is at this point that the concept of form factors enter the model. Because of the lack of a full understanding of the structure of hadrons, QED instead uses symmetry principles: Lorentz invariance and EM current conservation (gauge invariance) to obtain the general form of W µν [AH03]. So we use form factors to model the spacial extension of the particle whenever our knowledge about its internal structure is insufficient. From now on, we consider the proton as having structure in the finite volume. The hadronic tensor (given in Eq.(4.35)), which was extracted from the Feynman diagram (figure 4.4), does not take the structure of the proton into account; and it can also be rewritten as q 2 µν 1 µ 0ν ν 0µ µν p p2 + p2 p2 + g (4.52) W = 2m22 2 As we can see it takes three arguments p2 , p02 , and q. But they are not independent variable because p02 = p2 + q. Thus we notice that the hadronic tensor is the function of only two independent variables p2 and q. That is, W µν → W µν (p2 , q). Moreover, based on Eq.(4.33) we can deduce that the unknown quantity W µν (p2 , q) is certainly a second rank tensor. So in order to obtain the general form for the hadronic tensor, we must expand W µν in terms of the basis constructed from the two independent four-vector variables p2 and q and metric tensor g µν [AH03]. The available basis to build second-rank tensor are g µν , pµ2 pν2 , q µ q ν , pµ2 q ν + q µ pν2 , and pµ2 q ν − q µ pν2 . We note that the first four basis are symmetric and parity invariant, whereas, the last one is anti-symmetric (i.e. pν2 q µ − q ν pµ2 = −(pµ2 q ν − q µ pν2 ) swapping indices does not preserve the original expression). Thus the most general hadronic tensor may be written as W µν = −W1 g µν + W2 µ ν W3 µ ν W4 W5 p p + 2 (p2 q − q µ pν2 ) + 2 (pµ2 q ν + q µ pν2 ) + 2 q µ q ν 2 2 2 m2 m2 m2 m2 (4.53) where Wi are unknown functions of q 2 . For convenience, by collecting symmetric terms and anti-symmetric terms together, we can rewrite Eq.(4.53) as Thus W µν = WSµν + WAµν (4.54) Lµν W µν = Lµν (WSµν + WAµν ) (4.55) For an unpolarized electron, Lµν is completely symmetric. Since the contraction of symmetric tensor and anti-symmetric tensor is zero, we drop the second term at the right hand side of W3 eq.(4.55). In other words, the anti-symmetric term 2 (pµ2 q ν − q µ pν2 ) does not contribute to m2 |Mf i |2 [Gri87, HM84]. Then the possible form for W µν is W µν = −W1 g µν + W2 µ ν W4 µ ν W5 p p + 2 (p2 q + q µ pν2 ) + 2 q µ q ν 2 2 2 m2 m2 m2 (4.56) Section 4.5. Proton Form Factors and Electron-proton Elastic Scattering Page 26 These four functions are not independent [Gri87]. It is known that the W µν is the product of two proton currents; and consequently it satisfies the current conservation law at the hadronic vertex (or it is gauge invariant). That is, qµ W µν = qν W µν = 0 (4.57) We can find their relation by plugging Eq.(4.56) into Eq.(4.57). Hence qµ W µν = −W1 q ν + (p2 · q) W2 ν W4 ν 2 ν 2 W5 ν q + (p · q)q + q p q =0 2 2 +q 2 2 m2 m2 m22 Now by factoring pν and q ν out, we get W2 W4 2 W4 ν 2 W5 p2 · q 2 + q 2 p + −W1 + p2 · q 2 + q 2 q ν = 0 m2 m2 m2 m2 (4.58) (4.59) Eq.(4.59) holds if and only if the coefficients of pν and q ν vanish for pν and q ν are independent four-vectors. That is, W2 W4 p2 · q 2 + q 2 2 = 0 m2 m2 (4.60) W W4 5 2 −W1 + p2 · q 2 + q 2 = 0 m2 m2 which imply that W4 = − p2 · q W2 q2 and W5 = m22 W1 p2 · q + ( 2 )2 W2 2 q q Then substituting Eq.(4.61) back into Eq.(4.56) we get p2 · q µ ν W2 p2 · q 2 µ ν qµqν µ ν µ ν µν µν + 2 p2 p2 − 2 (p2 q + q p2 ) + ( 2 ) q q W = W1 −g + 2 q m2 q q (4.61) (4.62) The expression inside the parenthesis of the second term of Eq,(4.62) can be factorized into the form (pµ2 − Aq µ )(pν2 − Bq ν ). So after expanding this form and making a comparison with the p2 · q coefficient of W2 in Eq.(4.62) we find that A = B = 2 . Thus q p2 · q µ W2 p2 · q ν qµqν µ ν µν µν + 2 p2 − 2 q p2 − 2 q (4.63) W = W1 −g + 2 q m2 q q Had the anti-symmetric term not been dropped, we would have used qν W µν = 0 to show that W3 = 0 which would have supported our previous argument. By using the lepton tensor in Eq.(4.36), we contract the lepton tensor and the hadronic tensor, and by applying the current conservation law q µ Lµν = q ν Lµν = 0, the invariant amplitude becomes 2 2 e 2 Lµν Wµν |Mf i | = 2 q 2 2 (4.64) q2 2 W2 e 0 0 2 = W1 (p1 · p1 − 2m1 ) + 2 (p1 · p2 )(p2 · p1 ) + m2 q2 m2 4 Section 4.5. Proton Form Factors and Electron-proton Elastic Scattering Since we are working in the lab frame, by applying Eq.(4.38), we finally obtain 0 4 θ θ E E e 1 1 2 2 2W1 sin + W2 cos |Mf i |2 = 4 2 q m1 2 2 Page 27 (4.65) where E10 = E1 /(1 + (2E1 /m2 ) sin2 (θ/2)) and we neglect the electron rest-mass m1 in the numerators. Thus the differential cross-section can be expressed as dσe θ dσe E10 2 = W2 + 2W1 tan (4.66) dΩe lab dΩe M ott E1 2 This is known as Rosenbluth formula [GR03, Gri87]. Thus W1 (q 2 ) and W2 (q 2 ) are form factors. The measured quantity from the scattering experiment is the differential cross section. So by undertaking a series of e-p elastic experiment one can determine these form factors at various q 2 . This is done by varying E1 and θ for a fixed value of q 2 and out the same procedure for carrying θ (dσ /dΩ ) other values of q 2 . Then we plot (dσee/dΩee) lab against tan2 which is a straight line with slope M ott 2 (dσ /dΩ ) 2W1 and intersects the axis of (dσee/dΩee) lab at W2 . This may require extrapolation. Actually, M ott these two form factors can be combined to give the electric and magnetic form factors which, respectively, are used to describe the spacial distribution of the charge and the magnetic moment of the proton [BJ95]. 5. Simulation, Discussion and Conclusion 5.1 Simulation and Discussion For the e-p elastic scattering, we compare the derived differential cross-section to Rutherford and Mott cross-sections by regarding the target proton as the Dirac point particle. In fact, this treatment holds in the framework of Quantum Electrodynamics in which the structure of the hadron is unknown. Rutherford prediction does not take into account the spin of the particles. So in order to get the Rutherford cross-section for e-p scattering we have to assume that both are spin-0 particles and with relativistic correction it is given by [GR03, Ter06] dσe α2 = dΩe 4|p1 |2 β 2 sin4 ( θ2 ) (5.1) where E1 is the energy of incident electron in the lab frame. Eq.(5.1) only tells us that the scattering is from the point-charge target. The Mott cross-section is a modification of Rutherford cross-section for the scattering from the point-like target by replacing spin-0 by spin- 21 and was given by Eq.(4.51). 100 100 Our calculation at 10.0MeV Our calculation at 50.0MeV Our calculation at 90.0MeV Rutherford cross section at 10.0MeV Rutherford cross section at 50.0MeV Rutherford cross section at 90.0MeV 1 Differential cross section (barns/steradian) Differential cross section (barns/steradian) 1 0.01 1e-04 1e-06 1e-08 1e-10 0.01 1e-04 1e-06 1e-08 1e-10 1e-12 1e-12 (a) Our calculation at 10.0MeV Our calculation at 50.0MeV Our calculation at 90.0MeV Mott cross section at 10.0MeV Mott cross section at 50.0MeV Mott cross section at 90.0MeV 1e-14 0 20 40 60 80 100 120 Scattering angle (degree) 140 160 180 (b) 0 20 40 60 80 100 120 Scattering angle (degree) 140 160 Figure 5.1: Comparison of the calculated cross-section with (a) Rutherford differential crosssection and (b) Mott differential cross-section for 10 MeV, 50 MeV, and 90 MeV incident energies of electron. In our simulations, we use the final result of our calculation of differential cross-section written in Eq.(4.50). For the invariant amplitude we take the expression given in Eq.(4.37). Since we are working in the lab frame, we also use the values given for momenta in Eq.(4.38) as inputs in our Python piece of code. The study of the e-p elastic process is done in the energy range m1 < E1 < m2 of the incident electron. The mass of electron is m1 ≈ 0.50 MeV and the mass of the proton is m2 ≈ 939.0 MeV. Moreover, we choose to work in natural units so that mass, momentum and energy have the same units, MeV. The unit of space-time coordinates is (MeV)−1 . Note also that |Mf i |2 and dσ̄e /dΩe are evaluated at E10 = E1∗ given in Eq.(4.48) to 28 180 Section 5.1. Simulation and Discussion Page 29 obey the requirement of energy conservation. The units of differential cross-section and scattering angle are barns per steradian and degree, respectively. We fix the initial energy of electron, in order to observe the dependence of differential cross-section on the scattering angle and the vice versa. 1 1 Our calculation at 100.0MeV Our calculation at 150.0MeV Our calculation at 350.0MeV Our calculation at 700.0MeV Rutherford cross section at 100.0MeV Rutherford cross section at 150.0MeV Rutherford cross section at 350.0MeV Rutherford cross section at 700.0MeV 0.01 Differential cross section (barns/steradian) Differential cross section (barns/steradian) 0.01 0.0001 1e-006 1e-008 1e-010 0.0001 1e-006 1e-008 1e-010 1e-012 1e-012 1e-014 0 (a) Our calculation at 100.0MeV Our calculation at 150.0MeV Our calculation at 350.0MeV Our calculation at 700.0MeV Mott cross section at 100.0MeV Mott cross section at 150.0MeV Mott cross section at 350.0MeV Mott cross section at 700.0MeV 20 40 60 80 100 120 Scattering angle (degree) 140 160 180 0 (b) 20 40 60 80 100 120 Scattering angle (degree) 140 160 180 Figure 5.2: Comparison of the calculated cross-section with (a) Rutherford differential crosssection and (b) Mott differential cross-section for 100 MeV, 150 MeV, 350 MeV and 700 MeV incident energies of electron. As we can see from figure 5.1 (a) and 5.2 (a), the calculated cross-section shows deviation from the Rutherford cross-section, whereas, in figure 5.1 (b) and 5.2 (b) our result coincides with the Mott cross-section as expected. But for large E1 , our results separates from Mott cross section at larger angles closer to 180o . As θ increases dσ̄e /dΩe decreases. 9e-008 260 scattering angle= 40.0degrees scattering angle= 50.0degrees scattering angle= 60.0degrees scattering angle= 70.0degrees (a) 7e-008 220 6e-008 5e-008 4e-008 3e-008 200 180 160 140 2e-008 120 1e-008 100 0 100 for incident energy=100.0MeV for incident energy=150.0MeV for incident energy=200.0MeV for incident energy=250.0MeV 240 Energy of scattered electron (MeV) Differential cross section (barns/steradian) 8e-008 80 150 200 250 300 350 Incident energy (MeV) 400 450 500 (b) 0 20 40 60 80 100 120 Scattering angle (degree) 140 160 Figure 5.3: (a) Angular distribution of differential cross-section versus energy of the initial electron. (b) Simulation for energy the scattered electron versus the scattering angle at incident energies 100 MeV, 150 MeV, 200 MeV, and 250 MeV Figure 5.3 (a) is the illustration of the dependence of differential cross-section on the initial electron beam energy. The simulation is done for θ = 40, 50, 60, and 70 degrees. Figure 180 Section 5.2. Conclusion Page 30 5.3 (b) is the plot of the energy of scattered electrons as a function of scattering angle which was calculated before in chapter four (see Eq.(4.48)). For a particular incident electron energy, as the scattering angle increases, the scattered electron energy decreases. In other words, this means that whenever there is relatively high four-momentum transfer from electron to the proton during the moment of interaction, the electron scattered at relatively large angle. Note that the maximum scattering angle is 180o . 5.2 Conclusion In the relativistic scattering process, the Feynman diagram and the Feynman rules are the powerful tools used to evaluate the invariant amplitude which is proportional to the scattering cross-section. The result of this essay is found to be consistent with the Mott cross-section for the electron scattered from the spin- 12 point-like Dirac proton particularly when incident electron energy, E1 is much less than the rest-mass of the proton, m2 . It was also shown that the scattering crosssection depends on the two independent quantities: the energy of incident beam of electron and the scattering angle. In addition, in the lab frame, once the initial state four-vector momenta and the scattering angle are known, it is possible to determine the final state four-vector momenta. The other important outcome is that a relatively large number of electrons are scattered at small scattering angles. As the scattering angle increases at a fixed value of E1 , the number of scattered particles decreases. Also, at a particular scattering angle, as the incident beam energy increases, the differential cross-section decreases. Moreover, by applying symmetric principle, one can obtain more general way of describing the internal structure of the proton from the study the e-p elastic scattering. However, the detailed information of the structure of the proton can be found from the QCD treatment of e-p inelastic scattering. Therefore, this understanding can be extended to study scattering processes involving more than two particles in the final states. Finally, the writer wants to undertake further theoretical and experimental study of high energy electron-proton elastic and inelastic scattering in order to investigate the quark model of nucleons that describes the proton as the bound state of three quarks. Appendix A. Background Mathematics A.1 The Dirac Delta Function The Dirac delta function, δ(x), is an infinitely high, infinitesimally narrow spike at the origin [Gri87], with area of unity; mathematically it is defined as ∞ if x = xo (A.1) δ(x − xo ) = 0 if x 6= xo and Z ∞ −∞ δ(x − xo )dx = 1 (A.2) Properties 1). It follows from the definition that for any ordinary function f (x) Z ∞ Z ∞ f (x)δ(x − xo )dx = f (xo ) δ(x − xo )dx = f (xo ) −∞ and we can also have Z x2 x1 (A.3) −∞ f (x)δ(x − xo )dx = f (xo ) if x1 < xo < x2 0 if xo > x2 or xo < x1 (A.4) The three dimensional delta function δ 3 (~x) can be defined as δ 3 (~x) = δ(x)δ(y)δ(z), (A.5) and for any function f (~x) = f (x, y, z) integrated over a volume V gives Z f (~xo ) if ~xo lies inside V 3 (A.6) f (~x)δ (~x − ~xo ) = 0 if ~xo lies outside V V R R 1 2). If k is a non-zero real number, then δ(kx)dx = |k| δ(x)dx and for any arbitrary function f (x) we have Z Z f (xo ) f (xo ) f (x)δ(k(x − xo ))dx = δ(x − xo )dx = (A.7) |k| |k| 3). If function g(x) is an argument of delta function and if xo is the solution for g(x) = 0, then it follows Eq.(A.7) that Z Z 1 1 δ[g(x)]dx = (A.8) δ(x − x )dx = o |g 0(x)|x=xo |g 0(xo )| and for any continuous function f (x) at x = xo Z Z 1 f (xo ) f (x)δ[g(x)]dx = f (xo ) δ(x − x )dx = o |g 0 (x)|x=xo |g 0(xo )| 31 (A.9) Section A.2. Relativistic Notation and Four-vector Formalism A.2 Page 32 Relativistic Notation and Four-vector Formalism Any point or event in Minkowski space or space-time is specified by xµ , µ = 0,1,2,3 and explicitly written as [Bra05] xµ = (x0 , x1 , x2 , x3 ) = (ct, x, y, z) (A.10) = (ct, ~x) Any mathematical object aµ is called a four-vector. Four-vector, with upper index, is called contravariant and one with, lower index, is called covariant vector; and the contravariant is related to its covariant counterpart through the metric tensor gµν as X X xµ = gµν xν ; xµ = g µν xν (A.11) ν ν where the Minkowski metric is defined as gµν = g µν 1 0 0 0 0 −1 0 0 = 0 0 −1 0 0 0 0 −1 (A.12) Here we introduce Einstein summation that drops the summation symbol for expresP convention µ sions with repeated indices. That is µ gµν x → gµν xµ . Thus Eq.(A.11) rewritten as xµ = gµν xν ; xµ = g µν xν (A.13) Minkowski metric gµν is defined by the relation g µν gµα = δαν I (A.14) where δαν is the Kronecker delta and I is 4 × 4 identity metric. It is worth noting that here also we used the summation convention. If aµ and bµ are four-vectors, their scalar product is given by a · b = (a0 b0 − a1 b1 − a2 b2 − a3 b3 ) = gµν aν bµ = aµ bµ (A.15) Similarly, the scalar product of xµ with itself is x · x = gµν xν xµ = (c2 t2 − ~x · ~x) (A.16) The covariant gradient is obtained by differentiating with respect to contravariant four-vector xµ . That is, ∂ ∂ ∂ ∂ ∂ , , , ∂µ = µ = ∂x c∂t ∂x ∂y ∂z (A.17) ∂ ~ ,∇ = c∂t and ∂ µ µν = g ∂ν = ∂ ~ , −∇ c∂t (A.18) Section A.3. Dirac Algebra of Gamma Matrix and Trace Theorems Page 33 By taking their scalar product we get µ = ∂ ∂µ = ∂2 ~2 −∇ c2 ∂t2 (A.19) called the d’Alembertian operator. A.3 A.3.1 Dirac Algebra of Gamma Matrix and Trace Theorems Pauli Matrices These are the three Hermitian, unitary 2 × 2 matrices: 0 1 0 −i σ1 = , σ2 = , 1 0 i 0 σ3 = 1 0 0 −1 (A.20) product rules σ12 = σ22 = σ32 = I {σi , σj } = 2δij A.3.2 (anticommutation relation) (A.21) Dirac Matrices These are four unitary traceless 4 × 4 matrices I 0 0 σk 0 k , γ = , γ = −σ k 0 0 −I k = 1, 2, 3 (A.22) Here I is 2 × 2 identity matrix and 0 is 2 × 2 matrix of zero. and these matrices satisfy anticommutation relation {γ µ , γ ν } = 2g µν I (A.23) where g µν is the metric tensor. We notice that gµν g µν = 4. Now we define the fifth Dirac metric out of the four ones 0 I 5 0 1 2 3 (A.24) γ = iγ γ γ γ = I 0 (γ 0 )2 = (γ 5 )2 = I, and (γ k )2 = −I, where k = 1, 2, 3. And γ 5 γ µ = −γ µ γ 5 . Other properties of gamma matrices are γ µ = γ 0 (γ µ )† γ 0 = γ µ ; iγ 5 = γ 0 (γ 5 )† γ 0 = iγ 5 ; γ µ γ 5 = γ 0 (γ 5 )† (γ µ )† γ 0 = γ µ γ 5 (A.25) Section A.3. Dirac Algebra of Gamma Matrix and Trace Theorems A.3.3 Page 34 Trace Theorems Here we recall the slash notation γ · a = 6a we have introduced in chapter three. Theorem 1 Tr[γ µ γ ν ] = 4g µν (A.26) Proof: we use the anticommutation relation γ µ γ ν + γ ν γ µ = 2g µν . 1 Tr[γ µ γ ν ] = Tr[γ µ γ ν + γ ν γ µ ] 2 = g µν Tr[I] = 4g µν (A.27) And Tr[γ µ γ ν γ α γ β ] = 4(g µν g αβ − g µα g νβ + g µβ g να ) (A.28) Tr[γ µ γ ν γ α γ β ] = Tr[(2g µν I − γ ν γ µ )γ α γ β ] = 2g µν Tr[γ α γ β ] − Tr[γ ν γ µ γ α γ β ] = 8g µν g αβ − Tr[γ ν γ µ γ α γ β ] (A.29) µ ν ν µ µν Proof: We use the the anticommutation relation γ γ + γ γ = 2g . We do the same for Tr[γ ν γ µ γ α γ β ]. Now we apply the anticommutation relation for γ µ γ α also, then we get Tr[γ ν γ µ γ α γ β ] = Tr[γ ν (2g µαI − γ α γ µ )γ β ] (A.30) = 8g µα g νβ − Tr[γ ν γ α γ µ γ β ] Once again by using γ µ γ β + γ β γ µ = 2g µβ I, we get Tr[γ ν γ α γ µ γ β ] = 8g µβ g να − Tr[γ ν γ α γ β γ µ ] (A.31) By putting all these togather, we get Tr[γ µ γ ν γ α γ β ] = 8g µν g αβ − 8g µα g νβ + 8g µβ g να − Tr[γ ν γ α γ β γ µ ] (A.32) By recalling the cyclic property, we get Tr[γ µ γ ν γ α γ β ] = 4(g µν g αβ − g µα g νβ + g µβ g να ) (A.33) The trace of an odd number of γ matrices is zero. (A.34) Theorem 2 Proof: For n odd, Tr[6a1 · · · a 6 n ] = Tr[6a1 · · · a 6 nγ 5γ 5] = Tr[γ 5 6a1 · · · a 6 nγ 5] n = (−1) Tr[6a1 · · · a 6 nγ 5γ 5] =0 (A.35) Section A.3. Dirac Algebra of Gamma Matrix and Trace Theorems Page 35 Theorem 3 For any arbitrary four-vectors a and b, Proof: Tr[6a6b] = 4a · b (A.36) 1 Tr[6a6b] = Tr[6b6a] = Tr[6a6b + 6b6a] 2 1 = Tr[aµ bν (γ µ γ ν + γ ν γ µ )] 2 1 = Tr[aµ bν (2g µν )] 2 = Tr[a · b] = a · bTr[I] = 4a · b (A.37) Theorem 4 For any arbitrary four-vectors a, b, c and d, Tr[6a6b6c6d] = 4[(a · b)(c · d) − (a · c)(b · d) + (a · d)(b · c)] (A.38) Proof: From theorem 3 we have 6a6b + 6b6a = 2a · b. By using this and the cyclic property Tr[6a1 · · · a 6 n ] = Tr[6an a1 · · · a 6 n−1 ] Tr[6a6b6c6d] = 2a · b − Tr[6b6c6d6a] (A.39) Repeating this trick a number of times gives Tr[6a6b6c6d] = 2(a · b)4(c · d) − 2((a · b)4(a · c))Tr[6b6d] + Tr[6b6c6d6a] = 8(a · b)(c · d) − 8(a · c)(b · d) + 8(b · c)(a · d) − Tr[6b6c6d6a] (A.40) The rest it trivial. By using Eq.(A.28) we can also show that Tr[6aγ6bγ] = 4(aµ bν + aν bµ − a · bg µν ) (A.41) Appendix B. The Python Code #!usr/lib/python from Numeric import * import math import Gnuplot #************************************************************************* #This function returns the scalar product two four-vectors by using metric tensor #************************************************************************* def metric_tensor():#Eq.(A.12) m_tensor=zeros([4,4]) for i in range(len(m_tensor[0])): for j in range(len(m_tensor[0])): if i == j: if i==0: m_tensor[i][j]=1 else: m_tensor[i][j]=-1 else: m_tensor[i][j] return m_tensor def contractionM(x, y): y_transpos = reshape(y, (y.shape[1], y.shape[0])) y_con = matrixmultiply(metric_tensor(), y_transpos) result = matrixmultiply(x, y_con)[0][0] return result #************************************************************************** #This function returns the square of four-monentum transfer q^2 #************************************************************************** def Four_momentum_square(a,b): q = a-b q2 = contractionM(q,q) return q2 #*************************************************************************** #The Invariant Amplitude (matrix element) #*************************************************************************** def matrix_element():#Eq.(4.37) ga = Four_momentum_square(p1i, p1f) sta0 = ga**2 ; sta1 = (4*math.pi*alpha)**2 sta2 = 2*(m1*m2)**2; sta3 = (p3p4)*(p1p2) sta4 = (p2p3)*(p1p4); sta5 = m2**2*(p1p3) sta6 = m1**2*(p2p4) ME2 = sta1*(sta3 +sta4 - sta5 -sta6 + sta2)/(sta0*sta2) return ME2 #*************************************************************************** 36 Page 37 #The Differentail Cross-Section (our result) #*************************************************************************** def Scattering_cross_section():#Eq.(4.50) car0 = (m1**2*m2)/(2*math.pi)**2 car1 = math.sqrt(E1i**2 - m1**2) car2 = matrix_element() H = math.sqrt(E1f**2 - m1**2)/E2f Df = abs(-(2*E1f-(cc*E1f)/math.sqrt(E1f**2-m1**2)) /(2*math.sqrt(-cc*math.sqrt(E1f**2-m1**2)+E1f**2+bb))-1)#Eq.(4.49) Cross_section = (car0*car2*H)/(Df*car1) return Cross_section #************************************************************************** #This one returns the Rutherford’s cross-section #************************************************************************** def Rutherford_cs():#Eq.(5.1) rr = alpha=1.0/137; ad = arad/2 dd = 4*(E1i**2 - m1**2)**2*(math.sin(ad))**4 ruth = (rr**2*E1i**2)/dd return ruth #*************************************************************************** #This one returns the Mott cross-section #*************************************************************************** def Mott_sc():#Eq.(4.51) mo2 = Rutherford_cs()*(1-((E1i**2 - m1**2)/(E1i**2)) *(math.sin(arad/2))**2) return mo2 #**************************************************************************** #This one returns all the outputs such as differential cross-sections, energy #of the scattered electron, the four-momenta. Few modifiacation is required. #we choose to work in natural units #**************************************************************************** gn = Gnuplot.Gnuplot(persist=1) m1 = 0.50#MeV m2 = 939.0#MeV alpha = 1.0/137#fine structure constant Start = 10.0 grandO =[]; grandM =[]; grandR =[] for E1i in arange(Start, 1000.0, 80.0): p1ix = 0.0; p1iy = 0.0; p1iz = math.sqrt(E1i**2-m1**2) E2i = m2; p2ix = p2iy = p2iz = 0.0 p1i = array([[E1i, p1ix, p1iy, p1iz]]) p2i = array([[E2i, p2ix, p2iy, p2iz]]) list_scatO = []; list_scatR = []; list_scatM=[];list_scatE=[] for angle in arange(1.0,180.0, 2.25): arad = angle*math.pi/180.0 aa = E1i + E2i; bb = E1i**2 - 2*m1**2 + m2**2 cc = 2*math.sqrt(E1i**2 - m1**2)*math.cos(arad) Page 38 #E1f=(cc*math.sqrt((cc**2-4*aa**2)*m1**2+bb**2-2*aa**2*bb+aa**4) +2*aa*bb-2*aa**3)/(cc**2-4*aa**2)#this does’t obey conservation E1f = -(cc*math.sqrt((cc**2-4*aa**2)*m1**2+bb**2-2*aa**2*bb+aa**4) -2*aa*bb+2*aa**3)/(cc**2-4*aa**2)#Eq.(4.48) E2f = math.sqrt(bb + E1f**2 - cc*math.sqrt(E1f**2 - m1**2)) p1fx = math.sqrt(E1f**2 - m1**2)*math.sin(arad); p1fy = 0 p1fz = math.sqrt(E1f**2 - m1**2)*math.cos(arad) p1f = array([[E1f, p1fx, p1fy, p1fz]]) p2f = p1i+p2i-p1f#cosevation of for momentum at the calculated E1f. p2fx = p2f[0][1]; p2fy = p2f[0][2]; p2fz = p2f[0][3] p2f = array([[E2f, p2fx, p2fy, p2fz]]) p1p2 = contractionM(p1i, p2i); p1p3 = contractionM(p1i, p1f) p1p4 = contractionM(p1i, p2f); p2p3 = contractionM(p2i, p1f) p2p4 = contractionM(p2i, p2f); p3p4 = contractionM(p1f, p2f) our_cs = Scattering_cross_section()#Our calculation ruth = Rutherford_cs()#Rutherford cross section mott = Mott_sc()#Mott cross section list_scatO.append([angle, our_cs]) list_scatR.append([angle, ruth]) list_scatM.append([angle, mott]) pdataO = Gnuplot.Data(list_scatO, with ="lines", title = ’Our calculation at ’+str( E1i)+"MeV") pdataM = Gnuplot.Data(list_scatM, with ="points", title = ’Mott cross section at ’+str( E1i)+"MeV") pdataR = Gnuplot.Data(list_scatR, with ="points", title = ’Rutherford cross section at ’+str( E1i)+"MeV") grandO.append(pdataO) grandM.append(pdataM) grandR.append(pdataR) gn(’set xlabel "Scattering angle (degree)"’); gn(’set logscale y’) gn(’set grid’) gn(’set ylabel "Differential cross section (barns/steradian)"’) gn.plot(grandO[0],grandO[1],grandO[2],grandO[4]) gn.replot(grandM[0],grandM[1],grandM[2],grandM[4]) #gn.replot(grandR[0],grandR[1],grandR[2],grandR[4]) Acknowledgement First of all, I would like to lift the name of my lord on high. He is my shepherd. Jesus has led and comforted my life with his mighty hands and heavenly blesses and enabled me to bring this project to light. He is the prince of peace, the everlasting father, and the almighty God. His kingdom shall not end. I am grateful to the genuine supervision of Dr. Brandon Van der Ventel, who spent his invaluable time and effort guiding me in the proper direction as well as developed a strong confidence in me to do my level best. Without him all my effort would have been in vain. I would like to give many thanks to all AIMS staffs especially Henry Amuasi who committed himself to share my burdens throughout this project and gave me constructive comments and suggestions, Jan Groenewald who helped me a lot in fixing the problem with latex and also Ambrose Chongo who edited the earlier draft of this paper. I also appreciate Dr. Sam Webster who tutored me for more than a semester in a very friendly way. I would also like to express my heart felt gratitude to the founders of AIMS especially Prof. Neil Turok and Prof. Fritz Hahne through whom the widows of opportunity opened for African students. Finally, I am indebted to Girma Goro and all AIMS2006/7 students who have shown me the spirit of friendship, forgiveness and given me courage during my stay at AIMS. 39 Bibliography [AH03] I.J.R. Aitchison and A.J.G. Hey, Gauge theories in particle physics, Institute of Physics, London, 2003. [BJ95] W.E. Burcham and M. Jobes, Nuclear and particle physics, Longman, 1995. [Bra05] Andreas Brandhuber, Relativistic quantum ph.qmul.ac.uk/ andreas/RQM/rqm.html, 2005. [CR56] E.E. Chambert and R.Hofstadter, Phys. Rev. 103 (1956). [DD64] Bjorken J. D and Drell, Relativistic quantum mechanics, McGraw-Hill, 1964. [Gin04] Douglas M. Gingrich, Feynman rules for tree graphs, http://www.phys.ualberta.ca, 2004. [GR03] W. Greiner and J. Reinhardt, Quantum electrodynamics, Springer, 2003. [Gri87] David Griffiths, Introduction to elementary particles, Wiley, 1987. [HM84] Francis Halzen and Alan D. Martin, Quarks and leptons: an introductory course to modern particle physics, wiley, 1984. [Jaf98] Robert L. Jaffe, hep-ph/9811327 2 (1998). [Joa84] C. J. Joachain, Quantum collision theory, Elsevier Science, North Holland, 1984. mechanics, http://www.strings. [McM93] Sara McMury, Quantum mechanics, Addison-Wesely, 1993. [Ron94] William B. Ronlick, The fundamental particles and their interaction, 1994. [Ros50] M.N. Rosenbluth, Phys. Rev. 79 (1950). [Ryd96] Lewis H. Ryder, quantum field theory, Cambridge University Press, 1996. [Ter06] Dimitri Terryn, Scattering cross sections, blogspot.com/2006/06/scattering-cross-sections, 2006. 40 http://www.stringschool.