Date: 10/11/17-04:44:37 PM Z

I'm trying to calculate a diagram using three different approaches

1) OneLoop
2) TID -> OneLoop
3) OneLoopSimplify -> OneLoop

The first two give a correct result while the last one does not (Mathematica 11.1 and the latest stable version of FeynCalc).

Please note that the difference doesn't seem to be caused by the usual (D-4)/(D-4) terms. I guess that OneLoopSimplify treats chiral projections incorrectly. The example is appended below.

By the way, what is the ``official way" of calculating one-loop integrals, 1,2,3 or maybe none of them?

Best wishes,

In[2]:= InfParts := {B0[X__] -> FinB0[X] + 1/eps,
A0[Y_] -> FinA0[Y] + Y/eps};

R = (1 + GA[5])/2;
L = (1 - GA[5])/2;
pd[q_, m_] := PropagatorDenominator[Momentum[q, D], m];

In[9]:= DiagAnum = -I (SPD[k, k] MTD[mu, nu] -
mu].(Fal L + FalCC R).(GSD[k - p] + Mcc L + M R).GAD[
nu].(Fbe L + FbeCC R);

In[12]:= DiagA = ((4 Pi)^2/(2 Pi)^4) DiagAnum DiagAden ;

(*First way*)

In[16]:= ResA0 = OneLoop[k, DiagA];

In[17]:= ResA0 = ResA0 // DiracSimplify

In[18]:= ResA0 =
ResA0 /. InfParts /. D -> 4 - 2 eps // Series[#, {eps, 0, 0}] & //
Normal

(*Second way*)

In[30]:= ResA = OneLoopSimplify[DiagA, k];

In[20]:= ResA = OneLoop[k, ResA];

In[21]:= ResA = ResA // DiracSimplify;

In[22]:= ResA =
ResA /. InfParts /. D -> 4 - 2 eps // Series[#, {eps, 0, 0}] & //
Normal

(*Third way*)

In[23]:= ResA1 = TID[DiagA, k];

In[24]:= ResA1 = OneLoop[k, ResA1];

In[25]:= ResA1 = ResA1 // DiracSimplify;

In[26]:= ResA1 =
ResA1 /. InfParts /. D -> 4 - 2 eps // Series[#, {eps, 0, 0}] & //
Normal

(*Comparison*)

In[27]:= ResA1 - ResA0 // Simplify

Out[27]= 0

In[28]:= Limit[eps ( ResA1 - ResA), eps -> 0] // Simplify

Out[28]= -(3/
2) (FalCC M - Fal Mcc) (Fbe (-1 + DiracGamma[5]) +
FbeCC (1 + DiracGamma[5]))

This archive was generated by hypermail 2b29 : 02/16/19-08:40:02 AM Z CET