Semileptonic decays in the PQCD approach with the lattice QCD input

We study the semileptonic decays in the framework of the Standard Model (SM), by employing the perturbative QCD (PQCD) factorization formalism combined with the lattice QCD input for the relevant transition form factors. We calculate the branching ratios with , the ratios of the branching fractions and , and the physical observables , and . The “PQCD+Lattice” predictions for and agree with the available experimental measurements within errors. For the ratios and , the "PQCD+Lattice" predictions agree with the other predictions. For and , our theoretical predictions agree with the measured values within errors. Our theoretical predictions of the semileptonic decays considered could be tested in the near future by the LHCb and Belle II experiments.

The studies of semileptonic decays of meson play an essential role in testing the Standard Model (SM) as well as in the search for new physics (NP) beyond SM, since the lepton flavor universality (LFU) can be examined in this kind of decay modes. LFU implies that all electroweak gauge bosons ( , and ) have equivalent couplings to the three generations of leptons, and the only difference arises from the mass differences: . Therefore, if the experiments were to find signals of lepton flavor violation (LFV), that would be a true challenge for SM.
Besides the ratios , other physical observables, such as the longitudinal polarization of the tau lepton and the fraction of the longitudinal polarization , were measured recently by the Belle collaboration [20,21,30]: F L (D * ) = 0.60 ± 0.08(stat.) ± 0.04(syst.) [30]. They are compatible with the SM predictions [31], [32] and [33] for the decays. It is clear from the above points that although the anomalies of have become less serious recently, there is still a sizeable discrepancy with respect to the SM expectations, and that they should be investigated with complementary and more precise measurements in order to draw a conclusion. In addition to the decays, the decay mode is one of the best choices for cross examination. Systematic theoretical studies and experimental measurements of the semileptonic decays of meson are therefore very important and should be performed promptly. As illustrated by the Feynman diagrams in Fig. 1, the decays are closely related to the decays via the SU(3) flavor symmetry. They are all controlled by the same transitions at the quark level, but with a different spectator quark: or s quark. As a consequence, in the limit of the SU(3) flavor symmetry, these two sets of semileptonic decays must have very similar properties. If the current anomalies in decays are indeed induced by contributions from NP, they must appear in the decays. It is therefore necessary and very interesting to study the decays and to measure the corresponding ratios and other relevant physical observables, such as , and , in order to check if similar deviations exist.
As is well known, the central issue for calculating such semileptonic decays in the SM framework is the estimate of the values and shapes of the form factors for the transitions. However, calculations of these form factors are not easy and they cannot be estimated reliably in the whole range of momenta carried by the lepton pair using a single method, and their extrapolation is indispensable. There are many traditional methods and approaches to estimate the form factors and provide predictions of the ratios , for example, the heavy quark effective theory (HQET) [4,5], light cone sum rules (LCSR) [34][35][36], lattice QCD (LQCD) [37][38][39][40][41] and the perturbative QCD factorization approach (PQCD) [9,10].
In recent years, the PQCD factorization approach has been used to study various kinds of semileptonic decays of mesons, as given for example in Refs. [9,10,[42][43][44][45][46][47][48]. However, like in many other theoretical approaches, the form factors evaluated using the PQCD approach are only reliable in the low region. Therefore, extrapolations must be done in order to cover the whole range of values of the form factors . In order to improve the reliability of the values and shapes of the six form factors obtained in the PQCD approach, we include the lattice QCD result at the end point so that the extrapolation from the low region to the high region is more reliable.
In Refs. [9,10], the decays were studied by employing the PQCD approach only [9], or with the inclusion of the lattice QCD input [10]. In Ref. [42], the decays were studied by employing the PQCD approach without the lattice QCD input. In this paper, we study the semileptonic decays of B and mesons simultaneously and focus on the following three tasks: For the and transitions, we evaluate the relevant form factors in the region using the PQCD factorization approach, and then include the lattice QCD result at the end point in the fitting process and extrapolate the form factors in the entire momentum region by employing the Bourrely-Caprini-Lellouch (BCL) parametrization [49,50] instead of the pole models used in Refs. [9,10,42].
In addition to the calculations of the branching ratios and the ratios and , we also calculate three physical observables (which were not considered in [9,10,42]) for the decays of B and mesons: the longitudinal polarization of the tau lepton , the fraction of longitudinal polarization , and the forward-backward asymmetry of the tau lepton in the ordinary PQCD approach and the "PQCD +Lattice " approach (i.e. the PQCD approach with the inclusion of the lattice QCD input for the form factors).
(3) We present our theoretical predictions, compare them with the available experimental measurements, or the theoretical predictions obtained using other different theories or models.
The paper is organized as follows: In Sec. 2, we briefly review the kinematics of the decays. The calculations of the form factors for the transitions are then given. We use the lattice QCD result at given by the HPQCD group [41] as input in our extrapolation. Explicit expressions for the differential decay rates and additional physical observables are also given in Sec. 2. In Sec. 3, we present the theoretical predictions of the other physical observables obtained using the PQCD approach, the "PQCD+Lattice " approach and some other models. A short summary is given in the final section.

Kinematics and the wave functions
In the PQCD approach, the tree-level Feynman diagrams for decays 1) are shown in Fig. 1. We define the meson momentum as , the meson momentum as , and the polarization vectors in the at rest frames as in Ref. [51].
while and denote the longitudinal and transverse polarization of mesons, respectively. The parameter and r are defined as: where is the momentum of the lepton pair. The momenta of the spectator quarks in and mesons are parametrized as where and are the fraction of the momentum carried by the light spectator quark in the initial meson and the final state meson , respectively.
For the meson wave functions, we use the functions from Refs. [51,52].
The normalization factor depends on the values of the parameter and decay constant through the normalization relation: . In order to estimate the uncertainties of theoretical predictions, we set the shape parameter GeV and GeV.
For and mesons, we use the same wave functions as in Ref. [ with the distribution amplitudes where we set the parameters and . From the heavy quark limit, we assume that

Form factors in the PQCD approach
The form factors of the transition are defined in the same form as in Ref. [54] In order to cancel the poles at , should be equal to . For convenience, we also define the auxiliary form factors and , which are related to and by the following relations, For vector meson in the final state, the form factors for transitions are and , which are given in the following form [54] : The analytical expressions for the form factors mentioned above are the following: where the mass ratios are , .
Chinese Physics C Vol. 44, No. 5 (2020) 053102 053102-4 where the color factor is , the mass ratios are , and the functions are in the following form The explicit expressions for the hard functions , hard scales and the Sudakov factors are given in the Appendix.

Lattice input at
The lattice QCD has advantages in calculating the form factors for large , and it is generally believed that the lattice QCD predictions close to or at are reliable. In this work, we make use of the lattice QCD predictions of all relevant form factors at the endpoint as the additional input in the fitting process, so as to make the extrapolation of these form factors from the low region to more reliable. The form factors used in the lattice QCD are parametrized as follows [55,56]: where , the velocity transfer is , and is the polarization vector of the meson. With a simple transformation, we can relate them to the form factors used in our work: where , and where .
At the end point , we have Therefore, the relevant form factors at the endpoint become: Using the formulas above and including the lattice input [41, 57] with the relation between and near the endpoint evaluated in Ref. [57]: we find the following values of the relevant form factors at the endpoint : The uncertainty of the form factors comes from the errors of , given in Eq. (34), is around , while the uncertainty of is around 2%. We set conservatively a common error of 5% for the form factors in Eq. (36), which takes into account the small variations of the central values of and in recent years [41, 55-57].

Extrapolation and differential decay rates
We know that the form factors calculated in the PQCD approach are only reliable in the low region. In order to cover the whole momentum range we have to perform an extrapolation. In order to improve the reliability of the extrapolation, we also use the values of the form factors at the endpoint evaluated in lattice QCD.
In the previous works [9,42], we used the pole model parametrization [58] to perform the fit. In this work, we use the Bourrely-Caprini-Lellouch (BCL) parametrization [49] instead, and match it to the lattice input to improve the reliability of the extrapolation. We use our PQCD predictions of the form factors at sixteen points in the interval , and the lattice QCD value as an additional input at the endpoint . Analogous to Ref. [50], we consider only the first two terms of the series in parameter z: with , and are the masses of the low-lying resonances. The optimized value of is chosen to be the same as in Ref. [50]. Since the choice of depends on the charged currents of semileptonic decays, i.e. of the or the transitions, we use the same set of as in Ref. [50], where the decays have the same quark level transitions as in this paper.
With these extrapolations, we have access to the full where is the Fermi constant. The differential decay rates of the decay mode can be expressed as [60]: where is the mass of leptons e, or , and is the phase space factor.
For the decay mode , the differential decay widths can be written as [58]: with the phase space factor . The total differential decay width is defined as:

Additional physical observables , and
Besides the decay rates and the ratios , there are three additional physical observables for the semileptonic decays: the longitudinal polarization of the tau lepton , the fraction of longitudinal polarization and the forward-backward asymmetry of the tau lepton . These physical observables are sensitive to new physics beyond SM [2,61,62].
As for the definitions of these physical observables, we follow Refs. [31,63,64]. For the longitudinal polarization, for example, the authors of Refs. [31,63] used the rest frame, where the spatial components of the momentum transfer vanish. In this coordinate τ x − z system the direction of the B momenta is along the z axis, and the momentum lies in the plane. Here, we use the same definition: where denotes the decay rate of with the lepton helicity of . The explicit expressions for in SM can be found in the Appendix of Ref. [64]. For decays, we find where the functions can be written as combinations of the six form factors defined in Eqs. (18), (19), (22)-(27): λ The phase space factor is the same as in Eqs. (40)- (42).
The fraction of longitudinal polarization is defined using the secondary decay chain of the semileptonic decays. It is given in the form [64]: and the corresponding differential decay rates have the following form: The lepton forward-backward asymmetry is more complicated since it is defined in the rest frame. The explicit expression is [64]: where is the angle between the 3-momenta of and the initial B or in the rest frame. The function is the angular coefficient which can be written as [64]: With the above definitions and formulas, we are now ready to present our numerical results and a phenomenological analysis.

Numerical results and discussion
In the numerical calculations, we use the following parameters (masses and decay constants are in units of GeV) [28,29,65,66]:

Form factors
It is obvious that the theoretical predictions of the differential decay rates and other physical observables for the considered semileptonic meson decays strongly depend on the form factors , and . Specifically, controls the process , while and play key roles in the process . The value of these form factors at and their dependence in the range contain a lot of information about the relevant physical processes.
In Refs. [9,51,53], the applicability of the PQCD approach to transitions was examined, and it was shown that the PQCD approach with the inclusion of the Sudakov effects is applicable to the semileptonic decays in the low region. Therefore, we present our predictions of the form factors at in Table 1.
It can be seen from Table 1 that the form factors for  and  transitions are very similar at , which implies that the effect of SU(3) flavor symmetry breaking is small, less than 10%. In order to cover the whole region, an extrapolation of the form factors has to be made from up to . In this work, we make the extrapolation using two different methods.
The first method is analogous to that used in Refs. [9,42]. We first calculate the form factors at several points in the lower region using the expressions in Eqs. (22)- (27) and the definitions in Eqs. (18), (19). In the fitting process, we use the BCL parametrization formula given in Eqs. (37), (38) instead of the pole model parametrization used in Refs. [9,42].
The second method is the "PQCD+Lattice" method, similar to Ref. [10]. Since the lattice QCD results for the form factors are reliable and accurate at , we take the lattice QCD predictions of the form factors at the endpoint as an additional input in the fitting process. In order to match the lattice input, we use the BCL parametrization [49,50] for extrapolation.
In Figs. 2 and 3, we show the theoretical predictions of the dependence of the six form factors for and transitions, obtained by employing the PQCD approach (blue solid curves) and the "PQCD + Lattice " method (red dashed curves). The blue and red bands show the theoretical uncertainties. As a comparison, we also show the central values of the PQCD predictions of the form factors (orange dashed curves) with the pole model parametrization [9,42,58]. One can see from the theoretical predictions listed in Table 1 and illustrated in Figs. 2 and 3 that: The form factors for and transitions have very similar values in the whole range due to the SU(3) flavor symmetry. (2) The differences between the theoretical predictions of the form factors in the traditional PQCD approach and the "PQCD+Lattice" method are small for low and remain moderate for large , in particular for and . For , the difference is more pronounced for large . For the transition, we found ( ) in the PQCD approach ("PQCD + Lattice" method). For the transition, we found a very similar result.

Branching ratios and the ratio of Brs
Inserting the obtained form factors in the differential decay rates given in Eqs. (40)- (43), it is straightforward to integrate over the range . For the four semileptonic decays and with , the theoretical predictions of the branching ratios are the following:  Table 1. Theoretical predictions of the form factors at using the PQCD approach with BCL parametrization [49,50].
where the theoretical errors come mainly from the uncertainties of the input parameters GeV, and GeV. Errors from the uncertainties of the decay constants of the initial and the final mesons are small and have been neglected. In Table 2, we list our PQCD and "PQCD+ B 0 s Lattice " predictions of the branching ratios of the semileptonic decays of , where the total errors are obtained by adding the individual errors in quadrature. As a comparison, we also show predictions of some other theoretical approaches in the SM framework. Unfortunately, there are no experimental results available at present. In Table 3, we show the theoretical predictions of the ratios Fig. 2. (color online) Theoretical predictions of the dependence of the six form factors for transitions in the PQCD approach (blue solid curves) and the "PQCD+Lattice" method (red dashed curves) using the BCL parametrization [49,50]. The orange dashed curves denote the PQCD predictions using the traditional pole model [9,42,58].
Chinese Physics C Vol. 44, No. 5 (2020) 053102 of the branching fractions and defined in the same way as : where l denotes an electron or a muon.

B s
Following the same procedure as for decays, we calculate the branching ratios and the ratios for decays. The numerical results are listed in Table 4 for the branching ratios and in Table 5 for the ratios and . In the estimate of the errors, 0.04 GeV was used. From Tables 2 -5, the following points can be noted: For the semileptonic decays with , the "PQCD+Lattice " predictions of the branching ratios and the ratios agree with the available experimental data within errors. This can be viewed as evidence of the reliability of the "PQCD+Lattice " method for calculating semileptonic decays of B meson.
(2) B s For the four decay modes and the eight B decay modes, the "PQCD+Lattice" predictions of the Fig. 3. (color online) Theoretical predictions of the dependence of the six form factors for transitions in the PQCD approach (blue solid curves) and the "PQCD+Lattice" method (red dashed curves) using the BCL parametrization [49,50]. The orange dashed curves denote the PQCD predictions using the traditional pole model [9,42,58].
Chinese Physics C Vol. 44, No. 5 (2020) 053102 053102-10 branching ratios are generally smaller than the conventional PQCD predictions, but the differences are relatively small, less than 20%. The predictions also agree with the previous PQCD predictions given in Refs. [9,10,42] within the still large theoretical uncertainties.
For the ratios and , the theoretical errors largely cancel, and the PQCD and "PQCD+Lattice" predictions have a very small error, only about 5%. These predictions could be tested in the near future by the LH-Cb and Belle-II experiments. For the branching ratios, the predictions of the different theoretical approaches can be somewhat different for the same decay mode, but they still agree within errors, since the theoretical uncertainties are large. For the ratios and , our "PQCD+Lattice " predictions agree with the average of the SM predictions obtained by HQET and the lattice QCD input for the form factors [3][4][5][6]29].
Comparing the two sets of ratios and , listed in Tables 3 and 5, one can see that the SU(3) flavor symmetry is well preserved in the PQCD and the "PQCD+Lattaice" approaches. This point could also be tested by experiments. Besides the branching ratios and the ratios , we consider three additional physical observables: the longitudinal polarization of the tau lepton , the fraction of longitudinal polarization and the forwardbackward asymmetry of the tau lepton . These quantities could be measured by the LHCb and Belle experiments, and may be sensitive to new physics [2, 31, Table 2. Theoretical predictions (in units of ) of the branching ratios of the semileptonic decays of , with , obtained using various theoretical approaches [42,[67][68][69][70][71][72][73].   Table 4. PQCD and "PQCD+Lattice " predictions of the branching ratios (in units of ) of the eight semileptonic decays and with . As a comparison, the previous "PQCD+Lattice" predictions [10], the SM predictions based on HQET [11], and the world average of the measured values given in PDG 2018 [65] are also shown.  [61][62][63]. Investigation of these physical observables may provide new clues for understanding the puzzle, and it is therefore necessary and interesting to calculate them.
In Refs. [31][32][33], these three physical observables were calculated in the SM framework, and possible new physics effects were studied. We calculate these observables in the PQCD and "PQCD+Lattice" approaches using Eqs. (44)- (59), and show the numerical predictions of , and in Table 6. The measured values of and [20,21,30] given in Eqs. (4), (5) are also listed in Table 6. As a comparison, we also show the SM predictions of these physical observables given in Refs. [31,32] .
From the results listed in Eqs. (4), (5) and in Table 6, the following points can be noted: The uncertainties of the theoretical predictions of the physical observables , and are very small compared to the branching ratios, since the theoretical uncertainties largely cancel in the ratios.
The PQCD and "PQCD+Lattice " predictions of the physical observables , and for decays are very similar: the difference for a fixed decay mode is less than 5%. For of the decays, however, the difference is about 40%. The reason is the definition of the angular coeffi- where the term in Eq. (59) can be moderately changed by the high behavior of the form factors in the PQCD and "PQCD+Lattice" approaches. Clearly, a precise experimental measurement of the forward-backward asymmetry in the decays could be of great help to test and improve the factorization models.
For and , our theoretical predictions agree with the measured data within errors, partially due to the still large experimental errors. For all decay modes considered, the PQCD and "PQCD+Lattice " predictions of the three physical observables are consistent with the other approaches in the SM framework. We expect that these physical observables could be measured with high precision in the future by the LHCb and Belle-II experiments, which could help to test the theoretical models or approaches.

Summary
In this paper, we studied the semileptonic decays in the SM framework by employing the conventional PQCD factorization approach and the "PQCD+Lattice " approach. In the second approach, we took into account the lattice QCD results of the relevant form factors as input for the extrapolation from the low Table 5. PQCD and "PQCD+Lattice" predictions of the ratios and . As a comparison, the previous "PQCD+Lattice" predictions given in Ref. [10], the average of the SM predictions given in Ref. [29], measured values reported by the BaBar, Belle and LHCb collaborations [1,24,25], and the world averaged results from HFLAV group [29] are also shown.  Table 6. PQCD and "PQCD+Lattice" predictions of , and . The SM predictions and the measured values are also listed.

Ratios
Observable to the endpoint . We calculated the form factors , and of the transitions, provided the theoretical predictions of the branching ratios of the semileptonic decays and of the ratios and . In addition, we also gave our theoretical predictions of the additional physical observables: the longitudinal polarization of the tau lepton , the fraction of longitudinal polarization and the forward-backward asymmetry of the tau lepton . From the numerical calculations and phenomenological analysis we noted the following points: For the twelve semileptonic decay modes considered, the "PQCD+Lattice" predictions of the branching ratios are generally smaller than the conventional PQCD predictions, but the differences are relatively small, less than 20%. The "PQCD+Lattice" predictions of the branching ratios and the ratios agree with the available experimental measurements within errors.  They agree with the other SM predictions based on different approaches. These predictions could be tested by the LHCb and Belle-II experiments.
(3) P τ (D ( * ) (s) ) F L (D * (s) ) A FB (τ) For most of the observables , and , the PQCD and "PQCD+Lattice " predictions agree, and the difference is less than 5%. The relatively large difference of in the decays can be reasonably explained. For and , our theoretical predictions agree with the measurements within errors. For all decay modes considered, the PQCD and "PQCD+Lattice" predictions of these physical observables are consistent with the results of the other approaches in the SM framework. The LHCb and Belle-II experiments could help to test the above theoretical predictions.
We wish to thank Wen-Fei Wang and Ying-Ying Fan for valuable discussions. In this appendix, we present the explicit expressions for some functions which have already appeared in the previous sections. The hard functions and come form the Fourier transform and can be written as:

Appendix: Relevant functions
The threshold resummation factor is adopted from Ref. [74]: The factor contains the Sudakov logarithmic cor- where is defined in Eq. (7). The hard scale t and the quark anomalous dimension govern the aforementioned renormalization group evolution. The explicit expressions for the functions can be found in Appendix A of Ref. [75]. The hard scales in Eqs. (22)- (27) are chosen as the largest scales of virtuality of internal particles in the hard b quark decay diagram,