Effects of sfermion mixing induced by RGE running in the minimal flavor violating CMSSM

Within the Constrained Minimal Supersymmetric Standard Model (CMSSM) with Minimal Flavor Violation (MFV) for scalar quarks, we study numerically the effects of intergenerational squark mixing on B-physics observables, electroweak precision observables (EWPO), and the Higgs-boson mass predictions. In models with universal soft terms at the GUT scale, squark mixing is generated through the Renormalization Group Equations (RGEs) running from the GUT scale to the electroweak scale due to presence of non-diagonal Yukawa matrices in the RGEs, e.g. due to the CKM matrix. Our numerical analysis is based on the code Spheno for the RGE running and full one-loop calculations, supplemented by further higher-order corrections, at the electroweak scale of the precision observables as included in the code FeynHiggs. Taking the CMSSM as a concrete “realistic” example, we find that the B-physics observables as well as the Higgs mass predictions do not receive sizable corrections. On the other hand, in our numerical analysis we observe that the EWPO such as the W boson mass can receive relevant corrections. Such contributions could in principle be used to place new bounds on the CMSSM parameter space. We extend our numerical analysis to the CMSSM extended with a mechanism to explain neutrino masses (CMSSM-seesaw I), which induces flavor violation in the scalar lepton sector. The effects of slepton mixing on the analyzed observables are found to be, in general, smaller than those of squark mixing, but in our numerical analysis reach the level of the current experimental uncertainty for the EWPO.


Introduction
Supersymmetric (SUSY) extensions of the Standard Model (SM) are broadly considered as the most motivated and promising New Physics (NP) theories beyond the SM. The solution of the hierarchy problem, the gauge coupling unification and the possibility of having a natural cold dark matter candidate, constitute the most convincing arguments in favor of SUSY.
Within the Minimal Supersymmetric Standard Model (MSSM) [1][2][3], flavor mixing can occur in both scalar quark and scalar lepton sector. Here the possible presence of soft SUSY-breaking parameters in the squark and slepton sector, which are off-diagonal in flavor space (mass parameters as well as trilinear couplings) are the most general way to introduce flavor mixing within the MSSM. This, however, yields many new sources of flavor and CP-violation, which potentially lead to large non-standard effects in flavor processes, in conflict with the experimental bounds.
The SM has been very successfully tested by low-energy flavor observables both from the kaon and the B d sectors. In particular, the two B factories have established that B d flavor and CP-violating processes are well described by the SM up to an accuracy of the ∼10 % level [4]. This immediately implies a tension between the solution of the hierarchy problem, calling for a NP scale at or below the TeV scale, and the explanation of the flavor physics data, requiring a multi-TeV NP scale, if the new flavor-violating couplings are of generic size.
An elegant way to simultaneously solve the above problems is provided by the Minimal Flavor Violation (MFV) hypothesis [5][6][7][8], where flavor and CP-violation in the quark sector are assumed to be entirely described by the CKM matrix, even in theories beyond the SM. For example in MSSM, the off-diagonality in the sfermion mass matrix reflects the misalignment (in flavor space) between fermions and sfermions mass matrices, which cannot be diagonalized simultaneously. This misalignment can be produced from various origins. For instance, off-diagonal sfermion mass matrix entries can be generated by Renormalization Group Equations (RGE) running. Going from a high energy scale, where no flavor violation is assumed, down to the electroweak (EW) scale, such entries can be generated due to the presence of non-diagonal Yukawa matrices in the RGEs. For instance, in the Constrained Minimal Supersymmetric Standard Model (CMSSM, see [9] and references therein), the RGE effects on non-diagonal sfermion soft SUSY-breaking parameters are affected only by non-diagonal elements on the Yukawa couplings and the trilinear terms which are taken as proportional to the Yukawas at the GUT scale. We choose the following form of the Yukawa matrices (working in the Super-CKM basis [10]): Hence, all flavor violation in the quark and squark sector is generated by the RGEs and controlled by the CKM matrix, i.e. the Yukawa couplings have a strong impact on the size of the induced off-diagonal entries in the squark mass matrices.
The situation is somewhat different in the slepton sector where neutrinos are strictly massless (in the SM and the MSSM). Consequently, there is no slepton mixing, which would induce Lepton Flavor Violation (LFV) in the charged sector, allowing not yet observed processes like l i → l j γ (i > j; l 3,2,1 = τ, μ, e) [11]. However, in the neutral sector, we have strong experimental evidence that shows that the neutrinos are massive and mix among themselves [12][13][14][15][16][17][18][19][20]. In order to incorporate this, one needs to go beyond the MSSM to introduce a mechanism that generates neutrino masses. The simplest way would be to introduce Dirac masses, leaving, however, the extreme smallness of the neutrino masses unexplained. To overcome this problem, typically a seesaw mechanism is used to generate neutrino masses, and the PMNS matrix plays the role of the CKM matrix in the lepton sector. Extending the MFV hypothesis for leptons [21] we can assume that the flavor mixing in the lepton and slepton sector is induced and controlled by the seesaw mechanism.
In many analyses of the CMSSM, or extensions such as the NUHM1 or NUHM2 (see [9] and references therein), the hypothesis of MFV has been used, and it has been assumed that the contributions coming from MFV are negligible for other observables as well; see, e.g., [29][30][31][32]. In this paper we will perform a numerical analysis to see whether this assumption is justified, and whether including these MFV effects could in principle lead to additional constraints on the CMSSM parameter space. However, we do not attempt to find analytical solutions to analyze this, as they become extremely involved in the presence of intergenerational mixing in SUSY models. In this respect we numerically evaluate in the CMSSM and in the CMSSM-seesaw I the following set of observables: B physics observables (BPO), in particular BR(B → X s γ ), BR(B s → μ + μ − ) and M B s , electroweak precision observables (EWPO), in particular M W and the effective weak leptonic mixing angle, sin 2 θ eff , as well as the masses of the neutral and charged Higgs bosons in the MSSM.
In order to perform our calculations, we use the code SPheno [33,34] to generate the CMSSM (containing also the type I seesaw) particle spectrum by running RGE from the GUT down to the EW scale. The effects of the CKM matrix in the RGE running on the mixing in the scalar fermion sector thus fully relies on the SPheno implementation. The particle spectrum was then handed over in the form of an SLHA file [10,35] to FeynHiggs [36][37][38][39][40][41] to calculate EWPO and Higgs-boson masses. The B physics observables were calculated by the BPHYSICS subroutine included in the SuFla code [42,43] (see also [44][45][46] for the improved version used here).
Our setup provides an evaluation of flavor-violating effects in "realistic" MFV models (where flavor violation enters only via RGE running) using state-of-the-art tools, compared to state-of-the-art limits, where the size of the effects will also be compared to future sensitivities. Effects that may appear negligible now might be non-negligible in the future. Furthermore, in the case of lepton-flavor violation [case (ii) above], we are not aware of any analysis of this type.
The paper is organized as follows: First we review the main features of the MSSM with sfermion flavor mixing in MFV in Sect. 2. The computational setup is given in Sect. 3. The numerical results are presented in Sect. 4, where first we discuss the effect of squarks mixing in the CMSSM. In a second step we analyze the effects of slepton mixing i.e. the CMSSM-seesaw I. Our conclusions can be found in Sect. 5.

Model setup
In this section we will first review the CMSSM and the concept of MFV. Subsequently, we will discuss the MSSM, its seesaw extension and parameterization of sfermion mixing at low energy.

The CMSSM and MFV
The MSSM is the simplest supersymmetric structure we can build from the SM particle content. The general setup for the soft SUSY-breaking parameters is given by [1][2][3] Here m 2Q and m 2L are 3 × 3 matrices in family space (with i, j being the generation indices) for the soft masses of the left-handed squarkq L and sleptonl L SU (2) doublets, respectively. m 2 u , m 2d , and m 2 e contain the soft masses for right-handed up-type squarksũ R , down-type squarks d R , and charged sleptonẽ R SU (2) singlets, respectively. A u , A d , and A l are the 3 × 3 matrices for the trilinear couplings for up-type squarks, down-type squarks, and charged slepton, respectively;m 1 andm 2 are the soft masses of the Higgs sector. In the last line M 1 , M 2 , and M 3 define the bino, wino, and gluino mass terms, respectively.
Within the constrained MSSM the soft SUSY-breaking parameters are assumed to be universal at the Grand Unification scale M GUT ∼ 2 × 10 16 GeV, There is a common mass for all the scalars, m 2 0 , a single gaugino mass, m 1/2 , and all the trilinear soft-breaking terms are directly proportional to the corresponding Yukawa couplings in the superpotential with a proportionality constant A 0 e iφ A , containing a potential non-trivial complex phase.
With the use of the Renormalization Group Equations (RGE) of the MSSM, one can obtain the SUSY spectrum at the EW scale. All the SUSY masses and mixings are then given as a function of m 2 0 , m 1/2 , A 0 , and tan β = v 2 /v 1 , the ratio of the two vacuum expectation values (see below). We require radiative symmetry breaking to fix |μ| and |Bμ| [47,48] with the tree-level Higgs potential.
By definition, this model fulfills the MFV hypothesis, since the only flavor-violating terms stem from the CKM matrix. The important point is that, even in a model with universal soft SUSY-breaking terms at some high energy scale as the CMSSM, some off-diagonality in the squark mass matrices appears at the EW scale. Working in the basis where the squarks are rotated parallel to the quarks, the so-called Super CKM (SCKM) basis, the squark mass matrices are not flavor diagonal at the EW scale. This is due to the fact that at M GUT there exist two non-trivial flavor structures, namely the two Yukawa matrices for the up and down quarks, which are not simultaneously diagonalizable. This implies that through RGE evolution some flavor mixing leaks into the sfermion mass matrices. In a general SUSY model the presence of new flavor structures in the soft SUSY-breaking terms would generate large flavor mixing in the sfermion mass matrices. However, in the CMSSM, which we are investigating here, the two Yukawa matrices are the only source of flavor change. As always in the SCKM basis, any off-diagonal entry in the sfermion mass matrices at the EW scale will be necessarily proportional to a product of Yukawa couplings. The RGEs for the soft SUSY-breaking terms are sets of linear equations, and, thus, to match the correct chirality of the coupling, Yukawa couplings or trilinear soft terms must enter the RGE in pairs. (The same holds for the CMSSM-seesaw I; see below.)

MSSM and its seesaw extension
One can write the most general SU (3) C × SU (2) L × U (1) Y gauge invariant and renormalizable superpotential as where L i represents the chiral multiplet of a SU (2) L doublet lepton, E c i a SU (2) L singlet charged lepton, H 1 and H 2 two Higgs doublets with opposite hypercharge. Similarly Q, U , and D represent chiral multiplets of quarks of a SU (2) L doublet and two singlets with different U (1) Y charges. Three generations of leptons and quarks are assumed and thus the subscripts i and j run over 1 to 3. The symbol αβ is an anti-symmetric tensor with 12 = 1.
In order to provide an explanation for the (small) neutrino masses, the MSSM can be extended by the type-I seesaw mechanism [22][23][24][25][26][27][28]. The superpotential for CMSSMseesaw I can be written as where W MSSM is given in Eq. (4) and N c i is the additional superfield that contains the three right-handed neutrinos, ν Ri , and their scalar partners,ν Ri . M i j N denotes the 3×3 Majorana mass matrix for heavy right-handed neutrino. The full set of soft SUSY-breaking terms is given by By the seesaw mechanism three of the neutral fields acquire heavy masses and decouple at high energy scale that we will denote M N ; below this scale the effective theory contains the MSSM plus an operator that provides masses to the neutrinos.
This framework naturally explains neutrino oscillations in agreement with experimental data [12][13][14][15][16][17][18][19][20]. At the electroweak scale an effective Majorana mass matrix for light neutrinos, arises from Dirac neutrino Yukawa Y ν (which can be assumed of the same order as the charged-lepton and quark Yukawas), and heavy Majorana masses M N . The smallness of the neutrino masses implies that the scale M N is very high, O(10 14 GeV). From Eqs. (5) and (6) we can observe that one can choose a basis such that the Yukawa coupling matrix, Y i j l , and the mass matrix of the right-handed neutrinos, M i j N , are diagonalized as Y δ l and M δ R , respectively. In this case the neutrino Yukawa couplings Y i j ν are not generally diagonal, giving rise to LFV. Here it is important to note that the lepton-flavor conservation is not a consequence of the SM gauge symmetry, even in the absence of the righthanded neutrinos. Consequently, slepton-mass terms can violate the lepton-flavor conservation in a manner consistent with the gauge symmetry. Thus the scale of LFV can be identified with the EW scale, much lower than the righthanded neutrino scale M N , leading to potentially observable rates.
In the SM augmented by right-handed neutrinos, the flavor-violating processes such as μ → eγ , τ → μγ etc., whose rates are proportional to inverse powers of M δ R , would be highly suppressed with such a large M N scale, and hence are far beyond current experimental bounds. However, in SUSY theories, the neutrino Dirac couplings Y ν enter in the RGEs of the soft SUSY-breaking sneutrino and slepton masses, generating LFV. In the basis where the chargedlepton masses Y is diagonal, the soft slepton-mass matrix acquires corrections that contain off-diagonal contributions from the RGE running from M GUT down to the Majorana mass scale M N , of the following form (in the leading-log approximation) [49]: Consequently, even if the soft scalar masses were universal at the unification scale, quantum corrections between the GUT scale and the seesaw scale M N would modify this structure via renormalization-group running, which generates offdiagonal contributions [50][51][52][53][54][55] at M N in a basis such that Y is diagonal. Below this scale, the off-diagonal contributions remain almost unchanged. Therefore the seesaw mechanism induces non-trivial values for slepton δ F AB i j resulting in a prediction for LFV decays l i → l j γ , (i > j) that can be much larger than the non-SUSY case. These rates depend on the structure of Y ν at a seesaw scale M N in a basis where Y l and M N are diagonal. By using the approach of [55] a general form of Y ν containing all neutrino experimental information can be written as where R is a general orthogonal matrix and m δ ν denotes the diagonalized neutrino mass matrix. In this basis the matrix U can be identified with the U PMNS matrix obtained as In order to find values for the slepton generation mixing parameters we need a specific form of the product Y † ν Y ν as shown in Eq. (9). The simple consideration of direct hierarchical neutrinos with a common scale for right-handed neutrinos provides a representative reference value. In this case using Eq. (10) we find Here M R is the common mass assigned to the ν R . In the conditions considered here, LFV effects are independent of the matrix R.
For the numerical analysis the values of the Yukawa couplings etc. have to be set to yield values in agreement with the experimental data for neutrino masses and mixings. In our computation, by considering a normal hierarchy among the neutrino masses, we fix m ν 3 ∼ m 2 atm ∼ 0.05 eV and require m ν 2 /m ν 3 = 0.17, m ν 2 ∼ 100 · m ν 1 consistent with the measured values of m 2 sol and m 2 atm [56]. The matrix U is identified with U PMNS with the CP-phases set to zero and neutrino mixing angles set to the center of their experimental values.
One can observe that m eff remains unchanged by consistent changes on the scales of M N and Y ν . This is no longer correct for the off-diagonal entries in the slepton-mass matrices (parameterized by slepton δ F AB i j , see the next subsection). These quantities have a quadratic dependence on Y ν and a logarithmic dependence on M N ; see Eq. (9). Therefore larger values of M N imply larger LFV effects. By setting M N = 10 14 GeV, the largest values of Y ν are of about 0.29, this implies an important restriction on the parameters space arising from the BR(μ → eγ ) as will be discussed in Sects. 3 and 4. An example of models with almost degenerate ν R can be found in [50]. For our numerical analysis we tested several scenarios and we found that the one defined here is the simplest and also the one with larger LFV prediction.

Scalar fermion sector with flavor mixing
In this section we give a brief description about how we parameterize flavor mixing at the EW scale. We are using the same notation as in [44][45][46]57,58]. However, while in this section we give a general description, in our analysis below, contrary to our previous analyses [57], this time we concentrate on the origin of the flavor mixing as discussed in the previous sections.
The most general hypothesis for flavor mixing assumes a mass matrix that is not diagonal in flavor space, both for squarks and sleptons. In the squarks sector and charged slepton sector we have 6 × 6 mass matrices, based on the corresponding six electroweak interaction eigenstates,Ũ L ,R with U = u, c, t for up-type squarks,D L ,R with D = d, s, b for down-type squarks andL L ,R with L = e, μ, τ for charged sleptons. For the sneutrinos we have a 3×3 mass matrix, since within the MSSM even with type I seesaw (right-handed neutrinos decouple below their respective mass scale) we have only three electroweak interaction eigenstates,ν L with ν = ν e , ν μ , ν τ .
The non-diagonal entries in this 6 × 6 general matrix for sfermions can be described in terms of a set of dimensionless where F identifies the sfermion type, L , R refer to the "left-" and "right-handed" SUSY partners of the corresponding fermionic degrees of freedom, and i, j indices run over the three generations. (Non-zero values for the δ F AB i j are generated via the processes discussed in the previous subsections.) One usually writes the 6 × 6 non-diagonal mass matrices, M 2 u and M 2d , referred to the Super-CKM basis, being ordered, respectively, as referred to the Super-PMNS basis, being ordered as (ẽ L ,μ L ,τ L ,ẽ R ,μ R ,τ R ), and write them in terms of left-and right-handed blocks where where . It should be noted that the non-diagonality in flavor comes exclusively from the soft SUSY-breaking parameters, that could be non-vanishing for i = j, namely: the masses mQ i j and mL i j for the sfermion SU (2) doublets, the masses m 2Ũ , mẼ i j for the sfermion SU (2) singlets and the trilinear couplings, A f i j . In the sneutrino sector there is, correspondingly, a oneblock 3×3 mass matrix, that is referred to the (ν eL ,ν μL ,ν τ L ) electroweak interaction basis: where It is important to note that due to SU (2) L gauge invariance the same soft masses mQ i j enter in both up-type and down-type squarks mass matrices similarly mL i j enter in both the slepton and the sneutrino L L mass matrices. The soft SUSY-breaking parameters for the up-type squarks differ from corresponding ones for down-type squarks by a rotation with CKM matrix. The same would hold for sleptons i.e. the soft SUSY-breaking parameters of the sneutrinos would differ from the corresponding ones for charged sleptons by a rotation with the PMNS matrix. However, taking the neutrino masses and oscillations into account in the SM leads to LFV effects that are extremely small. For instance, in μ → eγ they are of O(10 −47 ) in the case of Dirac neutrinos with mass around 1 eV and maximal mixing [59][60][61][62], and of O(10 −40 ) in the case of Majorana neutrinos [59,62]. Consequently we do not expect large effects from the inclusion of neutrino mass effects here and neglect a rotation with the PMNS matrix. The sfermion mass matrices in terms of the δ F AB i j are given as In all this work, for simplicity, we are assuming that all δ F AB i j parameters are real and, therefore, hermiticity of M 2Q , . The next step is to rotate the squark states from the Super-CKM basis,q L ,R , to the physical basis. If we set the order in the Super-CKM basis as above, , and in the physical basis asũ 1,...,6 andd 1,...,6 , respectively, these last rotations are given by two 6 × 6 matrices, Rũ and Rd , yielding the diagonal mass-squared matrices for squarks as follows: Similarly we need to rotate the sleptons and sneutrinos from the electroweak interaction basis to the physical mass eigenstate basis, with Rl and Rν being the respective 6 × 6 and 3 × 3 unitary rotating matrices that yield the diagonal mass-squared matrices as follows:

Computational setup
Here we briefly describe our numerical setup. We first give some details on the running from the GUT to the EW scale, and subsequently describe the calculations of the observables evaluated at the EW scale.

From the GUT scale to the EW scale
The SUSY spectra have been generated with the code SPheno 3.2.4 [33,34] (for the CMSSM and the CMSSMseesaw I). We defined the SLHA [10,35] file at the GUT scale. In a first step within SPheno, gauge and Yukawa couplings at M Z scale are calculated using tree-level formulas. Fermion masses, the Z boson pole mass, the fine-structure constant α, the Fermi constant G F , and the strong coupling constant α s (M Z ) are used as input parameters. The gauge and Yukawa couplings, calculated at M Z , are then used as input for the one-loop RGEs to obtain the corresponding values at the GUT scale which is calculated from the requirement that g 1 = g 2 (where g 1,2 denote the gauge couplings of the U (1) and SU (2), respectively). The CMSSM boundary conditions are then applied to the complete set of two-loop RGEs and are evolved to the EW scale. At this point the SM and SUSY radiative corrections are applied to the gauge and Yukawa couplings, and the two-loop RGEs are again evolved to GUT scale. After applying the CMSSM boundary conditions again the two-loop RGEs are run down to EW scale to get SUSY spectrum. This procedure is iterated until the required precision is achieved. As stressed above, for the effects of the CKM matrix on the sfermion mixing we fully rely on Spheno. The output is then written in the form of an SLHA file, which is used as input to calculate low-energy observables discussed below.
For the CMSSM-seesaw I a similar procedure is applied, where the neutrino related input parameters are included in the respective SLHA input blocks (see [10,35] for details), the relevant numerical values are given in Sect. 2.2. For our scans of the CMSSM-seesaw I parameter space we use SPheno 3.2.4 [33,34] with the model "seesaw type-I". The value for Y ν is implemented as explained in Sect. 2.2, adjusting the matrix elements such that neutrino experimental parameters achieve the desired results after RGEs. The predictions for BR(l i → l j γ ) are also obtained with SPheno 3.2.4, see the discussion in Sect. 4.2. We checked that the use of this code produces results similar to the ones obtained by our private codes used in [50].

Calculations at the EW scale
Here we briefly review the various observables that we compute at the EW scale, either taking the non-zero δ F AB i j into account, or setting them to zero.

The MSSM Higgs sector
The MSSM Higgs sector consist of two Higgs doublets and predicts five physical Higgs bosons, the light and heavy CP-even h and H , the CP-odd A, and the charged Higgs boson, H ± . At tree level the Higgs sector is described with the help of two parameters: the mass of the A boson, M A , and tan β = v 2 /v 1 , the ratio of the two vacuum expectation values. The tree-level relations receive large higher-order corrections; see, e.g., [63,64]  Effects of sfermion mixing in the MSSM Higgs sector has already been calculated in a model independent way in the scalar quark sector [44][45][46]69], as well as independently in [70]. They have also been calculated in the scalar lepton sector in [57]. In both cases there are sizable corrections to the Higgs-boson masses, specially to the charged Higgsboson mass M H ± , assuming general NMFV in the squark and slepton sector.
In the Feynman diagrammatic approach that we are following here, the higher-order corrected CP-even Higgs-boson masses are derived by finding the poles of the (h, H )propagator matrix. The inverse of this matrix is given by Determining the poles of the matrix Higgs in Eq. (36) is equivalent to solving the equation Similarly, in the case of the charged Higgs sector, the corrected Higgs mass is derived by the position of the pole in the charged Higgs propagator, which is defined by The flavor-violating parameters enter into the one-loop prediction of the various (renormalized) Higgs-boson selfenergies, where details can be found in [44][45][46]57]. Numerically the results have been obtained using the code FeynHiggs [36][37][38][39][40][41], which contains the complete set of one-loop corrections from (flavor-violating) squark and slepton contributions (based on [44,45,57,69]). Those are supplemented with leading and sub-leading two-loop corrections as well as a resummation of leading and sub-leading logarithmic contributions from the t/t sector, all evaluated in the flavor conserving MSSM.

Electroweak precision observables
EWPO that are known with an accuracy at the per-mille level or better have the potential to allow for a discrimination between quantum effects of the SM and SUSY models; see [71] for a review. Examples are the W -boson mass M W and the Z -boson observables, such as the effective leptonic weak mixing angle sin 2 θ eff , whose present experimental uncertainties are [72] δ M exp,today W ∼ 15 MeV, δ sin 2 θ exp,today eff The experimental uncertainty will further be reduced [73,74] at the ILC and at the GigaZ option of the ILC, respectively. An even higher precision could be expected from the FCC-ee; see, e.g., [75].
The W -boson mass can be evaluated from where α is the fine-structure constant and G μ the Fermi constant. This relation arises from comparing the prediction for muon decay with the experimentally precisely known Fermi constant. The one-loop contributions to r can be written as where α is the shift in the fine-structure constant due to the light fermions of the SM, α ∝ log(M Z /m f ), and ρ is the leading contribution to the ρ parameter [76] from (certain) fermion and sfermion loops (see below). The remainder part ( r ) rem contains in particular the contributions from the Higgs sector. The effective leptonic weak mixing angle at the Z -boson resonance, sin 2 θ eff , is defined through the vector and axialvector couplings (g V and g A ) of leptons ( ) to the Z boson, measured at the Z -boson pole. If this vertex is written as Loop corrections enter through higher-order contributions to g V and g A . Both of these (pseudo-)observables are affected by shifts in the quantity ρ according to The quantity ρ is defined by the relation with the unrenormalized transverse parts of the Z -and Wboson self-energies at zero momentum, T Z ,W (0). It represents the leading universal corrections to the electroweak precision observables induced by mass splitting between partners in isospin doublets [76]. Consequently, it is sensitive to the mass-splitting effects induced by flavor mixing. The effects from flavor violation in the squark and slepton sector, entering via ρ have been evaluated in [57,69] and included in FeynHiggs. In particular, in [69] it has been shown that for the squark contributions ρ constitutes an excellent approximation to r . We use FeynHiggs for our numerical evaluation.
Concerning the expected effects in ρ some more detailed comments are in order. Within the SM the corrections to ρ stem from the splitting in one SU (2) doublet. Due to the mixing of various scalar fermion states the picture is slightly more involved in the MSSM. In MSSM without flavor violation the well-known results for the third generation squark contribution to ρ (without flavor mixing) can be written as with In the absence of intergenerational mixing there are only 2×2 mixing matrices to be taken into account, here parametrized by θt (θb) in the scalar top (bottom) case. Here one can see that squarks do not need to be the SU (2) partners to give a contribution to ρ. In particular the first two terms of Eq. (46) describe contributions from the same type (up type or down type) of scalar quarks. Going from this simple case to the one with generation mixing, one finds a contribution from all three generations, including two 6 × 6 mixing matrices (which are difficult to analyze analytically). For the sake of completeness, the two gauge boson self-energies are then given by (see also [69]) Here Rũ and Rd are the 6 × 6 rotation matrices for the upand down-type squarks, respectively; see Eq. (28). The finite part of the one point integral function is given by Here it is important to note that the corrections will come, as in Eq. (46), from states connected via SU (2) as well as from "same flavor" contributions stemming from the Z boson self-energy; see Eq. (45). Larger splitting between "same flavor" states due to the intergenerational mixing thus leads to the expectation of increasing contributions to ρ from flavor-violation effects.   included again what are known to be the dominant contributions to these three types of diagrams in scenarios with non-minimal flavor violation (for a review see, for instance, [79]): gluino contributions to box diagrams, chargino contributions to box and Z -penguin diagrams, and chargino and gluino contributions to double φ-penguin diagrams. More details about the calculations employed can be found in [44][45][46]. We perform our numerical calculation with the BPHYSICS subroutine taken from the SuFla code [42,43] (with some additions and improvements as detailed in [44][45][46]), which has been implemented as a subroutine into (a private version of) FeynHiggs. The present experimental status and SM prediction of these observables is given in Table 1 [80-87].

B-physics observables
with μ > 0. Primarily we are not interested in the absolute values for all these observables but the effects that come from flavor violation within the MFV framework, i.e. the effect from the off-diagonal entries in the sfermion mass matrices. We first calculate the low-energy observables by setting all δ F AB We then evaluate the "pure MFV effects", whereas δ U L R and δ DL R decrease with m 0 . This behavior can be understood for the RGEs of the non-diagonal SUSY breaking parameters (see, e.g., [88,89]), the δ QL L are defined as ratios of off-diagonal soft terms that grow with m 2 0 over diagonal soft masses that also grow with m 0 . However, the δ U L R and the δ DL R arise from the ratio of the RGE generated off-diagonal trilinear terms, which depend on the value of A 0 , which is considered fixed in our case, over diagonal soft masses growing with m 0 . As discussed above, these δ F AB i j = 0 are often neglected in phenomenological analyses of the CMSSM (see, e.g., [29][30][31][32]). We also emphasize that these effects are purely due to the presence of the CKM matrix on the RGEs; their contribution will vanish when the mixing of the two first generation with the third generation is neglected (as we have checked numerically).
In Figs. 4, 5 and 6 we analyze the effects of the nonzero δ F AB i j on the EWPO ρ MFV , M MFV W and sin 2 θ MFV eff , respectively. Here the same pattern is reflected for the EWPO, i.e. by increasing the value of tan β or A 0 , we find larger contributions to the EWPO. In particular, one can observe a nondecoupling effect for large values of m 0 . Larger soft SUSYbreaking parameters with the non-zero values in particular of δ QL L 23 , see above, lead to an enhanced splitting in masses belonging to an SU (2) doublet, and thus to an enhanced contribution to the ρ-parameter. The corresponding effects on M W and sin 2 θ eff , for m 0 3 TeV, exhibit corrections that are several times larger than the current experimental accuracy (whereas the SUSY corrections with all δ F AB i j = 0 decouple and go to zero). Consequently, including the non-zero values of the δ F AB i j and correctly taking these corrections into account, would yield an upper limit on m 0 , which in the known analyses so far is unconstrained from above [29][30][31][32]. A more detailed analysis within the CMSSM will be needed to determine the real upper bound on m 0 , which, however, is beyond the scope of this paper.
In order to gain more insight about the source of the large corrections to ρ (and thus to the EWPO), we show in Fig. 7 several relative mass (square) differences, (m 2 2 − m 2 1 )/(m 2 2 + m 2 1 ) in the m 0 -m 1/2 plane for fixed A 0 = 0 and tan β = 45. The left plot shows the mass difference for the two most stop-like squarks (i.e. in the limit of zero intergenerational mixing they coincide with the two scalar tops). The right plot shows the relative mass difference for the lightest most stoplike and most sbottom-like squark. (These results are simply the Spheno output in our scenario.) In both cases one can see that the relative mass differences increase (controlled by the non-zero δ F AB i j induced by the CKM matrix in the RGE running) in a fashion similar to the δ QL L discussed above, i.e. in particular for m 0 > m 1/2 > 1 TeV. These increasing mass differences lead (together with contributions from the mixing matrices) to the observed increase of ρ as in Fig. 4.
Our findings can be briefly compared to the existing literature. The EWPO in the context of flavor violation were evaluated first in [69], where correspondingly large corrections were found for large δ QL L 23 (in fact, that was the only parameter dependence analyzed in that paper, and only the mixing between the second and third generation of squarks was taken into account). Subsequently, the EWPO were also  [90]. The numerical analysis, however, was restricted to a degenerate and fixed SUSY mass scale. Correspondingly, no large effects with increasing SUSY mass scales were analyzed and only relative small corrections were found. Due to the different numerical setup, however, there is no contradiction with our results for ρ.
In Fig. 8  will increase as the scale M N increases but the parameter space excluded by the BR(μ → eγ ) bound will also increase. For example, by increasing M N by an order of magnitude, the largest entries in the matrix Y ν will become of O(1) and the bound on BR(μ → eγ ) will only be satisfied if m 0 ≈ 5 TeV.
Our numerical results in the CMSSM-seesaw I are shown in Figs. 9, 10, 11, 12, 13, 14 and 15. As in the CMSSM we present the results in the m 0 -m 1/2 plane for four combinations of tan β = 10, 45 (upper and lower row) and A 0 = 0, −3000 GeV (left and right column), again capturing the "extreme" cases. We start presenting the three most relevant δ F AB

Conclusions
In this paper we have investigated the CMSSM and the CMSSM-seesaw I (i.e. the CMSSM augmented by righthanded neutrinos to produce the observed neutrino mass pat- tern via the seesaw type I mechanism) under the hypothesis of Minimal Flavor Violation (MFV, i.e. the only source for flavor violation is the CKM matrix and/or the PMNS matrix in the case of the CMSSM-seesaw I). In many phenomenological analyses of the CMSSM the effects of intergenerational mixing in the squark and/or slepton sector are neglected. However, such mixings are naturally induced, assuming no flavor violation at the GUT scale, by the RGE running from the GUT to the EW scale exactly due to the presence of the CKM and/or the PMNS matrix. In this sense the CMSSM and the CMSSM-seesaw I represent two simple "realistic" GUT based models, in which flavor violation in induced solely by RGE running. The spectra of the CMSSM and CMSSM-seesaw I have been numerically evaluated with the help of the program SPheno by taking the GUT scale input run down via the appropriate RGEs to the EW scale.
We have evaluated the predictions for B-physics observables, MSSM Higgs-boson masses, electroweak precision observables in the CMSSM and CMSSM-seesaw I. In order to numerically analyze the effects of neglecting intergenerational mixing these observables have been evaluated with the full spectrum at the EW scale, as well as with the spectrum, but with all intergenerational mixing set numerically (artificially) to zero (as has been done in many phenomenological analyses). We did not attempt an analytical evaluation of the flavor-violation effects, as they become very involved in the case intergenerational sfermion mixing in SUSY models. The numerical difference in the various observables indicates the possible size of the effects neglected in those analyses. In this way it can be checked whether neglecting those mixing effects is a justified approximation.
Within the CMSSM we have taken a fixed grid of A 0 and tan β, while scanning the m 0 -m 1/2 plane. We found that the value of δ F AB i j increases with the increase of the A 0 or tan β values. The Higgs-boson masses receive corrections below current and future experimental uncertainties, where the shifts in M H ± were found largest at the level of O(100 MeV). Similarly for the B-physics observables the induced effects are at least one order of magnitude smaller than the current experimental uncertainty. For those two groups of observables the approximation of neglecting intergenerational mixing explicitly is a viable option.
The picture changes for the electroweak precision observables. We find that the masses of the squarks grow with m 0 , and so do the mixing terms, inducing a splitting between masses in an SU (2) doublet, leading (numerically) to a nondecoupling effect. For m 0 3 TeV the effects induced in M W and sin 2 θ eff are found to be several times larger than the current experimental uncertainties and could shift the CMSSM prediction outside the allowed experimental range. In this way, taking the intergenerational mixing into account could in principle set bounds on m 0 not present in recent phenomenological analyses. By investigating numerically squark mass differences, we have shown that this behavior can be traced back to the non-decoupling effects in the scalar quark mass matrices, provided by Spheno when taking into account the CKM matrix in the RGE running. However, we would like to point out that this bound only holds because of the particularly simple structure of the CMSSM and cannot be extended easily to other, more complicated model frameworks.
Going to the CMSSM-seesaw I the numerical results depend on the concrete model definition. We have chosen a set of parameter that reproduces correctly the observed neutrino data and simultaneously induces large LFV effects and induces relatively large corrections to the calculated observables. Consequently, parts of the parameter space are excluded by the experimental bounds on BR(μ → eγ ). Concerning the precision observables we find that B-physics observables are not affected, we also find that the additional effects induced by slepton-flavor violation on the Higgsboson masses are negligible. Again the EWPO are found to show the largest impact, where for M W numerical effects at the same level as the current experimental accuracy have been observed for very large values of m 0 . As above, we would like to point out that these effects are due to the relatively simple structure of the CMSSM-seesaw I.
To summarize: we have numerically analyzed two "realistic" GUT based models in which flavor violation in solely induced by the CKM matrix via RGE running (as evaluated using the Spheno code). We find that artificially setting all flavor-violating terms to zero in the CMSSM and CMSSMseesaw I is an acceptable approximation for B-physics observables, Higgs-boson masses (evaluated using a private version of FeynHiggs). However, in the electroweak precision observables (also evaluated with FeynHiggs) in our numerical analysis we find larger effects in the CMSSM and CMSSM-seesaw I. The numerical contributions are larger than the current experimental accuracy in M W and sin 2 θ eff . Taking those effects correctly into account could in principle place new bounds on m 0 not present in recent phenomenological analyses.