/**
* Split-SUSY Higgs self-energy 3-loop contribution from a gluino,
* \f$O(\alpha_t \alpha_s^2)\f$. Taken from arxiv:1312.5220,
* Eq. (4.8).
*
* @note The result contains the 3-loop tadpole diagrams. It is
* therefore not 1-particle irreducible (1PI).
*
* @note The sign of the result is opposite to arxiv:1312.5220,
* Eq. (4.8).
*
* @warning The result is in Landau gauge (\f$\xi = 0\f$).
*
* @param scale renormalization scale
* @param mt MS-bar top mass
* @param yt MS-bar Yukawa coupling
* @param g3 MS-bar strong gauge coupling
* @param mg MS-bar gluino mass
*
* @return real part of 3-loop gluino contribution to Higgs self-energy
*/
double self_energy_higgs_3loop_gluino_split(
double scale, double mt, double yt, double g3, double mg)
{
const double yt2 = Sqr(yt);
const double mt2 = Sqr(mt);
const double g34 = std::pow(g3, 4);
const double mg2 = Sqr(mg);
const double Q2 = Sqr(scale);
const double LogG = FiniteLog(mg2 / Q2);
const double LogG3 = std::pow(LogG, 3);
const double result = 64*g34*yt2*mt2*LogG3;
return - result * threeLoop;
}
/**
* Standard Model Higgs mass 4-loop contribution, \f$O(\alpha_t
* \alpha_s^3)\f$. Taken from arxiv:1508.00912, Eq. (5.5).
*
* @note The result contains the 4-loop tadpole diagrams. It is
* therefore not 1-particle irreducible (1PI).
*
* @param scale renormalization scale
* @param mt MS-bar top mass
* @param yt MS-bar Yukawa coupling
* @param g3 MS-bar strong gauge coupling
*
* @return real part of 4-loop correction \f$O(\alpha_t \alpha_s^3)\f$
*/
double delta_mh_4loop_at_as_as_as_sm(
double scale, double mt, double yt, double g3)
{
const double yt2 = Sqr(yt);
const double mt2 = Sqr(mt);
const double g36 = Power6(g3);
const double Q2 = Sqr(scale);
const double LogT = FiniteLog(mt2 / Q2);
const double LogT2 = Sqr(LogT);
const double LogT3 = Power3(LogT);
const double LogT4 = Power4(LogT);
const double result =
g36*yt2*mt2 * (23925.974863118638
+ 6435.327201984095*LogT
- 20675.746270332107*LogT2
- 3456*LogT3
+ 5520*LogT4);
return result * fourLoop;
}
