•The third diagram

$VerticesSpecifications = {{VertexFields -> {PseudoScalar[0][0], PhiMeson[0]}, PhiModel -> ChPT3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PseudoScalar[0][0], PhiMeson[0], Scalar[1][0]}, PhiModel -> ChPTW3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PhiMeson[0], PhiMeson[0], Scalar[1][0]}, PhiModel -> ChPTW3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PseudoScalar[0][0], PhiMeson[0], PhiMeson[0], PhiMeson[0]}, PhiModel -> ChPT3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PhiMeson[0], PhiMeson[0], PhiMeson[0], PhiMeson[0]}, PhiModel -> ChPT3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PhiMeson[0], PhiMeson[0], PhiMeson[0], Scalar[1][0]}, PhiModel -> ChPTW3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PhiMeson[0], PhiMeson[0], PhiMeson[0], PhiMeson[0], Scalar[1][0]}, PhiModel -> ChPTW3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> Automatic}, {VertexFields -> {PseudoScalar[0][0], PhiMeson[0], PhiMeson[0], PhiMeson[0], Scalar[1][0]}, PhiModel -> ChPTW3, PerturbationOrder -> {2}, CouplingSign -> 1, XFileName -> {"ChPTW3P00P10P10P10o2"}}} ;

InitializeModel["Automatic", GenericModel -> "Automatic", Reinitialize -> True] ;

tmp2 = M$CouplingMatrices[[-1, 2, 1, 1]] /. Wrap -> Identity /. {I3 -> I2, I1 -> i3, I4 -> I1} // Simplify

-4/9 i c _ 5^(  ) (3 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 1 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 1 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 6 i d _ (I _ 1 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) - 3 f _ (6 I _ 2 k2)^(3) f _ (I _ 1 i _ 3 k1)^(3) d _ (I _ 2 k1 k2)^(3) - 3 f _ (6 I _ 1 k2)^(3) f _ (I _ 2 i _ 3 k1)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) + 3 i d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) - 3 d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 1 k1 k2)^(3) - 3 d _ (I _ 2 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) f _ (I _ 1 k1 k2)^(3) - 3 d _ (I _ 1 k1 k2)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 2 i _ 3 k1)^(3) + 3 i f _ (6 I _ 2 k1)^(3) f _ (I _ 1 k1 k2)^(3) f _ (I _ 2 i _ 3 k2)^(3) + 3 d _ (6 i _ 3 k1)^(3) (d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k2)^(3) - i f _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k2)^(3) + 2 d _ (I _ 1 I _ 2 k2)^(3) (d _ (I _ 2 k1 k2)^(3) - i f _ (I _ 2 k1 k2)^(3))) - 3 i d _ (6 I _ 1 k1)^(3) d _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 3 d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 1 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 3 d _ (I _ 1 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 6 d _ (I _ 1 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 i f _ (6 I _ 2 k1)^(3) f _ (I _ 1 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 d _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 i f _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 6 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 I _ 2 k1)^(3) δ _ (0  i _ 3) + 6 i d _ (I _ 2 I _ 2 k1)^(3) f _ (6 I _ 1 k1)^(3) δ _ (0  i _ 3) + 12 i d _ (I _ 1 I _ 2 k1)^(3) f _ (6 I _ 2 k1)^(3) δ _ (0  i _ 3) + 3 d _ (6 I _ 2 k1)^(3) (i d _ (I _ 2 k1 k2)^(3) f _ (I _ 1 i _ 3 k2)^(3) + f _ (I _ 2 k1 k2)^(3) f _ (I _ 1 i _ 3 k2)^(3) - i d _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 1 k1 k2)^(3) + d _ (I _ 1 k1 k2)^(3) (d _ (I _ 2 i _ 3 k2)^(3) + i f _ (I _ 2 i _ 3 k2)^(3)) + f _ (I _ 1 k1 k2)^(3) f _ (I _ 2 i _ 3 k2)^(3) + d _ (I _ 1 i _ 3 k2)^(3) (d _ (I _ 2 k1 k2)^(3) - i f _ (I _ 2 k1 k2)^(3)) + 4 d _ (I _ 1 I _ 2 k1)^(3) δ _ (0  i _ 3)) + 2 d _ (I _ 2 I _ 2 i _ 3)^(3) δ _ (6  I _ 1) + 4 d _ (I _ 1 I _ 2 i _ 3)^(3) δ _ (6  I _ 2) + 6 d _ (I _ 1 I _ 2 I _ 2)^(3) δ _ (6  i _ 3) + 4 d _ (6 I _ 2 i _ 3)^(3) δ _ (I _ 1  I _ 2) - 4 i f _ (6 I _ 2 i _ 3)^(3) δ _ (I _ 1  I _ 2) + 8 δ _ (0  i _ 3) δ _ (6  I _ 2) δ _ (I _ 1  I _ 2) + 2 d _ (6 I _ 2 I _ 2)^(3) δ _ (I _ 1  i _ 3) + 2 d _ (6 I _ 1 i _ 3)^(3) δ _ (I _ 2  I _ 2) - 2 i f _ (6 I _ 1 i _ 3)^(3) δ _ (I _ 2  I _ 2) + 4 δ _ (0  i _ 3) δ _ (6  I _ 1) δ _ (I _ 2  I _ 2) + 4 d _ (6 I _ 1 I _ 2)^(3) δ _ (I _ 2  i _ 3))

tmp1 = -16 amplFC[[3]] /. {_DecayConstant -> 1, _SumOver -> 1, _QuarkCondensate -> 1, _FeynAmpDenominator -> 1, _PropagatorDenominator -> 1, Pi -> 1, SUNDelta[SUNIndex[i1], SUNIndex[I1]] -> 1} // Simplify

-4/9 i c _ 5^(  ) (3 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 1 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 1 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 6 i d _ (I _ 1 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) - 3 f _ (6 I _ 2 k2)^(3) f _ (I _ 1 i _ 3 k1)^(3) d _ (I _ 2 k1 k2)^(3) - 3 f _ (6 I _ 1 k2)^(3) f _ (I _ 2 i _ 3 k1)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) d _ (I _ 2 k1 k2)^(3) + 3 i d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) + 3 i d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) - 3 d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 1 k1 k2)^(3) - 3 d _ (I _ 2 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) f _ (I _ 1 k1 k2)^(3) - 3 d _ (I _ 1 k1 k2)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 2 i _ 3 k1)^(3) + 3 i f _ (6 I _ 2 k1)^(3) f _ (I _ 1 k1 k2)^(3) f _ (I _ 2 i _ 3 k2)^(3) + 3 d _ (6 i _ 3 k1)^(3) (d _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k2)^(3) - i f _ (I _ 1 k1 k2)^(3) d _ (I _ 2 I _ 2 k2)^(3) + 2 d _ (I _ 1 I _ 2 k2)^(3) (d _ (I _ 2 k1 k2)^(3) - i f _ (I _ 2 k1 k2)^(3))) - 3 i d _ (6 I _ 1 k1)^(3) d _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 3 d _ (I _ 2 i _ 3 k1)^(3) f _ (6 I _ 1 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 3 d _ (I _ 1 i _ 3 k1)^(3) f _ (6 I _ 2 k2)^(3) f _ (I _ 2 k1 k2)^(3) - 6 d _ (I _ 1 I _ 2 k1)^(3) f _ (6 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 i f _ (6 I _ 2 k1)^(3) f _ (I _ 1 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 d _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 3 i f _ (6 I _ 1 k1)^(3) f _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 2 k1 k2)^(3) + 6 d _ (6 I _ 1 k1)^(3) d _ (I _ 2 I _ 2 k1)^(3) δ _ (0  i _ 3) + 6 i d _ (I _ 2 I _ 2 k1)^(3) f _ (6 I _ 1 k1)^(3) δ _ (0  i _ 3) + 12 i d _ (I _ 1 I _ 2 k1)^(3) f _ (6 I _ 2 k1)^(3) δ _ (0  i _ 3) + 3 d _ (6 I _ 2 k1)^(3) (i d _ (I _ 2 k1 k2)^(3) f _ (I _ 1 i _ 3 k2)^(3) + f _ (I _ 2 k1 k2)^(3) f _ (I _ 1 i _ 3 k2)^(3) - i d _ (I _ 2 i _ 3 k2)^(3) f _ (I _ 1 k1 k2)^(3) + d _ (I _ 1 k1 k2)^(3) (d _ (I _ 2 i _ 3 k2)^(3) + i f _ (I _ 2 i _ 3 k2)^(3)) + f _ (I _ 1 k1 k2)^(3) f _ (I _ 2 i _ 3 k2)^(3) + d _ (I _ 1 i _ 3 k2)^(3) (d _ (I _ 2 k1 k2)^(3) - i f _ (I _ 2 k1 k2)^(3)) + 4 d _ (I _ 1 I _ 2 k1)^(3) δ _ (0  i _ 3)) + 2 d _ (I _ 2 I _ 2 i _ 3)^(3) δ _ (6  I _ 1) + 4 d _ (I _ 1 I _ 2 i _ 3)^(3) δ _ (6  I _ 2) + 6 d _ (I _ 1 I _ 2 I _ 2)^(3) δ _ (6  i _ 3) + 4 d _ (6 I _ 2 i _ 3)^(3) δ _ (I _ 1  I _ 2) - 4 i f _ (6 I _ 2 i _ 3)^(3) δ _ (I _ 1  I _ 2) + 8 δ _ (0  i _ 3) δ _ (6  I _ 2) δ _ (I _ 1  I _ 2) + 2 d _ (6 I _ 2 I _ 2)^(3) δ _ (I _ 1  i _ 3) + 2 d _ (6 I _ 1 i _ 3)^(3) δ _ (I _ 2  I _ 2) - 2 i f _ (6 I _ 1 i _ 3)^(3) δ _ (I _ 2  I _ 2) + 4 δ _ (0  i _ 3) δ _ (6  I _ 1) δ _ (I _ 2  I _ 2) + 4 d _ (6 I _ 1 I _ 2)^(3) δ _ (I _ 2  i _ 3))

tmp1 - tmp2 /. SUNIndex | ExplicitSUNIndex -> Identity // Simplify

0

Table[tmp1 /. {I1 -> 6, i3 -> 3} // SUNReduce // SUNReduce // SUNReduce, {I2, 8}] // Simplify

{4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  )}

Table[M$CouplingMatrices[[-1, 2, 1, 1]] /. Wrap -> Identity /. {I3 -> ii, I2 -> ii, I1 -> i3, I4 -> I1} /. ((a : (SUNDelta | SU3Delta | SUNF | SU3F | SUND | SU3D))[b__]) :> a[Sequence @@ (SUNIndex /@ {b})] /. {I1 -> 6, i3 -> 3} // SUNReduce // SUNReduce // SUNReduce, {ii, 8}] // Simplify

{4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  )}

Table[tmp2 /. ((a : (SUNDelta | SU3Delta | SUNF | SU3F | SUND | SU3D))[b__]) :> a[Sequence @@ (SUNIndex /@ {b})] /. {I1 -> 6, i3 -> 3} // SUNReduce // SUNReduce // SUNReduce, {I2, 8}] // Simplify

{4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  )}

mel = Table[(SUNReduce[#, Explicit -> True, HoldSums -> False] & /@ (tmp1 /. {I1 -> 6, i3 -> 3} // Expand)) // SUNReduce // SUNReduce // Simplify, {I2, 8}]

{4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4 i c _ 5^(  ), 4/3 i c _ 5^(  ), 4/3 i c _ 5^(  )}

mel2 = OneLoop[p3, #] & /@ Join[FeynAmpDenominator[PropagatorDenominator[Momentum[p3], ParticleMass[Pion]]] * mel[[{1, 2, 3}]], FeynAmpDenominator[PropagatorDenominator[Momentum[p3], ParticleMass[Kaon]]] * mel[[{4, 5, 6, 7}]], FeynAmpDenominator[PropagatorDenominator[Momentum[p3], ParticleMass[EtaMeson]]] * mel[[{8}]]]

{-4/3 π^2 A _ 0  ( (m _ π^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ π^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ π^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ K^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ K^(ó    ))^2 ) c _ 5^(  ), -4 π^2 A _ 0  ( (m _ K^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ K^(ó    ))^2 ) c _ 5^(  ), -4/3 π^2 A _ 0  ( (m _ η^(ó    ))^2 ) c _ 5^(  )}

((Plus @@ mel2 // VeltmanExpand[#, ExplicitLeutwylerJ0 -> True] &) /. udrules /. gellmannOkubo /. _LeutwylerLambda -> 0 // Simplify) /. toEtaRules /. _RenormalizationState -> Sequence[] // Simplify

4/9 π^2 c _ 5^(  ) ((9 log((m _ π^(ó    ))^2/μ^2) - log((m _ η^(ó    ))^2/μ^2)) (m _ π^(ó    ))^2 + 2 (9 log((m _ K^(ó    ))^2/μ^2) + 2 log((m _ η^(ó    ))^2/μ^2)) (m _ K^(ó    ))^2)

(ampinfinities[[3]] /. _LeutwylerLambda -> 0 /. p3 -> -p1 /. MomentaRules /. gellmannOkubo // Simplify) /. toEtaRules

(c _ 5^(  ) ((log((m _ η^(ó    ))^2/μ^2) - 9 log((m _ π^(ó    ))^2/μ^2)) (m _ π^(ó    ))^2 - 2 (9 log((m _ K^(ó    ))^2/μ^2) + 2 log((m _ η^(ó    ))^2/μ^2)) (m _ K^(ó    ))^2) (!, _ 0^(  ))^2)/(36 π^2 (f _ ϕ^(ó    ))^2 (p _ 1^2 - (m _ K^(ó    ))^2))


Converted by Mathematica  (July 10, 2003)