Next Article in Journal
Techniques for TPC Calibration: Application to Liquid Ar-TPCs
Previous Article in Journal
Acknowledgment to Reviewers of Particles in 2021
Previous Article in Special Issue
Implications of the Conformal Higgs Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

High-Precision Calculations of the Higgs Boson Mass

1
Departamento de Física, Universidad de Pamplona, Pamplona 543057, Colombia
2
Departamento de Física, Universidad Nacional de Colombia, Bogotá 111321, Colombia
*
Author to whom correspondence should be addressed.
Submission received: 31 December 2021 / Revised: 11 February 2022 / Accepted: 15 February 2022 / Published: 17 February 2022
(This article belongs to the Special Issue Higgs Physics)

Abstract

:
In this paper, we 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 even up to the 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 Model is included.

1. Introduction

The Standard Model (SM) of the electroweak force unification does not explain important phenomena of fundamental interactions such as the neutrino oscillation, the dark matter, the baryon asymmetry, the vacuum stability of the model, and the electroweak hierarchy in the mass scales, among others. The explanation of the above phenomena remains without a conclusive solution and requires new physics, 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 increasing discrepancy between the measurement of the muon anomalous magnetic moment and the corresponding SM theoretical prediction. The last report shows a 4.2 σ departure [1] from the reliable SM calculation, which was updated together with the inclusion of non-perturbative contributions in a 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 the 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 is 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 future colliders. A survey of the present status of analytic calculation methods of higher order Feynman diagrams with a 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. BSMs based on supersymmetry (SUSY) where the Higgs boson mass is a prediction of the theory, and their low-energy effective field theories (EFTs), 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 an SM-like Higgs at 125.09 ± 0.21 ( stat . ) ± 0.11 ( syst . ) GeV , but an agreement of the properties of the new particles with the present experimental bounds coming from the searches for dark matter candidates, the muon g 2 anomaly, the LHCb results, and the neutrino oscillations is also mandatory. All of these phenomenological restrictions impose important constraints on the parameter space of the BSMs and can serve as a guide to look for the parameter region where the new particles might have appreciable interactions. In addition, they imply a very accurate computation of the Higgs boson mass since the central experimental value of 125.09 GeV is only qualitative at the tree level in minimal supersymmetric models and the uncertainty associated with the theoretical predictions must reach the experimental accuracy of the LHC and future colliders, estimated to be of the order of 100 MeV [14]. The high-precision calculations of the radiative corrections to the Higgs boson mass are 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 the tree level and then in Section 3 we review 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 us to describe the recent experimental observations, such 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.

2. The Higgs Sector of the SM

In the SM, the Higgs sector contains just one scalar Higgs doublet,
H = ϕ + ϕ 0 ,
which is sufficient to give masses to the three generation of quarks, the three up-type 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
L H = D μ H D μ H V ( H ) .
The first component of the Higgs doublet, ϕ + = ϕ 1 + i ϕ 2 , is a complex field that has electric charge Q = 1 and a 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 equal to T 3 = 1 / 2 . As a consequence, both components have a U ( 1 ) hypercharge Y ( H ) = 2 ( Q T 3 ) = 1 . The covariant derivative, D μ , is defined as
D μ = μ i g 2 τ a W μ a i g 2 B μ ,
where τ a are the Pauli’s matrices and the parameters g and g are the interaction couplings associated with the S U ( 2 ) L and U ( 1 ) Y gauge symmetries, respectively. The electroweak gauge bosons are defined as
W μ ± = W μ 1 i W μ 2 2 ; Z μ 0 A μ = c o s θ W s i n θ W s i n θ W c o s θ W W μ 3 B μ ,
where the EW mixing angle, θ W , is related to the coupling constants g and g and the electric charge, e, by
g = e s i n θ W ; g = e c o s θ W .
Additionally, the classical potential V ( H ) is restricted by the gauge invariance and renormalizability to having the functional form:
V ( H ) = μ 2 H 2 + λ H 4 .
The Higgs quartic coupling λ is a dimensionless positive number so the potential is 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 S U ( 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
V ( H ) | H | 2 = 0 H 2 = μ 2 2 λ 1 2 v 2 .
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 S U ( 2 ) of weak isospin. The election of one of those minima spontaneously breaks the gauge symmetry, and the vacuum is not invariant under the symmetry group of the Lagrangian. 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 S U ( 2 ) symmetry, we can define a new real scalar field h ( x ) in the so-called unitary gauge, such that ϕ 1 = ϕ 2 = 0 , ϕ 0 = ( 1 / 2 ) ( v + h ) and, therefore,
H ( x ) = 1 2 0 v + h ( x ) .
The choice of ϕ 0 as the component of the Higgs doublet that has a vacuum expectation value ensures the conservation of the electric charge. In addition, since only one component has a vev when one chooses a direction, we have three broken continuous symmetries of the original S U ( 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:
V ( H ) = 1 2 2 λ v 2 h 2 + λ v h 3 + 1 4 λ h 4 .
It is possible to associate a physical particle with the quantum field h ( x ) with a tree-level mass proportional to the self-interaction coupling λ and to the vev of H:
M h = 2 λ v .
This neutral zero-spin particle is well known as the SM Higgs boson and, to date, its properties and all the SM predictions are in a great agreement with the Higgs boson properties measured at the LHC by the ATLAS and CMS experiments. It is worth mentioning 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 Equation (2) and after expanding H around its vev, we have:
D μ H D μ H = 1 2 μ h μ h + 1 4 g 2 v 2 W μ + W μ + 1 8 ( g 2 + g 2 ) v 2 Z μ Z μ 1 + h v 2 .
Whereas, when the Lagrangian develops additional mass terms for the EW gauge bosons, the photon remains massless and, therefore,
M W = 1 2 v g ; M Z = 1 2 v g 2 + g 2 ; M γ = 0 .
From the relation of the gauge boson mass M W with the Fermi coupling constant obtained from the muon decay process,
G μ 2 = g 2 8 M W 2 ( 1 + Δ r ) ,
where the Fermi constant G μ 1.17 × 10 5 G e V 2 and Δ r contains the radiative corrections to the relation, and we can estimate the tree-level value for the vev of H,
v = 1 2 G μ 246 GeV .
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, contrary to what happens in the Proca theory, where just the Lorentz covariance is used [18]. On the other hand, quark masses require the existence of left-handed S U ( 2 ) doublets and right-handed S U ( 2 ) singlets in such a way that
L q = y d Q ¯ L H d R + y u Q ¯ L H c u R + h . c .
where
Q ¯ L = u d L ; H c = 1 2 v + h ( x ) 0 ,
for quarks of 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,
m u , d = y u , d 2 v .
Similarly, masses of the leptons are generated from the interaction Lagrangian of the leptons with the Higgs field,
L L = y e L ¯ H e R + H e ¯ R L ,
where y e is the Yukawa coupling for the e-type lepton ( e = e , μ , τ ), L is an S U ( 2 ) doublet with components that contain a left-handed lepton, e L , and a left-handed neutrino, ν e L ,
L = ν e e L ,
and e R are right-handed leptons which are S U ( 2 ) singlets. The mass term for the lepton e is of the form:
m e = y e 2 v .
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 in 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 gaugeless limit approximation for two-loop order are 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 the leading three-loop contributions estimated using the 1PI effective potential approach in the limit where g 2 , g 2 , λ g s 2 , y t 2 are also 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 the 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 s 4 y t 2 ) + O( g s 2 y t 4 ) + O( y t 6 ), 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 from 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 the Q dependence flat. Additionally, 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 very 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], and the metastability is the preferred option with 99.3 % 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 larger extent, contain a theoretical error which reaches 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 coupling constants. This can be done with the SMDR code 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 are different from those used in SMDR. In mr, the renormalized vacuum, v t r e e , 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, as the tadpole diagrams involving the Higgs field include negative powers of the Higgs mass, leading to a parametrically slower convergence. The leading small- M h behavior is dominated by the one-loop correction, and the threshold corrections to M h and λ contain terms that diverge as 1 / M h 2 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 on the gauge-fixing condition [27]. As a consequence, in this framework the MS ¯ masses are gauge-dependent quantities. Nevertheless, an MS ¯ mass is not a physical quantity nor a parameter of the bare Lagrangian, moreover the physical observables are independent of 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 t r e e and the alternative computation implemented in SMDR deserve more attention in future studies.

3. Higgs Sector in Supersymmetry

We start by reviewing the Higgs sector of the minimal supersymmetric extension of the SM (MSSM) at the 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 consist 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 to 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:
H 1 = H 1 0 + 1 2 v 1 H 1 and H 2 = e i ψ 2 H 2 + H 2 0 + 1 2 v 2 ,
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 j 0 have neutral components which are vevless scalar fields, H j 0 = ϕ j 0 + i χ j 0 , 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 ) e m symmetry will not be spontaneously broken. The Higgs sector comes from the bare Lagrangian
L V H = ( D σ H 1 ) D σ H 1 + ( D σ H 2 ) D σ H 2 V ( H 1 , H 2 ) .
The kinetic term contains the covariant derivative
D σ = σ + i g τ a 2 W σ a + i g Y 2 B μ + i g s λ α 2 G σ α ,
where g , g, and g s are the U ( 1 ) Y , SU ( 2 ) L , and SU ( 3 ) gauge couplings, respectively. The SU ( 2 ) generator Pauli matrices are represented by τ a while λ α stands for the SU ( 3 ) Gell-Mann matrices. The classical Higgs potential, V ( H 1 , H 2 ) , is given by
V ( H 1 , H 2 ) = μ 2 + m H 1 2 H 1 2 + μ 2 + m H 2 2 H 2 2 + b ϵ a b H 1 a H 2 b + h . c . + 1 2 g 2 H 1 H 2 2 + 1 8 g 2 + g 2 H 2 2 H 1 2 2 .
The quadratic terms proportional to the higgsino mass parameter μ 2 in Equation (24) 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 also obtains 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 applications, one must introduce new terms that explicitly break supersymmetry and represent the low-energy effects of the unknown breaking mechanism. In Equation (24), the parameters m H 1 2 , m H 2 2 , 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 H 1 2 and m H 2 2 , which can of course have both signs, are needed in order to break the EW symmetry. The b-term is the only one that depends on the phases of the fields. The parameter b is real and positive and the fields H 1 0 and H 2 0 have equal and opposite phases which can both be reduced 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 super-renormalizable, with a 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, and the cancellation of quadratically divergent corrections to scalar masses is maintained for all perturbative orders. In addition, 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 Equation (24), we can derive the linear part of the Higgs potential on the basis ( ϕ j 0 , χ j 0 , H j ± ) where j = 1 , 2 . For the fields χ j 0 and H j ± , there are no contributions since the rMSSM Higgs potential is invariant under CP transformation. Thus, the linear terms of the effective potential T 1 ϕ 1 0 + T 2 ϕ 2 0 , with the coefficients T j , better known as Higgs tadpoles, have the expressions
T 1 2 v 1 = μ 2 + m H 1 2 b v 2 v 1 + 1 8 g 2 + g 2 v 1 2 v 2 2 , T 2 2 v 2 = μ 2 + m H 2 2 b v 1 v 2 + 1 8 g 2 + g 2 v 2 2 v 1 2 .
As the vevs v 1 and v 2 minimize the Higgs potential, the following stationary conditions are satisfied:
V | H j 0 | H j 0 = 0 ; H j ± = 0 = T j = 0 ; j = 1 , 2 .
From Equation (26), one can identify the conditions required for the stable minimum of the effective potential. First, note that along the direction | H 1 0 | = | H 2 0 | the potential will be unbounded from below and therefore it does not have a minimum unless
2 μ 2 + m H 1 2 + m H 2 2 > 2 b .
For μ 2 + m H 1 2 and μ 2 + m H 2 2 positive, the origin is not a minimum but a saddle point, and the minimum occurs at non-zero vevs of H j 0 , when
μ 2 + m H 1 2 μ 2 + m H 2 2 < b 2 .
The rMSSM Higgs potential develops a stable minimum if the conditions of the Equations (27) and (28) are met.
Turning to the bilinear part of the kinetic terms in the Lagrangian L V H (Equation (22)), the masses of the EW gauge bosons amount to
M W 2 = 1 4 g 2 v 1 2 + v 2 2 ; M Z 2 = 1 4 ( g 2 + g 2 ) v 1 2 + v 2 2 ; M γ 2 = 0 .
The Z boson mass determines the tree-level relation v 1 2 + v 2 2 174 GeV. The mass matrices of the rMSSM Higgs bosons can be identified from the bilinear part of the classical potential,
ϕ 1 0 ϕ 2 0 M ϕ 0 ϕ 1 0 ϕ 2 0 + χ 1 0 χ 2 0 M χ 0 χ 1 0 χ 2 0 + H 1 + H 2 + M H ± H 1 H 2 .
To derive the M ϕ 0 -matrix, the relations from Equations (25) and (26) and the definitions
M A 2 = b c o t β + t a n β , t a n β = v 2 v 1 ; 0 β π 2
are required. The tree-level mass matrix of the neutral ϕ 1 , 2 0 -bosons reads
M ϕ 0 2 = M Z 2 c β 2 + M A 2 s β 2 + T 1 2 v 1 M A 2 + M Z 2 s β c β M A 2 + M Z 2 s β c β M Z 2 s β 2 + M A 2 c β 2 + T 2 2 v 2 .
We have used the short notation s β = s i n ( β ) and c β = c o s ( β ) and we have explicitly written the contributions of the Higgs tadpoles, which vanish at the tree level according to Equation (26), because they will receive non-zero loop contributions when renormalization of the Higgs masses is 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 Equation (31). On the other hand, the tree-level mass matrices of the χ 1 , 2 0 and H 1 , 2 ± components are given by
M χ 0 2 = M A 2 s β 2 + T 1 2 v 1 M A 2 s β c β M A 2 s β c β M A 2 c β 2 + T 2 2 v 2
and
M H ± 2 = M A 2 s β 2 + T 1 2 v 1 + 1 2 g 2 v 2 2 M A 2 s β c β + 1 2 g 2 v 1 v 2 M A 2 s β c β + 1 2 g 2 v 1 v 2 M A 2 c β 2 + T 2 2 v 2 + 1 2 g 2 v 1 2 .
The potential can be given a physical basis, where the quadratic terms in the components of H j are diagonalized, through the rotations
ϕ 1 0 ϕ 2 0 = D α H h , χ 1 0 χ 2 0 = D β G 0 A , H 1 ± H 2 ± = D β G ± H ± ,
via the orthogonal transformation
D θ = c θ s θ s θ c θ .
On this basis, the Higgs sector has five physical Higgs bosons, three of them are neutral: The lightest (h) and heavy (H) CP-even Higgs bosons and the CP-odd Higgs boson (A). The other two, H ± , are charged and vevless. There are also three unphysical massless Goldstone bosons, G 0 and G ± , which are absorbed by the EW gauge fields to build up their longitudinal components just as in the SM. The angle β is linked to the vevs through Equation (31) while α can be determined from the rotation of Equation (32) in terms of the MSSM parameters,
t a n ( 2 α ) = t a n ( 2 β ) M A 2 + M Z 2 M A 2 M Z 2 ; π 2 < α < 0 .
Thus, after diagonalization, besides the EW boson masses of Equation (29), the rMSSM Higgs sector is parametrized in terms of two additional parameters: t a n β and the mass of the CP-odd Higgs boson m A . The masses of the charged Higgs bosons, m ± , are linearly dependent on m A , and they are usually used in the complex version of the MSSM. At the tree level, we have:
m A 2 = M A 2 + T 1 2 v 1 s β 2 + T 2 2 v 2 c β 2 ; m ± 2 = m A 2 + M W 2 .
The tree-level masses of the CP-even Higgs boson particles, h and H, follow as predictions
m h , H 2 = 1 2 m A 2 + M Z 2 m A 2 + M Z 2 2 4 m A 2 M Z 2 c o s 2 2 β .
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 is predicted in the MSSM. From the mass Formulas (38) and (39), the next important inequalities can be derived:
m h M Z ; m A m H ; M W m ± .
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, up to 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 the art of the quantum corrections to the lightest Higgs boson mass in the MSSM is very advanced and widely studied. The reported results include fixed-order diagrammatic calculations, resummation of the large logarithms in effective field theories (EFTs), and hybrid fixed-order and EFT calculations up to the 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 ˜ t 1 and m ˜ t 2 :
δ t ( 1 ) m h 2 = 3 G F 2 π 2 s β 2 M t 4 ln m ˜ t 1 m ˜ t 2 M t 2 .
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:
m h 2 = 1 2 m A 2 + δ t ( 1 ) m h 2 m A 2 + δ t ( 1 ) m h 2 2 4 m A 2 s β 2 δ t ( 1 ) m h 2 , m H 2 = m A 2 m h 2 + δ t ( 1 ) m h 2 ; m ± 2 = m A 2 + M W 2 .
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:
δ X t ( 1 ) m h 2 = 3 G F M t 4 X t 2 2 2 π 2 s β 2 2 f m ˜ t 1 2 , m ˜ t 2 2 + 2 m ˜ t 1 2 + m ˜ t 2 2 f m ˜ t 1 2 , m ˜ t 2 2 m ˜ t 1 2 m ˜ t 2 2 2 X t 2 ,
where the function f is defined as:
f x , y = 1 x y ln x y .
The dependence of the CP-even and charged Higgs boson masses on the parameter m A in Equation (42), including the dominant one-loop radiative corrections in a scenario where the squark masses are at the same supersymmetric scale M S U S Y , produces the curves in Figure 2. Note that, independently of t a n β , 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 an asymptotic value, showing a more regular dependence on m A for large t a n β , where m h m A for small m A while m h c o n s t 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 the one-loop level, yielding the expression:
m h 2 M Z 0 2 c 2 β 2 + δ t ( 1 ) m h 2 s β 2 + δ X t ( 1 ) m h 2 s β 2 .
If one considers large values of the parameter t a n β ( t a n β 30 ) in the scenario of maximal stop mixing, where the value of X t makes m h maximal (frequently referred to in literature as the m h m a x scenario), a general upper bound given by m h 140 GeV is found. Fortunately for the MSSM, the Higgs boson was discovered at the LHC within this energy region. The contributions of Equations (41) and (43) contain the bulk of the one-loop corrections. The subdominant contributions coming from higher-loop corrections can essentially be reduced to higher-order SQCD effects.
Figure 3—Left, shows the dependence of the tree-level, full one-loop, and leading QCD two-loop m h -predictions coming from FeynHiggs [43] as a function of the scale m A in the maximal stop mixing scenario for the same election of parameters as in Figure 2. In this plot, one can clearly see the good behavior of the perturbative expansion. If one looks at the numerical difference between the dashed and the dot-dashed curves as well as the difference between the dashed and the solid lines for large m A , one can easily conclude that the one-loop corrections (estimated to be about 45 GeV) are larger than the two-loop ones (roughly 12 GeV), as expected from a perturbative quantum field theory. Additionally, the two-loop corrections give negative contributions in contrast to the one-loop corrections which are positive. In the right panel, the three-loop predictions of m h at order α t α s 2 are depicted, firstly predicted in [44,45], as a function of m A (dot-dashed red line) in the m h m a x scenario. In general, the complete m h -prediction is built up as the sum of different contributions. The tree-level value accounts for about 60% of the renormalized mass, the one-loop correction represents about 32% of m h , the two-loop correction contributes about 6%, and finally the three-loop contribution represents about 1% of the total mass. The sizes of the three-loop quantum corrections are of the order of 1 GeV and have an opposite signs regarding the one-loop corrections for small values of M S U S Y of about 1 TeV. However, due to the lack of experimental evidence for SUSY particles at this energy, we have to consider benchmark scenarios with even higher SUSY energy scales. Considering the limit where all the soft SUSY-breaking masses as well as the CP-odd Higgs mass ( m A ) lie around the characteristic scale M S U S Y , the dependence of the m h -predictions on this scale can be studied.
In Figure 4—Left, we present a numerical comparison of our three-loop fixed-order predictions of m h to the two-loop predictions coming from FeynHiggs for large SUSY scales in the range 0.5 TeV M S U S Y 40 TeV . The dependence of m h on M S U S Y is evaluated for a vanishing stop mixing, X t = 0 , and for t a n β = 10 . The two- and three-loop fixed-order results are represented with the black dashed and the red dot-dashed lines, respectively. The brown band is the experimental Higgs boson mass measured at the LHC and its corresponding uncertainty. Note that for scales above 10 TeV the difference between the two-loop and three-loop results rapidly increases up to about 25 GeV when M S U S Y grows to up to 20 TeV and grows monotonically, reaching about 100 GeV at M S U S Y = 40 TeV. This pronounced behavior depends crucially on our election of the input parameters and is a consequence of the presence of n-loop logarithms of the form l o g n ( M S U S Y / M t ) in the three-loop fixed-order Higgs self-energy corrections. Those n-loop logarithms are known as large logarithms because they are the source of additional large contributions in the three-loop predictions of m h when there is a large hierarchy between the EW scale M t and the SUSY scale M S U S Y . The large logarithm contributions can spoil the convergence of the perturbative expansion and yield unreliable predictions of the Higgs boson mass. A fixed-order computation is thus recommended for low values of M S U S Y not separated too much from M t . 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 S U S Y can be evaluated up to superhigh energy scales of the order of the SUSY grand unified theory (GUT) 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 S U S Y between 10 7 GeV M S U S Y 10 16 GeV and for t a n β = 1 . As was also argued in [46], it is possible to fix upper bounds on the SUSY scale M S U S Y for scenarios with large values of t a n β ( tan β > 10 ) if we impose that the m h -predictions in the MSSM are compatible with the combined CMS/ATLAS measured Higgs boson mass within the actual combined uncertainties. However, for small values of t a n β ( t a n β < 10 ), such upper bounds coming from the LHC measurements cannot be derived and a strong dependence between M S U S Y and X t is observed. In fact, for large values of the stop mixing parameter, | X t | 3 M S U S Y , the m h -prediction is compatible with the LHC measurement for M S U S Y 10 7 GeV . It is possible, however, to derive upper bounds over M S U S Y for small t a n β values if we impose the positivity of the running Higgs self-coupling λ ( Q ) from the EW scale until the next threshold, Q = M S U S Y , which is equivalent to imposing the stability of the SM renormalization group improved Higgs effective potential,
V ( h ) R G I λ ( Q ) 4 h 4 ( Q ) ,
up to the scale where supersymmetry is supposed to appear. The SM is considered here as a low-energy EFT of the rMSSM in the limit where we integrate out all SUSY QCD particles at the scale M S U S Y , 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 therein. In the gray region of Figure 4—Right. we have pictured the points where the quartic Higgs self-coupling is negative at the SUSY scale, λ ( M S U S Y ) < 0 , leading to an unstable electroweak vacuum. If the stability of the EW vacuum is required, then the parameter space in the gray region is excluded and the intersection of the gray border line with the brown region could constrain the M S U S Y value as a function of m h . Note that the region of exclusion occurs, in general, around | X t / M S U S Y | = 0 and for large values of X t , | X t / M S U S Y | > 3 . Specifically, for a 125.09 GeV Higgs boson mass (dashed brown line) and a stop mixing of | X t / M S U S Y | = 1.2 , the intersection of the gray curve with the dashed brown line provides the upper bound M S U S Y 10 10 GeV . Assuming an uncertainty of ± 5 GeV on the Higgs boson mass, this bound can reach a maximum value of M S U S Y 10 11 GeV for | X t / M S U S Y | 1.4 as was also verified in [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 future colliders. An experimental test requires, therefore, the design of experiments able 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 caused a special interest in 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 these superhigh energies are not the only target where we could expect some SUSY evidence. There are recent phenomenological analyses 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].

4. Theoretical Uncertainties and Next to Minimal Supersymmetry

The computation of the Higgs boson mass in SUSY models requires imposing 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 a 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 TeV [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 allow a numerical estimation of the different theoretical errors associated with 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 improved perturbative computations in the future.
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 two-loop 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 also require 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 α s 2 ) while the calculation of the three-loop corrections that involve lower powers of the strong gauge coupling, e.g., O( α t 2 α s ), O( α t 3 ), 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 to date been performed only for the simplest heavy SUSY scenario where the MSSM is matched directly to the SM and all the SUSY particles are at the same scale M S U S Y in the limit M S U S Y M t . The next to next leading logarithms (NNLLs) 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 S U S Y which is varied from 0.4 TeV to 100 TeV. The stop mixing parameter is taken to be X t = 2.45 M S U S Y while t a n β = 20 . In addition, the squark masses and X t are running parameters evaluated at the scale Q = M S U S Y in the D R ¯ 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 S U S Y 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 the LHC (∼200 MeV) 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 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 yet been 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 with the effective potential approach which is equivalent to computing 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 allow 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 numerically compute the renormalized ϵ -finite master integrals for arbitrary external momentum invariants and, in principle, for arbitrary masses, the evaluation of non-planar self-energies is also possible. For a perturbative computation at order O( y t 6 ) 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 seen 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 h h t , G 0 G 0 t , and G ± t b . 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 ± t b is disregarded when it 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 also require diagrams with quartic vertices involving the stop fields and will allow an evaluation of the non-logarithmic terms at order O( y t 6 ). In fact, even if the MSSM supersymmetry is explicitly broken, the mechanism is soft, therefore the supersymmetric cancellation is 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 are still to be determined.
Although the MSSM is the BSM which has the most precise theoretical predictions of the Higgs boson mass and the 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 400 GeV [73,74,75,76]. The second one, taking 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], was confirmed by the LEP experiments [78] as well. In previous analyses focusing 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 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 also be 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, such 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 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 for 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 MSSMs contain more particles and additional large couplings. Thus, an estimation of the theoretical uncertainties is 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.

5. 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 tree-level expressions of the SM and MSSM Higgs potential together with a discussion about the effects of the dominant quantum corrections have been reviewed. The theoretical uncertainty associated with 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 the LHC is ∼0.2 GeV and is expected to reach a value of about 50 MeV in future colliders. We have described the main sources of theoretical uncertainties in both models. In particular, we have focused our attention on 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 it relevant to estimate the theory uncertainties for the next to MSSMs discussed in Section 4. We hope that this review will stimulate scientific discussion with groups interested in this research line.

Author Contributions

Conceptualization, E.R., R.F.; data curation, E.R. and R.F.; formal analysis, E.R., R.F.; investigation, E.R. and R.F.; methodology, E.R. and R.F.; project administration, E.R. and R.F.; resources, E.R. and R.F.; software, E.R. and R.F.; supervision, E.R. and R.F.; validation, E.R. and R.F.; visualization, E.R. and R.F.; writing—original draft, E.R. and R.F.; writing—review and editing, E.R. and R.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abi, B.; Albahri, T.; Al-Kilani, S.; Allspach, D.; Alonzi, L.P.; Anastasi, A.; Anisenkov, A.; Azfar, F.; Badgley, K.; Baeßler, S.; et al. Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm. Phys. Rev. Lett. 2021, 126, 141801. [Google Scholar] [PubMed]
  2. Aoyama, T.; Asmussen, N.; Benayoun, M.; Bijnens, J.; Blum, T.; Bruno, M.; Caprini, I.; Carloni Calame, C.M.; Cè, M.; Colangelo, G.; et al. The anomalous magnetic moment of the muon in the Standard Model. Phys. Rep. 2020, 887, 1–166. [Google Scholar] [CrossRef]
  3. Aaij, R.; Abellán Beteta, C.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H.; Aiola, S.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; et al. (LHCb). Test of lepton universality in beauty-quark decays. arXiv 2021, arXiv:2103.11769. [Google Scholar]
  4. Apollinari, G.; Béjar Alonso, I.; Brüning, O.; Fessia, P.; Lamont, M.; Rossi, L.; Tavian, L. High-Luminosity Large Hadron Collider (HL-LHC) Technical Design Report V. 0.1; CERN: Geneva, Switzerland, 2017. [Google Scholar] [CrossRef]
  5. Baer, H.; Barklow, T.; Fujii, K.; Gao, Y.; Hoang, A.; Kanemura, S.; List, J.; Logan, H.E.; Nomerotski, A.; Perelstein, M.; et al. The International Linear Collider Technical Design Report—Volume 2: Physics. Available online: https://arxiv.org/abs/1306.6352 (accessed on 30 December 2021).
  6. Barklow, T.; Brau, J.; Fujii, K.; Gao, J.; List, J.; Walker, N.; Yokoya, K. ILC Operating Scenarios. Available online: https://arxiv.org/abs/1506.07830 (accessed on 30 December 2021).
  7. Bicer, M.; Duran Yildiz, H.; Yildiz, I.; Coignet, G.; Delmastro, M.; Alexopoulos, T.; Grojean, C.; Antusch, S.; Sen, T.; He, H.-J.; et al. First Look at the Physics Case of TLEP. J. High Energy Phys. 2014, 1401, 164. [Google Scholar] [CrossRef] [Green Version]
  8. Blondel, A.; Freitas, A.; Gluza, J.; Heinemeyer, S.; Jadach, S.; Janot, P.; Riemann, T. Theory Requirements and Possibilities for the FCC-ee and other Future High Energy and Precision Frontier Lepton Colliders. Available online: https://arxiv.org/abs/1901.02648 (accessed on 30 December 2021).
  9. Blondel, A.; Gluza, J.; Jadach, S.; Janot, P.; Riemann, T.; Abreu, S.; Aguilera-Verdugo, J.J.; Arbuzov, A.B.; Baglio, J.; Bakshi, S.D.; et al. Theory for the FCC-ee: Report on the 11th FCC-ee Workshop; CERN: Geneva, Switzerland, 2020. [Google Scholar] [CrossRef]
  10. Ahmad, M.; Alves, D.; An, H.; An, Q.; Arhrib, A.; Arkani-Hamed, N.; Ahmed, I.; Bai, Y.; Ferroli, R.B.; Ban, Y.; et al. CEPC-SPPC Preliminary Conceptual Design Report. 2. Accelerator. Available online: https://inspirehep.net/literature/1395736 (accessed on 30 December 2021).
  11. Blümlein, J. Analytic Integration Methods in Quantum Field Theory: An Introduction. Available online: https://arxiv.org/abs/2103.10652 (accessed on 30 December 2021).
  12. ATLAS; CMS Collaboration. Combined Measurement of the Higgs Boson Mass in pp Collisions at s = 7 and 8 TeV with the ATLAS and CMS Experiments. Phys. Rev. Lett. 2015, 114, 191803. [Google Scholar] [CrossRef] [Green Version]
  13. ATLAS; CMS Collaboration. Measurements of the Higgs boson production and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the LHC pp collision data at s = 7 and 8 TeV. J. High Energy Phys. 2016, 8, 45. [Google Scholar] [CrossRef]
  14. Aguilar-Saavedra, J.A.; Alcaraz, J.; Ali, A.; Ambrosanio, S.; Andreazza, A.; Andruszkow, J.; Badelek, B.; Ballestrero, A.; Barklow, T.; Bartl, A.; et al. TESLA: The Superconducting Electron Positron Linear Collider with an Integrated X-ray Laser Laboratory. Technical Design Report. Part 3. Physics at an e+e Linear Collider. Available online: https://inspirehep.net/literature/559272 (accessed on 30 December 2021).
  15. Higgs, P.W. Broken Symmetries and the Masses of Gauge Bosons. Phys. Rev. Lett. 1964, 13, 508. [Google Scholar]
  16. Englert, F.; Brout, R. Broken Symmetry and the Mass of Gauge Vector Mesons. Phys. Rev. Lett. 1964, 13, 321. [Google Scholar]
  17. Iliopoulos, J.; Itzykson, C.; Martin, A. Functional methods and perturbation theory. Rev. Mod. Phys. 1975, 47, 165. [Google Scholar]
  18. Becchi, C. Introduction to Gauge Theories. Available online: https://arxiv.org/abs/hep-ph/9705211 (accessed on 30 December 2021).
  19. Bezrukov, F.; Kalmykov, M.Y.; Kniehl, B.A.; Shaposhnikov, M. Higgs Boson Mass and New Physics. Available online: https://arxiv.org/pdf/1205.2893.pdf (accessed on 30 December 2021).
  20. Degrassi, G.; Di Vita, S.; Elias-Miro, J.; Espinosa, J.R.; Giudice, G.F.; Isidori, G.; Strumia, A. Higgs mass and vacuum stability in the Standard Model at NNLO. J. High Energy Phys. 2012, 1208, 098. [Google Scholar] [CrossRef] [Green Version]
  21. Buttazzo, D.; Degrassi, G.; Giardino, P.P.; Giudice, G.F.; Sala, F.; Salvio, A.; Strumia, A. Investigating the near-criticality of the Higgs boson. J. High Energy Phys. 2013, 2013, 89. [Google Scholar] [CrossRef] [Green Version]
  22. Martin, S.P.; Robertson, D.G. Higgs boson mass in the Standard Model at two-loop order and beyond. Phys. Rev. D 2014, 90, 073010. [Google Scholar] [CrossRef] [Green Version]
  23. Martin, S.P.; Robertson, D.G. Standard Model parameters in the tadpole-free pure MS scheme. Phys. Rev. D 2019, 100, 073004. [Google Scholar] [CrossRef] [Green Version]
  24. Kniehl, B.A.; Pikelner, A.F.; Veretin, O.L. mr: A C++ library for the matching and running of the Standard Model parameters. Comput. Phys. Commun. 2016, 206, 84–96. [Google Scholar] [CrossRef] [Green Version]
  25. Kniehl, B.A.; Pikelner, A.F.; Veretin, O.L. Two-loop electroweak threshold corrections in the Standard Model. Nucl. Phys. B 2015, 896, 19–51. [Google Scholar] [CrossRef]
  26. Sirlin, A.; Zucchini, R. Dependence of the Higgs coupling hMS(M)onmH and the possible onset of new physics. Nucl. Phys. B 1986, 266, 389–409. [Google Scholar] [CrossRef]
  27. Fleischer, J.; Jegerlehner, F. Radiative corrections to Higgs-boson decays in the Weinberg-Salam model. Phys. Rev. D 1981, 23, 2001. [Google Scholar] [CrossRef]
  28. Coleman, S.R.; Mandula, J. All Possible Symmetries of the S Matrix. Phys. Rev. 1967, 159, 1251. [Google Scholar] [CrossRef]
  29. Haag, R.; Lopuszanski, J.T.; Sohnius, M. All Possible Generators of Supersymmetries of the S Matrix. Nucl. Phys. B 1975, 88, 257. [Google Scholar] [CrossRef]
  30. Langacker, P. The Standard Model and Beyond; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  31. Fayet, P.; Ferrara, S. Supersymmetry. Phys. Rep. 1977, 32, 249–334. [Google Scholar] [CrossRef]
  32. Nilles, H.P. Supersymmetry, supergravity and particle physics. Phys. Rep. 1984, 110, 1–162. [Google Scholar] [CrossRef]
  33. Haber, H.E.; Kane, G.L. The search for supersymmetry: Probing physics beyond the standard model. Phys. Rep. 1985, 117, 75–263. [Google Scholar] [CrossRef] [Green Version]
  34. Drees, M.; Godbole, R.; Roy, P. Theory and Phenomenology of Sparticles; World Scientific: Singapore, 2005. [Google Scholar]
  35. Kuroda, M. Complete Lagrangian of MSSM. Available online: https://arxiv.org/abs/hep-ph/9902340 (accessed on 30 December 2021).
  36. Terning, J. Modern Supersymmetry—Dinamics and Duality; CLARENDON PRESS: Oxford, UK, 2006. [Google Scholar]
  37. Baer, H.; Tata, X. Weak Scale Supersymmetry: From Superfields to Scattering Events; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
  38. Aitchison, I.J.R. Supersymmetry in Particle Physics: An Elementary Introduction; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  39. Witten, E. Dynamical breaking of supersymmetry. Nucl. Phys. B 1981, 188, 513–554. [Google Scholar] [CrossRef]
  40. Dimopoulos, S.; Georgi, H. Softly broken supersymmetry and SU(5). Nucl. Phys. B 1981, 193, 150–162. [Google Scholar] [CrossRef] [Green Version]
  41. Girardello, L.; Grisaru, M.T. Soft breaking of supersymmetry. Nucl. Phys. B 1982, 194, 65–76. [Google Scholar]
  42. Slavich, P.; Heinemeyer, S.; Bagnaschi, E.; Bahl, H.; Goodsell, M.; Haber, H.E.; Hahn, T.; Harlander, R.; Hollik, W.; Lee, G.; et al. Higgs-mass predictions in the MSSM and beyond. Eur. Phys. J. C 2021, 81, 450. [Google Scholar] [CrossRef]
  43. Bahl, H.; Hahn, T.; Heinemeyer, S.; Hollik, W.; Paßehr, S.; Rzehak, H.; Weiglein, G. Precision calculations in the MSSM Higgs-boson sector with FeynHiggs 2.14. Comput. Phys. Commun. 2020, 249, 107099. [Google Scholar] [CrossRef] [Green Version]
  44. Fazio, A.R.; Reyes R, E.A. The Lightest Higgs Boson Mass of the MSSM at Three-Loop Accuracy. Nucl. Phys. B 2019, 942, 164–183. [Google Scholar] [CrossRef]
  45. Reyes R, E.A. NNLO Vacuum Stability Analysis in the SM and the Higgs Bosons of the MSSM at Three-Loop Accuracy. Available online: https://repositorio.unal.edu.co/handle/unal/77990 (accessed on 15 June 2020).
  46. Reyes R, E.A.; Fazio, A.R. Comparison of the EFT Hybrid and Three-Loop Fixed-Order Calculations of the Lightest MSSM Higgs Boson Mass. Phys. Rev. D 2019, 100, 115017. [Google Scholar] [CrossRef] [Green Version]
  47. Bagnaschi, E.; Degrassi, G.; Paßehr, S.; Slavich, P. Full two-loop QCD corrections to the Higgs mass in the MSSM with heavy superpartners. Eur. Phys. J. C 2019, 79, 910. [Google Scholar] [CrossRef]
  48. Allanach, B.C.; Voigt, A. Uncertainties in the Lightest CP Even Higgs Boson Mass Prediction in the Minimal Supersymmetric Standard Model: Fixed Order Versus Effective Field Theory Prediction. Eur. Phys. J. C 2018, 78, 573. [Google Scholar] [CrossRef]
  49. Gelmini, G.B. Search for dark matter. Int. J. Mod. Phys. A 2008, 23, 4273–4288. [Google Scholar] [CrossRef]
  50. Feng, J.L. Dark Matter Candidates from Particle Physics and Methods of Detection. Ann. Rev. Astron. Astrophys. 2010, 48, 495–545. [Google Scholar] [CrossRef] [Green Version]
  51. Frenk, C.S.; White, S.D.M. Dark matter and cosmic structure. Ann. Phys. 2012, 524, 507–534. [Google Scholar] [CrossRef] [Green Version]
  52. Khlopov, M. Cosmological Probes for Supersymmetry. Symmetry 2015, 7, 815–842. [Google Scholar] [CrossRef] [Green Version]
  53. Ketov, S.V.; Khlopov, M.Y. Cosmological Probes of Supersymmetric Field Theory Models at Superhigh Energy Scales. Symmetry 2019, 11, 511. [Google Scholar] [CrossRef] [Green Version]
  54. Chakraborti, M.; Heinemeyer, S.; Saha, I. Improved (g-2)μ Measurements and Supersymmetry: Implications for e+e Colliders. Available online: https://arxiv.org/abs/2105.06408 (accessed on 30 December 2021).
  55. ATLAS Collaboration. Search for direct production of charginos, neutralinos and sleptons in final states with two leptons and missing transverse momentum in pp collisions at s = 8 TeV with the ATLAS detector. J. High Energy Phys. 2014, 5, 071. [Google Scholar] [CrossRef] [Green Version]
  56. ATLAS Collaboration. Searches for electroweak production of supersymmetric particles with compressed mass spectra in s = 13 TeV pp collisions with the ATLAS detector. Phys. Rev. D 2020, 101, 052005. [Google Scholar] [CrossRef] [Green Version]
  57. ATLAS Collaboration. Search for electroweak production of charginos and sleptons decaying into final states with two leptons and missing transverse momentum in s = 13 TeV pp collisions using the ATLAS detector. Eur. Phys. J. C 2020, 80, 123. [Google Scholar] [CrossRef] [Green Version]
  58. ATLAS Collaboration. Search for direct stau production in events with two hadronic -leptons in s = 13 TeV pp collisions with the ATLAS detector. Phys. Rev. D 2020, 101, 032009. [Google Scholar] [CrossRef] [Green Version]
  59. Arbey, A.; Battaglia, M.; Djouadi, A.; Mahmoudi, F.; Muhlleitner, M.; Spira, M. Higgs Properties and Supersymmetry: Constraints and Sensitivity from the LHC to an e+e Collider. Available online: https://arxiv.org/abs/2201.00070 (accessed on 30 December 2021).
  60. Redi, M.; Tesi, A. Implications of a Light Higgs in Composite Models. J. High Energy Phys. 2012, 10, 166. [Google Scholar] [CrossRef] [Green Version]
  61. Bellazzini, B.; Csáki, C.; Serra, J. Composite Higgses. Eur. Phys. J. C 2014, 74, 2766. [Google Scholar] [CrossRef] [Green Version]
  62. Dobado, A.; Espriu, D. Strongly coupled theories beyond the Standard Model. Prog. Part. Nucl. Phys. 2020, 115, 103813. [Google Scholar] [CrossRef]
  63. Branco, G.C.; Ferreira, P.M.; Lavoura, L.; Rebelo, M.N.; Sher, M.; Silva, J.P. Theory and phenomenology of two-Higgs-doublet models. Phys. Rep. 2012, 516, 1–102. [Google Scholar] [CrossRef] [Green Version]
  64. Bernon, J.; Gunion, J.F.; Haber, H.E.; Jiang, Y.; Kraml, S. Scrutinizing the Alignment Limit in Two-Higgs-Doublet Models. Part 1: mh = 125 GeV. Phys. Rev. D 2015, 92, 075004. [Google Scholar] [CrossRef] [Green Version]
  65. Goodsell, M.D.; Paßehr, S. All two-loop scalar self-energies and tadpoles in general renormalisable field theories. Eur. Phys. J. C 2020, 80, 417. [Google Scholar] [CrossRef]
  66. Bahl, H.; Heinemeyer, S.; Hollik, W.; Weiglein, G. Theoretical uncertainties in the MSSM Higgs boson mass calculation. Eur. Phys. J. C 2020, 80, 497. [Google Scholar] [CrossRef]
  67. Athron, P.; Park, J.-H.; Steudtner, T.; Stöckinger, D.; Voigt, A. Precise Higgs mass calculations in (non-)minimal supersymmetry at both high and low scales. J. High Energy Phys. 2017, 1, 079. [Google Scholar] [CrossRef] [Green Version]
  68. Kwasnitza, T.; Stöckinger, D.; Voigt, A. Improved MSSM Higgs mass calculation using the 3-loop FlexibleEFTHiggs approach including xt-resummation. J. High Energy Phys. 2020, 7, 197. [Google Scholar] [CrossRef]
  69. Stöckinger, D.; Unger, J. Three-loop MSSM Higgs-Boson Mass Predictions and Regularization by Dimensional Reduction. Nucl. Phys. B 2018, 935, 1–16. [Google Scholar] [CrossRef]
  70. Bauberger, S.; Freitas, A.; Wiegand, D. TVID 2: Evaluation of planar-type three-loop self-energy integrals with arbitrary masses. J. High Energy Phys. 2020, 2001, 024. [Google Scholar] [CrossRef] [Green Version]
  71. Martin, S.P. Renormalized ε-Finite Master Integrals and Their Virtues: The Three-Loop Self Energy Case. Available online: https://arxiv.org/abs/2112.07694 (accessed on 30 December 2021).
  72. Gnendiger, C.; Signer, A.; Stockinger, D.; Broggio, A.; Cherchiglia, A.L.; Driencourt-Mangin, F.; Fazio, A.R.; Hiller, B.; Mastrolia, P.; Peraro, T.; et al. To d, or not to d: Recent developments and comparisons of regularization schemes. Eur. Phys. J. C 2017, 77, 471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. ATLAS Collaboration. Search for heavy resonances decaying into a W or Z boson and a Higgs boson in final states with leptons and b-jets in 36 fb−1 of s = 13 TeV pp collisions with the ATLAS detector. J. High Energy Phys. 2018, 3, 174. [Google Scholar] [CrossRef] [Green Version]
  74. CMS Collaboration. Search for a heavy pseudoscalar boson decaying to a Z and a Higgs boson at s = 13 TeV. Eur. Phys. J. C 2019, 79, 564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. CMS Collaboration. Search for heavy Higgs bosons decaying to a top quark pair in proton-proton collisions at s = 13 TeV. J. High Energy Phys. 2020, 4, 171. [Google Scholar] [CrossRef]
  76. ATLAS Collaboration. Search for heavy Higgs bosons decaying into two tau leptons with the ATLAS detector using pp collisions at s = 13 TeV. Phys. Rev. Lett. 2020, 125, 051801. [Google Scholar] [CrossRef]
  77. CMS Collaboration. Search for a standard model-like Higgs boson in the mass range between 70 and 110 GeV in the diphoton final state in proton-proton collisions at s = 8 and 13 TeV. Phys. Lett. B 2019, 793, 320. [Google Scholar] [CrossRef]
  78. ALEPH Collaboration; DELPHI Collaboration; L3 Collaboration; OPAL Collaboration; The LEP Working Group for Higgs Boson Searches. Search for the Standard Model Higgs boson at LEP. Phys. Lett. B 2003, 565, 61–75. [Google Scholar] [CrossRef]
  79. Bechtle, P.; Haber, H.E.; Heinemeyer, S.; Stål, O.; Stefaniak, T.; Weiglein, G.; Zeune, L. The Light and Heavy Higgs Interpretation of the MSSM. Eur. Phys. J. C 2017, 77, 67. [Google Scholar] [CrossRef] [Green Version]
  80. Domingo, F.; Heinemeyer, S.; Paßehr, S.; Weiglein, G. Decays of the neutral Higgs bosons into SM fermions and gauge bosons in the CP-violating NMSSM. Eur. Phys. J. C 2018, 78, 942. [Google Scholar] [CrossRef] [Green Version]
  81. Choi, K.; Im, S.H.; Jeong, K.S.; Park, C.B. Light Higgs bosons in the general NMSSM. Eur. Phys. J. C 2019, 79, 956. [Google Scholar] [CrossRef] [Green Version]
  82. Bartl, A.; Hirsch, M.; Vicente, A.; Liebler, S.; Porod, W. LHC phenomenology of the uvSSM. J. High Energy Phys. 2009, 5, 120. [Google Scholar] [CrossRef]
  83. Fidalgo, J.; Lopez-Fogliani, D.E.; Muñoz, C.; Ruiz de Austri, R. Neutrino Physics and Spontaneous CP Violation in the uvSSM. J. High Energy Phys. 2009, 8, 105. [Google Scholar] [CrossRef] [Green Version]
  84. Biekötter, T.; Grohsjean, A.; Heinemeyer, S.; Schwanenberger, C.; Weiglein, G. Possible Indications for New Higgs Bosons in the Reach of the LHC: N2HDM and NMSSM Interpretations. Available online: https://arxiv.org/abs/2109.01128 (accessed on 30 December 2021).
  85. Cao, J.; Guo, X.; He, Y.; Wu, P.; Zhang, Y. The diphoton signal of the light Higgs boson in Natural NMSSM. Phys. Rev. D 2017, 95, 116001. [Google Scholar] [CrossRef] [Green Version]
  86. Cao, J.; Jia, X.; Yue, Y.; Zhou, H.; Zhu, P. 96 GeV diphoton excess in seesaw extensions of the natural NMSSM. Phys. Rev. D 2020, 101, 055008. [Google Scholar] [CrossRef] [Green Version]
  87. Altmannshofer, W.; Gadam, S.A.; Gori, S.; Hamer, N. Explaining (g-2)u with Multi-TeV Sleptons. J. High Energy Phys. 2021, 118. [Google Scholar] [CrossRef]
  88. Ke, W.; Slavich, P. Higgs-Mass Constraints on a Supersymmetric Solution of the Muon g-2 Anomaly. Available online: https://arxiv.org/abs/2109.15277 (accessed on 30 December 2021).
  89. Biekötter, T.; Heinemeyer, S.; Muñoz, C. Precise prediction for the Higgs-Boson Masses in the uvSSM with three right-handed neutrino superfields. Eur. Phys. J. C 2019, 79, 667. [Google Scholar] [CrossRef] [Green Version]
  90. Miller, D.J.; Nevzorov, R.; Zerwas, P.M. The Higgs Sector of the Next-to-Minimal Supersymmetric Standard Model. Nucl. Phys. B 2004, 681, 3–30. [Google Scholar] [CrossRef] [Green Version]
  91. Heinemeyer, S.; Kpatcha, E.; Lara, I.; López-Fogliani, D.E.; Muñoz, C.; Nagata, N. The new (g-2)u result and the uvSSM. Eur. Phys. J. C 2021, 81, 802. [Google Scholar] [CrossRef]
  92. Biekötter, T.; Heinemeyer, S.; Weiglein, G. Vacuum (Meta-)Stability in the uvSSM. Available online: https://arxiv.org/abs/2112.12132 (accessed on 30 December 2021).
Figure 1. Renormalized Higgs boson mass in the tadpole-free MS ¯ scheme at full one-loop (dashed), full two-loop (dotted), and three-loop leading logarithms (dot-dashed red) levels as a function of the renormalization scale Q.
Figure 1. Renormalized Higgs boson mass in the tadpole-free MS ¯ scheme at full one-loop (dashed), full two-loop (dotted), and three-loop leading logarithms (dot-dashed red) levels as a function of the renormalization scale Q.
Particles 05 00006 g001
Figure 2. One-loop predictions of the Higgs boson masses m h (solid), m H (dashed), and m ± (dotted) as a function of the CP-odd Higgs boson mass m A . We have considered t a n β = 30 in a scenario where M S U S Y = 1 TeV and μ = 200 GeV.
Figure 2. One-loop predictions of the Higgs boson masses m h (solid), m H (dashed), and m ± (dotted) as a function of the CP-odd Higgs boson mass m A . We have considered t a n β = 30 in a scenario where M S U S Y = 1 TeV and μ = 200 GeV.
Particles 05 00006 g002
Figure 3. (Left) Predictions of the lightest Higgs boson mass, m h , at tree (dotted), one-loop (solid), and two-loop (dashed) levels as a function of the CP-odd Higgs boson mass m A . We have considered the m h m a x scenario ( X t 2.4 M S U S Y ) at the kinematical point t a n β = 30 , M S U S Y = 1 TeV , μ = 200 GeV and with a gluino mass M g ˜ = 1.5 TeV . (Right) The three-loop m h -predictions (dot-dashed red line) are included for large m A in the maximal stop mixing scenario.
Figure 3. (Left) Predictions of the lightest Higgs boson mass, m h , at tree (dotted), one-loop (solid), and two-loop (dashed) levels as a function of the CP-odd Higgs boson mass m A . We have considered the m h m a x scenario ( X t 2.4 M S U S Y ) at the kinematical point t a n β = 30 , M S U S Y = 1 TeV , μ = 200 GeV and with a gluino mass M g ˜ = 1.5 TeV . (Right) The three-loop m h -predictions (dot-dashed red line) are included for large m A in the maximal stop mixing scenario.
Particles 05 00006 g003
Figure 4. (Left) Dependence of m h on the supersymmetric scale M S U S Y for a vanishing stop mixing ( X t = 0 ) at full two-loop (black dashed line) and O( α t α s 2 ) three-loop levels (red dot-dashed line). The brown band is the CMS/ATLAS Higgs boson mass: 125.09 ± 0.24 GeV . (Right) Higgs boson mass predictions as a function of M S U S Y and X t . The brown region corresponds to the points compatible with a Higgs boson mass of 125.09 ± 5.0 GeV for t a n β = 1 . The gray region represents the values of X t and M S U S Y where λ ( M S U S Y ) becomes negative.
Figure 4. (Left) Dependence of m h on the supersymmetric scale M S U S Y for a vanishing stop mixing ( X t = 0 ) at full two-loop (black dashed line) and O( α t α s 2 ) three-loop levels (red dot-dashed line). The brown band is the CMS/ATLAS Higgs boson mass: 125.09 ± 0.24 GeV . (Right) Higgs boson mass predictions as a function of M S U S Y and X t . The brown region corresponds to the points compatible with a Higgs boson mass of 125.09 ± 5.0 GeV for t a n β = 1 . The gray region represents the values of X t and M S U S Y where λ ( M S U S Y ) becomes negative.
Particles 05 00006 g004
Figure 5. Estimation of the theory uncertainty of the Higgs mass produced with FeynHiggs in the hybrid approach for an MSSM scenario with degenerate SUSY masses.
Figure 5. Estimation of the theory uncertainty of the Higgs mass produced with FeynHiggs in the hybrid approach for an MSSM scenario with degenerate SUSY masses.
Particles 05 00006 g005
Figure 6. Self-energy topologies contributing to the three-loop corrections to the Higgs boson mass at order O( y t 6 ) in the SM. Topology 6 is the only non-planar self-energy.
Figure 6. Self-energy topologies contributing to the three-loop corrections to the Higgs boson mass at order O( y t 6 ) in the SM. Topology 6 is the only non-planar self-energy.
Particles 05 00006 g006
Table 1. MSSM spectrum fields and their S U ( 3 ) × S U ( 2 ) × U ( 1 ) quantum numbers. In this table, only one generation of leptons and quarks is exhibited. For each lepton and quark supermultiplet, there is a corresponding multiplet of scalar superpartners, while for the Higgs supermultiplet an antiparticle multiplet of charge-conjugated fermions is associated.
Table 1. MSSM spectrum fields and their S U ( 3 ) × S U ( 2 ) × U ( 1 ) quantum numbers. In this table, only one generation of leptons and quarks is exhibited. For each lepton and quark supermultiplet, there is a corresponding multiplet of scalar superpartners, while for the Higgs supermultiplet an antiparticle multiplet of charge-conjugated fermions is associated.
SupermultipletSuperfieldBosonsFermions SU ( 3 ) × SU ( 2 ) × U ( 1 )
gauge/ino V ^ 8 V ^ V ^ g W ± , W 0 B g ˜ W ˜ ± , W ˜ 0 B ˜ 8 1 0 1 3 0 1 1 0
s/lepton L ^ E ^ c v ˜ L , e ˜ L e ˜ R + v , e L e L c 1 2 ( 1 ) 1 1 2
s/quark Q ^ U ^ c D ^ c u ˜ L , d ˜ L u ˜ R * d ˜ R * u , d L u L c d L c 3 2 1 / 3 3 ¯ 1 ( 4 / 3 ) 3 ¯ 1 2 / 3
Higgs/ino H ^ 1 H ^ 2 H 1 0 , H 1 H 2 + , H 2 0 H ˜ 1 0 , H ˜ 1 H ˜ 2 + , H ˜ 2 0 1 2 ( 1 ) 1 2 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Reyes, E.; Fazio, R. High-Precision Calculations of the Higgs Boson Mass. Particles 2022, 5, 53-73. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5010006

AMA Style

Reyes E, Fazio R. High-Precision Calculations of the Higgs Boson Mass. Particles. 2022; 5(1):53-73. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5010006

Chicago/Turabian Style

Reyes, Edilson, and Raffaele Fazio. 2022. "High-Precision Calculations of the Higgs Boson Mass" Particles 5, no. 1: 53-73. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5010006

Article Metrics

Back to TopTop