I am trying to evaluate the Diractrace using Feyncalc package. The expression involves projection operator 'ab' which comes from completeness relation of spin 3/2 particles (Delta-). I am having some issues while evaluating this. The newer version of mathematica i.e. 14.3 does not provide any output and kernel becomes dead when I try to run hadronic tensor (last line of the code). However, The older version 14.1 works fine. Can anyone help me with this?
<< FeynCalc`
pd = Momentum[Subscript[p, Δ]];
pl = Momentum[Subscript[p, Λ]];
pe = Momentum[Subscript[p, e]];
pnu = Momentum[Subscript[p, ν]];
MD = Symbol["MD"];
me = Symbol["me"];
ML = Symbol["ML"];
f1 = Symbol["f1"];
g1 = Symbol["g1"];
GF = Symbol["GF"];
mu = LorentzIndex[μ];
nu = LorentzIndex[ν];
ab = MT[mu,
nu] - (1/3) GA[mu] .
GA[nu] + (1/(3 MD)) (GA[nu]*FV[pd, mu] -
GA[mu]*FV[pd, nu]) - (2/(3 MD^2)) FV[pd, mu] . FV[pd, nu]
hadronictensor = (1/4)*
DiracTrace[(GS[pd] + MD) . ab .
GA[mu] . (f1 + g1 GA[5]) . (GS[pl] + ML) . GA[nu] . (f1 + g1 GA[5]),
DiracTraceEvaluate -> True]