High Precision Calculations of the Higgs Boson Mass

In this paper we are going to review the status of the computations of the perturbative quantum corrections to the Higgs boson mass in the Standard Model and in its supersymmetric extensions. In particular, supersymmetric theories require a very accurate computation of the Higgs boson mass, which includes corrections up to even three-loop level, since their predictions are limited by theoretical uncertainties. A discussion about these uncertainties in the context of the minimal and next to minimal supersymmetric standard models is included.


Introduction
The Standard Model (SM) of the electro-weak forces unification does not explain important phenomena of fundamental interactions like the neutrino oscillation, the dark matter, the baryon asymmetry, the vacuum stability of the model, the electro-weak hierarchy in the mass scales, among others. The explanation of the above phenomena remains without a conclusive solution and requires new physics, so called beyond the Standard Model (BSM). Moreover recent experimental results are moving away from their SM expectations. At the beginning of this year, the Fermilab Muon g-2 collaboration showed an increasingly discrepancy between the measurement of the muon anomalous magnetic momentum and the corresponding SM theoretical prediction. The last report shows a 4.2σ departure [1] from the reliable SM calculation, which are updated together with the inclusion of non-perturbative contributions in the so called white paper [2]. Another recent result carried out by the LHCb collaboration, the so-called R K anomaly which is also related to the muon, reported the measurement of the decay rate of B-mesons to a kaon and a pair of muons compared to the decay into a kaon and electrons, providing evidence of the violation of lepton-flavor universality at 3σ level [3]. Assuming those discrepancies are a consequence of new physical degrees of freedom, it is important to study their relation with the other SM problems and the possible consequences of the physics beyond the Standard Model. The agenda of new colliders, such as the high-luminosity (HL-) LHC [4], the International Linear Collider (ILC) [5,6], the Future Circular Collider (FCC-ee) [7,8,9], or the Circular Electron Positron Collider (CEPC) [10], includes the design of experiments with the necessary capabilities to measure with higher accuracy the properties of the precision observables in order to look for new physics signals coming from the results that deviate from the SM predictions. There will be also the necessity to interpret the coming experimental results and study the implications for specific UV models. This includes the development and improvement of precision calculation techniques. New tools for the evaluation of amplitudes at three, four and even five -loops together with the development of techniques for the tensorial reduction of multi-leg strong (QCD) and electroweak (EW) mixed processes are needed to study the Higgs properties, the EW physics, the top, the QCD and the Heavy Flavor physics, the possible BSM production and the global interpretation of all the expected results in the future colliders . A survey of the present status of analytic calculation methods of higher orders Feynman diagrams with zero or single mass scale is given in [11]. In particular, the recent discovery of a Higgs boson -like particle at the LHC [12] opened a window to the study of new phenomena. BSM models based on supersymmetry (SUSY) where the Higgs boson mass is a prediction of the theory, and their low-energy effective field theories (EFT), include additional neutral and charged particles in the Higgs sector with a mass above or below the energy scale where the LHC Higgs boson was found. The compatibility of this extended Higgs sector with the experimental results coming from the combined ATLAS/CMS measurements [13] will require a SM-like Higgs at 125.09 ± 0.21(stat.) ± 0.11(syst.) GeV , but also an agreement of the properties of the new particles with the present experimental bounds coming from the searches of dark matter candidates, the muon g − 2 anomaly, the LHCb results and the neutrino oscillations is mandatory. All of this phenomenological restrictions impose important constraints on the parameter space of the BSM models and can serve as a guide to look for the parameter region where the new particles might have appreciable interactions. In addition, it implies a very accurate computation of the Higgs boson mass since the central experimental value of 125.09 GeV is far to be accurate at tree-level in minimal supersymmetric models and besides the uncertainty associated to the theoretical predictions must reach the experimental accuracy at LHC and at future colliders, estimated to be of the order of 100 M eV [14]. The highprecision calculations of the radiative corrections to the Higgs boson mass is therefore of major importance to supersymmetric phenomenology. In this article we review the theoretical status of the Higgs boson mass in the SM and also in the main BSM, supersymmetric extensions. We start by introducing in Section 2 the Higgs sector of the SM at tree-level and then we review in Section 3 the Higgs sector in the context of softly broken supersymmetry. We focus our attention on the minimal supersymmetric extension of the SM (MSSM) with real parameters. The expressions of the tree-level Higgs boson masses and their leading quantum corrections in the above mentioned models together with a discussion about the conditions to guarantee the stability of the Higgs potential are also included in Section 2 and Section 3. In Section 4 a discussion about the theoretical uncertainties of the Higgs boson mass in supersymmetric models including some of the next to minimal extensions which enable to describe the recent experimental observations, as the one obtained in FERMILAB for the muon g − 2 anomaly, is included. Finally, in Section 5 we give a research outlook and conclude.

The Higgs Sector of the SM
In the SM, the Higgs sector contains just one scalar Higgs doublet which is sufficient to give masses to the three generation of quarks, the three uptype leptons and the EW gauge bosons Z 0 and W ± . Neutrinos instead remain massless particles because of chiral symmetry. The method to generate the masses is the famous Higgs mechanism of electroweak symmetry breaking [15,16] where the Higgs sector is introduced through the Lagrangian The first component of the Higgs doublet, φ + = φ 1 + iφ 2 , is a complex field that has electric charge Q = 1 and weak isospin third component T 3 = 1/2, while the component φ 0 = φ 3 + iφ 4 is a neutral complex field (Q = 0) and has weak isospin equals to T 3 = −1/2. As a consequence, both components have an where τ a are the Pauli's matrices, the parameters g and g are the interaction couplings associated to the SU (2) L and U (1) Y gauge symmetries respectively. The electroweak gauge bosons are defined as where the EW mixing angle, θ W , is related to the coupling constants g and g and the electric charge, e, by Additionally, the classical potential V (H) is restricted by the gauge invariance and renormalizability to have the functional form: The Higgs quartic coupling λ is a dimensionless positive number in order that the potential be bounded from below as H goes to infinity [17], while the quadratic parameter µ 2 can take either positive or negative values. In the case where µ 2 is positive, the classical potential has a minimum at H = 0, the SU (2) symmetry of the Lagrangian is manifest and µ 2 can be interpreted as the mass parameter of the field H. Note that the Lagrangian L H is also invariant under H → −H.
On the other hand, if µ 2 is negative the potential develops an infinite number of equivalent minima, determined by all the field configurations for which The new parameter v is called the vacuum expectation value (vev) of H and is defined up to a local gauge transformation under the gauge group SU (2) of weak isospin. The election of one of those minima breaks spontaneously the gauge symmetry, the vacuum is not invariant under the symmetry group of the Lagrangean. In order to determine the particle spectrum, we have to choose one of the possible minimum configurations and expand the field H around it. Due to the SU (2) symmetry, we can define a new real scalar field h(x) in the so called unitary gauge, such that The choice of φ 0 as the component of the Higgs doublet that gets a vacuum expectation value ensures the conservation of the electric charge. Besides, since only one component gets a vev when one chooses a direction, we have three broken continuous symmetries of the original SU (2) × U (1) continuous symmetry. As a consequence of the Goldstone theorem, the resulting spectrum contains three massless, zero-spin Goldstone bosons, that will become the zero spin component of the weak gauge bosons Z 0 and W ± . Expanding the classical potential around h = 0 we have: It is possible to associate a physical particle to the quantum field h(x) with a tree-level mass proportional to the self-interaction coupling λ and to the vev of H: This neutral zero-spin particle is well known as the SM Higgs boson and so far, its properties and all the SM predictions are in a great agreement with the Higgs boson properties measured at LHC by the ATLAS and CMS experiments. It is worth to mention that the mass of the Higgs boson in the SM is not a prediction of the theory, but a free input parameter whose numerical value must be specified in the theoretical computations of the other M h -dependent SM observables. From the kinetic term of the Lagrangian in eq. (2) and after expanding H around its vev, we have: Whereas the Lagrangian develops additional mass terms for the EW gauge bosons, the photon remains massless and therefore From the relation of the gauge boson mass M W with the Fermi coupling constant obtained from the muon decay process, where the Fermi constant G F ≈ 1.17 × 10 −5 GeV −2 and ∆r contains the radiative corrections to the relation, we can estimate the tree-level value for the vev Note that the mass generation for the EW bosons through the Higgs mechanism does not spoil the gauge invariance of the Lagrangian and consequently the theory is renormalizable, contrarily to what happens in the Proca theory, where the vector bosons mass is a given parameter of the bare Lagrangian [18]. On the other side, quark masses require the existence of left-handed SU (2) doublets and right-handed SU (2) singlets in such a way that where for quarks type up, u, and type down, d. The Yukawa parameters y u and y d are arbitrary, thus, the Higgs mechanism produces expressions for quark masses in terms of the Higgs vev, Similarly, masses of the leptons are generated from the interaction Lagrangian of the leptons with the Higgs field, where y e is the Yukawa coupling for the e-type lepton (e = e − , µ, τ ) L is an SU (2) doublet with components that contain a left-handed lepton, e L , and a left-handed neutrino, ν e L , and e R are right-handed leptons which are SU (2) singlets. The mass term for the lepton e is of the form: Due to the absence of a right-handed neutrino ν e R , there is no mass term for ν e . Experimentally, only left-handed neutrinos have been observed with masses that are very small but non-zero. The masses of the quarks and up-type leptons are not predicted in the model and therefore, just like the Higgs mass, they are input (running) parameters that must be measured at the experiments. The experimental measured values are related to the pole masses by the spectral representation of the full propagator. For the Higgs boson, the renormalized mass in the MS scheme has been computed to two-loop order including the mixed QCD/EW contributions to M h with the explicit inclusion of the tadpoles in the Higgs self-energy diagrams [19], also including the O(αα s ) corrections to the relation between the MS and pole masses of the Higgs boson. The radiative corrections in the gauge-less limit approximation at two-loop order is presented in [20], while the full two-loop correction in a hybrid MS/on-shell scheme can be found in [21]. In addition, the Higgs boson mass in the tadpole-free pure MS scheme including all the one-loop and two-loop contributions can be consulted in [22], where also the leading three-loop contributions estimated using the 1PI effective potential approach in the limit where g 2 , g 2 , λ g 2 s , y 2 t is included. The latter corrections together with the relations between the MS parameters and the on-shell masses of the SM are implemented in the public code SMDR [23]. With the help of SMDR the state of art of the Higgs boson mass in the SM can be easily visualized. Figure 1 shows the Higgs boson pole mass in the MS scheme as a function of the renormalization scale, Q, obtained from the SMDR code. The plot includes the full one-loop contribution (dashed black curve), full two-loop contribution (dotted black line) and the leading three-loop corrections (dot-dashed red line) at order O(g 4 s y 2 t )+O(g 2 s y 4 t )+O(y 6 t ), where g s is the strong coupling constant. The curves were drawn using the default input values in SMDR of the MS parameters at the initial scale Q 0 = 173.2 GeV . The running parameters of the model at Q are determined by the renormalization group evolution from the subtraction point defined at Q 0 . Note that the renormalization scale dependence makes the difference between the three loop and tree level result of the order of several tenths of MeV in Figure 1 when Q is close to 160 GeV . Nevertheless, as was argued in [22], this error can increase if we are far away of this scale. We must take into account that the three-loop correction to M h included in SMDR is in the 1PI effective potential approximation where the external momentum of the Higgs self-energies vanishes. The full determination of the complex pole squared mass requires, however, the inclusion of the external momentum dependence of the Higgs self energy. A complete Feynman diagrammatic computation of the three-loop Higgs self-energies is expected to make flat the Q dependence. Besides, even though the effective potential approximation might be expected to be the bulk of the three-loop correction, the complete two loop order with leading three-loop corrections of the pole mass in [22] may receive significant corrections because the ratio M h /M t ≈ 0.6 is not a really small expansion parameter. An additional important issue has to do with the stability of the SM vacuum. The LHC value of M h implies a near critical metastable vacuum [20,21], the metastability is the preferred option with 99.3% of probability, which means that the SM is at the border between stability and metastability regions and could be sitting on the stability region, i.e. it could reach and even cross the stability line, when its experimental and theoretical uncertainties have been taken into account. Including the two-loop threshold corrections of all SM parameters, the estimated overall theory error on M h is about 1.0 GeV . Although the dominant parametric uncertainties come from the top-quark mass and the QCD coupling, the limitations of the knowledge about the state of the SM vacuum should rely mainly on experimental errors. An unambiguously theoretical computation of M h should, to the possible extent, contains a theoretical error which reach at least the current experimental uncertainty. It is therefore important to perform the computation of the Higgs boson mass and its relation with the SM parameters as accurately as possible. Additionally, the information from the SM masses can be inverted to obtain the evolution of the interaction couplings constants. This can be done with the code SMDR as well, but there is also an alternative code named mr [24] which is based on the calculations reported in [25] that contains the relationships between the SM parameters of the on-shell scheme and their counterparts in the MS scheme at full two-loop order. It is worth mentioning that in mr the definition and renormalization of the vev is different from the used in SMDR. In mr, the renormalized vacuum, v tree , is defined as the minimum of the tree-level Higgs potential. The advantage of a scheme with this vacuum definition is that all the physical observables depending on it are gauge invariant due to the explicit insertion of the tadpole diagrams. Their insertion is done not only in the diagrams with counterterms but also in the 1PI Feynman diagrams as required to have the gauge independence of renormalized Green functions according to the gauge Ward identities. However, there is a disadvantage, the inclusion of tadpole diagrams involving the Higgs field include negative powers of the Higgs mass leading to a parametrically slower convergence. The leading small-M h behaviour is dominated by the one-loop correction, the threshold corrections to M h and λ contain terms that diverge as 1/M 2 h for M h → 0 [25,26]. In SMDR instead, the vev is defined to be the minimum of the full effective potential calculated in the Landau gauge. By this definition, the sum of all Higgs tadpole graphs, including the tree-level Higgs tadpole, vanishes identically, and therefore negative powers of λ and huge EW corrections are absent in the perturbative expansions of the pole masses and their relations with the MS parameters. This vev is in some sense a more faithful description of the true vacuum state. However, there is also a disadvantage. If the scheme is defined in terms of the self-energy diagrams without including the tadpole contribution, it gives rise to radiative corrections that are dependent of the gauge-fixing condition [27]. As a consequence, in this framework the MS masses are gauge-dependent quantities. Nevertheless, a MS mass is not a physical quantity nor a parameter of the bare Lagrangian, moreover the physical observable are independent on the running of the renormalized parameters (up to the error of the truncation of the perturbative series) and therefore the requirement of gauge-invariance for the running of the MS mass is not mandatory. The precise relationship between the threshold corrections that relate the MS masses to the parameters in the pole scheme with the three level vacuum v tree and the alternative computation implemented in SMDR deserve more attention for future studies.

Bosons
Fermions In this table only one generation of leptons and quarks is exhibited. For each lepton and quark super-multiplet there is associated a corresponding multiplet of scalar superpartners, while for the Higgs super-multiplet an anti-particle multiplet of chargeconjugated fermions is associated.

Higgs Sector in Supersymmetry
We start by reviewing the Higgs sector of the minimal supersymmetric extension of the SM (MSSM) at tree-level. Supersymmetry is the only fermionic generalization of the Poincaré symmetry of the space-time which transforms fermions into bosons and vice versa. The existence of such a non-trivial extension of the Poincaré symmetry is highly constrained by theoretical principles [28,29]. The MSSM contains the fields of the two-Higgs-doublet extension of the SM (2HDM) and their corresponding superpartners. The field content of the MSSM is summarized in Table 1 [30]. Specifically, the spectrum contains three gauge supermultiplets which consist of the gluons with their corresponding gluino fermionic superpartners and the EW gauge bosons with their gaugino superpartners. The matter supermultiplets consist of three generations of left-handed quarks and their scalar superpartners, the squarks, of three generations of left-handed leptons and their associated sleptons, and the corresponding antiparticles of all the fermions. The Higgs supermultiplets consists of two complex Higgs doublets, their higgsino fermionic superpartners and the corresponding antiparticles. In order to guarantee the cancellation of anomalies coming from the introduction of the higgsino superpartners and preserve supersymmetry and gauge invariance, the Higgs sector of the MSSM Lagrangian [31,32,33,34,35,36,37,38] requires the definition of at least two doublets with associated hypercharges Y 1 = −1 and Y 2 = +1, respectively. The vacuum expectation values v 1 and v 2 have a phase difference which is parametrized through the angle ψ 2 . In the real version of the MSSM (rMSSM) this angle is removed through a hypercharge rotation. We assume in the following that ψ 2 = 0. The second Higgs doublet is necessary, moreover, to generate mass for both up-and down-type quarks as well as charged leptons. The complex fields H 0 j have neutral components which are vevless scalar fields, H 0 j = φ 0 j + iχ 0 j , coming from the expansion around the minima of the Higgs effective potential, v 1,2 . The components H ± 1,2 are charged scalar fields with vanishing vacuum expectation values (vevs) to ensure that the U (1) em symmetry will not be spontaneously broken. The Higgs sector comes from the bare Lagrangian The kinetic term contains the covariant derivative where g , g and g s are the U(1) Y , SU(2) L and SU (3) gauge couplings respectively. The SU (2) generators Pauli matrices are represented by τ a while λ α stands for the SU (3) Gell-Mann matrices. The classical Higgs potential, The quadratic terms proportional to the higgsino mass parameter |µ| 2 in eq. (23) come from the F -contribution to the SUSY Lagrangian, while the terms with the EW gauge couplings (g, g ) come from the D-contribution; thus, from the D-terms, one obtains the four scalar vertices, which include the quartic Higgs self-interaction terms in the Higgs potential. In contrast to the SM, the Higgs quartic self-coupling is not a free parameter but is determined by the coupling (g 2 + g 2 )/8. From the F -terms, one obtains also another four scalar vertices, but they do not contribute to the Higgs potential since contain always at least two sfermions. The F and D fields are auxiliary fields that do not propagate in space-time and can be eliminated by applying the equations of motion. As supersymmetric particles have not been observed at the electroweak scale, supersymmetry is manifestly not an exact symmetry of the nature and therefore must be broken. Several SUSY breaking mechanisms have been considered [39,40,41] but in fact there is no consensus on how to break SUSY. In phenomenological application, one must introduce new terms that explicitly breaks supersymmetry and represent the low-energy effects of the unknown breaking mechanism. In eq. (23) the parameters m 2 H1 , m 2 H2 and b parameterize the soft SUSY breaking. The SUSY-invariant Lagrangian cannot accommodate electroweak symmetry breaking since the terms proportional to |µ| 2 are positive. Thus, the SUSY-violating parameters m 2 H1 and m 2 H2 , which can of course have either sign, are needed in order to break the EW symmetry. The b-term is the only that depends on the phases of the fields. The parameter b is real and positive and the fields H 0 1 and H 0 2 have equal and opposite phases which can be reduced both to zero through a U(1) Y gauge transformation since these fields have equal and opposite hypercharges. The vevs, v 1,2 , as well as the couplings are therefore all real, which means that CP invariance is not spontaneously broken by the classical Higgs potential of the rMSSM. The softness of the SUSY breaking terms means that they must be superrenormalizable, with mass dimension less than four, and therefore their couplings must have positive mass dimension. The reason is that a soft breaking term will not introduce additional divergences into the dimensionless coupling constants which guarantee the stability of the mass hierarchy, the cancellation of quadratically divergent corrections to scalar masses are maintained to all perturbative orders. Besides, it is important to emphasize that the mass terms which break SUSY and therefore the masses of the undiscovered SUSY particles do respect the SM gauge symmetries. The masses of the known SM particles all arise from the spontaneous breaking of the EW symmetry. From eq. (23) we can derive the linear part of the Higgs potential in the basis (φ 0 j , χ 0 j , H ± j ) where j = 1, 2. For the fields χ 0 j and H ± j there are no contributions since the rMSSM Higgs potential is invariant under CP-tranformation. Thus, the linear term of the effective potential is T 1 φ 0 1 + T 2 φ 0 2 , where the coefficients T j , better-known as Higgs tadpoles, have the expressions As the vevs v 1 and v 2 minimize the Higgs potential, the following stationary conditions are satisfied: From eq. (25) one can identify the conditions required for the stable minimum of V . First note that along the direction |H 0 1 | = |H 0 2 | the potential will be unbounded from below and therefore it does not have a minimum unless For |µ| 2 + m 2 H1 and |µ| 2 + m 2 H2 positive, the origin is not a minimum but a saddle point, and the minimum occur at non-zero vevs of H 0 j , when The rMSSM Higgs potential develops a stable minimum if the conditions of the eq. (26) and eq. (27) are met.
Turning to the bilinear part of the kinetic terms in the Lagrangian L V −H (eq. 21) the masses of the EW gauge bosons amount to The Z boson mass determines the tree-level relation v 2 1 + v 2 2 ≈ 174 GeV . The mass matrices of the rMSSM Higgs bosons can be identified from the bilinear part of the classical potential, To derive the M φ 0 -matrix the relations from eqs. (24) and (25) and the definitions are required. The tree-level mass matrix of the neutral φ 0 1,2 -bosons reads We have used the short notation s β = sin(β) and c β = cos(β) and we have written explicitly the contributions of the Higgs tadpoles, which vanish at treelevel according to eq. (25), because they will receive non-zero loop contributions when renormalization of the Higgs masses will be considered. Because both v 1 and v 2 are real and positive, the upper and lower bound on the angle β lies on the interval shown in eq. (29). By other side, the tree-level mass matrices of the χ 0 1,2 and H ± 1,2 components are given by and The potential can be brought into the physical basis, where the quadratic terms in the components of H j are diagonalized, through the rotations via the orthogonal transformation Thus, after diagonalization, besides of the EW boson masses of eq.  Higgs bosons, m ± , are linearly dependent on m A , they are usually used in the complex version of the MSSM. At tree-level we have: The tree-level masses of the CP-even Higgs boson particles, h and H, follow as predictions In most of the relevant phenomenology benchmark scenarios for MSSM Higgs boson searches, the LHC Higgs boson corresponds to the lightest CP-even Higgs boson with a mass m h which is not a free input parameter but it is predicted in the MSSM. From the mass formulas (35) and (36) the next important inequalities can be derived: As a consequence, the lightest Higgs boson mass is predicted to be bounded from above by the Z-boson mass, m h ≤ 91.2 GeV , modulo radiative corrections. This bound follows from the fact that the quartic coupling of the Higgs boson fields is determined by the size of the EW gauge couplings. The tree-level bound on m h has already been excluded by the experimental value found at the LHC. However, this tree-level prediction is strongly modified by higher-order quantum corrections making the MSSM compatible with the measured Higgs mass of 125 GeV and consistent with the similarities of the measured Higgs couplings to those in the SM [12]. The state of art of the quantum corrections to the lightest Higgs boson mass in the MSSM is really advanced and widely studied. The reported results include fixed-order diagrammatic calculations, resummation of the large logarithms in effective field theories (EFT) and hybrid fixed-order and EFT calculations up to three-loop level. All of these computations have been implemented in several public codes that provide a precise numerical evaluation of the Higgs masses in SUSY models. The relevant references and a comprehensive overview of the current status of Higgs-mass calculations in supersymmetry can be found in [42] and references therein. For the rMSSM, the dominant contributions come from the top-stop Yukawa corrections at zero external momentum in the EW gaugeless limit, which are proportional to the fourth power of the top mass, M t , and to the logarithm of the physical stop masses,m t1 and m t2 : The source of this correction is the soft breaking of supersymmetry which produces an incomplete cancellation between virtual top and stop loops. In this approximation the one-loop prediction of the Higgs boson masses can be expressed in the simple compact form: The Higgs mass can be also affected by a potentially large stop mixing X t due to the non-leading effects of the one-loop correction: where the function f is defined as: The dependence of the CP-even and charged Higgs boson masses on the parameter m A in eq. 39, including the dominant one-loop radiative corrections in a scenario where the squarks masses are put at the same supersymmetric scale M SU SY , produces the curves in Figure 2. Note that, independently of tanβ, the masses of the heavy Higgs particles, m H and m ± , grow linearly without boundary when the scale m A grows and have approximately the same magnitude, m H ≈ m ± . The lightest Higgs boson mass instead approaches to an asymptotic value, showing a more regular dependence on m A for large tanβ, where m h m A for small m A while m h const for large A-mass. The non-leading effects of the stop mixing X t allow a simple determination of an upper bound on the lightest Higgs boson mass at one-loop level, yielding the expression:  to higher-order SQCD effects.  For high SUSY energy scales, a resummation of the large logarithms with the help of the renormalization group techniques is required as was discussed in [46]. Once the resummation is performed, the dependence on M SU SY can be evaluated up to super-high energy scales of the order of the GUT (SUSY Grand Unified Theories) energy. In Figure 4 -Right we further explore the dependence of the Higgs boson mass on the SUSY scale and the stop mixing parameter. We have plotted the region of parameters which is compatible with a Higgs boson mass equal to 125.09 ± 5.0 GeV (brown region) for large values of M SU SY between 10 7 GeV M SU SY 10 16 GeV and for tanβ = 1. As was argued also in [46], it is possible to fix upper bounds on the SUSY scale M SU SY for scenarios with large values of tanβ (tan β > 10) if we impose that the m h -predictions in the MSSM is compatible with the combined CMS/ATLAS measured Higgs boson mass within the actual combined uncertainties. However, for small values of tanβ (tanβ < 10) such upper bounds coming from the LHC measurements cannot be derived and a strong dependence between M SU SY and X t is observed. In fact, for large values of the stop mixing parameter, |X t | 3M SU SY , the m h -prediction is compatible with the LHC measurement for M SU SY 10 7 GeV . It is possible, however, to derive upper bounds over M SU SY for small tanβ values if we impose the positivity of the running Higgs self-coupling λ(Q) from the EW scale until the next threshold, Q = M SU SY , which is equivalent to impose the stability of the SM renormalization group improved Higgs effective potential, up to the scale where supersymmetry is supposed to appear. The SM is considered here as a low-energy EFT of the rMSSM valid in the limit where we integrate out all SUSY QCD particles at the scale M SU SY , and therefore the RGEs correspond to the three-loop SM beta functions while the threshold corrections to the boundary condition for the Higgs quartic coupling contain two-loop corrections including the MSSM particles. The full two-loop threshold corrections to λ can be found in [47] and references there in. In the gray region of  [45,48]. However, such a region of superhigh SUSY energy scales cannot be tested at the current LHC or dark matter underground experiments and at the future colliders. An experimental test requires, therefore, the design of experiments capable to provide indirect signals of SUSY particles belonging to models for high energies of the order of the GUT scale, as is the case of supergravity. This has brought a special interest of the study of supersymmetry in the context of astroparticle physics and cosmology. The cosmological consequences and the possible applications to SUSY models can be found in [49,50,51,52,53]. It is worth mentioning that this super-high energies are not the only target where we could expect some SUSY evidence. There are recent phenomenological analysis that show the possibility to reproduce the 4.2σ discrepancy in the muon (g-2) anomaly together with the full Dark Matter (DM) relic density of the universe by taking the lightest supersymmetric particle (LSP) to be the DM candidate [54]. An upper limit of ∼ 600 GeV is obtained in this work for the LSP and next-to (N)LSP masses establishing clear search targets for the future HL-LHC EW searches and for future high-energy e + e − colliders. This bound is in accordance with the LHC exclusion bounds found for searches of direct production of charginos, neutralinos and sleptons [55,56,57,58]. We would like to remark that for the phenomenological MSSM (pMSSM) parameter ranges, where the SUSY masses are of the order of 1 to 10 TeV, the Higgs coupling measurements with the accuracies obtained on the LHC Run 2 data and those expected for the HL-LHC and future e + e − colliders can exclude up to 20% of the accepted pMSSM points, as was recently shown in [59]. Therefore, a large part of the pMSSM parameters are still to be probed, suggesting the conclusion that the properties of the observed Higgs boson are also MSSM-like. Finally, we stress that there are also well-motivated non-supersymmetric theories that naturally accommodate the experimental observation of a light 125 GeV Higgs boson and exhibit a large separation of scales able to match experimental observations, as is the case with the composite Higgs models [60,61,62] or the low-energy effective Higgs theories, such as the two Higgs doublet model (2HDM) [63,64].

Theoretical Uncertainties and Next to Minimal Supersymmetry
The computation of the Higgs boson mass in SUSY models requires to impose several approximations, including, in particular, the truncation of the involved perturbative expansions at some loop level, the limit where the effects of some subset of couplings within the considered loop-order are neglected or the approximation where some of the kinematics invariants vanish, as in the case of the effective potential approach where the limit of vanishing external momentum is adopted. Therefore, an estimation of the theoretical uncertainty, which measures the possible effects of the missing terms, must be included in a realistic computation of m h . For the case of the MSSM, the Higgs boson mass has been computed at fixed-order by using a Feynman diagrammatic approach, but also with the help of the renormalization group equations in an effective field theory (EFT) approach. Due to the well-known large logarithmic effects in the fixed-order calculation, an EFT computation where the large logarithms are resummed is considered the best option for high SUSY scales larger than about 1 T eV [46]. In fact, pure fixed-order calculations are more reliable when the SUSY masses are close to the EW scale, while pure EFT calculations are more reliable in heavy-SUSY scenarios. As a consequence, hybrid calculations where the fixed-order and the EFT results are combined in order to improve the theory uncertainty at each kinematical point have been implemented by different groups. Currently there are public codes which allows a numerical estimation of the different theoretical errors associated to the EFT, fixed-order and hybrid calculations of the Higgs boson mass. For an extensive and detailed explanation of these estimations, see [42] and references therein. Nonetheless, in the following we highlight some of the uncertainty sources which will require in the future improved perturbative computations. In the diagrammatic approach, a full two-loop fixed-order calculation of m h , of the lightest Higgs boson is still missing. The two-loop computation in the SQCD sector is complete, but the effects due to the EW gauge couplings are only available in the limit of vanishing external momentum. A complete twoloop calculation of scalar tadpoles and self-energies of a general renormalizable theory was presented recently in [65]; however, the relations between running parameters and on-shell observables require also the full two-loop gauge boson self-energies which are still missing. At three-loop level, contributions exist in the effective potential approach at order O(α t α 2 s ) while the calculation of the three-loop corrections that involve lower powers of the strong gauge coupling, e.g. O(α 2 t α s ), O(α 3 t ), contains only the logarithmic effects. On the other hand, in the EFT calculation of m h , beyond next to leading logarithmic (NLL) order, the estimation has so far been performed only for the simplest heavy-SUSY scenario where the MSSM is matched directly to the SM and all the SUSY particles are putted at the same scale M SU SY in the limit M SU SY M t . The next to next leading logarithms (NNLL) corrections neglect most of the effects that involve the EW gauge couplings, while at NNNLL they account only for the effects that involve the top Yukawa coupling combined with the highest powers of the strong gauge coupling. Taking into consideration these sources of uncertainties together with the additional theoretical and parametric uncertainties described in references [42,66,67,68], the codes FeynHiggs and FlexibleEFTHiggs allow for an estimation of the Higgs boson mass uncertainty as a function of the SUSY scale. In Figure 5 an estimation of the theory uncertainties in the hybrid approach of FeynHiggs as a function of the SUSY scale is drawn. The calculation was done in a scenario where all the squark masses (m q ), the gluino mass (m g ), the parameter µ and the CP-odd Higgs boson mass m A are set equal to a common scale M SU SY which is varied from 0.4 T eV to 100 T eV . The stop mixing parameter is taken to be X t = −2.45M SU SY while tanβ = 20. Besides, the squark masses and X t are running parameters evaluated at the scale Q = M SU SY in the DR renormalization scheme [69]. Assuming this MSSM scenario, the total theoretical uncertainty is in the range ∆m h ≈ 1 − 5 GeV , where, in particular, for large values of M SU SY the SUSY scale dependence is mild and the uncertainty decreases up to about 1.5 GeV . The size of ∆m h is significantly large compared with the current experimental error estimated at LHC (∼ 200 M eV ) which is about an order of magnitude less than the theoretical uncertainties estimated in Figure 5. Therefore, the theory uncertainty must be reduced by about a factor of 10 in order to reach the current and future experimental precision. This aim can be achieved by improving the accuracy in the measurement of the EW observables at the future colliders, but also through the explicit computations of the dominant missing terms in the theoretical prediction of m h , which additionally could tell us whether the current estimation of ∆m h is too optimistic or too pessimistic. If the accuracy of the Higgs mass calculation is improved, the existing uncertainty estimate must be adapted so that it simulates the dominant terms among those which have not been yet computed. Any future improvement in the accuracy of the fixed-order and EFT calculations of m h would also eventually improve the uncertainty in the hybrid calculations. There are several strategies to make an improvement. We just refer to one of them in this article which is related to a work in progress that our group is currently doing. The three-loop computation of the Higgs mass, in the SM as well as in the MSSM, was done in the effective potential approach which is equivalent to compute the self-energy corrections in the limit of vanishing external momentum. However, the inclusion of the external-momentum effects in the three-loop self-energies is actually attainable by numerical methods, based on the dispersion relations, implemented in TVID2 [70], that allows the evaluation of three-loop planar self-energies with arbitrary masses and external momentum. In addition, with the results presented recently in [71], where the differential equations method is used to compute numerically the renormalized -finite master integrals for arbitrary external momentum invariant and, in principle, for arbitrary masses, the evaluation of non-planar self-energies is also possible. For a perturbative computation at order O(y 6 t ) in the SM, which is the non-QCD dominant contribution, planar as well as non-planar self-energies with just cubic vertices are required, as you can see in Figure 6. The external lines of the topologies contain only the Higgs field (h). The internal lines (propagators) contain the top quark (t), bottom quark (b), Higgs and Goldstone bosons (G 0 and G ± ) fields. The cubic vertices involved in the computation are hht, G 0 G 0 t and G ± tb. In the Landau gauge, therefore, the Feynman integrals contain four different scales, namely, the Higgs boson mass, the top quark mass, the bottom quark mass and the external momentum of the self-energies. The mass of the bottom quark is considered here as a small variable which is maintained as a non-zero input parameter inside the propagators of the amplitudes, the latter with the aim of avoiding additional infrared divergences, but its contribution to the vertex G ± tb is disregarded when appears in the numerators of the integrals. In particular, Feynman diagrams involving the charge Goldstone boson G ± contain Dirac traces with the chiral matrix γ 5 and an arbitrary number of gamma matrices. Due to the presence of four independent momentum scales, the external momentum and the three loop momenta, the self-energy diagrams can contain contributions of Dirac traces with a single γ 5 and four γ matrices. The details of the regularization of these amplitudes and the consequent reduction and explicit evaluation of the resulting master integrals coming from the topologies in Figure 6 will be presented in a future publication. An extension of this calculation to the MSSM in a suitable regularization scheme that preserves supersymmetry would require also diagrams with quartic vertices involving the stop fields and will allow an evaluation of the non-logarithmic terms at order O(y 6 t ). In fact, even if in the MSSM supersymmetry is explicitly broken, the mechanism is soft therefore the supersymmetric cancellation are more efficiently realized in a regularization scheme, preserving supersymmetry [72]. The inclusion of the external-momentum effects in the self energies is relevant when going beyond the gaugeless limit in the MSSM, or in models with additional Higgs self-couplings. The estimations of those effects in the SQCD sector has to be still determined.
Although the MSSM is the BSM model which has the most precise theoretical predictions of the Higgs boson mass and their associated uncertainties, there are phenomenological indications that suggest the necessity of an extension of its Higgs sector. The searches of BSM Higgs bosons at the RUN2 of the LHC have shown several excesses in the data, from which two facts are remarkable. The first one is the excess showed in several events above the background expectation around the same mass scale of a hypothetical new Higgs boson with a mass of 400GeV [73,74,75,76]. The second one, take into account the possibility of a BSM Higgs boson with a mass lighter than the 125 GeV scale, due to the local excess observed around the scale of about 96 GeV by the CMS [77], confirmed by the LEP experiments [78] as well. In previous analyses focussing just on these two excesses it was shown that in the MSSM neither a 400 GeV Higgs boson nor the CMS excess at around 96 GeV can be realized [79]. Extended SUSY models instead, like the NMSSM [80,81] or the µνSSM [82,83], can account for the excesses at a level of roughly 1σ. An analysis of the possible interpretations of the observed excesses in the context of the N2HDM and the NMSSM models can be found in [84]. A complete explanation of the excesses found by LEP-II and CMS at 96 GeV in the natural NMSSM (nNMSSM) can be also found in [85,86]. Additionally, the minimal SUSY extension of the SM cannot explain the neutrino oscillations and is compatible with the 4.2σ anomaly in the muon g − 2 experiment just for a compressed region of EW parameters. The anomaly could be explained in terms of light EW SUSY particles, as the sleptons or the neutralino, with masses between ∼ 250 GeV and ∼ 700 GeV . If those BSM particles are discovered in this range, then the motivation for improving the theoretical accuracy of the MSSM Higgs-mass prediction will be even stronger. Otherwise, a next to MSSM model suitable to explain the muon g − 2 anomaly must be considered. This issue has been anticipated in recent works. There is an interesting proposal where the MSSM Higgs sector is extended in order to have four Higgs doublets, two of which couple just to quarks and leptons of the third generation, whereas the other two have much smaller vacuum expectation values and provide masses to the first and second generation. The model is called FSSM [87,88] and its region of parameters includes sleptons in the multi-TeV regime, beyond the current reach of the LHC, capable of solving the muon g − 2 anomaly. In particular, the extended Higgs sector of the FSSM can account the experimental (g − 2) µ results in scenarios where the stops need to be considerably heavier than smuons, higgsinos and EW gauginos. Another interesting alternative is the µνSSM which is a highly predictive model that can explain the neutrino oscillations via an EW seesaw mechanism [89], provide a solution of the so-called µ-problem [90] and solve the muon g − 2 anomaly [91] simultaneously. Recently, a vacuum stability analysis that takes into consideration the possibility of long-lived metastable vacua of the neutral scalar potential in the µνSSM was also presented [92]. The matter content of this model is enlarged to include right-handed neutrinos and their corresponding superpartners, the "right-handed" sneutrinos, that are gauge singlet scalar fields. Finally, we point out that, in order to accommodate the observed phenomenology, the MSSM parameters must be enlarged with additional degrees of freedom. This can go against the aim to improve the theoretical uncertainty of the Higgs boson mass since the next to MSSM models contains more particles and additional large couplings. Thus, an estimation of the theoretical uncertainties are also relevant in this context. The application of the existing techniques for the MSSM Higgs mass uncertainty calculations should be useful to accomplish this aim.

Conclusions
In this article, we have presented a general review of the precision calculations of the Higgs boson mass in the SM (M h ) as well as in the rMSSM (m h ). The treelevel expressions of the SM and MSSM Higgs potential together with a discussion about the effects of the dominant quantum corrections have been also reviewed. The theoretical uncertainty associated to M h in the SM can be as large as 1 GeV while for m h the uncertainty estimated in FeynHiggs and confirmed in FlexibleEFTHiggs is of the order of 1 − 5 GeV . In both cases, an improvement of the precision in the computation of the Higgs boson mass is required since the current experimental uncertainty obtained at LHC is ∼ 0.2 GeV and is expected to reach a value of about 50 M eV in future colliders. We have described the main sources of theoretical uncertainties in both models. In particular, we have focused our attention in the three-loop corrections to the Higgs boson mass where the fixed-order calculations based on the Feynman diagrammatic approach have been reported in the limit of zero external-momentum for the Higgs self-energies. We have proposed a way to improve the accuracy of the calculation based on the possibility to estimate the external-momentum effects that is reachable by the recent results presented in [70,71]. On the other hand, the SM and the MSSM predictions are moving away from the current observed phenomenology. We have discussed briefly the alternative models which can accommodate the recent results obtained at the RUN2 of the LHC and the muon g − 2 anomaly found at FERMILAB. In particular, these alternatives include some extensions of the minimal version of the supersymmetric standard model which enlarge the Higgs sector with additional large parameters. This could go against the intention to improve the theoretical uncertainty on the Higgs boson mass and therefore, as a research perspective, we consider relevant to estimate the theory uncertainties for the next to MSSM models discussed in Section 4. We hope that this review will stimulate the scientific discussion with groups interested in this research line.