Hey Vladyslav,
Thanks for the reply, though the problem persists. I actually tried your solution before but I ran into the same problem.
after doing:
res=((M12 // FermionSpinSum) /. DiracTrace -> Tr) // Contract
I still get a product of traces, because there are still dummy indexes inside them:
tr(...GA[alpha]...GA[alphalin]..) tr (... GA[alpha]...GA[alphalin]...)
Everything is multiplying, so there are no "free indexes" M12 is a Lorentz scalar.
Best,
Duarte
This archive was generated by hypermail 2b29 : 09/04/20-12:55:05 AM Z CEST