Study of Bc+ decays to the K+K−π+ final state by using Bs0, χc0 and D0 resonances and weak annihilation nonresonant topologys

In this research the weak decay of Bc+ decays to the K+K−π+ final state, which is being observed by LHCb collaboration for the first time, is calculated in the quasi-two-body decays which takes the Bs0, χc0 and D0 resonances and weak annihilation nonresonant contributions into account. In this process, the Bc+ meson decays first into Bs0π+, χc0π+ and D0π+ intermediate states, and then the Bs0, χc0 and D0 resonances decay into K+K− components, which undergo final state interaction. The mode of the Bc+ → D0(→K−π+)K+ is also associated with the calculation, in this mode the intermediate resonance D0 decays to the K−π+ final mesons. The resonances Bs0, χc0 and D0 effects in the Bc+ → Bs0(→K+K−)π+, Bc+ → χc0(→K+K−)π+ and Bc+ → D0(→K+K−)π+, D0(→ K−π+)K+ decays are described in terms of the quasi-two-body modes. There is a weak annihilation nonresonant contribution in which Bc+ decays to the K+K−π+ directly, so the point-like 3-body matrix element K+K−π+ud¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\langle {K}^{+}{K}^{-}{\pi}^{+}\left|u\overline{d}\right|0\right\rangle $$\end{document} is also considered. The decay mode of the Bc+→K¯∗0892K+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {B}_c^{+}\to {\overline{K}}^{\ast 0}(892){K}^{+} $$\end{document} is contributed to the annihilation contribution. The branching ratios of quasi-two-body decays expand in the range of (2.12 ± 0.61) × 10−6 to (7.56 ± 1.71) × 10−6.

c decays to the K + K − π + final state, which is being observed by LHCb collaboration for the first time, is calculated in the quasi-two-body decays which takes the B 0 s , χ c0 and D 0 resonances and weak annihilation nonresonant contributions into account. In this process, the B + c meson decays first into B 0 s π + , χ c0 π + and D 0 π + intermediate states, and then the B 0 s , χ c0 and D 0 resonances decay into K + K − components, which undergo final state interaction. The mode of the B + c → D 0 (→ K − π + )K + is also associated with the calculation, in this mode the intermediate resonance D 0 decays to the K − π + final mesons. The resonances B 0 s , χ c0 and D 0 effects in the B + c → B 0 s (→ K + K − )π + , B + c → χ c0 (→ K + K − )π + and B + c → D 0 (→ K + K − )π + , D 0 (→ K − π + )K + decays are described in terms of the quasitwo-body modes. There is a weak annihilation nonresonant contribution in which B +

Introduction
B c meson is one of the most interesting mesons that can be studied at the Tevatron, the discovery of the B c was reported by the CDF collaboration in the B c → J/ψl ±ν l process at Fermilab [1]. After that, the decay mode B ± c → J/ψπ ± has been observed by CDF and D0 collaboration significance of more than 8σ and 5σ respectively [2,3]. This decay mode addition B ± c → J/ψD ± s decay have also been observed by the LHCb collaboration at the LHC center-of-mass energy 7 TeV of proton-proton collisions [4,5]. Studies of B c properties are important, because it is made of two different heavy quarks, bottom-charm antiquark-quark pair. Each of the quarks can participate in a weak interaction in which other quark participates as a spectator. For this, B c is also the only meson in which decays of both heavy quarks compete with each other, therefore a wide range of decay channels are possible. However, a significant number of these channels has not been observed yet [6]. Unlike B 0 , B + and B 0 s mesons, more than 70% of the B + c width is due to c-quark decays, in which c → s transition has been observed with B + c → B 0 s π + decays [7]. Around 20% of its width is due to the b-quark decays [8]. In charmless final states, thebc → W + →qq annihilation amplitudes account for only 10% of the B + c width [9]. The mass of the B c meson has been predicted using a variety of theoretical techniques. Nonrelativistic potential models have been used to predict a mass of the B c in the range of 6247-6286 MeV/c2 [10][11][12], by using a perturbative QCD calculation slightly higher value is found [13] and recent a B c mass prediction of 6304 ± 12 +18 −0 MeV/c2 have been provided applying lattice QCD calculations [14], in which are heavier than that of other B mesons, this suggests an expected lifetime much shorter than those as τ B + c = 513.4±11.0±5.7 fs [15,16].

JHEP12(2017)019
The decay of B + c meson to three light charged hadrons like K + , K − and π + , which is observed by LHCb collaboration for the first time [9], provides a good way to study standard model for which has a large available phase space. In this work the decay mode of B + c → K + K − π + is being studied which includes other processes such as B + c → B 0 s (→ K + K − )π + decay mediated by c → s transition, charmonium mode B + c → χ c0 (→ K + K − )π + mediated by theb →c transition, B + c → D 0 (→ K + K − )π + mediated by theb →ū andb →d transitions and finally the decay of B + c → D 0 (→ K − π + )K + mediated by theb →ū andb →s transitions. In the standard model, there is another process that can be mediated via cb → W + → ud annihilation topology, in this mode the B + c decays with no charm and beauty particles in the final or intermediate states. In the B + c region 6.0 < m(K + K − π + ) < 6.5 GeV/c 2 , the signals were fitted by authors in ref. [9] separately for regions of the phase space corresponding to the different expected contributions: (a) the annihilation region, m(K − π + ) < 1.834 GeV/c 2 , they have claimed theK * 0 (892) meson can be found in this region. The contribution of the mode B + c →K * 0 (892)(→ K − π + )K + and direct annihilation are obtained separately which the estimation ofK * 0 (892) → K − π + is lower than the direct annihilation calculation, GeV/c 2 . A concentration of events was observed by [9] around m 2 (K + K − ) ∼ 11GeV 2 /c 4 , a one-dimensional projection of m(K + K − ) shows clustering near 3.41 GeV/c 2 , close to the mass of the charmonium state χ c0 where it has the highest branching fraction into the K + K − final state [17]. The accumulation of events was also observed near m 2 (K + K − ) ∼ 29GeV 2 /c 4 close to the mass of the B 0 s meson. In this research the study of the B + c → K + K − π + via quasi-two-body decay was considered, this decay mode was observed in the B 0 s → K + K − , χ c0 → K + K − , D 0 → K + K − and D 0 → K − π + channels. It is known that in the narrow width approximation, in the models we use to obtain the amplitudes of the decays, the 3-body decay rate obeys the factorization relation [18]

Decay amplitudes and branching fractions
In the standard factorization scheme, the decay amplitude is obtained by applying weak Hamiltonian which is given below:  where the weak current J µ is given by: of the d, s and b quarks, as given by the CKM matrix. Current matrix elements are defined as: It has been pointed out in the BSW2 model [19,20] that consistency with the heavy quark symmetry requires certain form factors such as F 0 and F 1 to have dipole q 2 dependence i.e. F P →P 1,2 0,1 In the following, we calculate the amplitudes of the all decay modes. Different with B u,d,s mesons, the B + c system consists of two heavy quarksb and c, which can decay individually. Here we will considerb decays while c acts as a spectator except B + c → B 0 s π + decay in which we have c → s transition. According to figure 1, at the quark level corresponding effective Hamiltonian is given by ..,10 c i Q i +c 7γ q 7γ +c 8g Q 8g ) where λ p is the CKM matrix elements, c i are the Wilson coefficients evaluated at the renor-JHEP12(2017)019 malization scale µ, Q p 1,2 are the left-handed current-current operators arising from W-boson exchange, Q 3,...,6 and Q 7,...,10 are QCD and electroweak penguin operators, and Q 7γ and Q 8g are the electromagnetic and chromomagnetic dipole operators. Because in the tree level diagrams (a 1 coefficient) bothb →ū andb →c transitions are available, so we have considered to current-current operators Q u 1 , Q u 2 , Q c 1 and Q c 2 as: The operators arise from the QCD-penguin diagrams (bothb →d andb →s transitions) which contribute in the a 4 coefficient we have considered: Here α and β are the SU(3) color indices and the subscript V − A represent the chiral projection 1 − γ 5 .
In the following, we calculate the amplitude of the B + c → B 0 s π + decay mediated by c → s transition. The amplitude of the B + c → B 0 s π + decay by using the color-allowed external W-emission tree diagram become: where F denote the factorized hadronic matrix element, which has the same definition as that in the "nonfactorizable" approach. The factorized matrix element of the and the decay constant of the π , after multiplying the form factor in the decay constant only the F [21]. Since the decay of B + c → B 0 s π + mediated by c → s transition so the m fit becomes: m D * + s . The expressions of decay amplitude for B + c → B 0 s π + within the factorization framework can be written as where a 1 = c 1 + c 2 /3. The branching fraction of B + c → B 0 s π + in B + c meson rest frame can be written as in which where τ B + c is the lifetime of the B + c meson and | p| is the absolute value of the 3-momentum of the B 0 s or π + mesons that can be calculated via: 2.2 B + c → χ c0 π + decay The same calculation, similar to the pervious calculation, can also be applied to the B + c → χ c0 π + mode, with the difference that the state of decay mediated byb →c transition so the factorized matrix element of the χ c0 (p 1 )|(cb) V −A |B + c (p) related to the form JHEP12(2017)019 forb →c transition. The amplitude and branching fraction of the B + c → χ c0 π + is the same with the B + c → B 0 s π + one which are given by and to obtain the absolute value of the 3-momentum of the χ c0 or π + mesons, the χ c0 meson mass should be replaced with the B 0 s meson mass in the pervious calculation. The branching ratios for other decays

B +
c → D 0 π + and B + c → D 0 K + decays These decays, in additionb →ū tree level a 1 coefficient, involve theb →d andb → s penguin amplitudes (a 4 coefficient) with the QCD penguins participating. The same amplitudes of the B + c → D 0 π + and B + c → D 0 K + decays are given by where a 4 = c 4 +c 3 /3 and forb →ū transition we use F B + 2.4 B 0 s → K + K − , D 0 → K + K − and D 0 → K − π + decays In the factorization approach the amplitudes of the B 0 s → K + K − , D 0 → K + K − and D 0 → K − π + decays have the following form JHEP12(2017)019

Annihilation topology
In this subsection we offer two ways for annihilation topology, the first is two body pure annihilation topology for χ c0 → K + K − decay, this decay mode proceeds only through the W-annihilation diagram and the second is nonresonant three body B + c → K + K − π + annihilation contribution. The contribution of the mode B + c →K * 0 (892)K + in the annihilation processes could be also prominent.

Pure annihilation
The χ c0 → K + K − decay is a pure annihilation decay channel, this process only occurs via annihilation between c and s quarks. In the factorization method, Feynman diagram for the χ c0 → K + K − decay is shown in figure 1. When all the basic building blocks equations are solved, for the case that both mesons are pseudoscalar, it is found that weak annihilation Kernels exhibit endpoint divergence. Divergence terms are determined by 1 0 dx/x and 1 0 dy/y. For the liberation of the divergence, a small ǫ of Λ QCD /m χ c0 order was added to the denominator. So the answer to the integral becomes ln(1 + ǫ)/ǫ form, which is shown with X A . Specifically, we treat X A = (1 + ρ A e iϕ A )ln(m χ c0 /Λ QCD ) as a arbitrary parameter obtained by using ρ A = 0.5 and a strong phase ϕ A = −55 • [22]. The factorization amplitude when both mesons are pseudoscalar is given by where b 1 is the building block of the non-singlet annihilation coefficient which is given by: The light-cone expansion implies that only leading-twist distribution amplitudes are needed in the heavy-quark limit. There exist however a number of subleading quark-antiquark distribution amplitudes of twist 3, which have large normalization factors for pseudoscalar mesons, e.g. r K χ for the kaon we use 2m 2 K / (m c − m u )(m u + m s ) 2.5.2 Effects of B + c →K * 0 (892)K + decay to the annihilation processes According to the panel (g) of the figure 1, the decay mode of the B + c →K * 0 (892)K + could be contributed to the annihilation processes. The amplitude of the B + c →K * 0 (892)K + annihilation decay reads when one of the final state mesons is vector and another is pseudoscalar, the building block of the non-singlet annihilation coefficient b 2 becomes: (8/9)πα s c 2 3(X A − 4 + π 2 /3) + rK * 0 , the ratio of the rK * 0 χ is calculated via: (2mK * 0/m b )(f ⊥ K * 0 /fK * 0).

Nonresonant annihilation contribution
for B + c → K + K − π + decay As for the three-body nonresonant annihilation matrix element K + (p 1 )K − (p 2 )π + (p 3 )| (ud) V −A |0 , we can show that it vanishes in the chiral limit owing to the helicity suppression.
To prove this claim, we first assume that the kaon and pion JHEP12(2017)019 mesons are soft.Tthe point-like 3-body matrix element is chiral-realization dependent, this realization dependence should be compensated by the pole contribution, so the three-body nonresonant annihilation matrix element becomes: (2i/f π ) p 2µ − p µ (p · p 2 )/(p 2 − m 2 K ) . It is worth stressing again that the above matrix element is valid only for low-momentum pseudoscalars. It is easily seen that in the chiral ) is needed to accommodate the fact that the final-state pseudoscalars are energetic rather than soft that is assumed to be 1/ MeV is the chiral-symmetry breaking scale. The direct three body weak annihilation amplitude reads (2.13) To obtain the p 1 · p 3 , we consider the decays of B + c meson (with 4-momentum of p and m mass) into three K + , K − and π + particles. Denote their masses m 1 , m 2 and m 3 and 4-momenta by p 1 , p 2 and p 3 , respectively. Energy-momentum conservation is expressed by p = p 1 + p 2 + p 3 . Define the following invariants s 12 = ( The three invariants s 12 , s 13 and s 23 are not independent, it follows from their definitions together with 4-momentum conservation that s 12 + s 13 + s 23 = m 2 + m 2 1 + m 2 2 + m 2 3 . We take s 12 = s and s 23 = t, so we have With these definitions, we obtain multiplying of the 4-momentum as: The decay width of a three-body process is given by where t min,max (s) = m 2 1 +m 2 3 −(1/(2s)) (m 2 −s−m 2 3 )(s−m 2 2 +m 2 1 )∓ λ(s, m 2 , m 2 3 ) λ(s, m 2 1 , m 2 2 ) , s min = (m 1 + m 2 ) 2 , s max = (m − m 3 ) 2 and λ(x, y, z) = x 2 + y 2 + z 2 − 2(xy + xz + yz).

Numerical results and conclusion
The theoretical predictions depend on many input parameters such as Wilson coefficients, the CKM matrix elements, masses, lifetimes, decay constants, form factors, and so on. We present all the relevant input parameters as follows: Wilson coefficients, the Wilson coefficients c 1 , c 2 , c 3 and c 4 in the effective weak Hamiltonian have been reliably evaluated by the next-to-leading logarithmic order. To proceed, we use the following numerical values at µ = m b scale, which have been obtained in the NDR scheme [22]: c 1 = 1.081, c 2 = −0.190, c 3 = 0.014 and c 4 = −0.036.
For the branching fraction of the B + c → D 0 K + and B + c → D 0 π + decays our calculations is obtained from B(B + c → D 0 K + ) = (5.69 + 0.74) × 10 −7 and B(B + c → D 0 π + ) = (10.53 + 1.47) × 10 −7 , as it has been obtained in [28]. The experimental result available for (f c /f u ) × B(B + c → D 0 K + ) is (9.3 +2.8 −2.5 ± 0.6) × 10 −7 and for (f c /f u ) × B(B + c → D 0 π + ) is less than 3.9 × 10 −7 [29], which implies f c /f u values in the range 0.004-0.012, the experimental branching fractions of them is obtained from 10 −4 order in which our result is a thousand times smaller than experimental one. Maybe for this reason for the mode B + c → D 0 (→ K − π + )K + , no significant deviation from the background-only hypothesis is observed [9].
For the B 0 s → K + K − mode, the decay amplitude is calculated at leading power in Λ QCD /m b and at next-to-leading order in α s using the QCD factorization approach. The calculation of the relevant hard-scattering kernels is completed. Important classes of power corrections, including "chirally-enhanced" terms and weak annihilation contributions, are estimated and included in the phenomenological analysis as given in refs. [30,31]. Applying other contributions from the Feynman graph to this decay (such as the contribution of tree level a 2 in addition to the a 1 mentioned in the text, penguin diagrams a 3 and a 4 and weak JHEP12(2017)019 annihilation of b 1 ) we obtain B(B 0 s → K + K − ) = (2.16 ± 0.35) × 10 −5 , this is while its experimental value is B(B 0 s → K + K − ) = (2.50 ± 0.17) × 10 −5 [17]. For the branching ratios of the D 0 → K + K − and D 0 → K − π + decays we obtain B(D 0 → K + K − ) = (3.78 ± 0.14) × 10 −3 and B(D 0 → K − π − ) = (4.56 ± 0.12)% where the experimental results of them are (4.01 ± 0.07) × 10 −3 and (3.93 ± 0.04)% [17].
For the pure annihilation mode of the χ c0 → K + K − theoretical calculation via factorization approach is obtained from 10 −11 order, while the experimental value of that is from 10 −3 order. This is a strange result because the χ c0 meson has very small lifetime of 10 −23 s. The fact is that the branching fraction of the pure annihilation decays by using the factorization approaches becomes smaller than experimental one [32,33], from the theoretical calculation it seems that in the decay of χ c0 → K + K − , before the K + and K − mesons are produced in the final states, the pare mesons such as D + s , D − s and D 0 ,D 0 are produced in the intermediate state so the final state interaction effects are needed.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.