description = "Q Qbar -> Ga Gl, QCD, matrix element squared, tree";
If[ $FrontEnd === Null,
$FeynCalcStartupMessages = False;
Print[description];
];
If[ $Notebooks === False,
$FeynCalcStartupMessages = False
];
$LoadAddOns = {"FeynArts"};
<< FeynCalc`
$FAVerbose = 0;
FCCheckVersion[9, 3, 1];\text{FeynCalc }\;\text{10.0.0 (dev version, 2023-12-20 22:40:59 +01:00, dff3b835). For help, use the }\underline{\text{online} \;\text{documentation}}\;\text{, check out the }\underline{\text{wiki}}\;\text{ or visit the }\underline{\text{forum}.}
\text{Please check our }\underline{\text{FAQ}}\;\text{ for answers to some common FeynCalc questions and have a look at the supplied }\underline{\text{examples}.}
\text{If you use FeynCalc in your research, please evaluate FeynCalcHowToCite[] to learn how to cite this software.}
\text{Please keep in mind that the proper academic attribution of our work is crucial to ensure the future development of this package!}
\text{FeynArts }\;\text{3.11 (3 Aug 2020) patched for use with FeynCalc, for documentation see the }\underline{\text{manual}}\;\text{ or visit }\underline{\text{www}.\text{feynarts}.\text{de}.}
\text{If you use FeynArts in your research, please cite}
\text{ $\bullet $ T. Hahn, Comput. Phys. Commun., 140, 418-431, 2001, arXiv:hep-ph/0012260}
Nicer typesetting
MakeBoxes[p1, TraditionalForm] := "\!\(\*SubscriptBox[\(p\), \(1\)]\)";
MakeBoxes[p2, TraditionalForm] := "\!\(\*SubscriptBox[\(p\), \(2\)]\)";
MakeBoxes[k1, TraditionalForm] := "\!\(\*SubscriptBox[\(k\), \(1\)]\)";
MakeBoxes[k2, TraditionalForm] := "\!\(\*SubscriptBox[\(k\), \(2\)]\)";diags = InsertFields[CreateTopologies[0, 2 -> 2], {F[3, {1}], -F[3, {1}]} ->
{V[1], V[5]}, InsertionLevel -> {Classes}, Model -> "SMQCD"];
Paint[diags, ColumnsXRows -> {2, 1}, Numbering -> Simple,
SheetHeader -> None, ImageSize -> {512, 256}];
amp[0] = FCFAConvert[CreateFeynAmp[diags], IncomingMomenta -> {p1, p2},
OutgoingMomenta -> {k1, k2}, UndoChiralSplittings -> True, ChangeDimension -> 4,
TransversePolarizationVectors -> {k2}, List -> False, SMP -> True,
Contract -> True, DropSumOver -> True, Prefactor -> 3/2 SMP["e_Q"]]-\frac{\text{e} e_Q g_s T_{\text{Col2}\;\text{Col1}}^{\text{Glu4}} \left(\varphi (-\overline{p_2},m_u)\right).\left(\bar{\gamma }\cdot \bar{\varepsilon }^*\left(k_1\right)\right).\left(\bar{\gamma }\cdot \left(\overline{k_1}-\overline{p_2}\right)+m_u\right).\left(\bar{\gamma }\cdot \bar{\varepsilon }^*\left(k_2\right)\right).\left(\varphi (\overline{p_1},m_u)\right)}{(\overline{p_2}-\overline{k_1}){}^2-m_u^2}-\frac{\text{e} e_Q g_s T_{\text{Col2}\;\text{Col1}}^{\text{Glu4}} \left(\varphi (-\overline{p_2},m_u)\right).\left(\bar{\gamma }\cdot \bar{\varepsilon }^*\left(k_2\right)\right).\left(\bar{\gamma }\cdot \left(\overline{k_2}-\overline{p_2}\right)+m_u\right).\left(\bar{\gamma }\cdot \bar{\varepsilon }^*\left(k_1\right)\right).\left(\varphi (\overline{p_1},m_u)\right)}{(\overline{p_2}-\overline{k_2}){}^2-m_u^2}
FCClearScalarProducts[];
SetMandelstam[s, t, u, p1, p2, -k1, -k2, SMP["m_u"], SMP["m_u"], M, 0];Now come the usual steps, but with some special features. We do not average over the polarizations of the virtual photon but use the gauge trick for the sum over its polarizations. In this case the sum goes over all 4 unphysical polarizations, not just 2.
ampSquared[0] = 1/(SUNN^2) (amp[0] (ComplexConjugate[amp[0]])) //
FeynAmpDenominatorExplicit // SUNSimplify[#, Explicit -> True,
SUNNToCACF -> False] & // FermionSpinSum[#, ExtraFactor -> 1/2^2] & //
DiracSimplify // DoPolarizationSums[#, k1, 0, VirtualBoson -> True] & //
DoPolarizationSums[#, k2, k1] & // TrickMandelstam[#, {s, t, u, 2 SMP["m_u"]^2 + M^2}] & //
Simplify-\frac{\text{e}^2 \left(N^2-1\right) e_Q^2 g_s^2 \left(-m_u^4 \left(2 M^4-2 M^2 (t+u)+3 t^2+14 t u+3 u^2\right)+m_u^2 \left(2 M^4 (t+u)-8 M^2 t u+t^3+7 t^2 u+7 t u^2+u^3\right)+6 m_u^8-t u \left(2 M^4-2 M^2 (t+u)+t^2+u^2\right)\right)}{N^2 \left(u-m_u^2\right){}^2 \left(t-m_u^2\right){}^2}
ampSquaredMassless[0] = ampSquared[0] // ReplaceAll[#, {SMP["m_u"] -> 0}] & //
TrickMandelstam[#, {s, t, u, M^2}] &-\frac{\text{e}^2 \left(1-N^2\right) e_Q^2 g_s^2 \left(2 M^4-2 M^2 t-2 M^2 u+t^2+u^2\right)}{N^2 t u}
ampSquaredMasslessSUNN3[0] = ampSquaredMassless[0] /. SUNN -> 3\frac{8 \;\text{e}^2 e_Q^2 g_s^2 \left(2 M^4-2 M^2 t-2 M^2 u+t^2+u^2\right)}{9 t u}
```mathematica knownResults = { (8/9) SMP[“e”]^2 SMP[“e_Q”]^2 SMP[“g_s”]^2 (u/t + t/u + 2 M^2 (-t - u + M^2)/(t u)) }; FCCompareResults[{ampSquaredMasslessSUNN3[0]}, {knownResults}, Text -> {“with R. Field, Applications of Perturbative QCD, Eq 5.2.3:”, “CORRECT.”, “WRONG!”}, Interrupt -> {Hold[Quit[1]], Automatic}]; Print[“Time used:”, Round[N[TimeUsed[], 3], 0.001], ” s.”];
```mathematica
\text{$\backslash $tCheck with R. Field, Applications of Perturbative QCD, Eq 5.2.3:} \;\text{CORRECT.}
\text{$\backslash $tCPU Time used: }26.974\text{ s.}