Muon $g-2$ in Split-Family SUSY in light of LHC Run II

TL;DR

This paper reexamines the Split-Family SUSY model, which explains muon g-2 deviation and Higgs mass with hierarchical sfermion masses, considering recent LHC constraints and flavor mixing effects.

Contribution

It provides an updated analysis of Split-Family SUSY in light of LHC data and explores flavor mixing constraints within the model.

Findings

The model remains viable under current LHC constraints.

Flavor changing neutral currents are not severely constrained with CKM as the sole source of mixing.

Large flavor mixing in SUSY breaking parameters is possible.

Abstract

The Split-Family supersymmetry is a model in which the sfermion masses of the first two generations are in $\mathcal{O}(100\text{-}1000)$ GeV while that of the third one is in $\mathcal{O}(10)$TeV. With such a hierarchical spectrum, the deviation of the muon $g-2$ and the observed Higgs boson mass are explained simultaneously. In this paper, we revisit the Split-Family SUSY model in light of the updated LHC constraints. We also study the flavor changing neutral current problems in the model. As we will show, the problems do not lead to stringent constraints when the Cabibbo-Kobayashi-Maskawa matrix is the only source of the flavor mixing. We also study how large flavor mixing in the supersymmetry breaking parameters is allowed.

Tables

Table 1. Table 1: A sample point in the universal gaugino mass case. m gluino subscript 𝑚 gluino m_{\rm gluino} , m Q ~ subscript 𝑚 ~ 𝑄 m_{\tilde{Q}} , m e ~ L ​ ( m μ ~ L ) subscript 𝑚 subscript ~ 𝑒 𝐿 subscript 𝑚 subscript ~ 𝜇 𝐿 m_{\tilde{e}_{L}}~{}(m_{\tilde{\mu}_{L}}) , m e ~ R ​ ( m μ ~ R ) subscript 𝑚 subscript ~ 𝑒 𝑅 subscript 𝑚 subscript ~ 𝜇 𝑅 m_{\tilde{e}_{R}}~{}(m_{\tilde{\mu}_{R}}) , m χ ~ 0 1 subscript 𝑚 subscript superscript ~ 𝜒 1 0 m_{\tilde{\chi}^{1}_{0}} , m χ ~ 0 2 subscript 𝑚 subscript superscript ~ 𝜒 2 0 m_{\tilde{\chi}^{2}_{0}} , m χ ~ 1 ± subscript 𝑚 superscript subscript ~ 𝜒 1 plus-or-minus m_{\tilde{\chi}_{1}^{\pm}} denote the masses of the gluino, the lightest squark, the lightest almost left handed selectron (smuon), the lightest almost right-handed selectron (smuon), the lightest neutralino, the next to the lightest neutralino, the lightest chargino, respectively. It should be noted that R 𝑅 R -parity violation is required for the decay of the LSP before the Big-Bang Nucleosynthesis (BBN). Then, we also need to introduce a dark matter candidate in the model. (See the next subsection.)

$m_0,m_3$	$0\mathrm{GeV},12\mathrm{TeV}$
$M_{1/2}$	$1650\mathrm{GeV}$
$\tan \beta$	$50$
$m_{\mathrm{Higgs}}$	$124.6\mathrm{GeV}$
$a_{\mu }$	$1.48\times {10}^{-9}$
LSP	charged slepton
$\Delta M_K$	$-3.0\times {10}^{-21}\mathrm{GeV}$
${ϵ}_K$	$-2.6\times {10}^{-7}$
$\Delta M_D$	$1.5\times {10}^{-17}\mathrm{GeV}$
$\mu$	$9751\mathrm{GeV}$
$m_{\mathrm{gluino}}$	$3619\mathrm{GeV}$
$m_{\tilde{Q}}$	$2730\mathrm{GeV}$
$m_{{\tilde{e}}_L}\left(m_{{\tilde{\mu }}_L}\right)$	$486\mathrm{GeV}$
$m_{{\tilde{e}}_R}\left(m_{{\tilde{\mu }}_R}\right)$	$942\mathrm{GeV}$
$m_{{\chi }_0^1}$	$766\mathrm{GeV}$
$m_{{\chi }_0^2}$	$1360\mathrm{GeV}$
$m_{{\chi }_1^{\pm }}$	$1360\mathrm{GeV}$

Table 2. Table 2: A sample parameter point in the case of the non-universal gaugino masses. The neutralino is the LSP. Ω ​ h 2 Ω superscript ℎ 2 \Omega h^{2} and σ SI superscript 𝜎 SI \sigma^{\rm SI} denote the current thermal relic abundance of the bino and the spin-independent bino-nucleon cross-section.

$m_0,m_3$	$650\mathrm{GeV},11\mathrm{TeV}$
$M_1,M_2,M_3$	$1.725\times 480\mathrm{GeV},480\mathrm{GeV},-2.6\mathrm{TeV}$
$\tan \beta$	$40$
$m_{\mathrm{Higgs}}$	$124.6\mathrm{GeV}$
$a_{\mu }$	$2.69\times {10}^{-9}$
LSP	bino
$\Omega h^2$	$0.119$
${\sigma }^{\mathrm{SI}}$	$2.4\times {10}^{-14}\mathrm{pb}$
$\Delta M_K$	$-2.0\times {10}^{-21}\mathrm{GeV}$
${ϵ}_K$	$1.2\times {10}^{-7}$
$\Delta M_D$	$1.5\times {10}^{-17}\mathrm{GeV}$
$\mu$	$9108\mathrm{GeV}$
$m_{\mathrm{gluino}}$	$5529\mathrm{GeV}$
$m_{\tilde{Q}}$	$4442\mathrm{GeV}$
$m_{{\tilde{e}}_L}\left(m_{{\tilde{\mu }}_L}\right)$	$481\mathrm{GeV}$
$m_{{\tilde{e}}_R}\left(m_{{\tilde{\mu }}_R}\right)$	$635\mathrm{GeV}$
$m_{{\chi }_0^1}$	$404\mathrm{GeV}$
$m_{{\chi }_0^2}$	$478\mathrm{GeV}$
$m_{{\chi }_1^{\pm }}$	$478\mathrm{GeV}$

Equations

$$\text{mit}\Delta a_{\mu }\equiv a_{\mu }^{\exp }-a_{\mu }^{\mathrm{SM}}=(27.06\pm 7.26)\times 10^{-10}.$$

$$\displaystyle m_{\rm soft}^{2}=\left(\begin{array}[]{ccc}m_{0}^{2}&0&0\\ 0&m_{0}^{2}&0\\ 0&0&m_{3}^{2}\end{array}\right),$$

$$W_{\mathrm{LEPTON}}$$

$$W_{\mathrm{QUARK}}$$

$$W_{\mathrm{LEPTON}}={\hat{f}}_E^{ii}{L}_i^{[e]}{\bar{E}}_i^{[e]}{H}_d,$$

$$W={\hat{f}}_U^{ii}{Q}_i^{[d]}{\bar{U}}_i^{[u]}{H}_u+(V_{\mathrm{CKM}}^{\ast }{\hat{f}}_D{)}^{ij}{Q}_i^{[d]}{\bar{D}}_j^{[d]}{H}_d.$$

$$W=(V_{\mathrm{CKM}}^T{\hat{f}}_U{)}^{ij}{Q}_i{\bar{U}}_j{H}_u+({\hat{f}}_D{)}^{ii}{Q}_i{\bar{D}}_i{H}_d,$$

$$\displaystyle m_{\rm soft}^{2}=V^{\dagger}_{\rm mix}\left(\begin{array}[]{ccc}m_{0}^{2}&0&0\\ 0&m_{0}^{2}&0\\ 0&0&m_{3}^{2}\end{array}\right)V_{\rm mix}.$$

$$m_0,m_3,m_{H_u}^2,m_{H_d}^2,\tan \beta ,M_1,M_2,M_3.$$

$$m_{H_u}^2=m_{H_d}^2=m_0^2,$$

$$|\mu {|}^2$$

$$M_{1/2}=M_1=M_2=M_3=\mathcal{O}(100-1000)\mathrm{GeV}.$$

$$M_1\simeq M_2=\mathcal{O}(100-1000)\mathrm{GeV},M_3\gtrsim -\mathcal{O}(1)\mathrm{TeV}.$$

$$\tan \beta \gtrsim 40.$$

$$M_1(M_{\mathrm{GUT}})\simeq 1.7\times M_2(M_{\mathrm{GUT}}),$$

$${\sigma }_{\mathrm{SI}}\lesssim 10^{-10}\left(\frac{m_{\mathrm{DM}}}{100\mathrm{GeV}}\right)\mathrm{pb}$$

$$|{ϵ}_K{|}^{\mathrm{ex}}=2.228(11)\times 10^{-3}(90\%CL).$$

$$|{ϵ}_K{|}^{\mathrm{SM},\mathrm{in}}=2.05(18)\times 10^{-3},$$

$$\text{mit}\Delta M_K^{\exp }=3.483(6)\times 10^{-15}\mathrm{GeV}.$$

$$|\text{mit}\Delta M_D^{\mathrm{ex}}|=0.63_{-0.29}^{+0.27}\times 10^{-14}\mathrm{GeV}.$$

$$\mathcal{B}({\mu }^{+}\to e^{+}\gamma )<4.2\times 10^{-13}(90\%\text{CL}).$$

$$\displaystyle{\small V_{\rm mix}=\left(\begin{array}[]{ccc}1&0&0\\ 0&{\rm cos}\,\theta_{23}&{\rm sin}\,\theta_{23}\\ 0&-{\rm sin}\,\theta_{23}&{\rm cos}\,\theta_{23}\end{array}\right)\left(\begin{array}[]{ccc}{\rm cos}\,\theta_{13}&0&{\rm sin}\,\theta_{13}\\ 0&1&0\\ -{\rm sin}\,\theta_{13}&0&{\rm cos}\,\theta_{13}\end{array}\right)\left(\begin{array}[]{ccc}{\rm cos}\,\theta_{12}&{\rm sin}\,\theta_{12}&0\\ -{\rm sin}\,\theta_{12}&{\rm cos}\,\theta_{12}&0\\ 0&0&1\end{array}\right).}$$

$$(\sin {\theta }_{12},\sin {\theta }_{13},\sin {\theta }_{23})=(\mathcal{O}(ϵ),\mathcal{O}({ϵ}^3),\mathcal{O}({ϵ}^2)),$$

$$(s_{12},s_{13},s_{23})=(R_1ϵ,R_2{ϵ}^3,R_3{ϵ}^2),$$

$$\frac{\text{mit}\Delta y_b}{y_b}\propto \frac{M_3\mu \tan \beta }{m_3^2},$$

$${ϵ}_{K_{\mathrm{SUSY}}}$$

$${\text{mit}\Delta }_{LL}^{ds}$$

$$\displaystyle m_{\rm soft}^{2}=\left(\begin{array}[]{ccc}m_{0}^{2}&0&0\\ 0&m_{2}^{2}&0\\ 0&0&m_{3}^{2}\end{array}\right),$$

$$W={\hat{f}}_E^{ij}{L}_i^{[e]}{\bar{E}}_j^{[e]}{H}_d+(U_{\mathrm{MNS}}^{\ast }{\hat{f}}_{\nu }{)}^{ij}{L}_i^{[e]}{N}_{Rj}{H}_u+\frac{1}{2}{M}_R{N}_{Ri}{N}_{Ri}.$$

$$({\text{mit}\Delta }_{12}^e{)}_{LL}\simeq \frac{1}{16{\pi }^2}{\left(\sqrt{\frac{M_R}{10^{15}\mathrm{GeV}}}\right)}^2\left(\frac{m_3^2}{m_0^2}\right)(U_{\mathrm{MNS}}{)}_{13}(U_{\mathrm{MNS}}^{\ast }{)}_{23}\ln (M_R/10^{16}\mathrm{GeV}),$$

$$\mathcal{B}\simeq \frac{48{\pi }^3{\alpha }_{\mathrm{em}}}{G_F^2}{\left(c_1\frac{{\alpha }_1}{4\pi }\frac{\mu m_{\mathrm{bino}}\tan \beta }{m_{\mathrm{slepton}}^4}({\text{mit}\Delta }_{12}^e{)}_{LL}\right)}^2.$$

Full text

††preprint: IPMU19-0044††preprint: TU-1076

I Introduction

The standard model (SM) of particle physics is now complete by the discovery of the Higgs boson with a mass around $125$ GeV Aad et al. (2012); Chatrchyan et al. (2012); Aad et al. (2015, 2016); collaboration (2018). In the minimal supersymmetric SM (MSSM) (see Nilles (1984) and references therein), the measured Higgs boson mass can be explained when the masses of the superpartners of the top quark (the stop) are in $\mathcal{O}(10\text{--}100)$ TeV Okada et al. (1991a, b); Ellis et al. (1991a); Haber and Hempfling (1991); Ellis et al. (1991b). No evidence of the super particles at the LHC experiments also suggests that their masses are in a multi-TeV range.

On the other hand, the measured value of the muon anomalous magnetic moment $a_{\mu}^{\rm exp}$ Hagiwara et al. (2011); Keshavarzi et al. (2018) deviates from the SM prediction $a_{\mu}^{\rm SM}$ Bennett et al. (2006); Roberts (2010) at about $3.7\sigma$,111Davier $et.al.$ have reported a deviation of $3.6~{}\sigma$ level Davier et al. (2011). The deviation has increased Jegerlehner (2018) by the recent more precise measurement of the fine-structure constant from atomic interferometry with Cesium133 et al. . See Lusiani (2018) for the future prospect of the muon $g-2$ measurement.

[TABLE]

This discrepancy can be explained in the MSSM when such as the smuons and the electroweakinos are in $\mathcal{O}(100)$ GeV Lopez et al. (1994); Chattopadhyay and Nath (1996); Moroi (1996).

In Ibe et al. (2013a), the Split-Family supersymmetry (SUSY) model has been proposed to explain the observed Higgs boson and the muon $g-2$ deviation simultaneously. There, the sfermions of the first two generations are in $\mathcal{O}(100\text{--}1000)$ GeV while that of the third one is in $\mathcal{O}(10)$ TeV.222For other simple possibilities and models in the MSSM, see Ref. $e.g.$  Evans et al. (2011, 2012); Ibe et al. (2013b); Sato et al. (2012); Akula and Nath (2013); Bhattacharyya et al. (2013, 2014); Chakrabortty et al. (2014); Evans et al. (2014); Mohanty et al. (2013); Babu et al. (2014); Gogoladze et al. (2014); Iwamoto et al. (2015); Adeel Ajaib et al. (2015); Chowdhury and Yokozaki (2015); Gogoladze et al. (2015); Harigaya et al. (2015a, b); Gogoladze and Un (2017); Yin and Yokozaki (2016); Yanagida and Yokozaki (2017); Wang et al. (2018). Such a hierarchical SUSY spectrum is motivated by the Yukawa hierarchy.

In this paper, we revisit the Split-Family SUSY model in light of the updated LHC constraints.333For model independent study, See Ajaib et al. (2015); Tran and Nguyen (2019). As we will see, almost the entire region which explains the muon $g-2$ within $2\sigma$ is excluded for the universal gaugino mass. We also show that the collider constraints can be evaded for the non-universal gauino masses while explaining the muon $g-2$.

We also study the FCNC problems in the Split-Family SUSY model. As we will explain, the precise construction of the model requires careful treatment of the family basis, which generically leads to sizable SUSY contributions to the flavor changing neutral currents (FCNC). To see such effects, we first discuss the case where the Cabibbo-Kobayashi-Maskawa (CKM) matrix is the only source of the flavor mixing. As we will see, the FCNC constraints are not so stringent in that case. We also discuss how large flavor mixing in the supersymmetry breaking parameters are allowed.

The organization of this paper is as follows. In Sec. II, we review the Split-Family SUSY model. In Sec. III, we update the constraints from the collider experiments. In Sec. IV, we discuss the FCNC constraints on the model. In Sec. V, we discuss the successful bottom-tau unification as a bonus feature of the model. The final section is devoted to our conclusion.

II Split-Family SUSY Model

The basic idea of the Split-Family model is to assume the hierarchical soft masses,

[TABLE]

where the masses of the first two generation sfermions, $m_{0}^{2}$, are in $m_{0}=\mathcal{O}(100-1000)$ GeV and that of the third generation sfermions, $m_{3}^{2}$, is in $m_{3}=\mathcal{O}(10)$ TeV. This structure may be related to the hierarchy of the Yukawa couplings. For example, an extra-dimensional setup can lead to the aligned hierarchy between the soft masses and the Yukawa couplings, where only the third generation resides on the brane close to those of the Higgs brane and the SUSY breaking brane (See $e.g.$ Gabella et al. (2007)). It is also possible to achieve the aligned hierarchy if the first two generations are pseudo-Nambu-Goldstone multiplets of some broken global symmetry (See $e.g.$ Kugo and Yanagida (1984); Yanagida and Yasui (1986); Evans et al. (2014)).

The precise construction of the Split-Family model requires careful treatment of the family basis. If the soft masses are universal, for example, the general Yukawa couplings in the superpotential do not lead to the SUSY FCNC contributions, since the soft breaking parameters are proportional to the unit matrix in any family basis. On the other hand, for the non-universal soft masses, the general Yukawa couplings result in the non-zero SUSY contributions to the flavor mixing. Since we have assumed the rather light sfermions, those mixing leads to unacceptably large FCNC processes.

In this paper, we put a phenomenological requirement that the soft masses in Eq. (5) are realized in the family basis where the Yukawa couplings in the superpotential are almost diagonal.444We also assume that the hierarchy of the soft masses and the Yukawa couplings are aligned. These assumptions are implicitly made, for example, in the models Gabella et al. (2007); Kugo and Yanagida (1984); Yanagida and Yasui (1986); Ibe et al. (2013a) to avoid too large FCNC. It should be noted, however, that there are unavoidable SUSY FCNC contributions even under these assumptions due to the effects of the CKM matrix. To demonstrate those effects, we study the following two scenarios.

• •

The CKM matrix is the only source for the flavor mixing (the minimal mixing scenario)

• •

Small flavor mixing comes from the supersymmetry breaking parameters (the small mixing scenario)

In the following, we call the first scenario the minimal mixing scenario and the second one the small mixing scenario. The minimal mixing scenario gives us a demonstration of how large SUSY FCNC contributions are expected by the effect of the non-universality of the soft masses in Eq. (5).

Minimal Mixing Scenario

Let us start from the minimal mixing scenario. The superpotential with general Yukawa matrices is given by,

[TABLE]

Here, $L,~{}\bar{E},~{}Q,~{}\bar{U},~{}\bar{D}$ are chiral superfields of the left-handed leptons, the left-handed anti- electrons, the left-handed quarks, the left-handed up-type anti-quarks, and the left-handed down type anti-quarks, respectively. The coefficients $f_{E},~{}f_{U},~{}f_{D}$ are the $3\times 3$ complex matrices. The subscripts $i,~{}j$ run from one to three and each one is contracted.

As we assume that the soft masses in Eq. (5) are defined in the family basis where the Yukawa couplings are diagonal, the lepton Yukawa matrix is given by the diagonal form,

[TABLE]

where the subscript $[e]$ shows the chiral superfields in the diagonalized basis, $\hat{f}_{E}$ is a diagonal and real-positive matrix.555Even if we include the PMNS effect, the lepton flavor mixing is not so large for the model with the degenerate right-handed masses $M_{R}\lesssim 10^{10}$ GeV (See the appendix C for more details). Thus, we ignore the finite neutrino masses in the superpotential. Note that the subscript $i$ in Eq. (8) runs from one to three, which corresponds to the charged lepton generation.

On the other hand, the up- and down-type quark Yukawa matrices are not diagonalized at the same time due to the CKM matrix. To demonstrate the minimal mixing scenario, we take a simple family basis, for example,

[TABLE]

Here, $\hat{f}_{U}$ and $\hat{f}_{D}$ are diagonal and real-positive matrices, $V_{CKM}$ is the CKM matrix, the superscript $*$ is a complex conjugate of a matrix, and the superscript $T$ denotes the transpose of a matrix. The subscripts $i,~{}j$ in Eq. (9), run from one to three, which corresponds to the quark generation. In this basis, the down-type Yukawa couplings are not diagonal due to the CKM matrix. We may consider another simple family basis where the CKM matrix appears in the up-type Yukawa couplings666That is, we may take the basis,

$\displaystyle W=(V_{\rm CKM}^{T}\hat{f}_{U})^{ij}Q_{i}{\bar{U}}_{j}H_{u}+(\hat{f}_{D})^{ii}Q_{i}{\bar{D}}_{i}H_{d},$

(10) while the soft masses are given by Eq. (5). (See Sec. IV for more details).

Small Mixing Scenario

The small mixing scenario is defined in the following way. We first take the family basis specified in Eq. (8) and Eq. (9). Then, we introduce the small flavor mixing matrix $V_{\rm mix}$ so that the mass matrices are given by

[TABLE]

In the small mixing scenario, we assume that $V_{\rm mix}$ is close to the unit matrix with the mixing angles of the order of the CKM angles. See Sec. IV for more details.

III Phenomenology of Split-Family SUSY model

In this section, we update the favored region for the muon $g-2$ and the Higgs boson mass in Ref. Ibe et al. (2013a) in light of the current LHC data. The SUSY contributions to the FCNC are discussed in the next section.

III.1 Parameter Choice At Input Scale

The free parameters in our analysis are,

[TABLE]

Here, $m^{2}_{H_{u}}~{}(m_{H_{d}}^{2})$ is the squared mass of the up-type (down-type) Higgs doublet, ${\rm tan}\beta$ is defined by the ratio of the vacuum expectation value of $H_{u}$ to that of $H_{d}$, and $M_{1},~{}M_{2}$, and $M_{3}$ are the gaugino mass parameters of the bino, the wino, and the gluino, respectively. We assume the same $m_{0}^{2}$ and $m_{3}^{2}$ for the squarks and the sleptons for simplicity.777These assumptions are motivated by the $SU(5)$ GUT.

We take the above free parameters at the GUT scale $M_{\rm GUT}=10^{16}\,{\rm GeV}$. By solving the renormalization group equations, we obtain the physical parameters. Notice that the Higgsino mass parameter $\mu$ is chosen to achieve electroweak symmetry breaking consistently (See Eq. (17)). It should be also noted that the following arguments are not sensitive to the choice of $m_{H_{u}}^{2}$ and $m_{H_{d}}^{2}$ as long as the electroweak symmetry breaking is successful.888When the required value of $\mu$ is as small as $m_{0}$, the Higgsino also contributes to the muon $g-2$ Ibe et al. (2013a). For simplicity, we take

[TABLE]

in the following analysis.

For $m^{2}_{H_{u}}=m_{H_{d}}^{2}=\mathcal{O}(m_{0}^{2})$, radiative corrections from the stop mass squared leads to a relatively large $|m^{2}_{H_{u}}|$ at the weak scale and it requires a large $|\mu|$. That is, the value of $|\mu|$ is determined by the minimization condition of Higgs potential,

[TABLE]

for ${\rm tan}\beta\gg 1$. Here, $m_{Z}$ is the mass of the $Z$ boson, $m_{Z}\simeq 91$ GeV. From the numerical analysis, $|m_{H_{u}}^{2}|$ at the electroweak scale is $\mathcal{O}(10)^{2}\,{\rm TeV}^{2}$ in the parameter space of our interest, and thus $|\mu|=\mathcal{O}(10)\,$TeV.

To explain the muon $g-2$, we require that the masses of the electroweakinos are also in $\mathcal{O}(100-1000)$ GeV. We consider two cases, the universal gaugino mass and the non-universal gaugino masses. In the case of the universal gaugino mass, we take

[TABLE]

For the non-universal gaugino masses, we take

[TABLE]

The reason of the choice of the sign of $M_{3}$ will be explained in Sec. V.

The SUSY contributions to the muon $g-2$ are proportional to ${\rm tan}\beta$ and we focus on

[TABLE]

The large ${\rm tan}\beta$ is also advantageous to explain the observed Higgs boson mass with $m_{3}=\mathcal{O}(10)$ TeV. As we will also discuss, the large ${\rm tan}\beta$ is important to achieve the bottom-tau unification (See Sec. V).

Throughout this paper, we assume that the SUSY breaking parameters do not have CP violating phases and the only source of the CP violation in the following analysis comes from the CP phase in the CKM matrix. Under these assumptions, we take the convention of $\mu>0$ without loss of generality (See the following section for more details).

III.2 SUSY Contributions To Muon $g-2$

Here, we show the parameter space to explain the muon $g-2$. In the MSSM, the relevant one loop contributions to the muon $g-2$ come from the diagrams with the smuon/neutralino or the muon-type sneutrino/chargino loops Cho et al. (2011). In our scenario, the one-loop bino-smuon diagram dominates the SUSY contributions to the muon $g-2$, the SUSY contribution is proportional to the $\mu$ parameter. The Higgsino contributions are, on the other hand, suppressed due to their heaviness.

In our analysis, we use the package SPheno-4.0.3999We slightly modify SPheno-4.0.3 to calculate the SUSY spectrum on the basis in Eq. (9). to calculate the low energy spectrum from the input parameter at the GUT scale. In the code, the renormalization group equations are solved including the flavor mixing parameters, which are relevant for analyses in the next section. We also use the package FeynHiggs2.14.3 Frank et al. (2007); Heinemeyer et al. (1999, 2000); Hahn et al. (2009); Degrassi et al. (2003) to compute the muon $g-2$ and the Higgs mass from the low energy spectrum obtained by SPheno.101010We added the two loop corrections for the large ${\rm tan}\beta$ given in Marchetti et al. (2009) for the SUSY contributions to the muon $g-2$. Although we use FeynHiggs to calculate the muon $g-2$, the following results are not changed even if we use SPheno.

III.2.1 Universal Gaugino Mass at the GUT Scale

For the universal gaugino mass, we search for the region to explain the observed muon $g-2$ on the $m_{0}\text{--}M_{1/2}$ plane for $m_{3}=12\,{\rm TeV}$ or $10\,{\rm TeV}$. In Fig. 1, the observed muon $g-2$ is explained in the orange (yellow) regions within $1\sigma$ ($2\sigma$). The gray shaded regions are excluded by the tachyonic masses of the sleptons in the first two generations due to the large two-loop renormalization group effects from the third generation masses. In most of the favored parameter space, the SM Higgs mass is in $124{\rm-}126\,{\rm GeV}$ for the central value of the measured top quark mass $m_{\rm top}=173.1\pm 0.9$ GeV Tanabashi et al. (2018). That is consistent with the observed Higgs boson mass $m_{H}=125.18\pm 0.16$ GeV Tanabashi et al. (2018) within the theoretical uncertainty.111111The theoretical uncertainty of the Higgs mass is about $2-3$ GeV Athron et al. (2017). We have also checked that the Higgs boson mass is consistent with the observed value within the theoretical error in Fig. 2, 3, 4, and Fig. 6. In the left figure, the measured muon $g-2$ is explained for slightly larger $M_{1/2}$ and $m_{0}$ compared with the right one. This is because the $\mu$ parameter is larger for the larger $m_{3}$, with which the SUSY contribution to the muon $g-2$ is slightly enhanced.

III.2.2 Non-Universal Gaugino Masses at the GUT Scale

For the non-universal gaugino masses, we also study the parameter space explaining the muon $g-2$ on the $m_{0}\text{--}M_{2}$ plane. In Fig. 2, we plot the orange (yellow) shaded regions predicting the observed muon $g-2$ within $1\sigma~{}(2\sigma)$. The gray shaded regions are excluded by the tachyonic masses of the sleptons. Here, the ratio of the gaugino masses are fixed to be $M_{1}=1.725\times M_{2}$ ($M_{1}=1.74\times M_{2}$) and $M_{3}=-2.6\,{\rm TeV}$ ($-2.5\,{\rm TeV}$) for $m_{3}=11$ TeV ($10$ TeV) at the GUT scale, respectively. The motivations of these choices will be explained in Sec. III.3.2. Again, the SM Higgs boson mass is consistent with the observed Higgs boson mass in the favored parameter space. As in the case of the universal gaugino mass, the muon $g-2$ is explained by slightly larger $m_{0}$ and $M_{2}$ for the larger $m_{3}$ when we compare both the figures.

III.3 Collider Constraints

Now, let us discuss the collider constraints. We will discuss cosmology in the next subsection.

III.3.1 Universal Gaugino Mass at the GUT Scale

The bino is the lightest SUSY particle (LSP) on the right sides of the black dashed lines in Fig. 3. In this case, the parameter space with the gluino and the squarks lighter than about $2.6$ TeV is excluded by the search for multi-jets plus missing transverse momentum at ATLAS $13$ TeV using $36$ fb*-1* Aaboud et al. (2018a). We show the excluded regions as the blue shaded ones in Fig. 3, which correspond to the $95\%$ CL limits.

A charged slepton is the LSP on the left side of the black dashed line in Fig. 3. Here, we assume that the LSP has a short lifetime by $R$-parity violation so that the scenario is consistent with cosmology. However, the size of the $R$-parity violation is limited from above not to wash out the baryon asymmetry made by baryogenesis (such as thermal leptogenesis). As a result, the charged slepton LSP is expected to be stable inside the detectors (See $e.g.$ Endo et al. (2010)). The heavy stable charged particles searches in Collaboration (2016) put upper limits on the production cross section of the SUSY particles, which is converted to the constraints on the mass parameters by using the cross-section given in Borschensky et al. (2014).

In Fig. 3, the blue shaded regions on the left of the black dashed lines are excluded by the constraints on the heavy stable charged particle ($95\%$ CL).121212There are two kinds of blue regions (the blue and the lighter blue regions). The blue plus light blue shaded regions are excluded here. We will explain this difference soon. In the figure, almost entire region favored by the muon $g-2$ is excluded except for a tiny region near $(m_{0},\,M_{1/2})=(0\,{\rm GeV},\,1.7\,{\rm TeV})$ for $m_{3}=12\,{\rm TeV}$ in the case of the universal gaugino mass.131313The muon $g-2$ can be explained for a larger $M_{1/2}$ for a larger ${\rm tan}\beta$. In this case, however, the CP-odd neutral Higgs scalar becomes tachyonic. In Tab. 1, we show a sample spectrum which evades the constraints while explains the muon $g-2$ within $2\sigma$ in the charged slepton LSP scenario.

In the analysis for the charged slepton LSP, the electroweakino productions are the dominant SUSY production modes, where we use the production cross section given in Fuks et al. (2012, 2013). For comparison, we also show the constraints assuming the SUSY production cross section for the degenerate squarks and gluino (the lighter blue shaded regions). In the actual spectrum, the gluino is heavier than the lightest squark in most of the favored parameter space, and hence, the production cross section of the colored SUSY particles Borschensky et al. (2014) is smaller than the degenerated case. As the figure shows, the constraints via the colored SUSY particle production are at most comparable or weaker than those from the electroweakino productions for the heavy stable charged particle searches.141414For the searches of multi-jets plus missing transverse energy which is relevant for the bino search, the colored SUSY production plays the dominant roles.

III.3.2 Non-Universal Gaugino Masses at the GUT Scale

In the case of the non-universal gaugino masses, the constraints from the collider searches can be weakened by several reasons. For the bino LSP cases, for example, the production cross section of the colored SUSY particles which are relevant for multi-jets plus missing transverse energy search is reduced if the gluino and squarks are heavy. The constraints from the heavy stable charged particle searches can be also evaded since the sneutrino can be lighter than the charged sleptons. In the following, we again take $M_{3}=-2.6\,{\rm TeV}$ or $-2.5\,{\rm TeV}$ at the GUT scale, which suppresses the colored SUSY particle production cross section. The choice of the sign of $M_{3}$ will be relevant in the discussion in Sec. V. We also take $M_{1}=1.725\times M_{2}$ or $M_{1}=1.74\times M_{2}$ with which the sneurinos are lighter than the charged sleptons in the favored parameter space.

On the right side of the black dashed line in Fig. 4, the neutralino is the LSP. As we have mentioned, there is no stringent collider constraint because the squarks and the gluinos become heavy by a rather large $|M_{3}|$. It should be also noted that the constraints on the missing transverse momentum from the electroweakino production are far less relevant due to the rather degenerate electroweakino spectrum, $(m_{\chi_{1}^{\pm}}-m_{\chi_{0}^{1}})/m_{\chi_{1}^{\pm}}\lesssim 30\%$ Aaboud et al. (2018b); Sirunyan et al. (2018).

Eventually, the searches for the missing transverse energy with the charged leptons put the most stringent constraints on the neutralino LSP region. In the figure, the blue shaded regions in Fig. 4 are excluded by the constraint in Aaboud et al. (2018c) ($95\%$ CL).151515The constraint from the latest result in collaboration (2019) is not stringent, where it is assumed that the slepton decays into a lepton and a neutralino with a $100$% branching ratio. In our scenario, a slepton also decays into a chargino and a neutrino. Furthermore, a charged lepton from a chargino decay becomes soft. Thus, the constraint is weakened. In Tab. 2, we show a sample spectrum which evades the constraints while explains the muon $g-2$ within $1\sigma$ in the neutralino LSP scenario.

On the left side of the black dashed line in the favored parameter space in Fig. 4, the sneturino is the LSP due to a relatively small $M_{2}$ compared with $M_{1}$. Because of the heavy colored SUSY particles and the degenerate electroweakinos, this region is less constrained by multi-jet plus missing energy searches Aaboud et al. (2018a, b); Sirunyan et al. (2018). The search for the missing transverse energy with the charged leptons are also less sensitive to this region. In the sneutrino LSP region, however, the constraints from the direct detection of dark matter are stringent if the sneutrino LSP is stable (See the following section for more detail discussion).

III.4 Cosmology

As we have seen above, most parameter region favored by the muon $g-2$ has been excluded for the universal gaugino mass by the LHC searches. In this subsection, we discuss cosmology focusing on the non-universal gaugino masses.

In most of the neutralino LSP region, the bino is the dominant component of the lightest neutralino. In general, the bino LSP is disfavored from a too large thermal relic abundance due to its small annihilation cross section.161616In some of the neutralino LSP region, the wino is the LSP. There, we have confirmed that the constraints from the direct detection are negligible due to the small thermal relic abundance of the wino. This problem can be evaded for our particular choice of the gaugino mass parameters,

[TABLE]

which leads to the rather degenerate physical bino and wino masses. In this case, the co-annihilation between the bino and the wino is efficient Birkedal-Hansen and Nelson (2003); Baer et al. (2005); Arkani-Hamed et al. (2006); Harigaya et al. (2014)171717See the other neutralino coannihilation with other light sparticles $e.g$ a stop Boehm et al. (2000), stau Ellis et al. (1998), gluino Profumo and Yaguna (2004)., which makes the bino abundance consistent with the observed dark matter density.

In Fig. 4, we draw the green line on which the thermal relic abundance of the bino corresponds to the observed dark matter density. To calculate the thermal relic abundance of dark matter (and its cross-section with nucleons), we use the package MicrOMEGAs_5.0.4 Bélanger et al. (2018). For the choice of the gaugino mass relation in Eq. (21), we find that the bino mass around $400$ GeV can explain the observed dark matter density. By tuning the ratio between $M_{2}$ and $M_{1}$ further, the dark matter density can be explained for the bino mass up to around $700$ GeV within the parameter region favored by the muon $g-2$.

As the bino is the LSP, the direct detection cross is highly suppressed. In fact, we find,

[TABLE]

where $m_{\rm DM}$ denotes the dark matter mass. This cross section is much smaller than the current constraint Akerib et al. (2017); Cui et al. (2017); Aprile et al. (2018).

For the sneutrino LSP, on the other hand, the relic abundance is much smaller than the observed dark matter density ($\Omega h^{2}=10^{-2}$) due to its large annihilation cross section in the parameter region favored by the muon $g-2$. Even with such a small relic abundance, however, the sneutrino LSP contribution to dark matter has been excluded by the direct detection experiments since it has a large scattering cross section with the nucleons, $\sigma_{\rm SI}\simeq\mathcal{O}(10^{-5})\,{\rm pb}$.

As a result, we find that tiny $R$-parity violation is required in the case of the sneutrino LSP as in the case of the charged slepton LSP (See the previous discussion in Sec. III.3.1). In those cases, we need dark matter candidate other than the LSP.

IV FCNC constraints in Split-Family SUSY model

.

In the Split-Family model, there is a non-trivial enhancement of the FCNC by the non-universality of the sfermion masses. In this section, we investigate the FCNC constraints on the model for the minimal and small mixing scenarios defined in Sec. II. In the minimal mixing scenario, the CKM matrix is the only source of the flavor mixing. We show that the FCNC constraints on the minimal scenario are not stringent. For the small mixing scenario, we demonstrate how large flavor mixing in the supersymmetry breaking parameters is allowed.

IV.1 Experimental FCNC limits

Let us first summarize the FCNC constraints relevant to the mixing parameters.

Meson mixing

The CP-violating parameter in the neutral kaon system $\epsilon_{K}$ gives the stringent constraint on the squark flavor mixing. The measured value of $\epsilon_{K}$ Tanabashi et al. (2018) is

[TABLE]

For the theoretical prediction of SM, we adopt Jang et al. (2018),

[TABLE]

which is derived by using the QCD sum-rule.181818We do not use the theoretical prediction on $\epsilon_{K}$ from the lattice QCD Jang et al. (2018).

The mass difference between the long-lived Kaon and the short-lived Kaon, ${\mit\Delta}M_{K}=m_{K^{0}_{L}}-m_{K^{0}_{S}}$, also puts the constraints on the mixing parameters. The measured ${\mit\Delta}M_{K}$ Tanabashi et al. (2018) is given by

[TABLE]

The SM prediction of ${\mit\Delta}M_{K}$, however, has a large uncertainty due to the unknown long-distance effects (see $e.g.$ Ref Grinstein (2016)). Thus, in this paper, we will just compare ${\mit\Delta}M_{K}^{\rm exp}$ with the order of the magnitude of SUSY contributions. The neutral $D$-meson oscillation may also put the constraints. The observed mass difference Tanabashi et al. (2018) is

[TABLE]

The theoretical calculation of ${\mit\Delta}M_{D}$ has large uncertainties as in the case of ${\mit\Delta}M_{K}$. Again, we only evaluate the order of the magnitude of the SUSY contribution with the $|{\mit\Delta}M_{D}^{\rm ex}|$.

Lepton Flavor Violation

For the slepton flavor mixings, the stringent constraint comes from the search for the decay $\mu^{+}\rightarrow e^{+}\gamma$ at the MEG experiment Mori (2017). The upper limit on the branching ratio of this process is

[TABLE]

IV.2 FCNC Constraints in the Minimal Mixing Scenario

Here, we investigate the FCNC constraints in the minimal mixing scenario. We search for the parameter space which is the same as the previous section in Fig. 1 and Fig. 2. In this scenario, the CKM matrix appearing in the Yukawa couplings in Eq. (9) leads to the flavor mixing in the squark mass matrices at the weak scale due to the non-universality of the split family structure at the GUT scale.

Such flavor mixing is constrained by $\epsilon_{K}$, where the CP violation comes solely from the CP phase of the CKM matrix. In Fig. 3, the red hatched regions are excluded. There, the region with the small $m_{0}$ and $M_{1/2}$ is excluded due to the light gluino and squarks. Notice that the SUSY FCNC contribution is larger for a larger $m_{3}^{2}/m_{0}^{2}$, and thus the constraint is severe for the left figure.191919In our analysis, we fully diagonalize the squark masses without using the mass insertion technique (See appendix A for more details). In Fig. 4, on the other hand, no constraint appears from $\epsilon_{K}$. This is because the SUSY contributions are suppressed due to the heavy gluino and squarks for $|M_{3}|\simeq 2.5-2.6$ TeV.

As a result, the Split-Family SUSY model is not stringently constrained by the flavor violation in the minimal mixing scenario. We have also confirmed that ${\mit\Delta}M_{K_{\rm SUSY}}$ and ${\mit\Delta}M_{D_{\rm SUSY}}$ become much smaller than the observed one (See $e.g.$ Tab. 1 and Tab. 2).202020In the minimal mixing scenario, we have also confirmed that the neutron electric dipole moment (EDM) from the CKM phase does not lead to the stringent constraint. See also Ref. Endo et al. (2011) for the EDM constraints.

Several comments are in order. As we have mentioned earlier, we may consider another simple family basis in Eq. (10). Since the flavor mixing is dominated by the renormalization group effects, we obtained similar constraint even in this case for the large ${\rm tan}\beta$.

One may also wonder how large SUSY FCNC contributions are expected if the first and the second generation soft SUSY breaking masses are not degenerate. We discuss this possibility in appendix B. There, the model is severely constrained by the FCNC even for the minimal mixing scenario.

IV.3 FCNC Constraints in Small Mixing Scenario

Now, let us discuss the FCNC problems in the small mixing scenario. As we will see, the flavor mixing in the slepton sector leads to sizable FCNC contributions because their masses are of $\mathcal{O}(100)$ GeV to explain the muon $g-2$.

IV.3.1 FCNC Constraints on Slepton Flavor Mixing

As we explained in Sec. II, we take the diagonal lepton Yukawa coupling while allowing small flavor mixing in the SUSY breaking soft masses as in Eq. (14). We use the following parametrization for the slepton mixing matrix $V_{\rm mix}$,

[TABLE]

Here, we again assume that no CP violation appears in the mixing matrix.

As the split-family structure is motivated by the Yukawa hierarchy, we also assume that the flavor mixing in the SUSY breaking sector has a structure similar to the CKM matrix. That is, we assume that the mixing angles are controlled by a small parameter $\epsilon$,

[TABLE]

where $\epsilon\simeq 0.2$ roughly mimics the mixing angles in the CKM matrix. In the following, we demonstrate how large flavor mixing is allowed in the slepton sector by using $\epsilon$.

In our numerical analysis, we take the following mixing angles,

[TABLE]

where $R_{i}$ shows the random numbers between $(0.5,~{}1.5)$. The random numbers $R_{i}$ for $m_{\tilde{L}}^{2}$ and $m_{\tilde{E}}^{2}$ are taken independently.

In Fig. 5, we show how large slepton mixing is tolerable from the current and future constraints on $\mathcal{B}(\mu^{+}\rightarrow e^{+}\gamma)$. There, we consider the input parameters in Tab. 1 and Tab. 2 for examples. The blue curve line denotes the model predictions for $\mathcal{B}(\mu^{+}\rightarrow e^{+}\gamma)$. The horizontal (dashed) lines are the current (future expected Baldini et al. (2018); Hew (2012)) experimental bounds. From the figures, the constraint from $\mathcal{B}(\mu^{+}\rightarrow e^{+}\gamma)$ requires $\epsilon\lesssim 0.06$.212121We checked that the constraints are not significantly changed even if we put additional CP phases in $V_{\rm mix}$ of the sleptons. Detail analysis will be done elsewhere.

Before closing this section, let us comment on the effect of the PMNS matrix. As in the case of the squark mixing, the PMNS matrix could lead to large flavor mixing in the slepton mass matrix at the weak scale even if we assume a diagonal soft masses (See Eq. (5)). However, the FCNC constraints due to the PMNS matrix are not stringent if the neutrino Yukawa matrix is small enough as long as the right-handed neutrino mass $M_{R}\lesssim 10^{10}$ GeV (See appendix C for more details).

IV.3.2 FCNC Constraints on Squark Flavor Mixing

In the small mixing scenario, we can also introduce $V_{\rm mix}$ to the squarks with the mixing angles given by Eq. (38). In this case, however, the FCNC constraints on the mixing angle $\epsilon$ are much weaker than the slepton case, and hence, we do not discuss them any further.

V Bottom-tau unification

In the Split-Family SUSY model, the muon $g-2$ can be explained by the small $|M_{1}|,~{}|M_{2}|,~{}m_{0}^{2}$, and the large ${\rm tan}\beta$. In addition, the stringent limits from the LHC experiments require the rather large $|M_{3}|$. Interestingly, these parameter sets are found to be appropriate to realize the bottom-tau unification Ajaib et al. (2014).

The bottom Yukawa coupling receives threshold corrections from the gluino loop diagrams Banks (1988); Hall et al. (1994); Hempfling (1994); Blazek et al. (1995),

[TABLE]

where $y_{b}$ is the bottom Yukawa coupling in the MSSM, and ${\mit\Delta}y_{b}$ is the threshold correction. The negative $M_{3}\,\mu\,{\rm tan}\beta$ gives the negative ${\mit\Delta}y_{b}$, which makes the bottom Yukawa coupling in the MSSM $y_{b}$ larger for a given bottom quark mass.222222In our notation, $y_{b}$ is real positive (See Eq. (9)). For the tau Yukawa coupling, on the other hand, the threshold corrections are small due to the small $M_{1}$ and $M_{2}$ as well as the small gauge coupling constants. As a result, the large negative contribution in Eq. (40) makes $y_{b}/y_{\tau}$ large at the weak scale, which is appropriate for the bottom-tau unification at the GUT scale. This is the reason that we take negative $M_{3}$ in the previous section.

In Fig. 6, we show the parameter space for the successful bottom-tau unification. There, we assume the input parameters in Tab. 2 other than $m_{0}$ and ${\rm tan}\beta$. The red lines are the contours of $|y_{b}/y_{\tau}-1|=~{}1.5\%,~{}3\%,~{}5\%$, where $y_{b}$ and $y_{\tau}$ are the bottom and tau Yukawa couplings at the GUT scale, respectively.232323The threshold correction from the gluino is sensitive to the matching scale, which requires careful treatments of the decoupling and the renormalization group effects Chigusa and Moroi (2016). The other color codes are the same as Fig. 4. We have checked that the Higgs mass is in $124{\rm-}126$ GeV in most of the favored region. We have also confirmed that the dark matter abundance is $\Omega h^{2}\simeq 0.12$ on the right side of the dotted line (bino LSP region) in the figure.242424On the left side of the dotted line, the sneutrino is the LSP as discussed in Sec. III. Thus, we find the parameter space favored by the muon $g-2$ is also appropriate to achieve the bottom-tau unification.

VI Conclusion

In this paper, we have revisited the Split-Family SUSY model. In the model, the sfermion masses of the first two generations are in the hundreds GeV range, while that of the third generation is in the tens TeV range. With this spectrum, the deviation of the muon $g-2$ and the observed Higgs boson mass are explained simultaneously.

In Sec. III, we have first shown the parameter space to explain the muon $g-2$ and the Higgs boson mass. In our analysis, we have searched for two cases of the universal gaugino mass and the non-universal gaugino masses. For the universal gaugino mass, almost the entire region to explain the muon $g-2$ within $2\sigma$ is excluded by the collider searches as shown in Fig. 3. This is due to the lightness of the squarks, the gluino, and the wino masses. For the non-universal gaugino masses, the gluino can be heavier with which the collider constraints can be easily evaded (See Fig. 4).252525See $e.g.$ Ref. Endo et al. (2013, 2014) for the discussion about the future sensitivity at the LHC and the ILC. We have also found the parameter space where the bino LSP can explain the observed dark matter density thanks to the coannihilation with the wino.

In Sec. IV, we have studied the FCNC problems in the Split-Family SUSY model. We have searched for two scenarios, $i.e.$ the minimal mixing scenario and the small mixing scenario (See Sec. IV for details). In the minimal scenario, we have assumed that the CKM matrix is the only source of the flavor mixing. There, we have shown that the SUSY FCNC contributions are small enough to evade the problem.

For the small mixing scenario, we have assumed the CKM like mixing matrix to the soft mass parameters (See Eq. (37) and around it). Then, we have demonstrated how large flavor mixing is allowed in the slepton sector. There, the most stringent constraint comes from the lepton flavor violation decay of $\mu^{+}\to e^{+}+\gamma$. We have shown that the mixing angles have to be relatively small, $\epsilon\lesssim 0.06$.

In Sec. V, we have discussed one bonus feature of the model, the bottom-tau unification. For the successful bottom-tau unification, the large threshold correction is required. Interestingly, such parameter space is compatible with the one favored by the muon $g-2$. In Fig. 6, we have shown that the bottom-tau unification is significantly improved for the large ${\rm tan}\beta$.

Several comments are in order. First, it is possible to achieve a small $\mu$ term in the Split-Family SUSY model. In fact, for $m_{H_{u,d}}^{2}=m_{3}^{2}$, the focus point mechanism Feng et al. (2000a, b) results in a small $\mu$ term. In such a case, the neutralino LSP can have a sizable Higgsino contribution, so that the dark matter-nucleon cross section becomes large.

Throughout this paper, we have assumed that the SUSY breaking parameters do not have any CP-violating phases. In fact, these are strong assumptions and it is highly non-trivial to achieve such soft SUSY breaking parameters from high energy theory. The CP violating phases in the SUSY breaking parameters are constrained by the measurements of the EDMs, which will be discussed elsewhere.

Acknowledgements.

This work is supported by JSPS KAKENHI Grant Numbers JP26104001 (T.T.Y), JP26104009 (T.T.Y), JP16H02176 (T.T.Y), JP17H02878 (M. I. and T. T. Y.), No. 15H05889 and No. 16H03991 (M. I.), JP15H05889 (N.Y.), JP15K21733 (N.Y.), JP17H05396 (N.Y.), JP17H02875 (N.Y.), and by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan (T.T.Y.). The work of M. S. is supported in part by a Research Fellowship for Young Scientists from the Japan Society for the Promotion of Science (JSPS).

Appendix A Validity of the Mass Insertion Approximation

In this appendix, we briefly discuss the validity of the mass insertion approximation for the SUSY FCNC contributions. For simplicity, we consider the case of the universal gaugino mass in the minimal mixing scenario.

In the mass insertion approximation, the gluino contribution to $\epsilon_{K}$ is estimated at the next leading order Ciuchini et al. (1998),

[TABLE]

Here, $\alpha_{s}=0.1184$, $m_{K}=0.498$ GeV, the Kaon decay constant $f_{K}=0.16$ GeV. The mixing parameter ${\mit\Delta}^{ds}_{LL}$ denotes an off-diagonal element for the left-handed down and strange squarks soft SUSY breaking mass squared matrix normalized by $1/m_{0}^{2}$,

[TABLE]

From Eq. (41), (42), $\epsilon_{K_{SUSY}}$ is enhanced by a factor of $m_{3}^{4}$.

The above approximation is, however, only valid when the sfermion mass splitting is small, $i.e.$ $m_{0}^{2}\simeq m_{3}^{2}$. For $m_{3}^{2}\gg m_{0}^{2}$, on the other hand, the approximation leads to the overestimation. Thus, we need to calculate the FCNC processes in the exact mass diagonalization. For this purpose, we use the package susy_flavor_2_54 Rosiek (2015); Crivellin et al. (2013); Rosiek et al. (2010).

The numerical result of $\epsilon_{K_{\rm SUSY}}$ as a function of $m_{3}$ is shown in Fig. 7. The blue dotted line corresponds to the result from the code, where we take the input parameter $m_{0}=1\,{\rm TeV},~{}M_{1/2}=1.5\,{\rm TeV},~{}{\rm tan}\beta=40,~{}m_{H_{u,d}}^{2}=(100\,{\rm GeV})^{2}$.262626Here, we only include the contribution from the $VLL$ four-quark operator Crivellin et al. (2013) to compare with the approximation in Eq. (41). We also switched off the resummation of chirally enhanced corrections Crivellin et al. (2013) to compare with the mass insertion approximation. The gray solid line is given by the mass insertion approximation in Eq. (41) for the same input parameters. Here, these input parameters are set at the SUSY scale to make comparison easier. From the figure, we confirm the overestimation of $\epsilon_{K}$ for $m_{3}^{2}\gg m_{0}^{2}$ in the mass insertion approximation.

Appendix B FCNC Constraints on First and Second Sfermion Soft Mass Splitting

In this appendix, we discuss how large mass splitting is allowed for the first and the second generation squarks from the FCNC constraints. To parametrize the splitting, we use the following non-universal soft masses for the squarks,

[TABLE]

where $m_{2}$ denotes the soft mass parameter of the second generation squarks.

In Fig. 8, we show the FCNC constraints for the case of the universal gaugino mass (the figure in the left). The gray shaded region is excluded by $\epsilon_{K}$ ($90\%$ CL). There, we use the input parameters in Tab. 1, but take $m_{0}$ and $m_{2}$ as variables. The parameter space with the ratio $m_{2}/m_{0}=\mathcal{O}(10)$ is excluded for $m_{0}\gtrsim 60$ GeV, where $m_{2}$ becomes comparable to the gluino mass contributions to the squark masses in the second generation.

In the right panel of Fig. 8, we also show the constraints for the non-universal gaugino masses. There, we use the input parameters in Tab. 2. The SUSY FCNC contributions are large with the splitting, $m_{2}/m_{1}\gtrsim 5$, where $m_{2}$ becomes comparable to the gluino mass contributions to the squark masses in the second generation.

Appendix C Lepton Flavor Violation from MNS Matrix Effect

In this appendix, we briefly discuss the constraint from the decay $\mu^{+}\to e^{+}+\gamma$, where the MNS matrix is taken into account as in the case of the minimal mixing scenario in the squark sector.

In the presence of the right-handed neutrinos, the superpotential of the lepton sector is given by,

[TABLE]

Here, $N_{Ri}\,(i=1{\rm-}3)$ are the right-handed neutrinos with three flavors. As in the case of the minimal mixing scenario in the squark sector, we assume the soft masses in Eq. (5) while the Yukawa couplings are diagonal matrices ($\hat{f}_{\nu}$) with the MNS matrix ($U_{\rm MNS}$).

In this setup, the flavor mixing parameter between the first and second generation of the charged selectrons is given by,

[TABLE]

where the parameter is defined by an off-diagonal element for the left-handed electron and muon sleptons soft SUSY breaking mass squared matrix normalized by $1/m_{0}^{2}$.

Then, the branching ratio of the process $\mu^{+}\to e^{+}+\gamma$ is given in the following Hisano et al. (1996); Arganda and Herrero (2006); Paradisi (2005),

[TABLE]

Here, $G_{F}=1.166\times 10^{-5}\,{\rm GeV}^{-2}$ is the Fermi coupling decay constant, $\alpha_{\rm em}\simeq 1/137$ is the fine-structure constant, $\alpha_{1}\simeq 1/60$ is the weak coupling constant. We also take $m_{\rm slepton}\simeq 480$ GeV, $\mu\simeq 9\times 10^{3}$ GeV, and $m_{\rm bino}\simeq 400$ GeV as in Tab. 2. The coefficient $c_{1}$ is about $0.1$ Hisano et al. (1996); Arganda and Herrero (2006); Paradisi (2005). The branching ratio is consistent with the current upper-limit $\mathcal{B}\lesssim 10^{-13}$ for $M_{R}\lesssim 10^{10}\,{\rm GeV}$. More detailed analysis remains for future work.