Lec 28

Topic 5
Interactions and Feynman Diagrams
Lecture 28
Basic Formulas for Cross Sections and Decay Rates
The Scattering cross section is defined by Peskin-Schroeder on page 100 by the formula
Number of scattering events of a given type
.
A`AB `B A
where A is the overlapping cross sectional area of the two beams. The differential cross section is
!

3
Y
X
1
d pf 1
2
4 (4)
d =
|M(pA, pB {pf })| (2)
pA + pB
pf ,
2EA2EB |vA vB |
(2)3 2Ef
f
Eq. (4.79). The matrix element M is simply the sum of Feynman diagrams to which the Feynman rules
are applied. The formula for the decay of an unstable particle is given in Eq. (4.86)
!

3
X
Y
d pf 1
1
2
4 (4)
d =
|M(mA {pf })| (2)
pA
pf .
2mA
(2)3 2Ef
f
PHY 510
10/30/2013
Topic 5
Lecture 28
FeynCalc: Tools for Field Theory Calculations

Download and install FeynCalc from the website http://www.feyncalc.org/ or directly from a Mathematica Notebook by invoking the Mathematica Import function:
Import["http://www.feyncalc.org/install.m"]
The second In definition statement Imports and evaluates the Mathematica code Bhabha.m from Dr. Romaos
web site http://porthos.ist.utl.pt/CTQFT/node3.html.
PHY 510
10/30/2013
Topic 5
Lecture 28
Bhabha Scattering Cross Section using FeynCalc

J.C. Romao, Computational Techniques in Quantum Field Theory has a good tutorial and example codes
for relativistic perturbation theory calculations using FeynCalc and other computer programs.
Bhabha scattering is the elastic electron-positron scattering process
e(p1) + e+(p2) e(p3) + e+(p4) .
The matrix element for the process is a sum of t-channel and s-channel contributions
M = M1 + M2
e2
e2
= u(p3) u(p1) v(p2)v(p4) v(p2) u(p1) u(p3)v(p4) ,

t
s
where the Mandelstam variables s, t, u are defined by
t = (p1 p3)2 ,
PHY 510
s = (p1 + p2) ,
u = (p1 p4)2 ,
3
s + t + u = 4m2e .
10/30/2013
Topic 5
Lecture 28
Feynman Rules for QED

From Peskin-Schroeder pages 118 and 123:
PHY 510
10/30/2013
Topic 5
Lecture 28
Mathematica Code Bhabha.m

Mathematica Input Syntax defines the meaning of Mathematica programs.
The code begins with a Mathematica Bracketed comment.
(* Program to Calculate the Traces of Eq. 4.43 de TC
Last Version: 21/04/2013
Author: Jorge C. Romao
email: [email protected]
*)
It is generally a good idea to clear Mathematicas memory of all previous calculations to avoid conflicting
definitions and other strange bugs. This is done by invoking Remove on all symbols in the global Standard
Namespace
Remove["Global*"]
The next statments invoke FeynCalc functions to Define Mathematica Functions to generate
, , 6 p = p = p , g
PHY 510
, p ,
, 1nn , p q , 6 p + m , R,L
1 5
=
.
2
10/30/2013
Topic 5
Lecture 28
dm[mu_]:=DiracMatrix[mu]
dm[5]:=DiracMatrix[5]
ds[p_]:=DiracSlash[p]
mt[mu_,nu_]:=MetricTensor[mu,nu]
fv[p_,mu_]:=FourVector[p,mu]
epsilon[a_,b_,c_,d_]:=LeviCivita[a,b,c,d]
id[n_]:=IdentityMatrix[n]
sp[p_,q_]:=ScalarProduct[p,q]
li[mu_]:=LorentzIndex[mu]
prop[p_,m_]:=ds[p]+m
PR=(1 + dm[5])/2
PL=(1 - dm[5])/2
Next, define and set variables Line1,...,Line5 to represent products of Dirac matrices
Line1:= ds[p3] . dm[mu] .
ds[p1] . dm[nu]
ds[p4] . dm[nu]
ds[p1] . dm[nu]
ds[p4] . dm[nu]
PHY 510
10/30/2013
Topic 5
Lecture 28
ds[p1] . dm[nu] . ds[p2] . dm[mu] . ds[p4] . dm[nu]
Next invoke the FeynCalc functions

Tr to perform the Trace of Dirac matrix products,
Contract to contract repeated Lorentz indices
and the Mathematica Simplify function to simplify the resulting expressions:
ans1= Simplify[Contract[Tr[Line1] Tr[Line2]]]
ans2= Simplify[Contract[Tr[Line3] Tr[Line4]]]
ans3= Simplify[Contract[Tr[Line5]]]
ans=ans1/4/t^2+ans2/4/s^2-2*ans3/4/s/t
Define a List of Replace Rules to make substitutions in an expressions
dot={sp[p1,p2]->s/2,sp[p3,p4]->s/2,sp[p1,p3]->-t/2,sp[p2,p4]->-t/2,
sp[p1,p4]->(s+t)/2,sp[p2,p3]->(s+t)/2}
PHY 510
10/30/2013
Topic 5
Lecture 28
ReplaceAll
res= e^4 ans /. dot
and Print the answer using LATEX formatting:
Print["\!\(\*FractionBox[\(1\), \(4\)]\) \!\(\*UnderscriptBox[\(\[Sum]\),
\(spins\)]\)|M\!\(\*SuperscriptBox[\(|\), \(2\)]\) = ",res]
(********************* Fim do Programa **********************)
PHY 510
10/30/2013

Lec 28

Uploaded by

Copyright:

Available Formats

Lec 28

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Lec 28

Uploaded by

Copyright:

Available Formats

Topic 5

Interactions and Feynman Diagrams

Basic Formulas for Cross Sections and Decay Rates

Interactions and Feynman Diagrams

FeynCalc: Tools for Field Theory Calculations

Interactions and Feynman Diagrams

Bhabha Scattering Cross Section using FeynCalc

= u(p3) u(p1) v(p2)v(p4) v(p2) u(p1) u(p3)v(p4) ,

Interactions and Feynman Diagrams

Feynman Rules for QED

Interactions and Feynman Diagrams

Mathematica Code Bhabha.m

Interactions and Feynman Diagrams

Line2:= ds[p2] . dm[mu] .

Line3:= ds[p2] . dm[mu] .

Line4:= ds[p3] . dm[mu] .

Line5:= ds[p3] . dm[mu] .

Interactions and Feynman Diagrams

ds[p1] . dm[nu] . ds[p2] . dm[mu] . ds[p4] . dm[nu]

Next invoke the FeynCalc functions

Interactions and Feynman Diagrams

You might also like