LAL16-011

Novel approach to measure the leptonic decays via charmed meson decays

N. T. Huong, E. Kou and B. Viaud

Faculty of Physics, VNU University of Science, Vietnam National University

334 Nguyen Trai, Thanh Xuan, Hanoi, Vietnam

Laboratoire de l’Accélérateur Linéaire, Univ. Paris-Sud, CNRS/IN2P3 (UMR 8607)

Université Paris-Saclay, 91898 Orsay Cédex, France

Abstract

In this article, we propose a novel approach to measure the branching ratios of the leptonic decays by using charmed meson decays, namely, and . We advocate that the data available at LHCb can already yield a new measurement of with accuracy competitive with the current world average. We also estimate that using the data collected by LHCb between 2015 and 2018 in proton-proton collisions at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 5.0 fb, the relative uncertainties to this branching ratio can be reduced down to %. We also show that the first observation of may be possible with the Upgrade of the LHCb experiment.

## 1 Introduction

In this article, we discuss a possible improved measurement of the rare and () decays. The pseudoscalar mesons , being , decay into a two-photon final state while they can also decay into two leptons. The leptonic decays are very rare because of being loop induced processes as shown in Fig. 1 (a). This loop diagrams are ultraviolet divergent, which requires a local counter term. These contributions are parameterized by the chiral coupling constants, and [1], and they have to be determined by e.g. using the experimental measurements of the branching ratios of the leptonic decays of .

The determination of the coupling constants has important consequences in various phenomenological studies. Let us give some examples.

Recently, two of the authors have computed the pion exchange contributions to the hyperfine splitting of the muonic hydrogen [2]. This contribution which comes from the diagram given in Fig. 1 (b) contains the same coupling as Fig. 1 (a). It was found that the theoretical uncertainty for the pion exchange contribution to the hyperfine splitting is indeed dominated by the uncertainties in the parameters, thus, it can be improved by measurements of with higher accuracy.

An important motivation is the theoretical prediction of the muon anomalous magnetic moment. The muon is one of the most precisely measured observable in particle physics: . The theoretical prediction within the Standard Model (SM), is also very precise while the comparison shows 3.1 deviation, according to the latest review [3]: . As the experimental uncertainties will be reduced by future experiments namely at Fermilab and J-PARC [4], the dominant uncertainty will be theoretical one, in particular, the so-called hadronic vacuum polarization effect and the light-by-light contributions. The dominant contributions to the hadronic light-by-light amplitude from pseudoscalar meson exchange are as shown in Fig. 1 (c). This contribution can be described by two counter-terms in the chiral perturbation theory [5, 6] and one of them is the coupling constants. The importance of determining these parameters from the leptonic decays has been also re-emphasized recently in [3, 7], since in particular, a tension between theory and experiment has been reported in the case of [8].

In this article, we propose to measure the coupling constants by using the semi-leptonic decays of charmed mesons, in particular, and . These charm decays are recently studied in detail at LHCb [15, 17] for searching new physics signals (see [9, 10], for examples of new physics scenarios).

In section 2, we first introduce the counter-term in terms of chiral Lagrangian and constrain them by using available experimental data and compare to the theoretical model computation by the lowest meson dominant model (LMD), which is often applied to compute the hadronic light-by-light contributions of muon . Then, we propose a measurement of at LHCb and assess the future sensitivity in section 3 and we conclude in section 4.

## 2 Determining the counter-term of leptonic pseudoscalar decays

In this section, we derive ratios of decay rates:

(1) |

where . The decay widths are:

(2) |

which comes from the Wess-Zumino term:

(3) | |||||

where is the electric charge of the and quarks and where are the Gell-mann matrices extended to the nonet symmetry with and where . For simplicity, we identify to for now. As we will see later-on, the different factors for cancel in the ratio so that in fact, we do not need to consider the pseudoscalar mixing angle as long as we use this ratio.

The loop diagram in Fig. 1 (a) has been computed in [6, 11]. The real part diverges and it requires a local counter term, which can be written to the lowest order in the chiral expansion, as [5]:

(4) |

where . The coefficients are renormalization scale and scheme dependent. Using this counter term, we can write the decay rates as:

(5) |

where the amplitude is given with a loop function as [6]:

(6) |

where

(7) |

and . This result implies that () can be expressed in a simple form as:

(8) |

As mentioned earlier, the overall factor cancels.

— | |||

— | |||

— | |||

In [6], the counter term is computed using the Lowest Meson Dominant model (LMD) as:

(9) |

where systematic theoretical error of 40 % is assumed. Using this value, we compute for all the channels. The results are given in Table 1. We can see that the theoretical prediction for shows a good agreement with the LMD prediction while in the channels a small deviation appears. Note that for , we use the latest KTeV measurement of [8] which subtracts the radiative process.

The counter term, in turn, can be extracted from the experimental values up to two hold ambiguity, namely those of and (including NLO QED corrections as in [12, 13]):

(10) | |||||

(11) |

The current situation is summarized in Fig. 2. We drop the second solution of for now as it is far from any other measurement. As we can see, the experimental extraction of the parameter has much larger errors than the theoretical prediction based on LMD model: if we choose the first solution for , the average with leads to the range of:

(12) |

Further improvements in determining the parameters through the experimental measurement of the leptonic pseudoscalar decays are desired. In the next section, we present a novel way to extract this parameter via charmed meson decays.

## 3 Measuring the decays via charmed meson decays at LHCb

The LHC produces a huge number of charmed hadrons, like the , or mesons, whose decays can be studied by the LHCb collaboration. While investigating decays in the data collected in 2011 [15], LHCb found a few dozens of events peaking around the meson mass in the dimuon mass spectrum. These events seem to come from the high branching ratio channels , followed by the rare decay . The branching ratio evaluated with this data has a larger uncertainty than the present world average: =(5.80.8)10 [14]. While LHCb accumulates more and more data, the collaboration will be able to measure these peaks with much higher precision. In this section, we evaluate LHCb’s potential for the measurement of and its impact on the determination of the parameter.

We also investigate a possible observation of the decay at LHCb. As shown in Table 1, if we treat in the same way as and , we expect its leptonic branching ratio to be . However, as is much heavier than the other mesons, it is quite possible that the mass correction is sizable. Therefore, a measurement of this branching ratio is very important to understand the applicability of our theoretical framework.

### 3.1 Sensitivity of the LHCb experiment to

The most precise measurement of the decay should be obtained by measuring . In data collected in proton-proton collisions at a centre-of-mass energy of 7 TeV and corresponding to an integrated luminosity of 1.0 fb, LHCb observed (297) decays in the region of the dimuon mass spectrum surrounding the mass of the meson [15]. An evaluation of the corresponding branching fraction made in [15] yields: =(2.20.6)10. The measurement of this branching fraction was not the goal of the analysis described in [15] and a more detailled analysis is necessary to measure it. On the other hand, by combining the presently available estimate with =(3.660.22)10 [14], we can obtain a reasonable precision: =(6.01.7)10.

The decay can also be used. Its branching fraction is also evaluated in Ref. [15]: =(6.82.1)10. Combined with =(1.70.09)10 [14], this yields =(4.01.3)10. The weighted average of these two evaluations is =(4.71.0)10. This corresponds to a relative uncertainty of about 20%. This value can be compared to the PDG average =(5.80.8)10 quoted in the previous sections which is dominated by the SATURNE II collaboration measurement [16]. The LHCb constraint on the parameter is shown in Fig. 2. Although the relative error in the branching ratio is 20 % for LHCb and 14 % for PDG value, since the branching fraction has parabolic dependence on the , the LHCb result, which has a smaller central value, leads to a much larger errors in for now. A simple average of the LHCb and PDG values, or , is shown in Figs. 2 and 3 (the gray bound in Fig. 3). By averaging with the , we obtain:

(13) |

Now let us discuss possible future LHCb measurements. First of all, adding the sample collected in 2012 (2 fb at a centre-of-mass energy of 8 TeV), three times more decay candidates would be already available. Therefore, with the data already available at LHCb, can be measured with a relative uncertainty of the order of 13%.

The uncertainty above is dominated by the statistical uncertainty on and . On the other hand, as more data becomes available, the share of uncertainties on and (presently known up 5 to 6%) will increase. Also, such measurements are in practice normalized to the decays : the value of is derived from the ratio between the yields of and decays measured in the sample, combined with the branching fractions . The latter can be determined with an uncertainty of about 5% by combining the measured values of and [14] and assuming no difference between the decays and .

In the 5 fb data sample collected between 2015 and 2018 at a centre-of-mass energy of 13 TeV (LHCb Run II), about 3 times more -meson decays should be produced (larger integrated luminosity and twice larger production cross-section). Using this data, can therefore be measured with a relative statistical uncertainty of about 6%. Accounting also for the uncertainty on the normalisation mode and on , we conclude that can be measured up to about a 10% uncertainty with LHCb Run II’s data.

The measurement of can be further improved by using the decay . In. [17], LHCb measured the branching fraction of the process in the dimuon mass range 675 875 MeV/. Focussing on events triggered thanks to the muons in the final state, roughly 15000 decays can be isolated in the 2 fb sample used in [17], where 5000 background candidates also lies in the signal region of the . In the dimuon mass range 538 558 MeV/, where almost all of decays should be reconstructed, we expect about ten times less background candidates. Using ()=0.01 and =(5.80.8)10 (Table 1) to predict 5.810, we rescale the yield of signal events found in 675 875 MeV/, where the branching fraction is =(4.20.4)10 [17]. We predict that about 200 decays can be found in the 2 fb sample collected in 2012. This corresponds to a significance of about 7 and suggests can be measured with a statistical precision of about 15%. Combined with the 10% uncertainty on in 675 875 MeV/, which would be used as the normalisation mode in this case, this would yield a 20% precise measurement of . This also assumes that () can be measured precisely at Flavor Factories in the coming years. The value assumed above is derived from =(1.090.16)10 [14] by betting that the branching fraction of , the corresponding Cabibbo favored mode, is an order of magnitude higher. A 20% precise measurement through this mode would provide a valuable cross-check of the value obtained through the decays . Combining both measurements, an overall 10% precision on is possible.

Looking further, the Upgrade of LHCb is planed to collect 50 fb of data at =14 TeV. We expect precisions of the order of 8% for the measurement using decays and of 10% for that using decays. The uncertainties do not decrease much due to the fact that we are not accounting for the improvements in the the uncertainties on the normalization modes nor on and (). These errors can be indeed reduced by LHCb or Belle II, BESIII experiments independently, which will further reduce the errors in .

### 3.2 Sensitivity of the LHCb experiment to

According to the prediction in Table 1, the branching ratio of the decay is 40 times lower than that of . It will therefore be far more difficult to measure. It should however be possible using the sample collected with the Upgraded detector, by measuring , which is predicted to be as low as about 5.510. We evaluate the yield of such decays based on the yield of decays measured in [15], corrected by this ratio :

(14) |

where all the branching fractions but are taken from [14]. We also assume that the Upgraded trigger system will double the signal efficiency. We derive the background yields that will also contribute to the measurement’s uncertainty from the yields observed below the peak of the distribution as shown in [15]. We count only candidates found in the dimuon mass range where the was looked for (850-1250 MeV/), and scale this number by the width ratio between this window and the 10 MeV/ window that would be used to search for a . For both the signal and background predictions, we also account for the larger integrated luminosity and production cross-section. With about 250 decays and about 1000 background events expected in these conditions, a signal significance of 7 is possible and is measurable with a precision of about 15%. If so, will be measured with a similar precision, due to the low impact of the uncertainties in =(3.940.25)10 and in =(4.540.16)10 [14]. The decay is less competitive due to a lower branching fraction, [14], which is not compensated for by the higher production cross section (3 times higher than that of the ) in proton-proton collisions. Therefore, a 25 % precise measurement is possible in this channel.

The decay can also be studied through the decay . The same computation as in Sect. 3.1, where 1.110 is substituted for , can be carried out to conclude that even LHCb’s Upgrade can’t provide a precise measurement of such a small branching ratio. Another experimental approach is possible: a measurement focussing on decays where the originated from . The mass difference can be used to reduce the background level by two orders of magnitude, at the price of about 10 times less signal decays due to a production cross-section about 2.5 times lower, to =23 and to the low efficiency with which soft pions can be reconstructed and selected. Assuming once more that the Upgrade of the experiment will allow to raise the trigger efficiency by a factor 2, we can hope to measure and therefore with a precision of about 30%.

We conclude that the Upgraded LHCb experiment can produce the first precise measurement of with an uncertainty lower than 15 %.

## 4 Conclusions

The process is known to provide crucial information for various low energy phenomena, including the light-by-light contribution to the muon anomalous magnetic moment. The theoretical model which is often used to estimate the light-by-light contribution predicts rather well the branching ratio of the process. However, the latest experimental results on the branching ratio by the KTeV collaboration shows a tension with the prediction from the same model. Since average is dominated by a single measurement, a confirmation of this result is desirable (see [3, 7] for recent discussions).

In this article, we propose a novel approach to determine the branching ratios of decays using the high static LHCb data of charmed meson decays. Based on the decays already observed by LHCb in the 1 sample collected in 2011, we find that this branching fraction could be measured with a precision of 20 %. This should be compared with the 13% uncertainty of the present world average: . Using the LHCb Run II data (, =13 TeV), we estimate that the statistical uncertainties can be reduced down to % level. The decay chain can be also used to further improve the precision. We estimate that using 2 from Run I data, can be determined at the 20 % level from this channel. This assumes that the branching fraction of the first stage of the decay chain, , can be determined at the few percent level using non leptonic decays of the , which seems in the reach of Belle II and BESIII. With this additional information, the knowledge of would be improved with respect to the present situation, both in precision (by about a factor two) and in robustness, due to the existence of two new independent analyses.

With the LHCb upgrade, another factor 10 in statistics will be gained and will make this extremely rare channel easier to observe. By then, the uncertainty on the will be limited by the knowledge of the normalization and the reference channels. The measurements of these channels by LHCb as well as the colliders such as Belle II and BESIII will become essential to obtain the required high precision measurement.

We have also studied a possible first observation of the decay which the theoretical model mentioned above predicts to be 40 times suppressed with respect to the decay. With the LHCb upgrade, this branching ratio could be measured with an uncertainty below 15 % via the channels. With an improved trigger efficiency, the channel could also be used to confirm an early observation in the previous mode, which would challenge the theoretical prediction.

## Acknowledgement

We would like to acknowledge B. Moussallam for collaborations at the beginning of this project and many advices. N. T. H. thanks LAL and IPNO for the hospitality and the financial support during her stay in Orsay.

## References

- [1] M. J. Savage, M. E. Luke and M. B. Wise, Phys. Lett. B 291 (1992) 481 doi:10.1016/0370-2693(92)91407-Z [hep-ph/9207233].
- [2] N. T. Huong, E. Kou and B. Moussallam, Phys. Rev. D 93 (2016) no.11, 114005 doi:10.1103/PhysRevD.93.114005 [arXiv:1511.06255 [hep-ph]].
- [3] A. Nyffeler, arXiv:1602.03398 [hep-ph].
- [4] T. P. Gorringe and D. W. Hertzog, Prog. Part. Nucl. Phys. 84 (2015) 73 doi:10.1016/j.ppnp.2015.06.001 [arXiv:1506.01465 [hep-ex]].
- [5] M. J. Ramsey-Musolf and M. B. Wise, Phys. Rev. Lett. 89 (2002) 041601 doi:10.1103/PhysRevLett.89.041601 [hep-ph/0201297].
- [6] M. Knecht, S. Peris, M. Perrottet and E. de Rafael, Phys. Rev. Lett. 83 (1999) 5230 doi:10.1103/PhysRevLett.83.5230 [hep-ph/9908283].
- [7] P. Masjuan and P. Sanchez-Puertas, arXiv:1512.09292 [hep-ph].
- [8] E. Abouzaid et al. [KTeV Collaboration], Phys. Rev. D 75 (2007) 012004 doi:10.1103/PhysRevD.75.012004 [hep-ex/0610072].
- [9] L. Cappiello, O. Cata and G. D’Ambrosio, JHEP 1304 (2013) 135 doi:10.1007/JHEP04(2013)135 [arXiv:1209.4235 [hep-ph]].
- [10] S. de Boer and G. Hiller, Phys. Rev. D 93 (2016) no.7, 074001 doi:10.1103/PhysRevD.93.074001 [arXiv:1510.00311 [hep-ph]].
- [11] C. Quigg and J. D. Jackson, UCRL-18487.
- [12] P. Vasko and J. Novotny, JHEP 1110 (2011) 122 doi:10.1007/JHEP10(2011)122 [arXiv:1106.5956 [hep-ph]].
- [13] T. Husek, K. Kampf and J. NovotnÃ½, Eur. Phys. J. C 74 (2014) 8, 3010 doi:10.1140/epjc/s10052-014-3010-4 [arXiv:1405.6927 [hep-ph]].
- [14] K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38 (2014) 090001. doi:10.1088/1674-1137/38/9/090001
- [15] R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 724 (2013) 203 doi:10.1016/j.physletb.2013.06.010 [arXiv:1304.6365 [hep-ex]].
- [16] R. Abegg et al., Phys. Rev. D 50 (1994) 92. doi:10.1103/PhysRevD.50.92
- [17] R. Aaij et al. [LHCb Collaboration], arXiv:1510.08367 [hep-ex].