Inspiration from Intersecting D-branes: General Supersymmetry Breaking Soft Terms in No-Scale ${\cal F}$-$SU(5)$

Motivated by D-brane model building, we evaluate the $\cal{F}$-$SU(5)$ model with additional vector-like particle multiplets, referred to as flippons, within the framework of No-Scale Supergravity with non-vanishing general supersymmetry breaking soft terms at the string scale. The viable phenomenology is uncovered by applying all current experimental constraints, including but not limited to the correct light Higgs boson mass, WMAP and Planck relic density measurements, and several LHC constraints on supersymmetric particle spectra. Four interesting regions of the parameter space arise, as well as mixed scenarios, given by: (i) light stop coannihilation; (ii) pure Higgsino dark matter; (iii) Higgs funnel; and (iv) light stau coannihilation. All regions can generate the observed value of the relic density commensurate with a 125 GeV light Higgs boson mass, with the exception of the relatively small relic density value for the pure Higgsino lightest supersymmetric particle (LSP). This work is concluded by gauging the model against present LHC search constraints and derivation of the final states observable at the LHC for each of these scenarios.


INTRODUCTION
Successful confirmation of the Standard Model (SM) was celebrated when the lightest CP-even Higgs boson with mass m h = 125.09 ± 0.24 GeV was discovered at the LHC in 2012 [1,2]. Despite the memorable occasion, severe anomalies persisted in the SM, for instance, the gauge hierarchy problem, conspicuous absence of gauge coupling unification, and lack of a plausible dark matter candidate, just to highlight a few. Beyond the SM (BSM) theories, and supersymmetry in particular, can rescue high-energy physics from these SM deficiencies. Supersymmetry (SUSY) can solve the gauge hierarchy problem and produce gauge coupling unification. In supersymmetric SMs (SSMs), the large top quark Yukawa coupling can radiatively break electroweak (EW) gauge symmetry. The lightest supersymmetric particle (LSP) neutralino ( χ 0 1 ) can serve as a viable dark matter candidate in SSMs with R-parity. Of particular significance, gauge coupling unification strongly suggests Grand Unified Theories (GUTs), and SUSY GUTs can be elegantly constructed from superstring theory. Supersymmetry thus provides rather auspicious new physics beyond the SM, beautifully merging low energy phenomenology with high-energy fundamental physics. Despite the successful intervention SUSY can inject into high-energy physics, it remains considerably challenging to obtain a light Higgs boson mass around 125 GeV in the Minimal SSM (MSSM) without embracing either multi-TeV top squarks with small mixing or TeV-scale top squarks with large mixing [3]. Compounding the effort are the strong constraints on the SSM viable parameter space established by the LHC SUSY searches. Principally among those LHC constraints are the masses for the gluino (g) and light stop (t 1 ), where the nominal exclusion curves imply the masses are heavier than about 1.9 TeV and 900 GeV, respectively [4], indicating that an electroweak fine-tuning problem may potentially lurk in SUSY.
The aforementioned threats to developing effective BSM constructions notwithstanding, string theory perseveres as one of the most promising theories for quantum gravity. As opposed to the conventional unification at the GUT scale realized in the typical SUSY GUT, we proposed the testable flipped SU (5)×U (1) X models [5][6][7] with additional TeV-scale vector-like particles [8], which we colorfully dubbed flippons, to obtain gauge coupling unification at the string scale. Subsequently, we further constructed this class of flipped SU (5) models from local F-theory model building [9,10]. These flipped SU (5) models with extra vector-like multiplets can be realized in free-fermionic string constructions also [11], hence we denoted them F -SU (5). Let's enumerate the "Miracles" [12] of the flippons in F -SU (5): (1) The lightest CP-even Higgs boson mass can be easily lifted to 125 GeV due to one-loop contributions from the Yukawa couplings between the flippons and Higgs fields [12,13].
(2) It is well-known that dimension-five proton decays mediated by colored Higgsinos are highly suppressed due to the missing partner mechanism and TeV-scale µ term in the flipped SU (5) × U (1) X models. In F -SU (5), the SU (3) C × SU (2) L gauge couplings do in fact unify at the traditional GUT scale while the two unified gauge couplings increase as a result of the vectorlike particle contributions [14,15]. Therefore, dimensionsix proton decays via heavy gauge boson exchanges are within the reach of the future proton decay experiments such as the Hyper-Kamiokande experiment. More concisely, F -SU (5) models differ from the minimal flipped SU (5) × U (1) X model since the proton lifetime in the minimal model is too lengthy for future proton decay experiments.
(3) The lightest neutralino serves as the LSP and is lighter than the light stau attributable to the longer running of the Renormalization Group Equations (RGEs) in No-Scale supergravity [16], allowing the LSP neutralino to prevail as a dark matter candidate [17][18][19]. More acutely, No-Scale F -SU (5) yields the uncommon mass hierarchy M (t 1 ) < M (g) < M (q) of a light stop and gluino both substantially lighter than all other squarks (q) [17][18][19]. The net effect of this rare SUSY spectrum is the production of four top quarks, leading to large multijet events at the LHC [20].
Several prior No-Scale F -SU (5) analyses have intimately examined vanishing SUSY breaking soft terms with the exception of a single unified gaugino parameter M 1/2 (see, for example, Refs. [12,24,25]). To complement those prior sweeping studies, we now consider in this paper non-zero general SUSY breaking soft terms in No-Scale F -SU (5), partially inspired by D-brane model building [26]. The low-energy particle spectra will be examined that are consistent with all the current experimental constraints, and interesting diverse regions of the parameter space that can generate the observed dark matter relic density and correct light Higgs boson mass, in addition to several other currently operating experiments, will be identified and analyzed discretely to derive low-energy phenomenology. Finally, benchmarks models will be classified and branching fractions computed to itemize the final states observable at the LHC. Here we only briefly review the minimal flipped SU (5) model [5][6][7], where the gauge group SU (5)×U (1) X can be embedded into the SO(10) model. More comprehensive discussions of the minimal flipped SU (5) model can be found in Refs. [12,20,23,24,27] and references therein. We define the generator U (1) Y ′ in SU (5) as and subsequently the hypercharge is given by We have three families of the SM fermions whose quantum numbers under SU (5) × U (1) X are respectively where i = 1, 2, 3. The SM particle assignments in F i ,f i andl i are where are the left-handed quark doublets, right-handed up-type quarks, downtype quarks, left-handed lepton doublets, right-handed charged leptons, and neutrinos, respectively. Generation of the heavy right-handed neutrino masses is accomplished by introduction of three SM singlets φ i .
Breaking the GUT and electroweak gauge symmetries is achieved via introduction of the two pairs of Higgs representations H = (10, 1), H = (10, −1), The states in the H multiplet are labeled by the same symbols as in the F multiplet, and for H we merely add "bar" above the fields. Explicitly, the Higgs particles are where H d and H u are one pair of Higgs doublets in the MSSM.
The SU (5) × U (1) X gauge symmetry is broken down to the SM gauge symmetry by the following Higgs superpotential at the GUT scale respectively with D h and D h , forming massive eigenstates with masses 2λ 1 < N c H > and 2λ 2 < N c H >. As a result, we naturally experience the doublet-triplet splitting due to the missing partner mechanism [7]. The triplets in h and h though only have a small mixing through the µ term, and as such, the colored Higgsinoexchange mediated proton decay is negligible, i.e., there is no dimension-5 proton decay problem.
To realize string-scale gauge coupling unification [8][9][10], we introduce the following vector-like particles (flippons) at the TeV scale XF = (10, 1) , XF = (10, −1) , The particle content from the decompositions of XF , XF , Xl, and Xl under the SM gauge symmetry are The quantum numbers for the extra vector-like particles under the SU The superpotential is and after the SU (5) × U (1) X gauge symmetry is broken down to the SM gauge symmetry, the above superpotential gives where y D ij , y Uν ij , y E ij , y N ij , y XF , and y XF are Yukawa couplings, µ is the bilinear Higgs mass term, and M φ ij , M XF and M Xl are masses for new particles. These new particles are of course our flippons, however, we shall not explicitly compute the masses M φ ij , M XF ,and M Xl here in this work, reserving that project for a future date. Nonetheless, we do implement a common mass decoupling scale M V for the flippon vector-like particles. Current LHC constraints on vector-like T and B quarks [28] provide lower limits of around 855 GeV for our (XQ, XQ c ) vector-like flippons and 735 GeV for our (XD, XD c ) vector-like flippons. Accordingly, we establish our lower M V limit at M V ≥ 855 GeV to ensure complete coverage of all experimentally viable flippon masses in our analysis.
As summarized in the prior section, the splitunification of flipped SU (5) [5][6][7] provides for fundamental GUT scale Higgs representations (not adjoints), natural doublet-triplet splitting, suppression of dimension-five proton decay [29], and a two-step see-saw mechanism for neutrino masses [30,31]. Additions to the one-loop gauge β-function coefficients b i to include contributions from the vector-like flippon multiplets induce the necessary flattening of the SU (3) Renormalization Group Equation (RGE) running (b 3 = 0) [17], translating into a clear separation between the primary SU (3) C × SU (2) L unification near 10 16 GeV and the secondary SU (5) × U (1) X unification at around 5×10 17 GeV, which we define as the M F scale, thus elevating unification to near the Planck mass. At the primary SU (3) C × SU (2) L unification near 10 16 GeV, the M 2 and M 3 gaugino mass terms are unified into a single term that we refer to as M 5 [32], where M 5 = M 2 = M 3 between the primary unification around 10 16 GeV and the secondary SU (5) × U (1) X unification near 5 × 10 17 GeV [17]. The M 1 gaugino mass term runs up to the secondary SU (5) × U (1) X unification at M F and unifies with M 5, by means of a slight shift due to U (1) X flux effects [32] between the primary unification around 10 16 GeV and the secondary SU (5) × U (1) X unification at M F [17]. Given this shift, we refer to the M 1 gaugino mass term above the primary unification around 10 16 GeV as M 1X [32]. The resulting baseline extension for logarithmic running of the No-Scale boundary conditions permits sufficient scale for natural dynamic evolution down to viable phenomenology at the electroweak scale. It is this associated flattening of the color-charged gaugino mass scale that produces the distinctive mass texture of M ( t 1 ) < M ( g) < M ( q), generating a light stop and gluino that are lighter than all other squarks [12]. Therefore, the most general supersymmetry breaking soft terms at the scale M F in our F -SU (5) model are and A b . General supersymmetry breaking soft terms of these type are partially inspired by D-brane model building [26], where F i , f i , l i , and h/h arise from intersections I: The F-SU (5) general SUSY breaking soft terms, in addition to the vector-like flippon decoupling scale MV , the low energy ratio of Higgs vacuum expectation values (VEVs) tanβ, and top quark mass Mt, for the four regions of the model space we study in this work. Each benchmark point is identified with an alphabetical label in order to link the data in TABLE I with the data in TABLES II -III. All masses are in GeV.
, rescaled dark matter directdetection cross sections (σSI , σSD), and p → e + π 0 proton decay rates (τp) for the F-SU (5) general SUSY breaking soft terms of TABLE I. The soft SUSY breaking terms that generate each of these values can be identified by the alphabetical label. The σSI and σSD cross-sections have been rescaled in accordance with Eq. (17). The numerical values given for ∆aµ are ×10 −10 , , rescaled spin-independent cross-sections σSI are ×10 −11 pb, rescaled spin-dependent cross-sections σSD are ×10 −9 pb, and proton decay rate τp are ×10 35 yrs.

Model
Benchmark GeV, 2 ≤ tanβ ≤ 60, and 855 ≤ M V ≤ 20, 000 GeV. A reasonable tolerance of the top quark mass is permitted around the world average [33], employing upper and lower limits of 171.8 ≤ M t ≤ 174.8 GeV. The WMAP 9-year [34] and 2015 Planck [35] relic density measurements are implemented, such that we constrain the model to satisfy both sets of data and allow the inclusion of multi-component dark matter beyond the neutralino, imposing limits of Ωh 2 ≤ 0.1221. We will further identify the subset of points that are consistent with the recently released 2018 Planck observation [36] of Ωh 2 = 0.120 ± 0.001. Regarding LHC gluino searches, a firm lower limit is imposed on the gluino mass of M ( g) ≥ 1.6 TeV, allowing for a reasonable lower boundary under all the LHC gluino searches.
The light Higgs boson mass theoretical calculation is allowed to float around the experimental central value of m h = 125.09 GeV, though we do enforce a larger range of 123 ≤ m h ≤ 128 GeV to account for a 2σ experimental uncertainty over and above a theoretical uncertainty of 1.5 GeV in our computations. The actual numerical value of the flippon Yukawa coupling remains an unknown, hence we permit the Yukawa coupling to range from a minimum value to its maximum in our light Higgs boson mass computations. A minimal flippon Yukawa coupling does not allow any vector-like flippon contributions, as the light Higgs boson mass in this case is then comprised of only the 1-loop and 2-loop SUSY contributions, chiefly from the coupling to the light stop. A maximum flippon Yukawa coupling implies the (XD, XD c ) flippon Yukawa coupling is fixed at Y XD = 0 and the (XU, XU c ) flippon Yukawa coupling is set at Y XU = 1, with the (XD, XD c ) flippon trilinear coupling A term set at A XD = 0 and the (XU, XU c ) A term fixed at A XU = A U = A 0 [12,13]. When necessary, we shall choose the maximum flippon Yukawa coupling in order to lift the light Higgs boson mass up to its observed value. However, some points in the model space can reach the Higgs boson observed value with no flippon contribution, thus we choose the minimal coupling in these instances. Regardless of our choice for the numerical value of the flippon Yukawa coupling, the calculation must return a value within 123 ≤ m h ≤ 128 GeV. For our chosen benchmarks points that are meant to be sufficiently representative of the model space, we will clearly annotate as to whether the given light Higgs boson mass consists of the minimum or maximum flippon Yukawa coupling.
A total of 110 million points were randomly sampled in scans implementing the SU (5) × U (1) X M F scale boundary conditions, with only about 37,000 of those points surviving the applied constraints. The SUSY mass spectra, relic density, rare decay processes, and direct dark matter detection cross-sections are calculated with MicrOMEGAs 2.1 [46] utilizing a proprietary mpi modification of the SuSpect 2.34 [47] codebase to run flippon and General No-Scale F -SU (5) enhanced RGEs, employing non-universal soft supersymmetry breaking parameters at the M F scale. Supersymmetric particle decays are computed with SUSY-HIT 1.5a [48]. The Particle Data Group [49] world average for the strong coupling constant is α S (M Z ) = 0.1181 ± 0.0011 at 1σ, and we assume a value in this work of α S = 0.1184. Notice that there is no intersection between the Light Stop Coannihilation and pure Higgsino regions, though a more detailed analysis relaxing the higgsino and light stop coannihilation requirements we applied here is in progress [50] to ascertain any potential union of these two significant regions. All points plot satisfy the light Higgs boson mass and relic density observations we outlined in this work, in addition to rare decay and current LHC SUSY search constraints. The round black dots represent those points that can also satisfy the recent 2018 Planck Collaboration satellite relic density measurements of Ωh 2 = 0.120 ± 0.001.
The data results are distributed across TABLES I -III. The benchmarks models of TABLES I -III originate via categorization of the viable model space into four dissimilar regions identified by LSP and NLSP characteristics, to be discussed shortly. The numerical relic density figures annotated in TA-BLE III consist only of the SUSY lightest neutralino χ 0 1 abundance, therefore regions with numerical values less than the WMAP9 1σ measurement lower bound of about Ωh 2 ≃ 0.1093 must welcome alternative additions to the total observed relic density by WMAP9 and Planck. Given the potential for multi-component dark matter in regions of low neutralino density, the spin-dependent and spin-independent cross-section calculations provided in

PHENOMENOLOGICAL RESULTS
The model space is constrained per the experimental constraints detailed in the prior section, with the exception of the LUX, PandaX-II, and XENON cross-sections, which we shall evaluate independently of all the other empirical measurements. Post application of the constraints, the surviving viable parameter space consists of four interesting regions segregated by LSP characteristics 1 , clearly indicating the light stop coannihilation strip, with the lower end severed by the most recent LHC constraints. Given that the small mass delta we require is less than M t , the light stop is forced into the decay channel t 1 → c + χ 0 1 and the gluino will always produce a light stop ( g → t 1 t), as highlighted in TABLE IV. The near degenerate LSP and light stop in this STOPC channel can evade the current ATLAS [51,52] and CMS [53][54][55][56]   The higgsino spectra are identified by the Higgs bilinear mass term µ driven via RGE running below the gaugino mass terms M 1 and M 2 at low-energy near electroweak symmetry breaking (EWSB), thus driving the χ 0 2 mass negative. This is despite the fact that the µ term can be larger or smaller than the gaugino mass terms M 5 and M 1X at the M F scale. The tight constraints applied within the higgsino region ensure an almost 100% higgsino LSP, though it apparently does prevent any combination of a pure higgsino LSP and light stop NLSP. This is evident in FIG. 1 by observing the gap between the HLSP and STOPC points, as no spectra with GeV can also generate a pure higgsino LSP. This notwithstanding, preliminary studies do show that a STOPC point with M ( t 1 )−M ( χ 0 1 ) ≤ 50 GeV can possess an LSP that is dominantly higgsino, for example, more than 60% and potentially as high as 80% or greater. While we do not integrate this more detailed HLSP+STOPC analysis into this work, it is currently in progress [50]. Such a combination of a dominantly higgsino LSP and light stop NLSP is rather difficult to produce naturally [57], though it has been uncovered in certain realistic intersecting D6-brane models [58]. As a result of our tight constraints on the higgsino in order to generate a pure higgsino LSP, most of the higgsino points in FIGs. 1 -4 have M ( t 1 ) > M ( g), thus we focus on these primary higgsino spectra here in this work. Unfortunately, as depicted in TABLE IV, there is no dominant decay channel to consistent final states for either the light stop or gluino in the HLSP region. As can be inferred from TABLE IV and our constraint of GeV, the chargino participates in the off-shell Standard Model boson mediated decays χ ± 1 → ud/cs+ χ 0 1 (36%/36%) or χ ± 1 → l ± +ν l + χ 0 1 (28%), and these pure higgsino points have an LHC established lower bound of about M ( χ ± 1 ) ∼ 140 GeV [59], thus we enforce M ( χ ± 1 ) > 140 GeV in FIGs. 1 -4. The chargino χ ± 1 is the NLSP for 99.6% of the Higgs Funnel points, with the remaining 0.4% possessing a stau NLSP. Nonetheless, only 0.07% of the Higgs Funnel also resides in the M ( τ ± 1 ) − M ( χ 0 1 ) ≤ 20 GeV stau coannihilation strip. We further observe that 11.5% of the Higgs Funnel has a dominant higgsino LSP, though only 0.5% are pure higgsino. The upper right plot space of FIG. 3 clearly shows the HF region with M (H 0 /A 0 ) plot as a function of the LSP χ 0 1 . The minimum heavy Higgs boson mass in the HF region is M (H 0 /A 0 ) ∼ 800 GeV, and these H 0 /A 0 near this lower boundary all have 30 ≤ tanβ ≤ 40, and given a branching ratio of more than 85% for H 0 /A 0 → bb, the lightest M (H 0 /A 0 ) in our model space persist comfortably beyond current LHC constraints for BSM Higgs searches [60].
In our light stau coannihilation strip, only 27.6% of the points have stau NLSP, with the remaining 72.4% having a chargino χ ± 1 NLSP. As one would expect though, four nearly degenerate sparticles ( χ 0 1 , χ ± 1 , χ 0 2 , τ ± 1 ) greatly suppresses the relic density, hence only the stau τ ± 1 NLSP points with larger χ ± 1 and χ 0 2 in the light stau coannihilation strip can generate the observed WMAP9 and Planck relic density. The light stau coannihilation strip IV: Dominant gluino ( g) and light stop ( t1) decay modes for each of the regions of the model space we study in this work, along with the associated branching ratios. Here, q = (u, d, c, s) and q = ( u, d, c, s).

Model
Dominant Decay Mode Branching Ratio is evident in the lower left panel of FIG. 2. We inspect the most recent ATLAS slepton constraints [59] on compressed spectra, such as we have in the STAUC region, therefore we attempt to approximately emulate the AT-LAS slepton exclusions by applying the following constraints on our model space, where both conditions for each must be TRUE for the point to be excluded: (1) M ( τ ± 1 ) < 100 GeV and M ( τ ± 1 ) − M ( χ 0 1 ) > 1 GeV; (2) 100 ≤ M ( τ ± 1 ) < 150 GeV and 2 ≤ M ( τ ± 1 ) − M ( χ 0 1 ) ≤ 14 GeV; and (3) 150 ≤ M ( τ ± 1 ) < 190 GeV and 3 ≤ M ( τ ± 1 )−M ( χ 0 1 ) ≤ 8 GeV. These restrictions are included in FIGs. 1 -4. Notice in constraint (1) just above that we allow for an off-shell tau lepton τ ± less than 1 GeV, an action that does capture a not insignificant number of viable points, in view of the fact that 10.5% of the STAUC region has M ( τ ± 1 ) − M ( χ 0 1 ) < 1 GeV. A quick review of TABLE IV indicates that the stau coannihilation spectra with light stau NLSP and M ( g) − M ( t 1 ) > M t produce exclusively a 4-top signature ( g → t 1 t → tt + χ 0 1 ), giving rise to large multijet events, as is typical in F -SU (5) [20]. The final states are not as clear when M ( g)−M ( t 1 ) < M t or M ( g) < M ( t 1 ) given that there is no dominant decay channel in either of these two cases. The model space is generally split between these two situations, where about half the STAUC region has M ( g) − M ( t 1 ) > M t and the other half of the STAUC region has M ( g) − M ( t 1 ) < M t or M ( g) < M ( t 1 ).

CONCLUSION
With our motivation partially inspired by D-brane model building, we presented an analysis of non-zero general SUSY breaking soft terms in F -SU (5). The methodology involved massive parallel computing given the large number of unknown parameters. While the resulting viable parameter space was extraordinarily constrained given the large number of computations, several interesting regions were uncovered that can simultaneously generate the WMAP and Planck observed relic density measurements and correct light Higgs boson mass, over and above satisfying many LHC search constraints. Four regions that accomplish these rely upon (i) light stop coannihilation, (ii) pure Higgsino dark matter, (iii) Higgs funnel, and (iv) light stau coannihilation, though the pure Higgsino LSP could not reach the observed relic density due to its large annihilation cross-section. Concluding the analysis was an effort to identify the decay modes to final states that could represent observable signatures of these four scenarios, discovering that the light stop coannihilation and light stau coannihilation decay channels possess very dominant branching fractions. The light stop coannihilation produces the typical top+charm quark final state via the gluino decay and the light stau coannihilation mostly leads to the characteristic F -SU (5) large multijet event, whereas the pure Higgsino LSP and Higgs funnel provide no dominant decay channel to final states. Though possibly furnishing the highest intrigue, our exploration tentatively revealed that the F -SU (5) model may indeed harbor a rather diminutive region exhibiting spectra with a mixed scenario of a dominant Higgsino LSP and light stop coannihilation, a rare yet rather natural SUSY spectrum, but that is the focus of our next endeavor.