Semileptonic decays of Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} meson to S-wave charmonium states in the perturbative QCD approach

Inspired by the recent measurement of the ratio of Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} branching fractions to J/ψπ+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi \pi ^+$$\end{document} and J/ψμ+νμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi \mu ^+\nu _{\mu }$$\end{document} final states at the LHCb detector, we study the semileptonic decays of Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} meson to the S-wave ground and radially excited 2S and 3S charmonium states with the perturbative QCD approach. After evaluating the form factors for the transitions Bc→P,V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c\rightarrow P,V$$\end{document}, where P and V denote pseudoscalar and vector S-wave charmonia, respectively, we calculate the branching ratios for all these semileptonic decays. The theoretical uncertainty of hadronic input parameters are reduced by utilizing the light-cone wave function for the Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} meson. It is found that the predicted branching ratios range from 10-7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-7}$$\end{document} up to 10-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-2}$$\end{document} and could be measured by the future LHCb experiment. Our prediction for the ratio of branching fractions BR(Bc+→J/Ψπ+)BR(Bc+→J/Ψμ+νμ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\mathcal {BR}(B_c^+\rightarrow J/\Psi \pi ^+)}{\mathcal {BR}(B_c^+\rightarrow J/\Psi \mu ^+\nu _{\mu })}$$\end{document} is in good agreement with the data. For Bc→Vlνl\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c\rightarrow V l \nu _l$$\end{document} decays, the relative contributions of the longitudinal and transverse polarization are discussed in different momentum transfer squared regions. These predictions will be tested on the ongoing and forthcoming experiments.


Introduction
Recently, the LHCb Collaboration has measured the semileptonic and hadronic decay rates of the B c meson and obtained BR(B + c →J/ π + ) BR(B + c →J/ μ + ν μ ) = 0.0469 ± 0.0028(stat) ± 0.0046(syst) [1]. It is a motivation to investigate the B c meson semileptonic decays to charmonium, which are easier to identify in experiment. Indeed, both the CDF and the D0 Collaboration have measured the lifetime of the B c meson through its semileptonic decays [2][3][4]. More recently, the LHCb Collaboration gave a more precise measurement of its lifetime using semileptonic B c → J/ψμν μ X decays [5], where X denotes any possible additional particles in the final states. At the a e-mail: jindui1127@126.com quark level, the semileptonic decays of the B c meson driven by a b → c transitions, where the effects of the strong interaction can be separated from the effects of the weak interaction into a set of Lorentz-invariant form factors. It may provide us with the information as regards the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V cb and the weak B c to charmonia transition form factors.
There are many theoretical approaches to the calculation of B c meson semileptonic decays to charmonium. Some of them are: the nonrelativistic QCD [6,7], the Bethe-Salpeter relativistic quark model [8], the relativistic quark model [9][10][11], the light-cone QCD sum rules approach [12,13], the covariant light-front model [14], the nonrelativistic quark model [15], the QCD potential model [16][17][18], and the lightfront quark model [19]. The perturbative QCD (pQCD) [20,21] is one of the recently developed theoretical tools based on QCD to deal with the nonleptonic and semileptonic B decays. So far the semileptonic B u,d,s,c decays have been studied systematically in the pQCD approach [22][23][24][25]. One may refer to the review paper [26] and the references therein.
In our previous work [27,28], we analyzed the two-body nonleptonic decays of the B c meson with the final states involving one S-wave charmonium using the perturbative QCD based on k T factorization. By using the harmonicosillator wave functions for the charmonium states, the obtained ratios of the branching fractions are consistent with the data and other studies. Especially some of our predictions were well tested by the recent experiments at ATLAS [29] and LHCb [30], which may indicate that the harmonic-oscillator wave functions for S-wave charmonium work well.
In this paper, we extend our previous pQCD analysis to the semileptonic B c decay such as B c → (η c (nS), ψ(nS))lν (here l stands for the leptons e, μ, and τ ) with the radial quantum number n = 1, 2, 3, while the higher 4S charmonia are not included here since their properties are still not understood well. The semileptonic decays B c → (J/ψ, η c )lν have been studied in pQCD [31], compared to which the new ingredients of this paper are the following.
(1) Instead of the traditional zero-point wave function for the B c meson, the light-cone wave function which was well developed in Ref. [32] is employed in order to reduce the uncertainties caused by the hadronic parameters. In addition, the charmonium distribution amplitudes are also extracted from the correspond Schrödinger states for the harmonic-oscillator potential. (2) Here, the momentum of the spectator charm quark is proportional to the corresponding meson momentum. In Ref. [31], the charm quark in the B c meson carries a momentum with only the minus component. That is, its invariant mass vanishes, while the charm quark in the final states is proportional to the charmonium meson momentum and its invariant mass does not vanish. This substantial revision will render our analysis more consistent. (3) We updated some input hadronic parameters according to the Particle Data Group 2014 [33]. (4) Besides including the B c → (J/ψ, η c )lν decays, the B c → P/V (2S, 3S)lν decays are also investigated, where it is theoretically easier compared with that of nonleptonic decays. Our goal is to provide a ready reference to the existing and forthcoming experiments to compare their data with the predictions in the pQCD approach.
The paper is organized as follows. In Sect. 2 we define kinematics and describe the wave functions of the initial and final states, while the analytic expressions for the transition form factors and the differential decay rate of the considered decay modes are given in Sect. 3. The numerical results and relevant discussions are given in Sect. 4. The final section is the conclusion. The evaluation of the 3S charmonium distribution amplitudes is relegated to the appendix.

Kinematics and the wave functions
It is convenient to work at the B c meson rest frame and the light cone coordinate. The B c meson momentum P 1 and the charmonium meson momentum P 2 are chosen as [34] with the ratio r = m/M and m(M) is the mass of the charmonium (B c ) meson. The factors η ± = η ± η 2 − 1 come with the definition of the η of the form [34] with the momentum transfer q = P 1 − P 2 . When the final state is a vector meson, the longitudinal and transverse polarization vector L ,T can be written as The momentum of the valence quarks k 1,2 , whose notation is displayed in Fig. 1, is parametrized as the k 1T,2T , x 1,2 represent the transverse momentum and longitudinal momentum fraction of the charm quark inside the meson, respectively. One should note that there is no endpoint singularity in the B c meson decays and the integral is still convergent without the parton transverse momentum k 1T of B c meson in the collinear factorization. However, we here still keep it to suppress some non-physical contributions near the singularity (for example the singularity at x 1 = 0.1923 for B c → J/ψ decay). Fig. 1 The leading-order Feynman diagrams for the semileptonic decays B + c → P/V l + ν l with l = (e, μ, τ) There are three typical scales of the B c to charmonium decays: M, m, and the heavy-meson and heavy-quark mass difference¯ . These three scales allow for a consistent power expansion in m/M and in¯ /m under the hierarchy of M m ¯ . In the heavy-quark and large-recoil limits, based on the k T factorization theorem, the corresponding form factors can be expressed as the convolution of the hard amplitude with B c and charmonium meson wave functions. The hard amplitude can be treated by perturbative QCD at the leading order in an α s expansion (single gluon exchange as depicted in Fig. 1). The higher-order radiative corrections generate the logarithm divergences, which can be absorbed into the meson wave functions. One also encounters double logarithm divergences when collinear and soft divergences overlap, which can be summed to all orders to give a Sudakov factor. After absorbing all the soft dynamics, the initial and final state meson wave functions can be treated as nonperturbative inputs, which are not calculable but universal. Similar to the situation of the B meson [35], under above hierarchy, at leading order in 1/M, the B c meson light-cone matrix element can be decomposed as [36][37][38][39][40] with the unit vectors v = (0, 1, 0 T ) on the light cone. Here, we only consider the contribution from B c , while the contribution of¯ B c starting from the next-to-leading-power¯ /M is numerically neglected [41,42]. In coordinate space B c can be expressed by The distribution amplitude φ B c is adopted in the form [32] with shape parameters ω = 0.5 ± 0.1 GeV and the normalization conditions N is the normalization constant. For the charmonium meson, because of its large mass, the higher-twist contributions are important. The light-cone wave functions are obtained in powers of m/E or¯ /E where E(≈ M) is the energy of the charmonium meson. In terms of the notation in Ref. [43], we decompose the nonlocal matrix elements for the longitudinally and transversely polarized vector mesons (V = J/ψ, ψ(2S), ψ(3S)) and pseudoscalar respectively. For the distribution amplitudes of the 1S and 2S states, the same form and parameters are adopted as in [27,28]. The distribution amplitudes of 3S states will be derived in the appendix.

Form factors and semileptonic differential decay rates
The two factorizable emission Feynman diagrams for the semileptonic B c decays are given in Fig. 1. The transition form factors, , and A 0,1,2 (q 2 ) are defined via the matrix element [44,45], with 0123 = +1. In the large-recoil limit (q 2 = 0), the following relations should hold to cancel the poles: In the pQCD framework, it is convenient to compute the other equivalent auxiliary form factors f 1 (q 2 ) and f 2 (q 2 ), which are related to F + (q 2 ) and F 0 (q 2 ) by [31] Following the derivation of the factorization formula for the B → P, B → V transitions [46], we obtain these form factors as follows: with . α e and β a,b are the virtuality of the internal gluon and quark, respectively. Their expressions are The explicit expressions of the functions E ab , the scales t a,b , and the hard functions h are referred to [27]. In fact, if we take q 2 → 0, these expressions are agree with the results in Ref. [28]. At the quark level, the charged current B c → P(V )lν decays occur via the b → clν l transition. The effective Hamiltonian for the b → clν l transition is written as [47] where G F = 1.16637 × 10 −5 GeV −2 is the Fermi coupling constant and V cb is one of the CKM matrix elements. The differential decay rate of B c → Plν reads [14] d where m l is the mass of the leptons and λ(q 2 ) = (M 2 + m 2 − q 2 ) 2 − 4M 2 m 2 . Since the electron and muon are very light compared with the charm quark, we can safely neglect the masses of these two kinds of leptons in the analysis. For the channel of B c → V lν, the decay rates in the transverse and longitudinal polarization of the vector charmonium can be formulated as [14] d The combined transverse and total differential decay widths are defined as

Numerical results and discussions
In our calculations, some parameters are used as inputs, which are listed in Table 1.
As known, the pQCD results of these form factors are reliable only in the small q 2 region. For the form factors in the large q 2 region, the fast rise of the pQCD results indicates that the perturbative calculation gradually becomes unreliable. In order to extend our results to the whole physical region, we first perform the pQCD calculations to these form factors in the lower q 2 region (q 2 ∈ (0, ξ(M − m) 2 ) with ξ = 0.2(0.5) for the B c → 1S(2S/3S) transition), and then we make an extrapolation for them to the larger q 2 region (q 2 ∈ (ξ(M − m) 2 , (M −m) 2 )). There exist in the literature several different approaches for extrapolating the form factors from the small q 2 region to the large q 2 region. The three-parameter form is one of the pervasive models, where the fit function is chosen as where F i denotes any of the form factors, and a, b are the fitted parameters.
Our results of the transition form factors at the scale q 2 = 0 together with the fitted parameters a, b are collected in Table  2, where the theoretical uncertainties are estimated including three aspects.
The first kind of uncertainties is from the shape parameters ω in the initial and final states and the charm-quark mass m c . In the evaluation, we vary the values of ω within a 20 % range and m c = 1.275 GeV by ±0.025 GeV. We find that, in this work, the form factors are less sensitive to these hadronic parameters than our previous studies [27,28]. For example, the error induced by m c is just a few percent here, while in Ref. [28] this can reach 10-20 %. This can be understood from the B c meson wave function. In Ref. [28], the δ function depends strongly on the mass of charm quark which results in a relative large uncertainty. The second error comes from the decay constants of the final charmonium meson, which are shown in Table 1. Due to the low accuracy measurement of the decay width of the double photons decay of the pseudoscalar charmonia, the relevant uncertainty of F 0.+ is large. The last one is caused by the variation of the hard scale from 0.75 to 1.25 t. Most of this uncertainty is less than 10%, which means the next-to-leading-order contributions can be safely neglected. The errors from the uncertainty of the CKM matrix elements are very small, and they have been neglected.
It shows that the B c → P/V (1S, 2S) transition form factors are a bit larger than our previous calculations [27,28]. It is because, here, instead of the traditional zero-point wave function, we have used the light-cone wave function for the B c meson [32]. The shape of the leading twist distribution amplitude of the B c meson together with the final S-wave charmonium states are displayed in Fig. 2. It is easy to see that Table 1 The values of the input parameters for numerical analysis. The tensor decay constant f T V are determined through the assumption f T V m V = 2 f V m c , which has been used in [48] Mass (GeV) M Bc = 6.277 [28] m b = 4.18 [28] m c = 1.275 [28] m τ = 1.777 [12] m J/ψ = 3.097 [27] m ηc = 2.981 [27] m ψ(2S) = 3.686 [28] m ηc(2S) = 3.639 [28] m ηc(3S) = 3.940 [12] m ψ(3S) = 4.040 [12] CKM V cb = 40.9 × 10 −3 [28] V ud = 0.97425 [28] Decay constants(MeV) f Bc = 489 [27] f π = 131 [27] f J/ψ = 405 ± 14 [27] f ηc = 420 ± 50 [27] f ψ(2S) = 296 +3 −2 [28] f ηc(2S) = 243 +79 −111 [28] f ψ(3S) = 187 ± 8 [12] f ηc(3S) = 180 +27 −32 [12] Table 2 The fit parameters a, b, and the pQCD predictions of F 0,+ (0), A 0,1,2 (0), and V (0) for B c → nS(n = 1, 2, 3) decays, where the uncertainties come from the hadronic parameters including shape parameters ω in the initial and final state wave functions and charm-quark mass m c , decay constants, and the hard scale t, respectively   Table 2. Since one of the peaks of the 2S charmonium states wave function is so close to the peak of the B c meson wave function, the overlaps between them are large, which enhances the values of the B c → P/V (2S) form factors. However, due to the presence of the nodes in the 3S states wave function and the smaller decay constants, the corresponding form factors of the B c decays to the 3S states are slightly suppressed. We plot the q 2 dependence of the weak form factors with center values without theoretical uncertainties in Fig. 3 for the six decay processes in their physical kinematic range. We can see the different q 2 dependence of the form factors among the B c decays to different S-wave charmonia clearly. For example, the form factors for the B c → P/V (1S) transition have a relatively strong q 2 dependence, but those of the B c → P/V (2S/3S) transition show a little weaker q 2 dependence. In addition, most of these form factors become larger with increasing q 2 . However, this behavior is not universal. For instance, from Fig. 3 some of the form factors for B c → P/V (2S/3S) decays decreases with the increasing q 2 in the large region. A similar situation also exists in the light-front quark model [19] and in the ISGW2 quark model [49]. This behavior of the difference for the corresponding final states is the consequence of their different nodal structure in the wave functions.
Integrating the expressions in Eqs. (23) and (24) over the variable q 2 in the physical kinematical region, one obtains the relevant decay widths. Then it is straightforward to calculate the branching ratios. The results of our evaluation of the branching ratios for all the considered decays appear in Table 3 in comparison with predictions of other approaches. For the B c → P/V (1S) decays, our results are comparable to those of [6,7] within the error bars, but larger than the results from other models due to the values of the weak form factors.
For the B c → P/V (2S, 3S) decays our predictions are generally close to the light-cone QCD sum rules results of [12]. However, the relativistic quark model predictions for the B c → P/V (3S) decays in Refs. [9,10] are typically smaller, which can be discriminated by the future LHC experiments.
From Table 3, we can see the former four processes have a relatively large branching ratio (10 −2 ), while the branching ratios of the last four processes are comparatively small (10 −7 ∼ 10 −3 ). They have the following hierarchy: This is due to the tighter phase space, smaller decay constants, and the less sensitive dependence of the form factors on the momentum transfer q 2 for the higher excited state, which can be seen in Fig. 3. The combined effect above suppresses the branching ratios of the semileptonic B c decays to radially excited charmonia. For decays to higher charmonium excitations such a suppression should be more pro-nounced. In order to reduce the theoretical uncertainties from the hadronic parameters and the decay constants, we defined six ratios between the electron and tau branching ratios, i.e. Branching ratios (in units of %) of B c → P/V lν l decays evaluated by pQCD and by other methods in the literature. The errors induced by the same sources as in Table 2 Modes This work [6,7] Table 3, we obtain where the errors correspond to the combined uncertainty in the hadronic parameters, decay constants, and the hard scale. Since these parameter dependences canceled out in Eq. (28), the total theoretical errors of these ratios are only a few percent, much smaller than those for the branching ratios. In general, these ratios are of the same order of magnitude in the different approaches except the light-cone QCD sum rules [12], where it is obtained the smallest values of R(η c (3S)) = 33.3. For a more direct comparison with the available experimental data [1], we need to recalculate some of the nonleptonic B c decays by using the same wave functions and input parameters as this paper, whose results are where the errors induced by the same sources as in Table 2.
The ratios among the branching fractions are shown explicitly in Table 4, from which we can see that the ratios BR(B + c →J/ψπ + ) BR(B + c →J/ψl + ν l ) and BR(B + c →ψ(2S)π + ) BR(B + c →J/ψπ + ) are well consistent with the recent data [1,30], and also comparable with the prediction of the NRQCD [6,7]. Furthermore the latter still agree with the previous pQCD calculations [28] 0.29, although both BR(B + c → ψ(2S)π + ) and BR(B + c → J/ψπ + )  [27,28]. We now investigate the relative importance of the longitudinal ( L ) and transverse ( T ) polarizations contributions to the branching ratios of B c → V lν l decays within Region (1), Region (2), and the whole physical region, whose results and the ratios L T are displayed separately in Table 5. For light electron and muon, the regions are defined as: Region (1) Table 5, all of L T are < 1 in Region (2), which means that the transverse polarization dominates the branching ratios in this region. It can be understood as follows. For the B c → 1S, 2S decays, the form factor V as shown in Fig. 3 increase as the q 2 increases, which enhances the transverse polarization contribution in the large q 2 region, while for the B c → 3S decay, although the value of V decreases gradually with increasing q 2 , the form factor A 1 , which gives a dominant contribution to L , is significantly suppressed in the large region, and as a results the dominant contributions to the branching ratios of B c → ψ(2S) decays come from Region (1).
For B c → ψ(2S, 3S)eν e decays L is comparable with T in the whole physical region. These results will be tested by LHCb and the forthcoming Super-B experiments.

Conclusion
We calculate the transition form factors and obtain the branching ratios of the semileptonic decays of B c meson to S-wave charmonium states by employing the pQCD factorization approach. By using the light-cone wave function for the B c meson, the theoretical uncertainties from the nonperturbative hadronic parameters are largely reduced. It is found that the processes of B c to the ground state charmonium have comparatively large branching ratios (10 −2 ), while the branching ratios of other processes are relatively small owing to the phase space suppression, smaller decay constants, and the weaker q 2 dependence of the form factors. The theoretically evaluated ratio BR(B + c →J/ π + ) BR(B + c →J/ μ + ν μ ) = 0.046 +0.003 −0.002 is consistent with the recent data from LHCb. In addition, some interesting ratios among these branching fractions are discussed and compared with other studies. In general, these ratios in the different approaches are of the same order of magnitude, while there are also large discrepancies for specific decay modes. The partial branching ratios for transverse and longitudinal polarizations were investigated separately in B c → V lν l decays. We found that the transverse polarization gives a large contribution in the large q 2 region. For the semileptonic B c → ψ(2S, 3S)eν e decays the longitudinal contribution is comparable with the transverse contribution in the whole physical region. These theoretical predictions could be tested at the ongoing and forthcoming experiments.