Revisiting R-parity violating interactions as an explanation of the B-physics anomalies

In the last few years, the ratios \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{D^{(*)}}$$\end{document}RD(∗) and of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{K^{(*)}}$$\end{document}RK(∗) have reportedly exhibited significant deviations from the relevant standard model predictions, hinting towards a possible violation of Lepton Flavor Universality and a window to New Physics. We investigate to what extent the inclusion of R-parity violating couplings in the minimal supersymmetric standard model can provide a better fit to the anomalies simultaneously. We perform this analysis employing an approximate, non-abelian \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {G}_f=U(2)_q \times U(2)_\ell $$\end{document}Gf=U(2)q×U(2)ℓ flavour symmetry, which features a natural explanation of the appropriate hierarchy of the R-parity violating couplings. We show that, under the requirement of a supersymmetric spectrum with much heavier left-handed doublet superpartners, our assumption favors a considerable enhancement in the tree-level charged-current \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B \rightarrow D^{(*)} \tau \overline{\nu }$$\end{document}B→D(∗)τν¯, while the the anomalies induced by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b \rightarrow s\ell ^+\ell ^-$$\end{document}b→sℓ+ℓ- receive up to an approximate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30 \%$$\end{document}30% improvement. The consistency with all relevant low-energy constraints is assessed.


Introduction
While exploring the limits of the standard model (SM), many experiments have spent decades looking for processes that do not respect Lepton Flavor Universality (LFU). Until 2014 this fundamental SM feature stood rather steady, but since then several B-physics experiments, started reporting results which conflict with it, to the surprise and excitement of the community. These results are encoded by ratios of branching ratios involving rare b decays and different lepton flavors.
The statistical significance of each anomaly does not exceed the 3σ level, but the overall set is very consistent and the pattern of deviations intriguing. Our analysis will be solely focused on the above mentioned LFU ratios; nevertheless, it is worth mentioning briefly the existence of additional data that exhibit tensions with the SM predictions. In particular, the most notable is a deviation of about 3σ reported [10,11] on the so-called P 5 differential observable of B → K * μμ decays. Even though, given the non-negligible SM uncertainties [12], the P 5 anomaly is also not an individually definite NP signal, it has been pointed out that it has a common model-independent solution with the R K ( * ) anomaly, namely the modification of a single amplitude, the one induced by the semi-leptonic di-muon vector and axial operators [13][14][15].
While the charged-current decays occur at tree-level, the FCNC decays appear at loop level in the SM, rendering a simultaneous explanation of both anomalies a notoriously difficult theoretical endeavor even in the most general effective field theory (EFT) scenarios without some degree of finetuning (see e.g. [16][17][18][19][20][21][22][23][24][25][26][27][28][29] for model-independent studies and  for attempts to cast specific NP models). Moreover, one observes, that the charged-current anomalies require an enhancement in the decay channel that involves the third generation SM fermions, i.e. b and τ , while, as already mentioned, the b → s ¯ anomalies are resolved by assuming purely muonic NP effects.
The above motivate NP scenarios in which the third generation SM fermions is to be treated specially. On the one hand, the special role of the third generation in the radiative corrections of the Higgs boson self-energy and thus the famous problem of naturalness, is evocative of theories that have traditionally addressed this problem, such as Supersymmetry (SUSY). 1 However, it can be easily checked [79] that the R-parity conserving Minimal Supersymmetric Standard Model (MSSM) introduces amplitudes which are orders of magnitude smaller than the ones required to accommodate the one-loop level anomalies, let alone the tree-level ones. The study of the phenomenology of R-parity violating (RPV) interactions has shown instead that either the charged-current or the FCNC anomalies can be individually resolved [61][62][63][64][65], but a united solution which accommodates all other relevant low-energy bounds is impossible.
On the other hand, from a flavour point of view, a non-Abelian, U (2) q × U (2) flavour symmetry, acting on the light generations of SM fermions, is one of the most interesting proposals [80]. Complemented with the dynamical assumption, that the NP sector is coupled preferentially to third generation, this setup can give a consistent picture for all low-energy flavour observables not only at Effective Field Theory (EFT) level [67-70] but also in UV complete models, in which U (2) q × U (2) flavour symmetry appears as a subgroup of a greater gauge sector and emerges at low energies [71,72]. Interestingly enough, the U (2) symmetries were initially proposed in the context of Supersymmetry [80,81] in order to solve the 'flavour' problem of the MSSM, i.e. the abundance of new parameters introduced at the explicit soft SUSY breaking sector and their conspicuous 'near-CKM' alignment that avoids unacceptably large flavor-changing and CP-violating effects.
By invoking an appropriate flavour symmetry, it is possible to link the RPV sector to the origin of masses and mixings, while naturally suppressing the RPV couplings within the experimental bounds [82]. In the current work, we employ the U (2) q ×U (2) flavour symmetry to give a natural justification to the phenomenologically preferable hierarchy of the the RPV couplings. Unlike all previous studies in this framework, we have also taken into account the leptonic current that can be generated by the RPV interactions, besides the usual leptoquark-like current, and how it could affect the relevant amplitudes. A final fit within the natural region of the parameter space reveals certain implications for the SUSY spectrum.
The paper is organised as follows. In Sect. 2 we review the RPV sector and the U (2) q × U (2) flavour symmetry. We suggest a suitable symmetry breaking pattern and rewrite the RPV couplings with the help of the resulting spurions. Subsequently, in Sect. 3 we examine the relevant latest, lowenergy constraints and in Sect. 4 we perform a χ 2 -fit of the anomalies including those constraints. Finally, the discussion of the results is summarized in the Conclusions.

R-parity violating interactions under the U(2) q × U(2) flavour symmetry
The most general renormalizable, R-parity odd superpotential consistent with the gauge symmetry and field content of the MSSM is where there is a summation over the generation indices i, j, k = 1, 2, 3, and summation over gauge indices is understood. One has for example, indices. Gauge invariance enforces antisymmetry of the λ i jk couplings with respect to their first two indices, Gauge invariance also enforces antisymmetry of the λ i jk couplings with respect to their last two indices, Equation (11) involves 48 parameters: 3 dimensionful parameters μ i mixing the charged lepton and down-type Higgs superfields, and 45 dimensionless Yukawa-type couplings divided into 9 λ i jk and 27 λ i jk couplings which break leptonnumber conservation, and 9 λ i jk couplings which break baryon-number conservation. For the sake of simplicity, we shall assume that these couplings are real numbers. The standard motivation for R-parity is that it leads to conserved baryon number and thus ensures proton stability. Nevertheless, according to modern theoretical developments [83], if the MSSM is an effective theory, rapid proton decay can also be induced by higher-dimensional, nonrenormalizable operators suppressed by a scale lower than 10 TeV. Consequently, one needs to impose a stand-alone baryon-number conservation symmetry rather than R-parity.
In this case, the trilinear terms controlled by the couplings λ and λ can be present, while λ = 0.
Regarding the leptonic terms with the coupling λ, they may seem at first irrelevant in the context of the B-physics anomalies, but as we discuss in Sect. 3 they can contribute to the b → sll decays at one-loop level.

Tree-level four-fermion operators
Expanded in standard four-component Dirac notation, the trilinear interaction terms associated with the λ and λ couplings of the RPV superpotential (11) read, respectively, Exchanging Sparticles, one obtains the following fourfermion operators at tree-level: We can simplify the low-energy spectrum by assuming a large mass splitting between the light third generation and the much heavier first two generations of Sfermions and Sleptons, which can be considered as effectively not supersymmetrized. We note, that this simplification is equally wellmotivated in theory [83,84]. In fact, most models of spon-taneous SUSY breaking predict a significantly lighter third generation at the electroweak scale due to the large RG effects proportional to the top Yukawa coupling y t [85]. What is more, as recently shown [65], one of the prominent attributes of SUSY, namely the gauge coupling unification is still preserved despite the decoupling of the first two generations and even in presence of RPV interactions.

Flavour structure
The flavour group we are considering is , under which the superfields transform as: The Higgs fields H u(d) are pure G f singlets. We introduce additional heavier 'flavon' fields, which are charged under the flavour symmetry [81]; in particular a doublet φ a i , a triplet (a 2-index symmetric tensor) S ab i and a singlet (a 2-index antisymmetric tensor) A ab i , where i = q, . The flavour groups U (2) i , are then broken by the vacuum expectation values (VEVs) of these flavon fields, such that, where M is the cut-off of the effective theory. A step-wise breaking is achieved, if we naturally assume , The mass matrices of the charged leptons and the down quarks assume the following form: where 005, the usual quark mass hierarchy is successfully reproduced. In the lepton sector, as it will become clear from the phenomenological analysis of Sect. 3, if one hopes to generate any considerable contribution to the b → s ¯ processes the effective coupling to muons must remain unsuppressed. This translates into a strong breaking of U (2) to U (1) by the doublet flavon VEV, i.e. = 1. The correct lepton mass matrix is then reproduced for S 0.06 and 0.004. The RPV bilinear and trilinear terms in the superpotential can be obtained by appropriately contracting the superfields appearing in Eq. (11) with the flavons. The order of magnitude of the RPV couplings is the governed by and , • λ i jk couplings: • λ i jk couplings: Generic RPV couplings λ i jk (or λ i jk ) can then be decomposed as products of O(1) parameters c i jk (or c i jk ) and the respective and suppression factors.

Constraints from low-energy observables
In this section, we analyse the main experimental constraints on the RPV interactions. The processes of interest are the ones that are affected by contributions of the O(1) couplings λ 323 , λ 233 and λ 333 or at least by the -suppressed couplings λ 223 , λ 232 , λ 323 and λ 332 . These include the R D ( * ) and R K ( * ) ratios, the B s −B s mixing, the B → K ( * ) νν and B → τν decays, the RGE effects in τ → νν and the Z coupling modification for the relevant λ couplings and the τ → νν decays for the only relevant λ coupling. We have explicitly checked that further processes that have been discussed in the bibliography, e.g. the decays B → πνν, and τ → 3μ, do not lead to any relevant constraints in our setup. Of course, all processes involving only the first two generations do not receive any contributions at all.

B → D ( * ) τν
Adding to the SM the RPV contribution generated by the respective operators in Eq. (16), one obtains the effective Lagrangian describing b → c semi-leptonic decays at treelevel, The ratios (1), (2) can then be easily written as: From the above definition and the weighted average of the R D and R D * central values and errors, one observes that the enhancement of R D ( * ) implies rather large λ 333 coupling and λ 323 λ 333 coupling combination (which is also enhanced by V cs /V cb ), while at the same time λ 233 and λ 223 λ 233 must be kept small.

B → K ( * ) ¯
Instead of using the ratios (5) and (6), we will instead regard the NP modification of the Wilson Coefficients C 9 , C 9 , C 10 and C 10 defined as: In our framework, the operatorμ i L γ μ μ i Ld k R γ μ d k R is generated at tree-level (see Eq.
In the following, the upper index μ is implied for the Wilson Coefficients.
Let us examine first the tree-level case, According to Eq. (22), we expect λ 232 λ 233 ∼ 2 q ≈ 0.025 and hence δC 9 ≈ (8 × 10 6 GeV 2 )/m 2 t L . Unless we introduce a significant deviation from the flavour symmetry expectation for the order of magnitude of λ 232 λ 233 , the left-handed Stop has to be very heavy m 2 t L 10 TeV in order to satisfy the experimental results in (29) and (30). We are choosing the latter and generalize the result for all left-handed Sparticles. Indeed, we confirm later on that the limit of decoupled lefthanded Sparticles is favored in our context. We also note, that for this reason the solutions of [64,66] for the R ( * ) K anomaly are not applicable.
Next, the NP effect at one-loop level is The first term corresponds to box diagrams with a W boson andb R in the loop, the second to box diagrams with twob R and the third to a box diagram withν L andt L . The final term is new in our analysis and arises from the diagram in Fig. 1.
The first two terms are unable to explain the R ( * ) K anomaly due to severe constraints from B → K ( * ) νν [63] and the third term is suppressed in our framework due to very heavy left-handed Sparticle mediators. It is left to see at what degree can the final term alleviate the tensions, once we have taken the rest of the constraints into consideration, especially those that concern the leptonic currents (see Sect. 3.7).

B s −B s mixing
The effective Hamiltonian relevant to B = 2 processes is where we have used the notation of Eq. (2.1) in [88] and expressed the Wilson Coefficient of the operator Q L R 1s , accordingly, by performing a Fierz transformation to the operator Q 5 in the so-called 'SUSY basis' [89]. This operator arises at tree-level by mediation of a left-handed Sneutrinō ν L , while the SM-like operator (b L γ μ d i ) 2 and the operator appear first only at one-loop level by various box diagrams [90]. We find which in the limit of decoupled left-handed Sparticles vanishes. Further, if we likewise choose to neglect the box diagrams with left-handed Sparticle mediators, we find The Wilson Coefficients of the rest of the operators do not receive any RPV contributions and finally, we may write the Eqs. (7.25), (7.27) in [88] as: bounds are set for all contributing λ coupling combinations.

Z → ¯ coupling
As investigated in [93,94], the leptonic Z coupling is modified via one-loop diagrams. In our setup, the Z boson decays to a tt pair, which is in turn connected by a virtualb R and eventually turns into a dilepton pair ¯ . One defines the following ratios of vector and axial-vector couplings v and a : In our context 2 and keeping only the term proportional to the top Yukawa, the Eq. (30) in [94] becomes: We denote the NP scale to be roughly mb R ≈ 1 TeV. 2 We match the effective Lagrangian of Eq. (19) with ours (16) by setting C 1 = 1 2 , C 3 = − 1 2 and λ u i j = λ e i j = λ i jk .
A comparison with the measured values [95], v τ v e exp = 0.959 ± 0.029 and a τ a e exp = 1.0019 ± 0.0015, (54) yields bounds on the coupling λ 333 . We note here, that the W coupling is also modified, but the constraints given by LFU violating τ decays on the same couplings are more stringent (see Sect. 3.7).

τ → νν
The purely leptonic operators resulting from the trilinear couplings in (14), affect the Lepton Flavor violating decays τ → eνν and τ → μνν at tree-level through the exchange of a third generation Sleptonτ R [96]. Additionally, the RGE effects driven by the top Yukawa y t interactions contribute also to LFU violation ( [94], via one-loop diagrams involving the leptoquark-like interactions in (15). The NP effects are probed by the ratio: where and At leading order, Eq. (58) becomes: Even though the experimental bounds are very stringent [97], i.e.
with appropriate fine-tuning of the couplings and the masses, one could still recover a non-negligible leptonic current interaction.

Numerical fit and discussion
After the decoupling of left-handed Sparticle related contributions, the low-energy observables discussed above depend solely on the RPV couplings λ 323 , λ 223 , λ 232 , λ 233 , λ 323 , λ 332 and λ 333 and on the masses mb R and mτ R . We have performed a combined fit of these parameters using as input the experimental data reported in Sect. 3 and various SM parameters [95]. The following, conservative, lower bounds for the right-handed Sparticles are imposed: mb R > 800 GeV and mτ R > 400 GeV [95]. We have also assumed that the O(1) parameters c i jk (or c i jk ), which together with the flavour suppression factors constitute the RPV couplings, are restricted by the unitarity bounds of √ 4π . For simplicity, we have assumed Gaussian errors for all the observables. The preferred region of the model parameters x has been determined minimizing the χ 2 distribution: where μ i and σ i are the central values and the 1σ uncertainties of the measured values [O i ] exp , respectively. The best-fit points for the c i jk (or c i jk ) parameters, as well as the flavour parametric scaling and the total value of the relevant RPV couplings are listed in Table 1.
The best-fit points for the masses are mb R = 900 GeV and mτ R = 9 TeV.
The improvement of the best-fit point of the total χ 2 with respect to the SM limit is χ 2 (x SM )χ 2 (x BF ) = 9.43. In Table  2, we show the contributions of the individual summands in the χ 2 .
In Fig. 2 we show the 68% CL and 95% CL regions of the r D ( * ) and δC 9 observables in the (λ 333 , λ 323 ), (λ 333 , mb R ), (λ 323 , mb R ) and (λ 333 , mτ R ) planes, after having fixed the rest of the parameters according to Table 1. The constraints of the other, relevant low-energy observables are presented as exclusion contours at 2σ . The degree of consistency of the best fit-region with the anomalies is illustrated in Fig. 3, Table 2 The individual χ 2 components evaluated at the best-fit point and compared with the respective SM limits Observables where we show the values of the two observables r D ( * ) and δC 9 within the 1σ preferred region ( χ 2 < 2.3).
First and foremost, we observe that our choice of flavor symmetry points towards an alleviation of tensions in the charged-current anomalies similar to the results in [65]. 3 As a matter of fact, it was commonly expected that new phenomena would not show up in tree-level processes, where the effect has to be comparable with the SM value, but rather in FCNC transitions where one has to simply compete against a SM loop-suppressed contribution. In this regard, it is interesting that the flavour structure favors an improvement of the tree-level anomalies rather than the loop-induced ones. Nevertheless, our model prediction can only exhibit a 1σ agreement with the present central value of R D ( * ) in case of a more than 2σ reduction of the ratio a τ /a e . This is, indeed, the main obstacle for obtaining a larger λ 333 , required for a perfect fit. Other bounds, including those from the B s −B s mixing and the B → τν decays, are found to be less significant.
Regarding R K ( * ) , most of the parameter space is excluded due to the bounds from the R τ/μ τ observable . In principle, an appropriate canceling of the second and the third term in (58) can lift this constraint, but considering the previous bounds on λ 333 , we can achieve at most a ∼ 10% cancellation at 2σ deviation from the central value. This still allows for a smaller mass for mτ R and as a result, a slightly better fit for δC 9 .
In total, we see that, by involving the leptonic RPV interactions, we can 'cure' approximately 30% of the b → s + − anomalies, when all other constraints are allowed to deviate by maximum 2σ from their experimental average. This is still an improvement over previous attempts to accommodate  both anomalies within the RPV framework. Please note, that the fit yields very small λ 233 and λ 233 couplings and thus both the second term in (32) and the NP muonic contribution in the denominator of (25) are negligible. This implies that there is no clash between the anomalies themselves, but rather the performance of the fit is limited by the other low-energy constraints.
On a final note, we would like to briefly comment on three, relevant issues. First, it is well-known, that RPV interactions that attempt to address LFU violating effects can also generate neutrino masses at one-loop level [85], The typical way-out is to postulate a direct cancellation between the trilinear coupling A b and the term μ tan β, which is though hard to justify theoretically. Other unrelated NP contributions to neutrino masses, e.g. the standard see-saw mechanism with heavy right-handed neutrinos, could be the cause of the suppression. We stress, here that the possibility of generating a canceling due to loop effects induced by the leptonic trilinear RPV interactions is ruled out due to the direct, tree-level bounds from τ decays on the λ couplings.
The next comment refers to the resulting SUSY mass spectrum. According to the fit, the mass of the right-handed superpartners spans the 1-10 TeV range, with theτ R being significantly heavier than theb R . In fact, only the mass scale ofb R appears to be within the reach of highp T searches in the near future. An even more striking assumption of our setup is the complete decoupling of the left-handed superparters. A theoretical motivation for this scenario lies in the theory of gauge-mediated supersymmetry breaking. As it turns out, LHC observable Squarks can only contribute enough quantum corrections to lift the Higgs mass if the lefthanded SU (2) doublets are much heavier than the SU (2) singlets [98].
Last but not least, as shown in [43], flavour symmetries similar to the one employed here, when gauged and broken at the TeVs scale, can provide a natural mechanism for generating additional LFU violating contributions to b → s + − transitions. A further discussion of this scenario, that would require an enlargement of the field content of the theory, is beyond the scope of this paper.

Conclusions
In this paper we have studied the consistency of the R D ( * ) and R K ( * ) anomalies with all relevant low-energy observables, in the context of RPV interactions controlled by a G f = U (2) q × U (2) flavour symmetry. This particular scenario favors a viable solution of the charged-current anomaly, at least as good as in the generic, effective RPV-SUSY scenario. However, as we have shown, a perfect fit of the anomalous R D ( * ) observable cannot be achieved without a significant modification of the Z boson coupling, occurring at one-loop level. What is more, upon inclusion of the leptonic, trilinear RPV interactions, we were able to generate, simultaneously, a contribution to the b → s + − transitions, which is also limited by the the same Z boson coupling bounds and the tree-level, Lepton Flavor violating τ decays. All in all, the flavour symmetry for natural values of the free parameters, as summarised in Table 1, can provide an explanation of the strength of the RPV interactions, required for a better fit of the B-physics anomalies.