Muon g-2 through a flavor structure on soft SUSY terms

In this work we analyze the possibility to explain the muon anomalous magnetic moment discrepancy within theory and experiment through lepton flavor violation processes. We propose a flavor extended MSSM by considering a hierarchical family structure for the trilinear scalar Soft-Supersymmetric terms of the Lagranagian, present at the SUSY breaking scale. We obtain analytical results for the rotation mass matrix, with the consequence of having non-universal slepton masses and the possibility of leptonic flavour mixing. The one-loop supersymmetric contributions to the leptonic flavour violating process $\tau \to \mu\gamma$ are calculated in the physical basis, with slepton flavour mixed states, instead of using the well known Mass Insertion Method. We present the regions in parameter space where the muon g-2 problem is either entirely solved or partially reduced through the contribution of these flavor violating processes.


Introduction
It is well known that in contrast to electric charge conservation, lepton number conservation is not associated with a gauge symmetry.In the Standard Model (SM), the spontaneous breaking of the electroweak symmetry produces eigenstates of the remaining gauge group that are not in general eigenstates of the mass matrix [1][2][3][4].But after diagonalization of the mass matrix, the electroweak coupling matrix is also diagonal in the mass basis, therefore there is no possibility for lepton flavor violation.Certainly this is now in contradiction with the experimental evidence on neutrino mixing [5][6][7][8] and also the possible LFV Higgs decay [9] which forces the structure of the models beyond the SM.
The original structure of the SM with massless, and thus degenerate neutrinos, implied separately τ, µ, e number conservation.In particular, the processes τ ± → l ± γ, (l = µ ± , e ± ) through gauge bosons loops are predicted to give 1 very low rates [10], even considering the experimental evidence on neutrino oscillations [5][6][7][8].Under this evidence the amplitudes for the Lepton Flavor Violation (LFV) processes at low energy are suppressed by an inverse power of the large Majorana mass scale M I used in the well-known seesaw model [11,12], which explains naturally the small masses for the active left-handed neutrinos.On the other hand, the experimental bounds for the branching ratio BR(τ ± → µ ± γ) [13] set strong restrictions on models of physics beyond the SM.
A realistic possibility of physics beyond the SM is offered by supersymmetry (SUSY), whose simplest realization containing the SM is the Minimal Supersymmetric Standard Model (MSSM) (see for instance [14]).In terms of supersymmetry, the SM is embedded in a larger symmetry which relates the fermionic with the bosonic degrees of freedom.As a consequence of this higher symmetry, the SUSY framework stabilizes the electroweak scale, provides us with dark matter candidates, as well as with a greater possibility of unification of all fundamental interactions and a solution to the hierarchy problem.
The discovery of the Higgs boson [15][16][17][18] and the search for sparticles at the LHC, have modified the parameter space of supersymmetry as a near electroweak (EW) scale model [19][20][21][22].The MSSM, as the first minimal supersymmetric extension of the SM, was conceived to be near to the electroweak scale, in order to set viable phenomenological scenarios to analyze with available experimental data.One important issue to be considered was the experimental absence of Flavor Changing Neutral Currents (FCNC), which lead to the simplifying assumption of universality in the sfermion scalar masses, keeping the desired good behavior of FCNC's (i.e.bounded) and in addition, reducing the number of free parameters.
The Constrained Minimal Supersymmetric Standard Model (CMSSM) was conceived under the assumption of Grand Unified Theories (GUT) structures.It considers in particular universal sfermion masses and alignment of the trilinear soft scalar terms, A f,ij to Yukawa couplings at the unification scale [23,24].Nevertheless, neutrino oscillations made it imperative to reconsider the flavour structure in the theoretical models.
The most recent LHC data points to a heavy spectrum for some of the SUSY particles in the case this constrained model were realized in nature.The relation between the Higgs mass and the fermions and sfermions masses in supersymmetric models indicate either higher stops masses or large mixture within stops [25].It is the squark sector, and particularly the stop and gluino, which tend to lift the mass scale of the MSSM [20-22, 26, 27].However, for the slepton sector the LHC data for the exclusion bounds are less restrictive and masses may still be below the TeV scale [28]. 2 .On the other hand, we could go beyond the constrained MSSM and explore other possibilities for the flavor structure.It is thus very relevant to search for SUSY effects to indirect electroweak precision processes through quantum corrections involving superparticles, as the phenomenologically viable parameter space is modified by experimental data, being this the main motivation of the present work.
In the MSSM the conventional mechanism to introduce LFV is through the Yukawa couplings of the right handed neutrinos, N i , which generate offdiagonal entries in the mass matrices for sleptons through renormalization effects [29,30], particularly in the LL block.Then the predicted rates for the τ → µγ and µ → eγ decays are not suppressed, and depend on the unknown Yukawa matrix elements, but they will not be detected in the future experiments if those elements are too small.In Ref. [31] the authors work also with these LFV processes, using the seesaw mechanism in the SM [32] and supersymmetric models to extended neutrino and sneutrino sectors, and perform the one-loop calculation through the Renormalization Group Equations (RGEs) based on leading-log approximation.In the SM they use the neutrino-gauge loops, while in the supersymmetric model they get the sneutrino-chargino loops.In Ref. [33] the authors noticed that there is another source of LFV, namely the left-right mixing terms in the slepton mass matrix, and that their contributions to the LFV processes can be large even when the off diagonal Yukawa couplings elements are small.Later, in a second paper [34], they incorporated the full mixing of the slepton masses and mixing in the neutralino and chargino sector and then performed a numerical diagonalization of the slepton mass matrices.An interesting result of their analysis is that the contribution from the left-right mixing is only important in the region where the mixing term is m τ µ tan β ∼ O( m2 S ) and they consider the trilinear soft terms A E,ij contribution negligible.In the above expression m τ is the tau mass, µ (µ susy throughout this paper3 ) is the Higgsino mass parameter;

is the ratio of Higgs vacuum expectation values (vevs) and m2
S is the supersymmetric scalar mass scale from the soft SUSY breaking.It is worth noting, though, that this analysis was done with very different considerations on experimental data than those we have now.
A more recent work on this relation between the seesaw mechanism for neutrino mixing and charged lepton flavor violation is done in Ref. [35], where a non-trivial neutrino Yukawa matrix at the GUT scale leads to non-vanishing LFV interactions at the seesaw scale through the RGEs.Another approach to the same problem has been done using high-scale supersymmetry in Ref. [36], where the Majorana mass matrix of right-handed neutrinos is taken to be diagonal and universal, while the neutrino Yukawa matrix is proportional to the neutrino Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix U P M N S , and the product of the left and right handed neutrino masses is . This neutrino Yukawa matrix, which would be present in low energy phenomenology, changes also with the RGE running of the soft SUSY breaking parameters.This scheme of FV was proposed in Ref. [37], where small offdiagonal elements of the slepton mass matrix are considered and, in the interaction basis, the FV processes are restricted by using these off-diagonal elements as free parameters; here the trilinear coupling is considered to be zero, A = 0.In Ref. [38] the trilinear coupling A 0 is considered only for the LR flavor mixing term, in the LR term of the corresponding slepton A 0 is set it to zero.There is also a more general phenomenological work considering non-diagonal LL, RR and EL blocks of sfermion mass matrices which are parameterized as a sfermion mass product and a free parameter for each matrix element in order to do a numerical evaluation of the processes in the mass basis [39], having all the elements of the 6 × 6 sfermion mass matrix as parameters that might be constrained by the LFV processes.Recent analysis of these general FV contributions are done in [40,41].This general sfermion mass matrix, although complete, implies a considerable increase in the number of parameters.Nevertheless, the authors found in seven different possible scenarios an upper bound for their off-diagonal parameter.We must say here that in most of the literature, although the calculation is done in a physical basis, what is done is a diagonalization of 2 × 2 blocks flavor sleptons and they still consider a flavour mixing parameter, which is off-diagonal on the mass matrix and is used as coupling in the MIA method, so their physical basis means that instead of using the interaction basis states li,L , li,R , they use li,1 , li,2 with i = 1, 2, 3 as flavors.
There is as well work on supersymmetric models where R − parity violation is considered in the allowed superpotential operators [42], with the consequence of having LFV couplings directly present in the model.
A very important issue to be considered when lepton flavor mixing is allowed is the extra contribution to the anomalous magnetic moment of the muon.The experimental value of the g − 2 is another element of the electroweak (EW) precision data which has not been completely explained by the SM [43][44][45], despite the efforts that have been made for improving the hadronic contribution calculations [46][47][48], the dominant source of uncertainty in the theoretical prediction.It is well known that the main MSSM contribution to g − 2 (we will call it a µ ), involves neutralino-slepton and chargino-sneutrino loops [37,49].Even the two-loop contribution in terms of tan β has been calculated in Ref. [50], where a reduction was found of the discrepancy coming from an extra contribution, within 14% to 6% of the one-loop MSSM contribution, depending on different scenarios of parameter space.
In Ref. [51] the supersymmetric calculation of a µ has been updated considering both the chargino-sneutrino loop and the neutralino-smuon loop.It was found that the chargino-sneutrino loop dominates, especially in the case where all the scalar masses are degenerate and, on the other hand, when the µ susy parameter is large, then χ0 − μ could be enhanced.There has also been work done relating the parameters for g-2 anomaly, flavour violation, and h → γγ in [52].
In this work we present an analysis of a flavor violating extension of the MSSM (FV-MSSM) one-loop contribution to a µ , which is driven by a LFV mechanism at tree level.The LFV process τ → µγ is used as an additional constraint of the parameter space of the FV-MSSM.Our strategy for the implementation of LFV consists in assuming that A E,ij -terms follow a particular structure in the context of textures.Furthermore, we take an ansatz for the mass matrix for sleptons, allowing an exact diagonalization [53] that results in a non-universal spectra for sfermion masses, providing a clear way for having flavour mixing within sleptons at tree level and the opportunity to work in the mass eigenstates basis.Concerning the extra contribution to the anomaly coming from the FV-MSSM, we assume that it comes mainly from the slepton-bino loop, a l B µ , and we compare with the usual MSSM contribution from this loop.
The paper is organized as follows: In Sect. 2 we present the flavor structure of sleptons from an ansatz for the trilinear scalar terms.Then in Sect.
3 we show the one-loop analytical calculation of BR(τ → µγ).In Sect. 4 we include the a µ calculation and present the combined results in Sect. 5. Finally, we discuss our conclusions in Sect.6.

Flavor structure in the soft SUSY breaking Lagrangian
If supersymmetry exists in Nature it has to be broken, since there is no evidence that these new particles exist at low energies [19].This symmetry breaking is achieved by the introduction of terms in the Lagrangian, which break SUSY in such a way as to decouple the SUSY partners from the SM particles, and at the same time stabilize the Higgs boson mass to solve the hierarchy problem (see for instance [23]).The soft SUSY breaking Lagrangian in general includes trilinear scalar couplings A H ij , as well as bilinear couplings B ij , scalar squared mass terms M 2 f , and mass terms for the gauginos M 2 i .
Specifically, for the scalar fermion part of the soft SUSY terms in absence of flavor mixing, as is considered in the MSSM, it will have the following structure: where f are the scalar fields in the supermultiplet.In the case of sfermions the L, R are just labels which point out to the fermionic SM partners, but as we are dealing with scalar fields they have no longer left and right SU (2) properties.In general they may mix in two physical states by means of a 2 × 2 rotation matrix, fL fR ↔ f1 f2 .
The first terms in (1) contribute to the diagonal terms of the 2 × 2 sfermion mass matrix, while the second ones are Higgs couplings with the different sfermions, and they contribute to the off-diagonal L − R terms of the mass matrix once the EW symmetry is spontaneously broken.As i is a flavour index we can see that Eq. ( 1) implies no flavor mixing.
In our case, where we do consider flavour mixing in the trilinear terms, A ij f would be a general 3 × 3 matrix, since we consider together the three flavours, with two scalar fields for each one.The complete fermionic trilinear terms are given as Here Q i is the squark doublet partner of the SM SU (2) left doublet and U j , D j are the corresponding squarks singlets, while Li is the slepton doublet and E j is the singlet.In this work in particular, we only analyze the sleptonic part.We will explain further in this paper the ansatz flavour structure we consider for this.Once the EW symmetry breaking is considered, the above Lagrangian (2) for the sleptonic sector takes the form The soft terms are not the only contributions to the sfermion mass elements, the supersymmetric auxiliary fields F and D coming from the superpotential also contribute to this mass matrix as we explain in the next section.

Mass matrix for sfermions
The contribution to the elements of the sfermion mass matrix come from the interaction of the Higgs scalars with the sfermions, which appear in different terms of the superpotential and soft-SUSY breaking terms as is fully explained in [54,55].In the case of the slepton mass matrix, as we said before, the contributions coming from mass soft terms are M 2 l,LL , M 2 l,RR , from trilinear couplings after EW symmetry breaking A l ij and from the F, D−terms.We arrange them in a block mass matrix as follows: The F f and D f are the auxiliary fields in the supermultiplets, which are introduced to have the same bosonic and fermionic degrees of freedom, but are dynamical spurious [14].The F -auxiliary field comes from the Higgs chiral superfields and contributes to the mass matrix as follows: From the D−auxiliary fields which come from the scalar superfields of fermions we have the following mass terms: where l = e, µ, τ .The elements of the sleptons mass matrix Eq. ( 3), for the different flavors given by i, j = e, µ, τ are

Soft trilinear terms ansatz
The lepton-flavor conservation is easily violated by taking non-vanishing offdiagonal elements for each matrix, the size of such elements is strongly constrained from the experiments.In the CMSSM, it is assumed that the soft sfermion 2 × 2 mass matrices m2 E , m2 L are proportional to the identity matrix, and A e,ij is proportional to the Yukawa matrix y e,ij .With these soft terms the lepton-flavor number is conserved exactly [33].The non-universality of scalar masses has been studied in supersymmetric models in the context of string theory [56].In Ref. [57], the authors assume a non-universality of scalar masses, through off-diagonal trilinear couplings at higher energies.In Refs.[58,59] a SU(3) flavor symmetry is introduced, then by means of the Froggat-Nielsen mechanism the associated flavon fields acquire vevs, which upon spontaneous symmetry breaking generate the couplings which mix flavours.
In the present work, we assume m2 RR,l ≈ m2 LL,l = m2 S but we propose that there is a mixing of two of the scalar lepton families in the LR mass terms.This mixing may come from a discrete flavor symmetry, as could be the extension of the SM with S 3 [60][61][62], or supersymmetric models with Q 6 [63][64][65][66], which have the fermions assigned to doublet and a third family in a singlet irreducible representations.In order to analyze the consequences of this flavor structure we construct an ansatz for the trilinear terms A t .Our procedure is similar to the work done in Ref. [67] for FCNC's in the quark sector through an ansatz of soft-SUSY terms.In our case we consider the whole two families contributions and of the same order of magnitude, having the following form for the trilinear term: In this case one could have at tree level the selectrons in a singlet irrep., decoupled from the other two families of sleptons.This would give rise to a 4 × 4 matrix, diagonalizable through a unitary matrix Z l, such that Since we assumed that the mixing is in the smuons and staus only and the selectrons are decoupled, the remaining 4 × 4 smuon-stau mass matrix will have the following form: where This way we will have physical non-degenerate slepton masses.4 where We may write the transformation which diagonalizes the mass matrix as in Ref. [53], as a 4 × 4 rotation matrix for sleptons Z l, which is in turn a 2 × 2 block matrix Z † l M 2 µ−τ Z l = M 2 l,diag , explicitly having the form where σ 3 is the Pauli matrix and The non-physical states are transformed to the physical eigenstates by where In the case of the MSSM without slepton mixing we would need to revert the similarity transformation performed as Z l M 2 l,diag (y = 0)Z † l = M μ,τ , vanishing also the mixing parameter, y = 0. Then we will get a diagonal by blocks matrix, where the two 2 × 2 bloques are the mass matrix for smuons and staus, respectively which can in turn be diagonalize separately as in the usual MSSM, obtaining the two sleptons physical states l1 , l2 for each flavor that we identify with the MSSM slepton eigenstates.The masses for the smuons would then be the usual ones, where

Neutralino-lepton-slepton interaction
We assume the usual MSSM form of neutralinos as a mixing of the fermionic part of vector superfields, i.e. gauginos and Higgsinos.The symmetric mass matrix for neutralinos is given by The diagonalization of the mass matrix implies transformation of the neutralinos as In the rotation matrix η is a diagonal matrix, whose elements η j are introduced in such a way as to change the phase of those neutralinos whose eigenvalues become negative after diagonalization, i.e. η j = 1 for m χ 0 j > 0 and η j = i for m χ 0 j < 0. The general interaction Lagrangian for neutralino-fermion-sfermion in the MSSM is given as follows [54] where the (l n ) and (r n ) are the left and right fermion-neutralino couplings, respectively.In this expression the P L,R are the ordinary chiral operators, and the labels for the corresponding scalar superpartners of fermions are L for sfermions X = 1, 2, 3 and R for X = 4, 5, 6 in the interaction basis and g is the U (1) coupling constant.
The neutralino-fermion-sfermion couplings in equation ( 18) are given by where ηΘ N is the rotation matrix which diagonalizes the neutralino mass matrix [68].Now, considering the sleptons mass eigenstates given in (15) we rewrite the neutralino-lepton-slepton interaction Lagrangian as where So, we can see here that we directly introduce the FV into the interaction Lagrangian avoiding the need of a mass insertion in the propagators of the loops.

BR(τ → µ + γ)
In general, the way lepton flavor violation is introduced in calculations in the supersymmetric loops is using the approximation method called Mass Insertion Approximation (MIA) [39,[69][70][71], which uses a Taylor expansion in a mass parameter [72] giving qualitative good results [73].Then the calculation is done in a non-mass eigenstate basis expanding around the universal squark masses [74].This method assumes that off-diagonal elements are small, which generates a strong restriction on the allowed SUSY parameters.On the other hand, working in the interaction basis the number of loops to be calculated is reduced to one, giving a simple analytical expression for the free parameters involved.Concerning flavour violation via neutrino and sneutrino mixing, including a right-handed neutrino [33], the MIA method is used to compute the one-loop amplitude for this process.
In this paper, rather than using the MIA method, we work in a physical basis by diagonalizing exactly the complete mass matrix obtaining mixed flavour sleptons, introducing only two free parameters, which we reduce to one by considering w = 1, assuming the soft trilinear term ansatz proposed in the previous section, Eq. (9).
We now use the couplings obtained to calculate FV processes to establish the feasibility of the ansatz.In particular, we calculate the supersymmetric sfermion-neutralino one-loop contribution to the leptonic flavor violation process τ → µ + γ, which corresponds to the Feynman diagram given in Fig. 1.The experimental bound to the branching ratio for this decay at 90% C.L. [13] is The loop diagrams shown in Fig. 1 are IR safe.A photon is radiated either by a slepton inside the loop or by the external lepton, all three diagrams are needed to achieve gauge invariance.To simplify the expressions, we have assumed that the lightest neutralino is mainly a Bino ( B), although the procedure can be generalized to any type of neutralino.
Figure 1: One Loop diagrams in the LFV process τ → µγ.The total amplitude is gauge invariant and finite in the UV region.
Considering the limit M 1 , M 2 , µ susy m Z [68], then the lightest neutralino is mostly Bino χ0 1 ≈ B then we take (Θ N ) 1i ≈ δ 1i in Eq. ( 22).The mass eigenvalue for the lightest neutralino is given by [23] m Then this would be a Bino-like neutralino in the limit for numerical values M 1 < M 2 |µ susy |.In this case the Bino-lepton-slepton coupling can be written as follows: where l runs over the eigenstates l1,2,3,4 given by Eq.( 15).For the τ → µ + γ decay the scalar and pseudoscalar couplings are given in Table 1.The total amplitude is gauge invariant and free from UV divergences, as it should be, and it can be written in the conventional form, where the one-loop functions E and F contain the sum of the contributions from sleptons l1,2,3,4 running inside the loop, .
The functions E l, F l are written in terms of Passarino-Veltman functions and can be evaluated either by LoopTools [75] or by Mathematica using the analytical expressions for C0 and B0 [76], where we have defined the ratio x = mµ mτ , and possible values of r = 1, 2 set by lr , and the η( l) function as follows: The differential decay width in the τ rest frame reads where p µ is the 3-vector of the muon.The branching ratio of the τ → µ + γ decay is given by the familiar expression, 4 The MSSM and the muon anomalous magnetic moment a µ The anomalous magnetic moment of the muon a µ ≡ g−2 2 is an important issue concerning electroweak precision tests of the SM.The gyromagnetic ratio g, whose value g = 2 is predicted at lowest order by the Dirac equation, will deviate from this value when quantum loop effects are considered.A significant difference between the next to leading order contributions computed within the SM and the experimental measurement would indicate the effects of new physics.
The experimental value for a µ from the Brookhaven experiment [77] differs from the SM prediction by about three standard deviations.In particular, in Ref. [43] it is found that the discrepancy is where a th µ is the theoretical anomalous magnetic moment of the muon coming only from the SM.
Three generic possible sources of this discrepancy have been pointed out [78].The first one is the measurement itself, although there is already an effort for measuring a µ to 0.14 ppp precision [79], and an improvement over this measurement is planned at the J-Parc muon g-2/EDM experiment [80] whose aim is to reach a precision of 0.1 ppm.
The second possible source of discrepancy are the uncertainties in the evaluation of the non-perturbative hadronic corrections that enter in the SM prediction for a µ .The hadronic contribution to a µ is separated in High Order (HO) and Leading Order (LO) contributions.The hadronic LO is under control, this piece is the dominant hadronic vacuum polarization contribution and can be calculated with a combination of the experimental cross section data involving e + e − annihilation to hadrons and perturbative QCD [48].The hadronic HO is made of a contribution at O(α 3 ) of diagrams containing vacuum polarization insertions [81,82] and the very well known hadronic Light by Light (LbL) contribution, which can only be determined from theory, with many models attempting its evaluation [83,84].The main source of the theoretical error for a µ comes from LO and LbL contributions.It is worth mentioning that the error in LO can be reduced by improving the measurements, whereas the error in LbL depends on the theoretical model.
The third possibility comes from loop corrections from new particles beyond the SM.There have already many analyses been done in this direction (see for instance [33,85,86]).
To calculate one-loop effects to g-2, for general contributions coming from different kind of particles Beyond the SM, there is a numerical code built using Mathematica [87].
The supersymmetry contribution to g-2, a SU SY µ , was first computed by Moroi Ref. [37] and recently updated in Ref. [88].In these works the large tan β scenario was studied, showing the dominance of the chargino-sneutrino loop over the neutralino-smuon loop, provided the scalar masses are degenerate, otherwise the µ susy parameter (Higgsino mass parameter) must be large allowing an enhancement of the muon-neutralino loop (χ 0 − μ).It was also shown that in the interaction basis the dominant contributions are proportional to µ susy M 1 tan β, then the sign and the size of the contribution to a SU SY µ depends on the nature of this product.Hence, the supersymmetric contributions to the anomaly are determined by how these elements are assumed (see for instance [37,88]).The results in the literature are usually obtained using the MIA approximation, however, there are some schemes where the work is done in the physical basis (e.g.[41]).The difference with the MIA method is not only the change in basis, but the restriction that is imposed a priori that some elements in the mass matrix are considered small compared to the diagonal ones.
There has been research toward an MSSM explanation to the g − 2 discrepancy related to LFV as in [89,90], since there is a correspondence between the diagrams in the MSSM that contribute to the anomalous magnetic moment of the muon and the diagrams that contribute to LFV processes.The process µ → e + γ have been used to constrain lepton flavor violation and as In this work we assume that there is room for an MSSM contribution to g − 2 through lepton flavor violation in the sleptonic sector.In particular, we search for the LFV process τ → µ + γ and calculate g − 2 through a mixing of smuon and stau families, a l B µ , Fig. 2. The ansatz proposed here avoids extra µ → e + γ contributions.To establish the restrictions on parameter space we consider a loose constraint, µ indicates that the lepton flavor violation supersymmetric loop through charged sleptons is not necessarily the only contribution to solve the discrepancy, Eq. (31).We also show the extreme case in parameter space where this loop contribution solves the discrepancy completely |a Exp µ −a T H µ | ≤ 1σ.
When taking into account the slepton-bino flavor violation contribution to g − 2, if the discrepancy is ≤ 1σ, it means that this contribution solves the whole g−2 problem.In the opposite scenario, |a Exp µ −a T H µ | ≈ 3.3σ means that the slepton-bino loop gives no significant contribution to the discrepancy.In here we will look at a possible contribution to g − 2 between both scenarios.
Using the LFV terms constructed previously we obtained the contribution to the anomalous magnetic moment of the muon a l B µ .Defining the ratio r = m µ /M 1 and taking the leading terms when r 2 → 0, and M 1 = m B as the Bino mass.
In order to compute the SUSY contribution to the g − 2 anomaly, we follow the method given in Ref. [91].All we have to do is to isolate the coefficient of the (p 1 + p 2 ) µ term, in other words, computing the one-loop contribution, we can write the result as follows: where the ellipsis indicates terms that are not proportional to (p 1 + p 2 ) µ .Then the anomaly can be defined as δa µ = g−2 2 = F 2 (0) with q = p 2 − p 1 .Keeping in mind that we require the magnetic interaction which is given by the terms in the loop process proportional to (p 1 + p 2 ) µ we write it as Considering only these terms in the interaction and gathering them, the contribution of the flavour violation loop to the g − 2 anomaly due to a given slepton l reads where , and ∆ l B = m 2 l − M 2 1 , having four contributions with l running from 1 to 4 with the values of the couplings S Bµ, l, P Bµ, l are given in Table 1.This expression is equivalent to the one presented in [92] and can be written using their notation as can be found in Appendix B.
The expression will be different from MIA because the off-diagonal elements LR are not explicit since we are in the physical basis.In the interaction basis, the LR terms appear with explicit SUSY free parameter dependence as they use directly the elements of the slepton mass matrix.Exact analytical expressions for the leading one-and two-loop contributions to g-2 in terms of interactions eigenstates can be found in Refs.[49,92], and references therein.By taking these expressions in the limit of large tan β and of the mass parameters in the smuon, chargino and neutralino mass matrices equal to a common scale M SU SY , the results calculated in the mass-insertion approximation in the same limit [37] are reproduced from the complete forms given in [92].We have explicitly checked that our one-loop results when no LFV terms are present coincide with the analytical expressions of ref. [92], and thus in the appropriate limits also with the MIA expressions.Our expressions for the contribution of the LFV terms to g-2 can be found in Appendix B.
Here we take a flavour structure with no a priori restrictions on the size of the mass matrix elements other than two family mixing, and the restrictions come directly from the comparison with experimental data.
The table shows the parameter space where the scan was performed.The values were taken at random for each variable within the bounds shown .

Results
We now analyze the region in parameter space allowed by the experimental bound on BR(τ → µγ), taking into account that the mixing parameters w, y represent at most a phase, i.e. the mixing terms in the LR term of the mass matrix are of the same order as A 0 , see Eq. ( 9), in contrast with the MIA method where this terms are considered small compared with the diagonal ones which is needed to apply the method.In the parameter space region comprised by Table 2, we are able to safely consider lepton flavour mixing in trilinear soft terms of the MSSM, and constrained it at the current experimental bounds BR exp (τ → µγ) < 4.4 × 10 −8 [93].Throughout parameter space we take M 1 < |µ susy |.We highlight the points where the g − 2 is solved completely, shown in black in all figures.In order to ensure that the lightest neutralino is mostly Bino, we further assume for these points M 1 1 3 |µ susy |.We found for the parameter values given in Table 2 that the BR(τ → µγ) is only partially restricted from experimental bound for mS 3200 GeV, also for M 1 4.5 TeV.Table 3 shows examples of different sets of values for random parameters given within the range in Table 2, consistent with the experimental bound on LFV and that also solve entirely the g-2 discrepancy, in all these points the Bino is considered as the LSP.From these sets of values it can be seen that the g − 2 discrepancy can be solved within the FV-MSSM by different possible combinations of the parameters.
The difference between the experimental value and the SM prediction for the anomalous magnetic moment, Eq.( 31), gives σ = 80 × 10 −11 .As we have already explained we distinguish between two possible ways the slepton contribution should be constrained, depending whether the loop is dominant in FV-MSSM or not: It is important to mention that we take the points that solve for "any contribution" as defined above (blue in graphs), because we are aware that this is only one of the possible supersymmetric contributions to g − 2. In a more general case we need to include the chargino-sneutrino contributions in order to have an entire picture of the parameter space.In the FV extension considered here this contribution will be the same as in the usual MSSM.For a more complete treatment right-handed neutrinos should be considered, together with LR mixing and the trilinear term.
In Fig. 3 we show the dependence of the BR(τ → µγ) on mS and on the Bino mass M 1 , and it is shown the stringent restrictions for these masses.In Fig. 4 we show the value of a l B µ for different values of the Bino and the SUSY scalar mass, the color code is clear from this figure.The blue points correspond to the mass scale for which there is any contribution to the discrepancy a µ Eq. (35).The black ones are those for which the discrepancy would be completely explained by the LFV contribution Eq. ( 36), for these points we take M 1 < 1 3 |µ susy | (otherwise we just take M 1 < |µ susy |).The red points are outside these ranges, i.e. are contributions non-compatible with experimental data of the muon g-2 anomaly to be solved.The green points show the results obtained by taking y = 0 in our ansatz, i.e. no FV, and calculating the smuon-Bino loops for g − 2 with the smuons masses as given in Eq. ( 17) and considering a trilinear coupling as A 0 = 0.
Figure 5 shows the relation of µ susy with tan β and trilinear coupling A 0 tan β M 1 (GeV) µ susy (GeV) mS (GeV) A 0 (GeV) for values for which the a µ discrepancy receives contributions from the LFV terms.We see that there is a quite symmetrical behavior for any sign of µ susy .
In order for the a µ problem to be solved entirely by LFV |µ susy | 4000 GeV and no restriction for tan β.For smaller values of A 0 there will be less restriction on µ susy .Although µ susy values could be restricted by other sectors of the MSSM, e.g. the radiative corrections to the lightest Higgs mass [25,94].
On the other hand, there are other SUSY models, where the value of µ susy could be naturally small [95].
Figure 6 shows the ratio of the Bino mass M 1 with SUSY scalar mass mS where the points showed are solutions to a µ discrepancy achieved up to 1σ by the LFV contribution.We see a highly restricted regions for 1.5 < M 1 / mS < 2.5, although we also have few points within 0.4 M 1 / mS < 0.85, but there are no points for 0.8 M 1 / mS < 1.2.We also see the behavior of these points the scalar mass is highly restricted to the range of values m S [500, 3400] GeV, reaching the top values for larger values of |µ susy | We consider that in the region of parameter space where the points that solve completely the g − 2 anomaly lie, the Bino-sleptons loop contribution will dominate over the chargino-sneutrino contribution.Under this consideration is possible that the allowed parameter space is different from the MSSM with no FV terms in the charged lepton sector, where the chargino-sneutrino contribution is the dominant one [96].

Summary and conclusions
We proposed an ansatz for the trilinear scalar couplings considering a two family flavour structure.We obtain a non-universal slepton spectrum and and slepton states are now flavor mixed.This specific family structure implies the possibility of lepton flavour violation although avoids extra LFV contributions to BR(µ → eγ).In the method we used the FV is absorbed into the Lagrangian couplings instead of introducing a mass-insertion term into the propagator as used commonly in the literature.This method does not require a priori approximations to reduce the loop amplitude integral expression.We analyzed the parameter space which gives values for these processes within experimental bounds.We considered that the lightest neutralino is mainly a Bino, specifically we consider the slepton-Bino loops.In order to have the Bino decoupled from Higgsino we take M 1 < |mu susy |.Under these assumptions we showed that this FV couplings will include a mixture of four types of sleptons running in the loop contributing to a µ , which in the interaction basis corresponds to the smuons and the staus, as can be seen in Fig. 2, and that for certain regions of parameter space it is possible to solve entirely the discrepancy between the experimental and theoretical values of a µ , in this case we specifically take a more restricted condition, M 1 < 1/3|mu susy |.The points that match with these conditions are given for the scalar SUSY mass scale mS involved in the LFV processes range between 450 mS < 3300 GeV, the upper bound in the scalar mass is reached for |µ susy | ∼ 14 T eV .The possible Bino mass needed in order to solve the a µ problem ranges from ∼ 350 GeV to ∼ 7.5 TeV, nevertheless the ratio of these masses is restricted to 0.4 M 1 mS < 3, although we have very few point for M 1 mS < 0.9, and the points around M 1 mS ∼ 1 are excluded.It is possible to contribute only partially to the a µ problem, in which case a much larger parameter space is allowed (blue points).This partial contribution to g − 2 will be important when the chargino-sneutrino contribution is included, since it might change the allowed parameter space.This complete analysis we leave to a forthcoming work.Nevertheless, is worth mentioning again that it is natural to have differences in the parameter space as compared to the usual MSSM, where the chargino-sneutrino contribution is the dominant one.
It is interesting to notice that considering off-diagonal elements in the LR of the mass matrix block to be as large as 1TeV does not necessarily blow up the BR(τ → µγ) process, instead, this assumption helps to reduce partially or completely the g − 2 discrepancy.In our case, we have considered offdiagonal terms in the soft trilinear couplings, of the order of 50 GeV < A 0 5TeV.We also compare our results with the no flavour violation g − 2 MSSM one-loop contribution, where we obtain the same expressions given in the literature for complete calculation and in the numerical results we obtain small positive contributions to g − 2 considering no contribution from the trilinear term A = 0.

A Loop amplitude for τ → µγ
We present here the expressions we obtain for the invariant amplitude of the processes given in Fig. 1.For clarity in the expressions we have defined g 2 c = tan 2 θwg 2 1

16
. For general leptons in external particles represented by i, j = e, µ, τ , the diagram in Fig. 1 (a) we have M a = −eg 2 c ū(p 1 ) (S i S j − P i P j ) + (S i P j − S j P i )γ 5 m B 1 (2π) 4 dq 4 2(p 2 + q) • D q D 1 D 2 −eg 2 c ū(p 1 ) (S i S j + P i P j ) + (S i P j + S j P i )γ 5 1 (2π) 4 dq 4 2(p 2 + q) • q D q D 1 D 2 , (37) where D q = q 2 − m 2 B , D 1 = (q + p 1 ) 2 − m 2 lr , D 2 = (q + p 2 ) 2 − m 2 lr , and is the photon polarization vector.For the τ → µγ decay, we have i = τ and j = µ and the S i,j , P i,j couplings are labeled as follows: S i = S Bτ l, S j = S Bµ l, P i = P Bτ l and P j = S Bµ l.All the possible sleptons running inside the loop are indicated by the index l = μ1 , μ2 , τ1 , τ2 .The corresponding values are given in Table 1.For the anomaly g − 2 we set i = j = µ.For the diagram Fig. 1(b) we have with Σ b = m B g 2 c (S i S j − P i P j ) + (S i P j − S j P i )γ 5 1 (2π) 4 dq 4 D q D 1 +g 2 c (S i S j + P i P j ) + (S i P j + S j P i )γ 5 1 (2π) 4 dq 4 q D q D 1 . ( The amplitude for Fig. 1(c) reads where Σ c = m B g 2 c (S i S j − P i P j ) + (S i P j − S j P i )γ 5 1 (2π) 4 dq 4 D q D 2 +g 2 c (S i S j + P i P j ) + (S i P j + S j P i )γ 5 1 (2π) 4 dq 4 q D q D 2 . ( The total amplitude which is the sum of Eqs.(37,38,40) is written as follows: In the case of i = τ and j = µ we would have the expressions for E ij and F ij as in Eqs.(27,28).

B The loop contribution to the muon anomaly
The loop amplitude 5 for the vertex correction is given by ū(p 1 )Γ µ u(p 2 ) = ıg c ū(p 1 ) S Bµ, l + P Bµ, lγ 5 where q 2 = (p 2 − p 1 ) 2 and the ellipsis means terms that are not involved in the determination of the anomaly contribution.The propagators are given by By setting q 2 = 0 and considering that the muon mass is negligible compared to the supersymmetric particle masses inside the loop, the contributions to the anomaly are found to be

Figure 3 :
Figure 3: The plots show the dependence on BR theo (τ → µγ) on the SUSY scalar mass mS (left) and the Bino mass M 1 (right) and on the ratio of them (down).The gray points are excluded by the experimental bound on BR(τ → µγ).The rest of the color code is shown explicit in Fig.4, which separates ranges of FV contributions to g − 2.

Figure 4 :
Figure 4: The plots show the dependence of the value of our calculation for a l B µ with the SUSY scalar mass (left) and the Bino mass (right).Here the color code used in Figs. 3, 5 and 6 is shown explicit as ranges of the a l Bµ .The green points correspond to no FV Bino-slepton loop, considering only the smuons in their mass eigenstates and A 0 = 0 the same as green points in previous figure (fig.7).

Figure 5 :
Figure 5: Values of tan β (left) and A 0 (right) dependence on µ SU SY for which the a µ discrepancy would get solved partially by the LFV contributions (blue), or completely up to 1σ with the restriction M 1 < 1 3 µ susy (black).

Figure 6 :
Figure 6: The values for which the LFV contribution would explain completely the a µ discrepancy within theory and experimental data up to 1σ, considering M 1 < 1 3 µ susy .We show ratio on susy mass parameters M 1 / mS (left) and mS (right), both with respect on µ susy values.

Table 3 :
Sample of parameter sets that solve entirely the muon g − 2 discrepancy, consistent with the experimental bound on BR(τ → µγ), calculated using random values of the parameters given in Table2.For all these sets the LSP is a Bino.