Next Article in Journal
Evaluation of Online Teaching Quality: An Extended Linguistic MAGDM Framework Based on Risk Preferences and Unknown Weight Information
Next Article in Special Issue
Chains of Interacting Solitons
Previous Article in Journal
Synchronous Steady State Solutions of a Symmetric Mixed Cubic-Superlinear Schrödinger System
Previous Article in Special Issue
On ξ-Process for DVCS-Amplitude
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flavour Anomalies in a U(1) SUSY Extension of the SM

by
Alexander Bednyakov
1,2,* and
Alfiia Mukhaeva
1
1
Joint Institute for Nuclear Research, 141980 Dubna, Russia
2
P.N. Lebedev Physical Institute of the Russian Academy of Sciences, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(2), 191; https://doi.org/10.3390/sym13020191
Submission received: 11 January 2021 / Revised: 20 January 2021 / Accepted: 21 January 2021 / Published: 26 January 2021
(This article belongs to the Special Issue Symmetry in Particle Physics II)

Abstract

:
Flavour anomalies have attracted a lot of attention over recent years as they provide unique hints for possible New Physics. Here, we consider a supersymmetric (SUSY) extension of the Standard Model (SM) with an additional anomaly-free gauge U ( 1 ) group. The key feature of our model is the particular choice of non-universal charges to the gauge boson Z , which not only allows a relaxation of the flavour discrepancies but, contrary to previous studies, can reproduce the SM mixing matrices both in the quark and lepton sectors. We pay special attention to the latter and explicitly enumerate all parameters relevant for our calculation in the low-energy effective theory. We find regions in the parameter space that satisfy experimental constraints on meson mixing and LHC Z searches and can alleviate the flavour anomalies. In addition, we also discuss the predictions for lepton-flavour violating decays B + K + μ τ and B + K + e τ .
PACS:
12.60.Jv; 13.30.Ce

1. Introduction

The Standard Model (SM) is a very successful theory, explaining most experimental results. However, there are experimental discrepancies with some of the SM predictions. Among interesting processes there are those that receive special attention in the literature. Most of them are Flavour-Changing-Neutral-Current (FCNC) transitions that in the SM are loop-suppressed and have enhanced sensitivity to the New Physics (NP) effect. The role of such kind of observable is two-fold. After being measured compatible with the SM, it poses severe constraints on the Beyond-the-SM (BSM) scenarios. On the other hand, if there is a tension between the SM and experiment, it stimulates various speculations on possible solutions in the context of BSM models.
In particular, one usually discusses the b s l + l transitions. The LHCb Collaboration has made measurements of B K * μ + μ [1] that deviate from the SM predictions [2]. The Belle Collaboration finds similar results [3]. The main discrepancy is in the angular observable P 5 [4], averaged over the invariant mass q 2 of the lepton pair in the ranges q 2 = [ 4.0 6.0 ] GeV 2 and q 2 = [ 6.0 8.0 ] GeV 2 . Recent LHCb measurement [5]
P 5 [ 4.0 6.0 ] = 0.439 ± 0.111 ± 0.036 , P 5 [ 6.0 8.0 ] = 0.583 ± 0.090 ± 0.030
reports the significance of deviation to be 2.5 σ and 3 σ , respectively (the significance of the discrepancy depends on the assumptions about the theoretical hadronic uncertainties [6]).
Other important observables are ratios [7,8,9] in the dilepton invariant mass-squared range 1 GeV 2 q 2 6 GeV 2   
R K B r ( B + K + μ + μ ) B r ( B + K + e + e ) = 0 . 846 0.054 + 0.060 ( s t a t ) ± 0.021 ( s y s t . ) ,
R K * B r ( B ± K * μ + μ ) B r ( B ± K * e + e ) = 0 . 69 0.07 + 0.11 ( s t a t ) ± 0.05 ( s y s t . )
that test lepton flavour non-universality (see also [10]). They also deviate from the SM predictions of R K S M = 1 and R K * S M = 0.996 . Even though the most recent data seems to be consistent with the SM at 2.5 σ [9], or has large error bars [11], the R K and R K * puzzles still provide very intriguing insights on possible NP.
Finally, let us mention constraints from the meson-mixing. The mass difference of the neutral B s B ¯ s meson system, Δ M s , provides a severe constraint for any NP model aiming at an explanation of the B-physics anomalies. For quite some time the SM value for Δ M s was in perfect agreement with experimental results, see e.g., [12]. Taking however, the most recent lattice inputs, in particular the new average provided by the Flavour Lattice Averaging Group (FLAG) one gets a SM value Δ M s S M = ( 18 . 4 1.2 + 0.7 ) ps 1 [13] considerably above the measurement [14]
Δ M s e x p = ( 17.757 ± 0.021 ) ps 1 .
One can also consider the constraints due to the mixing in the D 0 D ¯ 0 meson system. For example, for the dimensionless ratio of the D 0 mass difference Δ m over the averaged decay width Γ (see, e.g., reference [14,15])) we have (under assumption that there is no CP-violation)
x D Δ m Γ = 0 . 5 0.14 + 0.13 % .
The uncertainties of SM prediction is rather large and as SM value we take a rough estimate provided by Flavio (2.0.0) package [16] x D S M = 0.4 ± 0.4 % .
It is interesting that all the observables for b s l + l , B s B ¯ s , and D 0 D ¯ 0 can be explained by the weak effective theory (WET) with a Hamiltonian of the form (see, e.g., references [16,17]),
H e f f = C i O i + h . c . ,
where relevant operators are ( λ t V t b V t s * )
O 9 b s l l = 4 G F 2 α e 4 π λ t ( s ¯ L γ μ b L ) ( l ¯ γ μ l ) ,     O 10 b s l l = 4 G F 2 α e 4 π λ t ( s ¯ L γ μ b L ) ( l ¯ γ μ γ 5 l ) ,     O 9 b s l l = 4 G F 2 α e 4 π λ t ( s ¯ R γ μ b R ) ( l ¯ γ μ l ) , O 10 b s l l = 4 G F 2 α e 4 π λ t ( s ¯ R γ μ b R ) ( l ¯ γ μ γ 5 l )
for the b s l + l transitions, and
O B s B ¯ s L L = ( s ¯ L γ μ b L ) ( s ¯ L γ μ b L ) ,                     O B s B ¯ s L R = ( s ¯ R γ μ b R ) ( s ¯ L γ μ b L ) , O B s B ¯ s R R = ( s ¯ R γ μ b R ) ( s ¯ R γ μ b R )
together with
O D 0 D 0 ¯ L L = ( u ¯ L γ μ c L ) ( u ¯ L γ μ c L ) , O D 0 D 0 ¯ L R = ( u ¯ R γ μ c R ) ( u ¯ L γ μ c L ) , O D 0 D 0 ¯ R R = ( u ¯ R γ μ c R ) ( u ¯ R γ μ c R ) ,
for the above-mentioned meson mixing. For further convenience we also list here the operators that give rise to lepton-flavour violating decays of B mesons, which were studied experimentally by BaBar Collaboration [18,19]:   
O 9 b s μ τ = 4 G F 2 α e 4 π λ t ( s ¯ L γ μ b L ) ( τ ¯ γ μ μ ) , O 10 b s μ τ = 4 G F 2 α e 4 π λ t ( s ¯ L γ μ b L ) ( τ ¯ γ μ γ 5 μ ) , O 9 b s μ τ = 4 G F 2 α e 4 π λ t ( s ¯ R γ μ b R ) ( τ ¯ γ μ μ ) , O 10 b s μ τ = 4 G F 2 α e 4 π λ t ( s ¯ R γ μ b R ) ( τ ¯ γ μ γ 5 μ ) .
The analyses before Moriond 2019 [17] and after Moriond 2019 [20] show that several patterns of NP contributions explain the discrepancies significantly better than the SM. In all cases, there should be a sizable negative contribution to C 9 b s μ μ .
A way to address the b s l + l anomalies is to introduce a Z gauge boson which couples to muons and down-type quarks. For instance, U ( 1 ) μ τ gauge symmetry is employed to control the flavor dependent couplings of the Z boson [21]. It is shown in [22] that b s l + l anomalies can be successfully explained in models with a Z boson.
In this work, we consider a non-universal U ( 1 ) gauge extension of the Minimal Supersymmetric Standard Model (MSSM), with family dependent couplings to quarks and leptons [23,24]. Such an U ( 1 ) could emerge in GUT, superstring constructions or dynamical electroweak breaking theories. We take this U ( 1 ) extended supersymmetric model as a simple extension of MSSM, allowing more flexibility in model parameters.
We analyse this model in the context of all three SM fermion families. This allows us to explicitly discuss both the Cabibbo-Kobayashi-Maskawa (CKM) and Pontecorvo–Maki–Nakagawa–Sakata (PMNS) matrices together with additional mixing allowed in the model. The purpose of this paper is to demonstrate the existence of scenarios that can relax the tension in flavour anomalies and discuss some interesting manifestations of the model. A more detailed analysis of the expected phenomenology in a wider parameter space is delegated to future work.
The rest of this paper is organized as follows. The model is introduced in Section 2 and the relevant Wilson coefficients are computed in Section 3. In Section 4 we discuss the results of our fit and study model predictions. Section 5 is devoted to our conclusions.

2. Model Description

Let us briefly describe our model and set up the notation. We consider a U ( 1 ) extension of the MSSM similar to that of reference [24]. In addition to the chiral multiplets of the MSSM, we also introduce a singlet (strictly speaking, the field S is singlet only w.r.t. the SM gauge group) superfield S, which allows one to break U ( 1 ) spontaneously and generate mass for the corresponding Z boson. To account for the massive neutrinos three right-handed chiral superfields ν 1 , 2 , 3 c are also introduced.
The charges of the additional gauge group are not universal and, thus, potentially allow one to accommodate for the flavour anomalies discussed in the literature. As usual, the requirement that there should be no gauge anomalies in the model, imposes important restrictions on the charges. While in the non-supersymmetric U ( 1 ) extensions one usually takes into account only the SM fermions (see, e.g., [25,26]), with SUSY we have plenty of half-spin superpartners, which can also contribute to the gauge anomalies.
It is known that models with charge assignments ( B L ) i and ( L i L j ) , where B i ( L i ) are baryon (lepton) numbers of the i generations are free from anomalies. Due to this, in reference [24] the model based on Q = a ( B L ) 3 + b ( L μ L τ ) with a = 3 / 2 and b = 2 was considered. A key observation of reference [24] was the fact that the Higgs superfield H d and the chiral lepton superfields L i have the same SM charges, so one can, e.g., switch the U ( 1 ) charges of L 3 and H d without spoiling the anomaly cancellation. As a consequence, the contribution of the left-handed (LH) taus to the anomalies is replaced by that of higgsinos. Moreover, since the right-handed (RH) neutrinos also carry the corresponding lepton numbers, we can switch the U ( 1 ) charges of ν c and S.
Our initial motivation was to extend the study of reference [24] and analyse the effect of a more general mixing in both the quark and lepton sectors. However, one can show that the lepton mass matrices compatible with U ( 1 ) charges of [24] turn out to be block-diagonal. Due to this, one has unbroken global U ( 1 ) L 2 symmetry corresponding to the lepton number L 2 and can not accommodate for the PMNS mixing in the neutrino sector. Parameter counting based on broken flavour symmetries confirms this statement.
To circumvent this problem we modify the initial charge assignment to allow more general Yukawa textures in the lepton sector. We started with Q = a ( B L ) 3 + b ( L 2 L 3 ) + c ( L 1 L 2 ) and made the substitutions L 3 H d , ν 3 c S . Among possible solutions we have chosen the one with a = 3 , b = 2 , c = 1 :
Q = 3 B 3 L 1 L 2 H d + S + E 3 c ,
where B 3 ( L i ) assigns 1 / 3 (1) to the top-quark (i generation lepton) superfields (both LH and RH (we use LH charge-conjugated superfields to account for the RH particles, so in (11) one can write L i = E L i + ν L i E i c ν i c ), while H d , S and E 3 c are equal to one for the higgs H d , the singlet S and the right-handed tau E 3 c superfields, respectively, and zero otherwise. The quantum numbers corresponding to (11) can be found in Table 1.
One can see that (modulo 1 / 2 , which can be absorbed into redefinition of g E ) in the quark sector we have the same charges as in reference [24]. However, the charges of H d , S and E 3 c are flipped. In addition, the fields carrying L 1 are also coupled to Z . The corresponding R-parity conserving superpotential is given by
W = i , j = 1 , 2 Y u i j Q i H u U j c + Y u 33 Q 3 H u U 3 c ( Q 3 H d ) ( Y d 31 D 1 c + Y d 32 D 2 c ) + i , j = 1 , 2 Y ν i j L i H u ν j c + M 3 ν ν 3 c ν 3 c + Y ν 33 L 3 H u ν 3 c ( L 3 H d ) Y e 31 E 1 c + Y e 32 E 2 c + Y e 33 E c 3 + λ s S H u H d ,
where LH chiral quark (lepton) superfields are denoted by Q i ( L i ), and U c , D c , E c , and ν C correspond to up-quark, down-quark, charged-lepton and neutrino RH fields, respectively. Since ν 3 c is not charged we also add a Majorana mass M 3 ν . The two higgs superfields H u and H d are coupled to the singlet S, the vacuum expectation value (VEV) S of which gives rise to the effective μ parameter and provide a solution to the μ problem.
The gauge field Z couples to quarks and leptons as
L g E Z α b ¯ γ α b + t ¯ γ α t g E Z α i = 1 , 2 [ l ¯ i L γ α l i L + ν ¯ i L γ α ν i L ] + ν ¯ i R γ α ν i R i = 1 , 3 l ¯ i R γ α l i R .
Here g E corresponds to the U ( 1 ) gauge coupling and all fermions are written in the weak basis. In what follows we assume that the Z Z mixing is negligible [24]. The pattern (13) allows one to evade the constraints on the Z production cross-section from the LHC dilepton searches (see, e.g., reference [27]). The weak eigenstates in Equation (13) have to be rewritten in terms of mass eigenstates, which originate from the diagonalization of the mass matrices. The latter are related by spontaneous symmetry breaking to the (effective) Yukawa couplings. However, one can see that certain Yukawa interactions are not allowed at the tree level. Nevertheless, it is possible to consider the non-holomorphic soft SUSY-breaking terms [28]
L s o f t n h = i = 1 2 j = 1 3 C E i j ( H u * l ˜ i ) E ˜ j c + C D 33 H u * q ˜ 3 d ˜ 3 c + H u * i , j = 1 , 2 C D i j q ˜ i d ˜ j c + H d * q ˜ 1 C U 13 + q ˜ 2 C U 23 u ˜ 3 c + H d * l ˜ 1 C ν 13 + l ˜ 2 C ν 23 ν ˜ 3 c + h . c . .
which are not forbidden by the U ( 1 ) gauge symmetry and in which scalar superpartners of the SM fermions couple to the “wrong” Higgs doublets. Given (14), additional contribution to the fermion mass matrices are generated [28], which similar to reference [24] we denote by κ f i j v α / 2 , where v α is VEV of the “wrong” doublet for a fermion f, and κ i j incorporates the loop-induced correction.
It is convenient to combine the tree-level terms and the non-holomorphic contributions into effective Dirac mass matrices:
2 m u = Y u 11 v u Y u 12 v u κ u 13 v d Y u 21 v u Y u 22 v u κ u 23 v d 0 0 Y u 33 v u , 2 m d = κ d 11 v u κ d 12 v u 0 κ d 21 v u κ d 22 v u 0 Y d 31 v d Y d 32 v d κ d 33 v u ,
for quarks and
2 m e = κ e 11 v u κ e 12 v u κ e 13 v u κ e 21 v u κ e 22 v u κ e 23 v u Y e 31 v d Y e 32 v d Y e 33 v d , 2 m ν = Y ν 11 v u Y ν 12 v u κ ν 13 v d Y ν 21 v u Y ν 22 v u κ ν 23 v d 0 0 Y ν 33 v u
for leptons. Diagonalizing the matrices (in the case M 3 ν = 0 ) by (bi)unitary transformations and rewriting (13) in the mass basis we generate the tree-level FCNC transitions, governed by the mixing matrices. One can see that there is enough freedom to account for CKM and PMNS mixing. Indeed, let us consider an effective low-energy model with all the SUSY partners but the Z boson and Higgs fields integrated out.
We can count the number of independent “physical” parameters in the flavour sector of our effective Z model by the following reasoning. The SM gauge group respects the U ( 3 ) Q × U ( 3 ) L × U ( 3 ) D × U ( 3 ) U × U ( 3 ) E × U ( 3 ) ν flavour symmetry corresponding to independent rotations of quark (Q) and lepton (L) LH doublets, RH up-type (U) and down-type (D) quarks, and RH charged (E) and neutral ( ν ) leptons. The U ( 1 ) coupling to fermions (13) breaks this symmetry down to
U f l a v o u r = U ( 2 ) Q 12 · U ( 1 ) Q 3 × U ( 2 ) U 12 · U ( 1 ) U 3 × U ( 2 ) D 12 · U ( 1 ) D 3 × U ( 2 ) L 12 · U ( 1 ) L 3 · U ( 3 ) E × U ( 2 ) ν 12 · U ( 1 ) ν 3 ,
where, e.g., U ( 2 ) L 12 corresponds to U ( 2 ) rotations of the first two generations of LH doublets. In turn, the introduction of the effective Yukawa interactions, which give rise to mass matrices (15) and (16), breaks U f l a v o u r down to
U f l a v o u r U B ( 1 ) × U L ( 1 ) , M 3 ν = 0 U B ( 1 ) , M 3 ν 0
The broken generators of (17) can be used to get rid of the “unphysical” parameters of the low-energy model. Indeed, the latter is given by N p h y s = N t o t N b r o k e n , where N t o t is the total number of parameters in the effective Yukawa couplings, and N b r o k e n denotes the number of broken generators.
The counting goes as follows. In the quark sector we have 3 × 3 matrices m u and m d , which depend on 14 complex parameters. In the lepton sector there are 9 complex parameters in m e and m ν depends on 7 complex numbers. The is also a complex Majorana mass M 3 ν . One can write (each U ( n ) factor in (17) depends on n ( n 1 ) / 2 angles and n ( n + 1 ) / 2 phases)
N t o t q u a r k = 14 R e + 14 I m , N t o t l e p t o n = 16 R e + 16 I m + ( 1 R e + 1 I m ) M 3 ν 0 ,                    
N b r o k e n q u a r k = 3 a n g l e s + ( 12 1 ) p h a s e s , N b r o k e n l e p t o n = 5 a n g l e s + 14 p h a s e s ( 1 p h a s e s ) M 3 ν = 0
where we indicate the parameters corresponding to Majorana M 3 ν , and negative contribution accounts for unbroken U ( 1 ) B and U ( 1 ) L . As a consequence, we have
N p h y s q u a r k = 11 R e + 3 I m , N p h y s l e p t o n = 11 R e + 3 I m + ( 1 R e ) M 3 ν 0 .
One can see that in the quark sector there are four additional real parameters besides 6 quark masses, 3 CKM angles and 1 CKM phase. To simplify the analysis and have some symmetry between quarks and leptons, we assume that Majorana mass M 3 ν = 0 . In this case, 10 out of 14 parameters correspond to 6 lepton masses, 3 angles and 1 CP-violating phase in PMNS. Again, we are left with 4 additional parameters in the lepton sector.
It is convenient to incorporate the new parameters as angles and phases in the mixing matrices of quarks and leptons. In our study we relate the weak and the mass eigenstates by means of the mixing matrices (we use Dirac spinors here):
u L = V u L U L , u R = V u R U R , d L = V d L D L , d R = V d R D R ,
l L = V e L E L , l R = V e R E R , ν L = V v L N L , ν R = V v R N R .
where the left-hand side (LHS) corresponds to the weak basis, and the fields in the right-hand side (RHS) are in the mass basis. Since we want to reproduce the CKM and PMNS matrices, one should require that
V C K M = V u L V d L , U P M N S = V e L V ν L .
In addition, given the diagonal fermion mass matrices D u , D d , D e and D ν in the mass basis, we have to reproduce textures in the weak basis, i.e.,
m u = V u L D u V u R , m d = V d L D d V d R , m e = V e L D e V e R , m ν = V ν L D ν V ν R
should have the form (15) and (16). All new mixing parameters counted in (21) are introduced as four angles and four phases entering
V d L = c 13 d 0 e i ϕ 13 s 13 d e i ( ϕ 13 ϕ 23 ) s 13 d s 23 d c 23 d e i ϕ 23 c 13 d s 23 d e i ϕ 13 s 13 d c 23 d e i ϕ 23 s 23 d c 13 d c 23 d
and
V e L = c 13 e 0 e i χ 13 s 13 e e i ( χ 13 χ 23 ) s 13 e s 23 e c 23 e e i χ 23 c 13 e s 23 e e i χ 13 s 13 e c 23 e e i χ 23 s 23 e c 13 e c 23 e .
In (26) we use s 23 d sin ( α 23 ) , etc., while in (27) s 23 e sin ( β 23 ) , etc. Without loss of generality we can set V e R = 1 . The remaining mixing matrices can be parametrized in the same way as V d L and V e L , but all the corresponding angles and phases are determined from the conditions (24) and (25).
Plugging (22) and (23) into (13), we obtain the Z couplings to the mass eigenstates
L + g E Z α i , j = 1 , 3 V u R , i 3 V u R , j 3 * ( U ¯ i R γ α U j R ) + V u L , i 3 V u L , j 3 * ( U ¯ i L γ α U j L ) + g E Z α i , j = 1 , 3 V d R , i 3 V d R , j 3 * ( D ¯ i R γ α D j R ) + V d L , i 3 V d L , j 3 * ( D ¯ i L γ α D j L ) g E Z α i = 1 , 3 E ¯ i γ α E i + N ¯ i γ α N i V e L , 3 i * V e L , 3 j ( E ¯ i L γ α E j L ) + g E Z α i , j = 1 , 3 V ν L , 3 i * V ν L , 3 j ( N ¯ i L γ α N j L ) + V ν R , 3 i * V ν R , 3 j ( N ¯ i R γ α N j R ) .
One can see that all the Z couplings to the SM fermions are determined either by the third column of the quark mixing matrices
V d L , i 3 = s 13 d e i ϕ 13 , c 13 d s 23 d e i ϕ 23 , c 13 c 23 ,
V d R , i 3 = m b m s s 13 d e i ϕ 13 , m b m d c 13 d s 23 d e i ϕ 23 , m s m d c 13 d c 23 d m d 2 ( m b 2 s 23 2 + m s 2 c 23 2 ) c 13 2 + m b 2 m s 2 s 13 2 ,
V u L , i 3 = V ˜ u d , V ˜ c d , V ˜ t d , V u R , i 3 = m u V ˜ u d , m c V ˜ c d , m t V ˜ t d m t 2 | V ˜ t d | 2 + m c 2 | V ˜ c d | 2 + m u 2 | V ˜ u d | 2
or by the third row of the leptonic ones
V e L , 3 i = s 13 e e i χ 13 , c 13 e s 23 e e i χ 23 , c 13 e c 23 e ,
V ν L , 3 i = U ˜ l 1 , U ˜ l 2 , U ˜ l 3 , V ν R , 3 i = m ν 1 U ˜ l 1 , m ν 2 U ˜ l 2 , m ν 3 U ˜ l 3 m ν 3 2 | V ˜ l 3 | 2 + m ν 2 2 | V ˜ l 2 | 2 + m ν 1 2 | V ˜ l 1 | 2 .
For convenience, we introduce the following shorthand notation (the similarity of V ˜ q d and U ˜ l i is due to convenient parametrization of V d L and V e L )
V ˜ q d c d 13 ( V q b c 23 d V q s s 23 d e i ϕ 23 ) V q d s 13 d e i ϕ 13 , q = { u , c , t }
U ˜ l i c e 13 ( U τ i c 23 e U μ i s 23 e e i χ 23 ) U e i s 13 e e i χ 13 , i = { 1 , 2 , 3 } ,
where V u i d j and U e i j are matrix elements of CKM and PMNS, respectively.
Even in this effective model the total number of additional parameters is quite large. To simplify our analysis, we neglect the CP violation in the CKM and PMNS matrices and consider the CP-conserving cases with ϕ 13 = ϕ 23 = χ 13 = χ 23 = 0 . Moreover, it is clear that α 13 and β 13 induce FCNC involving first generation of the SM fermions and subject to tight constraints. Due to this, our main goal is to study the scenario with α 13 = β 13 = 0 . Nevertheless, we also analyse the case with β 23 = 0 , β 13 0 .
To summarize, we study a SUSY-motivated U ( 1 ) extension of the SM with the following set of parameters: the U ( 1 ) coupling ( g E ), the mass of the Z boson ( M Z ) and the two angles, either ( α 23 , β 23 ) or ( α 23 , β 13 ) . In the subsequent section various constraints on this parameter space are obtained and interesting signatures are considered.

3. Wilson Coefficients

Tree-level exchange of massive Z with M Z m b contributes to Wilson coefficients of relevant effective operators. The latter can be deduced from the Lagrangian (28). In our model we have:
C 9 ( ) b s l l = N b s · sin 2 α 23 · B b s L ( R ) B l l V , C 10 ( ) b s l l = N b s · sin 2 α 23 · B b s L ( R ) B l l A ,
where the normalization factor N b s 4 2 π 2 e 2 G F V t b V t s * g E M Z 2 depends on the Fermi constant G F , electric charge e, and the CKM matrix elements V t b V t s * . The factors
B b s L = 1 , B b s R = B b s L m b m s m b 2 sin 2 α 23 + m s 2 cos 2 α 23
correspond to LH and RH quark currents, respectively. One can see that for small α 23 m s / m b 0.02 the factor B b s R scales as m b / m s 40 . As a consequence, C 9 , 10 can be enhanced w.r.t C 9 , 10 . However, for larger values 0.15 α 23 π / 2 we have B b s R 1 .
The leptonic vector and axial-vector currents give rise to the factors
B l l V = 1 , l l = e e 3 + cos 2 β 23 4 , l l = μ μ 3 cos 2 β 23 4 , l l = τ τ sin 2 β 23 4 , l l = τ μ , B l l A = 0 , l l = e e 1 cos 2 β 23 4 , l l = μ μ 1 + cos 2 β 23 4 l l = τ τ sin 2 β 23 4 , l l = τ μ
for the case β 13 = 0 . It is worth noting the hierarchy
1 = B e e V B μ μ V 3 4 B τ τ V 1 2 B τ τ A 1 4 B μ μ A B e e A = 0 for π / 0 β 23 π / 4 , 1 = B e e V B τ τ V 3 4 B μ μ V 1 2 B μ μ A 1 4 B τ τ A B e e A = 0 for π / 4 β 23 π / 2 .
The leptonic factors for the case β 23 = 0 can be easily obtained from that of β 23 = 0 by the substitution β 23 β 13 , e μ .
Heavy Z in the spectrum also affects meson mixing via the Wilson coefficients
C B s B s ¯ L L = g E M Z 2 sin 2 2 α 23 8 ,
C B s B s ¯ L R = 2 C B s B s ¯ L L · B b s R ,
C B s B s ¯ R R = C B s B s ¯ L L · ( B b s R ) 2
in the case of B s B ¯ s , and for D 0 D ¯ 0 via
C D 0 D 0 ¯ L L = g E M Z 2 V ˜ c d 2 V ˜ u d 2 2
C D 0 D 0 ¯ L R = 2 C D 0 D 0 ¯ L L · B u c R
C D 0 D 0 ¯ R R = C D 0 D 0 ¯ L L · ( B u c R ) 2 .
In the considered limit V ˜ q d = V q b cos α 23 V q s sin α 23 , and for convenience we introduce
B u c R m c m u m t 2 V ˜ t d 2 + m c 2 V ˜ c d 2 + m u 2 V ˜ u d .
Due to the hierarchy in masses of up-type quarks and the structure of CKM, one can deduce that B u c R 10 3 , so the contribution of the RH currents are suppressed.
Wilson coefficients were evaluated by means of Flavio (2.0.0) [16] package at the scale 1 TeV scale. Given the Wilson coefficients we can calculate flavour observables. Let us now discuss the analysis of the predictions in more detail.

4. Results of the Study

Taking into account all the above-mentioned experimental constraints, we find the allowed parameter space in the U ( 1 ) model. We consider the log-likelihood function incorporating the experimental measurements, and fit the three parameters M Z / g E , α 23 , and either β 23 or β 13 . We assume that π / 2 α 23 0 , since negative angles give rise to the negative contribution to C 9 favoured by model-independent studies (see, e.g., [17]).
Figure 1 demonstrates how one- and three-sigma regions for different constraints overlap. The blue, green and purple regions depict bounds due to R K * , P 5 , and Δ M s , respectively. In spite of large uncertainty, the constraint on the D 0 D ¯ 0 mixing restricts | α 23 | 0.20 and is not shown.
The best-fit points are obtained by means of Iminuit package [29] that utilizes the MINOS algorithm [30]. One can see that the quark mixing angle is tightly bounded near α 23 0.05 ( β 13 = 0 ) and α 23 0.19 ( β 23 = 0 ). The minimum of the log-likelihood function correspond to our benchmark points (BMPs):
BMP 1 : α 23 = 0.054 ( 4 ) , β 23 = 0.66 ( 11 ) , M Z / g E = 10.55 ( 1.60 ) TeV ,
BMP 2 : α 23 = 0.190 ( 4 ) , β 13 = 1.04 ( 11 ) , M Z / g E = 17.0 ( 1.6 ) TeV ,
where we also indicate 1 σ intervals obtained from the fit.
It is clear that it is the bound on Δ M B s coming from the B s B ¯ s mixing that severely restricts α 23 . The recent lattice results [31] imply that the SM contribution to Δ M B s is slightly larger than the experimental central value. In our model we have operators involving RH currents that can alleviate the difference. Indeed, adopting the formula from [13] to our case
Δ M s S M + N P Δ M s S M 1 + 50 5 TeV M Z / g E 2 sin 2 ( 2 α 23 ) 1 + ( B b s R ) 2 9 B b s R
we see that for 1 + ( B b s R ) 2 9 B b s R < 0 it is possible to achieve Δ M s S M + N P Δ M s S M Δ M s e x p Δ M s S M 0.96 .
Given the formulas for R K and R K * (see Equations (14) and (17) of reference [32]), one can show that in the considered scenarios the differences ( B l l V A B l l V B l l A )
( B μ μ V A B e e V A ) β 13 = 0 = cos 2 β 23 1 2 ,
( B μ μ V A B e e V A ) β 23 = 0 = 1 cos 2 β 13 2 ,
control the sign of the corrections to R K and R K * . From Equation (38) we deduce that for α 23 < 0 either R K 1 and R K * R K ( β 13 = 0 ) or R K 1 , R K * R K ( β 23 = 0 ). Clearly, the first scenario is disfavored by current experimental bounds. Nevertheless, we analyse both of them and examine whether the flavour anomalies can be accounted for, or at least, relaxed in comparison with the SM.
To get separate bounds on g E and M Z we consider constraints on Z production at the LHC (see also, reference [33]). We take into account the results of LHC searches p p Z l l [27,34] and p p Z t t ¯ [35]. To estimate relevant cross-section we use the following simplified expression:
σ ( p p Z f ) = K · 2 J + 1 · q , q Γ ( X q q ) M Z · C q q ( s , M Z ) s · Br ( Z f ) ,
where J = 1 is the spin of Z , and the final states are either f = e e ¯ , μ μ ¯ or f = t t ¯ . The dimensionless factors C q q   
C q q = 4 π 9 M Z 2 / s 1 d x x q x q ¯ M Z 2 / ( s x ) + q M Z 2 / ( s x ) q ¯ x
are given in terms of quark distributions q ( x ) and are evaluated at the scale Q 2 = M Z 2 . To estimate the cross-section we use the MSTW2008NLO [36] set of PDFs and the code ManeParse [37]. For the dilepton final state we take into account the leading QCD corrections by introducing a K-factor K 1.26 (see e.g., reference [38]). In our study we assume that all SUSY particles coupled to Z are much heavier that the boson, so Z can only decay via the interactions given in (28).
In Figure 2 we present the regions in the ( g E , M Z ) plane for fixed values of the angles corresponding to our BMPs, which are excluded by recent searches [27] p p X l l ¯ . The constraints due to p p X t t ¯ turn out to be much weaker. It is worth noting that larger values of the quark mixing angle α 23 give rise to a more restrictive bound on the parameter space. On the contrary, non-zero values of β 23 ( 13 ) can reduce the coupling Z μ μ ( e e ) and, thus, relax the corresponding constraints.
Let us also mention that the value of the U ( 1 ) coupling at the weak scale can be bounded from above by the requirement that there should be no Landau pole up to Λ = 10 19 GeV (Planck scale). Given the one-loop beta-function
β ( g E ) = g E 3 16 π 2 i Q i 2 = 3 g E 2 2 π 2 g E ( μ ) 3 π 2 ln Λ μ 1 / 2 ,
one can deduce that (see also [24])
g E ( 1 TeV ) 0.3 ,
which motivates our choice of the g E upper limit in Figure 2.
One can see that BMP1 with M Z / g E 10.55 TeV lies just at the boundary of the excluded region, which make this point very unnatural. On the contrary, there is some freedom in the choice of M Z and g E for BMP2 with M Z / g E 17 TeV, we assume that our BMP2 has minimal possible value of M Z = 2300  GeV with the corresponding g E = 0.135 .
Table 2 presents the model predictions for our benchmark points. The uncertainties given in Table 2 for the observables are related to the variation of parameters and are calculated using the Flavio (2.0.0) [16] package. For the LHC cross-sections at 13 and 14 TeV we give our estimates of the upper bounds since we neglect possible decays of Z to non-SM particles.
The estimated uncertainties for Δ M B s are rather big, which indicates that the scenarios are fine-tuned. Nevertheless, we see that both our benchmark points can easily account for Δ M B s , and P 5 in the central q 2 region. As we discussed earlier, we have different predictions for R K and R K * in two considered cases. While the tension in R K * can be alleviated, R K > 1 suggests that BMP1 is excluded. For BMP2 the U ( 1 ) model can accommodate smaller values of R K . Yet we predict R K * > R K , so future measurements of R K and R K * can either favour or exclude the scenario.
In the Table 2 we also add our predictions for B R ( B + K + μ τ ) and B R ( B + K + e τ ) . The BaBar Collaboration give the following 90% C.L. upper limits [19]:
B r ( B + K + τ μ + ) < 4.5 × 10 5 , B r ( B + K + τ + μ ) < 2.8 × 10 5 ,
B r ( B + K + τ e + ) < 4.3 × 10 5 , B r ( B + K + τ + e ) < 1.5 × 10 5 .
The computed values are far below current experimental bounds (56). Unfortunately, the model does not produce B R ( B + K + μ ( e ) τ ) sufficiently large to be observed at Belle II [39] (5 × 10 7 ).

5. Conclusions

We investigated the possibility to accommodate for the well-known flavour anomalies in the context of a U ( 1 ) supersymmetric extension of the SM with non-holomorphic soft terms. Contrary to previous studies, we extend the analysis to account not only for the quark and lepton masses together with the CKM matrix but also for the PMNS mixing of neutrinos. Moreover, we enumerated all relevant mixing parameters in quark and lepton sectors and, in addition to encoded in CKM and PMNS matrices, introduce four CP-conserving angles and four CP-violating phases.
To simplify our phenomenological analysis we restricted ourselves to the CP-conserving scenario with two additional angles: α 23 in the quark sector and either β 23 or β 13 in the lepton sector. We considered four-dimensional parameter space ( g E , M Z , α 23 , β 23 ( 13 ) ) , computed relevant Wilson coefficients, and took into account experimental flavour data to fit the ratio M Z / g E of the Z mass and the U ( 1 ) coupling g E together with the additional fermion mixing angles.
Our analysis demonstrated that the most restrictive bound comes from the B s B ¯ s mixing. Nevertheless, due to the presence of right-handed operators, it is possible to relax the tension between the SM prediction and the experimental value. Another interesting observation is the hierarchy between R K and R K * . In the considered scenarios we either have R K > R K * ( β 13 = 0 ) or R K < R K * ( β 23 = 0 ). Unfortunately, the case β 13 = 0 give rise to R K 1 and it is hard to accommodate the anomaly. We took into account the results of LHC searches and found viable points with M Z 2300 GeV and g E 0.3 . We also considered lepton-flavour violating semileptonic decays of B mesons predicted in the model. Unfortunately, the branching ratios Br ( B + K + μ τ ) 10 8 and Br ( B + K + e τ ) 10 9 are far below current and future experimental limits.
Of course, our analysis is far from complete and we plan to extend it, e.g., by computing more observables and considering scenarios with CP-violation. We also expect that the case when both β 13 and β 23 are non-zero has richer phenomenology.
Finally, let us mention that the parametrization of the effective Yukawa matrices in terms of additional angles and phases proposed in the paper can be used in comprehensive analysis of the full model with the SARAH/SPheno toolkit [40].

Author Contributions

A.B. and A.M. contributed equally to the manuscript. Both authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Grant of the Russian Federation Government, Agreement No. 14.W03.31.0026 from 15.02.2018.

Acknowledgments

The authors are grateful to D. Straub for correspondence regarding the Flavio code.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aaij, R.; Adeva, B.; Adinolfi, M.; Adrover, C.; Affolder, A.; Ajaltouni, Z.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Measurement of Form-Factor-Independent Observables in the Decay B0K*0μ+μ. Phys. Rev. Lett. 2013, 111, 191801. [Google Scholar] [CrossRef] [Green Version]
  2. Egede, U.; Hurth, T.; Matias, J.; Ramon, M.; Reece, W. New observables in the decay mode B¯dK¯*0l+l. JHEP 2008, 11, 032. [Google Scholar] [CrossRef] [Green Version]
  3. Abdesselam, A.; Adachi, I.; Adamczyk, K.; Aihara, H.; Al Said, S.; Arinstein, K.; Arita, Y.; Asner, D.M.; Aso, T.; Atmacan, H.; et al. Angular analysis of B0K*(892)0+. In Proceedings of the LHC Ski 2016: A First Discussion of 13 TeV Results, Obergurgl, Austria, 10–15 April 2016. [Google Scholar]
  4. Descotes-Genon, S.; Hurth, T.; Matias, J.; Virto, J. Optimizing the basis of BK*ll observables in the full kinematic range. JHEP 2013, 5, 137. [Google Scholar] [CrossRef] [Green Version]
  5. Aaij, R.; Beteta, C.A.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H.; Aidala, C.A.; Aiola, S.; Ajaltouni, Z.; Akar, S.; et al. Measurement of CP-Averaged Observables in the B0K*0μ+μ Decay. Phys. Rev. Lett. 2020, 125, 011802. [Google Scholar] [CrossRef] [PubMed]
  6. Descotes-Genon, S.; Hofer, L.; Matias, J.; Virto, J. On the impact of power corrections in the prediction of BK*μ+μ observables. JHEP 2014, 12, 125. [Google Scholar] [CrossRef] [Green Version]
  7. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Test of lepton universality using B+K++ decays. Phys. Rev. Lett. 2014, 113, 151601. [Google Scholar] [CrossRef] [Green Version]
  8. Aaij, R.; Adeva, B.; Adinolfi, M.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; Alkhazov, G.; et al. Test of lepton universality with B0K*0+ decays. JHEP 2017, 8, 055. [Google Scholar] [CrossRef] [Green Version]
  9. Aaij, R.; Beteta, C.A.; Adeva, B.; Adinolfi, M.; Aidala, C.A.; Ajaltouni, Z.; Akar, S.; Albicocco, P.; Albrecht, J.; Alessio, F.; et al. Search for lepton-universality violation in B+K++ decays. Phys. Rev. Lett. 2019, 122, 191801. [Google Scholar] [CrossRef] [Green Version]
  10. Bifani, S.; Descotes-Genon, S.; Romero Vidal, A.; Schune, M.H. Review of Lepton Universality tests in B decays. J. Phys. G 2019, 46, 023001. [Google Scholar] [CrossRef] [Green Version]
  11. Abdesselam, A.; Adachi, I.; Adamczyk, K.; Ahn, J.K.; Aihara, H.; Said, S.A.; Arinstein, K.; Arita, Y.; Asner, D.M.; Atmacan, H.; et al. Test of lepton flavor universality in BK*+ decays at Belle. arXiv 2019, arXiv:1904.02440. [Google Scholar]
  12. Artuso, M.; Borissov, G.; Lenz, A. CP violation in the B s 0 system. Rev. Mod. Phys. 2016, 88, 045002. [Google Scholar] [CrossRef] [Green Version]
  13. Di Luzio, L.; Kirk, M.; Lenz, A.; Rauh, T. ΔMs theory precision confronts flavour anomalies. JHEP 2019, 12, 009. [Google Scholar] [CrossRef] [Green Version]
  14. Amhis, Y.S.; Banerjee, S.; Ben-Haim, E.; Bona, M.; Bozek, A.; Bozzi, C.; Brodzicka, J.; Chrzaszcz, M.J.; Dingfelder, J.C.; Egede, U.; et al. Averages of b-hadron, c-hadron, and τ-lepton properties as of 2018. arXiv 2019, arXiv:1909.12524. [Google Scholar]
  15. Zyla, P.; Barnett, R.M.; Beringer, J.; Dahl, O.; Dwyer, D.A.; Groom, D.E.; Lin, C.-J.; Lugovsky, K.S.; Pianori, E.; Robinson, D.J.; et al. Review of Particle Physics. PTEP 2020, 2020, 083C01. [Google Scholar] [CrossRef]
  16. Straub, D.M. flavio: A Python package for flavour and precision phenomenology in the Standard Model and beyond. arXiv 2018, arXiv:1810.08132. [Google Scholar]
  17. Altmannshofer, W.; Niehoff, C.; Stangl, P.; Straub, D.M. Status of the BK*μ+μ anomaly after Moriond 2017. Eur. Phys. J. C 2017, 77, 377. [Google Scholar] [CrossRef] [Green Version]
  18. Aubert, B.; Bona, M.; Boutigny, D.; Karyotakis, Y.; Lees, J.P.; Poireau, V.; Prudent, X.; Tisserand, V.; Zghiche, A.; Tico, J.G.; et al. Search for the decay B+K+τμ±. Phys. Rev. Lett. 2007, 99, 201801. [Google Scholar] [CrossRef] [Green Version]
  19. Lees, J.; Raven, H.G.; Collaboration BABAR. A search for the decay modes B+−h+−τ+−l. Phys. Rev. D 2012, 86, 012004. [Google Scholar] [CrossRef] [Green Version]
  20. Aebischer, J.; Altmannshofer, W.; Guadagnoli, D.; Reboud, M.; Stangl, P.; Straub, D.M. B-decay discrepancies after Moriond 2019. Eur. Phys. J. C 2020, 80, 252. [Google Scholar] [CrossRef] [Green Version]
  21. Altmannshofer, W.; Gori, S.; Pospelov, M.; Yavin, I. Quark flavor transitions in LμLτ models. Phys. Rev. D 2014, 89, 095033. [Google Scholar] [CrossRef] [Green Version]
  22. Allanach, B.; Queiroz, F.S.; Strumia, A.; Sun, S. Z models for the LHCb and g−2 muon anomalies. Phys. Rev. D 2016, 93, 055045. [Google Scholar] [CrossRef]
  23. Demir, D.A.; Kane, G.L.; Wang, T.T. The Minimal U(1)’ extension of the MSSM. Phys. Rev. D 2005, 72, 015012. [Google Scholar] [CrossRef] [Green Version]
  24. Duan, G.H.; Fan, X.; Frank, M.; Han, C.; Yang, J.M. A minimal U(1) extension of MSSM in light of the B decay anomaly. Phys. Lett. 2019, B789, 54–58. [Google Scholar] [CrossRef]
  25. Allanach, B.; Gripaios, B.; Tooby-Smith, J. Anomaly cancellation with an extra gauge boson. Phys. Rev. Lett. 2020, 125, 161601. [Google Scholar] [CrossRef]
  26. Allanach, B. U(1)B3-L2 Explanation of the Neutral Current B-Anomalies. arXiv 2020, arXiv:2009.02197. [Google Scholar]
  27. Aad, G.; Abbott, B.; Abbott, D.C.; Abdinov, O.; Abud, A.A.; Abeling, K.; Abhayasinghe, D.K.; Abidi, S.H.; Abouzeid, O.; Abraham, N.; et al. Search for high-mass dilepton resonances using 139 fb-1 of pp collision data collected at s =13 TeV with the ATLAS detector. Phys. Lett. 2019, B796, 68–87. [Google Scholar] [CrossRef]
  28. Borzumati, F.; Farrar, G.R.; Polonsky, N.; Thomas, S.D. Soft Yukawa couplings in supersymmetric theories. Nucl. Phys. B 1999, 555, 53–115. [Google Scholar] [CrossRef] [Green Version]
  29. Dembinski, H.; Ongmongkolkul, P.; Deil, C.; Hurtado, D.M.; Schreiner, H.; Feickert, M.; Andrew; Burr, C.; Rost, F.; Pearce, A. Scikit-Hep/Iminuit: v2.3.0. Available online: https://zenodo.org/record/4460551#.YA_XB1kzZhF (accessed on 11 January 2021).
  30. James, F.; Roos, M. Minuit—A System for Function Minimization and Analysis of the Parameter Errors and Correlations. Comput. Phys. Commun. 1975, 10, 343–367. [Google Scholar] [CrossRef]
  31. Aoki, S.; Aoki, Y.; Bečirević, D.; Blum, T.; Colangelo, G.; Collins, S.; Morte, M.D.; Dimopoulos, P.; Dürr, S.; Fukaya, H.; et al. FLAG Review 2019: Flavour Lattice Averaging Group (FLAG). Eur. Phys. J. C 2020, 80, 113. [Google Scholar] [CrossRef] [Green Version]
  32. D’Amico, G.; Nardecchia, M.; Panci, P.; Sannino, F.; Strumia, A.; Torre, R.; Urbano, A. Flavour anomalies after the RK* measurement. JHEP 2017, 9, 010. [Google Scholar] [CrossRef]
  33. Allanach, B.; Butterworth, J.; Corbett, T. Collider constraints on Z models for neutral current B-anomalies. JHEP 2019, 8, 106. [Google Scholar] [CrossRef] [Green Version]
  34. Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Ambrogi, F.; Asilar, E.; Bergauer, T.; Brandstetter, J.; Brondolin, E.; Dragicevic, M.; Erö, J.; et al. Search for high-mass resonances in dilepton final states in proton-proton collisions at s = 13 TeV. JHEP 2018, 6, 120. [Google Scholar] [CrossRef] [Green Version]
  35. Aad, G.; Abbott, B.; Abbott, D.C.; Abud, A.A.; Abeling, K.; Abhayasinghe, D.K.; SHAbouZeid, O.S.A.; Abraham, N.L.; Abramowicz, H.; Abreu, H.; et al. Search for t t ¯ resonances in fully hadronic final states in pp collisions at s = 13 TeV with the ATLAS detector. JHEP 2020, 10, 061. [Google Scholar] [CrossRef]
  36. Martin, A.D.; Stirling, W.J.; Thorne, R.S.; Watt, G. Parton distributions for the LHC. Eur. Phys. J. 2009, C63, 189–285. [Google Scholar] [CrossRef] [Green Version]
  37. Clark, D.B.; Godat, E.; Olness, F.I. ManeParse: A Mathematica reader for Parton Distribution Functions. Comput. Phys. Commun. 2017, 216, 126–137. [Google Scholar] [CrossRef] [Green Version]
  38. Fuks, B.; Klasen, M.; Ledroit, F.; Li, Q.; Morel, J. Precision predictions for Z-production at the CERN LHC: QCD matrix elements, parton showers, and joint resummation. Nucl. Phys. 2008, B797, 322–339. [Google Scholar] [CrossRef] [Green Version]
  39. Bhattacharya, B.; Datta, A.; Guévin, J.P.; London, D.; Watanabe, R. Simultaneous Explanation of the RK and RD(*) Puzzles: A Model Analysis. JHEP 2017, 1, 015. [Google Scholar] [CrossRef]
  40. Staub, F.; Ohl, T.; Porod, W.; Speckner, C. A Tool Box for Implementing Supersymmetric Models. Comput. Phys. Commun. 2012, 183, 2165–2206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Overlapping 1-3 σ regions due to experimental constraints on the indicated observables for β 13 = 0 (left) and β 23 = 0 (right). Our benchmark points (BMPs) are marked with an asterisk.
Figure 1. Overlapping 1-3 σ regions due to experimental constraints on the indicated observables for β 13 = 0 (left) and β 23 = 0 (right). Our benchmark points (BMPs) are marked with an asterisk.
Symmetry 13 00191 g001
Figure 2. Constraints on g E and M Z due to the searches p p Z l l [27]. The left figure corresponds to β 13 = 0 , while the right one to β 23 = 0 . The benchmark points are indicated and the corresponding parameters are shown. The shaded region is excluded. Straight lines correspond to different values of M Z / g E in TeV, and the green band reflects the uncertainty of the fitted value.
Figure 2. Constraints on g E and M Z due to the searches p p Z l l [27]. The left figure corresponds to β 13 = 0 , while the right one to β 23 = 0 . The benchmark points are indicated and the corresponding parameters are shown. The shaded region is excluded. Straight lines correspond to different values of M Z / g E in TeV, and the green band reflects the uncertainty of the fitted value.
Symmetry 13 00191 g002
Table 1. The anomaly-free U ( 1 ) charges considered in the paper.
Table 1. The anomaly-free U ( 1 ) charges considered in the paper.
Field Q Field Q Field Q
Q 1 , 2 0 U 1 , 2 c 0 D 1 , 2 c 0
Q 3 + 1 U 3 c 1 D 3 c 1
L 1 , 2 1 E 1 , 2 c + 1 ν 1 , 2 c + 1
L 3 0 E 3 c + 1 ν 3 c 0
H d 1 H u 0S + 1
Table 2. Model predictions for the best fit values. BMP1 corresponds to α 23 = 0.054 , β 23 = 0.66 , M Z = 3160 GeV, g E = 0.3 , while BMP2 to α 23 = 0.19 , β 13 = 1.04 , M Z = 2300 GeV, g E = 0.135 . We do not give experimental bounds [27] on σ l l ( s ) σ ( p p Z l l ) in the table since they are different for different M Z masses.
Table 2. Model predictions for the best fit values. BMP1 corresponds to α 23 = 0.054 , β 23 = 0.66 , M Z = 3160 GeV, g E = 0.3 , while BMP2 to α 23 = 0.19 , β 13 = 1.04 , M Z = 2300 GeV, g E = 0.135 . We do not give experimental bounds [27] on σ l l ( s ) σ ( p p Z l l ) in the table since they are different for different M Z masses.
Obs R K [ 1.0 6.0 ] R K * [ 1.1 6.0 ] P 5 [ 4.0 6.0 ] Δ M B s , ps 1 BR ( B K μ ( e ) τ ) σ ll ( 13 TeV ) , fb σ ll ( 14 TeV ) , fb
Exp 0.846 ( 63 ) 0.69 ( 11 ) 0.439 ( 111 ) 17.757 ( 21 ) <4.5 ×   10 5 --
SM [16] 1.0008 ( 3 ) 0.9964 ( 6 ) 0.76 ( 8 ) 18.9 ( 1.2 ) 0--
BMP1 1.21 ( 20 ) 0.76 ( 7 ) 0.43 ( 9 ) 17.4 ( 4.3 ) 7 (6) ×   10 8 0.017 0.029
BMP2 0.65 ( 6 ) 0.87 ( 2 ) 0.34 ( 9 ) 17.5 ( 8.3 ) 4 (2) ×   10 9 0.067 0.099
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bednyakov, A.; Mukhaeva, A. Flavour Anomalies in a U(1) SUSY Extension of the SM. Symmetry 2021, 13, 191. https://doi.org/10.3390/sym13020191

AMA Style

Bednyakov A, Mukhaeva A. Flavour Anomalies in a U(1) SUSY Extension of the SM. Symmetry. 2021; 13(2):191. https://doi.org/10.3390/sym13020191

Chicago/Turabian Style

Bednyakov, Alexander, and Alfiia Mukhaeva. 2021. "Flavour Anomalies in a U(1) SUSY Extension of the SM" Symmetry 13, no. 2: 191. https://doi.org/10.3390/sym13020191

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop