Exclusive semileptonic $B_c$-meson decays to radially excited charmonium and charm meson states

In the wake of recent measurements of ratios of semileptonic branching fractions: ${\cal R}_{J/\psi}, {\cal R}_D$ and ${\cal R}_{D^*}$ reported by the LHCb, BELLE and HFLAV Collaborations, we calculate invariant form factors for the exclusive semileptonic $B_c$-meson decays to radially excited charmonium and charm meson states in the full kinematical region within the framework of relativistic independent quark (RIQ)model. We evaluate the lepton mass effect in the decay processes induced by $b\to c,u$ transition at the quark level. Our predictions on branching fractions for $B_c\to \eta_c/\psi(nS)$ are found $\sim10^{-2}-10^{-4}$ and that for $B_c\to D^*(nS)$ are $\sim 10^{-4}$ in their $e^-$ decay modes, which lie within the detection accuracy of current experiment. Our predictions on branching fractions, forward backward asymmetry and asymmetry parameter are found in reasonable agreement with other model predictions; which can hopefully be tested in future experiments at LHC and Tevatron. Our predicted observable ${\cal R}$ in this sector are found comparable to other standard model (SM) predictions that violate the lepton flavor universality hinting at new physics beyond SM.


Introduction
One of the most engrossing puzzles in flavor physics in recent years has been the observed deviations of observables in the semileptonic B c -decays from the corresponding SM predictions. The measurements of observables R J/ψ : the ratio of semileptonic branching fractions reported by the LHCb Collaboration [1]: B(B c → J/ψ µν µ ) = 0.71 ± 0.17(stat) ± 0.18(syst) lie within 2σ above the range of existing SM predictions [2][3][4][5][6]. A series of measurements of identical decay modes B → a lopalmn95@gmail.com b skar09.sk@gmail.com D * lν l by BaBar [7,8], BELLE [9][10][11], and LHCb [12,13] also mark significant deviation of observables from the SM predictions. The average of the observed R D and R D * by the HFLAV Collaboration [14] which read reveal deviations from the SM predictions of R D 0.299 ± 0.011 from lattice calculation [15,16] and R D * 0.252 ± 0.003 from model calculation [17], respectively. Several experimental studies on B c -decays are cited in the literature, most of which pertain to B c decays to final mesons in their ground state only. This includes the observation of B c decays reported by the CDF Collaboration [18][19][20] and the observation of the decay modes: B c → J/ψ µν µ X [21], where X denotes any possible additional particle in the final state, yielding more precise measurement of B c lifetime: τ B c = 0.51 +0.18 −0.16 (stat)±0.03(syst)ps and its mass: M = 6.40±0.39± 0.13 GeV . The radially excited ψ(2S) and η c (2S) states in the heavy flavored charmonium sector and many charmoniumlike D 0 (2S) and D * 0 (2S) states have also been detected.
Several theoretical studies on semileptonic B c -decays to charmonium and charm meson states have been cited in the literature under different theoretical approaches. A few noteworthy among them are the potential model [22], non-relativistic quark model (NRQM) [5], the Baur-Stech-Wirbel (BSW) approach [23], QCD sum rules (QCDSR) [3,[24][25][26][27][28], nonrelativistic QCD (NRQCD) [29][30][31][32][33], studies based on the Bethe-Salpeter Equation [34,35], Bethe-Salpeter relativistic quark model [36], and improved Bethe-Salpeter method [37], relativistic quark model (RQM) [12, [38][39][40][41], relativistic constituent quark model (RQCM) on light front [42,43], arXiv:2204.04453v3 [hep-ph] 14 Jun 2022 Covariant confined quark model (CCQM) [44][45][46], Covariant light-front quark model (CLQM) [47], QCD potential model (QCDPM) [48][49][50], perturbative QCD (pQCD) [51][52][53][54][55][56][57][58][59][60][61], Lattice QCD [62], Light-cone QCD sum rules [63,64], and the light-front quark model (LFQM) [65]. The detection of the semileptonic B c -decay modes to radially excited charmonium and charm meson states, which are easier to identify in the experiment and the deviation of observables R from the SM predictions provide us necessary motivation to study these decay modes in the RIQ model framework. The RIQ model, developed by our group, has been applied extensively in the description of wide-ranging hadronic phenomena including the static properties of hadrons [66] and several decay properties such as the radiative; weak radiative, rare radiative; leptonic, weak leptonic, radiative leptonic; and non-leptonic decays. The semileptonic decays of heavy flavored mesons in their ground states have also been predicted within the framework of the RIQ model. It may be noted that our predictions of semileptonic decay modes [67][68][69][70] are based on vanishing lepton mass limit, where only the 3-vector (or space component) hadronic current form factors contribute to the decay amplitude. The scalar-(or time component) hadronic current form factor is not accessible in such descriptions. However, in the description of the semitauonic decay modes involving nonvanishing lepton mass, one needs to consider both the space-and time component of hadronic current form factors that contribute to the decay amplitude. In our recent analysis of exclusive semileptonic B c meson decays to the charmonium and charm meson in their ground states, we predict the lepton mass effect by evaluating the B c decays to the electronic as well as tauonic modes [71]. Inspired by our recent prediction on exclusive nonleptonic B c meson decays to radially excited charmonium and charmonium like states [72], we would like to extend the applicability of our model to analyse the semileptonic B c decays to the charmonium and charm mesons in their radially excited (nS) states, where n = 2, 3 and evaluate the lepton mass effect in those decay processes. In the process we intend to predict the branching fractions, ratio of branching fraction R in comparison to those obtained in our previous analysis [71]. We ignore the decay channels involving higher 4S final states since their properties are still not understood well.
The semileptonic decay modes:B c → η c (J/ψ)lν l and B C → D(D * )lν l are induced by b → clν l and b → ulν l transitions at the quark level, respectively. The kinematic range of q 2 for B c → η c (J/ψ)lν l is 0 ≤ q 2 ≤ 10 GeV 2 , and that for B c → D(D * )lν l is 0 ≤ q 2 ≤ 18 GeV 2 . In the B c -meson rest frame, the maximum recoil momenta of the final state charmonium (η c , J/ψ) and charm (D, D * ) meson states are estimated to be in the same range of magnitude as their masses. With this kinematic constraint, it is interesting to analyze here the decay modes B c → η c (J/ψ)lν l and B c → D(D * )lν l separately. Due to their simple theoretical description via a tree-level diagram in the SM, the effects of the strong interaction can be separated from the effects of the weak interaction into a set of Lorentz-invariant form factors. The study of semileptonic decays, therefore, reduces to a calculation of relevant weak form factors in a suitable phenomenological model framework.
It is worthwhile to note here following few points: (1) In some of the theoretical approaches like the pQCD and QCD sum rules, for example, the form factors are determined first, with an end point normalization at minimum q 2 (maximum recoil) or maximum q 2 (minimum recoil) and then extrapolated to the entire kinematical region using some monopole/ dipole /Gaussian ansatz. This makes the form factor estimation less reliable. In the present analysis, the relevant form factors are evaluated in the full kinematical range of q 2 , which makes the predictions of physical observables more accurate. (2) Since the decay modes considered here involve both the light and heavy leptons, we would evaluate both the vector-and scalar components of hadronic current form factors that contribute to the decay amplitude, study the lepton mass effects in such decays and predict the observables R in comparison with other SM predictions. (3) We update some input hadronic parameters according to the Particle Data Group [73] for our numerical analysis. The outline of the paper is as follows. In Sec. 2 we briefly discuss the general formalism and kinematics for semileptonic B c -meson decays. Section 3 describes the extraction of the transition form factors from the invariant transition amplitudes in the framework of the RIQ model. Section 4 is devoted to the numerical analysis of the form factors, decay rates, branching fractions, ratios of branching fractions, and in comparison with the results of other theoretical approaches. In Sec. 5 we summarize our results. In the Appendix, we briefly describe the RIQ model framework, the wave packet representation of the participating meson states and momentum probability amplitude of the constituent quark and antiquark in the meson bound state.

General formalism and kinematics
The general formalism and kinematics for exclusive semileptonic B c -decays: B c → η c /J/ψlν l and B c → D/D * lν l are described elaborately in Ref. [41,71], from which we quote here a few important expressions related to leptonic and hadronic parts of the invariant decay amplitude: where L µ and H µ are leptonic and hadronic parts, respectively: Here J h µ = V µ − A µ is the vector-axial vector current; q = c, u. We take (p, k) as four momenta of the parent(B c ) and daughter(X) meson with their respective spin state: S B c and S X and mass: M and m. k l and k ν are the four momenta of lepton and lepton neutrino, respectively. q = p − k = k l + k ν represents the four-momentum transfer. The hadronic amplitudes are covariantly expanded in terms of Lorentz invariant form factors. For (0 − → 0 − ) type transitions, the expansion is and for (0 − → 1 − ) type transitions, its covariant expansion is The angular decay distribution differential in the momentum transfer squared q 2 is obtained in the form Here L µσ and H µσ are the lepton and hadron tensor, respectively; m l is the mass of charged lepton. In our normalization, the lepton tensor L µσ is found to be It is convenient to express physical observables such as the hadron tensor on a helicity basis in which the helicity form factors are obtained in terms of the Lorentz invariant form factors. While doing the Lorentz contraction in Eq.(5) with the helicity amplitudes, one considers four covariant helicity projections ε µ (m) with m = +, −, 0 and m = t for spin 1 and spin 0 part, respectively, of the W o f f −shell involved in the decay process. In our attempt to study the lepton mass effects in the semileptonic decay modes, we need to consider the time component of the polarization ε µ (m = t) in addition to its other three components with m = +, −, 0. Using orthogonality and the completeness relations satisfied by helicity projections, the Lorentz contraction in Eq.(5) can be written as L µσ H µσ = L µ σ g µ µ g σ σ H µσ = L(m, n)g mm g nn H(m n ) with g mn = dia(+, −, −, −). Here the lepton and hadron tensors are introduced in the space of helicity components as: For the sake of convenience, we consider two frames of reference: i) thelν or lν center-of-mass frame and ii) the parent B c -meson rest frame. We evaluate lepton tensor L(m, n) in thelν or lν c.m. frame and hadron tensor H(m, n) in the B c -rest frame. The space-and-time components of the four momenta: p µ , k µ , q µ , and that of the polarization vectors: ε µ (±), ε µ (0), and ε µ (t) are also specified in the B c -rest frame from which the helicity components of the hadronic tensors are obtained in terms of Lorentz invariant form factors. For B c → (cc/ūc) s=0 transition, they are obtained as Using the space-and time components of the four momenta: (q µ , k µ l , k µ ν ) and helicity projections: ε µ (m = +, −, 0;t) in the (lν)-c.m. frame, it is straight forward to calculate the helicity representation L(m, n) of the lepton tensor (8). In the present analysis, we do not consider the azimuthal χ distribution of the lepton pair and therefore integrate over the azimuthal angle dependence of the lepton tensor. The differential (q 2 , cos θ ) distribution is then obtained in the form: dΓ dq 2 cos θ = The upper and lower signs associated with the parity-violating term in Eq.(11) refer to two cases: l −ν and l + ν, respectively. Out of seven terms in the r.h.s of Eq.(11), four terms identified as tilde ratesΓ i are linked with the lepton mass, and other terms identified as Γ i are lepton mass-independent.
Both are related via a flip factor m 2 l 2q 2 as: The tilde rates do not contribute to the decay rate in the vanishing lepton mass limit. They can be neglected for e and µ modes but they are expected to provide a sizeable contribution to the τ-modes. Therefore the tilde rates are crucial in evaluating the lepton mass effects in the semileptonic decay modes. The differential partial helicity rates dΓ i dq 2 are defined by Here H i (i = U, L, P, S, SL) represents a standard set of helicity structure function given by linear combinations of helicity components of hadron tensor H(m, n) = H m H † n : : Scalar Here we assume that the helicity amplitudes are real since the available q 2 -range: (q 2 ≤ (M − m) 2 ) is below the physical threshold q 2 = (M + m) 2 . Therefore we drop the angular terms that are multiplied by coefficients Im(H i H * j ), i = j * . Then integrating over cos θ one gets the differential q 2 distribution and finally integrating over q 2 , one obtains the total decay rate Γ as the sum of the partial decay rates : Γ i = (i = U, L, P) andΓ i (i = U, L, S, SL).
A quantity of interest is the forward-backward asymmetry A FB of the lepton in the (lν) c.m. frame which is defined as Another quantity of interest is the asymmetry parameter α * which is defined by rewriting Eq.(11) in terms of its cos 2 θ * dependence i.e dΓ ∝ 1 + α * cos 2 θ * . The asymmetry parameter α * which determines the transverse and longitudinal composition of the final state vector meson is given by: We list our predictions on helicity rates Γ i ,Γ i , A FB , and α * in Sec. IV.

Transition amplitude
The decay process physically occurs in the momentum eigenstate of the participating mesons. Therefore its field-theoretic description requires meson bound-states to be represented by appropriate momentum wave packets with suitable momentum and spin distribution between the constituent quarkantiquark pair in the corresponding meson core. In the relativistic independent particle picture of the RIQ model, the constituent quark antiquark are assumed to be in the state of independent and relativistic motion inside the meson boundstate: |B c (p b , p c ) , for example, with the momentum p b and p c , respectively. In this model, the momentum probability amplitude of the constituent quark and antiquark in the participating meson ground states are obtained by taking the momentum projections of corresponding quark orbitals. The corresponding expressions (A.6) are shown in the Appendix.
In the present context, one needs to extract equivalent model expressions for the momentum probability amplitude of constituent quark and antiquark inside the charmonium and charm mesons in their radially excited 2S and 3S states. For this, the quark orbitals for radially excited meson states are obtained by using the relevant Dirac quantum numbers and associated Laguerre polynomial in the general expression(A1) in the Appendix. Thereafter we take the momentum projection of the relevant quark orbitals and obtain the momentum probability amplitude of the quark and antiquark in the S-wave (2S and 3S) charmonium and charm meson states, respectively, in the form; and The binding energy of constituent quark and antiquark in the radially excited meson states are also obtained by solv-ing the appropriate cubic equations for the binding energy condition (A5) shown in the Appendix. After specifying the momentum probability amplitudes of constituent quarks and their binding energies, the radially excited final meson states, involved in the semileptonic B c -meson decays, are duly constructed which can be used to calculate the S-matrix elements for the decay process.
The internal dynamics at the constituent level, responsible for physically observable meson-level decay processes, are described by an otherwise unbound quark and antiquark using the usual Feynman technique. The constituent level Smatrix element S b→c,u f i , obtained from the Feynman diagram, when operated upon by the bag like operatorΛ used in the RIQ model yields meson-level S-matrix in the form Using the wave packet representation of the participating meson states, the S-matrix element for the decay process is obtained in the standard form In the B c -rest frame the hadronic amplitude H µ is (20) where E p b and E p b +k stand for the energy of the non-spectator quark of the parent and daughter meson, respectively, and S X |J h µ (0)|S B c represents symbolically the spin matrix elements of vector-axial vector current.

Weak decay form factors
For 0 − → 0 − transitions, the axial vector current does not contribute. The spin matrix elements corresponding to the non-vanishing vector current parts are obtained in the form: With the above spin matrix elements, the expressions for hadronic amplitudes (20) are compared with corresponding expressions in Eq. (3) yielding the form factors f + and f − as For (0 − → 1 − ) transitions, the spin matrix elements corresponding to the vector and axial-vector current are found separately as: With the spin matrix elements (24)(25)(26)(27), the expressions for hadronic amplitudes (20) are compared with corresponding expressions in Eq.(4). Then the model expressions of the form factors: V (q 2 ), A 0 (q 2 ), A + (q 2 ) and A − (q 2 ) are obtained in the form: With the form factors thus obtained in terms of model quantities, the helicity amplitudes and hence the decay rates for semileptonic B c -meson decays to radially excited (2S and 3S) charmonium and charm meson states are evaluated and our predictions are listed in Sec.4.

Numerical results and discussion
In this section, we present our numerical results on the semileptonic B c decay modes to S-wave charmonium and charm meson states. For numerical calculation, we need the model parameters (a, V 0 ), quark masses m q , and quark binding energies E q , which have already been fixed from hadron spectroscopy by fitting the data of heavy and heavy-light flavored mesons in their ground state as With this set of input parameters, wide-ranging hadronic phenomena have been described earlier in the framework of the RIQ model. However, for a description of the decay process, in the present context, involving radially excited final meson states, the constituent quarks in the meson bound states are expected to have higher binding energies compared to the binding energies in their ground states. We solve the cubic equation representing the binding energy condition (A5) for respective constituent quarks and obtain the binding energies of the quarks (c, u) in the radially excited 2S and 3S states of the (cc) and (ūc) systems as: For the CKM-parameters and B c -lifetime, we take their central values from the Particle Data Group (2020) [73] as: For the masses of participating mesons, taken as phenomenological inputs in the present calculation, we take their central values of the observed data from [73][74][75][76]. In the absence of observed data in the charmonium and charm meson sector, we take the corresponding predicted values from established theoretical approaches [77][78][79]. Accordingly, the updated meson masses used in our numerical analysis are listed in Table 1. With these input parameters, the Lorentz invariant form factors: (F + , F − ; A 0 , A + , A − , V ) representing decay amplitudes can be calculated from the overlapping integral of participating meson wave functions. We first study the q 2dependence of the invariant form factors in the allowed kinematic range. Since our main purpose is to study the lepton mass effect in semileptonic decay of B c meson, we plot q 2dependence of the form factors in Figs. 2 and 3 for decays in their e − , µ − and τ − -modes. We find that the behavior of form factors in e − and µ − -mode overlap in the entire kinematic range: 0 ≤ q 2 ≤ q 2 max . This is because of an insignificant change in the phase space boundary going from e − to µ − mode. As one can see here the maximal lepton energy shift is invisible at the usual scale of the plot. On the other hand for the decays in the τ − mode, the relevant form factors behave differently throughout the accessible kinematic range of q 2 min ≤ q 2 ≤ q 2 max , where q 2 min is +ve and away from q 2 → 0. The τ-phase space, as compared to the e − and µ − cases is considerably reduced and shifted to a large q 2 -region. In the present analysis, we shall therefore consider decays in their e − and τ − modes only for evaluating the lepton mass effects on the semileptonic B c meson decays. As mentioned earlier it is convenient to calculate decay amplitudes in the helicity basis in which the partial helicity rate and total decay rates are expressed in terms of helicity form factors. Using the helicity form factors it is straightforward to obtain partial helicity rates: dΓ i dq 2 (without flip) and dΓ i dq 2 (with flip) for decay processes considered here in their e − and τ − modes. For this we study the q 2 -dependence of the helicity form factors, partial helicity rates, and q 2 -spectra of semileptonic decay rates, separately in their e − and τ − modes before evaluating helicity rates. We first rescale the helicity form factors according to and for flip cases where and m 2 l 2q 2 : denotes a flip factor. In terms of the rescaled helicity form factors, the angle integrated differential q 2 -rate is expressed as: In Fig.4 we plot the q 2 −dependence of the rescaled helicity form factors h j andh j for B c → η c (2S, 3S) and B c → D(2S, 3S) decays, respectively in their e − and τ − modes. We find that the longitudinal no-flip amplitude h 0 is reduced in the low q 2 region going from e − to τ − mode. This reduction is due to the threshold like factor;  S−wave, where there is no pseudo-threshold factor to tamper the enhancement at large q 2 resulting from the time-like form factor in the helicity amplitude.
In Fig.5 we plot the q 2 -dependence of the rescaled helicity amplitudes for B c → ψ(2S, 3S) and B c → D * (2S, 3S) transitions in their e − and τ − mode. We find a significant reduction of the longitudinal no-flip amplitude h 0 in these decay modes. Contrary to B c → η c (2S, 3S)/D(2S, 3S) decays we find all the flip amplitudes are generally small compared to corresponding no-flip amplitudes. This is because of the partial wave structure of the scalar-current contribution. The effects of timelike (scalar)-current on semileptonic decay in the τ − mode are depicted in Figs.4 and 5. It may be noted here that the contributions of invariant form factors: F − (q 2 ) and A − (q 2 ) are dominant to determine the shape of the plot of flip helicity componenth t over the allowed kinematic range of q 2 . In B c → η c (2S, 3S)/D(2S, 3S) decay modes the contribution of F − (q 2 ) to the timelike helicity part is destructive for which the corresponding rescaled helicity amplitudeh t increases when F − (q 2 ) is switched off as shown in Fig.4. Since |h t | 2 contributes dominantly to the decay rate, an accurate determination of B c → η c (2S, 3S) and B c → D(2S, 3S) decays in their τ − -mode would help extract information on the sign and magnitude of the scalar- In Fig.6 we plot the q 2 -spectra for different helicity components in B c → η c (e − , τ − ) and B c → D(e − , τ − ) decay modes. As in the case of semileptonic B c decays to ground state charmonium and charm mesons [71], we also find here a considerable reduction of the longitudinal no-flip contribution in the τ − modes compared to the e − modes. Besides contribution of the longitudinal no-flip part the scalar flip components over the allowed kinematic range provides a sizable contribution. Fig.7 depicts the q 2 -spectra for different helicity components in B c → ψ(e − , τ − ) and B c → D * (e − , τ − ) modes. We also find here that the spin-flip parts are negligible compared to no-flip parts as seen in our earlier study [71]. Besides contribution of longitudinal no-flip part, the dominant contribution also comes here from the transverse no-flip parts. The helicity rates are found reduced uniformly going from e − to τ − modes, as expected. A significant contribution of the time like scalar-current viaS(τ) is obtained in the high q 2region.
In Figs. 8 and 9 we plot the total q 2 -spectra: dΓ dq 2 for B c → η c (ψ) and B c → D(D * ) decays, respectively, in their e − and τ − modes. For B c → η c (2S) decay in e − mode, the q 2 -spectra increases rapidly at q 2 → 0 to a peak and then flattens with a slow rise to another peak near q 2 2.5 GeV 2 . Thereafter it decreases slowly to zero value at q 2 7 GeV 2 . In the τ − -mode, the spectra, which begin at a positive q 2 3.3 GeV 2 away from q 2 → 0, increase slowly to a peak at q 2 5.5 GeV 2 and then decrease almost at the same rate to zero at q 2 7 GeV 2 . The τ − -mode spectra in these decays is found almost below e − -mode spectra throughout only to overcome slightly above the e − -spectra towards the end of A similar trend of q 2 -spectra is observed for B c → η c (3S) decays in their e − and τ − modes in the kinematic range of 0 ≤ q 2 ≤ 5 GeV 2 and 3.3 ≤ q 2 ≤ 5 GeV 2 , respectively. For B c → D(2S) decays, the e − mode spectra originating at q 2 → 0 rise slowly to a peak at about q 2 11.25 GeV 2 and then decrease to zero at q 2 13.3 GeV 2 . On the other hand, the τ − -mode spectra for these decays starts at q 2 3 GeV 2 , slowly rise to a peak at q 2 11.25 GeV 2 , and fall to zero thereafter at q 2 13.8 GeV 2 lying all along within its e −mode spectra. For B c → D(3S) decays the e − mode spectra rises fast near q 2 0, then rises moderately to a peak at q 2 6.75 GeV 2 and thereafter falls to zero at q 2 10.25 GeV 2 . Its τ-mode spectra originates at about q 2 3 GeV 2 away from q 2 → 0 and then rises to the peak at q 2 7.5 GeV 2 after which it decreases to zero at q 2 = 10.25 GeV 2 , lying all along within its e − -mode spectra. For B c → ψ(2S) decays we find a sharp rise of e − mode spectra to a peak at q 2 → 0 which rises thereafter to a prominent shoulder at q 2 3 GeV 2 beyond which it decreases in the region 3 ≤ q 2 ≤ 5.5 GeV 2 . Then the spectra flatten up to q 2 6 GeV 2 followed by its sharp fall to zero at about 6.7 GeV 2 . Its τ − -mode spectra originating at about q 2 3.2 GeV 2 attains a slow linear rise up to q 2 6 GeV 2 ; after which it slows down a bit only to hook up to a very high value at q 2 6.7 GeV 2 . On the other hand for B c → ψ(3S) decays in e − mode, the spectra, which begins with a sharper rise to a peak at q 2 → 0, flattens to attain a peak at q 2 2.3 GeV 2 and then decreases rapidly to zero at q 2 5 GeV 2 . In its τ − -mode, the spectra begins only at q 2 3.25 GeV 2 and increases slowly up to q 2 4.5 GeV 2 after which it attains a steep rise to a high value at about q 2 4.75 GeV 2 .
For B c → D * (2S) decays, the e − mode spectra attains a step rise to a peak at q 2 → 0. Then the spectra increases moderately up to q 2 5 GeV 2 beyond which it increases rapidly to develop a sharp shoulder only to fall rapidly to a very small value near q 2 13 GeV 2 . Its τ − mode spectra appears identical in shape through out the kinematic range lying all along within the e − -spectra as expected except at q 2 13 GeV 2 where it attains a steep rise to a very high value at q 2 13 GeV 2 . For B c → D * (3S) decays the e − - mode spectra rises swiftly at q 2 → 0 to a peak and then rises moderately till it attains a peak at about q 2 7 GeV 2 . Thereafter it decreases to a very low value at q 2 10 GeV 2 . Its τ − -mode spectra on the other hand begins also at about q 2 3 GeV 2 away from q 2 → 0 and rises afterward till it overtakes its e − -mode spectra at about q 2 9 GeV 2 . After overtaking the e − -mode spectra, it attains a steep rise to a very high value. It may be noted here that in B c → ψ(2S, 3S) and B c → D * (2S, 3S) decays, the shapes of the τ − -mode spectra are significantly reduced compared to that in the e −mode. However, towards the end of the kinematic region, the respective τ − -mode spectra are found to overtake the e − -mode spectra contrary to the naive phase space expectations.
Before evaluating the physical quantities of interest it is interesting to study the behavior of radial quark momentum distribution amplitude function of the decaying B c meson state together with those of the final S-wave charmonium and charm meson states. The behavior of momentum distribution amplitudes is shown in Fig.10. Note that the overlap regions between the momentum distribution amplitude function of the initial B c meson and final charmonium and charm meson states are obtained in the decreasing order from 1S to higher 2S and 3S final states. It is expected that contribution to decay rate should also be in decreasing order as on goes from ground to higher 2S and 3S final states.
In Tables 2 and 3 we list our results for the integrated partial helicity rates: Γ i (i = U, L, P) andΓ i (i = U, L, S, SL) as well as the total decay rates of semileptonic B c decays. The partial tilde decay rates for each decay process in the e − mode are found tiny as expected from Eq.(12) and can therefore be neglected. But corresponding rates in τ − -mode, obtained comparable to the untilde rates, can hardly be neglected. From the contributions of individual helicity rates, we predict the decay rates for each process, in its e − as well as τ − mode. Like all other model predictions, our predicted decay rates in τ − -modes are found, in general, smaller than that in e − -modes. For B c → ψ(2S) and B c → ψ(3S) transitions, our predicted decay rate in τ − mode is suppressed by a factor of ∼ 12 and ∼ 13, respectively. However, for B c → η c (2S) and B c → η c (3S) transitions the suppression is by a factor of ∼ 7 and ∼ 43, respectively compared to their respective decays in the e − mode. For the CKM suppressed transitions: B c → D(2S) and B c → D * (2S), the τ − modes are suppressed by a factor of ∼ 2 whereas for B c → D(3S) and B c → D * (3S) the suppression is by a factor of ∼ 2.5 over corresponding e − modes.
From the predicted decay rates shown in Tables 2 and 3 and B c -meson lifetime τ B c = 0.510 ps, we predict the branching fractions (BF) for semileptonic B c → η c (ψ) and B c → D(D * ) decays in their e − and τ − modes. Our results for such decays to radially excited 2S and 3S final meson states are shown in Table 4 and 5, respectively. Considering our predictions on semileptonic B c decays to charmonium and charm meson ground states [71], we find that our predicted branching fractions for decays to 1S, 2S, and 3S final states are obtained in the hierarchy: Our results for B c → η c (2S) and B c → η c (3S) decays in their e − mode are obtained ∼ 1.5 and ∼ 2.5 times down and in the τ − mode, the suppression is one and two orders of magnitude respectively, compared to our results in the final ground state charmonium. However, our predicted branching fractions for B c → ψ(2S) and B c → ψ(3S) decays in e − modes are ∼ 2 and ∼ 4 times down whereas, in τ-mode, they are 1/2 and 1 order of magnitude down, respectively compared to those obtained for decays to ground state charmonium (J/ψ). The node structure of the 2S wave function is responsible for small branching fractions. In the calculation of the overlap integral of meson wave functions, as there is no node for the initial wave functions, the contribution of positive and negative parts of the final wave function cancel each other out yielding a small branching fraction. About 3S final states, there are even more severe cancellations, which leads to still smaller branching fractions. As expected the tighter phase   Table 3 Helicity rates(in 10 −15 GeV) of semileptonic B c -meson decays into 3S charmonium and charm-meson states: space and weaker q 2 -dependence of form factors also lead to smaller branching fractions for transitions to higher excited 2S and 3S final states. Our predictions on branching fraction for decays to 2S final charmonium states in e and τ modes are in reasonable agreement with those of Refs. [26,55,67] and about 1 order of magnitude higher than the predictions in Refs. [39][40][41]. Although our predictions for decays to 3S charmonium states in their τ-modes are found over-estimated compared to other model predictions shown in Table. IV, our predictions for corresponding decay in e − -mode agree with those of [41,55]. As in all other model predictions we find our prediction for B c → ψ(2S) and B c → ψ(3S) is the largest of all decay modes considered in the present study.
Our results for B c → D(2S)/D * (2S) and B c → D(3S)/D * (3S) transitions are shown in Table 5. We find the branching fractions for decays in their e − as well as τ − modes in the decreasing order as one goes from 1S to higher 2S and 3S final charm meson states as expected. For B c → D decays in their e − as well as τ − mode, the predicted branching fractions decrease marginally. However, for B c → D * decays in e − mode our results decreases marginally from 1S to higher 2S and 3S states, but in τ − mode our predictions for 2S and 3S final states decrease ∼ 3 and ∼ 5 times, respectively, as compared to that for 1S final state. Using appropriate helicity rates for different decay modes, we evaluate two quantities of interest: the forward-backward asymmetry A FB (14) and asymmetry parameter α * (15). Our predicted A FB and α * for decays to 2S and 3S final meson states are shown in Table 6 and 7, respectively. For decays to spin 0 state, A FB , which is proportional to the helicity amplitudeSL, is found tiny in e − mode but non-negligible in τ − mode. Looking at the expression for A FB (14) it is easy to relate corresponding parameters in the e − and e + mode in the form: A FB (e − ) = −A FB (e + ) and A FB (τ − ) = −A FB (τ + ) for decays to spin 1 state. To determine the transverse and longitudinal components of the final state vector mesons in B c → ψ and B c → D * transitions in their e − and τ − modes, we calculate the asymmetry parameter α * . One can see from Tables 2 and 3 that the transverse and longitudinal pieces in B c → ψ(2S) and B c → ψ(3S) decays are found almost equal in e − mode. But the transverse component dominates over longitudinal parts in their τ − mode by a factor of ∼ 4. In the B c → D * (2S) and B c → D * (3S) decays, the transverse components dominate over longitudinal parts by a factor ∼ 2 in e − mode. However, in the τ − mode, transverse components dominate over longitudinal parts by a factor of ∼ 5 and ∼ 4 for B c → D * (2S) and B c → D * (3S), respectively. The asymmetry parameter is found close to −27% in e − mode and −13% in τ − mode for B c → ψ(2S) decays. For B c → D * (2S) decays it is found close to 17% and 12% in e − and τ − mode, respectively. In B c → ψ(3S) decays, α * is ob-  Table 5 Branching fractions(in %) of semileptonic B c decays into 2S and 3S charm meson states:  Table 6 Forward-backward asymmetry A FB and the asymmetry parameter α * for semileptonic B c -decays to 2S charmonium and charmed meson states.  Table 7 Forward-backward asymmetry A FB and the asymmetry parameter α * for semileptonic B c -decays to 3S charmonium and charmed meson states.

Decay process
tained close to −24% and −79% in their e − and τ − mode, respectively, whereas in B c → D * (3S) decays, it is found to be ∼ 4% in e − mode and −40% in τ − mode. This is because the transverse components of helicity amplitudes provide significant contributions compared to the tiny contributions of the longitudinal and scalar parts for decays in e − mode. On the other hand, the scalar flip partS of the helicity amplitude dominates over the rest part and contributes destructively in the decays in their τ − modes yielding a highly suppressed asymmetry parameter α * as low as −79% and −40% in B c → ψ(3S) and B c → D * (3S) decays, respectively.
Finally, we evaluate the observable: and our results are listed in Table 8 in comparison with other model predictions. The CKM matrix elements do not contribute to the ratio R. The uncertainties due to the model calculation and CKM parameter etc. are cancelled in the evaluation of the observable R. Therefore, contrary to other observables, the ratio R, in which the production of B c meson is cancelled totally, provides an essential test of the phenomenological model used in the description of decay processes. Our predicted R for B c decays to charmonium states are comparable to other standard model predictions as shown in Table 8

Summary and conclusion
In this paper, we study the exclusive B c -meson decays into radially excited 2S and 3S charmonium and charm meson states in the framework of the RIQ model based on an average flavor-independent confining potential in an equally mixed scalar-vector harmonic form. The invariant form factors representing decay amplitudes are extracted from the overlapping integral of meson wave functions derivable from the model dynamics. Since our main objective here is to evaluate the lepton mass effects on the semileptonic B c -meson decays, we analyze B c → η c /ψ(nS)lν and B c → D/D * (nS)lν in their e − and τ − modes separately. For this, we study the q 2 -dependence of relevant physical quantities such as the invariant form factors, helicity amplitudes, partial helicity rates, and the total q 2 -spectra for all decay processes analyzed in the present work. Before evaluating the decay rates we study the behavior of radial quark momentum distribution amplitude function of the decaying B c meson state together with those of the final S-wave charmonium and charm meson states. From the overlapping regions, one could assess the decay rates of different B c -meson decay modes involving daughter mesons from ground(1S) to 2S and 3S higher excited states. Considering contribution from different partial helicity rates: dΓ i dq 2 (i = U, L, P) and dΓ i dq 2 (i = U, L, S, SL), the total q 2 -spectrum dΓ dq 2 is obtained for each decay process, from which we predict the decay rates/branching fractions of the semileptonic B c -meson decays to their e − as well as τ −modes. Our predictions on the decay rates/ branching fractions for decay processes are found in overall agreement with other SM predictions. The decay rate for B c → ψ(2S) as well as B c → ψ(3S) transition is found largest of all. As expected, the decay rates for B c → η c /ψ(nS) decay induced by b → c transition at the quark level dominate over those for B c → D/D * (nS) decays induced by the quark level b → u transition in their respective e − as well as τ − modes. The contributions of longitudinal and transverse flip parts in the τ mode dominate over e − mode contribution as shown in Table 2. However, considering the contribution from no-flip parts, the total contribution to decay rate/branching fractions for all decays in the e − mode dominates over corresponding τ mode contribution. Like all other model predictions our predicted branching fractions for semileptonic B c decays to charmonium and charm meson ground [71] as well as 2S and 3S final states are obtained in the hierarchy: We also find that the τ − modes for all decays are suppressed in comparison with corresponding e − modes.
Using appropriate helicity rates we evaluate two quantities of interest:(1)Forward-backward asymmetry 'A FB ' and the asymmetry parameter 'α * '. A FB in the present analysis is found negligible for transitions into spin 0 state in their e − -mode but non-negligible in the τ − -mode, as expected. For transition into spin 1 state, we also find A FB (e − ) = −A FB (e + ) and A FB (τ − ) = −A FB (τ + ). We evaluate the asymmetry parameter α * , which determines the transverse and longitudinal components of the final vector meson state in B c → ψ(nS) and B c → D * (nS) transitions. We find the asymmetry parameter close to −27% in e − mode and −13% in τ − mode for B c → ψ(2S) decays. For B c → D * (2S) decays it is found close to 17% and 12% in e − and τ − mode, respectively. However, for B c → ψ(3S) decays, α * is obtained close to −24% and −79% in their e − and τ − mode, respectively, whereas in B c → D * (3S) decays, it is found to be ∼ 4% in e − mode and −40% in τ − mode. This is because the transverse components of helicity amplitudes provide significant contribution compared to the tiny contributions of the longitudinal and scalar parts for decays in e − mode. On the other hand, the scalar flip partS of the helicity amplitude dominates over the rest part and contributes destructively in the decays in their τ − modes yielding a highly suppressed asymmetry parameter α * as low as −79% and −40% in B c → ψ(3S) and B c → D * (3S) decays, respectively.
Finally, we evaluate the observable R, which corresponds to the ratios of branching fractions for transitions in e − modes to the corresponding value in τ − modes. Our predicted observable R for B c → η c (nS) and B c → ψ(nS) is found comparable to other SM predictions. The departure of the SM predictions of R from the experimental data highlights the puzzle in flavor physics and the failure of lepton flavor universality hinting at new physics beyond SM. However, in the absence of predicted data from established model approaches in the literature for B c decays to charm meson states, our predictions: R D (2S), R D (3S), R D * (2S) and R D * (3S) can be useful to identify the B c -channels characterized by clear experimental signature. Given the possibility of high statistics B c -events expected to yield up to 10 10 events per year at the Tevatron and LHC, semileptonic B c -meson decays into charmonium and charm mesons in their ground as well as excited states offer a fascinating area for future research. where, r nl = aω −1/4 q is a state independent length parameter, N nl is an overall normalization constant given by (ω q /a)(E q − m q ) = (4n + 2l − 1) (A.5) The solution of the cubic equation provides the zeroth-order binding energies of the confined quark and antiquark for all possible eigenmodes.
In the relativistic independent particle picture of this model, the constituent quark and antiquark are thought to move independently inside the B c -meson bound state with momentum p b and p c , respectively. Their momentum probability amplitudes are obtained in this model via momentum projection of respective quark orbitals (A1) in the following forms: For ground state mesons:(n = 1,l = 0) In terms of individual momentum probability amplitudes: G b (p b ) and G c (p c ) of the constituent quarks, G B c (p b , p − p b ), the effective momentum distribution function is taken in the form: in the straightforward extension of the ansatz of Margolis and Mendel in their bag model description [81]. In the RIQ model, the wavepacket representing a meson bound state |B c (p, S B c ) , for example, at a definite momentum p, and spin S B c is taken in the form [67][68][69][70][71][72] |B c (p, S B c ) =Λ (p, S B c )|(p b , λ b ); (p c , λ c ) (A.8) where |(p b , λ b ); (p c , λ c ) is the Fock space representation of the unbound quark and antiquark in a color-singlet configuration with their respective momentum and spin: (p b , λ b ) and (p c , λ c ). Hereb b † (p b , λ b ) andb c † (p c , λ c ) denote the quarkantiquark creation operator, respectively, andΛ (p, S B c ) is a bag like integral operator in the form: Here √ 3 is the effective color factor, ζ B c b c (λ b , λ c ) is the SU(6) spin-flavor coefficients for B c -meson and N(p) is the mesonstate normalization obtained in an integral form: by imposing the normalization condition B c (p)|B c (p ) = δ 3 (p − p ). The binding energy of constituent quark and antiquark in the parent and daughter meson in their ground as well as radially excited states are obtained by solving respective cubic equations representing appropriate bound state conditions(A5).