Impact of non-universal Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{Z}}^{{'}}$$\end{document} in the lepton flavour violating BBs→K∗φl1-l2+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{B}}\left( {{B}}_{{{s}}} \right) {\rightarrow }{{K}}^{{*}}\left( {\varphi } \right) {{l}}_{{{1}}}^{{-}}{{l}}_{{{2}}}^{{+}}$$\end{document} decays

In recent years, lepton flavour violating (LFV) decays are one of the most trending topics to probe new physics (NP). The latest results of LHCb have motivated us to study the LFV decays through b→s\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$$\end{document} transition. The branching ratios, forward backward asymmetries and longitudinal polarization fractions of B(Bs)→K∗(φ)l1l2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B(B_{s})\rightarrow K^{*}(\varphi )l_{1}l_{2}$$\end{document} decays are studied in non-universal Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^{'}$$\end{document} model. Here, we have structured the four-fold angular distribution of the decays in terms of transversity amplitudes. The variation of the observables in whole kinematic region shows the sensitivity of NP. The observables estimated in this work are very intriguing and might lead a new track towards NP in near future.


Introduction
Inspired by the various LHCb results of lepton flavour violation on  →  transition, we study the lepton flavour violating (  ) →  * () 1  2 decays in terms of transversity amplitudes in non-universal  ′ model.Recently, some B meson decay measurements disagree with the prediction of the standard model (SM) known as B flavour anomalies.Some of them are: (i) neutral-current anomalies involved by  →  +  − transition, (ii) charged-current anomalies involved by  →  − ̅  transition, (iii) lepton flavour violation (LFV).The LFV decays are highly suppressed in the SM as the predicted values in the SM lie far below recent experimental results.For example, in ref. [1], branching ratios for   →  ±  ∓ and   →  ±  ∓ decays are calculated in the SM of order 10 −54 whereas Babar and LHCb has constrained the values at the order of 10 −5 with 90% and 95% C. L in refs.[2,3].The SM theory conserves the generation lepton number of electroweak interactions but the lepton flavour violation has been traced by the observations of neutrino oscillation.The mismatch between weak and mass eigenstates of neutrinos causes the neutrino oscillation as well as the mixing between various generation of leptons via charged current interactions of W boson [4].It is well-known that the principle for the flavour changing neutral current (FCNC) transitions of lepton sector in LFV decays could be similar to that for the quark sector.The FCNC transitions in the lepton sector which are due to the mixing in the charged current (CC) interactions with the left-handed W boson and light neutrinos are very small because they are suppressed by powers of   2   2 ⁄ .There are many other observables to examine lepton flavour universality (LFU) with FCNC.After the discovery of anomaly of the observable  5 ′ , a discrepancy in the measurement of   ( * ) = (→ ( * )  +  − ) (→ ( * )  +  − ) was found.It is recently announced that (  )  = 0.846 −0.054−0.014+0.060+0.016, 1 ≤  2 ≤ 6 GeV 2 [5] (where  2 is the momentum transfer term) which is lesser than the SM value (  )  = 1.00 ± 0.01 [6] by 2.5 deviation.The value of   * is observed at LHCb [7].Belle experiment has also constrained   and   * values which are nearer to the SM [8] prediction but with greater uncertainties.The lepton flavour universality violation test is performed for the charged-current transitions in 2012 for the first time with the observable   and   * .The colliders of the B factories: Belle [9][10][11], BaBar [12] and LHCb [13] have set their values of   * .Recently, Belle has announced the values as [14]:   = 0.307 ± 0.37 ± 0.016 and   * = 0.283 ± 0.018 ± 0.14.All these values are more than the SM values given in ref. [15] and ref. [16] by 2.3 and 3.4 deviations respectively.
Although the SM has provided large number of predictions and most of them are experimentally confirmed, it provides many unanswered questions in terms of the above spoken anomalies.So it is clear that the SM is not complete and we need the physics beyond the SM.Various theoretical models are proposed to explain the experimental tensions of LFV processes within current sensitivity.The LFV decays, such as  → 3,  → 3,  →  ′  (where ,  ′ are leptons of different flavours and  reprsents meson) and radiative decays  →  etc, all are analysed in different NP models though they are not evidently seen experimentally but their experimental upper limits exist.The LFV decays are successfully illustrated with the effect of FCNC mediated  boson [4,17], non-universal  ′ boson [18], leptoquarks [19] and other NP models and also in model independent way [20].
In this work, we study the differential branching fractions, forward backward asymmetries and longitudinal polarization fractions of LFV decays  →  *  1  2 and   →  1  2 induced by the quark level transition b→  1  2 in  ′ model where  1 and  2 are charged leptons of different flavours.We constrain the NP couplings using several experimental upper limits.
Here, we consider non-universal  ′ model [21][22][23][24][25].In this model, the NP contributes at tree level via  ′ -mediated flavour changing  →  decays where  ′ associates to the quark sector ̅  as well as to the leptonic sector  1 −  2 + .Conventionally in different grand unified theories (GUTs) the  ′ mass is considered differently as it is not unveiled till now.So the  ′ mass is restricted by various experiments and detectors constraining its upper and lower limit.The model-dependent lower bound is set at 500 GeV [26][27][28].The mass of  ′ boson is predicted by Sahoo et al. in the range 1352-1665 GeV from   −  ̅  mixing [29] and the LHC also constrained the  ′ mass, the  −  ′ mixing angle and the coupling of extra (1) gauge group [30][31][32].Lower limit of  ′ mass is set as 2.42 TeV [30] and 4.1 TeV [33] for sequential standard model (SSM) and for  6 motivated   ′ by the ATLAS collaboration respectively.The CMS collaboration has also set the lower bound of   ′ as 4.5 TeV for sequential standard model (SSM) and 3.9 TeV for superstring-inspired model [31].Bandopadhyay et al. have constrained as   ′ > 4.4 TeV using the recent Drell-Yan data of the LHC [32].Recently [34], the mass difference of   meson is studied in extended standard model and the upper limit of  ′ mass is constrained as 9 TeV.The ATLAS collaboration has set the lower limit of  ′ mass at 4.5 TeV and 5.1 TeV for   ′ and   ′ at 95% C. L. [35].In this paper, we have taken the  ′ mass in TeV range.
The paper is arranged as follows: The effective Hamiltonian for LFV decays is discussed in Sec.II.The kinematics of the decays, detail of transversity amplitudes and the description of the observables are illustrated in Sec.III.The NP is introduced in the next section, i.e. the contribution of non-universal  ′ model is incorporated in Sec.IV.We have presented the numerical analysis in Sec.V.And finally we have concluded the findings in Sec.VI.
In the similar manner, we can write in the  * () rest frame where   ,  () and |p  | are defined by the similar type expressions.
The full angular distribution of the above decay can be read as, Proceeding with these considerations, we have found the differential decay rate of the above decay as, Another interesting tool for exploring new physics in these lepton flavour violating decays is the forward-backward asymmetry which is given below, Using the above expressions, we have formulated the following form as: The longitudinal polarization fraction can be extracted as, Here, angular coefficients are expressed in terms of transversity amplitudes  ⊥ ,() ,  ∥ ,() ,  0 ,() and   ,() .The form factors are taken from ref. [42] and the detail is discussed in Appendix A.

𝑅
are the lefthanded and right-handed couplings with the leptons respectively.The coupling parameter consists of a NP weak phase term, which is related as    = |   | −   .
Considering eq. ( 1) and eq.( 2), we have incorporated the NP terms as follows: Including the NP parts through the modified Wilson coefficients, we have studied the observables defined in eqs.( 5), ( 7) and ( 8) in the BSM physics.

V. Numerical Analysis
In this work, we study the differential branching ratio, forward backward asymmetry and longitudinal polarization for  →  transition.In our previous work [57], we have studied lepton flavour violating Λ  decays and constrained the LFV leptonic couplings from the recent bound of various LFV decays of same quark transition.The LFV quark couplings are taken from refs [39,58] and are recorded in Table-1.The form factor values are taken from Appendix B [42].The values of other input parameters are tabulated in the Appendix D. ) .( 16) The scenario  9 ′ =  10 ′ is strongly abandoned because it produces the result   ≳ 1 [62] which is not in agreement with the experiments.In ref. [63]  With all these numerical data we have varied the differential branching ratio within allowed kinematic region of  2 in Fig. 1a and Fig. 1b for  →  *  and  →  *  channels respectively and Fig. 2a and Fig. 2b for   →  and   →  channels respectively.Here, all the plots are drawn for the central values of the input parameters and form factors.The blue line represents the variation for scenario 2 and orange line represents the scenario 1.We observe that the values of differential branching ratios are greater at low kinematic region.The maximum value at scenario 2 infers the appreciable contribution of NP on the decays.Integrating over the whole  2 phase space we have estimated the average values of branching ratios for these decays which are shown in Table-3.We assume 10% uncertainties on the form factors [42] (as it is one of the main sources of vagueness in the evaluation of observable) and the error on the input values (which are indicated in Appendix D) in our calculation.Another promising tool to explore NP in these decays is the lepton side forward-backward asymmetries.The sensitiveness of the observable on the LFV Λ  baryon decays has been studied in the references [40,57].We have studied the variation of forward-backward asymmetries in the  ′ model in the allowed  2 region in Fig. 3a and Fig. 3b for  →  *  and  →  *  channels respectively, whereas Fig. 4a and Fig. 4b represent the variance for   →  and   →  channels respectively.The colour description of the figures is similar as previous.The plots interpret that the zero crossing positions are very responsive to NP.Here, Fig. 3a shows the position of zero crossing of  →  *  for scenario 2 at  2 = 10.58GeV 2 and for scenario 1 at  2 = 10.08 GeV 2 whereas in Fig. 3b the scenario 2 shows the zero crossing at  2 = 6.245GeV 2 for  →  *  but it is absent for scenario 1.The zero crossing is also absent for   →  in Fig. 4a and the observable decreases gradually, reaches at  minimum value and then again increases.At low  2 region the variation of both the scenarios are not considerable but at high  2 regime the 2 nd scenario attains the minimum value.And Fig. 4b shows the variation of the observable for   →  channel.Here, we observe the zero crossing for the 2 nd scenario only at  2 = 1.495GeV 2 and the 1 st scenario has increasing nature throughout the allowed  2 region.The values of the forward-backward asymmetries are calculated and represented in Table-4.

VI. Conclusion
Undoubtedly the SM achieved an enormous success in the theoretical community of high energy physics.But lepton flavour violation is still suppressed in the SM which infers its imperfection.So it can be said that there must be some physics beyond the SM.Recently the LHCb has set a few bounds on LFV  →  1 −  2 + decays [65].These results of LHCb accelerate the BSM physics towards a new horizon incorporating several advanced and remarkable analysis.Previously the LFV decays are studied in different NP models.In this paper, we have studied  →  *  1 −  2 + and   →  1 −  2 + decays in non-universal  ′ model.We have calculated the branching ratio values, forward backward asymmetries and longitudinal polarization asymmetries for these decays.In our last work, we have constrained the LFV leptonic couplings [57] which are applied here.Our plots show that the 2 nd scenario provides higher values of the observables which offer a clear cut signature of NP in these LFV decays.We notice the existence of zero crossings in the plots of   but  →  *  and   →  decays do not show the zero crossing positions.Actually zero crossing position is the point where the uncertainties due to form factor terms are cancelled out and the sensitivity due to NP can be identified.The zero position changes for two different scenarios of NP in Fig. 3a, Fig. 3b and Fig. 4a and this responsiveness proves the effect of NP.The   value attains greater negative value for   →  decay whereas (  ) →  * () have positive value throughout the whole kinematic region.These changes in nature of the observable are due to leptonic masses.
Another interesting observable we have included in this work is longitudinal polarization fraction.In the advance 2000's, this observable has gained attention for the first time through the well-known Polarization Puzzle.We observe that the observable serves almost similar variation for all the decays and also the higher value of NP couplings provide higher value of the observable.This character evolves the probing nature of the observable towards NP.
Previously the lepton violating B decays are studied with the effect of scalar leptoquarks [66], in SUSY [67], with the contribution of pseudo-scalar co-efficients (coupling to Higgs) [37], in generic  ′ models [38].These previous studies are mainly based upon the branching ratio except the reference [37].In scalar leptoquark model and SUSY the upper limit of the branching ratio for  →  *  was predicted as (3.4 × 10 With the quick expansion of LHCb energy and upcoming Belle II experiments the LFV decays can be investigated with high precision and accuracy.The study of these decays are actually smoking signal of NP which will be proved in near future with the developments of colliders.
The results obtained in Table-3, 4 and 5 will be very useful to the experimentalists as well as the theorists.

Appendix C
The new  ′ couplings help to explain the FCNC transitions at tree level.The FCNCs depend on the individual unitary transformations for the left and right chiral quarks (u and d), leptons and neutrinos diagonalizing their particular mass matrices [23].However, only the combination of left chiral quarks (u and d) that are occurred in the CKM matrix is known experimentally (with weak constraints on the leptonic analog from neutrino oscillations).
Mathematically, we have the current correspond to the additional (1) gauge group in eq.(11) where   ,  are the right handed and left handed chiral couplings.These chirals for up quark are flavour diagonal and family universal as Fig. 1aFig.1b

1 )
[54][55][56]:  ,  =  ,  ,  ,  =  ,   and    =    where  represents the 3 × 3 identity matrix and  ,  ,  ,  and    represent the chiral charges.The chirals for down quark are represented as: .Here, the family non-universality can be signalized by the term .Diagonalization of Yukawa matrices indicates the right-handed couplings as flavour-diagonal whereas the non-diagonal left handed part unveil the nonuniversal nature of the flavour couplings.So in our work, we have considered the left-handed  −  −  ′ couplings and neglected the couplings of the  ′ with right-handed quarks.Appendix D:

It is to be noted that the LFV decays always occur at tree level in this 𝑍 ′ model; therefore the NP should contribute in the fashion where 𝑡𝑡 ̅ loops get cancelled. Moreover, there is an electromagnetic operator 𝑂 7 in the SM for 𝑏 → 𝑠𝑙𝑙 transition. Non-universal 𝑍 ′ model is
basically sensitive to the semileptonic current operators  9 and  10 involving NP contributions [44][45][46]n the non-universal  ′ model, the newly announced massive gauge boson  ′ couples to quarks and the lepton pair.Here, an additional (1) ′ gauge group is introduced with the gauge group of the SM[23,43].The FCNC transitions are successfully illustrated at tree level by the characteristic of the  ′ couplings with fermions.Acording to refs.[44][45][46] ′ boson associates with the third generation of quark differently from the other two generations and it is similar for the lepton families also.We have considered different family non-universal  ′ couplings for different families in the model.I this work, we have studied with the diagonal and nonuniversal  ′ couplings.Here, this sum is done over all fermions and  , and   ,  are the chiral couplings of the new gauge boson. Atually the FCNCs originate in both left-handed and right-handed sectors at the ̅      +    ̅      ) (       ̅      , ,  and  ′ are the gauge couplings of  and  ′ bosons (here,  = .IV.  ′ contribution in  →  * ()  −   + decayswhere  ′ is the new gauge coupling linked with the extra (1) ′ group and the effective Hamiltonian has the form as, for  ′ of TeV-scale .The LEP experiments have also recommended the  ′ existence with the identical couplings as the SM  boson.If |   |~| (14) then    will be (10 −3 ) .With the above summarization the effective Hamiltonian for  →  * (̅ ) − ( ̅ 1  2 ) + ] + ℎ. ,(14)where    is the left-handed coupling of  ′ boson with quarks,   1  2and   1  2

Table - 1: Numerical values of coupling parameters
Here, we have considered the maximised values of the coupling parameters to intensify the NP in these decays.The leptonic couplings are composed in the following way:  1  2 = (  1  2  +   1  2  ) , and   1  2 = (  1  2  −   1  2

Table - 2: Maximum values of NP couplings
also this consideration  9 ′ =  10 ′ does not match the experimental values of branching ratios of   →  +  − and  →  +  − at high  2 region; whereas there are plethora of best results of the anomalies with the scenario  9 ′ = − 10 ′ [40, 59-64].Not only the model dependent analysis but also the model independent analysis support the fruitfulness of the well-considered scenario.Therefore, we have adapted the scenario  9 ′ = − 10 ′ in our investigation.Comprising all the above discussions, we have set two scenarios with maximum values of the NP couplings as Table-2 below.

Table - 3: Predicted values of branching ratios for LFV
→  * and   →  decays

decays in 1 st and 2 nd scenarios
Longitudinal polarization asymmetry, one of the most efficient tools to test NP, for  →  * (→  −  + ) 1  2 and   → (→  ̅ ) 1  2 decays are studied in non-universal  ′ model for the first time in this work.The variation of the observable is presented in the allowed kinematic