/** * 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; }