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_{(s)}$$\end{document} to light tensor meson form factors via LCSR in HQEFT with applications to semileptonic decays

In the present work, the form factors of 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_{(s)}$$\end{document} to light P-wave tensor mesons (a2(1320)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_2(1320)$$\end{document}, K2∗(1430)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^*_2(1430)$$\end{document}, f2(1270)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_2(1270)$$\end{document}, f2′(1525)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f^\prime _2(1525)$$\end{document}) are calculated via the light cone sum rules (LCSR) in the framework of heavy quark effective field theory (HQEFT). Firstly, the expressions of form factors in terms of the light cone distribution amplitudes (DAs) of tensor mesons are derived via the LCSR at the leading order of heavy quark expansion. It is found that the penguin type form factors can be obtained directly from the corresponding semileptonic ones, which is similar to the case of S-wave mesons. Considering the light tensor meson DAs to twist-3, we give the numerical results of form factors systematically. As applications, we investigate the branching ratios, longitudinal polarization fractions and forward-backward asymmetries of relevant semileptonic decays induced by charged current and flavor changing neutral current (FCNC) separately. Our results may be tested by more precise experiments in the future.


Introduction
B (s) exclusive decays are of great importance in testing the standard model (SM) and probing new physics beyond it. Up to now, the decays with only S-wave mesons (including pseudoscalar and vector mesons) in the final states have been analyzed extensively both in theory and in experiment. Yet, the studies on decays containing P-wave mesons (including scalar, axial vector and tensor mesons) are relatively fewer. In the last decades, a large amount of such decays have been established experimentally [1], promoting corresponding theoretical investigations on these decays.
For the decay modes involving light tensor mesons, a large isospin violation has already been experimentally detected in B → ωK * 2 (1430) mode [2]. Also, B → φ K * 2 (1430) is mainly dominated by the longitudinal polarization [3,4], in a e-mail: zuoyabing@lnnu.edu.cn (corresponding author) b e-mail: cxyue@lnnu.edu.cn contrast with B → φ K * where the transverse polarization is comparable with the longitudinal one. In the quark model, the low-lying tensor meson with J PC = 2 ++ can be modeled as a constituent quark-antiquark pair with total orbital angular momentum L = 1 and spin S = 1. The observed tensor mesons a 2 (1320), K * 2 (1430), f 2 (1270), f 2 (1525) form an light flavor SU (3) 1 3 P 2 nonet [1]. Theoretical calculations of exclusive decays require the knowledge of non-perturbative QCD, which are generally parameterized in terms of decay constants, form factors and some non-factorization contributions. The form factors of B (s) to tensor mesons were initially calculated in the Isgur-Scora-Grinstein-Wise (ISGW) quark model [5] and its improved version ISGW2 [6,7]. They were also investigated via the covariant light-front quark model (LFQM) and applied to the relevant radiative decays [8,9]. Based on the two-parton light cone distribution amplitudes (DAs) of tensor mesons given in Ref. [10], the form factors of B (s) to light tensor mesons were derived from light-cone sum rules (LCSR) [11] and perturbative QCD [12], respectively. In Ref. [12], the results of form factors were also used to explore the relevant semileptonic decays induced by charged current. The form factors of B to tensor mesons were also calculated using LCSR with B-meson DAs [13]. Additionally, B s (B) → K * 2 (a 2 , f 2 )lν and relevant form factors were investigated with the approach of three-point QCD sum rules [14].
As well known, the LCSR is a well established technique in hadron physics, which is suitable to deal with the heavy to light meson exclusive decays. The LCSR adopted in Ref. [11] is based on an expansion over m T /m b , with m T and m b denoting the mass of tensor meson and bottom quark, respectively. Therefore, the neglected higher terms in expansion may have larger contributions. In addition, this approach can not be applied to the case of D (s) to light tensor meson exclusive decays.
In the framework of HQEFT, the heavy quark expansion is actually based on a small parameterΛ/m Q , withΛ and m Q representing the characteristic bounding energy of heavy hadron and heavy quark mass, respectively [15,16]. Moreover, some interesting relations between the penguin type form factors and corresponding semileptonic ones in the whole region of momentum transferred have been found in the heavy to light S-wave meson decays at the leading order of heavy quark expansion [17,18]. Therefore, it is worthwhile to calculate the heavy to light tensor meson form factors via the LCSR in HQEFT and investigate if similar relations still hold in this case. In this paper, we intend to calculate the B (s) → a 2 (1320), K * 2 (1430), f 2 (1270), f 2 (1525) form factors systematically by using the LCSR in HQEFT and, as applications, investigate the relevant semileptonic decays.
The remaining part of this paper is organized as follows. In Sect. 2, we give the definitions of form factors and derive their expressions in terms of the DAs of tensor mesons by using the LCSR in HQEFT. Numerical analysis and results of the form factors are presented in Sect. 3. Then, based on the form factors given here, we investigate the branching ratios, longitudinal polarization fractions and forward-backward asymmetries of all the relevant semileptonic decays in Sect. 4. Section 5 is our summary.

Definitions of form factors
The polarization tensors of a tensor meson can be constructed in terms of the polarization vectors of a massive vector state moving along the z-axis as [10,19] with The polarization tensors ε μν (λ) are symmetric and traceless, satisfying the divergence free condition ε μν (λ) P ν = 0 and the orthogonal condition ε μν (λ) ε * (λ )μν = δ λλ . The form factors for B (s) to light tensor meson decays are defined via the corresponding hadronic matrix elements as follows [11], where q μ = ( p H − P) μ is the momentum transferred and e * μ . ε μναβ is the total antisymmetric Levi-Civita symbol, with ε 0123 = − 1. A 1 , A 2 , A 0 , V and T 1 , T 2 , T 3 are the semileptonic and penguin type form factors, respectively. The latter are only relevant for exclusive decays induced by FCNC.
To the leading order of heavy quark expansion, the hadronic matrix elements can be written as [15,16] T (P, λ)|qΓ b|H ( p H ) and |M v are the effective heavy quark field and heavy meson state in HQEFT, respectively. From heavy quark spin-flavor symmetry, the leading order hadronic matrix elements with Here,Λ = lim m Q →∞ΛH . L 1 , L 2 , L 3 , L 4 are the leading order wave functions describing the heavy to light meson transition matrix elements in HQEFT. M v are the heavy pseudoscalar meson wave function. According to Eqs. (4)-(9), we have where v · P = m 2 H +m 2 T −q 2 2m H is the energy of light tensor meson in the final state. L i (i = 1, 3, 4) are in general different from corresponding L i as they arise from different matrix elements.

LCSR in HQEFT
In the following, we intend to calculate the leading order wave functions describing heavy to light meson transition matrix elements in HQEFT, i.e., L 1 , L 1 , L 2 , L 3 , L 3 , L 4 , L 4 by using the LCSR. The procedure is similar to the case of B (s) decaying to light vector mesons [20]. Firstly, we define the following vacuum-tensor meson correlation functions, (20) where T denotes time ordered product of operators.
On one hand, inserting a complete set of states with B (s) meson quantum numbers between the currents in Eqs. (17)-(20) phenomenologically, we obtain The effective decay constant of heavy meson in HQEFT is defined as [15,16] 0|qΓ Substituting Eqs. (8), (9), (25) into (21)-(24), we obtain where k μ = (P+q) μ −m Q v μ , denoting the residual momentum in B (s) meson. The integral terms represent contributions coming from the higher resonances. The subtraction terms "subtr."are introduced to ensure the convergence of integral terms, which vanish automatically after Borel transformation and therefore do not affect the physics results.
On the other hand, we can calculate the vacuum-tensor meson correlation functions theoretically, According to the assumption of quark hadron duality, starting from some effective threshold s 0 , the theoretical spectral densities ρ th Defining y = v · P, ω = 2v · k and performing two sequential Borel transformation on Eqs. (30)-(33), we have Calculating the vacuum-tensor meson correlation functions at the leading order of heavy quark expansion in HQEFT and substituting into Eqs. (38)-(41), we obtain Here are the light cone DAs of tensor mesons, which are defined as [10] T (P, λ)|q 1 It is found that the wave functions L 1 , L 3 and L 4 are exactly the same as L 1 , L 2 and L 3 , respectively. With this consideration, from Eqs. (10) to (16), we obtain the following relations between the semileptonic and penguin type form factors, Therefore, the penguin type form factors can be obtained directly from the corresponding semileptonic ones, which is similar to the case of S-wave mesons [15,16].

Numerical analysis and results of the form factors
With Eqs. (10)- (13) and (50)-(57) given above, we are now in a position to calculate the B (s) → a 2 (1320), K * 2 (1430), f 2 (1270), f 2 (1525) form factors. For this purpose, it needs to know the light cone DAs of tensor mesons, which have been studied via QCD sum rules in Ref. [10]. In the present work, we shall neglect the contributions of g 3 (u) and h 3 (u) and consider the light cone DAs to twist-3.
The twist-2 and twist-3 light cone DAs of tensor mesons are given by [10] For the decay constants of light tensor mesons, we quote the values calculated via QCD sum rules [10], as shown in Table 1.
Other input parameters involved are as follows [1,15,16], As shown in Eqs. (50)-(53), two free parameters, i.e., s 0 , t are involved in our calculations. s 0 is related to the threshold energy of initial heavy meson and t is the Borel parameter. According to the requirements of LCSR, we choose the region of these two parameters so that the curves of form factors become most stable. In the evaluation, we adjust s 0 and t for all relevant B → T decays consistently and the same procedure is also performed for B s → T decays. The variation of B → a 2 (1320), f 2 (1270), K * 2 (1430) form factors as functions of t for different s 0 at the maximally recoiling point (q 2 = 0) are shown in Figs. 1, 2 and 3, respectively.
Note that we have only shown the curves of the semileptonic type form factors here in that the penguin type form factors can be obtained directly from the corresponding semileptonic ones as mentioned above. Based on the curves given here, we can see that choosing s 0 = 2.0 ± 0.2GeV, t = 2.0 ± 0.5GeV for B → T decays is reasonable. Similarly, we can determine the region of s 0 and t for B s → T decays. It is found that the region for B s → T is the same as B → T decays.
As well known, the LCSR is expected to break down at large momentum transfer, i.e., q 2 ∼ m 2 b − 2m bΛ . For this reason, it is needed to extrapolate our LCSR results to the whole physical region via certain parametrization. Various parametric forms exist such as the single pole form, modified pole form and series expansion [12,[21][22][23]. Here the series expansion parametrization is adopted due to the fact that it is based on a more physical analysis. In this parametrization, the form factors are written as [21][22][23]   where . A truncation has been made after the quadratic term in z(q 2 , t 0 ). R F are the resonances in the range t − < q 2 < t + with appropriate quantum numbers. The corresponding resonance masses can be found in Ref. [24], as listed in Table 2. n is the number of subtractions necessary to render the dispersion relation finite, with n = 1 for A 0 and n = 2 for A 1 , A 2 , V . Optimizing the form factor parameters by using the relation where Using the LCSR predictions in low q 2 region (concretely, we have chosen to be 0 ≤ q 2 ≤ 6GeV 2 ) to fit the parameters r F 1 , r F 2 in Eq. (67), we obtain the behavior of form factors in the whole physical region. The B → a 2 (1320) form factors as functions of q 2 in the whole physical region are shown in Fig. 4. We can see that the A 1 decreases with q 2 , while the behaviors of A 2 , A 0 , V are opposite. For other decay modes, the form factors as functions of q 2 have similar behaviors. Numerical results are presented in Table 3. The uncertainties of form factors come form the two free parameters s 0 and t. As can be seen from Table 3 Table 4, in which the penguin type form factors are also listed for a complete analysis. From Table 4, we can find that large differences exist between the results of quark model approaches (ISGW2 [6,7], LFQM [8,9]) and those of sum rules (this work, LCSRT [11], LCSRB [13],QCDSR [14]) and pQCD [12]. The form factors predicted by ISGW2 all have negative signs and the magnitudes of A 1 (0), A 2 (0) are typically larger. For the LFQM [8,9], The form factors obtained in this work are in general some larger, yet still compatible with those results of LCSRT [11], LCSRB [13], pQCD [12] and QCDSR [14] as a whole.

B (s) to light tensor meson semileptonic decays
In this section, we shall apply the form factors obtained above to investigate the corresponding semileptonic decays. For the charged current induced semileptonic decays B (s) → T lν l , the differential decay width with respect to q 2 can be written as [12] dΓ where dΓ L dq 2 and dΓ T dq 2 denote the longitudinal and transverse differential decay width, respectively. Their expressions have the following forms, with Here m l is the mass of charged lepton and q 2 represents the momentum square of lepton pair. The branching ratio can be obtained by integrating Eq. (68) over q 2 in the whole physical region and using the lifetimes of B (s) ( denoted as τ H in the following ) as inputs, In addition, it is also meaningful to define the longitudinal polarization fraction and forward-backward asymmetry (FBA) of lepton For the lifetimes of B (s) , masses of charged leptons, CKM matrix element |V ub | and Fermi coupling constant G F , we use the latest values given by the particle data group (PDG)  Table 5.
For comparison, we also list the values of Br and f L calculated via pQCD [12]. We can see that the results given in this work are compatible with the corresponding values obtained by pQCD on the whole. The Br are at the order of O(10 −4 ). The results of Br and f L are almost identical for l = e, μ, while those for l = τ are significantly smaller. Specially, the Br for l = τ are about one fifth of those for l = e, μ. The f L are around 70% for l = e, μ and 60% for l = τ . The first and second uncertainties of Br come from the form factors and Table 4 Comparison of the B (s) → T form factors at q 2 = 0 given in this work with other groups. LCSRT and LCSRB denotes LCSR with tensor meson and B DAs, respectively

Decays
References  CKM matrix element |V ub |, respectively. The uncertainties from these two sources are at the same order of magnitude, i.e., roughly (15)(16)(17)(18)(19)(20)%. In contrast, the uncertainties of f L and c 0 are much smaller, due to the fact that these two quantities are defined via the ratios of decay widths and therefore most of the uncertainties are canceled. For the FCNC induced semileptonic decays B (s) → T ll, the differential decay width can be written as [25] dΓ with Table 5 The Br, f L and c 0 for semileptonic B (s) → T lν l decays, where the values of Br, f L calculated by pQCD [12] The longitudinal polarization fraction where the longitudinal differential decay width In addition, the FBA with For the lifetimes of B (s) , masses of charged leptons and Fermi coupling constant G F , we use the same values as those adopted in charged current induced semileptonic decays, For other FCNC induced semileptonic decays, the behaviors of FBAs are similar. As can be seen from Fig. 6, the zero-crossing points of FBAs are the same for l = e, μ and do not exist for l = τ , which are different from the case of charged current induced semileptonic decays (c.f., Fig. 5). The numerical results of branching ratios, longitudinal polarization fractions and zero-crossing points of FBAs for all B (s) → T l − l + (l = e, μ, τ ) decay modes are collected in Table 7, where the corresponding values for B 0 →K * 0 2 l − l + (l = μ, τ ) andB s → f 2 l − l + (l = μ, τ ) decays calculated by pQCD [25] are also listed for comparison. We can see that our results as a whole are compatible with the corresponding pQCD values for these decay modes. The branching ratios decrease significantly with the mass of charged lepton. For the decay modes with ee or μμ in the final states, the branching ratios are at the order of O(10 −8 − 10 −7 ), while those for the decays with τ τ are smaller by three orders of magnitude. As in the case of charged current induced semileptonic decays, the first and second uncertainties of branching ratios comes from the form factors and CKM matrix elements, respectively. The uncertainties from these two sources are at the same order of magnitude and roughly (15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)%. The decay modes with μμ in the final states have maximum longitudinal polarization fractions, which are approximately 70%. In contrast, the longitu- Table 6 The values of Wilson coefficients up to leading logarithmic accuracy [26]

Summary
By using the LCSR in HQEFT, we systematically calculated the B (s) → a 2 (1320), K * 2 (1430), f 2 (1270), f 2 (1525) form factors. We derived the expressions of form factors in terms of the light cone DAs of tensor mesons and found that the penguin type form factors can be obtained directly from the cor-responding semileptonic ones at the leading order of heavy quark expansion. Considering the DAs of tensor mesons to twist-3, we give the numerical results of form factors systematically. Large differences exist between the results of quark model approaches(ISGW2 [6,7], LFQM [8,9]) and those of sum rules (this work, LCSRT [11], LCSRB [13],QCDSR [14]) and pQCD [12]. The form factors obtained in this work are in general some larger, yet still compatible with those results of LCSRT [11], LCSRB [13], pQCD [12] and QCDSR [14] as a whole.
With the form factors given here, we investigated the relevant charged current and FCNC induced semileptonic decays separately. For the decays induced by charged current, the zero-crossing points of FBAs increase significantly with the mass of charged lepton. The branching ratios are at the order Table 7 The Br, f L and c 0 for rare B (s) → T l − l + decays, where the values calculated by LCSRB [13] and pQCD [25] ). The results of branching ratios and longitudinal polarization fractions are almost identical for l = e, μ, while those for l = τ are significantly smaller. For the FCNC induced decays, the zero-crossing points of FBAs are the same for l = e, μ and do not exist for l = τ . The branching ratios decrease significantly with the mass of charged lepton. For the decay modes with ee or μμ in the final states, the branching ratios are at the order of O(10 −8 − 10 −7 ), while those for the decays with τ τ are smaller by three orders of magnitude. The decay modes with μμ in the final states have maximum longitudinal polarization fractions. For both the charged current and FCNC induced semileptonic decays, the uncertainties of branching ratios come from the form factors and CKM matrix elements, which are at the same order of magnitude, roughly (15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)%. In contrast, the uncertainties of longitudinal polarization fractions and zero-crossing points of FBAs are much smaller, due to the fact that these two quantities are defined via the ratios of decay widths and therefore most of the uncertainties are canceled. Our results may be tested by more precise experiments in the future.