Study of the Λb→N∗ℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _{b} \rightarrow N^*\ell ^+ \ell ^- $$\end{document} decay in light cone sum rules

The form factors of the Λb→N∗ℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _{b} \rightarrow N^*\ell ^+ \ell ^-$$\end{document} decay are calculated in the framework of the light cone QCD sum rules. In the calculations the contribution of the negative parity Λb∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _b^*$$\end{document} baryon is eliminated by constructing the sum rules for different Lorentz structures. Furthermore the branching ratio of the semileptonic Λb→N∗ℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _b \rightarrow N^*\ell ^+ \ell ^-$$\end{document} decay is calculated. The numerical study for the branching ratio of the Λb→N∗ℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _{b} \rightarrow N^*\ell ^+ \ell ^-$$\end{document} decay indicates that it is quite large and could be measurable at future planned experiments to be conducted at LHCb.


Introduction
Lately, exciting experimental results have been obtained in study of the rare decays of the heavy b baryon induced by the flavor changing neutral currents. The rare b → + − decays induced by the b → s transition were observed by the CDF [1] and LHCb collaborations [2]. Later the detailed analyses of the differential branching ratio and various symmetries of the b → + − decay have been performed performed at LHCb [3]. The LHCb collaboration firstly observed the rare b → pπ − μ + μ − decay induced by the b → d transition [4]. This observation motivated the theoretical study of the b → N + − decay, induced also by the b → d transition. This decay was studied within the framework of the light cone QCD sum rules method (LCSR) in [5]. The light cone QCD sum rules method (LCSR) [6,7] is hybrid of the traditional SVZ sum rules [8] and the methods used in hard exclusive processes. The other interesting decays induced by the b → d transition are the b → nucleon resonance decays. The analysis of these decays can provide complementary information about the properties of the nucleon resonances in principle which could experimentally be studa e-mail: taliev@metu.edu.tr b e-mail: tbarakat@ksu.edu.sa c e-mail: savci@metu.edu.tr ied at LHCb. It should be noted here that the comprehensive study of the nucleon resonance constitutes one of the main research directions of the research program that is planned for future study at Jefferson Laboratory [9]. The properties of the nucleon resonance N * in the b(c) → N * ν decay is investigated in framework of the LCSR in [10]. The present work is devoted to the study of the rare b → N * + − decay in the framework of the LCSR method. The paper is organized as follows: In Sect. 2 the LCSR for the relevant form factors appearing in the b ( * b ) → N * transitions are obtained. In Sect. 3 present the numerical analysis of the sum rules for the form factors. Using then the obtained results for the form factors we estimate the decay widths of the b ( * b ) → N * + − decays. This section ends with a conclusion.

Form factors of the b ( * b ) → N * + − decay in LCSR
In the present section we derive the LCSR for the transition form factors of the b ( * b ) → N * + − decay. Before giving the details of the calculations few words about the notation should be mentioned. In all further discussions the negative parity states of the b and N baryons are denoted as * b and N * , respectively. The b ( * b ) → N * + − decay at the quark level is described by the b → d transition. At the hadronic level b ( * b ) → N * + − decay is obtained by sandwiching the transition current between the b ( * b ) and N * states. The corresponding form factors of the vector, axial vector and tensor currents are defined as, The form factors responsible for * b → N * transition can be obtained from Eqs. (1)- (3) with the help of the following replacements: 5 . In order to derive the LCSR for these form factors we introduce the correlation function j μ ( p, q) = i d 4 where η b is the interpolating current of the b ( * b ) baryon, J j μ (x) is the heavy-light transition current which is set to, In the calculations we use the following most general form of the interpolating current for the b ( * b ) baryon, where a, b and c are the color indices, C is the charge conjugation operator, and β is an arbitrary parameter with β = −1 corresponding to the Ioffe current. The usual procedure in deriving the LCSR is to calculate the correlation function given in Eq. (4) in two different domains. On one side, insert a complete set of states with the quantum numbers of b between the two currents, and isolate the ground state contribution. On the other side use the operator product expansion (OPE) around the light cone where ( p = q) 2 , q 2 0. These two representations of the correlation function are then matched using the dispersion relations and quark-gluon duality ansatz. Finally, applying the Borel transformation in order to kill the the possible subtraction terms which could appear in the dispersion relations, and to suppress the contributions from higher states. It should be remarked here that the interpolating current η b has nonzero overlap not only with the J P = 1 2 + state but also with the J P = 1 2 − state. It is shown in [11] that the mass difference between between the J P = 1 2 states is about (200-300) MeV. For this reason the contribution of the negative parity b baryon should properly be taken into account.
After having mentioned these cautionary remarks we proceed to calculate the physical part of the correlation function. Saturating Eq. (4) with the ground and first excited b baryon we get, and summation is performed over the ground and first excited states of the b baryon. The decay constants of the positive and negative parity b baryons are determined as, Using the equation of motion and Eqs. (1)-(3) and (7), for the phenomenological part of the correlation function we get, Now we turn our attention to the calculation of the correlation function from the QCD side. At deep Eucledian domain ( p − q) 2 , q 2 0 the product of the two currents can be expanded around the light-cone x 2 0. After contracting the heavy quark fields which give the heavy quark propagator, the matrix element of the three quarks between the vacuum and the N * state is revealed. Decomposition of this matrix element in terms of the distribution amplitudes (DAs) with increasing twist is given in [12] (see Appendix A). After contracting the heavy b-quark fields, the correlation takes the form, where K i are the modified Bessel functions of the second kind, and G αβ is the gluon field strength tensor. Using the expression of the heavy quark propagator and definition of the matrix in terms of the DAs of the N * baryon, we can calculate the correlation function from the QCD side. Note that, using the equation of motion ( p − m N * )u N * ( p) = 0, the correlation function can be decomposed into six independent functions as follows, After performing Borel transformation, these invariant functions in the correlation function can in general be written as, We can write (13) as a dispersion integral in ( p − q) 2 as follows, Making the replacement 2 , the denominator of Eq. (13) takes the form, Using the quark-hadron ansatz the contribution of the hadronic states can be represented as, where s 0 is the continuum contribution. Finally Borel transformation can be performed on the hadronic and physical sides with the help of the replacement In implementing the Borel transformation and the continuum subtraction we use the following relations, where x 0 is the solution of the equation Equating the coefficients of the structures p μ γ 5 , p μ / qγ 5 , γ μ γ 5 , γ μ / qγ 5 , q μ γ 5 , and q μ / qγ 5 we get the following sum rules for the invariant functions of the transition current (bγ μ d), The results for the form factors induced by thebγ μ γ 5 d current can be obtained from Eq. (18) by making the following replacements: The sum rules of the (biσ μν q ν d) transition current can be obtained in similar manner, which are given below, where I I 1 , I I 2 , I I 3 , and I I 4 are the invariant functions for the structures p μ / qγ 5 , q μ / qγ 5 , q μ γ 5 , and γ μ γ 5 , respectively. The sum rules for the (biσ μν q ν γ 5 d) transition current can be obtained from Eq. (19) by making the replacements The explicit form of these invariant functions for the aforementioned structures are presented in Appendix-B.
Solving these equations we can eliminate the * pole from the sum rules. As the result we obtain the desired sum rules responsible for the b → N * transition. In the next section we present our numerical results on these form factors.
In full QCD for describing b → N * transition we have 10 independent form factors. The number of independent form factors for b → N * transitions reduces considerably in the heavy quark limit. Using the heavy quark symmetry and the soft collinear effective theory in [13][14][15] it is obtained that the matrix element for b → N * transition can be in terms of only one universal function ξ( pv).
to the leading order. Using this relation and the definitions of the matrix elements in full QCD presented in Eqs.
(1)-(3), we immediately get the following relations among the form factors: At the end of this section, we would like to make the following remark. For a complete description of the heavy-to-light baryon transition form factors, it is necessary to take into account non-factorizable contributions. The non-factorizable contribution is described by following non-local matrix element: The h μ generated by the including to the matrix element of the four-quark operators with B = 1 effective weak Hamiltonian together with electromagnetic current for relevant quarks. The matrix element h μ does not factorize into multiplication of leptonic current and form factors. The contributions can be described by "hard spectator scattering", "weak annihilation", "hard vertex corrections", "soft gluon radiations" (see Fig. 3. These figures are taken from  Table 2 The same as Table 1, but for LQSR-2 [16]). For calculations of these contributions some nonperturbative approaches are needed (More details about the non-factorizable effects, see [15]).
In the present work, we calculate only factorizable contributions and we plan to present a detailed analysis of the non-factorizable contributions elsewhere in future.

Numerical analysis
In this section we present our numerical results on the form factors that describe the b → N * transition. First let us specify the input parameters which are needed in performing the numerical calculations. The masses of the b and * b baryons which we use in our calculations are b = 5.62 GeV and * b = 5.85 GeV, and the mass of the nucleon is m N * = 1.52 GeV [17]. The residues λ b and λ * b of the relevant baryons are taken from [5] having the values λ b = (6.5 ± 1.5) × 10 −2 GeV 3 and λ * b = (7.5 ± 2.0) × 10 −2 GeV 3 . The mass of the b quark is assigned to its M S given as m b = (4.16 ± 0.03) GeV [17].  [12] (also in [18][19][20][21][22][23]), and for completeness they are presented in Appendix A.
The sum rules for the transition form factors contain three auxiliary parameters: The Borel mass parameter M 2 , the continuum threshold s 0 and the arbitrary parameter β. The Borel mass parameter and the continuum threshold s 0 are determined from the criteria that the sum rule dictates, i.e., the suppression of the contributions coming from the continuum states and the higher twist contributions should be satisfied. Our analysis shows that the working regions of M 2 and s 0 lie in the region M 2 = (10 ± 5) GeV 2 , s 0 = (40 ± 1) GeV 2 , when aforementioned conditions are fulfilled, and hence sum rules predictions are reliable. The final step of our analysis is is the determination of the working region of the parameter β. Our numerical study shows that when −1 ≤ cos θ ≤ −0.5, where tan θ = β the results for the residues and masses are rather stable with respect to the variation of β, and we choose β = −1.
The LCSR predictions are reliable up to the range In order to calculate the decay width the LCSR predictions for the form factors need to be extrap- Fig. 1 The dependence of the differential branching ratio for the b → N * μ + μ − transition on s, at s 0 = 40 GeV 2 , and M 2 = 25 GeV 2

Fig. 2
The same as in Fig. 1, but for the b → N * τ + τ − transition olated to the whole physical region. For this purpose we use the z-series parametrization that is proposed in [24], The parametrization that best reproduces the form factors predicted by the LCSR in the region q 2 ≤ 11 GeV 2 , is given as For the pole masses we use, In Tables 1 and 2 we present the fit parameters a 0 , a 1 and a 2 that results from our numerical analysis. From the results for the form factors (see Tables 1  and 2), we observe that the LQSR-2 predictions nicely reproduce the relations among form factors given in Eq. (23). Using the definition of the form factors the differential decay width is calculated in the standard manner whose result is given below: where v = 1 − 4m 2 /q 2 is the lepton velocity, λ(1, r, s) = α is the fine structural constant, and the expressions of 1 (s) and 2 (s) are given in the Appendix-C. The dependence of the differential branching ratios on q 2 for the b → N * μ + μ − and b → N * τ + τ − decays, at s 0 = 40 GeV 2 , M 2 = 25 GeV 2 and β = −1 are presented in Figs. 1 and 2, respectively. In these figures the graphical results predicted by the LCSR-1 and LCSR-2 cases are shown together (Fig. 3).
Performing integration over s in the range 4m 2 2 , we obtain the branching ratios for the Table 3 Branching ratios for the b → N * + − , ( = e, μ, τ) decays , μ, τ ) decays in the case when long distance effects are due to the J/ψ family, and these results are given in Table 3.
It follows from these results that, especially for the e and μ channels, the branching ratios are quite large and could potentially be measurable at LHCb. The discovery of these decays would provide useful information about the inner structure of the N * baryon.
Finally, a comparison of the b → N + − and b → N * + − decays shows that the central values of the branching ratio of the former is approximately two (eight) times larger than the b → N * e + e − ( b → N * τ + τ − ) decays.

Conclusion
In this work the transition form factors of the b → N * + − decay are estimated, which is an alternative approach to extract information about the inner properties of the N * baryon. The contribution coming from the * b baryon is eliminated by constructing the sum rules with the choice of several different Lorentz structures. Using this result, we also calculate the branching ratio of the b → N + − and decays. We see that the branching ratios of the b → N * e + e − and b → N * μ + μ − seem to be large enough to be detected at LHCb.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All the data that were used in the manuscript are available in the current literature.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix A: N * distribution amplitudes
In this Appendix, we present the N * DAs, which are necessary to calculate the → N * transition form factors. The DAs of the N * baryon are defined from the matrix element 0 abc u a α (a 1 x) . The general decomposition of this matrix in terms of the DAs of the N * baryon is given below (see [12]).
The functions labeled with calligraphic letters in the above expression do not possess definite twists but they can be written in terms of the N * distribution amplitudes (DAs) with definite and increasing twists via the scalar product p·x and the parameters a i , i = 1, 2, 3. The relations between the two sets of DAs for the N * , and for the scalar, pseudo-scalar, vector, axial vector and tensor DAs for nucleons are: where, x i with i = 1, 2 and 3 are longitudinal momentum fractions carried by the participating quarks.
The explicit expressions for the DAs up to twist 6 are given as: Twist-3 DAs: Twist-4 DAs: Twist-5 DAs: Twist-6: Finally the x 2 corrections to the corresponding expressions V M 1 , A M 1 , T M 1 for the leading twist DAs V 1 , A 1 and T 1 in the momentum fraction space are given as: The expression for the axial-vector function A M(u) 1 (x 2 ) is given as: Similarly, we get for the function T M(u) The following functions are encountered to the above amplitudes and they can be defined in terms of the 8 independent parameters, namely f N * , λ 1 , λ 2 and f u where the parameters A u 1 , V d 1 , f d 1 , f u 1 and f d 2 are defined as [17], The numerical values of the parameters ϕ 10 , ϕ 11 , ϕ 20 , ϕ 21 , ϕ 22 , η 10 , η 11 and f N * /λ N * 1 , and λ N * 1 /λ N 1 are presented in Table 4 (this table is taken from [10]).

Invariant function
. I  5 ( p, q) for the structure q μ γ 5

Invariant function
Invariant function I 6 ( p, q) for the structure q μ / qγ 5

Invariant functions of the transition current (biσ μν q ν d)
In this section we present the expressions for the invariant functions I I 1 , I I 2 , I I 3 , and I I 4 for coefficients of the structures p μ / qγ 5 , q μ / qγ 5 , q μ γ 5 , and γ μ γ 5 , respectively, for the transition current (biσ μν q ν d) which appear in the sum rules. Invariant function I I 1 ( p, q) for the structure p μ / qγ 5 Invariant function I I 2 ( p, q) for the structure q μ / qγ 5 Invariant function I I 3 ( p, q) for the structure q μ γ 5 Invariant function I I 4 ( p, q) for the structure γ μ γ 5 xm N * 3 A 6 (1 + β) + 3 V 6 (1 + β)

Appendix C: Differential widths for the b → N * + − decays
In this Appendix we present the differential widths for the b → N * + − , ( = e, μ, τ ) decays. After lengthy, but straightforward calculations, for the differential rate of the b → N * + − we get,  , and λ(a, b, c) = a 2 + b 2 + c 2 − 2ab − 2ac − 2bc is the usual triangle function. For the element of the CKM matrix V tb V * td = (8.2 ± 0.6) × 10 −3 has been used [13]. The functions T 1 (s) and T 2 (s) are given as: