•ϕϕϕϕ

The evaluated leading order Lagrangian:

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

Redundant terms are discarded:

$VeryVerbose = 2 ;

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

Using Method->Expand

Putting on dummy factors

Expanding NM products

Expanding DOT products

Expanding

Discarding terms

Applying CommutatorReduce

DeclareUScalar[QuantumField[aa___, Particle[bb__], cc___][x_]] ;

$VeryVerbose = 1 ;

lala = ($IsoIndicesCounter = 0 ; PhiToFC[IsoIndicesSupply[#]]) & /@ Expand[lll] ;

DeclareNonCommutative /@ {UMatrix[UGenerator[SUNIndex[j_], opt___], opts___], UMatrix[UGenerator[ExplicitSUNIndex[j_], opt___], opts___]} ;

Calculation of the Feynman rule:

fields = {QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I1]][p1], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I2]][p2], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I3]][p3], QuantumField[Particle[PhiMeson, RenormalizationState[0]], SUNIndex[I4]][p4]}

{ϕ^( )^I _ 1, ϕ^( )^I _ 2, ϕ^( )^I _ 3, ϕ^( )^I _ 4}

$VeryVerbose = 0 ;

lala // Length

19

ii = 0 ; res = (++ ii ; WriteString["stdout", ii, " "] ; FeynRule[#, fields]) & /@ lala ;

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

resu = Collect[res // IndicesCleanup, _UTrace1] ;

resul = Collect[resu, HoldPattern[Plus[__ ? ((! FreeQ[{##}, Momentum | ParticleMass, Infinity, Heads -> True]) &)]]] ;

Converted by Mathematica  (July 10, 2003)