Configuration:  "ChPTVirtualPhotons2"

$LoadPhi = True ; $LoadFeynArts = True ;

$Configuration = "ChPTVirtualPhotons2" ; $Lagrangians = {"ChPTVirtualPhotons2"[2], "ChPTVirtualPhotons2"[4]} ;

Needs["HighEnergyPhysics`FeynCalc`"] ;

FeynCalc 4.2.0
For help, type ?FeynCalc,
use the built-in help system
or visit www.feyncalc.org

Loading PHI

Loading FeynArts

FeynArts 3.1

by Hagen Eck, Sepp Kueblbeck, and Thomas Hahn

last revised 12 Feb 03

patched for use with FeynCalc by Frederik Orellana

•Preliminaries

Cosmetics:

VariableBoxes["k"] ; VariableBoxes["τ"] ; VariableBoxes["ρ"] ;

Keep things compact:

DeclareUScalar[UTrace1]

SetOptions[#, Explicit -> False] & /@ {MM, SMM, UChi, FieldStrengthTensorFull, FieldStrengthTensor, CovariantFieldDerivative, LeftComponent, RightComponent} ;

SetOptions[SUNReduce, FullReduce -> True] ; SetOptions[UReduce, FullReduce -> True] ; SetOptions[IndicesCleanup, LorentzDummys -> {"μ", "μ", "μ", "μ", "μ"}] ;

IsoVector[QuantumField[Particle[LeftComponent[0] | RightComponent[0]], LorentzIndex[_]]][_] := 0 ;

QuantumField[Particle[lr : (LeftComponent[0] | RightComponent[0])], r___, SUNIndex[0]][x_] := QuantumField[Particle[lr], r][x] ;

UTrace1[FieldStrengthTensorFull[LorentzIndex[μ1_], QuantumField[Particle[LeftComponent[0] | RightComponent[0]], LorentzIndex[μ2_]][x_], x_, ___]] := 0 ;

$UMatrices = Union[$UMatrices, {LeftComponent[0], RightComponent[0]}] ;

We use k1 and k2 for the indices with 7 values - 3 SU(2) and 4 Lorentz values.

projectionRules = {HoldPattern[KroneckerDelta[LorentzIndex[μ_], SUNIndex[j_]] UMatrix[UGenerator[SUNIndex[j_]]]] -> 0, HoldPattern[KroneckerDelta[LorentzIndex[μ_], SUNIndex[j_]] NM[___, UMatrix[UGenerator[SUNIndex[j_]]], ___]] -> 0, HoldPattern[KroneckerDelta[LorentzIndex[μ_], SUNIndex[j_]] UTrace1[NM[___, UMatrix[UGenerator[SUNIndex[j_]]], ___]]] -> 0, UMatrix[UGenerator[LorentzIndex[_]]] -> 0, KroneckerDelta[LorentzIndex[_], SUNIndex[_]]^2 -> 4} ;

SetOptions[CovariantFieldDerivative, Explicit -> False] ;

test1 = CovariantNabla[HRight[x], x, {μ}] /. $Substitutions /. {GRight[LorentzIndex[mu_]] -> GRight[mu], GLeft[LorentzIndex[mu_]] -> GLeft[mu]} // NMExpand // Expand // UReduce

1/2 (∂ _ μ(öÆ)^† '6 Q _ R '6 öÆ) + öÆ^† '6 ∂ _ μ(Q _ R) '6 öÆ + 1/2 (öÆ^† '6 Q _ R '6 ∂ _ μ(öÆ)) + 1/2 (∂ _ μ(öÆ) '6 Q _ L '6 öÆ^†) + öÆ '6 ∂ _ μ(Q _ L) '6 öÆ^† + 1/2 (öÆ '6 Q _ L '6 ∂ _ μ(öÆ)^†) - 1/2 i (öÆ^† '6 G _ μ^R '6 Q _ R '6 öÆ) + 1/2 i (öÆ^† '6 Q _ R '6 G _ μ^R '6 öÆ) - 1/2 i (öÆ '6 G _ μ^L '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ '6 Q _ L '6 G _ μ^L '6 öÆ^†) - 1/2 (∂ _ μ(öÆ)^† '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) - 1/2 (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 ∂ _ μ(öÆ)^†) - 1/2 (∂ _ μ(öÆ) '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) - 1/2 (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 ∂ _ μ(öÆ)) - 1/2 i (öÆ^† '6 G _ μ^R '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 G _ μ^L '6 öÆ^†) - 1/2 i (öÆ '6 G _ μ^L '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) + 1/2 i (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 G _ μ^R '6 öÆ)

test2 = Collect[ I/2 UCommutator[USmall[μ][x], HLeft[x]] + NM[Adjoint[SMM[x]], CQRight[μ][x], SMM[x]] + NM[SMM[x], CQLeft[μ][x], Adjoint[SMM[x]]] /. $Substitutions /. CovariantFieldDerivative[a_, x_, mu_] :> (FieldDerivative[a, x, {μ}] - I NM[GRight[μ][x], a] + I NM[a, GLeft[μ][x]]) /. MM[x] -> NM[SMM[x], SMM[x]] // NMExpand // Expand // UReduce, _NM]

1/2 (∂ _ μ(öÆ)^† '6 Q _ R '6 öÆ) + öÆ^† '6 ∂ _ μ(Q _ R) '6 öÆ + 1/2 (öÆ^† '6 Q _ R '6 ∂ _ μ(öÆ)) + 1/2 (∂ _ μ(öÆ) '6 Q _ L '6 öÆ^†) + öÆ '6 ∂ _ μ(Q _ L) '6 öÆ^† + 1/2 (öÆ '6 Q _ L '6 ∂ _ μ(öÆ)^†) - 1/2 i (öÆ^† '6 G _ μ^R '6 Q _ R '6 öÆ) + 1/2 i (öÆ^† '6 Q _ R '6 G _ μ^R '6 öÆ) - 1/2 i (öÆ '6 G _ μ^L '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ '6 Q _ L '6 G _ μ^L '6 öÆ^†) - 1/2 (∂ _ μ(öÆ)^† '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) - 1/2 (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 ∂ _ μ(öÆ)^†) - 1/2 (∂ _ μ(öÆ) '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) - 1/2 (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 ∂ _ μ(öÆ)) - 1/2 i (öÆ^† '6 G _ μ^R '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 G _ μ^L '6 öÆ^†) - 1/2 i (öÆ '6 G _ μ^L '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) + 1/2 i (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 G _ μ^R '6 öÆ)

res = Collect[test1 - test2, _NM]

0

test1 = CovariantNabla[HLeft[x], x, {μ}] /. $Substitutions /. {GRight[LorentzIndex[mu_]] -> GRight[mu], GLeft[LorentzIndex[mu_]] -> GLeft[mu]} // NMExpand // Expand // UReduce

1/2 (∂ _ μ(öÆ)^† '6 Q _ R '6 öÆ) + öÆ^† '6 ∂ _ μ(Q _ R) '6 öÆ + 1/2 (öÆ^† '6 Q _ R '6 ∂ _ μ(öÆ)) - 1/2 (∂ _ μ(öÆ) '6 Q _ L '6 öÆ^†) - öÆ '6 ∂ _ μ(Q _ L) '6 öÆ^† - 1/2 (öÆ '6 Q _ L '6 ∂ _ μ(öÆ)^†) - 1/2 i (öÆ^† '6 G _ μ^R '6 Q _ R '6 öÆ) + 1/2 i (öÆ^† '6 Q _ R '6 G _ μ^R '6 öÆ) + 1/2 i (öÆ '6 G _ μ^L '6 Q _ L '6 öÆ^†) - 1/2 i (öÆ '6 Q _ L '6 G _ μ^L '6 öÆ^†) + 1/2 (∂ _ μ(öÆ)^† '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) - 1/2 (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 ∂ _ μ(öÆ)^†) - 1/2 (∂ _ μ(öÆ) '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) + 1/2 (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 ∂ _ μ(öÆ)) + 1/2 i (öÆ^† '6 G _ μ^R '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 G _ μ^L '6 öÆ^†) - 1/2 i (öÆ '6 G _ μ^L '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) - 1/2 i (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 G _ μ^R '6 öÆ)

test2 = Collect[ I/2 UCommutator[USmall[μ][x], HRight[x]] + NM[Adjoint[SMM[x]], CQRight[μ][x], SMM[x]] - NM[SMM[x], CQLeft[μ][x], Adjoint[SMM[x]]] /. $Substitutions /. CovariantFieldDerivative[a_, x_, mu_] :> (FieldDerivative[a, x, {μ}] - I NM[GRight[μ][x], a] + I NM[a, GLeft[μ][x]]) /. MM[x] -> NM[SMM[x], SMM[x]] // NMExpand // Expand // UReduce, _NM]

1/2 (∂ _ μ(öÆ)^† '6 Q _ R '6 öÆ) + öÆ^† '6 ∂ _ μ(Q _ R) '6 öÆ + 1/2 (öÆ^† '6 Q _ R '6 ∂ _ μ(öÆ)) - 1/2 (∂ _ μ(öÆ) '6 Q _ L '6 öÆ^†) - öÆ '6 ∂ _ μ(Q _ L) '6 öÆ^† - 1/2 (öÆ '6 Q _ L '6 ∂ _ μ(öÆ)^†) - 1/2 i (öÆ^† '6 G _ μ^R '6 Q _ R '6 öÆ) + 1/2 i (öÆ^† '6 Q _ R '6 G _ μ^R '6 öÆ) + 1/2 i (öÆ '6 G _ μ^L '6 Q _ L '6 öÆ^†) - 1/2 i (öÆ '6 Q _ L '6 G _ μ^L '6 öÆ^†) + 1/2 (∂ _ μ(öÆ)^† '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) - 1/2 (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 ∂ _ μ(öÆ)^†) - 1/2 (∂ _ μ(öÆ) '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) + 1/2 (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 ∂ _ μ(öÆ)) + 1/2 i (öÆ^† '6 G _ μ^R '6 öÆ '6 öÆ '6 Q _ L '6 öÆ^†) + 1/2 i (öÆ^† '6 Q _ R '6 öÆ '6 öÆ '6 G _ μ^L '6 öÆ^†) - 1/2 i (öÆ '6 G _ μ^L '6 öÆ^† '6 öÆ^† '6 Q _ R '6 öÆ) - 1/2 i (öÆ '6 Q _ L '6 öÆ^† '6 öÆ^† '6 G _ μ^R '6 öÆ)

res = Collect[test1 - test2, _NM]

0

Converted by Mathematica  (July 10, 2003)