Probing the supersymmetric grand unified theories with gravity mediation at the future proton–proton colliders and hyper-Kamiokande experiment

With the grand desert hypothesis, we have proposed to probe the supersymmetric Grand Unified Theories (GUTs) at the future proton–proton (pp) colliders and Hyper-Kamiokande experiment previously. In this paper, we study the supersymmetric GUTs with gravity mediated supersymmetry breaking in details. First, considering the dimension-six proton decay via heavy gauge boson exchange, we point out that we can probe the supersymmetric GUTs with GUT scale MGUT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{GUT}$$\end{document} up to 1.778×1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.778\times 10^{16}$$\end{document} GeV at the Hyper-Kamiokande experiment. Second, for the supersymmetric GUTs with MGUT≥1.0×1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{GUT} \ge 1.0\times 10^{16}$$\end{document} GeV and MGUT≥1.2×1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{GUT} \ge 1.2\times 10^{16}$$\end{document} GeV, we show that the upper bounds on the universal gaugino mass are 7.2 TeV and 3.5 TeV, respectively, and thus the corresponding upper bounds on gluino mass are 15 TeV and 8 TeV, respectively. Also, we shall study the masses for charginos, neutralinos, squarks, sleptons, and Higgs particles in details. In particular, the supersymmetric GUTs with MGUT≤1.2×1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{GUT} \le 1.2\times 10^{16}$$\end{document} GeV can be probed at the Hyper-Kamiokande experiment, and the supersymmetric GUTs with MGUT≥1.2×1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{GUT}\ge 1.2\times 10^{16}$$\end{document} GeV can be probed at the future 100 TeV pp collider experiments such as the FCChh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{FCC}_\textrm{hh}$$\end{document} and SppC via gluino searches. Thus, the supersymmetric GUTs with gravity mediation can be probed by the FCChh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{FCC}_\textrm{hh}$$\end{document}, SppC, and Hyper-Kamiokande experiments. In our previous study, we have shown that the supersymmetric GUTs with anomaly and gauge mediated supersymmetry breakings are well within the reaches of these experiments. Therefore, our proposal provides the concrete scientific goal for the FCChh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{FCC}_\textrm{hh}$$\end{document}, SppC, and Hyper-Kamiokande experiments: probing the supersymmetric GUTs.


I. INTRODUCTION
It is well known that supersymmetry (SUSY) provides a natural solution to the gauge hierarchy problem in the Standard Model (SM).In the supersymmetric SMs (SSMs) with Z 2 R-parity, we can achieve the gauge coupling unification [1], have the Lightest Supersymmetric Particle (LSP) like neutralino as a dark matter (DM) candidate [2], and break the electroweak (EW) gauge symmetry radiatively due to the large top quark Yukawa coupling, etc.In particular, gauge coupling unification strongly suggests the Grand Unified Theories (GUTs) [3][4][5][6], and the SUSY GUTs can be constructed from superstring theory, which is the most competitive candidate for quantum gravity.Therefore, supersymmetry is not only the most promising new physics beyond the SM, but also a bridge between the low energy phenomenology and high-energy fundamental physics.
On the other hand, to probe the new physics beyond the SM, we have a few proposals for the future protonproton (pp) colliders, for example, the FCC hh [26] and SppC [27].The naive question is whether we can probe the supersymmetry at the FCC hh and SppC, but the answer is no since the sparticles can be very heavy.And then the interesting question is whether we can probe the supersymmetric GUTs at the future experiments even if there does exist the SUSY EWFT problem.In a previous study [28], with the grand desert hypothesis, we showed that the supersymmetric GUTs can be probed at the future proton-proton (pp) colliders and Hyper-Kamiokande experiment.For the GUTs with the GUT scale M GU T ≤ 1.0 × 10 16 GeV, the dimension-six proton decay via heavy gauge boson exchange can be probed at the Hyper-Kamiokande experiment.Also, for the GUTs with M GU T ≥ 1.0 × 10 16 GeV, we showed that the GUTs with anomaly and gauge mediated supersymmetry breakings are well within the reaches of the future 100 TeV pp colliders such as the FCC hh and SppC, and the supersymmetric GUTs with gravity mediated supersymmetry breaking can be probed at the future 160 TeV pp collider.Therefore, the remaining interesting question is whether we can probe the supersymmetric GUTs with gravity mediated supersymmetry breaking at the FCC hh and SppC.
In this paper, we shall study the supersymmetric GUTs with gravity mediated supersymmetry breaking [29][30][31] in details.First, considering the dimension-six proton decay via heavy gauge boson exchange, we point out that we can probe the supersymmetric GUTs with GUT scale M GU T up to 1.778×10 16 GeV at the Hyper-Kamiokande experiment.Second, for the supersymmetric GUTs with M GU T ≥ 1.0 × 10 16 GeV and M GU T ≥ 1.2 × 10 16 GeV, we show that the upper bounds on the universal gaugino mass are 7.2 TeV and 3.5 TeV, respectively, and thus the corresponding upper bounds on gluino mass are 15 TeV and 8 TeV, respectively.Also, we shall study the masses for charginos, neutralinos, squarks, sleptons, and Higgs particles in details.In particular, the supersymmetric GUTs with M GU T ≤ 1.2×10 16 GeV can be probed at the Hyper-Kamiokande experiment, and the supersymmetric GUTs with M GU T ≥ 1.2 × 10 16 GeV can be probed at the FCC hh and SppC experiments via gluino searches.Thus, the supersymmetric GUTs with gravity mediation can be probed by the FCC hh , SppC, and Hyper-Kamiokande experiments.In our previous study, we have shown that the supersymmetric GUTs with anomaly and gauge mediated supersymmetry breakings are well within the reaches of these experiments [28].Therefore, we propose the concrete scientific goal for the FCC hh , SppC, and Hyper-Kamiokande experiments: probing the supersymmetric GUTs.
This paper is organized as follows.In Section II, we discuss the supersymmetric GUT searches at the Hyper-Kamiokande experiment.In Section III, we study the supersymmetric GUTs with gravity mediated supersymmetry breaking in details, and their searches at the future proton-proton colliders.Our conclusion is given in Section IV.

II. PROBING THE SUPERSYMMETRIC GRAND UNIFIED THEORIES AT THE HYPER-KAMIOKANDE EXPERIMENT
In the GUTs, the well-know prediction is the dimension-six proton decay p → e + π 0 via heavy gauge boson exchange, and the proton lifetime is given by [32,33] where A R is the dimensionless one-loop renormalization factor associated with anomalous dimension of the relevant baryon-number violating operators, α GUT is the unified gauge coupling, and M X/Y is the mass for the heavy gauge bosons X µ /Y µ .The current lower limit on the proton lifetime from the Super-Kamiokande experiment is τ p > 1.6×10 34 years [34].Thus, we obtain M X/Y ≥ 1.0 × 10 16 GeV.At the future Hyper-Kamiokande experiment, we can probe the proton lifetime at least above 1.0 × 10 35 years [35].Thus, the GUTs with M X/Y ≤ 1.778 × 10 16 GeV is within the reach of the future Hyper-Kamiokande experiment.For more detail related GUTs and M GU T related bosons see [36][37][38][39][40][41][42][43][44].
To clarify the subtle point, we want to emphasize that the mass of the heavy gauge bosons X µ /Y µ is smaller than or equal to the GUT scale M GU T .Thus, the supersymmetric GUTs with GUT scale up to 1.778 × 10 16 GeV can be probed at the future Hyper-Kamiokande experiment.

III. PROBING THE SUPERSYMMETRIC GRAND UNIFIED THEORIES WITH GRAVITY MEDIATION AT THE FUTURE PROTON-PROTON COLLIDERS
The supersymmetry searches at the 100 TeV pp colliders have been studied previously [26,27,45,[47][48][49].For the integrated luminosity 30 ab −1 , Wino via Bino decay, gluino g via heavy flavor decay, gluino via light flavor decay, first-two generation squarks q, and stop can be discovered for their masses up to about 6.5 TeV, 11 TeV, 17 TeV, 14 TeV, and 11 TeV, respectively.Moreover, if the gluino and first-two generation squark masses are similar, they can be probed up to 20 TeV.
By the way, the correlations between the low energy SUSY spectra and the GUT scale have been studied via the one-loop renormalization group equations before, and it was found that the bound from dimension-six proton decay already excludes the gluinos and Winos heavier than about 120 TeV and 40 TeV, respectively, if their mass ratio M 3 /M 2 is about 3 [50].In our paper, we employ the ISAJET 7.85 package [51] to perform the scan, which will give us more precise results.To be concrete, for the supersymmetric GUTs with M GU T ≤ 1.0 × 10 16 GeV, we find that the current bound from dimension-six proton decay excludes the gluinos and Winos heavier than 15 TeV and 6 TeV, respectively.Because Winos might decay via Higgsinos as the benchmark point 4 given in the following subsection B, we are not sure whether Wino is within the reach of the FCC hh and SppC experiments, which will be studied elsewhere.Moreover, we shall study the masses for charginos, neutralinos, squarks, sleptons, and Higgs particles as well.

A. Phenomenological Constraints and Scanning Procedure
We summarize our scanning procedure and the experimental constraints in this part of our article.The package ISAJET 7.85 [51] is employed to perform the random scans over the parameter space of mSUGRA/CMSSM.In this code, the weak scale values of the gauge and third generation Yukawa couplings are evolved to M GUT via the MSSM RGEs in the DR regularization scheme.It is to be noted that we do not strictly enforce the unification condition g 3 = g 1 = g 2 at M GUT .It is because a few percent deviation from unification can be assigned to the unknown GUT-scale threshold corrections [52].Moreover we do not allow g 3 to deviate from the unification by more than about 3%.We define the boundary conditions of our model at M GUT and then all the Soft Supersymmetry Breaking (SSB) parameters, along with the gauge and Yukawa couplings, are evolved back to the weak scale M Z .
In evaluating Yukawa couplings, the SUSY threshold corrections [53] are taken into account at the common scale M SUSY = √ m tL m tR .The entire parameter set is iteratively run between M Z and M GUT using the full two-loop RGEs until a stable solution is obtained.To better account for the leading-log corrections, one-loop step-beta functions are adopted for the gauge and Yukawa couplings, and the SSB parameters m i are extracted from RGEs at appropriate scales as m i = m i (m i ).The RGE-improved one-loop effective potential is minimized at an optimized scale M SUSY , which effectively accounts for the leading two-loop corrections.The full one-loop radiative corrections are incorporated for all sparticles.
Note that we set µ > 0 and use m t = 173.
We would like to draw attention of the reader to the fact that the requirement of radiative electroweak symmetry breaking (REWSB) [55] puts an important theoretical constraint on the parameter space.Another important constraint comes from the limits on the cosmological abundance of stable charged particles [56].They exclude the parameter space where the charged SUSY particles, such as τ1 or t1 , become the LSP.The data points collected all satisfy the requirement of REWSB, with the neutralino being the LSP.In addition, after collecting the data, we impose the mass bounds on all the sparticles [57], and the constraints from rare decay processes B s → µ + µ − [58], b → sγ [59], and B u → τ ν τ [60].More explicitly, we set To be general, we do not require the relic abundance of the LSP neutralino to satisfy the Planck 2018 bound within 5σ: 0.114 ≤ Ω CDM h 2 (Planck) ≤ 0.126 [61].

B. Scan Results
We shall discuss results from the systematical scans.In given by two vertical lines, we obtain that the upper bounds on gluino masses are 8 TeV and 15 TeV respectively for the red and blue points.Therefore, we clearly show that SUSY GUTs with M GU T 1.2 × 10 16 GeV for gravity mediated SUSY breaking scenario [29][30][31], i.e. the red points, can be probed by the future 100 TeV pp colliders such as the FCC hh and SppC.Moreover, the blue points and orange points with M GU T ≤ 1.0 × 10 16 GeV can be explored by the Hyper-Kamiokande experiment.In the latter part of the paper we see the impact of these bounds on the fundamental parameters of the mSUGRA/CMSSM and sparticle spectrum.
In Fig. 2, we display plots the mSUGRA/CMSSM fundamental parameters as function of gluino mass.The color In Fig. 3, we present the gluino mass versus the left-handed slepton mass, light stau mass, light stop mass, and CP-odd Higgs boson mass.In the top upper panel, we display plot in m ẽL − m g plane.As we stated earlier, large m 0 term can give dominant contributions to the squark and slepton masses, so this plot is very much similar to m 0 − m g .We do not show the plot for m ẽR − m g because it is similar to m ẽL − m g plot.Similarly, we depict m τ1 − m g plot in the top right panel in which also have similar feature to m ẽL − m g .Plot in m t1 − m g plane is shown in lower left panel.Here we see the similar trend but the mass ranges are reduced.Plot in lower right panel is shown in m A − m g plane.Here we notice that m A can be as heavy as 64 TeV for both red and blue points.
In Fig. 4, we show the gluino mass versus the light chargino mass and lightest neutralino mass, the light chargino mass versus the lightest neutralino mass, and dark matter density versus the lightest neutralino mass.Plot in the top left panel shows a graph in m χ± 2 productions with l-mediated decays and the productions with the SM-boson-mediated decays, which require m χ± 1 700 GeV.This means that most of our points can satisfy these bounds.Moreover, we will show later that for some of the lighter solutions m χ± 1 and m χ0 1 are mass-degenerate and then can evade these constraints as well.A plot in m χ0 1 − m g plane is shown in the top right panel.This plot also has similar features to the previous plot.Here, the neutralino mass can reached up to 1.1 TeV for red points and 2 TeV for blue points.In the lower left panel, we show plot in m χ0 1 − m χ± 1 plane.The diagonal black line is just to show the region of chargino-neutralino coannihilation where chargino and neutralino masses are degenerate, and most of these points might have the Higgsino LSP.For red and blue points, the coannihilation region has the upper mass bounds on m χ0 1 about 1 TeV and 1.9 TeV, respectively.The coannihilation region with lower mass range as stated before is safe from collider constraints.In the right lower panel, we displays the plot with the LSP neutralino mass versus dark matter relic density Ωh 2 .Horizontal two black lines represent Planck 2018 5σ bounds on cold dark matter relic density as shown above.It can be seen that some of the red and blue points can satisfy the relic density bounds.In the following, we present four benchmark points in Table I for the red and blue points, which are consistent with the dark matter relic density bounds.

IV. CONCLUSION
We have studied the supersymmetric GUTs with gravity mediated supersymmetry breaking in details.First, considering the dimension-six proton decay via heavy gauge boson exchange, we pointed out that the supersymmetric GUTs with GUT scale M GU T up to 1.778 × 10 16 GeV can be probed at the Hyper-Kamiokande experiment.Second, for the supersymmetric GUTs with M GU T ≥ 1.0 × 10 16 GeV and M GU T ≥ 1.2 × 10 16 GeV, we showed that the upper bounds on the universal gaugino mass are 7.2 TeV and 3.5 TeV, respectively, and thus the corresponding upper bounds on gluino mass are 15 TeV and 8 TeV, respectively.In particular, the supersymmetric GUTs with M GU T ≤ 1.2 × 10 16 GeV can be probed at the Hyper-Kamiokande experiment, and the supersymmetric GUTs with M GU T ≥ 1.2 × 10 16 GeV can be probed at the FCC hh and SppC experiments via gluino searches.Thus, the supersymmetric GUTs with gravity mediation can be probed by the FCC hh , SppC, and Hyper-Kamiokande experiments.In our previous study, we have shown that the supersymmetric GUTs with anomaly and gauge mediated supersymmetry breakings are well within the reaches of these experiments.Therefore, we propose the concrete scientific goal for the FCC hh , SppC, and Hyper-Kamiokande experiments: probing the supersymmetric GUTs.

3
FIG. 1. Plot in M 1/2 − MGUT plane.Gray points are consistent with REWSB and LSP neutralino.Orange points satisfy the mass bounds including m h = 125 ± 3 GeV and the constraints from rare B− meson decays.Blue points form a subset of orange points and satisfy 1 MGUT 1×10 16 GeV, while red points form a subset of orange points and satisfy MGUT 1.2×10 16 GeV.Two horizontal blue and red lines represent MGUT = 1 × 10 16 GeV and MGUT = 1.2 × 10 16 GeV, respectively.The first vertical line shows the upper bound on M 1/2 for red points (M 1/2 = 3.5 TeV), and the second vertical line shows the upper bound on M 1/2 for blue points (M 1/2 = 7.2 TeV).

FIG. 2 .
FIG. 2. The gluino mass versus supersymmetry breaking soft terms and tan β.The color coding is the same as in Fig. 1.Horizontal red, black and blue lines represent the gluino mass upper bounds of 8 TeV, 11 TeV, and 15 TeV corresponding to the red points, gluino discovery via heavy flavor decay at 100 TeV pp collider [45], and blue points, respectively.The first vertical line shows the upper bound on M 1/2 for red points (M 1/2 = 3.5 TeV), and the second vertical line shows the upper bound on M 1/2 for blue points (M 1/2 = 7.2 TeV).

FIG. 3 .
FIG. 3. The gluino mass versus the left-handed slepton mass, light stau mass, light stop mass, and CP-odd Higgs boson mass.The color coding is the same as in Fig. 2.

Fig. 1 ,
we show plot M GU T as a function of M 1/2 .Gray points are consistent with REWSB and LSP neutralino.Orange points satisfy the mass bounds including m h = 125 ± 3 GeV and the constraints from rare B− meson decays.Blue points form a subset of orange points and satisfy 1 M U 1 × 10 16 GeV, while red points form a subset of orange points and satisfy M GU T 1.2 × 10 16 GeV.Two horizontal blue and red lines represent M GU T = 1 × 10 16 GeV and M GU T = 1.2 × 10 16 GeV, respectively.The first vertical line shows the upper bound on M 1/2 for red points (M 1/2 = 3.5 TeV), while the second vertical line shows the upper bound on M 1/2 for blue points (M 1/2 = 7.2 TeV).From the upper bounds on gaugino masses M 1/2

FIG. 4 .
FIG.4.The gluino mass versus the light chargino mass and lightest neutralino mass, the light chargino mass versus the lightest neutralino mass, and dark matter density versus the lightest neutralino mass.The color coding is the same as in Fig.2.

TABLE I .
Four benchmark points for the red and blue points, which are consistent with the dark matter relic density bounds.Sparticle and Higgs masses are in GeV units.