1. Introduction
The investigation of weak decays of mesons composed of a heavy quark and antiquark gives a very important insight in the heavy quark dynamics. The exclusive semileptonic decay processes of heavy mesons generated a great excitement not only in extracting the most accurate values of Cabbibo-Kobayashi Maskawa (CKM) matrix elements but also in testing diverse theoretical approaches to describe the internal structure of hadrons. The great virtue of semileptonic decay processes is that the effects of the strong interaction can be separated from the effects of the weak interaction into a set of Lorentzinvariant form factors, i.e., the essential informations of the strongly interacting quark/gluon structure inside hadrons. Thus, the theoretical problem associated with analyzing semileptonic decay processes is essentially that of calculating the weak form factors.
The decay properties of the Bc meson are of special interest, since it is the only heavy meson consisting of two heavy quarks with different flavor. This difference of quark flavors forbids annihilation into gluons. As a result, the excited Bc meson states lying below the BD meson threshold undergo pionic or radiative transitions to the pseudoscalar ground state which is considerably more stable than corresponding charmonium or bottomonium states and decays only weakly. The CDF Collaboration reported the discovery of the Bc ground state in collisions already more than ten years ago [1]. However, up till recently its mass was known with a very large error. Now it is measured with a good precision in the decay channel. More experimental data on masses and decays of the Bc mesons are expected to come in near future from the Tevatron at Fermilab and the Large Hadron Collider (LHC) at CERN. The estimates of the Bc decay rates indicate that the c quark transitions give the dominant contribution while the b quark transitions and weak annihilation contribute less. However, from the experimental point of view the Bc decays to charmonium are easier to identify. Indeed, CDF and D0 observed the Bc meson and measured its mass analyzing its semileptonic and nonleptonic decays.
There are many theoretical approaches to the calculation of exclusive Bc semileptonic decay modes. Some of them are: QCD sum rules [2-4], the relativistic quark model [5-7] based on an effective Lagrangian describing the coupling of hadrons to their constituent quarks, the quasipotential approach to the relativistic quark model [8-10], the instantaneous nonrelativistic approach to the Bethe-Salpeter (BS) equation [11], the relativistic quark model based on the BS equation [12,13], the QCD relativistic potential model [14], the relativistic quark-meson model [15], the nonrelativistic quark model [16], the covariant light-front quark model [17], and the constituent quark model [18-21] using BSW (Bauer, Stech, and Wirbel) model [22] and ISGW (Isgur, Scora, Grinstein, and Wise) model [23]. The purpose of this paper is to extend a QCD potential model [24] with coulombic part as perturbation to calculate the hadronic form factors and decay widths for the exclusive semileptonic decay of Bc meson.
Recently, we have reported the slope and curvature of I-W function for D and B mesons with the coulombic part of the potential as the perturbation in two particular renormalisation schemes and V scheme [24]. Instead of using a particular renormalisation schemes, here in this manuscript we use the strong coupling constant as a scale dependent parameter and compute the slope and curvature of I-W function for Bc meson with a different set of mass input parameters than that of reference [24]. We then use the I-W function to study the form factors and decay rates of Bc meson into its S wave charmonium states within the framework of the potential model.
The rest of the paper is organised as follows: Section 2 contains the formalism with its subsections containing the model wavefunction, Masses, hadronic form factors and decay widths. In Section 3 we place our results and conclusions.
2. Formalism
2.1. The Wavefunction
The wavefunction computed by Dalgarno method [25,26] with coulombic part of the potential as perturbation and linear part as parent has been reported in reference [24] and the alternate approach of choosing the linear part as perturbation has been reported earlier [27,28]. For the calculations on the hadronic matrix elements, only the long-distance behaviours of the wavefunctions are important or dominant, so that taking the Coulomb piece as a perturbation becomes legitimate. The wavefunction with linear part as parent becomes an Airy function, which in fact is a diverging function. So in this manuscript we consider the wave function upto the order. For completeness we summarise the main equations in this section.
The total wave function corrected upto first order with normalisation is [24]
(1)
(2)
where is the normalisation constant for the total wave function with subscript “coul” means coulombic potential as perturbation.
Where are given as [29,30]:
(3)
(4)
(5)
(6)
(7)
and
(8)
here b and c are the model input parameters as is used in our previous works [24,27,28]. “n” is the principal quantum no. (n = 1 for ground state), is the reduced mass of mesons and is the strong running coupling constant.
2.2. Pseudoscalar and Vector Form Factors
In the case of the final states corresponds to the, states, as the matrix element of any axial current between the two pseudoscalar mesons vanishes, only vector current contributes. Unlike in the case of electromagnetic current of the charged pions, here the vector current is not conserved as. So the matrix element of the hadronic current, between the two mesons is expressed in terms of two form factors as
(9)
where is the four momentum transfer and and are the dimensionless weak transition form factors corresponds to, which are functions of the invariant. Here varies within the range
.
The transition between the pseudoscalar Bc and the vector mesons depends on four independant form factors as,
(10)
(11)
In the present study we treat Bc system similar to D system as the ratio of the constituent quark masses in the Bc meson is very close to that of D meson. So we extend the HQET of the system for the study of Bc meson also. And on the basis of HQET, the most general form of the transition discussed by Equations (9) and (10) can be expressed as [31],
(12)
(13)
(14)
where is the Isgur-Wise function. For small, nonzero recoil, it is conventional to write the Isgur-Wise function as [32]:
(15)
where is given by,
(16)
The quantity is the slope of I-W function at and known as charge radius:
(17)
The second order derivative is the curvature of the I-W function known as convexity parameter:
(18)
For the heavy-light flavor mesons the I-W function can also be written as [27,33]:
(19)
where
(20)
In Equation (1), the strong coupling constant connected to the potential is a function of the momentum as
(21)
where is the number of flavour and is the renormlistion scale related to the constituent quark mass and is the QCD scale which is taken as 0.150 GeV by fixing at the Z boson mass (91 GeV). We use the most common renormalisation scale as
with [34] and then evaluate. For numerical calculation, we use the model input parameters as, GeV and.
Using Equations (1), (15) and (19), we compute the slope and curvature of I-W function as
and
Consequently, the form factors correspond to the final state are related to the Isgur-Wise function as [31]
(22)
and those related to the as the final hadronic state are given by
(23)
It is evident from Equation (23) that at
(24)
Thus knowing the masses and Isgur-Wise function of the transition, we will be able to compute respective form factors correspond to transitions.
It should be noted that by virtue of transversality of the lepton current in the limit, the probabilities of semileptonic decays into and are independent of. Thus, in calculation of these particular decay modes of Bc meson this form factor can be consistently neglected [3].
The differential semileptonic decay rates can be expressed in terms of these form factors by
1) decay ()
(25)
2) decay () The decays rate in transversally and longitudinally polarized vector mesons are defined by [35]
(26)
(27)
where helicity amplitudes are given by the following expressions
(28)
(29)
Thus the total decay rate is given by
(30)
where is the Fermi constant, is CKM matrix element (),
Integration over of these formulas gives the total rate of the corresponding semileptonic decay. The computed results for the form factors and the total decay rate is shown in Tables 1-3 respectively and compared with the other theoretical values.
3. Results and Discussion
In this paper, we have computed the slope and curvature of Isgur-Wise function for Bc meson, considering the coulombic part of the cornell potential as perturbation and compute the different formfactors. We summarise and comment the result in the present work as follows.
1) We relate the form factors to the Isgur-Wise function in equations [22,23] which are based on the heavy flavour symmetry and is broken in the case of mesons containing two heavy quarks [39]. Spin symmetry breaking effects can occur when the c-quarks recoil momentum is larger than mc. However we expects that the equations are applicable to other kinematic point since the recoil momentum state is small () due to its heavy mass [40]. We nave computed the hadronic form factors of Bc meson using HQET in analogy to D meson system and are compared with other available results which are shown in Tables 1 and 2.
2) The values of and C are found to be smaller than the other theoretical values for D and B mesons. The reason is presumably due to the truncating of higher terms in the Airy function (which is an infinite series).
3) Numerical solution of Schrodinger equation with the specific potential also gives more accurate result than the present one. However, such approach appears to lack of physical insight into the problem unlike the relatively crude potential model approach pursued here.
4) The computed decay width for in this approach is found to be in good agreement with other theoretical results but the corresponding one for overshoots the other results and is shown in table 3. The branching ratio obtained for this semileptonic transition becomes. This value is larger than the inclusive branching ratio reported in the PDG 2011 [41] , which in principle includes the exclusive channel under consideration.