Non-minimal flavour violation in A4 × SU(5) SUSY GUTs with smuon assisted dark matter

We study CP-conserving non-minimal flavour violation in A4 × SU(5) inspired Supersymmetric Grand Unified Theories (GUTs), focussing on the regions of parameter space where dark matter is successfully accommodated due to a light right-handed smuon a few GeV heavier than the lightest neutralino. In this region of parameter space we find that some of the flavour-violating parameters are constrained by the requirement of the dark matter relic density, due to the delicate interplay between the smuon and neutralino masses. By scanning over GUT scale flavour violating parameters, constrained by low-energy quark and lepton flavour violating observables, we find a striking difference in the results in which individual parameters are varied to those where multiple parameters are varied simultaneously, where the latter relaxes the constraints on flavour violating parameters due to cancellations and/or correlations. Since charged lepton-flavour violation provides the strongest constraints within a GUT framework, due to relations between quark and lepton flavour violation, we examine in detail a prominent correlation between some of the flavour violating parameters at the GUT scale consistent with the stringent lepton flavour violating process μ → eγ. We also examine the relation between GUT scale and low scale flavour violating parameters, for both quarks and leptons, and show how the usual expectations may be violated due to the correlations when multiple parameters are varied simultaneously.


Introduction
Despite the absence of any direct experimental evidence, supersymmetric extensions continue to provide attractive solutions to the shortcomings of the Standard Model (SM). Phenomenologically, they provide viable dark matter candidates and can readily accommodate seesaw explanations for neutrino masses. From a more theoretical point of view, they cure the hierarchy problem related to the Higgs mass and lead to a more precise gauge-coupling unification as compared to the Standard Model (SM). This last point can be seen as a hint towards Grand Unified Theories (GUTs), where gauge couplings, certain soft-breaking parameters, and representations of SM fields are unified.
The fact that the predicted superpartners have not been observed so far, may to some extent be moderated by the argument that current direct searches rely on specific assumptions, for example when all scalar masses are assumed to be degenerate at the Grand Unification scale. Moreover, as superpartner mass limits increase, the assumption of the Minimal Flavour Violation (MFV) paradigm postulating that the flavour structure of the theory is the same as in the SM, i.e. all flavour-violating interactions are related to the CKM-and PMNS-matrices only, may be relaxed without violating experimental limits. Relaxing this assumption and allowing for additional sources of flavour violation JHEP03(2019)067 leads to a modification of the assumed decay patterns, for example in squark searches. As a further consequence, the obtained mass limits may be considerably weakened [1][2][3]. In addition, it appears that a considerable region of the parameter space of the TeV-scale Minimal Supersymmetric Standard Model (MSSM) can accomodate such Non-Minimal Flavour Violation (NMFV) in the squark sector with respect to current experimental and theoretical constraints [4][5][6].
The link between NMFV terms at the TeV scale and the GUT scale is interesting from both the phenomenological and model-building point of view. First, although they may be numerically rather different, flavour violating interactions in the squark and slepton sectors are linked if NMFV is implemented in a unification framework at the high-scale. The same source of flavour violation may therefore be challenged by experimental data from both sectors. Second, gaining information on the flavour structure of a new physics framework at the TeV scale may give valuable hints towards Grand Unification and related flavour symmetries. The present paper is a first step in this direction.
In a recent paper by one of the authors [47], such a scenario was discussed in the framework of an SU(5) GUT combined with an A 4 family symmetry. 1 The idea was that the three 5 representations form a single triplet of the family symmetry with a unified soft mass m F , while the three 10 representations are singlets with independent soft masses m T 1 , m T 2 , m T 3 . Assuming MFV, it was shown that in order to account for the muon anomalous magnetic moment (g − 2) µ , dark matter and LHC data, non-universal gaugino masses M i (i = 1, 2, 3) at the high scale are required in the framework of the MSSM.
We focussed on a region of parameter space that has not been studied in detail before characterised by low higgsino mass µ ≈ −300 GeV, as required by the (g − 2) µ . The latter also required a right-handed smuonμ R with a mass around 100 GeV, and a neutralinõ χ 0 1 several GeV lighter which allows successful dark matter. The LHC will be able to fully test this scenario with the upgraded luminosity via muon-dominated tri-and dilepton signatures resulting from higgsino dominatedχ ± 1χ 0 2 andχ + 1χ − 1 production, as well as direct smuon production searches in the above region of parameter space.
The above study [47] was clearly concerned with the implications of the flavoured GUT model for the superpartner spectrum consistent with (g − 2) µ and the Dark Matter relic JHEP03(2019)067 density. However, for simplicity, it was assumed that there was no flavour violation at the GUT scale, whereas it is well known that such flavour violation is expected in these models [48][49][50]. The goal of the present study is to extend this work to the NMFV framework by introducing off-diagonal squark and slepton mass-squared terms in the Lagrangian at the GUT scale, motivated by the analyses in refs. [48][49][50] which show that such flavour violation is generically expected. Here, we take a phenomenological (or model independent) approach, and simply introduce flavour violating terms at high energy to explore their effect on low energy observables. To this end we consider two MFV reference parameter points, one of which is inspired by the findings of ref. [47] and involves a very light smuon capable of accounting for (g − 2) µ , and the other one with a heavier smuon, harder to discover at the LHC, but not able to account for (g − 2) µ . In both cases, we then perturb around these points, switching on off-diagonal mass terms, consistently with SU(5), arising from A 4 breaking effects. We find interesting correlations between the flavour violating parameters at the GUT scale consistent with the stringent lepton flavour violating processes µ → eγ.
The present paper is organised as follows: in section 2, we first introduce the A 4 ×SU(5) supersymmetric model under consideration and present the implementation of NMFV therein. Section 3 then presents the method and relevant tools for our analysis. Our results are presented and discussed in section 4. Finally, our conclusions are given in section 5.

Non-minimal flavour violation in SUSY GUTs
In this section, we present the model under consideration in this paper, namely the Minimal Supersymmetric Standard Model (MSSM) based on SU(5) Grand Unification, including an A 4 flavour symmetry entailing NMFV terms already at the GUT scale.

SUSY-breaking in the MSSM
Although the exact breaking mechanism is not completely understood, it is well known that Supersymmetry (SUSY) must be broken to some degree. The associated SUSY-breaking Lagrangian contains all terms which do not necessarily respect SUSY but hold to the tenets of gauge invariance and renormalisability. Considering the Minimal Supersymmetric Standard Model (MSSM), the SUSY-breaking Lagrangian reads While the soft mass and trilinear parameters appearing in eq. (2.1) are assumed to be diagonal matrices in flavour space within the MFV framework, they may comprise non-diagonal entries when relaxing this hypothesis and considering a NMFV scenario. It

JHEP03(2019)067
should be noted that generic SUSY models do not possess any symmetry preventing large off-diagonal elements in soft-SUSY parameters. The soft mass matrices are defined in the Super-CKM (SCKM) basis 2 as: which are associated to the left-handed squarks, right-handed up-and down-type squarks, left-handed sleptons and sneutrinos, and right-handed sleptons, respectively. In addition, there are the trilinear coupling matrices: for the up-and down-type squarks and the sleptons. Detailed expressions for the diagonal elements of the matrices given in eqs. (2.2) and (2.3) can be found in ref. [51]. Note that the soft mass matrices in eq. (2.2) are symmetric due to the requirement for hermiticity. It is convenient to parametrize the off-diagonal, i.e. flavour violating, elements of the above matrices in a dimensionless manner by normalizing them to the respective diagonal entries of the sfermion mass matrices. In the SCKM basis, this leads to the following parameters [39];

JHEP03(2019)067
with v u and v d being the vacuum expectation values of the up-and down-type Higgs doublets, respectively. Note that these definitions hold at any scale. In the following, the scales of interest will be the GUT and TeV scales. Moreover, the situation where all off-diagonal NMFV parameters defined in eq. (2.4) vanish corresponds to a scenario with quite minimal flavour violation.

SU(5) unification
The MSSM realization under consideration is based on the gauge group SU(5), which is the smallest group containing the SM gauge group, and can accomodate its matter fields in the F = 5 and T = 10 representations according to where r, b, g denote the quark colours, and c denotes CP -conjugated fermions. The Higgs doublets H u and H d , which break the electroweak symmetry, may arise from SU(5) multiplets H 5 and H 5 , provided the colour triplet components are heavy. The SU(5) gauge group may be broken by an additional Higgs multiplet in the 24 representation developing a vacuum expectation value Including the above arguments into a supersymmetric framework, SU(5) symmetry provides relationships between the soft terms belonging to the supermultiplets within a given representation. For the MSSM under consideration here, we can write down the soft-breaking Lagrangian in terms of SU(5) fields: where F and T are the superpartner fields of F and T given in eq. (2.5). Comparing this with eq. (2.1) leads to the relations that hold at the GUT scale. Note that renormalization group evolution towards lower scales will spoil these relations.

The A 4 × SU(5) model
In addition to the SU(5) grand unification, we impose an A 4 (alternating group of order 4) flavour symmetry on the model under consideration. To this end, we unify the three families of F =5 = (d c , L) into the triplet of A 4 leading to a unified soft mass parameter m F for the three generations. 3 The three families of T i = 10 i = (Q, u c , e c ) i are singlets of A 4 , which means that the three generations may have independent soft mass parameters m T 1 , m T 2 , m T 3 [48,[52][53][54][55].
Through breaking the discrete symmetry just below the GUT scale, we can induce flavour violation in our soft parameters. We express this primordial flavour violation as the matrices M 2 T , M 2 F , A F T , and A T T analogously to eq. (2.2) in the flavour basis of A 4 , that is, before rotation to the SCKM: Note that the breaking of A 4 enforces off-diagonal elements of the M 2 T and A F T matrices in eq. (2.10) to be smaller than diagonal entries, and we also assume that off-diagonal elements in the other matrices are small. 4 This provides a theoretical motivation for smallbut-non-zero flavour violation in such a class of models. SU(5) gives the following relationships between the dimensionless NMFV parameters in the basis before rotation to the SCKM (as denoted by the subscript '0'):

JHEP03(2019)067
These four matrices parameterise the flavour violation in the A 4 × SU(5) setup studied here. Note that δ T , δ F and δ T T are necessarily symmetric whereas δ F T is not (see eqs. (2.8) and (2.10)) leading to a total of 15 NMFV parameters at the GUT scale. It is apparent that we have flavour violation at phenomenological scales from two distinct sources: the presence of off-diagonal elements in various coupling matrices at the GUT scale due to A 4 breaking, and further effects on the off-diagonal elements induced by RGE running. We do not consider a specific breaking mechanism or pattern for the discrete symmetry.

Setup and tools
The aim of the present study is to assess the impact of flavour violating parameters introduced at the GUT scale on low-energy physics. In order to work with a concrete framework, we focus on the model based on A 4 × SU(5) presented in section 2. We test the high-scale model against the Dark Matter (DM) relic density along with leptonic and hadronic flavour changing observables and the mass of the Higgs boson.

MFV reference points
In order to focus on the impact of NMFV terms in the Lagrangian of our model, we start by choosing suitable reference scenarios respecting the MFV paradigm. From previous work [47] it is apparent that successfully imposing the dark matter relic density as well as the anomalous magnetic moment of the muon on the A 4 × SU(5) framework requires rather specific parameter configurations. More precisely, the corresponding parameter points feature a physical spectrum where the "right-handed" smuon is light and almost massdegenerate with the lightest neutralino, which is bino-like. This allows to simultaneously satisfy the (g − 2) µ and relic density constraints [56,57]. For our study, we choose two MFV reference scenarios, which are summarized in table 1.
The first reference point of our choice corresponds to the scenario labelled 'BP4' in ref. [47]. For practical reasons, mainly due to including NMFV terms at the GUT scale, we do not make use of the same version of the spectrum generator SPheno. In consequence, effects from renormalization group running differ slightly, and we have adapted the input parameters of the original BP4 reference scenario to the ones given in table 1. However, note that, although there is a small deviation for the TeV scale parameters as compared to scenario BP4 of ref. [47], the phenomenological aspects of our reference scenario at the TeV scale are unaffected. Let us recall that the rather low smuon mass parameter, m T 2 = 200 GeV, which leads to the physical mass mμ R = 102.1 GeV, is required in order to satisfy simultaneously the (g − 2) µ and relic density constraints as discussed in ref. [47]. This particular choice for m T 2 is an assumption in this work.
While current limits on "right-handed" smuons still allow masses as low as about 100 GeV [58], this first scenario is going to be severely challenged by ongoing LHC searches. For this reason, we choose to include a second reference point which is inspired by the first one but features larger smuon and neutralino masses. This still allows satisfaction of the relic density constraint due to efficient co-annihilation and avoids LHC limits to be JHEP03(2019)067 published in the near future. Note that, however, the higher smuon mass mμ R ∼ 250 GeV does not resolve the tension between the Standard Model and the experimental value of (g − 2) µ . Let us emphasize that both reference scenarios capture the essential results of ref. [47], namely almost mass-degenerate "right-handed" smuon and bino-like neutralino, while all other MSSM states are essentially decoupled.

JHEP03(2019)067
Annihilation channel Relative contribution to Ωχ0 Table 2. Dominant annihilation channels contributing to the annihilation cross-section and the neutralino relic density in the two MFV reference scenarios of table 1.
In both reference scenarios, the required neutralino relic density is met thanks to efficient co-annihilation with the smuon and even smuon pair annihilation. All (co)annihilation contributions are summarized in table 2. Neutralino pair annihilation mainly proceeds through t-and u-channel smuon exchange, while smuon pair annihilation proceeds through neutralino t-or u-channel exchange. Moreover, the relative importance of the co-annihilation and smuon pair annihilation with respect to the neutralino pair annihilation is governed by the Boltzmann factor involving the mass difference of the two particles [59]. The smuon mass therefore plays a central role in this context. Considering NMFV, the offdiagonal elements of the matrices in eqs. (2.2) and (2.3) not only violate flavour but can in addition have a significant impact on the smuon mass and thus on the relic density.

Introducing NMFV
Starting from the two MFV reference points, we study the impact of flavour violating soft terms by perturbing around this scenario. Keeping the MFV parameters fixed at the values given in table 1, we perform a random scan on the flavour violating parameters introduced in eq. (2.11) at the GUT scale using flat prior distributions. In practice, we vary the NMFV parameters both independently and as part of a multi-dimensional scan over all parameters simultaneously. We subsequently study the impact of the constraints detailed in table 3.
More precisely, we require the Higgs-boson mass to be reasonably close to the observed value of about 125 GeV, where we account for a theory uncertainty of 2.5 GeV from the SPheno calculation. For the B s -meson oscillation, we consider the experimental value ∆M Bs = (17.757 ± 0.021) ps −1 [56] and add a theory uncertainty of 1.35 ps −1 [63] which dominates over the experimental error. For the neutralino relic density, we require that the lightest neutralino accounts for the totality of observed cold dark matter. The error given by the Planck collaboration is augmented in order to take into account the 1% accuracy of the theoretical calculation of the relic density by micrOMEGAs. For further details on experimental constraints we refer the reader to table 3 and the references therein.
Finally, note that although the reference scenarios defined in table 1 have in part been obtained considering the anomalous magnetic moment of the muon as a key observable [47], we do not take into account this constraint here. Since (g − 2) µ is a flavour-conserving   [66][67][68][69]. From the resulting code SPhenoMSSM we obtain through two-loop renormalization group equations for the soft-breaking parameters and the physical mass spectrum at the TeV scale, as well as numerical predictions for flavour observables listed in table 3. The neutralino relic density Ωχ0 1 h 2 is computed using the public package micrOMEGAs 4.3.5 [70][71][72]. Again, we have used SARAH to obtain the CalcHEP model files necessary to accomodate NFMV effects in the calculation. Our computational setup is summarized in figure 1. The mass spectrum obtained from SPhenoMSSM is handed to micrOMEGAs by making use of the SUSY Les Houches Accord 2 [73]. Note that, since the spin-independent scattering cross-sections related to direct dark matter detection given in table 1 are relatively low as compared to the corresponding experimental limits, we do not explicitly evaluate these cross-section in our NMFV scan. Before running SPheno, we first need to perform a CKM transformation to certain GUT scale matrices to comply with the basis that SPheno requires for the input parameters (see appendix A). Let us note that, for typical values of Yukawa parameters inserted into our MFV reference points, CKM matrix running between the GUT and TeV scales has been found to be negligible. We therefore assume that the CKM matrix is identical across all scales.
In the full multi-dimensional scan, the studied range for each parameter is set empirically to give reasonable computational efficiency as informed by one-dimensional scans over individual parameters. The obtained ranges have been increased slightly to be able to study whether correlations between the different NMFV parameters may result in larger allowed ranges as compared to the one-dimensional scan. The applied limiting values for each MFV scenario under consideration and for each NMFV parameter are given in table 4.
Already from the individual scans, it becomes apparent that for certain NMFV parameters, especially in the case of Scenario 1, small deviations from the MFV case can induce either a charged dark matter candidate (the smuon in this case) or tachyonic mass spectra.  Table 4. Ranges of the NMFV parameters defined at the GUT scale (see eq. (2.11)) for our multi-dimensional scans around the reference scenarios. Those parameters given as 0.0 have been switched off, since even small variations lead to tachyonic mass spectra and/or a charged LSP.

JHEP03(2019)067
We therefore set throughout the analysis of Scenario 1, and scan over the remaining 12 NMFV parameters according to table 4. This situation does not occur for Scenario 2, where we vary all 15 NMFV parameters.
Starting from parameters at the GUT scale, we test each point against the observables listed in table 3. Points which do not satisfy all the imposed constraints within the associated uncertainties are collected in the prior distribution only, while those which comply with all constraints are in addition recorded as part of the posterior distribution. In examining the latter, we obtain the allowed ranges for each of the NMFV parameters. In addition, by comparing the prior and posterior distributions, and taking into account posterior distributions based on a single constraint, we identify the most important constraints among those listed in  Table 5. Estimated allowed GUT scale flavour violation for both reference scenarios and impactful constraints ordered from the most to the least constraining. Where square brackets are shown open, we scan up to these values but even if we noticed some impact from constraints, it seems that the allowed region can be larger and extrapolation to concrete limits is not straightforward. * denotes parameters fixed to 0 in order to satisfy LSP and physical mass spectrum requirements. † stands for extrapolated ranges, meaning that the posterior does not actually drop to 0 but extrapolation to a limit is reasonable. A parameter that is constrained by 'prior' is limited by LSP and physical mass requirement.

Results and discussion
In this section we present the results of our analysis. Before coming to a more detailed discussion, we start by presenting the general aspects and the obtained limits on the NMFV parameters, presented in table 5. Ultimately, we perform two different kinds of scan on the parameter space: "individual" scans, where only a single δ is varied and all others are set to zero, and "simultaneous" scans where all of the NMFV parameters are varied at the same time according to the ranges given in table 4. From the multi-dimensional scan, we conclude that for the majority of the considered NMFV parameters, the most sensitive observables are the branching ratios of µ → eγ and µ → 3e, as well as the neutralino relic density Ωχ0 1 h 2 . As discussed in section 3, the impact of the relic density can be attributed to the small mass difference between the neutralino and the smuon, which depends strongly on the off-diagonal elements in the slepton mass matrix. Since both our reference scenario exhibit a relatively small value of (m T ) 22 , already rather tiny flavour violating elements can be excluded by current data.  Figure 2. Feynman diagrams that contribute to µ → eγ, dashed line represents a slepton and δ denotes mass insertion parameters. Photon should be taken to be emitted from any particle charged under QED.
Although the experimental limit is more stringent (by about a factor of two) for the decay µ → eγ, the µ → 3e decay has about the same constraining power and is in certain cases even the dominant constraint. This is explained as follows: the amplitude of µ → eγ is helicity-suppressed, and therefore contains a suppression factor m e /m µ . While this is also the case for µ → 3e diagrams related to those of µ → eγ, there are additional fourpoint diagrams, where the helicity suppression is lifted since no photon is involved. Despite the additional gauge coupling and the greater degree of loop suppression, these diagrams are numerically competitive to those of µ → eγ.
One can see that NMFV parameters mixing the first or second generation with the third generation are also mainly constrained by the decays µ → eγ and µ → 3e rather than by the corresponding τ decays such as τ → µγ or τ → eγ. This can be traced to the better experimental precision of the muonic decay measurements with respect to the analogous tau decays. Even though NMFV parameters mediating e − τ or µ − τ transitions lead to the dominant contributions of the tau decays, these parameters also can enter into the muon decay amplitudes. For example, if the µ → eγ process includes a stau in the loop, the corresponding amplitude is proportional to terms including products of the type (δ) 23 (δ) 13 . See figure 2 for a diagrammatic representation. Since the muon decay limits are stronger than the tau decay limits by four to five orders of magnitude, the e − τ and µ − τ mixing parameters are constrained by the e − µ processes first. We have explicitly checked this by artificially lowering the bounds on tau decays. In this case, the tau decay becomes the dominant constraint for the (δ) 13 and (δ) 23 parameters. Finally, we observe that the constraints coming from the hadronic sector, such as the decays B → X s γ or B s → µ + µ − , which are dominant in the case of NMFV in the squark sector alone [6], are not competitive as compared to the leptonic constraints mentioned above. This can be traced to the greater experimental precision of dedicated leptonic measurements compared to meson decays.

Scan around Scenario 1
We discuss here in detail the results obtained for the full NMFV scan around the reference scenario 1. The MFV parameters are fixed at the values given in table 1, while we scan over the NMFV parameters according to the ranges given in table 4, either individually (i.e. keeping all but one parameter to zero), or simultaneously. For each performed scan, we JHEP03(2019)067 record the prior distribution containing all points featuring a physical mass spectrum and neutralino dark matter candidate (see also figure 1) as well as the posterior distribution obtained when imposing either one or all constraints summarized in table 3. Figure 3 shows the obtained prior and posterior distributions for the NMFV parameter (δ F ) 12 . The viable region for this parameter with respect to the imposed constraints is much larger for the case of the simultaneous scan as compared to the individual scan result. Indeed, it is possible that more than one of the NMFV parameters enters the calculation of one or more observables. In such a case, interferences and/or cancellations between the contributions induced by different NMFV parameters can occur. As a consequence, they give rise to viable regions of parameter space that would not be fully explored when varying each parameter in isolation. This is seen quantitatively as a broadening of posterior distributions when comparing a simultaneous scan result against a histogram from an individual scan. Let us emphasize that this feature is present for several of the flavour violating parameters under consideration in our study. Figure 4, panel b) shows the action of a single observable, BR(µ → 3e), on the same parameter (δ F ) 12 for simulultaneous scan, and can thus be directly compared to figure 3. 5 Since the shape of the single-constraint posterior almost matches the posterior obtained imposing all constraints, we conclude that this parameter is mainly limited by the µ → 3e lepton decay bound. The µ → eγ observable is less important in this case (see table 5, corresponding posterior not shown).
Coming to the parameter (δ T ) 12 shown in figure 5 including all constraints, note that the obtained viable interval is again broadened when comparing the individual scan, leading to |(δ T ) 12 | 0.2×10 −2 , with the simultaneous one yielding the range |(δ T ) 12 | 1.6×10 −2 . For the same NMFV parameter (δ T ) 12 , we detail in figure 6 the effect of the three most important experimental constraints in the simultaneous scan. The µ → eγ constraint can be seen to admit the entirety of the scanned region of parameter space in the simultaneous scan, whereas it is far the most stringent constraint in the individual scan (see figure 5). Indeed, µ → 3e is the most constraining observable for this parameter when varied along with other flavour violating entries of mass matrices. In addition, figure 6 illustrates how the obtained shape of the posterior distribution is due to the influence of three experimental constraints imposed on the parameter space.
We now discuss the parameter (δ T ) 13 shown in figure 7. We can notice that it is constrained only by the neutralino relic density and that the flavour constraints have no effect. This gives insight on the unexpected shape of the posterior distribution: as we have seen for two examples above, other NMFV parameters are allowed under flavour constraints to shift significantly away from zero. This has a marked effect in reducing superpartner masses which are determined by diagonalising the mass-squared matrices from eq. (2.2). This applies in particular to the "right-handed" smuon mass, as the initial smallness of m T 2 means that small NMFV parameters can slightly lower the smuon mass. As a further consequence, the relic density is then reduced due to the smaller mass difference between smuon and neutralino, which increases the importance of co-annihilation and smuon pair annihilation. However, the smuon mass also is influenced by (δ T ) 13 , which by virtue of being unconstrained by flavour observables, may be non-zero. Moreover, this particular parameter increases the lightest smuon mass due to the specific hierarchies in the mass matrix. The smuon mass being decreased by other non-zero NMFV parameters, (δ T ) 13 being non-zero then re-establishes the initial mass difference between the smuon and neutralino allowing the relic density to stay within the Planck limits. If one relaxes the assumption that the neutralinoχ 0 1 is the only dark matter candidate, i.e. relax the lower limit on the relic density, then the caracteristic shape observed for (δ T ) 13  Any NMFV parameters among those listed in table 4 whose distributions are not detailed here do not have any interesting phenomenona associated with the imposed constraints, therefore the reader can deduce the full effect and resulting ranges from table 5. Recall that for this scenario, the parameters (δ T ) 23 , (δ F T ) 21 , and (δ F T ) 32 have been set to zero due to requirements for a physical spectrum and neutral LSP. For all δ T T parameters, the main requirements are for a physically relevant spectrum and uncharged LSP, hence we conclude that the prior distribution dominantes over flavour observables that we test against here. Finally, we do not discuss the δ F T parameters as the corresponding results are much the same as for the scan around Scenario 2 presented in the following.
From the discussed results related to reference Scenario 1, it is clear that varying the NMFV parameters individually is not sufficient to properly explore the entirety of parameter space. For this reason, we do not discuss individual variations any further.

Scan around Scenario 2
Here, we discuss selected results of the simultaneous scan of all 15 NMFV parameters around Scenario 2. NMFV parameters are varied according to the ranges given in table 4, while the MFV parameters are fixed to the values given in table 1. Note that the change of the MFV parameters as compared to Scenario 1 allows the variation of all 15 NMFV parameters, while three of them were set to zero for Scenario 1. This yields limits on the full range of flavour violation allowed in Scenario 2.
Starting the discussion with the parameter (δ T ) 13 for which we present the resulting prior and posterior distributions in figure 8, we observe the same feature as for Scenario 1 (see figure 7), but more pronounced. Again, slightly positive or negative values for (δ T ) 13 counteract the effects of other NMFV parameters on the neutralino relic density as explained previously. Coming to the parameter (δ F ) 13 , figure 9 shows that, rather than a single observable having a clear effect, cumulatively µ → eγ, µ → 3e, and Ωχ0 1 h 2 constrain the parameter together with each having a similar effect. Here, we see particularly the effect of flavour violating muon decays on (δ) 13 parameters as elaborated upon in the beginning of section 4.
In the same way as for Scenario 1, all δ T T parameters are constrained by the "prior" requirement of a physical mass spectrum and a neutralino dark matter candidate. Flavour observables have a negligible effect (see table 5).
An example of the posterior distribution for δ F T parameters is shown in figure 10, namely for (δ F T ) 13 . This parameter is constrained almost entirely by the relic density bound, as can be seen in the similarity of the two panels. Let us recall that complete information on limits and dominant constraints of all NMFV parameters associated with Scenario 2 is summarized in table 5.

SUSY scale NMFV parameters for Scenario 2
While from the model-building point of view it is useful to explore the allowed level of flavour violation at the GUT scale, it is equally important to explore the resulting physics at the SUSY scale. Renormalization group running from the GUT scale to the SUSY scale will break the unification conditions given in eq. (2.9) and consequently in eq. (2.11). The fact that these relations are not valid any more below the GUT scale is an essential and intrinsic part of Grand Unification. The present section is devoted to highlighting selected results related to the NMFV parameters obtained at the SUSY scale. More precisely, we study the behaviour of different SUSY scale NMFV parameters which stem from a single NMFV parameter at the GUT scale.
In figure 11 we show the example of (δ F ) 12 , defined at the GUT scale, and the two resulting SUSY scale parameters (δ L LL ) 12 and (δ D RR ) 12 , which belong to the slepton and down-type squark sectors, respectively. First, we see that the prior distribution is altered

JHEP03(2019)067
by the renormalization group effects between the GUT scale in panel a) and the SUSY scale distributions in panels b) and c). The imposed flat priors at the GUT scale are transformed into almost Gaussian-like distributions at the SUSY scale. Looking at the corresponding posteriors, the SUSY scale distributions look even more peaked than the corresponding GUT scale histrograms. Second, it is interesting to note that, at the SUSY scale, the allowed range for the hadronic parameter (δ D RR ) 12 is wider than that for the related leptonic parameter (δ L LL ) 12 in the simultaneous scan. This behaviour is somewhat unexpected, since the gluino running, which is blind to flavour, drives the diagonal squark mass parameters higher, while it leaves the leptonic ones unaffected. In turn, this is expected to reduce the squark NMFV parameters once normalized as per eq. (2.4) [39]. We find that this behaviour is confirmed for all NMFV parameters stemming from individual scans (see examples in figure 11 panels d) and e)), agreeing with the results presented in ref. [39]. However, for the δ F parameters, the reverse is true when considering the simultaneous scan. We suspect that strong renormalization group effects are the cause of this feature, due to the fact that multiple NMFV parameters interact with each other during the evolution from the GUT scale to the SUSY scale.

Parameter correlations
In this section, we examine more closely the correlation between certain NMFV parameters, mentioned already several times in the above discussion, and being the reason that scanning over all parameters simultaneously is ultimately required. The key is that cancellations may exist between the contributions from certain parameters in the calculation of a given observable. However, dealing with analytical results for the different experimental constraints is difficult and beyond the scope of this work. Instead, we choose to take advantage of the numerical results, showing posterior distributions of more than one NMFV parameter together.
The first panel in figure 12 shows viable parameter points that seem to follow a "golden line", with an increased density of points concentrated around a linear relationship between the GUT scale parameters (δ F ) 12 and (δ F T ) 12 . Indeed, the impact of BR(µ → eγ) is suppressed in this line due to cancellation between the two parameters in the analytic expression for this observable. One can also see this in the right panel, only those points lying close to or along said correlation line are consistent with the experimental limits. Said correlation could provide an interesting hint for future SUSY GUT model building.
The analytic expression for the decay rate of µ → eγ can be written as [39] BR( i → j γ) where the branching ratio of the decay i → j ν i ν j is a constant with respect to the NMFV parameters under consideration in the present work. For real NMFV parameters, the form factors F L,R are related to the flavour violating parameters at the SUSY scale according to   The coefficients c i (i = 1, . . . , 4) are combinations of loop factors, masses, and other numerical inputs which can be assumed to be constant in our analysis. Minimizing the form factors F L,R in eq. (4.2) to yield small µ → eγ branching ratios and hence satisfy the experimental constraint leads to relations of the form corresponding to the observed lines in figures 12 and 13. As such, the "golden line" that we recover purely from our numerical analysis is consistent with the analytic formulae for this lepton flavour-violating decay.

Conclusion
In this paper we have considered CP-conserving non-minimal flavour violation in A 4 × SU(5) inspired Supersymmetric Grand Unified Theories (GUTs), focussing on the regions of parameter space where Dark Matter is successfully accommodated due to a light righthanded smuon a few GeV heavier than the lightest neutralino dark matter candidate. Such regions of parameter space are obtained by choosing the second generation T 2 to have a light soft mass, while the heavy gluino mass ensured that all squarks in this multiplet are heavy after RG running to low energy. We have considered two scenarios along those lines, one with a very light right-handed smuon, which is capable of being discovered or excluded by the LHC very soon, but which can account for the (g − 2) µ results, and another scenario with a somewhat heavier smuon. In such regions of parameter space we have found that some of the flavour violating parameters, in particular (δ T ) 13 and (δ F T ) 32 , are constrained by the requirement of dark matter relic density, due to the delicate interplay between the smuon and neutralino masses. By scanning over many of the GUT scale flavour violating parameters, constrained by low energy quark and lepton flavour violating observables, we have discovered a striking difference between the results in which individual parameters are varied to those where multiple parameters are varied simultaneously, where the latter relaxes the constraints on flavour violating parameters due to cancellations and/or correlations. Since charged lepton flavour violation provides the strongest constraints within a GUT framework, due to relations between quark and lepton flavour violation, we have examined in detail a prominent correlation between the flavour violating parameters (δ F ) 12 and (δ F T ) 12 at the GUT scale consistent with the stringent lepton flavour violating process µ → eγ. By switching on both flavour violating parameters together, we have seen that much larger flavour violation is allowed than if only one of them were permitted separately. We have examined this correlation also in terms of the resulting low energy flavour violating parameters in the quark and lepton sectors, and have provided some analytic estimates to understand the origin of the observed correlation.
Precision flavour physics measurements could present challenges to this work and warrant further attention. Particularly, situations such as this often predict small-but-non-zero branching ratios for the LFV decays µ → eγ and µ → 3e, hence stricter bounds on such processes will further limit the amount of NMFV allowed in such scenarios. Figures 12  and 13 are purely data-driven and shows the regions that experimental data prefers; a model which predicts such a correlation could allow reasonable flavour violation and still be preferred over other such models.
In general, we have examined the relation between GUT scale and low scale flavour violating parameters, for both quarks and leptons, and shown how the usual expectations may be violated due to the correlations when multiple parameters are varied simultaneously. We have presented results in the framework of non-minimal flavour violation in A 4 × SU(5) inspired Supersymmetric Grand Unified Theories, with smuon assisted dark matter. Such a framework is interesting since it allows both successful dark matter and contributions to (g − 2) µ , as well as providing the smoking gun prediction of a light right-handed smuon accessible at LHC energies.

A SPheno and the SCKM basis
The CKM basis is the one in which the up-and down-type quark Yukawa matrices are diagonal. The Super-CKM basis (SCKM) is obtained analogously, i.e. the squarks undergo the same rotations as their SM partners. This basis is convenient for phenomenological studies, and allows for a consistent expression of flavour vioation throughout the literature. The different rotations for the SM quark and lepton fields are: where the primed fields are in the flavour basis and the bare fields are in the basis of diagonal Yukawa couplings. The misalignment between up-and down-type quarks leads to the usual CKM matrix: In order to account for the change to the SCKM basis, the numerical programme SPheno assumes diagonal down-type Yukawa matrices. In this case the CKM matrix is the following: In SU(5)-like models, the choice of the representations F =5 and T = 10 forces relationships between Yukawa couplings to hold at the unification scale: As a consequence, we have V u L = V u R in this case, meaning both lepton and down-type Yukawas are simultaneously diagonal. For consistency, we then have to perform a systematic CKM rotation for all terms involving V u L and V u R . The soft-breaking terms of the

JHEP03(2019)067
Lagrangian transform as follows when switching to the SCKM basis: Consequently, in the SCKM basis, where the down-type Yukawa matrix is diagonal following the SPheno requirements and assuming the SU(5) relations, the 6 ×6 soft mass matrices in the MSSM (once trilinear couplings have been taken into account) are up to the D-terms and SM masses. Since SPheno automatically ensures the CKM rotation for the left-left block of M 2 U , we enforce the rotation for the other blocks of the up-type squark mass matrix by hand before running SPheno.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.