SUSY threshold corrections to quark and lepton mixing inspired by SO(10) GUT models

,


Introduction
Observing Higgs boson [1,2] establishes the standard model (SM) completely. However, some topics such as non-zero neutrino masses and existence of the dark matter are not explained in the SM, and then, many physics beyond the SM (bSMs) are established to explain these topics. So far, many experiments have tried to find the signals of the new physics and they reported the results about upper and/or lower limits or observed values for some observables. We can now test the bSMs by using these limits and observables.
One of the most important tools to test the bSMs is flavor phenomena. Flavor observables can probe effects of new particles with masses far above the energy scale of current collider experiments. In the SM, the only origin of flavor violations (FVs) is the Cabibbo-Kobayashi-Maskawa (CKM) matrix [3,4] which is defined by a product of the mixing matrices of left-handed quarks. Small mixing angles of the CKM matrix cause small FVs for the SM predictions, and therefore, flavor observables show us clear signal of the bSMs.
In general, there are additional origins of the FVs in the bSMs. Especially, the diagonalizing matrices play a crucial role among these origins. Well-known examples are supersymmetric (SUSY) models. In these models, FVs can be induced by mediating the SUSY particles. To predict flavor observables via contributions from SUSY particles, mass insertion parameters (MIPs) in a super-CKM basis [5] are often used. The super-CKM basis is a basis in which the rotations for obtaining diagonal Yukawa matrices are applied to sfermions as well as the SM fermions. Because of this rotation, MIPs include information JHEP02(2020)103 of the corresponding diagonalizing matrices. However, most of the diagonalizing matrices are unknown ones because what we know about these matrices from the experiments is only the CKM matrix. 1 One of the possibilities to know the information of these unknown diagonalizing matrices is considering the model with grand unifications. In grand unified theories (GUTs), the three gauge interactions in the SM are unified into a single one, and moreover, the particles in the SM are unified into fewer multiplets [6]. The latter unification reduces degrees of freedom of the diagonalizing matrices. In the GUTs, new particles are predicted, and they also induce FVs through unknown diagonalizing matrices. Although these contributions are negligible for most processes because new particles which induce FVs have masses around the GUT scale, proton decays have a sensitivity for these contributions. Recently, the Super-Kamiokande Collaboration reported some events which can be candidates for proton decay, and these include a candidate for FV decay of the proton such as p → πµ decay mode [7,8]. It has been confirmed that in some models, this FV decay of the proton can be comparable to main decay mode of the proton [9][10][11][12].
In the minimal SO(10) GUT models [13,14], 2 all quarks and leptons are unified into only one multiplet 16 for one family. Because of this unification, all diagonalizing matrices become the CKM-like matrix to realize the CKM matrix except for a light neutrino mixing matrix. Note that the mixing angles of the light neutrino mixing matrix can be large because of new degrees of freedom from seesaw mechanism [15][16][17][18][19][20][21][22][23][24]. However because of the unification, it is not easy to realize measured quark and lepton masses and mixings. One useful correction to solve this problem is a SUSY threshold correction [5,[25][26][27].
In this paper, we discuss how large diagonalizing matrices originated from the SUSY threshold corrections are allowed in the minimal SUSY SO(10) GUT model. In this setup, we may obtain some specific signal for the model from the effect of large diagonalizing matrices. Note that large off-diagonal parts of fermion mass matrices which generate large diagonalizing matrices come from large off-diagonal parts of the A-terms through the SUSY threshold corrections. However, these off-diagonal parts are strongly constrained by flavor phenomena like SUSY flavor-changing neutral currents (FCNCs). 3 Furthermore, these also 1 Although the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix is also the experimental observable, we focus only on the correction to the CKM matrix in this paper. This is because the corrections considered in this paper are expected to be small, and then, such small corrections are irrelevant for the PMNS matrix which has large mixing angles. 2 In this paper, we define the models in which only 16-dimensional representation is used for quarks and leptons as minimal SO(10) GUT models. In addition, we use only one 10-dimensional representation for two Higgses in SUSY models. 3 The effects of these chirally enhanced corrections on FCNC processes have been studied, for example in ref. [44].

JHEP02(2020)103
cause the negative mass squared for sfermions through the renormalization group equation (RGE) effects. Then, we calculate the maximal size of diagonalizing matrices numerically without conflicting with these concerns. This paper is organized as follow. In section 2, we show the notation used in this paper. We give the method of our calculation in section 3, and our main results are shown in section 4. The conclusion is in section 5. In addition, we show the detail of the SUSY threshold corrections in appendix A. We also summarize the relevant constraints to our calculation in appendix B. The discussion about the maximal SUSY threshold correction which satisfy the constraints are shown in appendix C.

Notation
In this section, we summarize our notation in this paper. The superpotential for the Yukawa interactions and µ term in the interaction basis is and H d are superfields of the left-handed quark doublet, left-handed lepton doublet, right-handed up-type quark, right-handed down-type quark, right-handed charged-lepton, up-type Higgs doublet, and down-type Higgs doublet, respectively. In this paper, we use the capital letters for superfields and the small letters for the corresponding SM components. Here, Y u , Y d , and Y e are the up-type quark, downtype quark, and charged lepton Yukawa matrices, a, b = 1, 2 are weak isospin indices, i, j = 1, 2, 3 are generation indices, and 12 ≡ +1 is an antisymmetric tensor. To obtain fermion masses, neutral components of CP-even Higgses develop vacuum expectation values (VEVs) h 0 u(d) ≡ v u(d) whose ratio is defined as tan β ≡ v u /v d . The mass basis of fermions are obtained by diagonalizing Yukawa matrices as where ψ and ψ indicate the SM fermion in the interaction basis and mass basis, L ψ and R ψ are diagonalizing matrices for the left-handed and right-handed fermion, and (Ŷ ψ ) ii is the Yukawa coupling of i-th generation fermion ψ. The same rotations are applied to the corresponding sfermions because of the SUSY, and this basis is called the super-CKM basis. In this paper,X refers to the parameter of X in the super-CKM basis. We calculate mixing angles for the diagonalizing matrices, defined as

JHEP02(2020)103
To calculate the SUSY threshold corrections, we use following soft-breaking parameters in a Lagrangian L soft , where M i (i = 1, 2, 3) are gaugino masses, andψ is the SUSY partner for the corresponding SM particle ψ. Moreover, we define following 6×6 sfermion mass squared matrices as (2.9) The sfermion mass squared matrixM 2 ψ ij is diagonalized by a unitary matrix Uψ aŝ The SUSY threshold corrections are induced by chirality-flipping self-energies. In ref. [45], the details of the calculation for these corrections are shown in a decoupling limit M SUSY v where M SUSY and v are the SUSY scale and the VEV of the SM Higgs. We also apply to this decoupling limit in this paper. In the appendix A, we summarize the SUSY threshold corrections in our notation.
The SUSY threshold corrections affect on the diagonalizing matrices if chirality-flipping self-energies have off-diagonal parts which come from the off-diagonal parts of A-terms and sfermion mass squared matrices. The size of these off-diagonal parts are important to obtain the large SUSY threshold corrections. However, these off-diagonal parts are also related to the SUSY FCNCs, and therefore, these sizes are constrained from the experimental measurements. Moreover, large off-diagonal parts of A-terms make the sfermion squared mass negative and cause the a charge and color breaking (CCB) [46][47][48][49][50][51]. In the appendix B, we summarize the SUSY FCNC and CCB constraints, used in this paper. which is needed to realize the electroweak symmetry breaking (EWSB). Note that for simplicity and for reducing the parameters, we assume ∆m 2 h ∼ M 2 SUSY . We discuss the Yukawa unification at the beginning 5 since we can fix part of the input parameters from it.

SUSY threshold correction for realizing Yukawa unification
The SUSY threshold correction plays an important role for realizing the Yukawa unification in the SUSY SO(10) GUT models. When we consider this Yukawa unification, we can fix the unified gauge coupling g GUT and (3,3) component of unified Yukawa matrix (Y GUT ) 33 . The latter one is obtained by setting the values of the SUSY threshold corrections and tan β appropriately. Note that some of the other parameters can be fixed by the other assumptions, as we explain later.
First of all, we find the required SUSY threshold correction and tan β for realizing the Yukawa unification. In this calculation, we start from the central values of the experimental measurements summarized in table 1. There are three scales which are relevant to our calculation: M Z , M SUSY , and the GUT scale M GUT . First, we calculate the parameters at the M Z scale from the input parameters in table 1. For calculating the running quark masses, we use the Mathematica package RunDec [54], while for the running lepton masses, we translate the lepton pole masses to running masses at the M Z scale, following ref. [55]. Next, we derive the couplings at the SUSY scale from those at the M Z scale, using the SM RGE at the two-loop level [56][57][58][59]. 6 At the SUSY scale, we translate the MS scheme into the DR scheme according to ref. [64]. After applying the matching condition by using the tan β, we add the SUSY threshold corrections as where ∆y ψ shows the SUSY threshold correction for the fermion ψ. However in this section, we only consider ∆y b because this contribution is dominant over the other ones. Here (0) means parameters below the SUSY scale, namely, SM parameters. Finally, we derive the couplings at the GUT scale from those at the SUSY scale, using the minimal SUSY SM (MSSM) RGE at the two-loop level [65]. Note that in our calculation, M GUT and g GUT are defined as where g 1 is the GUT normalized gauge coupling for the U(1) Y , and g 2 is the gauge coupling for the SU(2) L . By means of this procedure, we can obtain the required value of ∆y b and tan β for realizing the Yukawa unification: y t = y b = y τ .
In figure 1, we show the dependence of ∆y b and tan β on M SUSY . The blue line and axis show the required SUSY threshold correction, and the red ones show the tan β. From this JHEP02(2020)103 figure, importance of large tan β is clear which has been already discussed in refs. [66][67][68][69][70]. Interestingly, it is shown that a maximal SUSY threshold correction is larger than the required one because of the large µ contribution when the SUSY scale is below 100 TeV. In appendix C, we discuss this issue by considering the CCB and EWSB constraints. Note that the SUSY threshold correction for the diagonal part is dominated by the contributions from µ parameter in the case of large tan β. Therefore in this paper, the µ parameter is fixed to realize the required SUSY threshold correction, and we use a maximal A-term to maximize the effect to the diagonalizing matrices.

SUSY threshold correction for the diagonalizing matrices
In this section, we explain how to calculate the SUSY threshold correction for the diagonalizing matrices. In this calculation, we use the eigenvalues of Yukawa matrices at the GUT scale. Therefore, when the SUSY threshold corrections are absent, the Yukawa matrices are diagonal matrices at any energy scales, and hence, the corresponding diagonalizing matrices are unit matrices. There are some advantages in this setup. First one is that contributions from the SUSY threshold corrections to the diagonalizing matrices become clear. Origins of the SUSY threshold corrections to the diagonalizing matrices are off-diagonal parts of SUSY parameters A ψ and m 2 ψ in the super-CKM basis. If the Yukawa matrices have the off-diagonal elements, these parts become the combination of those in the interaction basis and the diagonalizing matrices, which leads to the contributions are unclear. The second advantage is that we can exclude some unfavorable contributions to the SUSY threshold corrections. In many cases, family symmetries which show us how to realize the CKM matrix are introduced with the SO(10) symmetry. However for this calculation, these family symmetries are unfavorable because these affect and fix the diagonalizing matrices. By using diagonal Yukawa matrices, we can remove these unfavorable contributions.
In our calculation, we use the following assumptions at the GUT scale: where δA indicates the off-diagonal part of the A-terms. Note that if δA = 0, each diagonalizing matrix becomes the unit matrix. Since we set ∆m 2 h ∼ M 2 SUSY , there are 5 parameters in this setup, M 1/2 , A 0 , δA, m 16 , and m 10 . In order to reduce the parameters, we take the following assumptions: (i) the gluino mass at the SUSY scale is around M SUSY ; (ii) the SUSY scale is defined by two stop masses as Therefore, once we fix A 0 , δA, and m 10 in each M SUSY case, we can determine M 1/2 and m 16 by (i) and (ii), respectively. The remaining parameters can be determined so that the SUSY FCNC and CCB constraints are satisfied and the required SUSY threshold correction are obtained. Note that we also check whether the EWSB can be realized. Among such parameter sets, we take the set which maximizes the A-terms and calculate the diagonalizing matrices.

JHEP02(2020)103
After we calculate the diagonalizing matrices, we extract mixing angles by following eqs. (2.3) and (2.4). Moreover we compare these mixing angles and those of the measured CKM matrices [53], which are where the CKM matrix V CKM is defined as If the diagonalizing matrices are as large as or larger than the CKM matrix, the SUSY threshold corrections affect on the CKM matrix. In that case, the CKM matrix is obtained by the combination of the diagonalizing matrices from the SUSY threshold corrections and Yukawa matrices, or some cancellations are needed between the diagonalizing matrices to realize the CKM matrix. Since s CKM 13 is smaller than the other elements, we calculate s ψL,R 13 by setting δA as In this case, we should take care about the SUSY FCNCs, especially B d -B d mixing and τ → eγ. The details of calculations and constraints for these processes are summarized in the appendix B. In the other case of δA, the contributions to the corresponding diagonalizing matrices are expected to be small due to the largeness of s CKM ij and/or the SUSY FCNC constraints.

Results
The input parameters which satisfies the constraints and realizes the required SUSY threshold correction simultaneously are shown in table 2. Here, A max 0 is the maximal value of A 0 which satisfies the CCB constraints. Note that the value of ∆m 2 h for M SUSY = 10 TeV should be smaller than M 2 SUSY . We will explain this issue later. By using these input parameters, we can calculate the mixing angles for L ψ and R ψ . We show the results in table 3. From these results, we found that s qL,R

13
(q = u, d) can be as large as s CKM 13 when M SUSY becomes large. This is because the SUSY FCNC constraints become weak in large M SUSY case, and therefore, we can take large A off . In this case, however, a strong cancellation between |m 2 hu | and m 2 h d is needed to realize the EWSB. In order to specify that our input parameters realize the maximal size of s ψL,R 13 , 7 we show the parameter spaces for s dL 13 in each M SUSY case in figures 2 and 3. We checked that 7 The "maximal size" of the mixing angle in this paper means the mixing angle which is calculated by the maximal value of the off-diagonal element of A-terms A off through the SUSY threshold correction. Obviously, the maximal value of A off can be calculated by the SUSY FCNC and the CCB constraints.    We found that such parameter space can be obtained when ∆m 2 h 0.48M 2 SUSY , and then, we set ∆m 2 h = 0.48M 2 SUSY . On the other hand, we can obtain the parameter space in M SUSY = 2, 100 TeV case even when ∆m 2 h = M 2 SUSY since one of the EWSB and SUSY FCNC constraints is irrelevant in these cases. constraints become severe because of the large values of |m 2 hu | and |m 2 h d |. This means that a strong cancellation between |m 2 hu | and m 2 h d is needed, and therefore, the allowed parameter space is constrained by the EWSB constraints. Interestingly, these constraints determine the size of s dL 13 9.63 × 10 −3 in our setup when the required SUSY threshold correction is obtained. In the middle M SUSY case, like M SUSY = 10 TeV, both constraints relate to the maximal value of s dL 13 , as shown in figure 3. It is notable that in our setup, the obtained mixing angles through the SUSY threshold corrections are related with each other due to the SO(10) unification. This means that some of the predictions of FV process can be calculated once one of the predictions of the process are given. Therefore, we may be able to obtain the specific predictions and can test the models by the experiments.
Finally, we comment on the SM fermion masses. Since we do not consider the SUSY threshold corrections for y t and y τ when we obtain the required SUSY threshold correction for y b and tan β, the predictions of y t and y τ are expected to slightly deviate from the experimental values. We found that the deviation of the top Yukawa coupling is O(1) % or less, while that of τ Yukawa coupling is about 3 ∼ 4 %. Note that the deviation of bottom Yukawa coupling should be small since we find the parameters for realizing the required SUSY threshold correction. The result for its deviation is less than 0.1 %, as expected.

Conclusion
In this paper, we considered the model with the Yukawa unification inspired by the SO (10) GUT. In this model, we investigated the maximal size of the mixing angles of diagonalizing matrices for the Yukawa matrices, originated from the off-diagonal parts of SUSY threshold corrections. We need large off-diagonal parts of A-terms to obtain the sizable SUSY threshold corrections. However, these parts are constrained by the SUSY FCNC processes. Moreover, the CCB constraints also relate to the maximal value of A-terms. Therefore, we can estimate the maximal contributions to the mixing angles from the SUSY threshold corrections by using these constraints.
We firstly found the required SUSY threshold correction and tan β for realizing the Yukawa unification. As reported in the other papers, the large tan β ∼ 50 and enough size of ∆y b are needed. We emphasize that required size of ∆y b is smaller than the maximal one when M SUSY ≤ 100 TeV. Therefore, we can obtain required ∆y b by choosing the input parameters at the GUT scale.
Next, we showed the results for the s ψL,R

13
because this element has a possibility to reach the size of corresponding element of the CKM matrix, s CKM

13
. As a result, we found that s qL,R 13 , and therefore, we have to find the appropriate value of ∆m 2 h to satisfy these constraints.

JHEP02(2020)103
Importantly, the mixing angles obtained from the SUSY threshold corrections are related with each other when one considers the SO(10) unification as in our setup. Therefore, if one of the FV processes is confirmed, the contributions to the other processes can be calculated, and we can test the models by the (future) experiments. In addition, if some experimental results of the FV process can be explained by these contributions, we may be able to obtain the signature of the SO(10) unification. Furthermore, if the CKM matrix is realized by the other contributions in GUT models, we may find the contribution of the unification according to our results, although it deviates by the RGE effects.

Acknowledgments
We thank B. Bajc for advise for starting this project and on present situation of Yukawa unification in SUSY SO(10) GUT. We also thank J. Hisano for fruitful discussion. Y.M.

A SUSY threshold corrections
To calculate the SUSY threshold corrections, chirality-flipping self-energies Σ ψ LR are important, which are calculated in ref. [45] in the decoupling limit M SUSY v. We can divide Σ ψ LR for ψ = u, d, e as For this calculation, we consider the effective Lagrangian within a decoupling limit M SU SY v, and hence, we neglect the chirality-flipping elements ∆ ψ LR and ∆ ψ RL . Then the diagonalizing matrices for sfermion mass squared matrices, denoted as Uψ, can be defined by where Γ = L, R is the chirality. We summarize the chirality-flipping self-energies by using these parameters.
The chirality-flipping self-energies from neutralino contribution are where M W is the W boson mass and g i is the gauge coupling with GUT notation.

JHEP02(2020)103
The chirality-flipping self-energies from chargino contribution are Note that we neglect the chargino contributions to up-quark self-energies (Σ uχ ± LR ) ij since these contributions are not tan β-enhanced ones.
At the SUSY scale, the above corrections are added to the fermion mass matrices as Then, the chirality-flipping self-energies Σ ψ LR and the SUSY threshold corrections ∆y ψ defined in eq. (3.1) have following relations:

B Constraints on the SUSY contributions
In this appendix, we summarize the constraints on the SUSY contributions from the SUSY FCNC and CCB constraints. First, we will show the SUSY FCNC constraints. For the squark sector, the upper bounds on SUSY contributions are mainly come from the constraints of the neutral meson mixing. This can be calculated by the general forms which can be found in ref. [71], using the MIPs. In our calculation, B d -B d mixing is important. In order to check the experimental bounds, it is useful to define the following variables: where H SM,NP eff is the effective Hamiltonian of the SM and the MSSM, respectively. In this notation, C B d and φ B d are related to the observables as = 0.738. The upper bounds on the SUSY contributions in slepton sector are mainly come from the constraints of l i → l j γ process. The branching ratio of this process can be calculated by chargino and neutralino contributions as and (A R ij ) can be estimated by the mass insertion approximation, and these expressions can be found in refs. [74,75]. Note that the current experimental upper bounds of BR(l i → l j γ) exp are summarized in ref. [53]: BR(µ → eγ) < 4.2 × 10 −13 , BR(τ → eγ) < 3.3 × 10 −8 , and BR(τ → µγ) < 4.4 × 10 −8 . We emphasize that the bound of τ → eγ is satisfied in our calculation, even when M SUSY = 2 TeV.
In addition to above constraints, there are CCB constraints, which can be read as [46][47][48][49][50][51] where k = max(i, j). Note that because of the sfermion mass spectrum, the bound for (Â ψ ) 33 is stronger than the other ones.

C Maximal SUSY threshold correction
In figure 1, we showed the required SUSY threshold correction to realize SO(10) Yukawa unification. Another important issue is whether the maximal SUSY threshold correction can be larger than the required one. The SUSY threshold corrections for the bottom Yukawa coupling mainly come from the gluino contribution: . Similar notations for m 2 t L,R and m 2 τ L,R are used. In this paper, the SUSY scale is defined as

JHEP02(2020)103
However, |a b | and µ have upper limits which come from the CCB and EWSB constraints. The maximal value of |a b | can be estimated as y , according to eq. (B.5). On the other hand, that of µ is complicated because µ depends on the MSSM Higgs mass parameters, m 2 hu and m 2 h d , through the tree-level condition to have the EWSB minimum in the Higgs potential [76] Note that the EWSB can be occurred when the following conditions are satisfied: Therefore, in the following discussion, we assume the case with m 2 hu < 0 to satisfy the first inequality of eq. (C.4).
From eq. (C.1), ∆y b is dominated by the contribution from µ because of the large tan β enhancement. Therefore, we first discuss the maximal value of µ which satisfies the CCB and EWSB constraints. It is clear that m 2 hu is important to obtain the maximal value of µ. Note that the effect of a b is discussed later.
The maximal value of |m 2 hu | is m 2 h d is needed to maximize µ, the EWSB constraints cannot be satisfied with |m 2 hu | = 2M 2 SUSY in the large M SUSY when m 2 h d > 0. This can be understood from eqs. (C.3) and (C.4). By combining these equation and inequalities, we can obtain following inequalities: Then, a large |m 2 hu | needs a small m 2 h d in order to satisfy the right part of the inequalities. This means that m 2 h d becomes negative when M SUSY M th where M th ≡ m Z 2 tan 2 β −1 ≈ 2.25 TeV when tan β ≈ 50. As a result, the values of |m 2 hu | and m 2 h d which maximize µ are and the maximal µ can be obtained from eq. (C.3) as In this case, the CCB constraint ofÂ u is |(Â u ) ij | 2 ≤ 0, which leads to |(Â u ) 33 | = |a t | = 0. However, that of |a b | is where is realized at the SUSY scale through the MSSM RGE. We found that δm 2 bt ∼ 0.08M 2 SUSY in our calculation. Therefore, the maximal value of |a b | is We emphasize that eq. (C.9) is obtained by assuming that the contribution from µ to ∆y b is dominated. In order to check this and to see the dependence of a b , we also consider the case where |a b | max is the following expression: Here, M 2 ≥ 0 is a parameter for taking |a b | max to be large. This modification can be justified by changing m 2 hu and m 2 h d as without conflicting with the CCB and EWSB constraints. Note that since m 2 hu < 0 is needed for the correct EWSB, the range for M 2 is 0 ≤ M 2 < 2M 2 SUSY . This modification affects the value of µ max through m 2 hu : µ max = |m 2 hu | = 2M 2 SUSY − M 2 . By using above expressions for m 2 hu , m 2 h d , and a b , we can calculate the maximal value of ∆y b in each M SUSY . We show the results in figure 4. The blue line shows the required SUSY threshold correction for realizing the SO(10) Yukawa unification, while the red ones JHEP02(2020)103 show the maximal SUSY threshold corrections with different values of M 2 , calculated from the above discussion. The thick red line which is the case with M 2 = 0 corresponds to the maximal SUSY threshold correction in the range of 0 ≤ M 2 < 2M 2 SUSY . Here, we assume that M 3 = M SUSY as in our setup, C 0 (|M 3  is not important to obtain the maximal value of ∆y b unless this is enough large to compete the other contributions. We found that δm 2 bt is less than 0.1M 2 SUSY in our calculation since we use the universal sfermion mass as the input at the GUT scale.
We can find important results from figure 4. First, the maximal SUSY threshold correction is larger than the required one when M SUSY ≤ 100 TeV, and then, we can find the parameter set which realizes the required SUSY threshold correction. Second, the maximal SUSY threshold correction decreases when M 2 becomes large. This means that ∆y b is dominated by the contribution from µ, as one can understand from the plot with M 2 = 2M 2 SUSY case, namely, µ max = 0 case. Therefore, we can obtain the required SUSY threshold correction from the appropriate values of m 2 hu and m 2 h d . Third, as the SUSY scale increases, the required SUSY threshold correction increases, while the maximal SUSY threshold correction decreases. This is because the rough expression for the maximal value of ∆y b is and this combination becomes small when the SUSY scale become high. Therefore, we can expect that there is an intersection point between the required and maximal SUSY threshold corrections. This means that there is an upper bound on the SUSY scale when we consider the SO(10) Yukawa unification, although its bound is expected to be far above M SUSY = 100 TeV. Finally, we comment on the situation in which we do not consider in the above discussion. If m 2 τ L + m 2 τ R < 2M 2 SUSY − m 2 Z 2 tan 2 β − 1 due to the small m 2 τ R through the MSSM RGE, the maximal value of |m 2 h d | is m 2 τ L + m 2 τ R from the CCB constraint of (Â e ) 33 when m 2 h d < 0. In this case, the maximal value of |m 2 hu | is modified as 13) in order to satisfy the EWSB constraints. This is smaller than 2M 2 SUSY , and therefore, the maximal value of µ is also smaller than the above case. Note that in this case, since , we can safely ignore the contribution from a b , according to the above results. Therefore, we can obtain the maximal ∆y b from figure 4 straightforwardly. For example, when |m 2 hu | 1.5M 2 SUSY , predicted value of ∆y b is around the red line with M 2 = 0.5M 2 SUSY . However, m 2 τ L + m 2 τ R is not far below 2M 2 SUSY because of the input parameter at the GUT scale. Therefore, we can conclude that the maximal SUSY threshold correction is always larger than the required one when M SUSY ≤ 100 TeV.

JHEP02(2020)103
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.