Contents lists available at ScienceDirect

Physics Letters B

www.elsevier.com/locate/physletb

Sparticle spectroscopy of the minimal S0(10) model

Takeshi Fukuyamaa, Nobuchika 0kadab, Hieu Minh Tranc *

a Research Center for Nuclear Physics (RCNP), Osaka University, Ibaraki, Osaka, 567-0047, Japan b Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA c Hanoi University of Science and Technology, 1 Dai Co Viet Road, Hanoi, Viet Nam

CrossMark

A R T I C L E I N F 0

Article history:

Received 2 December 2016

Received in revised form 5 February 2017

Accepted 6 February 2017

Available online 14 February 2017

Editor: J. Hisano

A B S T R A C T

The supersymmetric (SUSY) minimal S0(10) model is a well-motivated grand unified theory, where the Standard Model (SM) fermions have Yukawa couplings with only one 10-plet and one 126-plet Higgs fields and it is highly non-trivial if the realistic quark and lepton mass matrices can be reproduced in this context. It has been known that the best fit for all the SM fermion mass matrices is achieved by a vacuum expectation value of the 126-plet Higgs field being at the intermediate scale of around 0(1013) GeV. Under the presence of the S0(10) symmetry breaking at the intermediate scale, the successful SM gauge coupling unification is at risk and likely to be spoiled. Recently, it has been shown that the low-energy fermion mass matrices, except for the down-quark mass predicted to be too low, are very well-fitted without the intermediate scale. In order to resolve the too-low down quark mass while keeping the other fittings intact, we consider SUSY threshold corrections to reproduce the right down quark mass. It turns out that this requires flavor-dependent soft parameters. Motivated by this fact, we calculate particle mass spectra at low energies with flavor-dependent sfermion masses at the grand unification scale. We present a benchmark particle mass spectrum which satisfies a variety of phenomenological constraints, in particular, the observed SM-like Higgs boson mass of around 125 GeV and the relic abundance of the neutralino dark matter as well as the experimental result of the muon anomalous magnetic moment. In the resultant mass spectrum, sleptons in the first and second generations, bino and winos are all light, and this scenario can be tested at the LHC Run-2 in the near future.

© 2017 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license

(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

1. Introduction

Although the supersymmetric (SUSY) extension of the Standard Model (SM) is one of the most promising candidates for new physics beyond the SM, the results of the search for SUSY particles at the Large Hadron Collider (LHC) [1] has been increasing tensions to the minimal SUSY SM (MSSM). However, in most of MSSM analysis, the universal boundary conditions are taken into account for sfermion masses and gaugino masses at the grand unified theory (GUT) scale (MGUT — 1016 GeV), which seems to make the tensions stronger. It is worth considering a more general spectrum for sparticles [2].

Despite the tensions to the MSSM, the SUSY GUT is still a very attractive paradigm for physics beyond the SM at high energies, supported by the successful gauge coupling unification at MGUT — 1016 GeV in the context of the MSSM. Among various

* Corresponding author.

E-mail addresses: fukuyama@se.ritsumei.ac.jp (T. Fukuyama), okadan@ua.edu (N. Okada), hieu.tranminh@hust.edu.vn (H.M. Tran).

SUSY GUT models, the renormalizable minimal SUSY S0(10) model (minimal S0(10) model) has attracted a lot of attention, where one 10-plet and one 126-plet Higgs fields couple with 16-plet MSSM matters. The guiding principles for this GUT model are the minimality of the Higgs multiplets coupled with the MSSM matters and the renormalizablity. The minimal S0(10) model can not only reproduce the quark and lepton mass matrices at low energies but also provide the symmetry breaking pattern from S0(10) to the SM gauge group almost unambiguously [3-9]. However, there has been a very serious tension to the model, namely, in order to fit all the fermion mass matrix data at low energies the 126-plet Higgs field needs to develop a vacuum expectation value (VEV) at the intermediate scale of around 1013 GeV [9,10]. This indicates that many exotic states appear at the intermediate scale, so that the success of the SM gauge coupling unification at the GUT scale is spoiled.

Recently, one of the authors (T.F.), together with K. Ichikawa and Y. Mimura, has revisited the minimal S0(10) model to fit the fermion mass matrix data with a new strategy [11]. Here, the authors of this paper have noticed that thanks to the recent progress

http://dx.doi.org/10.1016/j.physletb.2017.02.021

0370-2693/© 2017 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.

in the neutrino oscillation experiments, the mass matrix data in the lepton sector are more accurate than those in the quark sector. Instead of the quark mass matrix data, they have employed the mass matrix data in the lepton sector as the inputs to fit the fermion mass matrices. Through the x2-analysis for the general parameter space, they have found parameter regions to nicely fit the fermion mass matrix data. Although the best fit is obtained with the intermediate scale VEV of the 126-plet Higgs field, which is consistent with the result in [9,10], it has been found that a parameter set without the intermediate scale can results in a very good fit for the fermion mass matrix data, except for the down quark mass being too small.

The fitting result with no intermediate scale is particularly interesting since in this case the successful gauge coupling unification is kept intact. As discussed in [11], the too small down quark mass can be resolved by SUSY threshold corrections through ¿4-terms [14]. However, not to alter the other fittings which are very good, the ¿4-term corrections must be sizable only for the down quark Yukawa coupling, in other words, the soft SUSY breaking parameters are flavor-dependent. Motivated by this fact, we investigate a flavor-dependent sparticle mass spectrum in the context of the minimal S0(10) model and a variety of phenomeno-logical constraints. Considering very severe SUSY FCNC constraints, we impose the universal mass for the sfermions in the matter 16-plets of the first and second generations. Among other phe-nomenological constraints, we focus on a parameter set which can not only reproduce the measured SM-like Higgs boson of around 125 GeV and the observed relic density of the neutralino dark matter, but also fill the discrepancy between the experimental result and the SM prediction for the muon anomalous magnetic dipole moment a/x = 1 (g— 2). In order for a/x to receive sizable spar-ticle contributions, the universal mass inputs for sfermions in the first and second generation 16-plets as well as the gauginos must be < 1 TeV, while the 125 GeV Higgs mass requires the soft mass input for the 3rd generation 16-plet to be much larger. Hence, we have a large mass splitting between sfermions in the first two generations and the third generation. A similar direction was previously addressed in [15], though our arguments have the unique reasoning for them.

2. Minimal SO(10) model and fermion mass matrices

In this section, we briefly review the minimal S0(10) model and its fitting for the fermion mass matrices with the new strategy introduced in [11], focusing on the points which are used in this paper.

One of the most important features of the minimal S0(10) model is that all SM matter particles are set in a single 16-plet (f; i = 1, 2, 3) for each generation. Renormalizability indicates that its coupling with the Higgs multiplets is limited in the Yukawa coupling: 16 ® 16 = 10 e 120 e 126. Therefore, the Higgs multiplet must belong to 10, 120 and 126 to make the Yukawa couplings invariant under S0(10). One 10-plet Higgs (H) is indispensable but not sufficient since we have Cabbibo-Kobayashi-Maskwa (CKM) and Maki-Nakagawa-Sakata (MNS) mixing matrices. Minimally, we add only one 126-plet Higgs (A). This 126-plet Higgs (not 120-plet Higgs) is required from the following reasons. Under the Pati-Salam subgroup of S0(10), the 126-plet Higgs includes both SU(2)R and SU(2)l triplets ((10, 1, 3) and (10, 3, 1)) which generate the neutrino mass matrix through the type I [16] and type II [17] seesaw mechanism, respectively. Also the VEVs of the 126 + 126-plet Higgs fields in the B — L = ±2 direction reduce the rank of the S0(10) gauge symmetry in the renormalizable model and guarantees the conservation of R-parity. Thus, in the minimal S0(10) model the Yukawa coupling has the form of

Wy = 2 hyf fjH + 1 j fj A.

Due to the S0(10) algebra, the coupling matrices are symmetric, hjj = hji and fj = fji.

After the S0(10) symmetry is broken down to the SM gauge symmetry, the fermion Yukawa matrices are given as follows:

Yu = h + T2 f, Yd = n(h + f), Ye = n(h — 3 f),

Y v = h — 3r2 f, (2)

where ri and r2 depend on the Higgs mixing (doublet Higgs mixing in 10 and 126), and h and f are the original Yukawa matrices h and f multiplied by Higgs mixings,

U12 U11

h = V 11h, f = —^ f, r1 = —. 11 ' V3r1 ' 1 V11

r2 = r1—, (3) U12

where U and V are the unitary matrices to diagonalize the mass matrices of the MSSM Higgs doublets in the 10-plet and 126-plet. See [18] for details.

The charged fermion mass matrix and the Dirac neutrino mass matrix are obtained as

Mu = YuVu, Md = YdVd, Me = YeVd, M^ = Yv Vu,

where Vu and Vd are the VEVs of up- and down-type Higgs fields. From the mass formula Eq. (2), we obtain the relation of the mass matrices as

Md = Me + --F = rMu + F,

1 — r2

rMD = Me + 3 F, where

r = r1-= r1 cot P,

and the matrix F is F = r(1 — r2) fVu .

Roughly, we obtain r ~ m^/mt.

The right-handed Majorana neutrino mass matrix is obtained

MR = V2fvR,

where vR is a VEV of 126. Practically, we denote Mr = CrVrf,

cr = (11)

U12 V13

for the current notation. Because U12 and V13 are components of the diagonalization unitary matrix for doublet Higgs fields, cR has a minimal value. The size of cR is related to the size of original f coupling, which will be important to derive the GUT scale threshold correction for flavor violation. The seesaw neutrino mass matrix can be written as [16,17]

Mv = Ml — MDM R 1(MD )t

where ML is the left-handed neutrino Majorana mass which comes from SU(2)l triplet coupling,

Table 1

The reference parameters used in the fit.

Fixed Values

mu 3.961 x 10-4 GeV [19]

mt 71.0883 GeV [19]

me 3.585 x 10-4 GeV [19]

mß 7.5639 x 10-2 GeV [19]

mz 1.3146 GeV [19]

A<! 7.54 x 10-5 eV2 [21]

Am2trn 2.4 x 10-3 eV2 [21]

012 0.583 [21]

»23 0.710, n/4 [21]

013 0.156 [21]

Ml = cLvLf. (13)

We note here that the S0(10) breaking vacuum can be specified by one complex parameter for the minimal Higgs contents: 10 +126 +126 + 210 [6,7]. The 210-plet is required to connect the 10-plet and the 126-plet Higgs fields. The coefficient r2 is determined by the complex parameter. This is one of the reason why the model with the minimal number of Higgs multiplets in the Higgs superpotential does not provide a solution consistent with the gauge coupling unification in the minimal SUSY version of the minimal S0(10) [10].

In Ref. [11], using the parameters in lepton sector as inputs, the X2 analysis has been performed for fitting the quark masses and mixings under the current observables. The fit provides information about the behavior of the fitting parameters with respect to the scale of 126 + 126 VEVs (vR). The scale plays a key role for the breaking the rank-5 S0(10) symmetry down to the rank-4 SM gauge symmetry and is also important to see how the light neutrino mass is generated by the seesaw mechanism [16,17].

The values and the uncertainties used in the fit at the GUT scale are summarized in Table 1, where the quark mass values are taken from Table III in Ref. [19] (MSSM, tan ß = 10), the quark mixing angles and the CP-phase at the GUT scale are evaluated from the PDG data [20] by DR scheme, and the neutrino oscillation data are taken from the global fit results in Ref. [21].

In Ref. [11], the best fit is obtained when vR ~ 1013 GeV as was reported in Refs. [9,10]. However, such a low VEV scale spoils the successful gauge coupling unification in the context of the MSSM [9,10]. Therefore, we focus on solutions for vR being around the GUT scale. In Table 2, we cite the results in Ref. [11] for the type II seesaw dominated case for the neutrino mass matrix. Although they are not the best fit found in Ref. [11], the fits are very good except for the down quark mass.

In the results, the main pull comes from the down quark mass which is too low. However, this small down quark mass is advantageous for making the proton lifetime longer under the relation of (Mv)n = (M v) 12 = 0 [18,22]. In order to resolve the too-low down quark mass, we consider SUSY threshold corrections from the gluino loop diagrams [14] as

&md =---2 (mdß tan ß + ^^ , (14)

m2 V V2 J

where mg and mq are the gluino and squark masses, respectively, Vd is the VEV of the down-type Higgs doublet, and ad is the A-term parameter defined as

Aoft 3 addt QHd + h.c. (15)

in the trilinear scalar coupling among the squarks and the down-type Higgs doublet. For example, when we take mg ~ m^ ~ 1 TeV, ad ~ 10 GeV is sufficient to yield the right size of the down quark

Parameters to fit

mc 0.1930 ± 0.025 GeV [19]

md (9.316 ± 3.8) x 10-4 GeV [19]

ms (1.76702 ± 0.5) x 10-2 GeV [19]

mb (0.9898 ± 0.03) GeV [19]

Vus 0.224 ± 0.002 [20]

Vcb (3.7 ± 0.13) x 10-2 [20]

Vub (3.7 ± 0.45) x 10-3 [20]

^KM 1.18 ± 0.2 [20]

Table 2

The fit result for type II (which is called Type I+II in many literatures) cited from Ref. [11].

Type II

CrVr (GeV) 8.86 x 1016 9.22 x 1016

023 0.710 n/4

ae -0.66648 -0.25301

aß -2.8148 2.8177

ax -0.53961 -0.84287

«2 -2.8709 -3.1146

«3 -1.9809 -2.8604

^PMNS -2.3550 -3.1131

log10(mt/GeV) -11.207 -11.173

15.545 16.156

Arg(S) 0.43912 0.51567

mc (GeV) 0.1978 0.1989

md (GeV) 0.0004138 0.0003936

ms (GeV) 0.01980 0.02028

mb (GeV) 0.9903 0.9901

Vuc 0.2240 0.2241

Vsb 0.003765 0.003724

Vub 0.03694 0.03695

SKM 1.195 1.160

mc 0.191 0.236

md -1.363 -1.416

ms 0.426 0.522

mb 0.017 0.010

Vuc -0.000 0.027

Vsb 0.144 0.054

Vub -0.044 -0.039

Skm 0.075 -0.099

r 0.0230 0.0233

r2 2.15 + 0.227i 2.22 - 0.16

X2 2.10 2.35

mass from the SUSY threshold corrections. In analysis of sparticle mass spectrum, a universal A-term parameter ( A0) is often used, which has a relation to ad as YdA0 = ad with Yd being a down-type quark Yukawa coupling. We can see that if we apply the universal A-term parameter for the SUSY threshold corrections, the strange and bottom quarks also receive sizable corrections which ruin the good mass fitting. Hence, non-universal A-term parameters are necessary.

3. Sparticle mass spectroscopy and current experiments

In the minimal S0(10) model, the quark and lepton supermultiplets in a same generation are unified into a 16-representation. Therefore, their corresponding soft masses and trilinear couplings are also unified in the same way, significantly reducing the total number of free parameters compared to the general MSSM case. For simplicity, we assume that at the GUT scale the soft masses

of the first and the second generations are universal while that of the third generation is independent. As a result, in our analysis, the boundary conditions for the soft parameters at the GUT scale consist of five free parameters {M1/2,m0, m3, A0, tan/)} which are, respectively, the universal gaugino masses, the universal masses among the first two generations sfermions and the Higgs doublets, the third generation sfermion mass, the universal trilinear coupling, and the ratio between the VEVs of two MSSM Higgs doublets.1

For numerical analysis, we employ SOFTSUSY version 3.7.3 [23] to solve the RG equations with inputs at the GUT scale, and calculate the corresponding sparticle mass spectrum at the low energies. The constrained observables are computed by micrOMEGAs version 4.3.1 [24-26]. The Higgs sector is tested with the package HiggsBounds version 4.3.1 [27-30] by considering the exclusion limits from Higgs searches at the LEP, the Tevatron and the LHC experiments. We also check the validity of the benchmark with the package SModelS version 1.0.4 [31] which decomposes the model signal into simplified model spectrum topologies and compares with the LHC bounds.

We consider a variety of phenomenological constraints from collider physics. While the LEP limits is important to constrain the sleptons, neutralino and chargino sectors, the LHC experiments are more sensitive to colored sparticles. The SM-like Higgs boson mass obtained by the combined ATLAS and CMS analysis [32], the sparticle contributions (Ap) to the electroweak p -parameter [20], and the anomalous muon magnetic moment [33-36] are taken into account. Other constraints are from the branching ratios of rare de-

cay processes: b -D± — t±vT [20],

>s + y [37], B0 — D± — |±v„ [20],

and a constraint on the Kaon

[38], B± — t±Vt [37],

system characterized by the decay parameter [39]

Rl23 =

Vus(Kl2) Vud(0+

Vus( Kb)

V ud(n 12)

where the CKM matrix elements, Vus and Vud, are measured from the corresponding 3-body semileptonic Kaon decay (Ki3), 2-body leptonic Kaon and pion decay (Ki2, n2), and super-allowed nuclear beta decay (0+ ^ 0+). The limits for these constraints are as follows

mh = 125.09 ± 0.21(stat.) ± 0.11(syst.), (17)

Ap < 8.8 x 10-4, (18)

7.73 x 10-10 < Aa^ < 42.14 x 10-10, (19)

3.11 x 10-4 < BR(b — s + y)< 3.87 x 10-4, [2a ] (20)

1.7 x 10-9 < BR(B0 — ) < 4.3 x 10-9, [2a ] (21)

BR p(b± — t±vt) 0.70 < SM( x-< 1.82, BR (b ± — t ±vt) [2a ] (22)

5.07 x 10-2 < BR(D± — t±vt) < 6.03 x 10-2, [2a ] (23)

5.06 x 10-3 < BR(D± — |±v|) < 6.06 x 10-3, Rl23 = 1.004 ± 0.007. [2a ] (24) (25)

Assuming the R-parity conservation, the model predicts the lightest neutralino to be a dark matter (DM) candidate that can

1 As we have discussed, non-universal trilinear couplings are essential to resolve the too-low down quark mass while keeping the other fermion masses intact. Nevertheless, we can employ the universal trilinear coupling in our analysis for the

renormalization group (RG) evolutions of the soft parameters, since the trilinear couplings for the first two generations contribute to the RG evolutions only through their Yukawa couplings and their effects are negligibly small.

Table 3

Benchmark particle mass spectrum for a set of inputs (in units of GeV). M1/2, m0, m3, A0 are, respectively, the universal gaugino mass, the universal sfermion soft mass of the 1st and 2nd generations and Higgs doublets, the soft mass of the 3rd generation sfermions, and the common trilinear coupling at the GUT scale.

M1/2 900

m0 393.515

m3 9000

A0 -13920

tan p 10

h 125.09

H0, A0 8682

H ± 8683

g 2070

X 0 399.7, 768.4

X 0 8739, 8740

*Í2 768.6, 8740

u, c L 1558

u, c R 1524

d, ~sl 1560

d, sr 1514

?1,2 5832,7649

b1,2 7644, 9080

~e,i vl 585.3, 583.8

e,iL 590.6, 589.1

e,iR 404.0, 399.8

VI 8976

T 1,2 8903, 8978

ap 3.30 x 10-7

aai 1.39 x 10-9

BR(b — sy) 3.33 x 10-4

BR(B0 — |+ i ) 3.08 x 10-9

BRMSSM(B± — T BRi"(B± — T± ±Vt) V t) 1.00

BR(D± — T± Vt) 5.17 x 10-2

BR(D± — |±V|) 5.33 x 10-3

Rl23 1.000

^CDMh2 0.1188

aX -P (pb) 1.507 x 10-

asXD-p (pb) 4.589 x 10-i

explain a wide range of cosmological observations such as gravitational lensing and galaxy rotation curves. Therefore, we furthermore take into account constraints from DM searches which play an essential role in the SUSY phenomenology and DM physics. The DM relic density is precisely measured by the Planck experiment [40]:

^CDMh2 = 0.1188 ± 0.0010. (26)

Regarding to DM direct detections, the most severe constraint on spin-independent cross section between DM and nuclei is set by the LUX experiment [41]:

a* p < 8 x 10

pb, mwiMP ~ 400 GeV. (27)

On the other hand, the DM indirect detection at the IceCube experiment [42] with the vV channel imposes the most stringent upper limits on the spin-dependent cross section between DM and nuclei for a DM mass of about a few hundreds GeV:

ctsD"P ~ 10-5 pb, mwiMP ~ 400 GeV. (28)

in Table 3, a benchmark mass spectrum is shown for a given set of inputs at the GUT scale with tan /) = 10 corresponding to

2 x 10"8 1 x 10"' 5 x 1(T'

< 2 x 10-i 1x10 5x10

.-10 .

2 x 10"10 L

600 m-[ GeV]

Fig. 1. The rough estimate of Aa/ as a function of smuon mass, along with the experimental result, 7.73 x 10-10 < Aa/ < 42.14 x 10-10 (dashed horizontal lines).

the x2 fitting results in Ref. [11]. The benchmark satisfies all the above phenomenological constraints. Here, the heavy squark masses and large trilinear couplings of the third generation at the order of 0(10) TeV play a crucial role to reproduce the Higgs boson mass of around 125 GeV via quantum corrections.2 Since the third generation sfermions are all heavy, their visible effect is only the contribution to the Higgs boson mass. We find that the SUSY contribution to the muon anomalous magnetic moment is at the order of 0(10-9) which satisfies the experimental limits at the level of 2a. It is due to a large mass splitting between the first two generations and the third generation sfermions and relatively light smuons and gauginos. In order to illustrate a response of Aa/ to the sparticle mass, we show Aa/ as a function of smuon mass in Fig. 1. Here, we have used a very rough formula for the SUSY contributions as [43]

a2 ( m

4n \ m

Aa„ — tan p

where a2 is the SM SU(2) gauge coupling.

For this benchmark, the lightest sparticle (LSP) is a bino-like neutralino playing the role of the DM, and the next-to-the LSP (NLSP) is the right-handed smuon. A mass degeneracy between these LSP and NLSP is required to yield the correct DM relic density via neutralino-slepton co-annihilation processes. Since the scattering cross sections between the neutralino LSP and nuclei are a few orders of magnitude smaller than the current sensitivity of the LUX experiment, drastic improvements on new generations of DM direct and indirect searches are necessary to capture signatures of the bino-like neutralino. In spite of that, the light sleptons, neutralinos and charginos can be explored at the LHC Run-2 in the near future. The projected high energy linear collider with the center of mass energy around 1 TeV can shed more lights on the physics of these relatively low mass sparticles due to its high precision.

In the minimal S0(10) model, after the parameter fitting, the neutrino Dirac Yukawa coupling (Yv) is unambiguously deter-

2 The large A0 in Table 3 may generate a true minimum away from the SM vacuum and make the SM vacuum unstable. We have checked that the absolute stability conditions for the SM vacuum are satisfied in the direction of the 3rd generation sfermions, namely, A2 < 3(m? + m? + m2Hu + |/|2) [12]. On the other hand, the stability conditions in the directions of the 1st and 2nd generation sfermions are not satisfied, since they are relatively light compared to the corresponding A-term parameters. However, such charge/color breaking vacua appear far away from the SM vacuum because of the very small Yukawa couplings of the 1st and 2nd generations and the tunneling rate is extremely small [13]. Hence, we conclude that our benchmark spectrum is cosmologically safe.

m-[ GeV]

Fig. 2. The branching ratio BR(/ ^ ey) as a function of the left-handed slepton mass. The solid lines from left to right correspond to a variety of choices of the GUT scale, MGUT = 1.20, 1.23, 1.30, and 2.0 in units of 1016 GeV, respectively. The current experimental upper bound by the MEG experiment [45] is denoted as the horizontal dashed line.

mined. Using this neutrino Dirac Yukawa couplings, we calculate the lepton flavor violating (LFV) processes. The LFV effect most directly emerges in the left-handed slepton mass matrix generated through the RG equations such as [44]

d/ K) i,=№ d/ K)

= / — \m%

'ij a/

I MSSM

+ jn (m?Y?Yv + YVYvm? + 2YVm? Y

+ 2mHuY V Y ?)

where the first term in the right hand side denotes the normal MSSM term with no LFV. We employ the data obtained in Ref. [11] for the heavy Majorana neutrino mass eigenvalues,

MR1 = 6.9 x 108 GeV, MR2 = 3.3 x 1014 GeV

MR3 = 1.2 x 1016 GeV, and the neutrino Dirac Yukawa coupling,

/0.000111 0.000203 + 0.000217 i 0.00888 + 0.00372 i\ Yv = 0.000440 0.0308 - 0.0248 i 0.0426 + 0.0013 i , \ 0.00607 -0.0276 - 0.0069 i 0.990 - 0.278 i J

for the case of cRvR = 8.86 x 1016 GeV (the second column in Table 2). We solve Eq. (30) from the GUT scale to low energies with the inputs of the Majorana neutrino masses, the Yv data, the slepton mass matrix m£ = m? = diag(m0,m0, m3) and mHu = m0. Using the generated off-diagonal elements of the (left-handed) slepton mass matrix squared m2 at low energies, we roughly estimate the LFV decay rate of the charged leptons by [44]

a2 |(Am2)ij|2

r(£i ^ £jy) m|( x

tan2 p.

Our results for the branching ratio of the / ^ ey process, which provides the most severe constraint, are shown in Fig. 2 for a various values of the GUT scale, along with the final results of the MEG experiment [45], BR(/+ ^ e+y) < 4.2 x 10-13 at 90% C.L. (horizontal dashed line). We find that the resultant branching ratio is sensitive to the choice of MGUT, and hence a definite prediction

requires a precise determination of the GUT scale. We leave this issue for future work.

Finally we consider proton decay in our scenario. in SUSY GUT models, the proton decay process induced by the dimension-5 operators via color-triplet Higgsino exchanges [46] is severely constrained. in evaluating the proton decay amplitude through the dimension-five operators, the information of the fermion Yukawa matrices is necessary. Again, we employ the data of the Yukawa matrices obtained in Ref. [11]. Recently, in Ref. [22] the proton decay widths have been calculated in detail with the same Yukawa matrix data and it has been found that the proton lifetime is about 6 times larger than the current experimental limit. Here, a special Yukawa structure, namely f11 and /12 are extremely small, plays a crucial role to suppress the proton decay amplitude. Since the same Yukawa matrix data are used, we simply follow Ref. [22] for our calculation of the proton decay amplitudes. The only difference from the calculation in Ref. [22] is the choice of the particle mass spectrum in our benchmark on Table 3, while mq = ^ = 2 TeV has been taken in Ref. [22]. With our benchmark particle mass spectrum, we have estimated sparticle loop corrections using the formulas given, for example, in Ref. [47] and found that the proton decay width is enhanced roughly by a factor of 1.5 from the one obtained in Ref. [22]. Therefore, the proton lifetime is about 4 times larger than the current experimental bound. This proton lifetime can be tested by the projected Hyper-Kamiokande experiment [48].

4. Conclusions

The SUSY minimal S0(10) model is a well motivated GUT scenario for new physics at very high energies, where the SM fermions have Yukawa couplings with only one 10-plet and one 126-plet Higgs fields. With such a limited number of free parameters, it is highly non-trivial if this model can fit all the data of the SM fermion mass matrices at low energies. it has been known that the best fit for all the SM fermion mass matrices is achieved by a 126-plet Higgs field VEV being at the intermediate scale of around 0(1013) GeV. in the minimal S0(10) model with a renor-malizable Higgs sector, many exotic states emerge in association with the intermediate-scale S0(10) breaking, so that the success of the SM gauge coupling unification is spoiled.

The fitting of the fermion mass matrices in the minimal S0(10) model has been revisited in Ref. [11] and it has been found that the low-energy fermion mass matrices, except for the down-quark mass predicted to be too low, are very well-fitted without the intermediate scale. This result is particularly interesting since the successful gauge coupling unification at the GUT scale is kept intact. We have discussed that the too-small down quark mass can be resolved by the SUSY threshold corrections via the gluino loop diagrams with a suitable choice of soft parameters. We have found that the flavor-dependent soft parameters are essential in order for the SUSY threshold corrections to work only for resolving the too-small down quark mass.

Motivated by the need of the flavor-dependent soft parameters, we have considered flavor-dependent inputs for the soft parameters at the GUT scale and calculated particle mass spectra at low energies. For particle mass spectra, we have taken a variety of phe-nomenological constraints into account. The most important constraints are (i) the SM-like Higgs boson mass of around 125 GeV, (ii) the observed DM relic density for the neutralino LSP and (iii) the discrepancy between the measured value and the SM prediction for the muon anomalous magnetic moment, which is filled by sparticle contributions. The resultant sparticle mass spectra are found to be very characteristic such that sleptons in the first and second generations, bino and winos are all light and have masses

< 1 TeV. We present a benchmark particle mass spectrum which satisfies all phenomenological constraints. Our scenario with the light sparticles can be tested at the LHC Run-2 in the near future. With the benchmark particle mass spectrum, we have estimated the proton lifetime and found it to be about 4 times larger than the current experimental limit, which can be tested in the the projected Hyper-Kamiokande experiment.

Acknowledgements

We would like to thank Yukihiro Mimura for useful discussions and comments. H.M.T. would like to thank the Department of Physics and Astronomy at the University of Alabama for hospitality during his visit. The work of T.F is supported in part by Grant-in-Aid for Science Research from Japan Ministry of Education, Science and Culture (No. 26247036). The work of N.O. is supported in part by the United States Department of Energy Grant (DE-SC0013680). The work of H.M.T. is supported in part by Vietnam Education Foundation, and Vietnam National Foundation for Science and Technology Development (NAFOSTED) under the grant No. 103.01-2014.22.

References

[1] https://twiki.cern.ch/twiki/bin/view/AtlasPublic/SupersymmetryPublicResults; https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsSUS.

[2] See, for example C.F. Berger, J.S. Gainer, J.L. Hewett, T.G. Rizzo, Supersymme-try without prejudice, J. High Energy Phys. 0902 (2009) 023, arXiv:0812.0980 [hep-ph];

J.A. Conley, J.S. Gainer, J.L. Hewett, M.P. Le, T.G. Rizzo, Supersymmetry without prejudice at the LHC, Eur. Phys. J. C 71 (2011) 1697, arXiv:1009.2539 [hep-ph]; J.A. Conley, J.S. Gainer, J.L. Hewett, M.P. Le, T.G. Rizzo, Supersymmetry without prejudice at the 7 TeV LHC, arXiv:1103.1697 [hep-ph].

[3] K.S. Babu, R.N. Mohapatra, Predictive neutrino spectrum in minimal SO(10) grand unification, Phys. Rev. Lett. 70 (1993) 2845, arXiv:hep-ph/9209215.

[4] K. Matsuda, Y. Koide, T. Fukuyama, Can the SO(10) model with two Higgs doublets reproduce the observed fermion masses?, Phys. Rev. D 64 (2001) 053015, arXiv:hep-ph/0010026;

K. Matsuda, Y. Koide, T. Fukuyama, H. Nishiura, How far can the SO(10) two Higgs model describe the observed neutrino masses and mixings?, Phys. Rev. D 65 (2002) 033008, Erratum: Phys. Rev. D 65 (2002) 079904, arXiv:hep-ph/0108202.

[5] T. Fukuyama, N. Okada, Neutrino oscillation data versus minimal supersymmet-ric SO(10) model, J. High Energy Phys. 0211 (2002) 011, arXiv:hep-ph/0205066.

[6] T.E. Clark, T.K. Kuo, N. Nakagawa, A SO(10) supersymmetric grand unified theory, Phys. Lett. B 115 (1982) 26;

C.S. Aulakh, R.N. Mohapatra, Implications of supersymmetric SO(10) grand unification, Phys. Rev. D 28 (1983) 217;

D.G. Lee, Symmetry breaking and mass spectra in the minimal supersymmetric SO(10) grand unified theory, Phys. Rev. D 49 (1994) 1417;

C.S. Aulakh, B. Bajc, A. Melfo, G. Senjanovic, F. Vissani, The minimal super-symmetric grand unified theory, Phys. Lett. B 588 (2004) 196, arXiv:hep-ph/0306242.

[7] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, N. Okada, General formulation for proton decay rate in minimal supersymmetric SO(10) GUT, Eur. Phys. J. C 42 (2005) 191, arXiv:hep-ph/0401213;

T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, N. Okada, SO(10) group theory for the unified model building, J. Math. Phys. 46 (2005) 033505, arXiv:hep-ph/0405300;

T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, N. Okada, Higgs masses in the minimal SUSY SO(10) GUT, Phys. Rev. D 72 (2005) 051701, arXiv:hep-ph/0412348;

B. Bajc, A. Melfo, G. Senjanovic, F. Vissani, The minimal supersymmetric grand unified theory. 1. Symmetry breaking and the particle spectrum, Phys. Rev. D 70 (2004) 035007, arXiv:hep-ph/0402122;

B. Bajc, A. Melfo, G. Senjanovic, F. Vissani, Fermion mass relations in a super-symmetric SO(10) theory, Phys. Lett. B 634 (2006) 272, arXiv:hep-ph/0511352;

C.S. Aulakh, A. Girdhar, SO(10) MSGUT: spectra, couplings and threshold effects, Nucl. Phys. B 711 (2005) 275, arXiv:hep-ph/0405074.

[8] B. Bajc, G. Senjanovic, F. Vissani, b-tau unification and large atmospheric mixing: a case for noncanonical seesaw, Phys. Rev. Lett. 90 (2003) 051802, arXiv:hep-ph/0210207;

B. Bajc, G. Senjanovic, F. Vissani, Probing the nature of the seesaw in renormal-izable SO(10), Phys. Rev. D 70 (2004) 093002, arXiv:hep-ph/0402140;

H.S. Goh, R.N. Mohapatra, S.-P. Ng, Minimal SUSY S0(10), b tau unification and large neutrino mixings, Phys. Lett. B 570 (2003) 215, arXiv:hep-ph/0303055;

H.S. Goh, R.N. Mohapatra, S.-P. Ng, Minimal SUSY S0(10) model and predictions for neutrino mixings and leptonic CP violation, Phys. Rev. D 68 (2003) 115008, arXiv:hep-ph/0308197;

B. Dutta, Y. Mimura, R.N. Mohapatra, Phys. Rev. D 69 (2004) 115014, arXiv:hep-ph/0402113;

B. Dutta, Y. Mimura, R.N. Mohapatra, Neutrino masses and mixings in a predictive S0(10) model with CKM CP violation, Phys. Lett. B 603 (2004) 35, arXiv:hep-ph/0406262. [9] K.S. Babu, C. Macesanu, Neutrino masses and mixings in a minimal S0(10) model, Phys. Rev. D 72 (2005) 115003, arXiv:hep-ph/0505200.

[10] S. Bertolini, T. Schwetz, M. Malinsky, Fermion masses and mixings in S0(10) models and the neutrino challenge to SUSY GUTs, Phys. Rev. D 73 (2006) 115012, arXiv:hep-ph/0605006.

[11] T. Fukuyama, K. Ichikawa, Y. Mimura, Revisiting fermion mass and mixing fits in the minimal SUSY S0(10) GUT, Phys. Rev. D 94 (7) (2016) 075018, arXiv:1508.07078 [hep-ph].

[12] J.F. Gunion, H.E. Haber, M. Sher, Charge/color breaking minima and Aparameter bounds in supersymmetric models, Nucl. Phys. B 306 (1988) 1;

J.E. Camargo-Molina, B. 0'Leary, W. Porod, F. Staub, Stability of the CMSSM against sfermion VEVs, J. High Energy Phys. 1312 (2013) 103, arXiv:1309.7212 [hep-ph].

[13] A. Kusenko, P. Langacker, G. Segre, Phase transitions and vacuum tunneling into charge and color breaking minima in the MSSM, Phys. Rev. D 54 (1996) 5824, arXiv:hep-ph/9602414.

[14] L.J. Hall, R. Rattazzi, U. Sarid, The top quark mass in supersymmetric S0(10) unification, Phys. Rev. D 50 (1994) 7048, arXiv:hep-ph/9306309.

[15] M. Ibe, T.T. Yanagida, N. Yokozaki, Muon g-2 and 125 GeV Higgs in split-family supersymmetry, J. High Energy Phys. 1308 (2013) 067, arXiv:1303.6995 [hep-ph];

S. Akula, P. Nath, Gluino-driven radiative breaking, Higgs boson mass, muon g-2, and the Higgs diphoton decay in supergravity unification, Phys. Rev. D 87 (11) (2013) 115022, arXiv:1304.5526 [hep-ph];

M.A. Ajaib, I. Gogoladze, Q. Shafi, C.S. Un, Split sfermion families, Yukawa unification and muon g — 2, J. High Energy Phys. 1405 (2014) 079, arXiv:1402.4918 [hep-ph];

K.S. Babu, I. Gogoladze, Q. Shafi, C.S. Un, Muon g-2, 125 GeV Higgs boson, and neutralino dark matter in a flavor symmetry-based MSSM, Phys. Rev. D 90 (11) (2014) 116002, arXiv:1406.6965 [hep-ph];

M. Adeel Ajaib, I. Gogoladze, Q. Shafi, GUT-inspired supersymmetric model for h — YY and the muon g-2, Phys. Rev. D 91 (9) (2015) 095005, arXiv:1501.04125 [hep-ph];

D. Chowdhury, N. Yokozaki, Muon g — 2 in anomaly mediated SUSY breaking, J. High Energy Phys. 1508 (2015) 111, arXiv:1505.05153 [hep-ph];

I. Gogoladze, Q. Shafi, C.S. Un, Reconciling the muon g — 2, a 125 GeV Higgs boson, and dark matter in gauge mediation models, Phys. Rev. D 92 (11) (2015) 115014, arXiv:1509.07906 [hep-ph];

N. 0kada, H.M. Tran, 125 GeV Higgs boson mass and muon g — 2 in 5D MSSM, Phys. Rev. D 94 (7) (2016) 075016, arXiv:1606.05329 [hep-ph]; W. Yin, N. Yokozaki, Splitting mass spectra and muon g — 2 in Higgs-anomaly mediation, Phys. Lett. B 762 (2016) 72, arXiv:1607.05705 [hep-ph].

[16] P. Minkowski, / -— eY at a rate of one out of 109 muon decays?, Phys. Lett. B 67 (1977) 421;

T. Yanagida, Horizontal symmetry and masses of neutrinos, Conf. Proc. C 7902131 (1979) 95;

M. Gell-Mann, P. Ramond, R. Slansky, Complex spinors and unified theories, Conf. Proc. C 790927 (1979) 315, arXiv:1306.4669 [hep-th]; S.L. Glashow, The future of elementary particle physics, NAT0 Sci. Ser. B 61 (1980) 687;

R.N. Mohapatra, G. Senjanovic, Neutrino masses and mixings in gauge models with spontaneous parity violation, Phys. Rev. D 23 (1981) 165.

[17] G. Lazarides, Q. Shafi, C. Wetterich, Proton lifetime and fermion masses in an S0(10) model, Nucl. Phys. B 181 (1981) 287;

R.N. Mohapatra, G. Senjanovic, Neutrino masses and mixings in gauge models with spontaneous parity violation, Phys. Rev. D 23 (1981) 165; M. Magg, C. Wetterich, Neutrino mass problem and gauge hierarchy, Phys. Lett. B 94 (1980) 61;

J. Schechter, J.W.F. Valle, Neutrino masses in SU(2) x U(1) theories, Phys. Rev. D 22 (1980) 2227.

[18] B. Dutta, Y. Mimura, R.N. Mohapatra, Neutrino mixing predictions of a minimal S0(10) model with suppressed proton decay, Phys. Rev. D 72 (2005) 075009, arXiv:hep-ph/0507319.

[19] K. Bora, Updated values of running quark and lepton masses at GUT scale in SM, 2HDM and MSSM, Horizon, J. Phys. 2 (2012) 112, arXiv:1206.5909 [hep-ph].

[20] K.A. Olive, et al., Particle Data Group Collaboration, Review of particle physics, Chin. Phys. C 38 (2014) 090001.

[21] F. Capozzi, G.L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo, Phys. Rev. D 89 (2014) 093018, arXiv:1312.2878 [hep-ph].

[22] T. Fukuyama, K. Ichikawa, Y. Mimura, Relation between proton decay and PMNS phase in the minimal SUSY SO (10) GUT, Phys. Lett. B 764 (2017) 114, arXiv:1609.08640 [hep-ph].

[23] B.C. Allanach, SOFTSUSY: a program for calculating supersymmetric spectra, Comput. Phys. Commun. 143 (2002) 305, arXiv:hep-ph/0104145.

[24] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, micrOMEGAs4.1: two dark matter candidates, Comput. Phys. Commun. 192 (2015) 322, arXiv:1407.6129 [hep-ph].

[25] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, micrOMEGAs: version 1.3, Comput. Phys. Commun. 174 (2006) 577, arXiv:hep-ph/0405253.

[26] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, MicrOMEGAs: a program for calculating the relic density in the MSSM, Comput. Phys. Commun. 149 (2002) 103, arXiv:hep-ph/0112278.

[27] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, K.E. Williams, HiggsBounds: confronting arbitrary Higgs sectors with exclusion bounds from LEP and the Tevatron, Comput. Phys. Commun. 181 (2010) 138, arXiv:0811.4169 [hep-ph].

[28] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, K.E. Williams, HiggsBounds 2.0.0: confronting neutral and charged Higgs sector predictions with exclusion bounds from LEP and the tevatron, Comput. Phys. Commun. 182 (2011) 2605, arXiv:1102.1898 [hep-ph].

[29] P. Bechtle, O. Brein, S. Heinemeyer, O. Stal, T. Stefaniak, G. Weiglein, K. Williams, Recent developments in HiggsBounds and a preview of HiggsSignals, PoS CHARGED 2012 (2012) 024, arXiv:1301.2345 [hep-ph].

[30] P. Bechtle, O. Brein, S. Heinemeyer, O. Stal, T. Stefaniak, G. Weiglein, K.E. Williams, HiggsBounds — 4: improved tests of extended Higgs sectors against exclusion bounds from LEP, the tevatron and the LHC, Eur. Phys. J. C 74 (3) (2014) 2693, arXiv:1311.0055 [hep-ph].

[31] S. Kraml, S. Kulkarni, U. Laa, A. Lessa, W. Magerl, D. Proschofsky-Spindler, W. Waltenberger, SModelS: a tool for interpreting simplified-model results from the LHC and its application to supersymmetry, Eur. Phys. J. C 74 (2014) 2868, arXiv:1312.4175 [hep-ph].

[32] G. Aad, et al., ATLAS Collaboration, CMS Collaboration, Combined measurement of the Higgs boson mass in pp collisions at -Js = 7 and 8 TeV with the ATLAS and CMS experiments, Phys. Rev. Lett. 114 (2015) 191803, arXiv:1503.07589 [hep-ex].

[33] G.W. Bennett, et al., Muon g-2 Collaboration, Final report of the muon E821 anomalous magnetic moment measurement at BNL, Phys. Rev. D 73 ( 2006) 072003, arXiv:hep-ex/0602035.

[34] B.L. Roberts, Status of the fermilab muon (g — 2) experiment, Chin. Phys. C 34 (2010) 741, arXiv:1001.2898 [hep-ex].

[35] K. Hagiwara, R. Liao, A.D. Martin, D. Nomura, T. Teubner, (g — 2)ß and a(M2Z) re-evaluated using new precise data, J. Phys. G 38 (2011) 085003, arXiv:1105.3149 [hep-ph].

[36] T. Aoyama, M. Hayakawa, T. Kinoshita, M. Nio, Complete tenth-order QED contribution to the muon g-2, Phys. Rev. Lett. 109 (2012) 111808, arXiv:1205.5370 [hep-ph].

[37] Y. Amhis, et al., Heavy Flavor Averaging Group Collaboration, Averages of B-hadron, C-hadron, and tau-lepton properties as of early 2012, arXiv:1207.1158 [hep-ex].

[38] C. Patrignani, et al., Particle Data Group Collaboration, Chin. Phys. C 40 (10) (2016) 100001.

[39] M. Antonelli, et al., FlaviaNet Working Group on Kaon Decays Collaboration, Precision tests of the standard model with leptonic and semileptonic kaon decays, arXiv:0801.1817 [hep-ph].

[40] P.A.R. Ade, et al., Planck Collaboration, Planck 2015 results. XIII. Cosmological parameters, arXiv:1502.01589 [astro-ph.CO].

[41] D.S. Akerib, et al., LUX Collaboration, Improved limits on scattering of weakly interacting massive particles from reanalysis of 2013 LUX data, Phys. Rev. Lett. 116 (16) (2016) 161301, arXiv:1512.03506 [astro-ph.CO];

D.S. Akerib, et al., Results from a search for dark matter in the complete LUX exposure, arXiv:1608.07648 [astro-ph.CO].

[42] M.G. Aartsen, et al., IceCube Collaboration, Improved limits on dark matter annihilation in the Sun with the 79-string IceCube detector and implications for supersymmetry, J. Cosmol. Astropart. Phys. 1604 (04) (2016) 022, arXiv:1601. 00653 [hep-ph].

[43] T. Moroi, The muon anomalous magnetic dipole moment in the minimal super-symmetric standard model, Phys. Rev. D 53 (1996) 6565, Erratum: Phys. Rev. D 56 (1997) 4424, arXiv:hep-ph/9512396.

[44] J. Hisano, T. Moroi, K. Tobe, M. Yamaguchi, T. Yanagida, Lepton flavor violation in the supersymmetric standard model with seesaw induced neutrino masses, Phys. Lett. B 357 (1995) 579, arXiv:hep-ph/9501407;

J. Hisano, T. Moroi, K. Tobe, M. Yamaguchi, Lepton flavor violation via right-handed neutrino Yukawa couplings in supersymmetric standard model, Phys. Rev. D 53 (1996 ) 2442, arXiv:hep-ph/9510309.

[45] T. Mori, MEG Collaboration, Final results of the MEG experiment, arXiv:1606. 08168 [hep-ex].

[46] S. Weinberg, Supersymmetry at ordinary energies. 1. Masses and conservation laws, Phys. Rev. D 26 (1982) 287;

N. Sakai, T. Yanagida, Proton decay in a class of supersymmetric grand unified models, Nucl. Phys. B 197 (1982) 533.

[47] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, N. 0kada, General formulation for proton decay rate in minimal supersymmetric S0(10) GUT, Eur. Phys. J. C 42 (2005) 191, arXiv:hep-ph/0401213.

[48] K. Abe, et al., Hyper-Kamiokande Working Group, A long baseline neutrino oscillation experiment using J-PARC neutrino beam and Hyper-Kamiokande, arXiv:1412.4673 [physics.ins-det];

K. Abe, et al., Hyper-Kamiokande Proto-Collaboration, Physics potential of a long-baseline neutrino oscillation experiment using a J-PARC neutrino beam and Hyper-Kamiokande, PTEP 2015 (2015) 053C02, http://dx.doi.org/10.1093/ ptep/ptv061, arXiv:1502.05199 [hep-ex];

K. Abe, et al., Hyper-Kamiokande Proto-Collaboration, Physics potentials with the second Hyper-Kamiokande detector in Korea, arXiv:1611.06118 [hep-ex].