Model independent analysis of B∗→Pℓν¯ℓ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ B^* \rightarrow P \ell {\bar{\nu }}_\ell $$\end{document} decay processes

Very compelling deviations in the recently observed lepton nonuniversality observables (RD(∗),RK(∗),RJ/ψ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ (R_{D^{(*)}}, R_{K^{(*)}}, R_{J/\psi } )$$\end{document} of semileptonic B meson decays from their Standard Model predictions hint towards the presence of some kind of new physics beyond it. In this regard, we investigate the effect of new physics in the semileptonic B¯d(s)∗→Pℓν¯ℓ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{B}}_{d(s)}^* \rightarrow P \ell {\bar{\nu }}_\ell $$\end{document} decay processes, where P=D,π(Ds,K\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=D,\pi (D_s,K$$\end{document}), in a model independent way. We consider the presence of additional vector and scalar type interactions and constrain the corresponding new couplings by fitting Br(Bu+→τ+ντ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{Br(B_{u}^+ \rightarrow \tau ^+ \nu _\tau )}$$\end{document}, Br(B→πτν¯τ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{\mathrm Br}(B \rightarrow \pi \tau {{\bar{\nu }}}_\tau )}$$\end{document}, Br(Bc+→τ+ντ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{Br(B_{c}^+ \rightarrow \tau ^+ \nu _\tau )}$$\end{document}, Rπl\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_\pi ^l$$\end{document}, RD(∗)\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} and RJ/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{J/\psi }$$\end{document} data. Using the constrained new parameters, we estimate the branching ratios, forward–backward asymmetry, lepton-spin asymmetry and lepton non-universality observables of B¯d,s∗→Pτν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{B}}_{d,s}^{*} \rightarrow P \tau {{\bar{\nu }}}_\tau $$\end{document} processes. We find that the branching ratios of these decay modes are sizeable and deviate significantly (for vector-type couplings) from their corresponding standard model values, which are expected to be within the reach of Run III of Large Hadron Collider experiment.

Being motivated by these observed anomalies in various B meson decays, in this work we would like to investigate the impact of new physics on the differential decay rate and various other observables like forward-backward asymmetry, lepton-spin asymmetry and lepton nonuniversality (LNU) observable of weakly decaying vector B * d, (s) meson to a pseudoscalar P(= D(D s ), π(K )) meson mediated through the quark level transitions b → (c, u) ν . Although such hadrons decay primarily through the electromagnetic process B * d,s → B d,s γ , and their weak decay channels are expected to be quite suppressed, the situation has improved considerably with the advent of the high luminosity Belle II experiment. For instance, as discussed in Ref. [21], using the production cross section of ϒ(5S) in e − e + collision as σ (e + e − → ϒ(5S)) = 0.301nb and Br(ϒ(5S) → B * B * ) = (38.1 ± 3.4)% [20], about 4 × 10 9 B * meson pairs (B * u,d +B * u,d ) are expected to be produced per year. This in turn implies that the rare B * decay modes with branching fraction > O(10 −9 ) are likely to be observed at Belle II. Hence, Belle II experiment would be quite instrumental in search for the rare decay modes of the excited B mesons. In addition the LHC experiment will also play a pivotal role in the search for B * decay channels, as the production cross section of ϒ(5S) is much larger in pp collision compared to e + e − collision. On the other hand, the study of B * meson decays has also received considerable attention in recent times. In the literature [22][23][24][25], the leptonic decay modes of B * s,d mesons are investigated in SM and in the context of various new physics models. The analysis of semileptonic weak decays B * → P ν both in the SM and in the presence of NP are discussed in the Refs. [21,26,27].
The layout of the paper is as follows. In Sect. 2, we illustrate the theoretical framework required to analyse the decay processes B * → P ν in the effective theory formalism. The expressions for the differential decay rate and other observables like forward-backward asymmetry, lepton nonuniversality (R * P ) and the lepton-spin asymmetry are presented in this section. The constraints on the new couplings using χ 2 fit from R D ( * ) , R J/ψ , R l π , Br(B u,c → τ ν), Br(B → πτν) observables are obtained in Sect. 3. Our results are discussed in Sect. 4 followed by the summary of our work in Sect. 5.

Theoretical framework
The most general effective Lagrangian for B * → P ν processes mediated by b → q −ν (q = u, c), in the effective field theory approach can be expressed as [28], where P is any pseudoscalar meson, G F is the Fermi constant, V qb is the CKM matrix element, V L ,R , S L ,R , T L are the new vector, scalar, and tensor type new physics couplings, which are zero in the standard model. All these new physics couplings are considered to be complex. Furthermore, we consider the neutrinos as left handed. We assume the NP effect is mainly through the third generation leptons and do not consider the effect of tensor operators in our analysis for simplicity. Here (q, ) L ,R = P L ,R (q, ), where P L ,R = (1 ∓ γ 5 )/2 are the chiral projection operators.
We consider the kinematics of the decay process B * → P ν using helicity amplitudes. In this formalism, the decay process B * → P ν is considered to proceed throughB * → PW * − , where the off-shell W * − decays to −ν . One can write the amplitude from Eq. (5) as where C k (μ) represents the Wilson coefficient with values C k (μ) = 1 for SM, V L ,R , S L ,R for NP beyond SM, k denotes the product of gamma matrices, which gives rise to different Lorentz structure of hadronic and leptonic currents of Eq. (5) i.e., k = γ μ (1±γ 5 ), and (1±γ 5 ). Hence, the square of the matrix element can be expressed as the product of leptonic (L μν ) and hadronic (H μν ) tensors (related to the corresponding helicity amplitudes) where the superscripts i, j represent the combination of four operators (V ∓ A), (S ∓ P) in the effective Lagrangian (5), C i j (μ) denotes the product of Wilson coefficients C i and C j . We omit these superscripts in the following discussion for convenience. It should be noted that, the polarization vector of the off-shell particle W * (¯ μ (m)), satisfies the following orthonormality and completeness relations: * μ (m)¯ μ (m ) = g mm , where g mm = diag(+, −, −, −) and m, m = ±, 0, t represent the transverse, longitudinal and time-like polarization components. Now inserting the completeness relation from Eq. (8) into (7), the product of L μν and H μν can be expressed as where L(m, n) = L μν¯ μ (m)¯ * ν (n) and H (m, n) = H μν¯ * μ (m)¯ ν (n) are the Lorentz invariant parameters, and hence their values are independent of any specific reference frame. So for calculational convenience, we will evaluate H (m, n) in the B * rest frame and L(m, n) in −ν center of mass frame as discussed in [21,26].

Hadronic helicity amplitudes
In the rest frame of B * meson, we consider the pseudoscalar meson P to be moving along the positive z-direction. The polarization vector of the virtual W * boson are chosen to bē 2 , is the momentum transferred square and λ(a, b, c) = a 2 + b 2 + c 2 − 2(ab + bc + ca). The polarization vector of the on-shell B * meson ε μ (m = 0, ±), takes the form In order to calculate the hadronic helicity amplitudes, we use the following matrix elements of B * → P transition where V (q 2 ), A 0,1,2 (q 2 ) are the various form factors. The matrix elements for the scalar and pseudoscalar currents can be obtained by using the equation of motion as where the m b,q represent the current quark masses evaluated at the b-quark mass scale. The helicity amplitudes are defined as where for convenience, we use the notations λ B * = 0, ± and λ W * = 0, ±, t to represent the helicity states of the B * and W * boson. Thus, with Eqs. (12), (14) and (15), one obtains the following non-vanishing helicity amplitudes

Leptonic helicty amplitudes
The leptonic helicity amplitudes are defined as where In the center of mass frame of −ν , the four momenta of andν pair are expressed as where E = (q 2 + m 2 )/2 q 2 , | p | = (q 2 − m 2 )/2 q 2 and θ is the angle between the three momenta of of P and . The polarization vector of the virtual W * boson in this frame is Thus, with Eqs. (17) and (19), one obtains the following nonvanishing contributions

Decay distribution and other observables
The double differential decay rate of B * → P ν decay process can be expressed as Now, with Eqs. (16) and (20), one can obtain L μν H μν in terms of Wigner d J -functions as [26] L μν H μν where J and J take the values 0 and 1 and the various helicity components run over their allowed values. Thus, one can obtain the differential decay rate to particular leptonic helicity state (λ = ± 1 2 ) as From Eqs. (24) and (25), one can obtain the differential decay rate as where the values of the helicity amplitudes are given in Eq. (16). Apart from the differential decay rate, the other NP sensitive observables, considered here are • Lepton nonuniversality observable: where l denotes the light leptons l = e, μ. • Forward-backward asymmetry: which can be expressed in terms of the helicity amplitudes as where the parameters X and Y are given as • Lepton-spin asymmetry:

Form factors and their q 2 dependence
The main inputs required for the numerical analysis are the values of the form factors. As the first principle lattice calculation results of the form factors for B * d,s → D, D s (π, K ) transitions are not yet available, we use their values evaluated in the BSW model [29,30]. Their values at zero-momentum transfer are listed below The q 2 dependence of the form factors can be written as, where m B q (0 ± ) and m B q (1 ± ) are the pole masses whose values are presented in Table 1. In our analysis, we consider 10% uncertainty in the values of hadronic form factors at q 2 = 0.

Constraints on new couplings
In this analysis the new couplings are considered to be complex. Considering the contribution of only one coefficient at a time with all others set to zero, we perform the chi-square fitting for the individual complex couplings. The χ 2 is defined as where O th i represent the theoretical predictions of the observables, O  2 contain the 1σ errors from theory and experiment. We constrain the real and imaginary parts of new coefficients related to b → cτν τ quark level transitions from the χ 2 fit of R D ( * ) , R J/ψ and Br(B + c → τ + ν τ ) observables and the couplings associated with b → uτν τ processes are constrained from the fit of R l π , Br(B + u → τ + ν) and Br(B 0 → π + τ −ν ) data. The updated values of all the observables used for fitting are taken from [20] and are listed in Table 2. The upper limit on the branching ratio of B + c → τ + ν τ decay mode with the present world average of the B c lifetime is [31] Br(B + c → τ + ν τ ) 30%.
We use the theoretical expressions of these observables and their SM predictions from [32] and have listed them in Table 2. Br(B 0 → π + τ ν) < 2.5 × 10 −4 (9.40 ± 0.75) × 10 −5  Fig. 2. It should be noted that, the best-fit values are degenerate in the presence of V L coupling (V L , S L and S R couplings) for b → c (b → u) processes and for each of these couplings, we have considered only benchmark values. The best-fit values and the corresponding 1σ ranges, which are obtained from the joint confidence regions of the real and imaginary planes of these new couplings, are presented in Table 3. The χ 2 /d.o.f, as well as the pull χ 2 SM − χ 2 best−fit , for all the coefficients are also listed in this Table. One can notice that, the Wilson coefficient corresponding to b → c scalar operators have χ 2 /d.o.f > 1, which implies that the fit is not robust. However, the pull values of V L ,R coefficients of b → c implicit that the measured data are consistent with our model in the presence of either V L or V R and can be a viable candidate for explaining the b → cτν τ anomalies.

Effect of new coefficients on B * d,s → ( D, D s , π, K )τν τ decay modes
After collecting all the theoretical expressions of required observables and getting knowledge on the allowed ranges of new parameters, we now proceed towards numerical analysis. The particles masses and the values of the CKM elements and the Fermi constant G F are taken from PDG [20]. (36) Fig. 2 Constraints on individual new complex coefficients associated with b → uτν τ processes from the χ 2 fit of R l π , Br(B + u → τ + ν τ ) and upper limit on Br(B 0 → π + τ −ν τ )  )), so the branching fractions of B * + → P ν processes are roughly one-third of B * d → P ν . Hence, those results are not presented in this work. Furthermore, we assume that the new physics will couple only to third generation leptons, so the B * d,s → Pμν μ processes will not be affected by the presence of new physics operators, and their standard model branching fractions are listed in Table 4, which are expected to be within the reach of LHC experiment. TheB * d,s → (D, D s )τ −ν τ processes proceed through b → c quark level transitions, so we use the constrained values of the new couplings obtained for b → cτν τ in order to calculate the associated observables of these processes. Similarly we use the allowed parameter space obtained for b → uτν τ process to compute the observables associated with B * d,s → (π, K )τ −ν τ decay process as they are mediated by b → u quark level transitions. In the following subsections, we discuss the effect of the presence of one Wilson coefficient at a time on various observables of B * d,s → (D, D s , π, K )τ ν τ decay modes.

Effect of V L only
Here we consider the case, where the additional contribution to the SM Lagrangian arising only from V L coefficient and all other new coefficients are set to zero i.e., (S L = S R = V R = 0). Using the best-fit values and 1σ allowed parameter space of V L , obtained from the χ 2 fit of R D ( * ) , R J/ψ , , we then calculate the differential decay rate, LNU observable, lepton spin asymmetry and forward-backward asymmetry of B * 0 → D + τ ν and B * s → D + s τ ν (B * 0 → π + τ ν and B * s → K + τ ν) decay processes. In the left panel of Fig. 3, we show the q 2 variation of decay rate (top) and R * D observable (bottom) of B * 0 → D + τ ν process and the corresponding plots for B * 0 → π + τ ν channel are presented in the right panel of this figure. Here the blue dashed lines correspond to the SM prediction and the cyan bands represent the 1σ uncertainty, arising due to the errors in CKM matrix elements, hadronic form factors and the lifetime of B * meson. The solid black lines are obtained by using the best-fit values of the left handed vectorial new V L coupling and the orange bands represent the 1σ allowed ranges, which includes the SM uncertainties as well as the uncertainties due to the new couplings. From the plots, one can notice significant deviation in the branching ratios and LNU observables from their corresponding SM predictions due to presence of additional V L coefficient. To quantify these deviations, we define the pull metric at the observable level as  1.239) for B * → πτ ν process. The Pull value for R * D and R * π are found to be large as the SM uncertainties cancel out in these observ-ables, thus providing significantly large pull value. The plots for B * s → D + s τ ν (B * s → K + τ ν) process follow the same form as B * 0 → D + τ ν (B * 0 → π + τ ν), and hence, are not included in this article. The numerical values of these observables are presented in Table 5. Furthermore, no deviation has been observed in the forward-backward asymmetry and lepton-spin asymmetry observables from their SM results, so we don't provide the corresponding plots. The values of q 2 at which the forward-backward asymmetry vanishes are provided in Table 7.

Effect of V R only
In this scenario, we explore the effect of only V R coefficient on the decay rate and angular observables of B * → (D + , π + )τ ν τ processes. Using the best-fit values and corresponding 1σ allowed ranges of V R coefficients associated with b → (c, u)τν τ transitions, we present the plots for the decay rate (left-top panel), R * D (left-middle) and forward-backward asymmetry (left-bottom panel) of B * → D + τ ν decay modes in Fig. 4. The corresponding plots for B * → π + τ ν process are depicted in the right panel of Fig. 4. Here the solid black lines are obtained by using the best-fit values of new V R couplings and the gray bands by including 1σ uncertainties of all input values. Reasonable deviation in all the observables (except the leptonspin asymmetry) from their SM results are found due to the presence of additional V R coefficient, with Pull values Pull Br/R * D /A FB = 0.429/3.21/3.391 for B * → D + τ ν process and Pull Br/R * π /A FB = 0.368/1.203/1.323 for B * → πτ ν. In Table 5, we present the numerical values of decay rates and all these parameters. Due to the additional contribution from V R coefficient, we notice deviation in the zero crossing of the forward-backward asymmetry towards high q 2 and the q 2 values of the zero crossing point are given in Table 7.

Effect of S L only
In this subsection, we consider the contribution of S L new coefficient by assuming that all other new Wilson coefficients have vanishing values. As seen from Figs. 1 and 2, the S L parameters are severely constrained by the current data. Within the allowed parameter space for S L coefficient presented in Table 3, we show the q 2 variation of lepton-spin asymmetry (top) and forward-backward asymmetry (bottom) of B * → D + τν (B * → π + τν) process on the left panel (right panel) of Fig. 5. Here the plots obtained from the best-fit values (1σ range) of S L coupling are represented by dashed black lines (red bands). The numerical values of these observables are given in Table 6. With the additional S L contribution, the deviation in the branching ratios and LNU observables from their SM predictions are found to be Fig. 3 The q 2 variation of differential decay rates and LNU observ-  minimal. Though the lepton spin asymmetry and forwardbackward asymmetry observables of B * → D + τν channel provide slight deviation from their SM results, the deviation is negligible in the B * → π + τν modes. The zero crossing point of the forward-backward asymmetry of B * → D + τν process shifted sightly towards the low q 2 region. The A P FB vanishing values of q 2 predicted from the best-fit values and 1σ range of new S L coefficient are presented in Table 7.

Effect of S R only
Here we investigate the observables of B * → (D + , π + )τν decay modes by considering the presence of only additional S R coefficient. Using the available experimental data on b → (u, c)τν transitions, we fit the corresponding S R coefficients, which is already discussed in section II. In the left panel of Fig. 6, we present the q 2 variation of decay rate (top), R * D (second from top), lepton spin asymmetry (third from top) and forward-backward asymmetry (bottom) of B * → D + τν and the corresponding plots for B * → π + τν are shown in the right panel. Here the black dashed lines (magenta bands) are obtained from the best-fit values (1σ range) of S R coupling and other input parameters. In this case also, the deviation in the lepton spin asymmetry and forwardbackward asymmetry observables are comparatively large, whereas the deviations in the branching ratios and LNU observables are nominal. The numerical values are presented in Table 6. From Fig. 6, one can notice that the zero crossing point of the forward-backward asymmetry deviates significantly towards left (low q 2 region) and the corresponding q 2 values of the crossings are shown in Table 7.

Summary and conclusion
The rare decay modes of B mesons have been extensively studied both theoretically and experimentally in order to crit-  ically test the standard model prediction and to look for new physics beyond it. In this regard, the rare decay channels of the corresponding vector mesons i.e., the B * decay modes are essential as they can provide complementary ways to go beyond the standard model. However, the weak decay channels of B * vector mesons are not much explored experimen-tally as they decay dominantly through electromagnetic process B * → Bγ . Recently, with the advent of high luminosity LHCb experiment the sensitivity for the branching fractions of various rare decay modes is expected to reach the level ∼ O(10 −9 ). Thus, the LHCb would be an ideal platform to explore the rare decay modes of B * mesons.  ered the generalized effective Lagrangian in the presence of vector and scalar type new physics operators. Considering only one new coefficient to be present at a time, and assuming the new couplings as complex, we constrained the new parameters associated with b → cτν τ processes by performing χ 2 fit from R D ( * ) , R J/ψ parameters and the upper limit on B + c → τ + ν τ branching fraction. The new couplings of b → uτν τ processes are constrained by using experimental data on the branching ratios of B u → τ ν τ and B → πτ ν τ and R l π parameter. Using the best-fit values and the corresponding 1σ ranges of new individual complex Wilson coefficients, we computed the branching ratios, forwardbackward asymmetry, lepton spin asymmetry and lepton nonuniversality observables of B * d,(s) → D + (D + s )τ −ν τ and B * d,(s) → π + (K + )τν τ decay processes. We have also shown the values of q 2 at which the forward-backward asymmetry vanishes. The branching fractions and LNU observables of these decay modes in the presence of additional V L coefficient have significant deviations from their corresponding standard model predictions, whereas no deviations have been found in the lepton spin asymmetry and forward-backward asymmetry observables. Due to the additional contributions from V R new coefficient, profound deviations have observed in the decay rates, lepton nonuniversality observable and the forward-backward asymmetry of bothB * → (D, π)τν τ processes. Due to the presence of V R coupling, the zero crossing of forward-backward asymmetry has shifted towards high q 2 region for all decay modes. In the presence of S L coefficient, none of the observables are affected and there is practically no deviation from SM results. Only the leptonspin asymmetry and forward-backward asymmetry observables ofB * → Dτν τ show slight deviation due to additional S L coupling. On the other hand, in the presence of S R coupling, the lepton spin asymmetry and the forward-backward asymmetry show reasonable deviations from their SM predictions and the decay rate lepton nonuniversality observables remain unchanged. The zero crossing of forwardbackward asymmetry of all decay modes in the presence of S R coefficient is found to be shifted towards low q 2 region. To conclude, we noticed significant deviations in some of the observables from their standard model predictions in pres-ence of new couplings. The observation of these decay modes of vector B * mesons at LHC experiment will definitely shed light on the nature of new physics.