The LHC Higgs Boson Discovery: Updated implications for Finite Unified Theories and the SUSY breaking scale

Finite Unified Theories (FUTs) are N = 1 supersymmetric Grand Unified Theories which can be made finite to all orders in perturbation theory, based on the principle of reduction of couplings. The latter consists in searching for renormalization group invariant relations among parameters of a renormalizable theory holding to all orders in perturbation theory. FUTs have proven very successful so far. In particular, they predicted the top quark mass one and half years before its experimental discovery, while around five years before the Higgs boson discovery a particular FUT was predicting the light Higgs boson in the mass range ~ 121 - 126 GeV, in striking agreement with the discovery at LHC. Here we review the basic properties of the supersymmetric theories and in particular finite theories resulting from the application of the method of reduction of couplings in their dimensionless and dimensionful sectors. Then we analyse the phenomenologically favoured FUT, based on SU(5). This particular FUT leads to a finiteness constrained version of the MSSM, which naturally predicts a relatively heavy spectrum with coloured supersymmetric particles above 2.7 TeV, consistent with the non-observation of those particles at the LHC. The electroweak supersymmetric spectrum starts below 1 TeV and large parts of the allowed spectrum of the lighter might be accessible at CLIC. The FCC-hh will be able to fully test the predicted spectrum.


Introduction
In 2012 the discovery of a new particle at the LHC was announced [1,2]. Within theoretical and experimental uncertainties the new particle is compatible with predictions for the Higgs boson of the Standard Model (SM) [3,4], constituting a milestone in the quest for understanding the physics of electroweak symmetry breaking (EWSB). However, taking the experimental results and the respective uncertainties into account, also many models beyond the SM can accomodate the data. Furthermore, the hierarchy problem, the neutrino masses, the Dark Matter, the over twenty free parameters of the model, just to name some questions, ask for a more fundamental theory to answer some, if not all, of those.
Therefore, one of the main aims of this fundamental theory is to relate these free parameters, or rephrasing it, to achieve a reduction of these parameters in favour of a smaller number (or ideally only one). This reduction is usually based in the introduction of a larger symmetry, rendering the theory more predictive. Very good examples are the Grand Unified Theories (GUTs) [5][6][7][8][9] and their supersymmetric extensions [10,11]. The case of minimal SU (5) is one example, where the number of couplings is reduced to one due to the corresponding unification. Data from LEP [12] suggested that a N = 1 global supersymmetry (SUSY) [10,11] is required in order for the prediction to be viable. Relations among the Yukawa couplings are also suggested in GUTs. For example, the SU (5) predicts the ratio of the tau to the bottom mass M τ /M b [13] in the SM. GUTs introduce, however, new complications such as the different ways of breaking this larger symmetry as well as new degrees of freedom.
A way to relate the Yukawa and the gauge sector, in other words achieving Gauge -Yukawa Unification (GYU) [14][15][16] seems to be a natural extension of the GUTs. The possibility that N = 2 SUSY [17] plays such a role is highly limited due to the prediction of light mirror fermions. Other phenomenological drawbacks appear in composite models and superstring theories.
The above approach seems to need SUSY as an essential ingredient. However the breaking of SUSY has to be understood too, since it provides the SM with several predictions for its free parameters. Actually, the RGI relation searches have been extended to the soft SUSY breaking (SSB) sector [25,[30][31][32] relating parameters of mass dimension one and two. This is indeed possible to be done on the RGI surface which is defined by the solution of the reduction equations.

Theoretical basis
In this section we outline the idea of reduction of couplings. Any RGI relation among couplings (which does not depend on the renormalization scale µ explicitly) can be expressed, in the implicit form Φ(g 1 , · · · , g A ) = const., which has to satisfy the partial differential equation (PDE) where β a is the β-function of g a . This PDE is equivalent to a set of ordinary differential equations, the so-called reduction equations (REs) [18,19,49], β g dg a dg = β a , a = 1, · · · , A , where g and β g are the primary coupling and its β-function, and the counting on a does not include g. Since maximally (A − 1) independent RGI "constraints" in the A-dimensional space of couplings can be imposed by the Φ a 's, one could in principle express all the couplings in terms of a single coupling g. However, a closer look to the set of Eqs. (2) reveals that their general solutions contain as many integration constants as the number of equations themselves. Thus, using such integration constants we have just traded an integration constant for each ordinary renormalized coupling, and consequently, these general solutions cannot be considered as reduced ones. The crucial requirement in the search for RGE relations is to demand power series solutions to the REs, which preserve perturbative renormalizability. Such an ansatz fixes the corresponding integration constant in each of the REs and picks up a special solution out of the general one. Remarkably, the uniqueness of such power series solutions can be decided already at the one-loop level [18,19,49]. To illustrate this, let us assume that the β-functions have the form where · · · stands for higher order terms, and β (1) bcd a 's are symmetric in b, c, d. We then assume that the ρ (n) a 's with n ≤ r have been uniquely determined. To obtain ρ (r+1) a 's, we insert the power series (3) into the REs (2) and collect terms of O(g 2r+3 ) and find d =g where the r.h.s. is known by assumption, and Therefore, the ρ (n) a 's for all n > 1 for a given set of ρ a 's can be uniquely determined if det M (n) d a = 0 for all n ≥ 0. As it will be clear later by examining specific examples, the various couplings in SUSY theories have the same asymptotic behaviour. Therefore, searching for a power series solution of the form (3) to the REs (2) is justified. This is not the case in non-SUSY theories, although the deeper reason for this fact is not fully understood.
The possibility of coupling unification described in this section is without any doubt attractive, because the "completely reduced" theory contains only one independent coupling, but it can be unrealistic. Therefore, one often would like to impose fewer RGI constraints, and this is the idea of partial reduction [50,51].

Reduction of dimension one and two parameters
The reduction of couplings was originally formulated for massless theories on the basis of the Callan-Symanzik equation [18,19]. The extension to theories with massive parameters is not straightforward if one wants to keep the generality and the rigor on the same level as for the massless case; one has to fulfill a set of requirements coming from the renormalization group equations, the Callan-Symanzik equations etc. along with the normalization conditions imposed on irreducible Green's functions [52]. There has been a lot of progress in this direction starting from ref. [25], where it was assumed that a massindependent renormalization scheme could be employed so that all the RG functions have only trivial dependencies on dimensional parameters and then the mass parameters were introduced similarly to couplings (i.e. as a power series in the couplings). This choice was justified later in [31,53] where the scheme independence of the reduction principle has been proven generally, i.e it was shown that apart from dimensionless couplings, pole masses and gauge parameters, the model may also involve coupling parameters carrying a dimension and masses. Therefore here, to simplify the analysis, we follow Ref. [25] and we, too, use a mass-independent renormalization scheme.
We start by considering a renormalizable theory which contain a set of (N + 1) dimension zero couplings, (ĝ 0 ,ĝ 1 , ...,ĝ N ), a set of L parameters with mass-dimension one, ĥ 1 , ...,ĥ L , and a set of M parameters with mass-dimension two, m 2 1 , ...,m 2 M . The renormalized irreducible vertex function Γ satisfies the RG equation where where µ is the energy scale, while β i are the β-functions of the various dimensionless couplings g i , Φ I are the various matter fields and γ m 2 α , γ h a and γ φI J are the mass, trilinear coupling and wave function anomalous dimensions, respectively (where I enumerates the matter fields). In a mass independent renormalization scheme, the γ's are given by γ m 2 ,ab α (g 0 , g 1 , ..., g N )ĥ aĥb , where γ h,b a , γ m 2 ,β α and γ m 2 ,ab α are power series of the g's (which are dimensionless) in perturbation theory.
The assumptions that the reduced theory is perturbatively renormalizable means that the functionsĝ i , f b a , e β α and k ab α , defined in (10), should be expressed as a power series in the primary coupling g:ĝ The above expansion coefficients can be found by inserting these power series into Eqs. (11), (12) and requiring the equations to be satisfied at each order of g. It should be noted that the existence of a unique power series solution is a non-trivial matter: It depends on the theory as well as on the choice of the set of independent parameters. It should also be noted that in the case that there are no independent mass-dimension one parameters (ĥ) the reduction of these terms take naturally the form where M is a mass-dimension one parameter which could be a gaugino mass which corresponds to the independent (gauge) coupling. In case, on top of that, there are no independent mass-dimension two parameters (m 2 ), the corresponding reduction takes analogous formm

Finiteness in N = 1 Supersymmetric Gauge Theories
Let us consider a chiral, anomaly free, N = 1 globally SUSY gauge theory based on a group G with gauge coupling constant g. The superpotential of the theory is given by where m ij and C ijk are gauge invariant tensors and the matter field φ i transforms according to the irreducible representation R i of the gauge group G. The renormalization constants associated with the superpotential (15), assuming that SUSY is preserved, are The N = 1 non-renormalization theorem [54][55][56] ensures that there are no mass and cubic-interaction-term infinities and therefore As a result the only surviving possible infinities are the wave-function renormalization constants Z j i , i.e. one infinity for each field. The one-loop β-function of the gauge coupling g is given by [57] where l(R i ) is the Dynkin index of R i and C 2 (G) is the quadratic Casimir invariant of the adjoint representation of the gauge group G. The β-functions of C ijk , by virtue of the non-renormalization theorem, are related to the anomalous dimension matrix γ ij of the matter fields φ i as: At one-loop level γ ij is [57] γ where C 2 (R) is the quadratic Casimir invariant of the representation R i , and C ijk = C * ijk . Since dimensional coupling parameters such as masses and couplings of cubic scalar field terms do not influence the asymptotic properties of a theory on which we are interested here, it is sufficient to take into account only the dimensionless SUSY couplings such as g and C ijk . So we neglect the existence of dimensional parameters and assume furthermore that C ijk are real so that C 2 ijk are always positive numbers. As one can see from Eqs. (20) and (22), all the one-loop β-functions of the theory vanish if β (1) The conditions for finiteness for N = 1 field theories with SU (N ) gauge symmetry are discussed in [58], and the analysis of the anomaly-free and no-charge renormalization requirements for these theories can be found in [59]. A very interesting result is that the conditions (23,24) are necessary and sufficient for finiteness at the two-loop level [57,[60][61][62][63].
In case SUSY is broken by soft terms, the requirement of finiteness in the one-loop soft breaking terms imposes further constraints among themselves [34]. In addition, the same set of conditions that are sufficient for one-loop finiteness of the soft breaking terms renders the soft sector of the theory two-loop finite [35].
The one-and two-loop finiteness conditions (23,24) restrict considerably the possible choices of the irreducible representations (irreps) R i for a given group G as well as the Yukawa couplings in the superpotential (15). Note in particular that the finiteness conditions cannot be applied to the minimal SUSY standard model (MSSM), since the presence of a U (1) gauge group is incompatible with the condition (23), due to C 2 [U (1)] = 0. This naturally leads to the expectation that finiteness should be attained at the grand unified level only, the MSSM being just the corresponding, low-energy, effective theory.
Another important consequence of one-and two-loop finiteness is that SUSY (most probably) can only be broken due to the soft breaking terms. Indeed, due to the unacceptability of gauge singlets, F-type spontaneous symmetry breaking [64] terms are incompatible with finiteness, as well as D-type [65] spontaneous breaking which requires the existence of a U (1) gauge group.
A natural question to ask is what happens at higher-loop orders. The answer is contained in a theorem [28,66] which states the necessary and sufficient conditions to achieve finiteness at all orders. Before we discuss the theorem let us make some introductory remarks. The finiteness conditions impose relations between gauge and Yukawa couplings. To require such relations which render the couplings mutually dependent at a given renormalization point is trivial. What is not trivial is to guarantee that relations leading to a reduction of the couplings hold at any renormalization point. As we have seen, the necessary and also sufficient condition for this to happen is to require that such relations are solutions to the REs and hold at all orders. Remarkably, the existence of all-order power series solutions to (25) can be decided at one-loop level, as already mentioned.
Let us now turn to the all-order finiteness theorem [28,66], which states the conditions under which an N = 1 SUSY gauge theory can become finite to all orders in the sense of vanishing β-functions, that is of physical scale invariance. It is based on (a) the structure of the supercurrent in N = 1 SUSY gauge theory [67][68][69], and on (b) the non-renormalization properties of N = 1 chiral anomalies [27,28,66,70,71]. Details on the proof can be found in refs. [28,66] and further discussion in Refs. [27,29,[70][71][72]. Here, following mostly Ref. [72] we present a comprehensible sketch of the proof.
Consider an N = 1 SUSY gauge theory, with simple Lie group G. The content of this theory is given at the classical level by the matter supermultiplets S i , which contain a scalar field φ i and a Weyl spinor ψ ia , and the vector supermultiplet V a , which contains a gauge vector field A a µ and a gaugino Weyl spinor λ a α . Let us first recall certain facts about the theory: (1) A massless N = 1 SUSY theory is invariant under a U (1) chiral transformation R under which the various fields transform as follows The corresponding axial Noether current is conserved classically, while in the quantum case is violated by the axial anomaly From its known topological origin in ordinary gauge theories [73][74][75], one would expect the axial vector current J µ R to satisfy the Adler-Bardeen theorem and receive corrections only at the one-loop level. Indeed it has been shown that the same non-renormalization theorem holds also in SUSY theories [27,70,71]. Therefore (2) The massless theory we consider is scale invariant at the classical level and, in general, there is a scale anomaly due to radiative corrections. The scale anomaly appears in the trace of the energy momentum tensor T µν , which is traceless classically. It has the form (3) Massless, N = 1 SUSY gauge theories are classically invariant under the SUSY extension of the conformal group -the superconformal group. Examining the superconformal algebra, it can be seen that the subset of superconformal transformations consisting of translations, SUSY transformations, and axial R transformations is closed under SUSY, i.e. these transformations form a representation of SUSY. It follows that the conserved currents corresponding to these transformations make up a supermultiplet represented by an axial vector superfield called the supercurrent J, where J µ R is the current associated to R invariance, Q µ α is the one associated to SUSY invariance, and T µ ν the one associated to translational invariance (energy-momentum tensor).
The anomalies of the R current J µ R , the trace anomalies of the SUSY current, and the energy-momentum tensor, form also a second supermultiplet, called the supertrace anomaly where T µ µ is given in Eq. (30) and (4) It is very important to note that the Noether current defined in (27) is not the same as the current associated to R invariance that appears in the supercurrent J in (31), but they coincide in the tree approximation. So starting from a unique classical Noether current J µ R(class) , the Noether current J µ R is defined as the quantum extension of J µ

R(class)
which allows for the validity of the non-renormalization theorem. On the other hand J µ R , is defined to belong to the supercurrent J, together with the energy-momentum tensor. The two requirements cannot be fulfilled by a single current operator at the same time.
Although the Noether current J µ R which obeys (28) and the current J µ R belonging to the supercurrent multiplet J are not the same, there is a relation [28,66] between quantities associated with them where r was given in Eq. (29). The r A are the non-renormalized coefficients of the anomalies of the Noether currents associated to the chiral invariances of the superpotential, and -like r-are strictly one-loop quantities. The γ A 's are linear combinations of the anomalous dimensions of the matter fields, and x g , and x ijk are radiative correction quantities. The structure of equality (34) is independent of the renormalization scheme. One-loop finiteness, i.e. vanishing of the β-functions at one loop, implies that the Yukawa couplings C ijk must be functions of the gauge coupling g. To find a similar condition to all orders it is necessary and sufficient for the Yukawa couplings to be a formal power series in g, which is solution of the REs (25).
We can now state the theorem for all-order vanishing β-functions.

Theorem:
Consider an N = 1 SUSY Yang-Mills theory, with simple gauge group. If the following conditions are satisfied: 1. There is no gauge anomaly.
2. The gauge β-function vanishes at one loop 3. There exist solutions of the form to the conditions of vanishing one-loop matter fields anomalous dimensions 4. These solutions are isolated and non-degenerate when considered as solutions of vanishing one-loop Yukawa β-functions: Then, each of the solutions (36) can be uniquely extended to a formal power series in g, and the associated super Yang-Mills models depend on the single coupling constant g with a β function which vanishes at all orders.
It is important to note a few things: The requirement of isolated and non-degenerate solutions guarantees the existence of a unique formal power series solution to the reduction equations. The vanishing of the gauge β function at one loop, β (1) g , is equivalent to the vanishing of the R current anomaly (28). The vanishing of the anomalous dimensions at one loop implies the vanishing of the Yukawa couplings β functions at that order. It also implies the vanishing of the chiral anomaly coefficients r A . This last property is a necessary condition for having β functions vanishing at all orders. 1

Proof:
Insert β ijk as given by the REs into the relationship (34) between the axial anomalies coefficients and the β-functions. Since these chiral anomalies vanish, we get for β g an homogeneous equation of the form The solution of this equation in the sense of a formal power series in is β g = 0, order by order. Therefore, due to the REs (25), β ijk = 0 too. Thus we see that finiteness and reduction of couplings are intimately related. Since an equation like Eq.(34) is lacking in non-SUSY theories, one cannot extend the validity of a similar theorem in such theories.

The SSB sector of reduced N = 1 SUSY and Finite Theories
As we have seen in subsection 2.1, the method of reducing the dimensionless couplings has been extended [25], to the soft SUSY breaking (SSB) dimensionful parameters of N = 1 SUSY theories. In addition it was found [40] that RGI SSB scalar masses in GYU models satisfy a universal sum rule.
Consider the superpotential given by (15) along with the Lagrangian for SSB terms where the φ i are the scalar parts of the chiral superfields Φ i , λ are the gauginos and M their unified mass. We assume that the reduction equations admit power series solutions of the form If we knew higher-loop β-functions explicitly, we could follow the same procedure and find higher-loop RGI relations among SSB terms. However, the β-functions of the soft scalar masses are explicitly known only up to two loops. In order to obtain higher-loop results some relations among β-functions are needed.
In the case of finite theories we assume that the gauge group is a simple group and the one-loop β-function of the gauge coupling g vanishes. According to the finiteness theorem of Refs. [28,66] and the assumption given in (40), the theory is then finite to all orders in perturbation theory, if, among others, the one-loop anomalous dimensions γ j(1) i vanish. The one-and two-loop finiteness for h ijk can be achieved by [35] where . . . stand for higher order terms. Now, to obtain the two-loop sum rule for soft scalar masses, we assume that the lowest order coefficients ρ ijk (0) and also (m 2 ) i j satisfy the diagonality relations ρ ipq(0) ρ jpq (0) ∝ δ j i for all p and q and (m 2 ) i j = m 2 j δ i j , respectively. Then we find the following soft scalar-mass sum rule [16,42,77] ( for i, j, k with ρ ijk (0) = 0, where ∆ (2) is the two-loop correction which vanishes for the universal choice in accordance with the previous findings of Ref. [35] (in the above relation T (R l ) is the Dynkin index of the R l irrep). Making use of the spurion technique [56,[78][79][80][81], it is possible to find the following all-loop relations among SSB β-functions, [82][83][84][85][86][87] where (γ 1 ) i j = Oγ i j , C lmn = (C lmn ) * , and Assuming, following [83], that the relation among couplings is all-loop RGI and using the all-loop gauge β-function of Novikov et al. [43][44][45] given by the all-loop RGI sum rule [41] was found: In addition the exact-β-function for m 2 in the NSVZ scheme has been obtained [41] for the first time and is given by Surprisingly enough, the all-loop result (53) coincides with the superstring result for the finite case in a certain class of orbifold models [42] if d ln C ijk /d ln g = 1.

All-loop RGI Relations in the SSB Sector
Let us now see how the all-loop results on the SSB β-functions, Eqs. (45)-(50), lead to all-loop RGI relations. We make two assumptions: (a) the existence of an RGI surface on which C = C(g), or equivalently that holds, i.e. reduction of couplings is possible, and (b) the existence of an RGI surface on which holds, too, in all orders. Then one can prove, [88,89], that the following relations are RGI to all loops (note that in both (a) and (b) assumptions above we do not rely on specific solutions of these equations) where M 0 is an arbitrary reference mass scale to be specified shortly.
The assumption that for an RGI surface F (g, C ijk , C * ijk ) leads to where Eq.(55) has been used. Now let us consider the partial differential operator O in Eq.(48) which, assuming Eq.(51) becomes In turn, β M given in Eq. (45) becomes which by integration provides us [88,90] where M 0 is the integration constant and can be associated to the unification scale M U in GUTs or to the gravitino mass m 3/2 in a supergravity framework. Therefore, Eq. (65) becomes the all-loop RGE Eq. (57). Note that β M using Eqs. (64) and (65) can be written as Similarly Next, form Eq. (51) and (65) we obtain while β ijk h , given in Eq.(46) and using Eq.(67), becomes [88] which shows that Eq.(68) is all-loop RGI. In a similar way Eq. (59) can be shown to be all-loop RGI. Finally, we would like to emphasize that under the same assumptions (a) and (b) the sum rule given in Eq. (53) has been proven [41] to be all-loop RGI, which (using Eq.(65)) gives us a generalization of Eq.(60) to be applied in considerations of non-universal soft scalar masses, which are necessary in many cases. Moreover, the sum rule holds also in the more general cases, discussed in subsection 2.1, according to which exact relations among the squared scalar and gaugino masses can be found.

A successful Finite Unified Theory
We review an all-loop FUT with SU (5) as gauge group, where the reduction of couplings has been applied to the third generation of quarks and leptons. This FUT was selected previously on the basis of agreement with the known experimental data at the time [33] and was predicting the Higgs mass to be in the range 121 − 126 GeV almost five years before the discovery. The particle content of the model we will study, which we denote SU (5)-FUT consists of the following supermultiplets: three (5 + 10), needed for each of the three generations of quarks and leptons, four (5 + 5) and one 24 considered as Higgs supermultiplets. When the gauge group of the finite GUT is broken the theory is no longer finite, and we will assume that we are left with the MSSM [15,[20][21][22][23][24].
A predictive GYU SU (5) model which is finite to all orders, in addition to the requirements mentioned already, should also have the following properties: 1. One-loop anomalous dimensions are diagonal, i.e., γ 3. The two Higgs doublets of the MSSM should mostly be made out of a pair of Higgs quintet and anti-quintet, which couple to the third generation.
After the reduction of couplings the symmetry is enhanced, leading to the following superpotential [42,91] The non-degenerate and isolated solutions to γ (1) i = 0 give us: and from the sum rule we obtain: i.e., in this case we have only two free parameters m 10 and M for the dimensionful sector.
As already mentioned, after the SU (5) gauge symmetry breaking we assume we have the MSSM, i.e. only two Higgs doublets. This can be achieved by introducing appropriate mass terms that allow to perform a rotation of the Higgs sector [20,24,[92][93][94], in such a way that only one pair of Higgs doublets, coupled mostly to the third family, remains light and acquires vacuum expectation values. To avoid fast proton decay the usual fine tuning to achieve doublet-triplet splitting is performed, although the mechanism is not identical to minimal SU (5), since we have an extended Higgs sector.
Thus, after the gauge symmetry of the GUT theory is broken, we are left with the MSSM, with the boundary conditions for the third family given by the finiteness conditions, while the other two families are not restricted.

Phenomenological constraints
In this section we briefly review the relevant experimental constraints that we apply in our phenomenological analysis.

Flavour Constraints
We consider four types of flavour constraints to apply to the SU(5)-FUT, where SUSY is known to have significant impact. More specifically, we consider the flavour observables BR(b → sγ), BR(B s → µ + µ − ), BR(B u → τ ν) and ∆M Bs . 2 The uncertainties are the linear combination of the experimental error and twice the theoretical uncertainty in the MSSM. 3 In the case that no specific estimate is available, we use the SM uncertainty.
For the branching ratio BR(b → sγ), we take a value from the Heavy Flavour Averaging Group (HFAG) [95,96] BR For the branching ratio BR(B s → µ + µ − ), a combination of CMS and LHCb data [97][98][99][100][101] is used For the B u decay to τ ν we use the limit [96,102,103] while for ∆M Bs [104,105] ∆M exp At the end of the phenomenological discussion we also comment on the cold dark matter (CDM) density. It is well known that the lightest neutralino, being the lightest SUSY particle (LSP), is an excellent candidate for CDM [106]. Consequently one can in principle demand that the lightest neutralino is indeed the LSP and parameters leading to a different LSP could be discarded.
The current bound, favoured by a joint analysis of WMAP/Planck and other astrophysical and cosmological data, is at the 2 σ level given by the range [107,108],

The light Higgs boson mass
The quartic couplings in the Higgs potential are given by the SM gauge couplings. As a consequence, the lightest Higgs mass is not a free parameter, but rather predicted in terms of the other parameters of the model. Higher order corrections are crucial for a precise predictions of M h , see Refs. [109][110][111] for reviews. The discovery of a Higgs boson at ATLAS and CMS in July 2012 [112,113] can be interpreted as the discovery of the light CP-even Higgs boson of the MSSM Higgs spectrum [114][115][116]. The experimental average for the (SM) Higgs boson mass is taken to be [117] M exp and adding in quadrature a 3 (2) GeV theory uncertainty [118][119][120] for the Higgs mass caclulation in the MSSM we arrive at an allowing range of We used the code FeynHiggs

Numerical analysis
In this section we will analyse the particle spectrum predicted in the SU (5) The bottom mass is calculated at M Z to avoid uncertainties that come from running down to the pole mass; the leading SUSY radiative corrections to the bottom and tau masses have been taken into account [123]. We use the following value for the bottom mass at The bounds on the m b (M Z ) and the m t mass clearly single out µ < 0, as the solution most compatible with these experimental constraints.
As was already mentioned, for the lightest Higgs boson mass we used the code FeynHiggs (2.14.0 beta). The prediction for M h of SU (5)-FUT with µ < 0 is shown in Fig.2, in a range where the unified gaugino mass varies from 0.5 TeV M 9 TeV. The green points include the B-physics constraints. One should keep in mind that these predictions are subject to a theory uncertainty of 3(2) GeV [118]. Older analysis, including in particular less refined evaluations of the light Higgs boson mass, are given in Refs. [46,124,125].
The allowed values of the Higgs mas put a limit on the allowed values of the SUSY masses, as can be seen in Fig.3. In the left (right) plot we impose M h = 125.1 ± 3.1(2.1) GeV as discussed above. In particular very heavy coloured SUSY particles are favoured (nearly independent of the M h uncertainty), in agreement with the nonobservation of those particles at the LHC [126]. Overall, the allowed coloured SUSY masses would remain unobservable at the (HL-)LHC, the ILC or CLIC. However, the coloured spectrum would be accessible at the FCC-hh [127], as could the full heavy Higgs boson spectrum. On the other hand, the lightest observable SUSY particle (LSOP) is the Figure 3: The left (right) plot shows the spectrum of the SU (5)-FUT (with µ < 0) model after imposing the constraint M h = 125.1 ± 3.1(2.1) GeV. The light (green) points are the various Higgs boson masses, the dark (blue) points following are the two scalar top and bottom masses, the gray ones are the gluino masses, then come the scalar tau masses in orange (light gray), the darker (red) points to the right are the two chargino masses followed by the lighter shaded (pink) points indicating the neutralino masses. scalar tau. Some parts of the allowed spectrum of the lighter scalar tau or the lighter charginos/neutralinos might be accessible at CLIC with √ s = 3 TeV. In Table 1 we show two example spectra of the SU (5)-FUT (with µ < 0) which span the mass range of the parameter space that is in agreement with the B-physics observables and the Higgs-boson mass measurement. We give the lightest and the heaviest spectrum for δM h = 2.1 and δM h = 3. We find that no point of SU (5)-FUT (with µ < 0) fulfills the strict bound of Eq. (77). (For our evaluation we have used the code MicroMegas [128][129][130].) Consequently, on a more general basis a mechanism is needed in our model to reduce the CDM abundance in the early universe. This issue could, for instance, be related to another problem, that of neutrino masses. This type of masses cannot be generated naturally within the class of finite unified theories that we are considering in this paper, although a non-zero value for neutrino masses has clearly been established [103]. However, the class of FUTs discussed  Table 1: Two example spectra of the SU (5)-FUT (with µ < 0) . All masses are in GeV and rounded to 1 (0.1) GeV (for the light Higgs mass).
here can, in principle, be easily extended by introducing bilinear R-parity violating terms that preserve finiteness and introduce the desired neutrino masses [131][132][133]. R-parity violation [134] would have a small impact on the collider phenomenology presented here (apart from the fact that SUSY search strategies could not rely on a 'missing energy' signature), but remove the CDM bound of Eq. (77) completely. The details of such a possibility in the present framework attempting to provide the models with realistic neutrino masses will be discussed elsewhere. Other mechanisms, not involving R-parity violation (and keeping the 'missing energy' signature), that could be invoked if the amount of CDM appears to be too large, concern the cosmology of the early universe. For instance, "thermal inflation" [135] or "late time entropy injection" [136] could bring the CDM density into agreement with the WMAP measurements. This kind of modifications of the physics scenario neither concerns the theory basis nor the collider phenomenology, but could have a strong impact on the CDM bounds. (Lower values than those permitted by Eq. (77) are naturally allowed if another particle than the lightest neutralino constitutes CDM.)

Conclusions
The MSSM is considered a very attractive candidate for describing physics beyond the SM. However, the serious problem of the SM having too many free parameters is further proliferated in the MSSM. Assuming a GUT beyond the scale of gauge coupling unification, based on the idea that a Particle Physics Theory should be more symmetric at higher scales, seems to fit to the MSSM. On the other hand, the unification scenario seems to be unable to further reduce the number of free parameters. Attempting to reduce the free parameters of a theory, a new approach was proposed in Refs. [18,19] based on the possible existence of RGI relations among couplings. Although this approach could uncover further symmetries, its application opens new horizons too. At least the Finite Unified Theories seem to be a very promising field for applying the reduction approach. In the FUT case, the discovery of RGI relations among couplings above the unification scale ensures at the same time finiteness to all orders.
The discussion in the previous sections of this paper shows that the predictions of the particular FUT discussed here are impressive. In addition, one could add some comments on a successful FUT from the theoretical side, too. The developments on treating the problem of divergencies include string and non-commutative theories, as well as N = 4 SUSY theories [137,138], N = 8 supergravity [139][140][141][142][143] and the AdS/CFT correspondence [144]. It is very interesting that the N = 1 FUT discussed here includes many ideas which survived phenomenological and theoretical tests as well as the ultraviolet divergence problem. It is actually solving that problem in a minimal way.
We concentrated our examination on the predictions of one particular SU (5) Finite Unified Theory, including the restrictions of third generation quark masses and B-physics observables. The model, SU (5)-FUT (with µ < 0), is consistent with all the phenomenological constraints. Compared to our previous analyses [46,47,124,125], the improved evaluation of M h prefers a heavier (Higgs) spectrum and thus in general allows only a very heavy SUSY spectrum. The coloured spectrum could easily escape the (HL-)LHC searches, but can likely be tested at the FCC-hh. The lower part of the electroweak spectrum could be accessible at CLIC.