Name: Michael Park Date: 12/19/16-05:48:18 AM Z


Hi. I have been having a lot of trouble trying to compute the scalar C0 integral for general masses using the Package-X routines.

Using the FeynCalc interface

in: PaXEvaluate[C0[0,0,mH^2,ma^2,ma^2,ma^2]]

out: log^2((-mH^2+Sqrt[mH^4-4 mH^2 MT3^2]+2 MT3^2)/(2 MT3^2))/(2 mH^2)

…however, when I try more general masses it seems to give up

in: PaXEvaluate[C0[0,0,mH^2,ma^2,mb^2,ma^2]]
out: C0[0,0,mH^2,ma^2,mb^2,ma^2]

Similarly, if I use the straight Package-X routines it gives me a reasonable answer for massless final states and one flavor in the loop

in: LoopIntegrate[Spur[\[DoubleStruckOne]],
     q, {q - k1, ma}, {q + k2, ma}, {q, ma}] /. {k1.k1 -> 0,
     k2.k2 -> 0, q.q -> mH^2, q -> k1 + k2,
     k1.k2 -> (mH^2 - 2 mW^2)/2} /. {LTensor[k1, \[Mu]] -> 0,
    LTensor[k2, \[Nu]] -> 0} // LoopRefine // Simplify
out: (2 Log[(2 ma^2 - mH^2 + 2 mW^2 +
   Sqrt[(mH^2 - 2 mW^2) (-4 ma^2 + mH^2 - 2 mW^2)])/(
  2 ma^2)]^2)/(mH^2 - 2 mW^2)

..however if I make the final state particles massive or if I enter more general loop masses, the computation seems to quit

I was wondering if there was any way to handle these more complicated C0 integrals? Thanks very much!