# Effects of supersymmetric CP violating phases on 𝐵→𝑋_𝑠⁢𝑙⁺⁢𝑙⁻ and ϵ_𝐾

## Abstract

We consider effects of CP violating phases in $`\mu `$ and $`A_t`$ parameters in the effective supersymmetric standard model on $`BX_sl^+l^{}`$ and $`ϵ_K`$. Scanning over the MSSM parameter space with experimental constraints including edm constraints from Chang-Keung-Pilaftsis (CKP) mechanism, we find that the $`\mathrm{Br}(BX_sl^+l^{})`$ can be enhanced by upto $`85\%`$ compared to the standard model (SM) prediction, and its correlation with $`\mathrm{Br}(BX_s\gamma )`$ is distinctly different from the minimal supergravity scenario. Also we find $`1ϵ_K/ϵ_K^{SM}1.4`$, and fully supersymmetric CP violation in $`K_L\pi \pi `$ is not possible. Namely, $`|ϵ_K^{\mathrm{SUSY}}|O(10^5)`$ if the phases of $`\mu `$ and $`A_t`$ are the sole origin of CP violation.

preprint:     April, 1999 KAIST-TH 99/1  hep-ph/9904283 SNUTP 99-003      

1. In the minimal supersymmetric standard model (MSSM), there can be many new CP violating (CPV) phases beyond the KM phase in the standard model (SM) both in the flavor conserving and flavor violating sectors. The flavor conserving CPV phases in the MSSM are strongly constrained by electron/neutron electric dipole moment (edm) and believed to be very small ($`\delta 10^2`$ for $`M_{\mathrm{SUSY}}O(100)`$ GeV ) . Or, one can imagine that the 1st/2nd generation scalar fermions are very heavy so that edm constraints are evaded via decoupling even for CPV phases of order $`O(1)`$ . Also it is possible that various contributions to electron/neutron EDM cancel with each other in substantial parts of the MSSM parameter space even if SUSY CPV phases are $`O(1)`$ and SUSY particles are relatively light  . In the last two cases where SUSY CPV phases are of $`O(1)`$, these phases may affect $`B`$ and $`K`$ physics in various manners. In the previous letter , we presented effects of these SUSY CPV phases on $`B`$ physics : the $`B^0\overline{B^0}`$ mixing and the direct asymmetry in $`BX_s\gamma `$, assuming that EDM constraints and SUSY FCNC problems are evaded by heavy 1st/2nd generation scalar fermions. In this letter, we extend our previous work to $`BX_sl^+l^{}`$ and $`ϵ_K`$ (see also Ref. .) within the same assumptions.

An important ingredient for large $`\mathrm{tan}\beta `$ in our model is the constraint on the $`\mu `$ and $`A_t`$ phases coming from electron/neutron edm’s through Chang-Keung-Pilaftsis (CKP) mechanism . Two loop diagrams with CP-odd higgs and photon (gluon) exchanges between the fermion line and the sfermion loop (mainly stops and sbottoms) can contribute significantly to electron/neutron edm’s in the large $`\mathrm{tan}\beta `$ region. The authors of Ref.  find that

$$(\frac{d_f}{e})_{\mathrm{CKP}}=Q_f\frac{3\alpha _{\mathrm{em}}}{64\pi ^2}\frac{R_fm_f}{M_A^2}\underset{q=t,b}{}\xi _qQ_q^2\left[F\left(\frac{M_{\stackrel{~}{q}_1}^2}{M_A^2}\right)F\left(\frac{M_{\stackrel{~}{q}_2}^2}{M_A^2}\right)\right],$$
(1)   
where $`R_f=\mathrm{cot}\beta (\mathrm{tan}\beta )`$ for $`I_{3f}=1/2(1/2)`$, and

$$\xi _t=\frac{\mathrm{sin}2\theta _{\stackrel{~}{t}}m_t\mathrm{Im}(\mu e^{i\delta _t})}{\mathrm{sin}^2\beta v^2},\xi _b=\frac{\mathrm{sin}2\theta _{\stackrel{~}{b}}m_b\mathrm{Im}(A_be^{i\delta _b})}{\mathrm{sin}\beta \mathrm{cos}\beta v^2},$$
(2)   
with $`\delta _q=\mathrm{Arg}(A_q+R_q\mu ^{})`$, and $`F(z)`$ is a two-loop function given in Ref. . This new contribution is independent of the 1st/2nd generation scalar fermion masses, so that it does not decouple for heavy 1st/2nd generation scalar fermions. Therefore it can be important for the electron or down quark edm for the large $`\mathrm{tan}\beta `$ case. This is in sharp contrast with the usual one-loop contributions to edm’s, for which 

$$\left(\frac{d_f}{e}\right)10^{25}\mathrm{cm}\times \frac{\{\mathrm{Im}\mu ,\mathrm{Im}A_f\}}{\mathrm{max}(M_{\stackrel{~}{f}},M_\lambda )}\left(\frac{1\mathrm{TeV}}{\mathrm{max}(M_{\stackrel{~}{f}},M_\lambda )}\right)^2\left(\frac{m_f}{10\mathrm{MeV}}\right),$$
(3)   
and one can evade the edm constraints by having small phases for $`\mu ,A_{e,u,d}`$, or heavy 1st/2nd generation scalar fermions. However, this would involve enlargement of our model parameter space, since one has to consider the sbottom sector as well as the stop sector. Therefore, more parameters have to be introduced in principle : $`m_{\stackrel{~}{b}}^2`$ and $`A_b`$ where $`A_b`$ may be complex like $`A_t`$. In order to avoid such enlargement, we will assume that there is no accidental cancellation between the stop and sbottom loop contributions.

This CKP edm constraint has not been included in the recent paper by Demir et al. , who made claims that there could be a large new phase shift in the $`B^0\overline{B^0}`$ mixing and it is possible to have a fully supersymmetric $`ϵ_K`$ from the phases of $`\mu `$ and $`A_t`$ only. However, if $`\mathrm{tan}\beta `$ is large ($`\mathrm{tan}\beta 60`$) as in Ref. , the CKP edm constraints via the CKP mechanism have to be properly included. This constraint reduces the possible new phase shift in the $`B^0\overline{B^0}`$ mixing to a very small number, $`2|\theta _d|1^{}`$, as demonstrated in Fig. 1 (a) of Ref. . On the other hand, the CKP edm constraint does not affect too much the direct CP asymmetry in $`BX_s\gamma `$ .

In this work, we continue studying the effects of the phases of $`\mu `$ and $`A_t`$ on $`BX_sl^+l^{}`$ and $`ϵ_K`$. We also reconsider a possibility of fully supersymmetric CP violation, namely generating $`ϵ_K`$ entirely from the phases of $`\mu `$ and $`A_t`$ with vanishing KM phase ($`\delta _{\mathrm{KM}}=0`$). Our conclusion is at variance with the claim made in Ref. .

2. As in Refs.  , we assume that the 1st and the 2nd family squarks are degenerate and very heavy in order to solve the SUSY FCNC/CP problems. Only the third family squarks can be light enough to affect $`B`$ and $`K`$ physics. We also ignore possible flavor changing squark mass matrix elements that could generate gluino-mediated flavor changing neutral current (FCNC) processes in addition to those effects we consider below, relegating the details to the existing literature -. Therefore the only source of the FCNC in our case is the CKM matrix, whereas there are new CPV phases coming from the phases of $`\mu `$ and $`A_t`$ parameters (see below), in addition to the KM phase $`\delta _{KM}`$. Definitions for the chargino and stop mass matrices are the same as Ref. . There are two new flavor conserving CPV phases in our model, $`\mathrm{Arg}(\mu )`$ and $`\mathrm{Arg}(A_t)`$ in the basis where $`M_2`$ is real.

We scan over the MSSM parameter space as in Ref.  indicated below (including that relevant to the EWBGEN scenario in the MSSM) :

$`80\mathrm{GeV}<|\mu |<1\mathrm{TeV},`$ $`80\mathrm{GeV}<M_2<1\mathrm{TeV},`$ (4)   
$`60\mathrm{GeV}<M_A<1\mathrm{TeV},`$ $`2<\mathrm{tan}\beta <70,`$ (5)   
$`(130\mathrm{GeV})^2<M_Q^2`$ $`<`$ $`(1\mathrm{TeV})^2,`$ (6)   
$`(80\mathrm{GeV})^2<M_U^2`$ $`<`$ $`(500\mathrm{GeV})^2,`$ (7)   
$`0<\varphi _\mu ,\varphi _{A_t}<2\pi ,`$ $`0<|A_t|<1.5\mathrm{TeV},`$ (8)   
with the following experimental constraints : $`M_{\stackrel{~}{t}_1}>80`$ GeV independent of the mixing angle $`\theta _{\stackrel{~}{t}}`$, $`M_{\stackrel{~}{\chi ^\pm }}>83`$ GeV, and $`0.77R_\gamma 1.15`$ , where $`R_\gamma `$ is defined as $`R_\gamma =BR(BX_s\gamma )^{expt}/BR(BX_s\gamma )^{SM}`$ and $`BR(BX_s\gamma )^{SM}=(3.29\pm 0.44)\times 10^4`$. We also impose $`\mathrm{Br}(BX_{sg})<6.8\%`$ , and vary $`\mathrm{tan}\beta `$ from 2 to 70 This may be too large for perturbation theory to be valid, but we did extend to $`\mathrm{tan}\beta 70`$ in order to check the claims made in Ref. .. This parameter space is larger than that in the constrained MSSM (CMSSM) where the universality of soft terms at the GUT scale is assumed. Especially, our parameter space includes the electroweak baryogenesis scenario in the MSSM . In the numerical analysis, we used the following numbers for the input parameters (running masses in the $`\overline{MS}`$ scheme are used for the quark masses) : $`\overline{m_c}(m_c(pole))=1.25`$ GeV, $`\overline{m_b}(m_b(pole))=4.3`$ GeV, $`\overline{m_t}(m_t(pole))=165`$ GeV, and $`|V_{cb}|=0.0410,|V_{tb}|=1,|V_{ts}|=0.0400`$ and $`\gamma (\varphi _3)=90^{}`$ in the CKM matrix elements.

3. Let us first consider the branching ratio for $`BX_sl^+l^{}`$. The SM and the MSSM contributions to this decay were considered by several groups  and , respectively. We use the standard notation for the effective Hamiltonian for this decay as described in Refs.  and . The new CPV phases in $`C_{7,9,10}`$ can affect the branching ratio and other observables in $`BX_sl^+l^{}`$ as discussed in the first half of Ref.  in a model independent way. In the second half of Ref. , specific supersymmetric models were presented where new CPV phases reside in flavor changing squark mass matrices. In the present work, new CPV phases lie in flavor conserving sector, namely in $`A_t`$ and $`\mu `$ parameters. Although these new phases are flavor conserving, they affect the branching ratio of $`BX_sl^+l^{}`$ and its correlation with $`Br(BX_s\gamma )`$, as discussed in the first half of Ref. . Note that $`C_{9,10}`$ depend on the sneutrino mass, and we have scanned over $`60\mathrm{GeV}<m_{\stackrel{~}{\nu }}<200\mathrm{GeV}`$. In the numerical evaluation for $`R_{ll}\mathrm{Br}(BX_sl^+l^{})/\mathrm{Br}(BX_sl^+l^{})_{\mathrm{SM}}`$, we considered the nonresonant contributions only for simplicity, neglecting the contributions from $`J/\psi ,\psi ^{^{}},etc`$.. It would be straightforward to incorporate these resonance effects. In Figs. 1 (a) and (b), we plot the correlations of $`R_{\mu \mu }`$ with $`\mathrm{Br}(BX_s\gamma )`$ and $`\mathrm{tan}\beta `$, respectively. Those points that (do not) satisfy the CKP edm constraints are denoted by the squares (crosses). Some points are denoted by both the square and the cross. This means that there are two classes of points in the MSSM parameter space, and for one class the CKP edm constraints are satisfied but for another class the CKP edm constraints are not satisfied, and these two classes happen to lead to the same branching ratios for $`BX_s\gamma `$ and $`R_{ll}`$. In the presence of the new phases $`\varphi _\mu `$ and $`\varphi _{A_t}`$, $`R_{\mu \mu }`$ can be as large as 1.85, and the deviations from the SM prediction can be large, if $`\mathrm{tan}\beta >8`$. As noticed in Ref. , the correlation between the $`\mathrm{Br}(BX_s\gamma )`$ and $`R_{ll}`$ is distinctly different from that in the minimal supergravisty case . In the latter case, only the envelop of Fig. 1 (a) is allowed, whereas everywhere in between is allowed in the presence of new CPV phases in the MSSM. Even if one introduces the phases of $`\mu `$ and $`A_0`$ at GUT scale in the minimal supergravity scenario, this correlation does not change very much from the case of the minimal supergravity scenario with real $`\mu `$ and $`A_0`$, since the $`A_0`$ phase becomes very small at the electroweak scale because of the renormalization effects . Only $`\mu `$ phase can affect the electroweak scale physics, but this phase is strongly constrained by the usual edm constraints so that $`\mu `$ should be essentially real parameter. Therefore the correlation between $`BX_s\gamma `$ and $`R_{ll}`$ can be a clean distinction between the minimal supergravity scenario and our model (or some other models with new CPV phases in the flavor changing ).

4. The new complex phases in $`\mu `$ and $`A_t`$ will also affect the $`K^0\overline{K^0}`$ mixing. The relevant $`\mathrm{\Delta }S=2`$ effective Hamiltonian is given by

$$H_{\mathrm{eff}}^{\mathrm{\Delta }S=2}=\frac{G_F^2M_W^2}{(2\pi )^2}\underset{i=1}{\overset{3}{}}C_iQ_i,$$
(9)   
where

$`C_1(\mu _0)`$ $`=`$ $`\left(V_{td}^{}V_{ts}\right)^2\left[F_V^W(3;3)+F_V^H(3;3)+A_V^C\right]`$ (10)   
$`+`$ $`\left(V_{cd}^{}V_{cs}\right)^2\left[F_V^W(2;2)+F_V^H(2;2)\right]`$ (11)   
$`+`$ $`2\left(V_{td}^{}V_{ts}V_{cd}^{}V_{cs}\right)\left[F_V^W(3;2)+F_V^H(3;2)\right],`$ (12)   
$`C_2(\mu _0)`$ $`=`$ $`\left(V_{td}^{}V_{ts}\right)^2F_S^H(3;3)+\left(V_{cd}^{}V_{cs}\right)^2F_S^H(2;2)`$ (13)   
$`+`$ $`2\left(V_{td}^{}V_{ts}V_{cd}^{}V_{cs}\right)F_S^H(3;2),`$ (14)   
$`C_3(\mu _0)`$ $`=`$ $`\left(V_{td}^{}V_{ts}\right)^2A_S^C,`$ (15)   
where the charm quark contributions have been kept. The superscripts $`W,H,C`$ denote the $`W^\pm ,H^\pm `$ and chargino contributions respectively, and

$`A_V^C`$ $`=`$ $`{\displaystyle \underset{i,j=1}{\overset{2}{}}}{\displaystyle \underset{k,l=1}{\overset{2}{}}}{\displaystyle \frac{1}{4}}G^{(3,k)i}G^{(3,k)j}G^{(3,l)i}G^{(3,l)j}Y_1(r_k,r_l,s_i,s_j),`$ (16)   
$`A_S^C`$ $`=`$ $`{\displaystyle \underset{i,j=1}{\overset{2}{}}}{\displaystyle \underset{k,l=1}{\overset{2}{}}}H^{(3,k)i}G^{(3,k)j}G^{(3,l)i}H^{(3,l)j}Y_2(r_k,r_l,s_i,s_j).`$ (17)   
Here $`G^{(3,k)i}`$ and $`H^{(3,k)i}`$ are the couplings of $`k`$th stop and $`i`$th chargino with left-handed and right-handed quarks, respectively :

$`G^{(3,k)i}`$ $`=`$ $`\sqrt{2}C_{R1i}^{}S_{tk1}{\displaystyle \frac{C_{R2i}^{}S_{tk2}}{\mathrm{sin}\beta }}{\displaystyle \frac{m_t}{M_W}},`$ (18)   
$`H^{(3,k)i}`$ $`=`$ $`{\displaystyle \frac{C_{L2i}^{}S_{tk1}}{\mathrm{cos}\beta }}{\displaystyle \frac{m_s}{M_W}},`$ (19)   
and $`C_{L,R}`$ and $`S_t`$ are unitary matrices that diagonalize the chargino and stop mass matrices . : $`C_R^{}M_\chi ^{}C_L=\mathrm{diag}(M_{\stackrel{~}{\chi _1}},M_{\stackrel{~}{\chi _2}})`$ and $`S_tM_{\stackrel{~}{t}}^2S_t^{}=\mathrm{diag}(M_{\stackrel{~}{t}_1}^2,M_{\stackrel{~}{t}_2}^2)`$. Explicit forms for functions $`Y_{1,2}`$ and $`F`$’s can be found in Ref. , and $`r_k=M_{\stackrel{~}{t}_k}^2/M_W^2`$ and $`s_i=M_{\stackrel{~}{\chi ^\pm }_i}/M_W^2`$. It should be noted that $`C_2(\mu _0)`$ was misidentified as $`C_3^H(\mu _0)`$ in Ref. . The gluino and neutralino contributions are negligible in our model. The Wilson coefficients at lower scales are obtained by renomalization group running. The relevant formulae with the NLO QCD corrections at $`\mu =2`$ GeV are given in Ref. . It is important to note that $`C_1(\mu _0)`$ and $`C_2(\mu _0)`$ are real relative to the SM contribution in our model. On the other hand, the chargino exchange contributions to $`C_3(\mu _0)`$ (namely $`A_S^C`$) are generically complex relative to the SM contributions, and can generate a new phase shift in the $`K^0\overline{K^0}`$ mixing relative to the SM value. This effect is in fact significant for large $`\mathrm{tan}\beta (1/\mathrm{cos}\beta )`$ , since $`C_3(\mu _0)`$ is proportional to $`(m_s/M_W\mathrm{cos}\beta )^2`$.

The CP violating parameter $`ϵ_K`$ can be calculated from

$$ϵ_K\frac{e^{i\pi /4}\mathrm{Im}M_{12}}{\sqrt{2}\mathrm{\Delta }M_K},$$
(20)   
where $`M_{12}`$ can be obtained from the $`\mathrm{\Delta }S=2`$ effective Hamiltonian through $`2M_KM_{12}=K^0|H_{\mathrm{eff}}^{\mathrm{\Delta }S=2}|\overline{K^0}`$. For $`\mathrm{\Delta }M_K`$, we use the experimental value $`\mathrm{\Delta }M_K=(3.489\pm 0.009)\times 10^{12}`$ MeV, instead of theoretical relation $`\mathrm{\Delta }M_K=2\mathrm{R}\mathrm{e}M_{12}`$, since the long distance contributions to $`M_{12}`$ is hard to calculate reliably unlike the $`\mathrm{\Delta }S=2`$ box diagrams. For the strange quark mass, we use the $`\overline{\mathrm{MS}}`$ mass at $`\mu =2`$ GeV scale : $`m_s(\mu =2\mathrm{G}\mathrm{e}\mathrm{V})=125`$ MeV. In Figs. 2 (a) and (b), we plot the results of scanning the MSSM parameter space : the correlations between $`ϵ_K/ϵ_K^{\mathrm{SM}}`$ and (a) $`\mathrm{tan}\beta `$ and (b) the lighter stop mass. We note that $`ϵ_K/ϵ_K^{\mathrm{SM}}`$ can be as large as $`1.4`$ for $`\delta _{KM}=90^{}`$ if $`\mathrm{tan}\beta `$ is small. This is a factor 2 larger deviation from the SM compared to the minimal supergravity case . The dependence on the lighter stop is close to the case of the minimal supergravity case, but we can have a larger deviations. Such deviation is reasonably close to the experimental value, and will affect the CKM phenomenology at a certain level.

In the MSSM with new CPV phases, there is an intriguing possibility that the observed CP violation in $`K_L\pi \pi `$ is fully due to the complex parameters $`\mu `$ and $`A_t`$ in the soft SUSY breaking terms which also break CP softly. This possibility was recently considered by Demir et al. . Their claim was that it was possible to generate $`ϵ_K`$ entirely from SUSY CPV phases for large $`\mathrm{tan}\beta 60`$ with certain choice of soft parameters <sup>§</sup><sup>§</sup>§Their choice of parameters leads to $`M_{\chi ^\pm }=80`$ GeV and $`M_{\stackrel{~}{t}}=85`$ GeV, which are very close to the recent lower limits set by LEP2 experiments.. In such a scenario, only $`\mathrm{Im}(A_S^C)`$ in Eq. (6) can contribute to $`ϵ_K`$, if we ignore a possible mixing between $`C_2`$ and $`C_3`$ under QCD renormalization. In actual numerical analysis we have included this effect using the results in Ref. . We repeated their calculations using the same set of parameters, but could not confirm their claim. For $`\delta _{KM}=0^{}`$, we found that the supersymmetric $`ϵ_K`$ is less than $`2\times 10^5`$, which is too small compared to the observed value : $`|ϵ_K|=(2.280\pm 0.019)\times 10^3`$ determined from $`K_{L,S}\pi ^+\pi ^{}`$ .

Let us give a simple estimate for fully supersymmetric $`ϵ_K`$, in which case only $`C_3(\mu _0)`$ develops imaginary part and can contribute to $`ϵ_K`$. For $`m_{\stackrel{~}{t}_1}m_{\chi ^\pm }M_W`$, we would get $`Y_2Y_2(1,1,1,1)=1/6`$, and

$`|G^{(3,k)i}|O(1),\mathrm{and}|H^{(3,k)i}|{\displaystyle \frac{m_s\mathrm{tan}\beta }{M_W}},`$

because any components of unitary matrices $`C_R`$ and $`S_t`$ are $`O(1)`$. Therefore $`\mathrm{Im}(A_S^C)O(10^3)`$. Now using

$$\mathrm{Im}(M_{12})=\frac{G_F^2M_W^2}{(2\pi )^2}f_K^2M_K\left(\frac{M_K}{m_s}\right)^2\frac{1}{24}B_3(\mu )\mathrm{Im}(C_3(\mu )),$$
(21)   
and Eq. (9), we get $`|ϵ_K|2\times 10^5`$.

7. In conclusion, we extended our previous studies of SUSY CPV phases to $`BX_sl^+l^{}`$ and $`ϵ_K`$. Our results can be summarized as follows :

* The branching ratio for $`BX_sl^+l^{}`$ can be enhanced upto $`85\%`$ compared to the SM prediction, and the correlation between $`\mathrm{Br}(BX_s\gamma )`$ and $`\mathrm{Br}(BX_sl^+l^{})`$ is distinctly different from the minimal supergravity scenario (CMSSM) (even with new CP violating phases)  in the presence of new CP violating phases in $`C_{7,8,9}`$ as demonstared in model-independent analysis by Kim, Ko and Lee .
* $`ϵ_K/ϵ_K^{SM}`$ can be as large as 1.4 for $`\delta _{KM}=90^{}`$. This is the extent to which the new phases in $`\mu `$ and $`A_t`$ can affect the construction of the unitarity triangle through $`ϵ_K`$.
* Fully supersymmetric CP violation is not possible even for large $`\mathrm{tan}\beta 60`$ and light enough chargino and stop, contrary to the claim made in Ref. . With real CKM matrix elements, we get very small $`|ϵ_K|O(10^5)`$, which is two orders of magnitude smaller than the experimental value.

Before closing this paper, we’d like to emphasize that all of our results are based on the assumption that there are no new CPV phases in the flavor changing sector. Once this assumption is relaxed, then gluino-mediated FCNC with additional new CPV phases may play important roles, and many of our results may change . Within our assumption, the results presented here and in Ref.  are conservative since we did not impose any conditions on the soft SUSY breaking terms except that the resulting mass spectra for chargino, stop and other sparticles satisfy the current lower bounds from LEP and Tevatron. More detailed analysis of phenomenological implications of our works on $`B_{d(s)}^0\overline{B_{d(s)}^0}`$ mixing, $`BX_{s(d)}\gamma ,X_{s(d)}l^+l^{},B_{s(d)}^0l^+l^{}`$ and their direct CP asymmetries will be presented elsewhere.

###### Acknowledgements.

The authors wich to thank G.C. Cho for clarifying $`O_2`$ and $`O_3`$ in Ref. , and A. Ali, A. Grant, A. Pilaftsis and O. Vives for useful communications. A part of this work was done while one of the authors (PK) was visiting Harvard University under the Distinguished Scholar Exchange Program of Korea Research Foundation. This work is supported in part by KOSEF Contract No. 971-0201-002-2, KOSEF through Center for Theoretical Physics at Seoul National University, Korea Rsearch Foundation Program 1998-015-D00054 (PK), and by KOSEF Postdoctoral Fellowship Program (SB).