Yukawa-unified natural supersymmetry

Previous work on t − b − τ Yukawa-unified supersymmetry, as expected from SUSY GUT theories based on the gauge group SO(10), tended to have exceedingly large electroweak fine-tuning (EWFT). Here, we examine supersymmetric models where we simultaneously require low EWFT (“natural SUSY”) and a high degree of Yukawa coupling unification, along with a light Higgs scalar with mh ~ 125 GeV. As Yukawa unification requires large tan β ~ 50, while EWFT requires rather light third generation squarks and low μ ≈ 100 − 250 GeV, B-physics constraints from BR(B → Xsγ) and BR(Bs → μ+μ−) can be severe. We are able to find models with EWFT Δ ≲ 50 − 100 (better than 1–2% EWFT) and with Yukawa unification as low as Ryuk ~ 1.2 (20% unification). The unification is lessened to Ryuk ~ 1.3 when B-physics constraints are imposed. We present several Yukawa-unified natural SUSY (YUNS) benchmark points. LHC searches will be able to access gluinos in the lower 1 – 2 TeV portion of their predicted mass range although much of YUNS parameter space may lie beyond LHC14 reach. If heavy Higgs bosons can be accessed at a high rate, then the rare H, A → μ+μ− decay might allow a determination of tan β ~ 50 as predicted by YUNS models. Finally, the predicted light higgsinos should be accessible to a linear e+e− collider with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$ \sqrt{s}\sim 0.5 $\end{document} TeV.


Introduction
A striking feature in nature is that all the fermions of each generation fill out a complete 16-dimensional spinor multiplet of the gauge group SO(10) [1][2][3][4][5]. While ordinary grand unified theories (GUTs) suffer from the notorious gauge hierarchy problem, supersymmetric (SUSY) GUTs not only tame this hierarchy problem [6,7], but they also receive support from the well-known unification of gauge couplings [8][9][10][11][12][13]. In the simplest SO(10) SUSY GUT theories, where both MSSM Higgs doublets H u and H d occupy the same 10-dimensional representation, one also expects unification of third generation Yukawa couplings f t , f b and f τ at M GUT ≃ 2 × 10 16 GeV . The t − b − τ Yukawa coupling unification is highly sensitive to both 2-loop renormalization group running (RGEs) and to threshold corrections when transitioning between MSSM and SM effective theories at the SUSY particle mass scale. Thus, the entire SUSY mass spectrum enters into a precise computation of Yukawa coupling unification. Many groups have explored t−b−τ Yukawa unification (YU) in SUSY theories [38][39][40][41][42][43][44][45][46][47][48][49][50]. It has been found that, for µ > 0, YU can occur at the few percent level in either the Higgs splitting (HS) model or the DR3 model (D-term splitting, right-hand neutrino effects and third generation splitting), provided that the GUT scale soft SUSY breaking (SSB) terms are related as A 2 0 ≃ 2m 2 10 ≃ 4m 2 16 , (1.1) with A 0 < 0. Moreover, the GUT scale Higgs splitting m 2 Hu < m 2 H d is needed to allow for an appropriate radiative breakdown of electroweak symmetry. In these models, first generation squarks and sleptons are required to be in the multi-TeV range [41,46] JHEP12(2012)066 while third generation sfermions are driven to much lighter TeV-scale masses. A benefit of these models is that the light SUSY Higgs boson mass m h tends naturally to be in the 125 GeV range, 1 as required by the recent LHC discovery [53,54]. The gaugino masses from YU SUSY are expected to be quite light, with typically mg 500 GeV (now excluded by LHC searches for gluino pair production along with cascade decay into states containing b-quarks [55]), although solutions with heavier gluinos ∼ 1-2 TeV can also be found [56]. In all these cases, the superpotential µ parameter, which is extracted from the electroweak minimization conditions, occupies values in the 1-10 TeV regime, leading to severe electroweak fine-tuning (EWFT).
In this work, we examine to what extent it is possible to reconcile YU with low EWFT. Minimization of the SUSY scalar potential allows one to relate the Z mass scale to the superpartner mass scale via the well-known relation (1. 2) The radiative corrections Σ u and Σ d are given in the 1-loop approximation of the Higgs effective potential by: where ∆V is the one-loop correction to the tree-level potential, and the derivative is evaluated in the physical vacuum: i.e. the fields are set to their vacuum expectation values after evaluating the derivative. At the one-loop level and in the limit of setting first/second generation Yukawa couplings to zero, Σ u contains 18 and Σ d contains 19 separate contributions from various particles/sparticles [57]. We include contributions from W ± , Z,t 1,2 , b 1,2 ,τ 1,2 , W 1,2 , Z 1,2,3,4 , t, b and τ , h, H and H ± . We adopt a scale choice Q 2 = mt 1 mt 2 to minimize the largest of the logarithms. The dominant contribution to the terms Σ u,d arise from superpotential Yukawa interactions of third generation squarks involving the top quark Yukawa coupling. For instance, the dominant contribution to Σ u is given by This expression thus grows quadratically with the stop mass. We adopt the fine-tuning measure from [57], which requires that each of the 40 terms on the right-hand-side (r.h.s. ) of eq. (1.2) should be of order ∼ m 2 where Λ max ∼ 100 − 300 GeV, depending on how much EWFT one is willing to tolerate. This measure of fine-tuning is similar to (but not exactly the same as) Kitano-Nomura [58,59] but different from Barbieri-Giudice [60] beyond the tree-level. In the JHEP12(2012)066 following, we will use the fine-tuning parameter where lower values of ∆ correspond to less fine-tuning, and e.g. ∆ = 20 would correspond to ∆ −1 = 5% fine-tuning. Our goal in this paper is to search for parameter choices which 1. Maximize the degree of Yukawa coupling unification, i.e. minimize with each Yukawa coupling evaluated at the GUT scale. Thus, a value of R yuk = 1 would give perfect Yukawa coupling unification.
2. Have as low EWFT as possible (in practice, we will require ∆ 100, or better than 1% EWFT).

3.
Have m h ≈ 125 GeV in accord with the recent LHC discovery of a Higgs-like resonance. In practice, we will require 122 GeV < m h < 128 GeV to allow for a roughly 2-3 GeV error in the RG-improved one-loop effective potential calculation of the Higgs mass m h .
We recognize that in addition to EWFT, there also exists a fine-tuning associated with generating particular weak scale SUSY spectra from distinct GUT scale parameters [60], (see also [61,62] for related discussions). Here we adopt the less restrictive weak scale finetuning condition, which nonetheless turns out to be indeed very restrictive. In this vein, we regard particular GUT scale parameters as merely a parametrization of our ignorance of the mechanism of SUSY breaking and soft term generation. 2 For instance, in this paper we will require rather low values of superpotential µ parameter to avoid excessive EWFT in eq. (1.2). In generic SUSY models, the value of µ is expected to be of order M Planck since it is a dimensionful SUSY-preserving parameter. Excessively large µ can be avoided ala Giudice-Masiero [65] where the µ superpotential term is forbidden by some high scale symmetry, but then is regenerated as a soft SUSY breaking term via a Higgs-Higgs coupling to the hidden sector. For the remainder of this paper, in section 2 we present details of our scan over SUSY parameter space, and which constraints are invoked in our analysis. In section 3, we present the results of our parameter space scans. We will find that requiring ∆ 100 and m h ∼ 125 GeV allows for R yuk 1.3 once all applicable constraints are taken into account. While this degree of Yukawa unification is not optimal, we feel it is still useful in that it might guide model builders towards models including additional GUT scale threshold corrections or extra matter or above-GUT-scale running which may ameliorate the situation. In section 5, we discuss observable consequences of YUNS for LHC, ILC and dark matter searches. We pay some attention to methods which might allow one to distinguish YUNS from generic NS models at lower tan β values. In section 6 we present a summary and conclusions.

JHEP12(2012)066
2 Parameter space and Yukawa unification For our calculations, we adopt the ISAJET 7.83 [66] SUSY spectrum generator ISAS-UGRA [67,68]. ISASUGRA begins the calculation of the sparticle mass spectrum with input DR gauge couplings and f b , f τ Yukawa couplings at the scale Q = M Z (f t running begins at Q = m t ) and evolves the 6 couplings up in energy to scale Q = M GUT (defined as the value Q where g 1 = g 2 ) using two-loop RGEs. At Q = M GUT , we input the soft SUSY breaking parameters as boundary conditions, and evolve the set of 26 coupled MSSM RGEs [69] back down in scale to Q = M Z . Full two-loop MSSM RGEs are used for soft term evolution, while the gauge and Yukawa coupling evolution includes threshold effects in the one-loop beta-functions, so the gauge and Yukawa couplings transition smoothly from the MSSM to SM effective theories as different mass thresholds are passed. In ISASUGRA, the values of SSB terms of sparticles which mix are frozen out at the scale Q ≡ M SUSY = √ mt L mt R , while non-mixing SSB terms are frozen out at their own mass scale [67,68]. The scalar potential is minimized using the RG-improved one-loop MSSM effective potential evaluated at an optimized scale Q = M SUSY which accounts for leading two-loop effects [70]. Once the tree-level sparticle mass spectrum is computed, full one-loop radiative corrections are calculated for all sparticle and Higgs boson masses, including complete one-loop weak scale threshold corrections for the top, bottom and tau masses at scale Q = M SUSY [71][72][73][74]. Since the GUT scale Yukawa couplings are modified by the threshold corrections, the ISAJET RGE solution must be imposed iteratively with successive updown running until a convergent sparticle mass solution is found. Since ISASUGRA uses a "tower of effective theories" approach to RG evolution, we expect a more accurate evaluation of the sparticle mass spectrum for models with split spectra (this procedure sums the logarithms of potentially large ratios of sparticle masses) than with programs which make an all-at-once transition from the MSSM to SM effective theories. The fine-tuning measure ∆ described in section 1 has been implemented in ISAJET 7.83 [66].
In models of "natural SUSY" (NS) [57][58][59][75][76][77][78][79][80][81], the first requirement to gain a low EWFT is that the µ parameter be of the order of ∼ M Z , while in Yukawa-unified SUSY, it is necessary to invoke some manner of Higgs soft term splitting at the GUT scale in order to obtain raditive EWSB. For these reasons, the model parameter space chosen in this study is the two-parameter non-universal Higgs model (NUHM2) where the weak scale values of µ and m A are input in lieu of the GUT scale values of m 2 Hu and m 2 H d . In addition, to allow for (sub)TeV-scale third generation masses (as required by EWFT from eq. (1.2) along with at least a partial decoupling solution to the SUSY flavor and CP problems, we allow for split first/second and third generations at the GUT scale. Thus, the parameter space we choose is given by Here, m 16 (1, 2) and m 16 (3) are the first/second and third generation sfermion soft masses, respectively; m 1/2 ≡ M 1 = M 2 = M 3 is the universal gaugino mass parameter; and We search for mass spectra with low EWFT ∆ and low R yuk by performing a vast random scan over the following parameter ranges (masses in GeV units): 100 GeV < µ < 250 GeV, 500 GeV < m A < 5000 GeV.
The lower limit on m 1/2 comes from the approximate LHC bound of mg 900 GeV (for mq ≫ mg) [82,83], while the lower bound on m A comes from LHC searches for A, H → τ + τ − which require m A > 500 GeV at tan β ∼ 50 [84]. We require of our solutions that 1. electroweak symmetry be radiatively broken (REWSB), 2. the neutralino Z 1 is the lightest MSSM particle, 3. the light chargino mass obeys the rather model independent LEP2 limit that m W 1 > 103.5 GeV [85], and 4. the light Higgs mass falls within the window m h = 122-128 GeV, where we adopt ±3 GeV as theoretical error on the Higgs mass calculation.
Regarding B-physics constraints, we consider BR(B → X s γ) = (3.55 ± 0.34) × 10 −4 , where experimental and SM theoretical uncertainties [86] have been added in quadrature. We combine the above uncertainty with a ∼ 5% uncertainty from new physics (SUSY) contributions, 3 and then impose the BR(B → X s γ) constraint at the 2σ level. For the BR(B s → µ + µ − ) constraint, we assume a theoretical uncertainty of 10% [87]. We combine this uncertainty linearly with the 95% CL limit BR(B s → µ + µ − ) < 4.2 × 10 −9 . This leads to the the following limits which we will use throughout the numerical analysis.
Regarding neutralino relic density, we remark here that models of natural SUSY contain a higgsino-like lightest neutralino with thermal abundance of typically Ω TP  3) The light blue triangles obey the BR(B s → µ + µ − ) constraint, but deviate from the measured BR(B → X s γ) by more than 2σ. Finally, the pink squares satisfy both the BR(B s → µ + µ − ) and BR(B → X s γ) constraints. Only points with ∆ < 1000 and R yuk < 1.5 are shown.
of NS models in the sense that if one invokes the axion solution to the strong CP problem, then one expects mixed axion-higgsino dark matter, where the higgsino portion is typically enhanced by thermal axino production and decay to higgsinos in the early universe. Thus, a thermal under-abundance leaves room for additional non-thermal higgsino production plus an axion component to the dark matter [88].

Scan results
As our first result, we show in figure 1 points from our parameter space scan in the R yuk vs. ∆ plane. All points have m h = 122-128 GeV and obey the current LEP and LHC SUSY mass limits. The dark blue triangles however violate the B-physics constraints of eq. (2.3). The light blue triangles obey the BR(B s → µ + µ − ) constraint, but deviate from the measured BR(B → X s γ) by more than 2σ. Finally, the pink squares satisfy both the BR(B s → µ + µ − ) and BR(B → X s γ) constraints. This color scheme is used throughout the remainder of the paper. We see already from this plot that flavor physics constraints significantly affect the parameter space of Yukawa-unified natural SUSY. This is to be expected, since naturalness requires lighter third generation squarks while Yukawa unification requires tan β ∼ 50: both these effects bolster SUSY contributions to B-physics observables. (Analogous observations were made in [38,39,41,46,89] in the context of generic YU models.) Without flavor constrains, we can obtain R yuk as low as ∼ 1.2. The BR(B s → µ + µ − ) constraint pushes this up to R yuk 1.27, and the BR(B → X s γ) constraint to R yuk 1.3. Aside from the B-physics constraints, it is intriguing that points with lowest ∆ also have lowest values of R yuk . This is because low ∆ requires light third generation squark masses, while at the same time Yukawa unification requires large SUSY threshold corrections which also require lighter third generation squarks.  In figure 2, we show the value of R yuk versus various input parameters for ∆ < 100 (less than 1% fine-tuning). Here, the behavior deviates considerably from t−b−τ unified models with large µ > 0 and a spectrum derived from the radiatively driven inverted scalar mass hierarchy [40,41,46]. From frame a), we see that YUNS models actually prefer large m 1/2 whereas generic Yukawa unified SUSY models (with R yuk 1.1 but arbitrary EWFT), YUS for short, prefer low m 1/2 . The preference for large m 1/2 helps to avoid LHC constraints on the gluino mass. If we impose the B-physics constraints, then the distribution flattens out with some preference for lower m 1/2 values. In frame b), where R yuk is plotted vs.  In figure 4 we show R yuk versus various sparticle and Higgs masses. In frame a), the distribution versus m h is shown, and we see that low R yuk prefers the lower range of m h . This is understandable since low EWFT prefers lower top and bottom squark masses, which may not feed a sufficient radiative correction into the m h computation. These low third generation squark mass solutions also tend to give large contributions to B-constraint. If we impose B-constraints, then m h can live in the 122 − 128 GeV range, at the cost of larger R yuk . In frame b), we show the distribution versus mg. While all solutions -especially JHEP12(2012)066  low R yuk ones -favor the heavier range of mg (mg ∼ 2 − 5 TeV, likely beyond LHC reach), the solutions obeying B-constraints tend to slightly favor lower mg, possibly within range of LHC14 searches. In frame c), the distribution versus mt 1 is shown. Here, we see a clear demarcation: B-constraints favor mt 1 1.5 TeV to suppress SUSY loop contributions to BR(B → X s γ). This constraint forces the minimum R yuk to move from ∼ 1.18 to about 1.3. (Note that additional small flavor-violating contributions to the MSSM Lagrangian could alter the predicted BR(B → X s γ) and/or BR(B s → µ + µ − ) rates [101,102].) In frame d ), we show the distribution versus pseudoscalar. Higgs mass m A . We see low R yuk favors the lower range of m A , although values up to and beyond 5 TeV are also possible (at the cost of R yuk ∼ 1.4, however). Since tan β ∼ 50, LHC searches for A, H → τ + τ − will access a significant range of m A in this case. It is also possible for some range of m A 1 TeV for LHC to access bA → bµ + µ − production [105].
In figure 5, we show scatter plots of YUNS points with R yuk < 1.4 and ∆ < 100 in a) mt 1 vs. mb 1 and b) mt 2 vs. mb 2 space. From frame a), we see as expected that mt 1 and mb 1 are correlated, with some solutions reaching well below mt 1 ∼ 500 GeV, which should be accessible to LHC searches. However, in this case these points all violate B-constraints, so  that requiring B-constraints within measured range requires instead mt 1 , mb 1 1.5 TeV, likely beyond the 14 TeV LHC reach. In frame b), we see also that mb 2 and mt 2 are correlated. This is different from usual NS, where mb 2 can be far above mt 1,2 and mb 1 . The reason here is that f b is large and so there is also a non-negligible contribution to Σ u fromb 1,2 . Imposing B-constraints, we find mt 2 and mb 2 both 2 TeV (the former aids in lifting m h into its measured range).

Benchmark points
In this section we present several YUNS benchmark points, see table 1, and compare to one YUS benchmark point (HSb, from the "just-so" HS model) from ref. [90].
For HSb, R yuk ∼ 1.02, which is nearly perfect Yukawa coupling unification. The Higgs mass m h ≃ 127.8 GeV is also sufficiently heavy. Unfortunately, the point is now excluded by LHC7 searches for multi-jet+E miss T plus one b-tag searches, since mg is only 351 GeV. We also list in table 1 the EWFT measure; for HSb we have ∆ = 2489, indicating exceptionally high level of fine-tuning. Much of this comes from the µ parameter which turns out to be of order 3 TeV, and thus requires a large value of m 2 Hu at the weak scale to cancel against. In contrast, point YUNS1 in column 3 has low fine-tuning of ∆ = 39, at the cost of relaxing R yuk to 1.37. For YUNS1, mg ≃ 1.3 TeV, with first/second generation squarks at ∼ 19 TeV. The point is likely beyond LHC8 reach, but should be accessible to LHC14 with 10-100 fb −1 . The point is slightly in violation of our listed condition on BR(b → sγ). However, by lowering m 0 (1, 2) slightly, then mt 1 is increased and also weak scale flavor violation is decreased [103,104] (the GUT scale squark mass matrices become more nearly proportional to the identity matrix). For instance, lowering m0(1, 2) to 15 TeV increases BR(b → sγ) to 3.1 × 10 −4 , while R increases to 1.41 and ∆ increases to 54. In column 4, we list YUNS2 with R yuk = 1.3, as low as allowed by B-constraints, but with ∆ ∼ 100. This point has mg ∼ 4.3 TeV and mq ∼ 17.5 TeV, so it is likely beyond LHC reach, including a high-luminosity upgrade. While the higgsino-like chargino is only 112.5 GeV, it decays via 3-body mode into a higgsino-like Z 1 with m Z 1 = 106.   Table 1. Parameters and masses in GeV units for HSb [90] and three Yukawa-unified natural SUSY (YUNS) benchmark points. We also show B-decay constraints and dark matter relic density and (in)direct detection cross sections.
visible decay products are very soft, and likely impossible to observe above SM backgrounds at the LHC. Point YUNS3, listed in column 5, features R yuk as low as 1. are beyond LHC reach, the rather light top and bottom squarks may be accessible to LHC searches. All these points have Ω Z 1 h 2 ≪ 0.11, leaving room for non-thermal higgsino production and axions. This contrasts point HSb, which has a much too thermal abundance and so would need an extremely light axino or huge late-time entropy production to tame this over-abundance [88].
In figure 6, we show the Yukawa coupling evolution of f t , f b and f τ versus renormalization group scale Q, from m weak to M GUT for benchmark points YUNS1 and YUNS3. These can be compared to similar plots for YUS, as in e.g. figure 6 of ref. [91]. The SUSY threshold corrections implemented at the scale Q = √ mt 1 mt 2 show up as jumps in the curves. In the case of YUS models, the m b threshold correction is positive due to a largẽ t i W j loop contribution which goes like δ b ∼ (f 2 t /32π 2 )(µA t /m 2 t ) tan β, where both µ and A t are extremely large. For YUNS models, with rather low µ, these loops are suppressed and in the case of YUNS1, thegb i loops actually dominate, and are of opposite sign to thẽ t i W j loops, leading to the slight downward jump of f b and thus bad Yukawa coupling unification. For YUNS3, thet i W j loops are larger, and the jump goes upwards, thus providing better Yukawa coupling unification.

Yukawa-unified natural SUSY: LHC, ILC and DM searches
In this section, we discuss the observable consequences of YUNS for LHC, ILC and direct and indirect WIMP and also axion searches.

YUNS at LHC
In natural SUSY models, it is favorable to have multi-TeV first/second generation squarks and sleptons because 1. they are safely beyond current LHC searches, 2. they provide at least a partial solution to the SUSY flavor and CP problems and 3. they provide additional suppression of third generation scalar masses via large 2-loop RGE effects. However, this means they are likely beyond any conceivable LHC reach. Third generation squarks may JHEP12(2012)066 be much lighter, and in generic NS models are naively expected to be below the TeV scale (but see ref. [62] where 1-4 TeV third generation squarks work just fine, and lift the value of m h into its measured range). In the case of YUNS, with tan β ∼ 50, the combined light squarks and large tan β usually imply violation of B-constraints, and if these are imposed, then top and bottom squarks are beyond 1.5 TeV, and likely inaccessible to LHC searches.
For YUNS, the gluino mass may lie anywhere in the 1-5 TeV range. It has been estimated in ref. [92][93][94][95][96][97][98][99] that LHC14 with 100 fb −1 should be able to access gluino pair production in the case of heavy squarks for mg up to 1.8 TeV. In the case of YUNS with heavier top and bottom squarks, theg is expected to dominantly decay via 3-body modes into tt Z i and tb W i final states. The gluino pair events will thus contain multi-jets plus missing energy plus isolated leptons plus several identifiable b-jets [100]. While the higgsinolike chargino and neutralino production cross sections can be large, their decays to soft visible particles, arising from the small energy release in their 3-body decays, will be difficult to detect at LHC above SM backgrounds [78].
In the case where m A 1 TeV, then it may be possible to detect A, H → τ + τ − , especially if these are produced in association with b-jets, e.g. pp → bA, bH production. It may also be possible to detect bA, bH production with A, H → µ + µ − [105], since the production and decay are enhanced at large tan β. In this case, the A, H mass and width may be determined by reconstructing m(µ + µ − ). At large tan β, this width is typically in the tends of GeV range and is very sensitive to tan β. This reaction offers a method to distinguish YUNS from NS, in that the former is expected to occur at tan β ∼ 50.

YUNS at ILC
A linear e + e − collider operating at √ s ∼ 250 − 500 GeV would in many ways be an optimal discovery machine for YUNS. The reason is that by construction µ 250 GeV, so chargino and neutralino pair production should always be available. While the small energy release in W 1 and Z 2 decay is problematic at LHC, it should be much more easily observable in the clean environment of an e + e − collider. In this sense, an e + e − collider operating at √ s ∼ 250 − 500 GeV would be a higgsino in addition to a Higgs factory. It is also possible that some lighter third generation squarks are accessible to ILC with √ s ∼ 1 TeV or CLIC with √ s = 3 TeV, depending if one avoids B-constraints and accepts the low mass, low R yuk solutions.

Higgsino-like WIMPs
A generic prediction of both NS and YUNS models is that the LSP is a higgsino-like WIMP with a typical under-abundance of thermally produced (TP) neutralinos Ω TP Z 1 h 2 ∼ 0.002−0.01. However, in cases where m 1/2 is as low as ∼ 300 GeV and µ is as large as 200− 250 GeV, then there can be substantial bino-higgsino mixing, boosting Ω TP which contains a pseudoscalar axion a as well as a spin-1/2 axinoã and a spin-0 saxion s. In such models, one expects both saxion and axino masses at or around the SUSY breaking scale, so that dark matter is comprised of an axion-WIMP admixture. In this case, thermal production of axinos and subsequent decay to states such agg bolster the WIMP abundance beyond its TP-value. In addition, axions are produced as usual via coherent oscillations. It is also possible to suppress the WIMP abundance in the mixed a Z 1 cosmology via latetime entropy production from saxion production and decay, although this case tends to be highly constrained by maintaining successful Big Bang nucleosynthesis (BBN). The upshot is that in the mixed a Z 1 dark matter scenario, it may be possible to detect both a WIMP and an axion.
In the case of the YUNS model, the higgsino-like neutralinos have a substantial spinindependent (SI) direct detection cross section, as illustrated in table 1, where σ SI ( Z 1 p) ∼ 10 −8 pb. While this level of direct detection cross section is now highly constrained by recent XENON100 [106] results, one must bear in mind that the higgsino-like WIMPs would constitute only a portion of the dark matter, so their local abundance might be up to a factor of 30 lower than is commonly assumed. In table 1, we also list σ SD ( Z 1 p) (relevant for WIMP detection at IceCube) and σv | v→0 , relevant for detection of dark matter annihilation into gamma rays or anti-matter throughout the cosmos. While these cross sections are also at potentially observable levels, again one must take into account that the overall WIMP abundance may be up to a factor of about 30 below what is commonly assumed. Plots of σ SI ( Z 1 p), σ SD ( Z 1 p) and σv | v→0 have been presented in the case of higgsino-like WIMPs in refs. [78] and [81] and so similar plots will not be reproduced here.

Summary and conclusions
Previous analyses of t − b − τ Yukawa-unified models suffer from two problems: 1. they tend to predict a light gluino mg 500 GeV (alhough solutions are possible for much heavier gluinos) which is now excluded by LHC searches, and 2. they suffer from extreme JHEP12(2012)066 fine-tuning in the electroweak sector. In this paper, we examined how well the Yukawa couplings could unify in the natural SUSY context, where µ ∼ 100 − 250 GeV, while at the same time requiring the light Higgs mass m h ∼ 122 − 128 GeV. The small value of µ suppresses the larget i W j loop contributions to the b-quark Yukawa coupling which seem to be needed for precision Yukawa coupling unification. Nonetheless, by scanning over NUHM2 parameters with split third generation, we are able to find solutions with R yuk as low as 1.18. These solutions, with very light third generations squarks and tan β ∼ 50 tend to violate B-physics constraints. If B-physics constraints are imposed, then only R yuk ∼ 1.3 can be achieved.
The Yukawa-unified natural SUSY spectra have important differences from previous YU spectra. The gluino mass can easily be in the 1-4 TeV range, thus avoiding LHC constraints from SUSY searches. The light higgsino-like charginos and neutralinos decay to soft particles, also avoiding LHC searches. However, light higgsinos should be easily accessible to an ILC with √ s ∼ 0.25 − 05 TeV, as is typical of all NS models. In addition, as in all NS models, the lightest neutralino is higgsino-like with a typical thermal underabundance of WIMP dark matter. We regard this as a positive feature in that the WIMP abundance is typically increased in non-standard (but more attractive) cosmologies such as those conatining mixed axion-neutralino cold dark matter.
The question arises as to how to distinguish YUNS from ordinary NS. The YUNS model requires tan β ∼ 50, which leads to large production cross sections for heavy Higgs bosons A and H at LHC, and large widths for these particles. If the rare decays A, H → µ + µ − can be identified with suficiently high statistics (perhaps at a luminosity upgraded LHC), then the widths may be measured with precision, allowing one to highly constrain tan β, and perhaps verify that it is consistent with YUNS models.
Note added. After completion of this study, a first measurement of BR(B s → µ + µ − ) = 3.2 +1.5 −1.2 × 10 −9 became available [107]. Requiring consistency with this measurement at the 2σ level actually gives a less severe constraint than the upper bound in eq. (2.3). Our conclusions regarding the minimum achievable R yuk remain unchanged.