-- Verify formulas for electron-positron annihilation E = sqrt(p^2 + m^2) p1 = (E, 0, 0, p) p2 = (E, 0, 0, -p) p3 = (E, E sin(theta) cos(phi), E sin(theta) sin(phi), E cos(theta)) p4 = (E, -E sin(theta) cos(phi), -E sin(theta) sin(phi), -E cos(theta)) u11 = (E + m, 0, p1[4], p1[2] + i p1[3]) / sqrt(E + m) u12 = (0, E + m, p1[2] - i p1[3], -p1[4]) / sqrt(E + m) v21 = (p2[4], p2[2] + i p2[3], E + m, 0) / sqrt(E + m) v22 = (p2[2] - i p2[3], -p2[4], 0, E + m) / sqrt(E + m) I = ((1,0,0,0),(0,1,0,0),(0,0,1,0),(0,0,0,1)) gmunu = ((1,0,0,0),(0,-1,0,0),(0,0,-1,0),(0,0,0,-1)) gamma0 = ((1,0,0,0),(0,1,0,0),(0,0,-1,0),(0,0,0,-1)) gamma1 = ((0,0,0,1),(0,0,1,0),(0,-1,0,0),(-1,0,0,0)) gamma2 = ((0,0,0,-i),(0,0,i,0),(0,i,0,0),(-i,0,0,0)) gamma3 = ((0,0,1,0),(0,0,0,-1),(-1,0,0,0),(0,1,0,0)) gamma = (gamma0,gamma1,gamma2,gamma3) gammaT = transpose(gamma) gammaL = transpose(dot(gmunu,gamma)) q1 = p1 - p3 q2 = p1 - p4 qslash1 = dot(q1,gmunu,gamma) qslash2 = dot(q2,gmunu,gamma) "Verify Casimir trick" v21bar = dot(conj(v21),gamma0) -- adjoint of v21 v22bar = dot(conj(v22),gamma0) -- adjoint of v22 M111 = dot(v21bar, -i e gammaT, qslash1 + m I, -i e gammaT, u11) M112 = dot(v21bar, -i e gammaT, qslash1 + m I, -i e gammaT, u12) M121 = dot(v22bar, -i e gammaT, qslash1 + m I, -i e gammaT, u11) M122 = dot(v22bar, -i e gammaT, qslash1 + m I, -i e gammaT, u12) M211 = dot(v21bar, -i e gammaT, qslash2 + m I, -i e gammaT, u11) M212 = dot(v21bar, -i e gammaT, qslash2 + m I, -i e gammaT, u12) M221 = dot(v22bar, -i e gammaT, qslash2 + m I, -i e gammaT, u11) M222 = dot(v22bar, -i e gammaT, qslash2 + m I, -i e gammaT, u12) M11 = contract(dot(M111, gmunu, transpose(conj(M111)), gmunu)) + contract(dot(M112, gmunu, transpose(conj(M112)), gmunu)) + contract(dot(M121, gmunu, transpose(conj(M121)), gmunu)) + contract(dot(M122, gmunu, transpose(conj(M122)), gmunu)) M12 = contract(dot(M111, gmunu, conj(M211), gmunu)) + contract(dot(M112, gmunu, conj(M212), gmunu)) + contract(dot(M121, gmunu, conj(M221), gmunu)) + contract(dot(M122, gmunu, conj(M222), gmunu)) M21 = contract(dot(M211, gmunu, conj(M111), gmunu)) + contract(dot(M212, gmunu, conj(M112), gmunu)) + contract(dot(M221, gmunu, conj(M121), gmunu)) + contract(dot(M222, gmunu, conj(M122), gmunu)) M22 = contract(dot(M211, gmunu, transpose(conj(M211)), gmunu)) + contract(dot(M212, gmunu, transpose(conj(M212)), gmunu)) + contract(dot(M221, gmunu, transpose(conj(M221)), gmunu)) + contract(dot(M222, gmunu, transpose(conj(M222)), gmunu)) check(M12 == M21) pslash1 = dot(p1,gmunu,gamma) pslash2 = dot(p2,gmunu,gamma) P1 = pslash1 + m I P2 = pslash2 - m I Q1 = qslash1 + m I Q2 = qslash2 + m I T = dot(P1,gammaT,Q1,gammaT,P2,gammaL,Q1,gammaL) f11 = contract(T,3,4,2,3,1,2) T = dot(P1,gammaT,Q2,gammaT,P2,gammaL,Q1,gammaL) f12 = contract(T,3,5,2,3,1,2) T = dot(P1,gammaT,Q2,gammaT,P2,gammaL,Q2,gammaL) f22 = contract(T,3,4,2,3,1,2) check(e^4 f11 == M11) check(2 e^4 f12 == M12 + M21) check(e^4 f22 == M22) "ok" "Verify probability density" p12 = dot(p1,gmunu,p2) p13 = dot(p1,gmunu,p3) p14 = dot(p1,gmunu,p4) check(f11 == 32 p13 p14 - 32 m^2 p12 + 64 m^2 p13 + 32 m^2 p14 - 64 m^4) check(f12 == 16 m^2 p13 + 16 m^2 p14 - 32 m^4) check(f22 == 32 p13 p14 - 32 m^2 p12 + 32 m^2 p13 + 64 m^2 p14 - 64 m^4) s = dot(p1 + p2, gmunu, p1 + p2) t = dot(p1 - p3, gmunu, p1 - p3) u = dot(p1 - p4, gmunu, p1 - p4) check(f11 == 8 t u - 24 t m^2 - 8 u m^2 - 8 m^4) check(f12 == 8 s m^2 - 32 m^4) check(f22 == 8 t u - 8 t m^2 - 24 u m^2 - 8 m^4) d11 = (t - m^2)^2 d12 = (t - m^2) (u - m^2) d22 = (u - m^2)^2 f = e^4 (f11/d11 + 2 f12/d12 + f22/d22) / 4 -- high energy approximation m = 0 check(s == 4 E^2) check(t == -2 E^2 (1 - cos(theta))) check(u == -2 E^2 (1 + cos(theta))) a = 1 + cos(theta) b = 1 - cos(theta) check(f == 2 e^4 (a/b + b/a)) m = quote(m) -- verify integral a = 1 + cos(theta) b = 1 - cos(theta) f = a/b + b/a I = 2 cos(theta) + 2 log(1 - cos(theta)) - 2 log(1 + cos(theta)) check(f sin(theta) == d(I,theta)) "ok" Run