Purely leptonic decays of the ground charged vector mesons
Abstract
The study of the purely leptonic decays of the ground charged vector mesons is very interesting and significant in determining the CKM matrix elements, obtaining the decay constant of vector mesons, examining the lepton flavor universality, and searching for new physics beyond the standard model. These purely leptonic decays of the ground charged vector mesons are induced by the weak interactions within the standard model, and usually have very small branching ratios, , , , , and . Inspired by the potential prospects of LHCb, Belle-II, STCF, CEPC and FCC-ee experiments, we discussed the probabilities of experimental investigation on these purely leptonic decays. It is found that the measurements of these decays might be possible and feasible with the improvement of data statistics, analytical technique, and measurement precision in the future. (1) With the hadron-hadron collisions, the purely leptonic decays of , , and mesons might be accessible at LHC experiments. (2) With the collisions, the purely leptonic decays of and mesons might be measurable with over bosons available at CEPC and FCC-ee experiments. In addition, the decays could also be studied at Belle-II and SCTF experiments.
Published : Eur. Phys. J. C 81, 1110 (2021).
I Introduction
In the quark model pdg2020 ; pl.8.214 ; zweig , mesons are generally regarded as bound states of the valence quark and antiquark . The classifications of mesons are usually based on the spin-parity quantum number of the system. The spin of meson is given by the relation . The orbital angular momentum and total spin of the system are respectively and , where for antiparallel quark spins, and for parallel quark spins. By convention, quarks have a positive parity and antiquarks have a negative parity. Hence, the parity of meson is . The states are the ground-state pseudoscalars with and vectors with . Both quarks and leptons are fermions with spin . Mesons are composed of a pair of fermions — quark and antiquark, therefore, they could in principle decay into a pair of fermions, for example, lepton and antilepton. The experimental observation of the two-body purely leptonic decays of mesons could be a clear and characteristic manifestation of the quark model. These leptonic decays provide us with valuable opportunities to fully investigate the microstructure and properties of mesons. The study of two-body purely leptonic decays of mesons is very interesting and significant.
The valence quarks of the electrically charged mesons must have different flavors. Within the standard model (SM) of elementary particles, the purely leptonic decays of the charged mesons (PLDCM) are typically induced by the tree-level exchange of the gauge bosons , the quanta of the weak interaction fields. Up to today, the masses of all the experimentally observed mesons are much less than those of bosons. Consequently, the massive bosons are virtual propagators rather than physical particles in the true picture of PLDCM. Phenomenologically, by integrating out the contributions from heavy dynamical degrees of freedom such as the fields, PLDCM can be properly described by the low-energy effective theory in analogy with the Fermi theory for decays. Considering the fact that leptons are free from the strong interactions, the corresponding effective Hamiltonian RevModPhys.68.1125 for PLDCM could be written as the product of quark currents and leptonic currents,
(1) |
where the contributions of the bosons are embodied in the Fermi coupling constant pdg2020 , and is the Cabibbo-Kobayashi-Maskawa (CKM) Cabibbo ; Kobayashi matrix element between the quarks in the charged mesons. The decay amplitudes can be written as,
(2) |
The leptonic part of amplitudes can be calculated reliably with the perturbative theory. The hadronic matrix elements (HMEs) interpolating the diquark currents between the mesons concerned and the vacuum can be parameterized by the decay constants.
With the conventions of Refs. JHEP.0605.004 ; JHEP.0703.069 , the HMEs of diquark currents are defined as,
(3) | |||||
(4) | |||||
(5) | |||||
(6) |
where the nonperturbative parameters of and are the decay constants of pseudoscalar and vector mesons, respectively; and and are the mass and polarization vector, respectively. To the lowest order, the decay widths are written as,
(7) | |||||
(8) |
where and are the masses of the charged pseudoscalar meson and lepton, respectively.
It is clearly seen from the above formula that the highly precise measurements of PLDCM will allow the relatively accurate determinations of the product of the decay constants and CKM elements, . Theoretically, the decay constants are nonperturbative parameters, and they are closely related with the wave functions at the origin which cannot be computed from first principles. There still exist some discrepancies among theoretical results of the decay constants with different methods, such as the potential model, QCD sum rules, lattice QCD, and so on. If the magnitudes of CKM element are fixed to the values of Ref. pdg2020 , the decay constants will be experimentally measured, and be used to seriously examine the different calculations on the decay constants with various theoretical models. Likewise, if the decay constants are well known to sufficient precision, the magnitudes of the corresponding CKM element will be experimentally determined, and provide complementary information to those from other processes. Within SM, the and decays are induecd by the axial-vector current of Eq.(4) and vector current of Eq.(5), respectively; and the electroweak interactions assign the vector-minus-axial-vector () currents to the bosons. The CKM elements determined from two different and complementary parts of the electroweak interactions, charged vector and axial-vector currents, could be independently examined. The latest CKM elements determined by PLDCM, such as , and , differ somewhat from those by exclusive and inclusive semileptonic meson decays pdg2020 . The CKM elements extracted from various processes can be combined to test the electroweak characteristic charged-current interactions.
Within SM, the lepton-gauge-boson electroweak gauge couplings are generally believed to be universal and process independent, which is called lepton flavor universality (LFU). However, there are some hints of LFU discrepancies between SM predictions and experimental measurements, such as the ratios of branching fractions of semileptonic decays with pdg2020 . The LFU validity can be carefully investigated through the PLDCM processes. Beyond SM, some possible new heavy particles accompanied with novel interactions, such as the charged higgs bosons, would affect PLDCM and LFU, and might lead to detectable effects. So PLDCM provide good arenas to search for the smoking gun of new physics (NP) beyond SM.
By considering the angular momentum conservation and the final states including a left-handed neutrino or right-handed antineutrino, the purely leptonic decay width of charged pseudoscalar meson, Eq.(7), is proportional to the square of the lepton mass. This is called helicity suppression. While there is no helicity suppression for the purely leptonic decay of charged vector meson (PLDCV). From the analytical expressions of Eq.(7) and Eq.(8), the decay width of pseudoscalar meson is suppressed by the factor compared with that of vector meson. What’s more, both the masses and the decay constants of vector mesons are relatively larger than those of corresponding pseudoscalar mesons, which would result in an enhancement of the decay widths for vector mesons. Of course, the vector mesons decay dominantly through the strong and/or electromagnetic interactions. The branching ratios for the PLDCV weak decays are usually very small, sometimes even close to the accessible limits of the existing and the coming experiments.
Inspired by the potential prospects of the future high-intensity and high-energy frontiers, along with the noticeable increase of experimental data statistics, the remarkable improvement of analytical technique and the continuous enhancement of measurement precision, the carefully experimental study of PLDCV might be possible and feasible. In this paper, we will focus on the PLDCV within SM to just provide a ready reference. The review of the purely leptonic decays of charged pseudoscalar mesons can be found in Ref. pdg2020 .
II decays
The mass of the meson, MeV pdg2020 , is much larger than that of two-pion pair. The rate of the meson decay into two pions via the strong interactions is almost 100%, which results in the very short lifetime s pdg2020 . The direct measurements of the electroweak properties of the meson would definitely be very challenging. It is evident from Eq.(8) that the parameter of could be experimentally determined from the observations of decay widths for the decays (if it is not specified, the corresponding charge-conjugation processes are included in this paper), with the coupling constant , the masses of lepton and meson .
The precise values of the CKM element in ascending order of measurement accuracy mainly come from transitions between the super-allowed nuclear analog states with quantum number of both and isospin , between mirror nuclei with , between neutron and proton, between charged and neutral pions CKM2016-vud . These four results for are basically consistent with one another. The result of the super-allowed nuclear transitions has an uncertainty a factor of about 10 smaller than the other results, and thus dominates the weighted average value CKM2016-vud . The best value from super-allowed nuclear transitions is pdg2020 , which is smaller compared with the 2018 value pdg2018 , as illustrated Fig. 1. This reduction of the value of leads to a slight deviation from the first row unitarity requirement . The current precision of the CKM element is about 0.01%. The latest value from the global fit in SM, pdg2020 , will be used in our calculation.

DVP | PhysRevD.69.094013 111With light-cone wave functions and parameters of Ref. epjc.22.655 . | PhysRevD.69.094013 222With light-cone wave functions and parameters of Ref. PhysRevD.60.114023 . | PhysRevD.69.094013 333With Gaussian wave functions and parameters of Ref. epjc.22.655 . | PhysRevD.69.094013 444With Gaussian wave functions and parameters of Ref. PhysRevD.60.114023 . | PhysRevD.94.074018 |
---|---|---|---|---|---|
SR | plb.436.351 555With nonlocal condensates like functions. | plb.436.351 666With Ball-Braun wave functions PhysRevD.54.2182 . | plb.436.351 777With Chernyak-Zhitnitsky wave functions PhysRept.112.173 . | PhysRevD.58.094016 | JHEP.0604.046 |
RQM | epja.24.411 888With Gaussian spatial wave functions and adjusted parameters of Ref. epja.24.411 . | epja.24.411 999With Gaussian spatial wave functions and parameters of Ref. plb.602.212 . | epja.24.411 101010With rational spatial wave functions and adjusted parameters of Ref. epja.24.411 . | epja.24.411 111111With rational spatial wave functions and parameters of Ref. plb.602.212 . | plb.635.93 |
LFQM | PhysRevD.75.034019 121212With Coulomb plus linear potential model. | PhysRevD.75.034019 131313With Coulomb plus harmonic oscillator potential model. | PhysRevC.92.055203 121212With Coulomb plus linear potential model. | cpc.42.073102 141414With a dilation parameter GeV. | cpc.42.073102 151515With a dilation parameter GeV. |
LFQM | jpg.39.025005 | PhysRevD.98.114018 | PhysRevD.100.014026 | ||
LQCD | PhysRevD.65.054505 | PTP.119.599 | PhysRevD.80.054510 | JHEP.1704.082 | cpc.42.063102 |
other | plb.635.93 | PhysRevC.60.055214 | pan.79.444 |
The decay constant is an very important characteristics of the meson. Compared with the CKM element , the present precision of decay constant is still not very high and needs to be improved. Theoretically, the estimations from different methods are more or less different from each other and even calculations with the same method sometimes give the diverse results. Some theoretical estimations on the decay constant are presented in Table 1. Experimentally, the decay constant can be obtained from the 1-prong hadronic decay. The partial width for the decay is given by Ref. JHEP.1504.101 ,
(9) |
where the factor includes the electroweak corrections JHEP.1504.101 ; PhysRevLett.61.1815 ; PhysRevD.42.3888 . With the mass MeV and lifetime fs pdg2020 , and branching ratio % pdg2020 , one can easily extract the decay constant MeV, which agrees well with the latest numerical simulation result from lattice QCD MeV cpc.42.063102 . The more accurate decay constant will be used in our calculation.
For the decays, one can obtain the PLDCV partial decay widths with Eq.(8) and the corresponding branching ratios with the full width MeV pdg2020 ,
(10) | |||||
(11) | |||||
(12) | |||||
(13) |
where the uncertainties come from the uncertainties of mass , decay constant and CKM element , and additional decay width for branching ratios. Clearly, the branching ratios are very small. Given the identification efficiency and pollution from background, the decays might be measured only with more than events available.
There are at least three possible ways to experimentally produce the charged mesons in the electron-position collisions, (a) the prompt pair production , (b) the pair production via decay , and (c) the single production via decay . The cross section has been determined by the BaBar group to be fb near the center-of-mass energy GeV PhysRevD.78.071103 . Assuming the production cross section PhysRevC.75.065202 ; plb.763.87 , it could be speculated that fb near GeV. There would be only about pairs with a data sample of PTEP.2019.123C01 near at the Belle-II detector or a data sample of epjconf.212.01010 near with the future super-tau-charm factory like STCF or SCTF PhysRevLett.127.012003 ; STCF ; SCTF . The charge mesons can in principle be produced from the , and decays. The branching ratios are
(14) | |||||
(15) | |||||
(16) | |||||
(17) |
where the brancing ratio is assumed to be the same order of magnitude as pdg2020 from the phenomenological analysis based on the flavor- symmetry PhysRevD.32.2961 . Now, there are events at Belle epjc74.3026 , events PhysRevLett.127.012003 at BES-III, and events at KLOE/KLOE-2 KLOE2.2020 available. It is expected that only about events PTEP.2019.123C01 and events at SCTF or STCF PhysRevLett.127.012003 could be accumulated. It is clearly seen that unless a very significant enhancement to branching ratios from some NP, the experimental data on the meson are too scarce to search for the decays at the electron-position collisions in the near future, which result in the natural difficulties to understand the meson.
The production cross sections of prompt and -from-b mesons in proton–proton collisions at TeV are measured by LHCb to be b and b, respectively, JHEP.2015.10.172 . It is expected that some events could be accumulated at TeV with an integrated luminosity of at LHCb 1808.08865 . There are only about mesons from decays available for prying into the PLDCV decays. At the same time, the inclusive cross-sections for prompt charm production at LHCb at TeV are measured to be JHEP.2016.03.159 . Analogically assuming the inclusive cross section of prompt meson production at LHCb at TeV is , some events would be accumulated with an integrated luminosity of at LHCb 1808.08865 . Optimistically assuming the reconstruction efficiency is about , there would be about events of the decays at LHCb, and more events with the enhanced branching ratios from NP contributions. Even through it will be very challenging for experimental analysis due to the complex background in hadron-hadron collisions, there is still a strong presumption that the decays could be explored and studied at LHC in the future. In addition, it is expected that an integrated luminosity exceeding would be reached at the future HE-LHC experiments epjst.228.1109 . More experimental data at HE-LHC would make the study of the decays indeed feasible in hadron-hadron collisions.
III decays
The parameter product could be experimentally determined from the decays using Eq.(8). Like the meson, the mass of the meson, MeV, is above the threshold of pair, and the partial branching ratio of the meson decay into pair via the strong interactions is almost 100% pdg2020 . It is not hard to imagine that the very short lifetime s would enable the measurements of the electroweak properties of the meson to be very challenging or nearly impossible.
The CKM element up to the order of , where is a Wolfenstein parameter. The current precision of the CKM element from purely leptonic and semileptonic meson decays and hadronic decays are 0.2%, 0.3% and 0.6%, respectively. It is seen from Fig. 2 that these three results for are not very consistent with one another. So if the decays could be measured, they would provide another determination and constraint to . Probably due to the reduction of the value of , the latest value from the global fit in SM, pdg2020 , is slightly larger than the 2018 value, to satisfy the first row unitarity requirement.

RQM | epja.24.411 111With Gaussian spatial wave functions and adjusted parameters of Ref. epja.24.411 . | epja.24.411 222With Gaussian spatial wave functions and parameters of Ref. plb.602.212 . | epja.24.411 333With rational spatial wave functions and adjusted parameters of Ref. epja.24.411 . | epja.24.411 444With rational spatial wave functions and parameters of Ref. plb.602.212 . | plb.635.93 |
---|---|---|---|---|---|
LFQM | PhysRevD.75.034019 555With Coulomb plus linear potential model. | PhysRevD.75.034019 666With Coulomb plus harmonic oscillator potential model. | PhysRevC.92.055203 555With Coulomb plus linear potential model. | cpc.42.073102 777With a dilation parameter GeV. | cpc.42.073102 888With a dilation parameter GeV. |
LFQM | jpg.39.025005 | PhysRevD.98.114018 | PhysRevD.100.014026 | ||
LQCD | PhysRevD.65.054505 999The meson mass is used as input. | PhysRevD.65.054505 101010The meson mass is used as input. | PhysRevD.84.014505 | ||
other | PhysRevD.58.094016 | plb.635.93 | PhysRevC.60.055214 |
Some theoretical results on the decay constant are presented in Table 2. Like the case of the decay constant , the model dependence of theoretical estimations on the decay constant is also obvious. Experimentally, the decay constant can be obtained from the hadronic decays. Using Eq.(9) and experimental data on branching ratio % pdg2020 , one can obtain the decay constant MeV. The value of is much less than that of LQCD results, and will be used in our calculation.
For the decays, the SM expectations on the partial decay widths and branching ratios are,
(18) | |||||
(19) | |||||
(20) | |||||
(21) |
The decay width MeV pdg2020 is used in our calculation. It is apparent that more than events are the minimum requirement for experimentally studying the decays.
Based on the -spin symmetry, the production mechanism of the mesons in electron-position collisions is similar to that of the mesons. An educated guess is that the cross section fb and fb near and , respectively. The branching ratios of decays are pdg2020 ,
(22) | |||||
(23) | |||||
(24) | |||||
(25) | |||||
(26) | |||||
(27) |
It is approximately estimated that %. Hence, the experiemtal data on the mesons at the collisions, which would be available by either the prompt pair production at SuperKEKB and SCTF experiments or the production via decay at SCTF, are far from sufficient for investigating the decays. If we assume that the inclusive cross section of prompt meson production in collisions at the center-of-mass energy of TeV is similar to that of mesons, about , there would be some events to be available with an integrated luminosity of at LHCb, which correspond to about events of the decays. It should be some glimmer of hope for observation and scrutinies of the decays at hadron-hadron collisions in the future, particularly at the planning HE-HLC.
IV decays
The mass of mesons, MeV, are just above the threshold of pair. The meson decays via the strong interactions are dominant, and the ratio of branching ratios pdg2020 , , basically agrees with the relations of isospin symmetry. It should be pointed out that the strong decays are highly suppressed by the compact phase spaces because of 6 MeV. The branching ratio of the magnetic dipole transition is small, pdg2020 . Hence, the decay width of mesons is narrow, keV pdg2020 . From the decays, the parameter is expected to be experimentally determined.

NRQM | PhysRevD.71.113006 | epl.115.21002 111Without QCD radiative corrections. | epl.115.21002 222With QCD radiative corrections. | epjp.132.80 111Without QCD radiative corrections. | epjp.132.80 222With QCD radiative corrections. |
---|---|---|---|---|---|
NRQM | plb.635.93 | ahep.2018.7032041 111Without QCD radiative corrections. | ahep.2018.7032041 222With QCD radiative corrections. | mpla.17.803 | |
RQM | plb.635.93 | mpla.17.803 | PhysRevD.55.6944 333With constituent quark masses for the light quarks and . | PhysRevD.55.6944 444With current quark masses for the light quarks and . | |
LFQM | PhysRevD.75.073016 555With Coulomb plus linear potential model. | PhysRevD.75.073016 666With Coulomb plus harmonic oscillator potential model. | PhysRevD.81.114024 777With the Gaussian type wave functions. | PhysRevD.81.114024 888With the power-law type wave functions. | PhysRevD.98.114018 |
LFQM | PhysRevC.92.055203 999With Coulomb plus linear potential model. | cpc.42.073102 101010With a dilation parameter GeV. | cpc.42.073102 111111With a dilation parameter GeV. | jpg.39.025005 | PhysRevD.100.014026 |
LFQM | epjc.76.313 777With the Gaussian type wave functions. | epjc.76.313 888With the power-law type wave functions. | epjp.133.134 121212With Martin potential model plb.93.338 . | epjp.133.134 131313With Cornell potential model zpc.33.135 . | epjp.133.134 141414With logarithmic potential model plb.71.153 . |
LQCD | PhysRevD.60.074501 | npb.619.507 151515With npb.619.507 and MeV pdg2020 . | jhep.1202.042 | PhysRevD.96.034524 | cpc.45.023109 |
SR | PhysRevD.88.014015 | epjc.75.427 | plb.735.12 | ijmpa.30.1550116 | |
other | epjp.133.134 161616With harmonic plus Yukawa potential model epjp.132.80 . | PhysRevD.75.116001 | plb.633.492 | PhysRevD.58.014007 |
Currently, the precise values of the CKM element comes mainly from the leptonic and semileptonic meson decays pdg2020 , as illustrated in Fig. 3. Because of the decay width of Eq.(7) being proportional to , the decay is helicity suppressed. And the decay suffers from the complications caused by the additional neutrino in decays. The decay is the most favorable mode for experimental measurement. For the values of from the purely leptonic decay , the experimentally statistical uncertainties are dominant uncertainties. For the values of from the semileptonic meson decays, the theoretical uncertainties from the form factor controlled by nonperturbative dynamics are dominant uncertainties. It is clearly seen from Fig. 3 that the experimental uncertainties have not decreased significantly recently. Besides, can also be determined from the neutrino-induced charm production data pdg2020 , but the relevant experimental data have not been updated after the measurements given by the CHARM-II Collaboration in 1999 epjc.11.19 . According to the Wolfenstein parameterization of the CKM matrix, there is an approximate relation between its elements up to . However, the measurement precision of the CKM element from both leptonic and semileptonic meson decays is generally about an order of magnitude smaller than that of from leptonic and semileptonic meson decays for the moment. The most precise values are from the global fit in SM, pdg2020 with uncertainties .
The information about the decay constant has not yet been obtained experimentally by now. Some theoretical results on are listed in Table 3. The theoretical discrepancies among various methods are obvious. In our calculation, as a conservative estimate, we will take the recent value MeV PhysRevD.100.014026 from the light front quark model, which agrees basically with the values MeV cpc.45.023109 from the recent lattice QCD simulation.
After some simple computation with Eq.(8), we obtain the partial decay widths and branching ratios for the decays as follows.
(28) | |||||
(29) | |||||
(30) | |||||
(31) |
These branching ratios are consistent with those of Ref. ctp.67.655 if the different values of decay constants are considered. The relatively large uncertainties of branching ratios come from the uncertainties of mass , width , decay constant and the CKM element . To experimentally study the decays, more than events are needed. Due to the short lifetime of lepton and the lepton number conservation in decays, additional neutrinos will make the measurement of the decay to have a poor reconstruction efficiency and to be very challenging. Perhaps some or more events are necessarily required to study the decay.


Above the open charm production threshold, there are several charmonium resonances and charmonium-like structures decaying predominantly into pairs of charmed meson final states. The studies of Belle PhysRevD.97.012002 , BaBar PhysRevD.79.092001 and CLEO-c PhysRevD.80.072001 ; cpc.42.043002 collaborations have shown that there is a sharply peaked structure and a broad plateau just above threshold, as illustrated in Fig. 4. Assuming the exclusive cross sections near threshold and , there will be about events corresponding to the total integrated luminosity of at future STCF, and about events corresponding to a data sample of at SuperKEKB. In addition, about bosons will be produced on the on the schedule of the large international scientific project of Circular Electron Positron Collider (CEPC) cepc and bosons at Future Circular Collider (FCC-ee) fcc . Considering the branching ratio pdg2020 , the boson decays will yield more than events at the tera-Z factories. So the decays could be investigated at Belle-II, SCTF, CEPC and FCC-ee experiments.
In hadron-hadron collisions, the inclusive cross sections for the pair and meson production are measured to be b and b at the center-of-mass energy of TeV by the LHCb group, with the transverse momentum within the range of JHEP.2016.03.159 . Some events could be accumulated with the integrated luminosity at LHCb. The total cross sections of charm and production measured at TeV by the ALICE group are mb and mb, respectively JHEP.2012.07.191 . The total cross sections of charm production measured at TeV by the ATLAS group are mb npb.907.717 . The production cross section at ATLAS should be very close to that at ALICE based on an educated guess. In addition, the meson can also produced from decays with the fragmentation fraction about epjc.75.19 . The -quark production cross sections at 13 TeV determined by LHCb and ALICE are about b JHEP.2015.10.172 and b 2108.02523 , respectively. So more than events from decays could be accumulated with the integrated luminosity at LHCb. All in all, the large cross section of meson plus the high luminosity at hadron-hadron collisions result in the abundant events, and make the carefully experimental study of the , decays, even the decay, to be possible and practicable.
V decays
The mesons have explicitly nonzero quantum number of electric charges, charm and strange, . Considering the conservation of the charm and strange quantum number in the strong and electromagnetic interactions, and the mass of mesons, MeV pdg2020 , being just above the threshold of pair but below the threshold of pair, the decays are the only allowable hadronic decay modes. However, the decays are highly suppressed due to four factors: (1) from the dynamical view, the decays are induced by the the electromagnetic interactions rather than the strong interactions because of the isospin non-conservation between the initial and final states, (2) from the perspective of the conservation of angular momentum, the orbital angular momentum of final states should be , so the decays are induced by the contributions of the -wave amplitudes, (3) from the phenomenological view, the decays are suppressed by the the Okubo-Zweig-Iizuka rules zweig ; ozi-o ; ozi-i because the quark lines of pion disconnect from those of the system, (4) from the kinematic view, the phase spaces of final states are very compact because of 9 MeV. Hence, the branching ratio for the hadronic decay is very small pdg2020 . And the branching ratio of the electromagnetic radiative decay is dominant, pdg2020 . Except for the , and final states, other decay modes of the mesons have not yet been observed pdg2020 . The weak decays are favored by the CKM element . The information about the can be experimentally obtained from the decays.

NRQM | PhysRevD.71.113006 | epl.115.21002 111Without QCD radiative corrections. | epl.115.21002 222With QCD radiative corrections. | epjp.132.80 111Without QCD radiative corrections. | epjp.132.80 222With QCD radiative corrections. |
---|---|---|---|---|---|
NRQM | plb.635.93 | ahep.2018.7032041 111Without QCD radiative corrections. | ahep.2018.7032041 222With QCD radiative corrections. | mpla.17.803 | |
RQM | plb.635.93 | mpla.17.803 | PhysRevD.55.6944 333With constituent quark masses for the light quarks and . | PhysRevD.55.6944 444With current quark masses for the light quarks and . | PhysRevD.96.016017 |
LFQM | PhysRevD.75.073016 555With Coulomb plus linear potential model. | PhysRevD.75.073016 666With Coulomb plus harmonic oscillator potential model. | PhysRevD.81.114024 777With the Gaussian type wave functions. | PhysRevD.81.114024 888With the power-law type wave functions. | PhysRevD.98.114018 |
LFQM | PhysRevC.92.055203 999With Coulomb plus linear potential model. | cpc.42.073102 101010With a dilation parameter GeV. | cpc.42.073102 111111With a dilation parameter GeV. | jpg.39.025005 | PhysRevD.100.014026 |
LFQM | epjc.76.313 777With the Gaussian type wave functions. | epjc.76.313 888With the power-law type wave functions. | epjp.133.134 121212With Martin potential model plb.93.338 . | epjp.133.134 131313With Cornell potential model zpc.33.135 . | epjp.133.134 141414With logarithmic potential model plb.71.153 . |
LQCD | PhysRevD.60.074501 | npb.619.507 151515With npb.619.507 and MeV pdg2020 . | jhep.1202.042 | ||
LQCD | PhysRevD.96.034524 | PhysRevLett.112.212002 | cpc.45.023109 | ||
SR | PhysRevD.88.014015 | epjc.75.427 | plb.735.12 | ijmpa.30.1550116 | |
other | epjp.133.134 161616With harmonic plus Yukawa potential model epjp.132.80 . | PhysRevD.75.116001 | plb.633.492 | PhysRevD.58.014007 |
The direct determinations of the CKM element come mainly from leptonic decays and semileptonic decays, as shown in Fig. 5. The uncertainties of from the leptonic decays, about , are dominated by the experimental uncertainties. The uncertainties of from the semileptonic decays, about , are dominated by the theoretical calculations of the form factors. It is wroth noting that the recent CKM element determined by the BES-III group from the and decays based on available data is PhysRevD.104.052009 , where the systematic (second) uncertainties has outweighed the statistical (first) one. This value is very close to the precise result from the global fit, pdg2020 that will be used in this paper.
By now, a relatively little information about the properties of the mesons is available. For example, the quantum number of , the decay constant , and the width have not yet been determined or confirmed explicitly by experiments. It is generally thought that the of the mesons is consistent with from decay modes PhysRevLett.75.3232 . Some theoretical results on the decay constant are listed in Table 4. It can be seen that the theoretical results are various. The recent LQCD results on the decay constant from ETM PhysRevD.96.034524 , HPQCD PhysRevLett.112.212002 and QCD cpc.45.023109 groups are in reasonable agreement with each other within an error range. The latest decay constant MeV from LQCD calculation cpc.45.023109 will be used for an estimation for PLDCV of the mesons in this paper. The experimental upper limit of the decay width is MeV at the 90 % confidence level set by the CLEO collaboration in 1995 PhysRevLett.75.3232 . An approximate relation for the decay width, , is often used in theoretical calculation. The radiative transition process, , is a parity conserving decay. The parity and angular momentum conservation implies that the orbital angular momentum of final states . There are many theoretical calculation on the decay width , for example, Refs. PhysRevLett.112.212002 ; PhysRevD.18.2537 ; PhysRevD.21.203 ; PhysRevD.31.1081 ; PhysRevD.37.2564 ; plb.284.421 ; plb.334.169 ; PhysRevD.47.1030 ; plb.316.555 ; plb.334.175 ; plb.336.113 ; PhysRevD.49.299 ; zpc.67.633 ; PhysRevD.52.6383 ; mpla.12.3027 ; npa.658.249 ; npa.671.380 ; jpg.27.1519 ; PhysRevD.64.094007 ; epja.13.363 ; plb.537.241 ; PhysRevD.68.054024 ; PhysRevD.72.094004 ; ijmpa.25.2063 ; epjc.75.243 ; epja.52.90 ; ijmpa.31.1650109 ; PhysRevD.94.113011 ; PhysRevD.100.016019 ; PhysRevD.101.054019 ; jhep.2020.04.023 ; 2106.13617 . The partial decay width for the magnetic dipole transition is generally written as fayyazuddin ,
(32) |
with the definition of the magnetic dipole moment and the momentum of photon in the rest frame of the vector meson,
(33) |
(34) |
where and are the electric charge in the unit of and mass of the constituent quark, respectively. With MeV, MeV, MeV and the
(35) |
(36) |
one can obtain keV and keV fayyazuddin . The theoretical value of partial decay width is roughly consistent with the corresponding experimental data keV within regions pdg2020 . For the moment, we will use keV in the calculation to give an estimate of branching ratios for the decays.
(37) | |||||
(38) | |||||
(39) | |||||
(40) |
If considering the experimental measurement efficiency, there are at least more than events to experimentally study the decays. And more than events might be needed to explore the decay.


In the electron-positron collisions, the cross sections of and production have been experimentally studied by the Belle PhysRevD.83.011101 , BaBar PhysRevD.82.052004 and CLEO-c PhysRevD.80.072001 ; cpc.42.043002 groups, as illustrated in Fig. 6. Assuming the exclusive cross sections near threshold and , there will be about events corresponding to a data sample of at STCF, and about events corresponding to a data sample of at SuperKEKB. In addition, considering the branching ratio pdg2020 and the fragmentation fraction epjc.76.397 , there will be more than (and ) events corresponding to cepc (and fcc ) bosons at the future CEPC (and FCC-ee). So the decays (with , and ) could be measured at Belle-II, SCTF, CEPC and FCC-ee experiments.
In hadron-hadron collisions, the inclusive cross sections for the pair production are mb at the center-of-mass energy of TeV at LHCb JHEP.2016.03.159 , mb and mb at TeV at ALICE JHEP.2012.07.191 and ATLAS npb.907.717 , respectively. With the fragmentation fraction epjc.76.397 , there will be about events corresponding a data sample of at LHCb, and more events available at ALICE and ATLAS. So the , , decays could be measured precisely at LHCb, ALICE and ATLAS experiments.
VI decays
The experimental information about the mesons are very scarce. The already known information about the mesons are their quark composition with the quark model assignment, the isospin , the spin-parity quantum number and the mass MeV pdg2020 . Due to the mass difference MeV , the electromagnetic radiative transition certainly will be the important and dominant decay mode. The photon in the decay is very soft, with the momentum MeV in the center-of-mass of the mesons. No signal event of the decay has yet been found. The decays offer a complementary decay modes of the meson. It can be seen from Eq.(8) that the information about could be obtained, however, the partial width for the purely leptonic decays are highly suppressed by the CKM element of .
The precise determinations of the CKM element are very central and important to verify the CKM picture of SM, where is the angle of the unitarity triangle of . The experimental determination of from the inclusive decay is complicated mainly by the large backgrounds from the CKM-favored decay. The experimental extraction of from exclusive decay is subject to the form factors calculated with the lattice QCD or QCD sum rules. The latest values obtained from inclusive and exclusive determinations are pdg2020
(41) |
(42) |
It is clearly seen that (1) the difference between inclusive and exclusive determinations of is obvious. (2) The best determinations of are from exclusive semileptonic decays, with a precision of about . The experimental errors for the exclusive semileptonic decays are expected to decrease from the current to based on a dataset of at Belle-II experiments PTEP.2019.123C01 . (3) The theoretical uncertainties are larger than experimental ones. Moreover, can also be experimentally determined from the leptonic decay and the semileptonic hyperon decay . The constraint from the global fit gives pdg2020 , which will be used in our calculation.
NRQM | plb.635.93 | PhysRevD.71.113006 | epjp.132.80 111Without QCD radiative corrections. | epjp.132.80 222With QCD radiative corrections. | |
---|---|---|---|---|---|
NRQM | ahep.2018.7032041 111Without QCD radiative corrections. | ahep.2018.7032041 222With QCD radiative corrections. | mpla.17.803 | epl.116.31004 111Without QCD radiative corrections. | epl.116.31004 222With QCD radiative corrections. |
RQM | plb.635.93 | mpla.17.803 | PhysRevD.55.6944 333With , and the constituent quark masses for the light quarks and . | PhysRevD.55.6944 444With , and the current quark masses for the light quarks and . | |
LFQM | PhysRevC.92.055203 555With Coulomb plus linear potential model. | cpc.42.073102 666With a dilation parameter GeV. | cpc.42.073102 777With a dilation parameter GeV. | jpg.39.025005 | PhysRevD.100.014026 |
LFQM | PhysRevD.75.073016 888With Coulomb plus linear potential model. | PhysRevD.75.073016 999With Coulomb plus harmonic oscillator potential model. | PhysRevD.81.114024 101010With the Gaussian type wave functions. | PhysRevD.81.114024 111111With the power-law type wave functions. | PhysRevD.98.114018 |
LFQM | epjc.76.313 101010With the Gaussian type wave functions. | epjc.76.313 111111With the power-law type wave functions. | epjp.133.134 121212With Martin potential model plb.93.338 . | epjp.133.134 131313With Cornell potential model zpc.33.135 . | epjp.133.134 141414With logarithmic potential model plb.71.153 . |
LQCD | PhysRevD.60.074501 | npb.619.507 151515With npb.619.507 and MeV pdg2020 . | PhysRevD.96.034524 | PhysRevD.91.114509 | |
SR | PhysRevD.88.014015 | epjc.75.427 | ijmpa.30.1550116 | 2106.13617 | PhysRevD.91.116009 |
other | epjp.133.134 161616With harmonic plus Yukawa potential model epjp.132.80 . | PhysRevD.75.116001 | plb.633.492 | PhysRevD.58.014007 |
Some theoretical results on the decay constant are collected in Table 5. The recent result MeV from LQCD calculation PhysRevD.96.034524 will be used in this paper. There are many theoretical calculation on the decay width , for example, Refs. PhysRevD.47.1030 ; plb.316.555 ; plb.334.175 ; PhysRevD.49.299 ; zpc.67.633 ; mpla.12.3027 ; npa.658.249 ; npa.671.380 ; jpg.27.1519 ; PhysRevD.64.094007 ; epja.13.363 ; plb.537.241 ; PhysRevD.68.054024 ; ijmpa.25.2063 ; epja.52.90 ; PhysRevD.94.113011 ; PhysRevD.100.016019 ; jhep.2020.04.023 ; 2106.13617 . With the formula of Eq.(32), the quark mass MeV fayyazuddin and GeV pdg2020 , and the magnetic dipole momentum
(43) |
one can obtain eV. The partial decay width and branching ratios for the decays are
(44) | |||||
(45) | |||||
(46) | |||||
(47) |
It is expected that there should be at least more than events available for experimental study of the decays.
The experimental study has shown that the exclusive cross sections for the final states of , and will have a large share of the total cross sections above the open bottom threshold, for example pdg2020 ,
(48) | |||
(49) | |||
(50) | |||
(51) | |||
(52) | |||
(53) |
There are about events corresponding to the dataset of at Belle experiments at the disposal PTEP.2019.123C01 . About events with a dataset at Belle-II are an outside estimate. Assuming the inclusive branching ratio , there will be some events at most at Belle-II. And it is more important that the vast majority of the data will be taken at resonance rather than mesons at Belle-II experiments, and lies below the threshold. So the probability of direct observation of the decays at Belle-II experiments should be very tiny. Considering about bosons at FCC-ee fcc and branching ratio pdg2020 , and assuming the fragmentation fraction plb.576.29 , there will be more than events to search for the decays. The -quark production cross sections at the center-of-mass energy 13 TeV is about b at LHCb JHEP.2015.10.172 . There will be more than events with a dateset of at LHCb and fragmentation fraction . Hence, the , , decays could be investigated at FCC-ee and LHCb experiments in the future.
VII decays
According to the conventional quark-model assignments, the mesons consist of two heavy quarks with different flavor numbers . Up to today, the experimental information of the meson is still very limited. For example, the potential candidate of the meson has not yet been determined. It is generally believed that the mass of the meson should be in the region between MeV recently measured by LHCb JHEP.2020.07.123 and MeV obtained by LHCb PhysRevLett.122.232001 (or MeV given by CMS PhysRevLett.122.132001 ), where the , and particles correspond to the sibling isoscalar states with quantum numbers of , and , respectively. So the branching ratios for the strong decays are zero, because meson is below the pair threshold. The experimental particle physicists are earnestly looking for and identifying the meson, a long-expected charming beauty. For the moment, almost all of the information available about the properties of meson (such as the mass, decay constant, lifetime, decay modes and so on) come from theoretical estimates. There are too many estimations on the meson mass with various theoretical models, for example, in Refs. zpc.3.165 ; PhysRevD.21.3180 ; PhysRevD.23.2724 ; PhysRevD.24.132 ; zpc.12.63 ; PhysRevD.32.189 ; ijmpa.6.2309 ; PhysRevD.44.212 ; PhysRevD.46.1165 ; PhysRevD.49.5845 ; zpc.64.57 ; Phys.Usp.38.1 ; PhysRevD.51.1248 ; PhysRevD.51.3613 ; PhysRevD.52.5229 ; PhysRevD.53.312 ; plb.382.131 ; epjc.4.107 ; PhysRevD.59.094001 ; PhysRevD.60.074006 ; PhysRevD.62.114024 ; PhysRevD.70.054017 ; ijmpa.19.1771 ; epjc.37.323 ; PhysRevD.71.034006 ; Pramana.66.953 ; plb.651.171 ; PhysRevC.78.055202 ; jpg.36.035003 ; npa.848.299 ; PhysRevD.81.076005 ; PhysRevLett.104.022001 ; PhysRevD.81.071502 ; epjc.71.1825 ; pan.74.631 ; PhysRevD.86.094510 ; epja.49.131 ; epja.50.154 ; ijmpa.32.1750021 ; PhysRevLett.121.202002 ; epjc.78.592 ; PhysRevD.99.054025 ; PhysRevD.99.096020 ; PhysRevD.100.096002 ; PhysRevD.100.114032 ; epja.55.82 ; npb.947.114726 ; epjc.80.223 ; plb.807.135522 ; PhysRevD.101.056002 ; PhysRevD.102.034002 ; pan.83.634 ; plb.816.136277 ; fewbody.62.39 ; epjc.81.327 . The recent result from lattice QCD calculation, MeV PhysRevLett.121.202002 , which are basically consistent with other estimations, will be used in this paper. Clearly, it is foreseeable that the isospin violating decay is explicitly forbidden by the law of energy conservation, because of MeV . Hence, the electromagnetic radiative transition should be the dominant decay mode. In addition, the photon in the magnetic dipole transition is very soft in the rest frame of the meson. This might be one main reason why the unambiguously experimental identification of the meson is very challenging. As an important complementary decay modes, the meson has very rich weak decay channels, which could be approximately classified into three classes: (1) the valence quark weak decay accompanied by the spectator quark, (2) the valence quark weak decay accompanied by the spectator quark, and (3) the and quarks annihilation into a virtual boson. The purely leptonic decays belong to the third case, which are favored by the CKM element . And the information of could be obtained from the decays.
NRQM | mpla.17.803 | epjc.78.592 | pan.83.634 | ||
---|---|---|---|---|---|
RQM | mpla.17.803 | PhysRevD.51.3613 111With Martin potential model plb.93.338 . | PhysRevD.51.3613 222With Coulomb plus linear potential model. | PhysRevD.51.3613 333Obtained from the scaling relation. | |
LFQM | PhysRevC.92.055203 | jpg.39.025005 | PhysRevD.81.114024 444With the Gaussian type wave functions. | PhysRevD.81.114024 555With the power-law type wave functions. | PhysRevD.98.114018 |
LFQM | epjc.76.313 444With the Gaussian type wave functions. | epjc.76.313 555With the power-law type wave functions. | PhysRevD.80.054016 666With Coulomb plus linear potential model. | PhysRevD.80.054016 777With Coulomb plus harmonic oscillator potential model. | PhysRevD.98.114038 |
LQCD | PhysRevD.91.114509 | lattice2018.273 | |||
SR | epja.49.131 888With the current quark mass. | epja.49.131 999With the pole quark mass. | npb.947.114726 | plb.807.135522 101010With inputs from the inverse Laplace-type model. | plb.807.135522 111111With inputs from Heavy Quark Symmetry. |
other | PhysRevD.75.116001 | plb.633.492 | PhysRevD.99.054025 |
The current values of the CKM element come mainly from inclusive and exclusive semileptonic decays of meson to charm pdg2020 . The average values obtained from inclusive decays and exclusieve decays are and , respectively pdg2020 . The lepton flavor non-universality in the ratio complicate the determination of . In addition, can also be obtained from the PLDCM decays, although none of the measurements has reached a competitive level of precision due to either the serious helicity suppression for , decays or other additional neutrinos from decay for decay. The global SM fit value is pdg2020 , which will be used in this paper.
Both valence quarks of the mesons are regarded as heavy quarks. Their Compton wave lengths are much shorter than a typical hadron size. The spin-flavor symmetry in the heavy quark limit would lead to an approximation between decay constants . Some theoretical results on the decay constant are collected in Table 6. The recent lattice QCD calculation MeV lattice2018.273 will be used in this paper. As it is well known that the magnetic momentum of both and quarks are inversely proportional to their mass. The magnetic dipole momentum
(54) |
should be very small. With the quark mass GeV and GeV, one can obtain eV using the formula of Eq.(32). The partial decay width and branching ratios for the decays are estimated to be,
(55) | |||||
(56) | |||||
(57) | |||||
(58) |
To experimentally investigate the decays, there should be at least more than events available.
More than bosons are expected at the future colliders of CEPC cepc and FCC-ee fcc . Considering the branching ratio pdg2020 and fragmentation fraction npa.953.21 ; cpc.43.083101 ; PhysRevD.100.034004 , there will be more than events to search for the , , decays. In addition, the production cross sections at LHC are estimated to be about nb for collisions at TeV, about mb for -Pb collisions at TeV and some mb for Pb-Pb collisions at TeV, respectively PhysRevD.97.114022 . There will be more than events corresponding to a dataset of at LHCb for collisions. Hence, the , , decays are expected to be carefully measured at LHCb experiments in the future.
VIII summary
The mass of the charged vector mesons are generally larger than that of the corresponding ground pseudoscalar mesons. The vector mesons decay mainly through the strong or/and electromagnetic interactions. These facts will inevitably result in that the branching ratios of the vector meson weak decays are often very tiny. Inspired by the potential prospects of existing and coming high-luminosity experiments, more and more experimental data will be accumulated, and higher measurement precision level will be reached. The probabilities of experimental investigation on the purely leptonic decays of charged vector mesons are discussed in this paper. We found that (1) for both and mesons, their widths are large due to the dominance of strong decay. Their PLDCV branching ratios are estimated at the order of . Although extremely complicated and difficult, the PLDCV decays , , might be measurable due to the huge data of the and mesons at LHCb. (2) The PLDCV decays are favored by the CKM element . Their branching ratios are about . The PLDCV decays , , could be carefully studied at the Belle-II, SCTF or STCF, CEPC, FCC-ee, LHCb experiments. (3) For the mesons below the thresholds and the mesons below both and thresholds, they decay predominantly through the magnetic dipole transitions. The branching ratios of the PLDCV decays favored by the CKM element could reach up to . The PLDCV decays , , might be searched for at the CEPC, FCC-ee, LHCb experiments. Our rough estimations and findings are summed in Table 7. We wish that our investigation could provoke physicists’ researching interest in PLDCV and offer a ready reference for the future experimental analysis.
decay modes | Belle-II | SCTF/STCF | CEPC | FCC-ee | LHCb | |
---|---|---|---|---|---|---|
, | ||||||
, | ||||||
, , | ||||||
, , | ||||||
, , | ||||||
, , |
Acknowledgments
The work is supported by the National Natural Science Foundation of China (Grant Nos. 11705047, 11981240403, U1632109, 11547014), the Chinese Academy of Sciences Large-Scale Scientific Facility Program (1G2017IHEPKFYJ01) and the Program for Innovative Research Team in University of Henan Province (19IRTSTHN018). We thank Prof. Haibo Li (IHEP@CAS), Prof. Shuangshi Fang (IHEP@CAS), Prof. Frank Porter (Caltech), Prof. Antimo Palano (INFN), Prof. Chengping Shen (Fudan University), Dr. Xiao Han (Fudan University), Prof. Xiaolin Kang (China University of Geosciences), Ms. Qingping Ji (Henan Normal University), Ms. Huijing Li (Henan Normal University) for their kindly help and valuable discussion.
References
- (1) P. Zyla et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2020, 083C01 (2020).
- (2) M. Gell-Mann, Phys. Lett. 8, 214 (1964).
- (3) G. Zweig, CERN-TH-401, 402, 412 (1964).
- (4) G. Buchalla, A. Buras, M. Lautenbacher, Rev. Mod. Phys. 68, 1125, (1996).
- (5) N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963).
- (6) M. Kobayashi, T. Maskawa, Prog. Theor. Phys. 49, 652 (1973).
- (7) P. Ball, V. Braun, A. Lenz, JHEP 05, 004 (2006).
- (8) P. Ball, J. Jones, JHEP 03, 069 (2007).
- (9) J. Hardy, I. Towner, PoS CKM2016, 028 (2016).
- (10) M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98, 030001 (2018).
- (11) J. Forshaw, R. Sandapen, G. Shaw, Phys. Rev. D 69, 094013 (2004).
- (12) M. McDermott, R. Sandapen, G. Shaw, Eur. Phys. J. C 22, 655 (2002).
- (13) K. Golec-Biernat, M. Wusthoff, Phys. Rev. D 60, 114023 (1999).
- (14) M. Ahmady, R. Sandapen, N. Sharma, Phys. Rev. D 94, 074018 (2016).
- (15) A. Bakulev, S. Mikhailov, Phys. Lett. B 436, 351 (1998).
- (16) P. Ball, V. Braun, Phys. Rev. D 54, 2182 (1996).
- (17) V. Chernyak, A. Zhitnitsky, Phys. Rept. 112, 173 (1984).
- (18) P. Ball, V. Braun, Phys. Rev. D 58, 094016 (1998).
- (19) P. Ball, R. Zwicky, JHEP 04, 046 (2006).
- (20) J. He, B. Juliá-Díaz, Y. Dong, Eur. Phys. J. A 24, 411 (2005).
- (21) J. He, B. Juliá-Díaz, Y. Dong, Phys. Lett. B 602, 212 (2004).
- (22) D. Ebert, R. Faustov, V. Galkin, Phys. Lett. B 635, 93 (2006).
- (23) H. Choi, C. Ji, Phys. Rev. D 75, 034019 (2007).
- (24) H. Choi, C. Ji, Z. Li, H. Ryu, Phys. Rev. C 92, 055203 (2015).
- (25) Q. Chang, X. Li, X. Li, F. Su, Chin. Phys. C 42, 073102 (2018).
- (26) R. Verma, J. Phys. G 39, 025005 (2012).
- (27) Q. Chang, X. Li, X. Li, et al., Phys. Rev. D 98, 114018 (2018).
- (28) N. Dhiman, H. Dahiya, C. Ji et al., Phys. Rev. D 100, 014026 (2019).
- (29) A. Khan et al. (CP-PACS Collaboration), Phys. Rev. D 65, 054505 (2002); 67, 059901 (2003)(E).
- (30) K. Hashimoto, T. Izubuchi, Prog. Theor. Phys. 119, 599 (2008).
- (31) K. Jansen, C. McNeile, C. Michael et al., Phys. Rev. D 80, 054510 (2009).
- (32) V. Braun, P. Bruns, S. Collins et al., JHEP 04, 082 (2017).
- (33) W. Sun et al. (QCD Collaboration), Chin. Phys. C 42, 063102 (2018).
- (34) P. Maris, P. Tandy, Phys. Rev. C 60, 055214 (1999).
- (35) Y. Simonov, Phys. Atom. Nucl. 79, 444 (2016).
- (36) Y. Grossman, M. König, M. Neubert, JHEP 04, 101 (2015).
- (37) W. Marciano, A. Sirlin, Phys. Rev. Lett. 61, 1815 (1988).
- (38) E. Braaten, C. Li, Phys. Rev. D 42, 3888 (1990).
- (39) B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 78, 071103 (2008).
- (40) C. Adamuščín, G. Gakh, E. Tomasi-Gustafsson, Phys. Rev. C 75, 065202 (2007).
- (41) J. Melo, C. Ji, T. Frederico, Phys. Lett. B 763, 87 (2016).
- (42) E. Kou et al., Prog. Theor. Exp. Phys. 123C01 (2019); 029201 (2020)(E).
- (43) T. Uglov, EPJ Web of Conf. 212, 01010 (2019).
- (44) C. Yuan, M. Karliner, Phys. Rev. Lett. 127, 0120031 (2021).
- (45) X. Shi, X. Zhou, X. Qin, H. Peng, JINST 16, P03029 (2021).
- (46) A. Barnyakov et al., Nucl. Inst. Meth. Phys. Res. A 958, 162352 (2020).
- (47) H. Haber, J. Perrier, Phys. Rev. D 32, 2961 (1985).
- (48) A. Bevan et al., Eur. Phys. J. C 74, 3026 (2014).
- (49) X. Kang, PoS EPS-HEP2019, 511 (2020).
- (50) R. Aaij et al. (LHCb Collaboration), JHEP 10, 172 (2015); 05, 063 (2017)(E).
- (51) I. Bediaga et al. (LHCb Collaboration), arXiv:1808.08865.
- (52) R. Aaij et al. (LHCb Collaboration), JHEP 03, 159 (2016); 09, 013 (2016)(E); 05, 074 (2017)(E).
- (53) A. Abada et al., Eur. Phys. J. Special Topics 228, 1109 (2019).
- (54) P. Dimopoulos, F. Mescia, A. Vladikas, Phys. Rev. D 84, 014505 (2011).
- (55) P. Vilain et al. (CHARM-II Collaboration), Eur. Phys. J. C 11, 19 (1999).
- (56) C. Albertus, E. Hernández, J. Nieves et al., Phys. Rev. D 71, 113006 (2005).
- (57) B. Yazarloo, H. Mehraban, EPL. 115, 21002 (2016).
- (58) B. Yazarloo, H. Mehraban, Eur. Phys. J. Plus 132, 80 (2017).
- (59) M. Abu-Shady, E. Khokha, Advances in High Energy Phys. 7032041 (2018).
- (60) D. Ebert, R. Faustov, V. Galkin, Mod. Phys. Lett. A 17, 803 (2002); Phys. Rev. D 67, 014027 (2003).
- (61) D. Hwang, G. Kim, Phys. Rev. D 55, 6944 (1997).
- (62) H. Choi, Phys. Rev. D 75, 073016 (2007).
- (63) C. Hwang, Phys. Rev. D 81, 114024 (2010).
- (64) C. Geng, C. Lih, C. Xia, Eur. Phys. J. C 76, 313 (2016).
- (65) N. Dhiman, H. Dahiya, Eur. Phys. J. Plus 133, 134 (2017).
- (66) A. Martin, Phys. Lett. B 93, 338 (1980).
- (67) K. Hagiwara, A. Martin, A. Peacock, Z. Phys. C 33, 135 (1986).
- (68) C. Quigg, J. Rosner, Phys. Lett. B 71, 153 (1977).
- (69) D. Becirevic, Ph. Boucaud, J. Leroy et al., Phys. Rev. D 60, 074501 (1999).
- (70) K. Bowler et al. (UKQCD collaboration), Nucl. Phys. B 619, 507 (2001).
- (71) D. Bečirević, V. Lubicz, F. Sanfilippo et al., JHEP 02, 042 (2012).
- (72) V. Lubicz, A. Melis, S. Simula (ETM Collaboration), Phys. Rev. D 96, 034524 (2017).
- (73) Y. Chen (QCD Collaboration), Chin. Phys. C 45, 023109 (2021).
- (74) P. Gelhausen, A. Khodjamirian, A. Pivovarov et al., Phys. Rev. D 88, 014015 (2013); 89, 099901 (2013)(E); 91, 099901 (2015)(E).
- (75) Z. Wang, Eur. Phys. J. C 75, 427 (2015).
- (76) W. Lucha, D. Melikhov, S. Simula, Phys. Lett. B 735, 12 (2014).
- (77) S. Narison, Int. J. Mod. Phys. A 30, 1550116 (2015).
- (78) C. Albertus, E. Hernández, J. Nieves et al., Phys. Rev. D 75, 116001 (2007).
- (79) G. Wang, Phys. Lett. B 633, 492 (2006).
- (80) A. El-Hady, A. Datta, J. Vary, Phys. Rev. D 58, 014007 (1998).
- (81) B. Zhou, J. Sun, Y. Zhang, Commun. Theor. Phys. 67, 655 (2017).
- (82) V. Zhukova et al (Belle Collaboration), Phys. Rev. D 97, 012002 (2018).
- (83) B. Aubert et al (BaBar Collaboration), Phys. Rev. D 79, 092001 (2009).
- (84) D. Cronin-Hennessy et al (CLEO Collaboration), Phys. Rev. D 80, 072001 (2009).
- (85) X. Dong, L. Wang, C. Yuan, Chin. Phys. C 42, 043002 (2018).
- (86) J. Costa et al., IHEP-CEPC-DR-2018-02, arXiv:1811.10545.
- (87) A. Abada et al., Eur. Phys. J. C 79, 474 (2019).
- (88) B. Abelev et al. (ALICE Collaboration), JHEP 07, 191 (2012).
- (89) G. Aad et al. (ATLAS Collaboration), Nucl. Phys. B 907, 117 (2016).
- (90) L. Gladilin, Eur. Phys. J. C 75, 19 (2015).
- (91) S. Acharya et al. (ALICE Collaboration), arXiv:2108.02523.
- (92) S. Okubo, Phys. Lett. 5, 165 (1963).
- (93) J. Iizuka, Prog. Theor. Phys. Suppl. 37-38, 21 (1966).
- (94) M. Ablikim et al. (BESIII Collaboration), Phys. Rev. D 104, 052009 (2021).
- (95) N. Soni, J. Pandya, Phys. Rev. D 96, 016017 (2017); 99, 059901 (2019)(E).
- (96) J. Gronberg et al. (CLEO Collaboration), Phys. Rev. Lett. 75, 3232 (1995).
- (97) G. Donald, C. Davies, J. Koponen, G. Lepage, Phys. Rev. Lett. 112, 212002 (2014).
- (98) R. Hackman, N. Deshpande, D. Dicus et al., Phys. Rev. D 18, 2537 (1978).
- (99) E. Eichten, K. Gottfried, T. Kinoshita et al., Phys. Rev. D 21, 203 (1980).
- (100) W. Wilcox, O. Maxwell, K Milton, Phys. Rev. D 31, 1081 (1985).
- (101) G. Miller, P. Singer, Phys. Rev. D 37, 2564 (1988).
- (102) A. Kamal, P. Xu, Phys. Lett. B 284, 421 (1992).
- (103) T. Aliev, E. Iltan, N. Pak, Phys. Lett. B 334, 169 (1994).
- (104) P. O’Donnell, P. Xu, Phys. Lett. B 336, 113 (1994).
- (105) G. Burdman, E. Golowich, J. Hewett et al., Phys. Rev. D 52, 6383 (1995).
- (106) F. Close, E. Swanson, Phys. Rev. D 72, 094004 (2005).
- (107) G. Yu, Z. Li, Z. Wang, Eur. Phys. J. C 75, 243 (2015).
- (108) N. Li, Y. Wu, Int. J. Mod. Phys. A 31, 1650109 (2016).
- (109) B. Yang, B. Wang, L. Meng et al., Phys. Rev. D 101, 054019 (2020).
- (110) H. Cheng, C. Cheung, G. Lin et al., Phys. Rev. D 47, 1030 (1993).
- (111) P. Colangelo, F. Fazio, G. Nardulli, Phys. Lett. B 316, 555 (1993).
- (112) P. Colangelo, F. Fazio, G. Nardulli, Phys. Lett. B 334, 175 (1994).
- (113) N. Barik, P. Dash, Phys. Rev. D 49, 299 (1994); 53, 4110 (1996)(E).
- (114) M. Ivanov, Y, Valit, Z. Phys. C 67, 633 (1995).
- (115) S. Zhu, Z. Yang, Mod. Phys. Lett. A 12, 3027 (1997).
- (116) S. Jena, S. Panda, T. Tripathy, Nucl. Phys. A 658, 249 (1999).
- (117) Y. Dong, A. Faessler, K. Shimizu, Nucl. Phys. A 671, 380 (2000).
- (118) S. Jena, P. Panda, K. Sahu, J. Phys. G 27, 1519 (2001).
- (119) J. Goity, W. Roberts, Phys. Rev. D 64, 094007 (2001).
- (120) R. Bonnaz, B. Silvestre-Brac, C Gignoux, Eur. Phys. J. A 13, 363 (2002).
- (121) D. Ebert, R. Faustov, V. Galkin, Phys. Lett. B 537, 241 (2002).
- (122) W. Bardeen, E. Eichten, C. Hill, Phys. Rev. D 68, 054024 (2003).
- (123) S. Jena, M. Muni, H. Pattnaik et al., Int. J. Mod. Phys. A 25, 2063 (2010).
- (124) V. Šimonis, Eur. Phys. J. A 52, 90 (2016).
- (125) M. Priyadarsini, P. Dash, S. Kar et al., Phys. Rev. D 94, 113011 (2016).
- (126) B. Wang, B. Yang, L. Meng et al., Phys. Rev. D 100, 016019 (2019).
- (127) H. Li, C. Lü, C. Wang et al., JHEP 04, 023 (2020).
- (128) B. Pullin, R. Zwicky, JHEP 09, 023 (2021).
- (129) Fayyazuddin, Riazuddin, A modern introduction to particle physics, 3rd edition, 2011.
- (130) G. Pakhlova et al. (Belle Collaboration), Phys. Rev. D 83, 011101 (2011).
- (131) P. Sanchez et al. (BaBar Collaboration), Phys. Rev. D 82, 052004 (2010).
- (132) M. Lisovyi, A. Verbytskyi, O. Zenaiev, Eur. Phys. J. C 76, 397 (2016).
- (133) B. Yazarloo, H. Mehraban, EPL 116, 31004 (2016).
- (134) B. Colquhoun et al. (HPQCD Collaboration), Phys. Rev. D 91, 114509 (2015).
- (135) W. Lucha, D. Melikhov, S. Simula, Phys. Rev. D 91, 116009 (2015).
- (136) J. Abdallah et al. (DELPHI Collaboration), Phys. Lett. B 576, 29 (2003).
- (137) R. Aaij et al. (LHCb Collaboration), JHEP 07, 123 (2020).
- (138) R. Aaij et al. (LHCb Collaboration), Phys. Rev. Lett. 122, 232001 (2019).
- (139) A. Sirunyan et al. (CMS Collaboration), Phys. Rev. Lett. 122, 132001 (2019).
- (140) S. Nussinov, Z. Phys. C 3, 165 (1979).
- (141) D. Stanley, D. Robson, Phys. Rev. D 21, 3180 (1980).
- (142) E. Eichten, F. Feinberg, Phys. Rev. D 23, 2724 (1981).
- (143) W. Buchmüller, S. Tye, Phys. Rev. D 24, 132 (1981).
- (144) A. Kaidalov, Z. Phys. C 12, 63 (1982).
- (145) S. Godfrey, N. Isgur, Phys. Rev. D 32, 189 (1985).
- (146) S. Gershtein, A. Likhoded, S.Slabospitsky, Int. J. Mod. Phys. A 6, 2309 (1991).
- (147) W. Kwong, J. Rosner, Phys. Rev. D 44, 212 (1991).
- (148) Y. Chen, Y. Kuang, Phys. Rev. D 46, 1165 (1992); 47, 350 (1993)(E).
- (149) E. Eichten, C. Quigg, Phys. Rev. D 49, 5845 (1994).
- (150) E. Bagan, H. Dosch, P. Gosdzinsky et al., Z. Phys. C 64, 57 (1994).
- (151) S. Gershtein, V. Kiselev, A. Likhoded et al. Phys. Usp. 38, 1 (1995).
- (152) R. Roncaglia, A. Dzierba, D. Lichtenberg et al. Phys. Rev. D 51, 1248 (1995).
- (153) S. Gershtein, V. Kiselev, A. Likhoded et al., Phys. Rev. D 51, 3613 (1995).
- (154) J. Zeng, J. Orden, W. Roberts, Phys. Rev. D 52, 5229 (1995).
- (155) S. Gupta, J. Johnson, Phys. Rev. D 53, 312 (1996).
- (156) C. Davies, K. Hornbostel, G. Lepage et al., Phys. Lett. B 382, 131 (1996).
- (157) L. Motyka, K. Zalewski, Eur. Phys. J. C 4, 107 (1998).
- (158) A. El-Hady, M. Lodhi, J. Vary, Phys. Rev. D 59, 094001 (1999).
- (159) L. Fulcher, Phys. Rev. D 60, 074006 (1999).
- (160) M. Baldicchi, G. Prosperi, Phys. Rev. D 62, 114024 (2000).
- (161) S. Godfrey, Phys. Rev. D 70, 054017 (2004).
- (162) S. Ikhdair, R. Sever, Int. J. Mod. Phys. A 19, 1771 (2004).
- (163) D. Li, B. Ma, Y. Li et al., Eur. Phys. J. C 37, 323 (2004).
- (164) A. El-Hady, J. Spence, J. Vary, Phys. Rev. D 71, 034006 (2005).
- (165) A. Rai, P. Vinodkumar, Pramana 66, 953 (2006).
- (166) T. Chiu et al. (TWQCD Collaboration), Phys. Lett. B 651, 171 (2007).
- (167) A. Rai, B. Patel, P. Vinodkumar, Phys. Rev. C 78, 055202 (2008).
- (168) B. Patel, P. Vinodkumar, J. Phys. G 36, 035003 (2009).
- (169) A. Parmar, B. Patel, P. Vinodkumar, Nucl. Phys. B 848, 299 (2010).
- (170) K. Wei, X. Guo, Phys. Rev. D 81, 076005 (2010).
- (171) E. Gregory et al. (HPQCD Collaboration), Phys. Rev. Lett. 104, 022001 (2010).
- (172) A. Badalian, B. Bakker, I. Danilkin, Phys. Rev. D 81, 071502 (2010); 81, 099902 (2010)(E).
- (173) D. Ebert, R. Faustov, V. Galkin, Eur. Phys. J. C 71, 1825 (2011).
- (174) A. Badalian, B. Bakker, I. Danilkin, Phys. Atom. Nucl. 74, 631 (2011).
- (175) R. Dowdall, C. Davies, T. Hammant et al., Phys. Rev. D 86, 094510 (2012).
- (176) Z. Wang, Eur. Phys. J. A 49, 131 (2013).
- (177) N. Devlani, V. Kher, A. Rai, Eur. Phys. J. A 50, 154 (2014).
- (178) A. Monteiro, M. Bhat, and K. Kumar, Int. J. Mod. Phys. A 32, 1750021 (2017); Phys. Rev. D 95, 054016 (2017).
- (179) N. Mathur, M. Padmanath, S. Mondal, Phys. Rev. Lett. 121, 202002 (2018).
- (180) N. Soni, B. Joshi, R. Shah et al., Eur. Phys. J. C 78, 592 (2018).
- (181) E. Eichten, C. Quigg, Phys. Rev. D 99, 054025 (2019).
- (182) Q. Li, M. Liu, L. Lu et al., Phys. Rev. D 99, 096020 (2019).
- (183) I. Asghar, F. Akram, B. Masud et al., Phys. Rev. D 100, 096002 (2019).
- (184) L. Gutiérrez-Guerrero, A. Bashir, A. Bedolla et al., Phys. Rev. D 100, 114032 (2019).
- (185) N. Akbar, F. Akram, B. Masud et al., Eur. Phys. J. A 55, 82 (2019).
- (186) T. Aliev, T. Barakat, S. Bilmis, Nucl. Phys. B 947, 114726 (2019).
- (187) P. Ortega, J. Segovia, D. Entem et al., Eur. Phys. J. C 80, 223 (2020).
- (188) S. Narison, Phys. Lett. B 807, 135522 (2020).
- (189) M. Chen, L. Chang, Y. Liu, Phys. Rev. D 101, 056002 (2020).
- (190) L. Abreu, F. Júnior, A. Favero, Phys. Rev. D 102, 034002 (2020).
- (191) N. Akbar, Phys. Atom. Nucl. 83, 634 (2020).
- (192) R. Ding, B. Wan, Z. Chen et al, Phys. Lett. B 816, 136277 (2021).
- (193) J. Ahmed, R. Manzoor, L. Chang et al., Few-Body Syst. 62, 39 (2021).
- (194) P. Yin, Z. Cui, C. Roberts et al, Eur. Phys. J. C 81, 327 (2021).
- (195) H. Choi, C. Ji, Phys. Rev. D 80, 054016 (2009).
- (196) S. Tang, Y. Li, P. Maris et al., Phys. Rev. D 98, 114038 (2018).
- (197) D. Becirevic et al. (ETM Collaboration), PoS LATTICE2018, 273 (2019)
- (198) G. Boroun, S. Zarrin, S. Dadfar, Nucl. Phys. A 953, 21 (2016).
- (199) D. Yang, W. zhang, Chin. Phys. C 43, 083101 (2019).
- (200) X. Zheng, C. Chang, T. Feng et al., Phys. Rev. D 100, 034004 (2019).
- (201) G. Chen, C. Chang, X. Wu, Phys. Rev. D 97, 114022 (2018).