Semileptonic decays of $B_c$ meson to S-wave charmonium states in the perturbative QCD approach

Inspired by the recent measurement of the ratio of $B_c$ branching fractions to $J/\psi \pi^+$ and $J/\psi \mu^+\nu_{\mu}$ final states at the LHCb detector, we study the semileptonic decays of $B_c$ 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 $B_c\rightarrow P,V$, 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 $B_c$ meson. It is found that the predicted branching ratios range from $10^{-6}$ up to $10^{-2}$ and could be measured by the future LHCb experiment. Our prediction for the ratio of branching fractions $\frac{\mathcal {BR}(B_c^+\rightarrow J/\Psi \pi^+)}{\mathcal {BR}(B_c^+\rightarrow J/\Psi \mu^+\nu_{\mu})}$ is in good agreement with the data. For $B_c\rightarrow V l \nu_l$ 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.

Inspired by the recent measurement of the ratio of Bc branching fractions to J/ψπ + and J/ψµ + νµ final states at the LHCb detector, we study the semileptonic decays of Bc 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 , 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 lightcone wave function for the Bc meson. It is found that the predicted branching ratios range from 10 −7 up to 10 −2 and could be measured by the future LHCb experiment. Our prediction for the ratio of branching fractions is in good agreement with the data. For Bc → V lν l 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.

I. 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 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], the Bethe-Salpeter relativistic quark model [7], the relativistic quark model [8,9], the light-cone QCD sum rules approach [10,11], the covariant light-front model [12], the nonrelativistic quark model [13], the QCD potential model [14,15], and the light front quark model [16]. The perturbative QCD (pQCD) [17] 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 [18][19][20][21]. One may refer to the review paper [22] and the references therein.
In our previous work [23,24], 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 [25] and LHCb [26], 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 [27], 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. [28] 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. [27], 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 [29]. (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. II 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. III. The numerical results and relevant discussions are given in Sect. IV. The final section is theconclusion. The evaluation of the 3S charmonium distribution amplitudes is relegated to the appendix.

II. 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 [30] 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 [30] 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 end-point 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 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 [31], under above hierarchy, at leading order in 1/M , the B c meson light-cone matrix element can be decomposed as [32] with the unit vectors v = (0, 1, 0 T ) on the light cone. Here, we only consider the contribution from Φ Bc , while the contribution ofΦ Bc starting from the next-to-leading-powerΛ/M is numerically neglected [33,34]. In coordinate space Φ Bc can be expressed by The distribution amplitude φ Bc is adopted in the form as [28] φ 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. [35], we decompose the nonlocal matrix elements for the longitudinally and transversely polarized vector mesons (V = J/ψ, ψ(2S), ψ(3S)) and pseudoscalar mesons (P = η c , η c (2S), η c (3S)) into respectively. For the distribution amplitudes of the 1S and 2S states, the same form and parameters are adopted as in [23,24]. The distribution amplitudes of the 3S states will be derived in the appendix.

III. 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, F + (q 2 ), F 0 (q 2 ), V (q 2 ), and A 0,1,2 (q 2 ) are defined via the matrix element [36], 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 [27] Following the derivation of the factorization formula for the B → P , B → V transitions [37], 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 [23]. In fact, if we take q 2 → 0, these expressions are agree with the results in Ref. [24]. 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 [38] 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 → P lν reads [12] 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 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 [12] dΓ Lifetime τB c = 0.453 × 10 −12 s [23] The combined transverse and total differential decay widths are defined as

IV. NUMERICAL RESULTS AND DISCUSSIONS
In our calculations, some parameters are used as inputs, which are listed in Table I. 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 II, 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 [23,24]. For example, the error induced by m c is just a few percent here, while in Ref. [24] this can reach 10 − 20%. This can be understood from the B c meson wave function. In Ref. [24], the δ function depend 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 I. 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.25t. 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 [23,24]. 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 [28]. 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 the dashed line is broader in shape than that of the zero-point wave function (φ Bc (x) ∝ δ(x − r c )). The overlap between the initial and final state wave functions becomes larger in this work, which certainly induces larger form factors. We also can see that the form factors of the B c weak transitions to the 2S charmonium states at zero momentum transfer are comparable with the corresponding values of B c decays to the 1S charmonium states in Table II. 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 II: The fit parameters a, b, and the pQCD predictions of F0,+(0), A0,1,2(0), and V (0) for Bc → 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 mc, decay constants, and the hard scale t, respectively.  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 [16] and in the ISGW2 quark model [40]. This behavior of the difference for the corresponding final states is the consequence of their different   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 III in comparison with predictions of other approaches. For the B c → P/V (1S) decays, our results are comparable to those of [6] 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 [10]. However, the relativistic quark model predictions for the B c → P/V (3S) decays in Refs. [8] are typically smaller, which can be discriminated by the future LHC experiments.
From Table III, 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 pronounced. 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.
From our numerical values listed in Table III, 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 [10], 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 II. The ratios among the branching fractions are shown explicitly in Table IV, 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,26], and also comparable with the prediction of the NRQCD [6]. Furthermore the latter still agree with the previous pQCD calculations [24] 0.29, although both BR(B + c → ψ(2S)π + ) and BR(B + c → J/ψπ + ) are enhanced compared with the corresponding values of [23,24].
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 V. For light electron and muon, the regions are defined as: Region (1) Table V, 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 increase, 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.

V. CONCLUSION
We calculate the transition form factors and obtain the branching ratios of the semileptonic decays of B c meson to Swave 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.
We then propose the 3S states distribution amplitudes inferred from Eq. (A6), with the Φ asy (x) being the asymptotic models, which are given in [41]. Therefore, we have the distribution amplitudes for the radially excited charmonium mesons η c (3S) and ψ(3S): with the normalization conditions N c above is the color number, N i (i = L, t, V, s) are the normalization constants. f 3S and f T 3S are vector and tensor decay constants, respectively. Since the energy spectrum of a three-dimensional harmonic oscillator is given by E nl = [2(n − 1) + l + 3 2 ]ω, the value of the frequency ω can be determined by the difference between the two adjacent energy states. Here, the parameter ω ≈ (m 4S − m 3S )/2 ≈ 0.1 GeV.