Name: Marco V Date: 04/19/18-12:54:01 AM Z


Hello,

I am using FeynCalc 9.2.0 and I am trying to compute the amplitude for a box diagram in QED, but when I use TID to decompose the loop integral, it seems that FC does not convert a denominator into a PaVe scalar function.

Setting the integrand

numBOX = e^4 (SpinorUBarD[q1,
       m\[Mu]].GAD[\[Alpha]].(GSD[p1 - k] +
        m\[Mu]).GAD[\[Beta]].SpinorUD[p1, m\[Mu]] SpinorUBarD[q2,
       me].GAD[\[Alpha]].(GSD[k + p2] + me).GAD[\[Beta]].SpinorUD[p2,
       me]);

integrandBOX =
  numBOX*FAD[{p1 - k, m\[Mu]}, {k + p2, me}, {k + p2 - q2,
     SmallVariable[\[Lambda]]}, {k, SmallVariable[\[Lambda]]}];

and putting it in TID

MBOX = ChangeDimension[(1/(16 Pi^4)) TID[integrandBOX, k,
     UsePaVeBasis -> True], 4] // Simplify

the result I obtain is mostly in terms of D coefficient functions and independent from the loop momentum k (as I expect), but there is also a term

-2 (me^2 + m\[Mu]^2 - s) Spinor[Momentum[q1], m\[Mu],
   1].GA[\[Alpha]].Spinor[Momentum[p1], m\[Mu], 1] Spinor[
   Momentum[q2], me, 1].GA[\[Alpha]].Spinor[Momentum[p2], me,
   1] FAD[{k, SmallVariable[\[Lambda]]}, {k + p1, m\[Mu]}, {k - p2,
   me}, {k - p2 + q2, SmallVariable[\[Lambda]]}, Dimension -> 4]

which I understand should correspond to a D0 but it is not transformed in it.

If I try to use TID on the last expression, it remains unchanged.
If instead I use OneLoop it gives me exactly the D0 which I was looking for.

I suppose the problem is related to the choice of dimension, but I don’t understand why the rest of the result is correctly in terms of PaVes. Can you explain what I am doing wrong?

Many Thanks