•ϕAS

IsoVector[QuantumField[___, Particle[Vector[0], ___], ___], ___][_] := 0 ;  QuantumField[___, Particle[Vector[0], ___], ___][_] := 0 ;  QuantumField[___, Particle[Scalar[2], ___], ___, ExplicitSUNIndex[_], ___][_] := 0 ;

The evaluated leading order Lagrangian:

LoadLagrangian[ChPTWS3[2]] ;

ll = ArgumentsSupply[Lagrangian[ChPTWS3[2]], x, RenormalizationState[0], ExpansionOrder -> 1, DropOrder -> 1, DiagonalToU -> True] ;

Redundant terms are discarded:

lll = DiscardTerms[ll, Retain -> {Particle[PhiMeson , RenormalizationState[0]] -> 1, Particle[AxialVector[0] , RenormalizationState[0]] -> 1, Particle[Scalar[1], RenormalizationState[0]] -> 1}, CommutatorReduce -> True, Method -> Expand] // Simplify

-(c _ 2^(  ) (< σ^6 '6 ∂ _ μ(Overscript[ϕ^( ), ->]) · Overscript[σ, ->] '6 Overscript[A^( ) _ μ, ->] · Overscript[σ, ->] > + < σ^6 '6 Overscript[A^( ) _ μ, ->] · Overscript[σ, ->] '6 ∂ _ μ(Overscript[ϕ^( ), ->]) · Overscript[σ, ->] >) s^( ))/f _ ϕ^(ó    )

Generator matrices are traced:

llle = ExpandU[ExpandU[lll, CommutatorReduce -> True], CommutatorReduce -> True] ;

Indices are supplied:

$IsoIndicesCounter = 0 ;

llll = llle // IsoIndicesSupply // IndicesCleanup ;

Calculation of the Feynman rule:

fields = {QuantumField[Particle[AxialVector[0], RenormalizationState[0]], LorentzIndex[μ1], SUNIndex[I1]][p1], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I2]][p2], QuantumField[Particle[Scalar[1], RenormalizationState[0]]][p3]}

{A^( ) _ μ _ 1^I _ 1, ϕ^( )^I _ 2, s^( )}

lal = Expand[llll] ;

melsimplified = If[Head[lal] == Plus, Plus @@ (IndicesCleanup[SUNReduce[FeynRule[#, fields], FullReduce -> True]] & /@ (List @@ lal))] ;

Another check that two different evaluations with specific components give the same result:

(SUNReduce[#, Explicit -> True, HoldSums -> False] & /@ (melsimplified /. {I1 -> 6, I2 -> 3} // Expand)) // Simplify

(2 c _ 2^(  ) p _ 2^μ _ 1)/f _ ϕ^(ó    )

(SUNReduce /@ SUNReduce /@ SUNReduce /@ SUNReduce /@ (melsimplified /. {I1 -> 6, I2 -> 3} // Expand)) // Simplify

(2 c _ 2^(  ) p _ 2^μ _ 1)/f _ ϕ^(ó    )


Converted by Mathematica  (July 10, 2003)