-
According to the
qˉq quark model assignments for mesons, bottom-charmed mesons are unique particles consisting of two heavy quarks with different flavors, because it is generally assumed that the top quark being the heaviest element fermion of the standard model (SM) has a very short lifetime and decays before hadronization. The bottom-charmed mesons are isospin singlets, and have nonzero additive quantum numbers, Q= B= C= ±1 , where Q, B and C are respectively electric charge, bottom and charm. The bottom-charmed mesons, analogical with hidden-flavored charmoniumcˉc and bottomoniumbˉb , are often considered as nonrelativistic bound states. Hyperfine splitting interactions divide the ground bottom-charmed states into the spin singlet11S0 and triplet13S1 , corresponding to the pseudoscalarBc meson with the spin-parityJP = 0− and vectorB∗c meson withJP = 1− , respectively.The pseudoscalar
Bc meson was first discovered viaBc → J/ψπ andJ/ψππ+π− decay modes with hadronicZ0 decay sample data collected by the DELPHI detector at thee+e− collider LEP in 1997 [1]. The natural properties of theBc meson including its mass, lifetime, and spin-parity have now been well determined [2]. In the meantime, the vectorB∗c meson has not been definitively identified by experimental physicists yet. All the information on theB∗c meson comes only from theoretical calculation for the moment. Undoubtedly, it is a matter of great urgency to find and identify theB∗c meson in experiments, which is the basis and prerequisite for obtaining a deeper insight into the desiredB∗c meson, and distinguishing by different theoretical models.The production probability of bottom-charmed mesons is far less than that of heavy-light
Bu,d mesons and heavy-heavy charmonium and bottomonium. First, two heavy quark pairs of bothbˉb andcˉc should be almost simultaneously produced in a high energy process, and then a b (orˉb ) quark from abˉb pair andˉc (or c) quark from acˉc pair should be lucky enough to combine and finally form bottom-charmed mesons [3]. In positron-electron collisions, bottom-charmed mesons can be produced viaZ0 boson decay,e++e−→Z0→(bˉc)+ˉb+c,
(1) for example, the observation of
Bc mesons by DELPHI [1], ALEPH [4] and OPAL [5] detectors in LEP experiments. In hadron-hadron collisions, the bottom-charmed mesons can be produced via gluon-gluon fusion and quark-antiquark annihilation,g+g→Z0→(bˉc)+ˉb+c,
(2) q+ˉq→Z0→(bˉc)+ˉb+c,
(3) for example, the observation of
Bc meson by CDF [6] and D0 [7] at Fermilab Tevatron, by LHCb [8] in LHC experiments. More and moreB∗c mesons are expected to be accessible in the future. Given the branching ratio ofZ0 boson decayBr(Z0→bˉb) = 15.12(5)% [2] and the bottom quark fragmentation fractionf(b→B∗c) ∼ 6×10−4 [9–11], there will be more than108 B∗c mesons from1012 Z0 boson decays at the Circular Electron Positron Collider (CEPC) [12], and more than109 B∗c mesons from1013 Z0 boson decays at the Future Circular Collider (FCC-ee) [13]. Assuming theB∗c production cross section is about100 nb forpp collisions at√s = 13 TeV [14], more than3×1010 B∗c mesons will be available with a total integrated luminosity of 300fb−1 at LHCb of HL-LHC [15]. The huge amount of experimental data provides an excellent opportunity and solid foundation to carefully study theB∗c mesons.Because of its nonrelativistic nature, the mass of the
B∗c meson is approximately equal to the sum of the mass of bottom and charm quarks, i.e.,mB∗c ≃ mb + mc . A recent lattice calculation givesmB∗c = 6331(7) MeV [16], which agrees basically with many other theoretical estimations with various theoretical models (e.g., references [150–207] of Ref. [17]). Owing to the hierarchical relationships among mesonic masses,mB∗c < mB + mD andmB∗c < mBc + mπ , theB∗c meson can not decay through the strong interaction. The dominant decay of theB∗c meson is the radiative transition process,B∗c → Bc + γ . This partial decay width can be written as [18],Γ(B∗c→Bcγ)=43αemk3γ|μB∗cBc|2,
(4) where the center-of-mass momentum of a photon in the rest frame of the
B∗c meson and the magnetic dipole (M1) momentum are respectively defined askγ=m2B∗c−m2Bc2mB∗c≈56MeV,
(5) μB∗cBc=⟨Bc|∑i=b,cQi2miσiz|B∗c⟩=16(2mc−1mb).
(6) It is clear that, kinematically, the parity and angular momentum conservation in electromagnetic interactions require the orbital angular momentum between the
Bc meson and photon to be L= 1 . In addition, the photon is very soft, resulting in a very small phase space. Dynamically, the decay width is proportional to the fine-structure constantαem of the electromagnetic interactions and the module square of the magnetic dipole momentum, while the magnetic dipole momentum is inversely proportional to the mass of heavy quarks. The combined effects of both kinematical and dynamical factors produce a very narrow decay width,Γ(B∗c→Bcγ) ≈ 60 eV [17]. A good approximation is the full decay width of theB∗c mesonΓB∗c ≈ Γ(B∗c→Bcγ) . Experimentally, the electromagnetic processB∗c → Bcγ with an occurrence probability of almost100% should be easily detected at thee+e− CEPC and FCC-ee colliders, thanks to the excellent photon resolution of the electromagnetic calorimeter and the fine performance in reconstruction technology and methods of the charged particle tracks. Because the masses ofB∗c andBc mesons are very close, the detection of the photon from theB∗c → Bcγ process plays a critical role in distinguishing between theB∗c andBc mesons. The identification of the soft photon is very challenging. Moreover, the photon from the M1 transitionB∗c → Bc + γ decay is bound to be seriously affected by those from bremsstrahlung radiation and chaotic electromagnetic backgrounds, which results in identification inefficiencies. Besides the electromagnetic interactions, theB∗c meson can also decay through weak interactions in the SM of elementary particles. The narrowness of the full width of theB∗c meson affords great potential for experimental investigations onB∗c weak decays. TheB∗c weak decays provide a good opportunity and an important and useful alternative to find and explore the foreseeableB∗c mesons with some novel methods in future high-luminosity and high-precision experiments.Compared with electromagnetic
B∗c decays, there are plenty ofB∗c meson weak decay processes. Based on the weak interaction couplings between particles, both component quarks of theB∗c meson, the heavy bottom and charm quarks can transmit into first and second family quarks lighter than themselves. TheseB∗c weak processes can be classified into three types, similar to those for the pseudoscalarBc meson: (1) the bottom quark decays while the charm quark remain quiescent as a spectator; (2) the charm quark decays while the bottom quark remain inactive as a spectator; (3) the bottom and charm quarks annihilate into a virtual charged W boson. The purely leptonicB∗c decays belonging to type (3) suffer additional complications from the final neutrinos. The charged hadrons fromB∗c weak decays, such as pions and kaons, are relatively easy to identify in sensitive particle detectors. In this paper, we estimate the branching ratio forB∗+c → BP ,BV ,B∗P ,ψP ,ηcP ,ηcV decays arising from external W emission with the factorization approach, where P (V) denotes the positively charged pion and kaon (ρ+ andK∗+ ). Some studies [19–22] have shown that these processes in question have relatively large branching ratios betweenB∗c meson weak decays and should have priority to be investigated experimentally. Here, we hope to provide a feasibility analysis of searching for theB∗c meson via some particular nonleptonic weak decays in future high-energy and high-luminosity experiments.The remainder of this paper is organized as follows. The effective Hamiltonian for the nonleptonic
B∗c weak decays in question is given in Section II. Branching ratios and our comments are presented in Section III. Section IV is devoted to a brief summary. The decay amplitudes are listed in Appendix. -
The effective Hamiltonian for the concerned nonleptonic
B∗c decays is written as [23]Heff=GF√2V∗cbVuq1{C1Ob1+C2Ob2}
(7) +GF√2V∗cq2Vuq3{C1Oc1+C2Oc2},
(8) where Fermi constant
GF ≈ 1.166×10−5 GeV−2 [2] and Wilson coefficientsC1,2 are process-independent couplings. Eq. (7) and Eq. (8) correspond to type (1) and (2)B∗c decays, respectively.Vcq andVuq are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [24, 25], and their magnitudes have been determined experimentally [2].|Vud|=0.97370(14),
(9) |Vus|=0.2245(8),
(10) |Vcd|=0.221(4),
(11) |Vcs|=0.987(11),
(12) |Vcb|=0.0410(14).
(13) The expressions of effective tetra-quark operators are written as
Ob1=[ˉbαγμ(1−γ5)cα][ˉuβγμ(1−γ5)q1β],
(14) Ob2=[ˉbαγμ(1−γ5)cβ][ˉuβγμ(1−γ5)q1α],
(15) Oc1=[ˉq2αγμ(1−γ5)cα][ˉuβγμ(1−γ5)q3β],
(16) Oc2=[ˉq2αγμ(1−γ5)cβ][ˉuβγμ(1−γ5)q3α],
(17) where
α andβ are the color indices,q1,2,3 = d and s.It is clear that the calculation of hadronic matrix elements (HMEs) is the remaining central missions of decay amplitudes. HMEs trap the tetra-quark operators between initial and final mesons, and involve long- and short-distance contributions. Recently, several attractive phenomenological methods based on power counting rules and QCD perturbative calculations have been developed to deal with HMEs, such as the QCD factorization (QCDF) approach built upon the collinear approximation [26–36], and the perturbative QCD (pQCD) approach where contributions from the transverse momentum of valence quarks and Sudakov factors for hadronic wave functions are taken into consideration [37–43]. The theme of both the QCDF and pQCD approaches is to properly factorize the perturbative and nonperturbative contributions contained in HMEs and appropriately evaluate their shares. The main objective of this paper is to determine whether the
B∗c meson can be explored in future experiments, so an approximate estimation of branching ratio is completely sufficient. In this sense, the naive factorization (NF) approach [44] is used and applied to nonleptonicB∗c meson weak decay in our calculation. The physics picture of the NF approach is very simple and clear. The results from the NF approach can be regarded as the leading order approximation of those from the QCDF approach. What matters is that the NF approach often has good performance for nonleptonic B and D weak decays induced by external W emission interactions. Using the NF approximation as a working hypothesis, it is usually assumed that the final state interactions and annihilation contributions may be disregarded. The product of quark currents in effective tetra-quark operators could be replaced by the product of the corresponding hadronic currents formed by the physical hadrons involved. HMEs of tetra-quark operators are separated into two HMEs of hadronic currents, which are further parameterized by hadronic decay constants and transition form factors.Due to their relatively large mass, vector mesons decay dominantly through strong and/or electromagnetic interactions within the SM. The weak decays of the vector mesons are rare processes and usually draw less attention. Following the conventions in Refs. [45–49], the hadronic decay constants and transition form factors are defined by HMEs of color-singlet diquark currents.
jVμ=ˉq1γμq2,
(18) jAμ=ˉq1γμγ5q2,
(19) ⟨P(k)|jVμ|0⟩=0,
(20) ⟨P(k)|jAμ|0⟩=−ifPkμ,
(21) ⟨V(k,ϵ)|jVμ|0⟩=fVmVϵμ,
(22) ⟨V(k,ϵ)|jAμ|0⟩=0,
(23) ⟨P(p2)|jVμ|B∗c(ϵ,p1)⟩=εμναβϵνPαqβVB∗c→P(q2)mB∗c+mP,
(24) ⟨P(p2)|jAμ|B∗c(ϵ,p1)⟩=+i2mB∗cϵ⋅qq2qμAB∗c→P0(q2)+iϵμ(mB∗c+mP)AB∗c→P1(q2)+iϵ⋅qmB∗c+mPPμAB∗c→P2(q2)−i2mB∗cϵ⋅qq2qμAB∗c→P3(q2),
(25) ⟨V(ϵ2,p2)|jVμ|B∗c(ϵ1,p1)⟩=−(ϵ1⋅ϵ∗2){PμVB∗c→V1(q2)−qμVB∗c→V2(q2)}+(ϵ1⋅q)(ϵ∗2⋅q)m2B∗c−m2V{[Pμ−m2B∗c−m2Vq2qμ]VB∗c→V3(q2)+m2B∗c−m2Vq2qμVB∗c→V4(q2)}−(ϵ1⋅q)ϵ∗2,μVB∗c→V5(q2)+(ϵ∗2⋅q)ϵ1,μVB∗c→V6(q2),
(26) ⟨V(ϵ2,p2)|jAμ|B∗c(ϵ1,p1)⟩=−iεμναβϵα1ϵ∗β2{[Pν−m2B∗c−m2Vq2qν]AB∗c→V1(q2)+m2B∗c−m2Vq2qνAB∗c→V2(q2)}−iεμναβPαqβm2B∗c−m2V{(ϵ∗2⋅q)ϵν1AB∗c→V3(q2)−(ϵ1⋅q)ϵ∗,ν2AB∗c→V4(q2)},
(27) where
mP (mV ) andfP (fV ) are the mass and decay constant of final pseudoscalar (vector) mesons, respectively;ϵi is the polarization vector;Ai andVi are mesonic transition form factors; the momentum P= p1 + p2 and q= p1 − p2 . There are some relationships among form factors,AB∗c→P0(0)=AB∗c→P3(0),
(28) AB∗c→V1(0)=AB∗c→V2(0),
(29) VB∗c→V3(0)=VB∗c→V4(0).
(30) The decay constants and form factors are nonperturbative hadronic parameters. These parameters are universal, and can be obtained from data and some nonperturbative methods. The experimental data on the
B∗c meson are still unavailable. Phenomenologically, form factors at the poleq2 = 0 are expressed as the overlap integrals of mesonic wave functions [45], where the mesonic wave functions are nonperturbative but process-independent physical quantities. For example, the form factors for theB∗c → Bd,s transitions have been calculated with the Wirbel-Stech-Bauer model in Ref. [20]. Recently, the form factors for theB∗c → J/ψ ,B∗d,s transitions have also been investigated with the light front quark model (LFQM) in Ref. [48]. Additionally, as noted, the latest decay constants obtained with LFQM [50, 51] are generally consistent with the data as well as those from lattice QCD simulations and QCD sum rule approaches. Following the calculation in Ref. [48], we obtain the form factors for theB∗c → ψ(2S) ,ηc(1S) ,ηc(2S) ,Bd,s transitions with LFQM. The numerical values of the form factors are listed in Table 1. The theoretical uncertainties are not particularly important or worthy of our attention for the moment. After all, the magnitude order estimation of branching ratios for the nonleptonicB∗c weak decays is quite sufficient for the present purposes.fπ = 130.2±1.2 MeVfK = 155.7±0.3 MeVfρ = 207.7±1.6 MeVfK∗ = 202.5+6.5−6.7 MeVVB∗c→ψ(1S)1 = 0.56 VB∗c→ψ(1S)2 = 0.33 VB∗c→ψ(1S)3,4 = 0.20 VB∗c→ψ(1S)5 = 1.17 VB∗c→ψ(1S)6 = 0.65 AB∗c→ψ(1S)1,2 = 0.54 AB∗c→ψ(1S)3 = 0.13 AB∗c→ψ(1S)4 = 0.14 VB∗c→ψ(2S)1 = 0.39 VB∗c→ψ(2S)2 = 0.32 VB∗c→ψ(2S)3,4 = 0.12 VB∗c→ψ(2S)5 = 0.79 VB∗c→ψ(2S)6 = 0.48 AB∗c→ψ(2S)1,2 = 0.37 AB∗c→ψ(2S)3 = 0.13 AB∗c→ψ(2S)4 = 0.07 VB∗c→B∗s1 = 0.63 VB∗c→B∗s2 = 1.06 VB∗c→B∗s3,4 = 0.40 VB∗c→B∗s5 = 3.52 VB∗c→B∗s6 = 3.02 AB∗c→B∗s1,2 = 0.53 AB∗c→B∗s3 = 0.73 AB∗c→B∗s4 = 0.85 VB∗c→B∗1 = 0.52 VB∗c→B∗2 = 1.18 VB∗c→B∗3.4 = 0.40 VB∗c→B∗5 = 3.15 VB∗c→B∗6 = 2.66 AB∗c→B∗1,2 = 0.43 AB∗c→B∗3 = 0.81 AB∗c→B∗4 = 0.89 VB∗c→ηc(1S) = 0.91 AB∗c→ηc(1S)0 = 0.66 AB∗c→ηc(1S)1 = 0.69 AB∗c→ηc(1S)2 = 0.59 VB∗c→ηc(2S) = 0.59 AB∗c→ηc(2S)0 = 0.43 AB∗c→ηc(2S)1 = 0.41 AB∗c→ηc(2S)2 = 0.51 VB∗c→Bs = 3.40 AB∗c→Bs0 = 0.69 AB∗c→Bs1 = 0.75 AB∗c→Bs2 = 0.96 VB∗c→B = 3.08 AB∗c→B0 = 0.60 AB∗c→B1 = 0.65 AB∗c→B2 = 0.91 -
The branching ratio of nonleptonic
B∗c decays is defined asBr=pcm24πm2B∗cΓB∗c¯|A|2,
(31) where
pcm is the center-of-mass momentum of final states in the rest frame of theB∗c meson.A denotes the decay amplitudes, which are collected in the Appendix. With the input parameters in Table 1 and Table 2, we obtain the branching ratios, which are listed in Table 3.mπ = 139.57 mρ = 775.26±0.25 mBd = 5279.65±0.12 mB∗d = 5324.70±0.21 mK = 493.677±0.016 mK∗ = 891.66±0.26 mηc(1S) = 2983.9±0.5 mψ(1S) = 3096.9 mBs = 5366.88±0.14 mB∗s = 5415.4+1.8−1.5 mηc(2S) = 3637.5±1.1 mψ(2S) = 3686.10±0.06 Table 2. Mesonic mass (in units of MeV) [2], where their central values are regarded as the default inputs unless otherwise specified.
decay mode CKM factor branching ratio event numbers unit here previous CEPC FCC-ee LHCb Bsπ+ V∗csVud ∼ O(1) 10−7 4.4 4.0 [20]7.3 [20]9.8 [20]44 442 13256 Bsρ+ V∗csVud ∼ O(1) 10−6 1.9 0.6 [20]1.2 [20]1.7 [20]193 1930 57887 B∗sπ+ V∗csVud ∼ O(1) 10−6 1.4 143 1428 42848 BsK+ V∗csVus ∼ O(λ) 10−8 2.2 2.0 [20]3.6 [20]4.9 [20]2 22 658 BsK∗+ V∗csVus ∼ O(λ) 10−8 7.2 2.2 [20]4.1 [20]6.0 [20]7 72 2169 B∗sK+ V∗csVus ∼ O(λ) 10−8 8.4 8 84 2533 Bdπ+ V∗cdVud ∼ O(λ) 10−8 2.1 2.1 [20]4.5 [20]6.5 [20]2 21 640 Bdρ+ V∗cdVud ∼ O(λ) 10−8 9.7 3.5 [20]7.5 [20]11.5 [20]10 97 2924 B∗dπ+ V∗cdVud ∼ O(λ) 10−8 6.3 6 63 1884 BdK+ V∗cdVus ∼ O(λ2) 10−9 1.1 1.2 [20]2.4 [20]3.5 [20]0 1 34 BdK∗+ V∗cdVus ∼ O(λ2) 10−9 4.4 1.5 [20]3.1 [20]4.8 [20]0 4 133 B∗dK+ V∗cdVus ∼ O(λ2) 10−9 3.9 0 4 117 ηc(1S)π+ V∗cbVud ∼ O(λ2) 10−8 1.5 2.2 [19]2 15 455 ηc(1S)ρ+ V∗cbVud ∼ O(λ2) 10−8 4.3 3.0 [21]2.5 [21]4 43 1300 ψ(1S)π+ V∗cbVud ∼ O(λ2) 10−8 5.1 9.2 [19]2.4 [22]5 51 1519 ηc(2S)π+ V∗cbVud ∼ O(λ2) 10−9 4.1 2.4 [19]0 4 123 ηc(2S)ρ+ V∗cbVud ∼ O(λ2) 10−8 1.2 1 12 351 ψ(2S)π+ V∗cbVud ∼ O(λ2) 10−8 1.6 3.2 [19]2 16 471 ηc(1S)K+ V∗cbVus ∼ O(λ3) 10−9 1.1 1.7 [19]0 1 33 ηc(1S)K∗+ V∗cbVus ∼ O(λ3) 10−9 2.3 1.7 [21]1.4 [21]0 2 68 ψ(1S)K+ V∗cbVus ∼ O(λ3) 10−9 3.8 7.3 [19]0 4 113 ηc(2S)K+ V∗cbVus ∼ O(λ3) 10−10 2.9 3.4 [19]0 0 9 ηc(2S)K∗+ V∗cbVus ∼ O(λ3) 10−10 6.1 0 1 18 ψ(2S)K+ V∗cbVus ∼ O(λ3) 10−9 1.2 2.4 [19]0 1 35 Table 3. Branching ratios and event numbers of nonleptonic
B∗c weak decays, assuming that about108 ,109 and3×1010 B∗c mesons will be available in the future CEPC, FCC-ee and LHCb@HL-LHC experiments, respectively;λ ≈ 0.2 is the phenomenological Wolfenstein parameter.Our comments on branching ratios are as follows.
(1) For the nonleptonic
B∗c weak decays induced by the external W emission interactions, branching ratios using the NF approach are generally of the same order of magnitude as previous estimations with the QCD-inspired QCDF and pQCD approach [19–21]. The consistency of the results between different approaches indicates that our results in this paper might be reasonable. In addition, the branching ratios for theB∗c → B∗s,dπ ,B∗s,dK ,ηc(2S)ρ ,ηc(2S)K∗ decays are estimated for the first time. It is found that theB∗c → B∗sπ decay has a relatively large branching ratio among the two-mesonB∗c weak decays. In addition, both the initialB∗c meson and one of the final mesons of the processes in question contain one and/or two of the heavy quarks. Due to the fact that the light quarks and gluon clouds are almost blind to the spin of the heavy quark with the approximation of the heavy quark limit, the heavy quark spin symmetry should be expected to relate the initial and final mesons [52]. In the extreme non-relativistic limit, it is expected to have [52]r=Br(1−→1−π)Br(1−→0−π) ≈ dBr(1−→1−ℓν)dq2|q2=m2πBr(1−→0−ℓν)dq2|q2=m2π ≈ 3.
(32) The mass difference in the phase factors furnishes the above relations in Eq. (32). Here, some relative branching ratios are
r1=Br(B∗c→B∗sπ)Br(B∗c→Bsπ)=3.2,
(33) r2=Br(B∗c→B∗dπ)Br(B∗c→Bdπ)=2.9,
(34) r3=Br(B∗c→B∗sK)Br(B∗c→BsK)=3.8,
(35) r4=Br(B∗c→B∗dK)Br(B∗c→BdK)=3.4,
(36) r5=Br(B∗c→ψ(1S)π)Br(B∗c→ηc(1S)π)=3.3,
(37) r6=Br(B∗c→ψ(1S)K)Br(B∗c→ηc(1S)K)=3.4,
(38) r7=Br(B∗c→ψ(2S)π)Br(B∗c→ηc(2S)π)=3.8,
(39) r8=Br(B∗c→ψ(2S)K)Br(B∗c→ηc(2S)K)=3.9.
(40) For the charm quark decay, the ratio
r1 of the Cabibbo-favoredB∗c → B∗sπ ,Bsπ decays is generally consistent with both the expectation Eq. (32) from the heavy quark spin symmetry and the ratio of semileptonicJ/ψ weak decaysBr(J/ψ→D∗se+νe)/Br(J/ψ→Dse+νe) ≈ 3.1 obtained with the QCD sum rules [49].(2) For the
B∗c decays induced either by the bottom quark decay [type (1)] or by the charm quark decay [type (2)], there are some clear hierarchical relationships among branching ratios according to the CKM factors of decay amplitudes. The Cabibbo-favoredB∗c → B∗sπ andBsρ decays have relatively maximum branching ratios which can reach up toO(10−6) . The CKM-suppressedB∗c → ηc(1S,2S)K ,ψ(1S,2S)K , andηc(1S,2S)K∗ decays, whose amplitudes are proportional to Wolfenstein parameterλ3 , have relatively minimal branching ratios. In addition, there are three partial wave amplitudes for final states containing one pseudoscalar plus vector mesons, while there is only the p-wave amplitude for two final pseudoscalar mesons. Hence, there are some relationships,Br(B∗c→B∗sπ) > Br(B∗c→Bsπ) andBr(B∗c→Bsρ) > Br(B∗c→Bsπ) , and these hierarchical relationships are also true for other similar cases where final states have the same valence quark compositions but different orbit-spin couplings between quarks.(3) A more careful theoretical investigation of these decays in question is necessary and helpful to explore the hadronic
B∗c weak decays experimentally. It should be pointed out that many influences can affect the final numerical results on nonleptonicB∗c weak decays, such as the nonfactorizable contributions. It is clear from Eq. (9)–(13) that|Vcb| has the largest uncertainty, about3% , among the CKM elements involved. Except for the decay constant offK∗ , the uncertainties from other decay constants in Table 1 is less than1% . For the nonleptonic decays induced by the external W boson emission interactions, the nonfactorizable contributions to the coefficienta1 is about15% for charm quark decays [47], i.e., theB∗c → B(∗)d,s transitions, and about5% for bottom quark decays [35, 36], i.e., theB∗c → ψ andηc transitions. Large theoretical uncertainties come mainly from the hadronic transition form factors, and some of them can reach about40% [48]. In addition, the undetermined decay width of theB∗c meson (for example,Γ(B∗c→γBc) = 20 − 135 eV, see Table 2 of Ref. [21]) will bring very large theoretical uncertainties to branching ratios. An accurate theoretical prediction seems to be temporarily unavailable. All numbers in Table 3 are rough estimates and only indicative of the expected order of magnitude. Here, what matters to us is whether it is possible to study the unacquaintedB∗c mesons in future experiments. Hence, a rough estimate is sufficient.(4) As it is well known that the vector
ρ andK∗ mesons are resonances, they will decay immediately via the strong interaction, with branching ratiosBr(ρ→ππ) ∼ 100% andBr(K∗→Kπ) ∼ 100% [2]. The vectorρ (K∗ ) meson is reconstructed experimentally by the final pseudoscalar mesons. An educated guess is that the branching ratios for the three-body decay modes,B∗c → Bs,dππ ,Bs,dπK ,ηcππ ,ηcπK , should be of a similar order of magnitude to those for theB∗c → Bs,dρ ,Bs,dK∗ ,ηcρ ,ηcK∗ , decays, respectively. If the contributions from other possible resonances reconstructed from theππ ,πK ,Bs,dπ ,Bs,dK ,ηcπ andηcK states are taken into consideration, the above hadronic three-bodyB∗c decays are likely to have even larger branching ratios. All in all, it is not utopian to expect that the Cabibbo-favoredB∗c → Bsπ ,B∗sπ ,Bsρ decays, the singly-Cabibbo-suppressedB∗c → BsK ,B∗sK ,BsK∗ ,Bdπ ,B∗dπ ,Bdρ decays, and even the CKM-suppressedB∗c → ηc(1S)π ,ηc(1S,2S)ρ ,ψ(1S,2S)π decays might be observable in the future CEPC, FCC-ee and LHCb@HL-LHC experiments. -
It has been established that the charged ground vector
B∗c meson carrying explicit flavor numbers should exist according to the quark model, but to date this has been merely on the theoretical calculation level rather than the experimental measurement level. The identification of theB∗c meson in experiments is necessary and of important significance to the quark model and SM. The signal reconstruction of theB∗c meson from the electromagnetic decayB∗c → Bcγ is severely hindered by electromagnetic background pollution. Inspired by the prospects of huge numbers ofB∗c mesons in future high-energy and high-luminosity colliders, an attractive and competitive choice might be searching for theB∗c meson from its nonleptonic weak decays, where the charged final hadrons are comparatively easily and effectively identified in experiments. In this paper, we study two kinds of nonleptonicB∗c meson weak decays resulting from external W boson emission interactions by using the factorization approximation and form factors from the light front quark model, one originating from bottom quark decays and the other from charm quark decays. It is found that the branching ratios for the Cabibbo-favoredB∗c → B∗sπ ,Bsρ decays can reach up toO(10−6) , and have the first priority to be studied in experiments. For the singly-Cabibbo-suppressedB∗c → BsK ,B∗sK ,BsK∗ ,Bdπ ,B∗dπ ,Bdρ decays and the CKM-suppressedB∗c → ηc(1S,2S)π ,ηc(1S,2S)ρ ,ψ(1S,2S)π decays, several hundred or even thousands of events might be observable in CEPC, FCC-ee and LHCb@HL-LHC experiments. This study provides a ready and helpful reference for experimental discovery and investigation ofB∗c mesons in the future. -
With the conventions of Ref. [46], the amplitudes for nonleptonic
B∗c weak decays can be written as follows.Based on the conservation of angular momentum, there are only p-wave amplitudes contributing to
B∗c meson decay into two pseudoscalar mesons.A(B∗+c→B0sπ+)=V∗csVudfπMBsp(ϵB∗c⋅pBs),
{\cal A}(B_{c}^{{\ast}+}{\to}B_{s}^{0}K^+) = V_{cs}^{\ast}\,V_{us}\, f_{K}\, {\cal M}_{p}^{B_{s}}\, ({\epsilon}_{B_{c}^{\ast}}{\cdot}p_{B_{s}}) , \tag{A2}
{\cal A}(B_{c}^{{\ast}+}{\to}B_{d}^{0}{\pi}^+) = V_{cd}^{\ast}\,V_{ud}\, f_{\pi}\, {\cal M}_{p}^{B}\, ({\epsilon}_{B_{c}^{\ast}}{\cdot}p_{B_{d}}) , \tag{A3}
{\cal A}(B_{c}^{{\ast}+}{\to}B_{d}^{0}K^+) = V_{cd}^{\ast}\,V_{us}\, f_{K}\, {\cal M}_{p}^{B}\, ({\epsilon}_{B_{c}^{\ast}}{\cdot}p_{B_{d}}) , \tag{A4}
{\cal A}(B_{c}^{{\ast}+}{\to}{\eta}_{c}{\pi}^+) = V_{cb}^{\ast}\,V_{ud}\, f_{\pi}\, {\cal M}_{p}^{{\eta}_{c}}\, ({\epsilon}_{B_{c}^{\ast}}{\cdot}p_{{\eta}_{c}}) , \tag{A5}
{\cal A}(B_{c}^{{\ast}+}{\to}{\eta}_{c}K^+) = V_{cd}^{\ast}\,V_{us}\, f_{K}\, {\cal M}_{p}^{{\eta}_{c}}\, ({\epsilon}_{B_{c}^{\ast}}{\cdot}p_{{\eta}_{c}}) , \tag{A6}
with the common factor of p-wave partial amplitudes,
{\cal M}_{p}^{B_{s}} = \sqrt{2}\, G_{\rm F}\, a_{1}\, m_{B_{c}^{\ast}}\, A_{0}^{B_{c}^{\ast}{\to}B_{s}} , \tag{A7}
{\cal M}_{p}^{B} = \sqrt{2}\, G_{\rm F}\, a_{1} \, m_{B_{c}^{\ast}}\, A_{0}^{B_{c}^{\ast}{\to}B} , \tag{A8}
{\cal M}_{p}^{{\eta}_{c}} = \sqrt{2}\, G_{\rm F}\, a_{1} \, m_{B_{c}^{\ast}}\, A_{0}^{B_{c}^{\ast}{\to}{\eta}_{c}} , \tag{A9}
where coefficient
a_{1} = C_{1} + C_{2}/N_{c} is generally influenced by nonfactorizable contributions and can be regarded as a phenomenological parameter, especially for charm quark decays. The value ofa_{1} {\approx} 1.2 is used in our calculation.There are three partial wave amplitudes contributing to
B_{c}^{\ast} meson decay into one pseudoscalar meson plus one vector meson. The general decay amplitude is written as\begin{aligned}[b] {\cal A}(B_{c}^{\ast}{\to}VP) =& {\cal M}_{s}\, ( {\epsilon}_{B_{c}^{\ast}}{\cdot}{\epsilon}_{V}^{\ast} ) \\ &+ \frac{ {\cal M}_{d} }{ m_{V_{B_{c}^{\ast}}}\,m_{V} }\, ( {\epsilon}_{B_{c}^{\ast}}{\cdot}p_{V} ) \, ( {\epsilon}_{V}^{\ast}{\cdot}p_{B_{c}^{\ast}} ) \\ & + \frac{ {\cal M}_{p} }{ m_{B_{c}^{\ast}}\,m_{V} }\, {\varepsilon}_{{\mu}{\nu}{\alpha}{\beta}}\, {\epsilon}_{B_{c}^{\ast}}^{\mu}\, {\epsilon}_{V}^{{\ast}{\nu}}\, p_{B_{c}^{\ast}}^{\alpha}\, p_{V}^{\beta} , \end{aligned}\tag{A10}
where
{\cal M}_{s,p,d} correspond to the s-, p-, and d-wave partial amplitudes.For
B_{c}^{{\ast}+} {\to} B_{s}^{{\ast}0}{\pi}^{+} ,B_{s}^{{\ast}0}K^{+} decays, one has{\cal M}_{s}^{ B_{s}^{\ast}{\pi} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, \big( m_{B_{c}^{\ast}}^{2}-m_{B_{s}^{\ast}}^{2} \big)\, V_{1}^{ B_{c}^{\ast} {\to} B_{s}^{\ast} } , \tag{A11}
\begin{aligned}[b] {\cal M}_{d}^{ B_{s}^{\ast}{\pi} } =& -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{B_{s}^{\ast}}\, \big( V_{4}^{ B_{c}^{\ast} {\to} B_{s}^{\ast} } \\& - V_{5}^{ B_{c}^{\ast} {\to} B_{s}^{\ast} } + V_{6}^{ B_{c}^{\ast} {\to} B_{s}^{\ast} } \big) , \end{aligned}\tag{A12}
{\cal M}_{p}^{ B_{s}^{\ast}{\pi} } = -2\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{B_{s}^{\ast}}\, A_{1}^{ B_{c}^{\ast} {\to} B_{s}^{\ast} } , \tag{A13}
{\cal M}_{i}^{ B_{s}^{\ast}K } = {\cal M}_{i}^{ B_{s}^{\ast}{\pi} } \big( V_{ud}{\to}V_{us},\,f_{\pi}{\to}f_{K} \big), \quad {\rm{for }}~~ i=s,p,d . \tag{A14}
For
B_{c}^{{\ast}+} {\to} B_{d}^{{\ast}0}{\pi}^{+} ,B_{d}^{{\ast}0}K^{+} decays, one has{\cal M}_{s}^{ B_{d}^{\ast}{\pi} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, \big( m_{B_{c}^{\ast}}^{2}-m_{B_{d}^{\ast}}^{2} \big)\, V_{1}^{ B_{c}^{\ast} {\to} B_{d}^{\ast} } , \tag{A15}
\begin{aligned}[b] {\cal M}_{d}^{ B_{d}^{\ast}{\pi} }=& -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{B_{d}^{\ast}}\, \big( V_{4}^{ B_{c}^{\ast} {\to} B_{d}^{\ast} } \\& - V_{5}^{ B_{c}^{\ast} {\to} B_{d}^{\ast} } + V_{6}^{ B_{c}^{\ast} {\to} B_{d}^{\ast} } \big) , \end{aligned}\tag{A16}
{\cal M}_{p}^{ B_{d}^{\ast}{\pi} } = -2\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{B_{d}^{\ast}}\, A_{1}^{ B_{c}^{\ast} {\to} B_{d}^{\ast} } , \tag{A17}
{\cal M}_{i}^{ B_{d}^{\ast}K } = {\cal M}_{i}^{ B_{d}^{\ast}{\pi} } \big( V_{ud}{\to}V_{us},\,f_{\pi}{\to}f_{K} \big), \quad {\rm{for }}~~ i=s,p,d . \tag{A18}
For
B_{c}^{{\ast}+} {\to} {\psi}{\pi}^{+} ,{\psi}K^{+} decays, one has{\cal M}_{s}^{ {\psi}{\pi} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, \big( m_{B_{c}^{\ast}}^{2}-m_{\psi}^{2} \big)\, V_{1}^{ B_{c}^{\ast} {\to} {\psi} } ,\tag{A19}
\begin{aligned}[b] {\cal M}_{d}^{ {\psi}{\pi} } =& -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{\psi}\, \big( V_{4}^{ B_{c}^{\ast} {\to} {\psi} } \\& - V_{5}^{ B_{c}^{\ast} {\to} {\psi} } + V_{6}^{ B_{c}^{\ast} {\to} {\psi} } \big) , \end{aligned}\tag{A20}
{\cal M}_{p}^{ {\psi}{\pi} } = -2\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{ \pi }\, a_{1}\, m_{B_{c}^{\ast}}\, m_{\psi}\, A_{1}^{ B_{c}^{\ast} {\to} {\psi} } , \tag{A21}
{\cal M}_{i}^{ {\psi}K } = {\cal M}_{i}^{ {\psi}{\pi} } \big( V_{ud}{\to}V_{us},\,f_{\pi}{\to}f_{K} \big), \quad {\rm{for }}~~ i=s,p,d . \tag{A22}
For
B_{c}^{{\ast}+} {\to} B_{s}^{0}{\rho}^{+} ,B_{s}^{0}K^{{\ast}+} decays, one has{\cal M}_{s}^{ B_{s}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \big( m_{B_{c}^{\ast}}+m_{B_{s}} \big)\, A_{1}^{ B_{c}^{\ast} {\to} B_{s} } ,\tag{A23}
{\cal M}_{d}^{ B_{s}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{B_{s}} }\, A_{2}^{ B_{c}^{\ast} {\to} B_{s} } , \tag{A24}
{\cal M}_{p}^{ B_{s}{\rho} } = -\;\,\frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cs}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{B_{s}} }\, V^{ B_{c}^{\ast} {\to} B_{s} } , \tag{A25}
\begin{equation} {\cal M}_{i}^{ B_{s}K^{\ast} } = {\cal M}_{i}^{ B_{s}{\rho} } \big( V_{ud}{\to}V_{us}, \, f_{\rho}{\to}f_{K^{\ast}}, \, m_{\rho}{\to}m_{K^{\ast}} \big), \quad {\rm{for }}~~ i=s,p,d . \end{equation} \tag{A26}
For
B_{c}^{{\ast}+} {\to} B_{d}^{0}{\rho}^{+} ,B_{d}^{0}K^{{\ast}+} decays, one has{\cal M}_{s}^{ B_{d}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \big( m_{B_{c}^{\ast}}+m_{B_{d}} \big)\, A_{1}^{ B_{c}^{\ast} {\to} B_{d} } , \tag{A27}
{\cal M}_{d}^{ B_{d}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{B_{d}} }\, A_{2}^{ B_{c}^{\ast} {\to} B_{d} } , \tag{A28}
{\cal M}_{p}^{ B_{d}{\rho} } = -\;\,\frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cd}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{B_{d}} }\, V^{ B_{c}^{\ast} {\to} B_{d} } , \tag{A29}
\begin{equation} {\cal M}_{i}^{ B_{d}K^{\ast} } = {\cal M}_{i}^{ B_{d}{\rho} } \big( V_{ud}{\to}V_{us}, \, f_{\rho}{\to}f_{K^{\ast}}, \, m_{\rho}{\to}m_{K^{\ast}} \big), \quad {\rm{for }}~~ i=s,p,d . \end{equation} \tag{A30}
For
B_{c}^{{\ast}+} {\to} {\eta}_{c}{\rho}^{+} ,{\eta}_{c}K^{{\ast}+} decays, one has{\cal M}_{s}^{ {\eta}_{c}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \big( m_{B_{c}^{\ast}}+m_{{\eta}_{c}} \big)\, A_{1}^{ B_{c}^{\ast} {\to} {\eta}_{c} } , \tag{A31}
{\cal M}_{d}^{ {\eta}_{c}{\rho} } = -{\rm i}\, \frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{{\eta}_{c}} }\, A_{2}^{ B_{c}^{\ast} {\to} {\eta}_{c} } , \tag{A32}
{\cal M}_{p}^{ {\eta}_{c}{\rho} } = -\;\,\frac{ G_{\rm F} }{ \sqrt{2} }\, V_{cb}^{\ast}\,V_{ud}\, f_{\rho}\, m_{\rho}\, a_{1}\, \frac{ 2\,m_{B_{c}^{\ast}}\,m_{\rho} } { m_{B_{c}^{\ast}}+m_{{\eta}_{c}} }\, V^{ B_{c}^{\ast} {\to} {\eta}_{c} } , \tag{A33}
\begin{equation} {\cal M}_{i}^{ {\eta}_{c}K^{\ast} } = {\cal M}_{i}^{ {\eta}_{c}{\rho} } \big( V_{ud}{\to}V_{us}, \, f_{\rho}{\to}f_{K^{\ast}}, \, m_{\rho}{\to}m_{K^{\ast}} \big), \quad {\rm{for }}~~ i=s,p,d . \end{equation}\tag{A34}
Possibility of experimental study on nonleptonic B∗c weak decays
- Received Date: 2022-09-23
- Available Online: 2023-01-15
Abstract: The ground vector