Form factors and differential distributions in rare radiative leptonic B-decays
We study rare radiative leptonic decays and within relativistic quark model. In addition to previous analysis we give the estimations of the branching ratios for four values of the minimal photon energy, which correspond to photon selection criteria of the Belle and LHCb detectors. We find out that the branching ratios only slightly change. The highest values corresponding to are and . We present the distribution of the forward-backward asymmetry.
Form factors and differential distributions in rare radiative leptonic B-decays
Anastasiia Kozachukaaae-mail: firstname.lastname@example.org \@andlastauthor Nikolai Nikitin
1D. V. Skobeltsyn Institute of Nuclear Physics, M. V. Lomonosov
Moscow State University, 119991, Moscow, Russia
\@textsuperscript2M. V. Lomonosov Moscow State University, Faculty of Physics, 119991 Moscow, Russia
\@textsuperscript3A. I. Alikhanov Institute for Theoretical and Experimental Physics, 117218 Moscow, Russia
\@textsuperscript4National Research Nuclear University MEPhI, 115409 Moscow, Russia
At the quark level rare radiative leptonic decays correspond to transitions, which proceed through flavor-changing neutral currents (FCNCs) forbidden at tree level in the Standard Model (SM). These currents are described by penguin and box diagrams containing loops and thus lead to small values of the branching ratios. The typical order is – (see e.g. ). Possible contributions of new particles to the loops make these decays sensitive to New Physics. Therefore there have been a lot of studies in this sector of flavor physics and as a result at the moment we have several deviations from the SM of the order of (see discussion in [2, 3, 4]). We list some of them here:
same was observed for , in the range the discrepancy for the branching ratio is more than .
Finally, the value of the branching ratio of is 25% lower than the SM predicts, but only at . For the branching ratios of and decays the following relation takes place
were the squared ratio of masses means that the radiative decay does not have the chirality constraint, comes from the photon emission and in the denominator is the difference between three- and two-particle phase space. For muons one can easily get the estimate , which means that the branching ratios are approximately equal . In fact, is a little bit larger due to additional dynamical effects, such as resonant contributions.
The paper is organized as follows: In Section 2 we discuss the contributions to the decay amplitude . In Section 3 we calculate the transition form-factors via dispersion approach based on constituent quark picture. In Section 4 we give numerical predictions for the branching ratios, differential distributions of the decay rates and forward-backward asymmetry.
2 The effective Hamiltonian and the amplitude
where is the Fermi constant, are the scale-dependent set of Wilson coefficients, and are the basis operators. For decays the scale parameter is approximately equal to GeV. The amplitudes of the basis operators between the initial and the final states may be parameterized in terms of Lorentz-invariant form factors. The calculation of these form factors is the main challenge of this work because they contain non-perturbative QCD effects.
2.1 Direct emission of the real photon from B-meson valence quarks
The most important part which contains non-perturbative QCD contributions corresponds to the cases, when the real photon is directly emitted from the valence or quarks, and the pair is coupled to the penguin. The effective Hamiltonian in this case takes the form111Our notations and conventions are: , , , , .
and the corresponding diagrams are shown in Fig. 1. The coefficient includes long-distance effects related to resonances in the -channel [26, 27, 28]. The transition form factors of the basis operators in (2.1) are defined according to 
The penguin form factors are defined as functions of two variables: is the momentum of the photon emitted from the penguin, and is the momentum of the photon emitted from the valence quark of the meson. We calculate the form factors in the framework of the dispersion approach based on constituent quark picture, the details are presented in Section 3.
2.2 Direct emission of the virtual photon from B-meson valence quarks
Another process contributing to the amplitude is that with the real photon emitted from the penguin, whereas one of the valence quarks directly emits the virtual photon which then goes into the final pair.
This process is described by the diagrams of Fig. 2. The corresponding amplitude has the same structure as the part of the amplitude in 2.1 with replaced by . The form factors at the necessary timelike momentum transfers are not known. The difficulty is connected with appearance of neutral light vector meson resonances, and for -decays and for -decays, in the physical decay region. We calculate the form factors for with the use of gauge-invariant version [29, 30] of the vector meson dominance [31, 32, 33]
where and are the mass and the width of the vector meson resonance, are the transition form factors, defined according to the relations
Fig. 3 represents diagrams describing bremsstrahlung. The corresponding contribution to the amplitude reads
2.4 Weak annihilation contribution
3 Transition form factors
We calculate the transition form factors in the framework of relativistic quark model, which is a dispersion approach based on constituent quark picture [34, 35]. All hadron observables are given by dispersion representations in terms of hadronic relativistic wave functions and spectral densities of corresponding feynman diagrams with constituent quarks in the loops. For the wave functions we use a Gaussian parametrization . The simplest relation can be obtained for a pseudoscalar or vector meson decay constant
were is the meson relativistic wave function and is the spectral density. The latter is obtained as a direct result of feynman rules for the corresponding feynman diagram.
The example of the diagram for a B-meson decay constant is given in Fig. 5. In this work we consider meson-to-photon transitions; the corresponding form factors may be obtained in the form of the spectral representation
were and are momenta of the emitted photons. The form factors were calculated in . We now perform the calculation of the form factors . Each of these form factors contains two contributions corresponding to the cases when the photon is emitted from or quarks of the B-meson:
The spectral representations of the form factors in (3) have the form
where is the mass of the quark, which emits the photon, is the mass of the spectator, and
The model contains only a few parameters such as the constituent quark masses and the parameter of the wave function . These parameters were fixed in  using relations (9) for meson decay constants so that our results reproduce the predictions from QCD sum rules and lattice QCD.
4 Numerical results
4.1 Branching ratios
For numerical estimates we use the following values of Wilson coefficients at GeV:
, , , . The evolution including cc-resonances is taken from [26, 27, 28]. We obtained several distributions of the differential branching ratios, they are shown in Fig. 6 and 7.
|80 MeV||100 MeV||500 MeV||1000 MeV|
In Table 1 we present the results for different values of the minimal photon energy . Numerical values of are given in the -meson rest frame and are different in others because energy is not Lorentz-invariant. The values in the range correspond to different photon selection criteria at the Belle II detector , while the interval is relevant for those at the LHCb detector [38, 39]. We only take typical values as in LHCb reference frame B-mesons have a sufficiently wide energy distribution. However, from Table 1 it is clear that the branching ratios only slightly depend on the particular choice of , and the related error is smaller then that of the model. We also use the value to compare the results with those of the previous work .
4.2 Forward-backward asymmetry
We obtained distribution of the forward-backward asymmetry, defined by the relation
where , is the angle between and . The distribution is presented in Fig. 8. The measurement of the forward-backward asymmetry seems to be a hard if not impossible task, because the final state does not carry any information about the flavor of the decaying B-meson. In addition, the signs of the asymmetries corresponding to and mesons are opposite. In the absence of flavor cut the total asymmetry equals zero aside from CP-violating effects. It appears that such selection is impossible at LHCb, but at Belle II one can use the fact that neutral B-mesons are produced in an entangled state. Thus, if one of the B-mesons decays to a state with a certain flavor, the other B-meson which decays to has the opposite flavor. Now, if the interval between the decays is less than half of the of oscillation period, one can claim that the flavor of the second B-meson is known with sufficient probability. For each selection procedure one can also account for the oscillations contribution and therefore improve the prediction accuracy. A method for such calculations was developed in .
The decay rates and the forward-backward asymmetry were previously calculated in several works [23, 24, 25, 41, 43, 44]. In  and  not all the contributions were taken into account, and in [23, 24, 25, 41] the transition form factors were estimated from symmetry considerations coming from LEET. We made direct calculation of the form factors in the framework of the relativistic quark model. Our results agree nicely with . The results [23, 43] are based on not fully consistent models for the form factors and therefore do not seem to us convincing.
We obtained predictions for the differential distributions and the branching ratios for the and decays taking into account the following contributions to the amplitude of the process: photon emission from the and the valence quarks of the -meson, weak annihilation, and bremsstrahlung. The corresponding form factors were calculated in the framework of the relativistic quark model. We represented the numerical estimates for several values of the minimal photon energy corresponding to photon selection criteria of the Belle and LHCb detectors and found out that the branching ratios almost do not change with the photon energy cut.
The work was supported by grant 16-12-10280 of the Russian Science Foundation. The authors are thankful to Dmitri Melikhov for useful discussions.
-  A. Ali, hep-ph/9606324
-  S. L. Glashow, D. Guadagnoli and K. Lane, Phys. Rev. Lett. 114, 091801 (2015)
-  D. Guadagnoli and K. Lane, Phys. Lett. B 751, 54 (2015)
-  D. Guadagnoli, D. Melikhov and M. Reboud, Phys. Lett. B 760, 442 (2016)
-  R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 113, 151601 (2014)
-  C. Bobeth, G. Hiller and G. Piranishvili, JHEP 0712, 040 (2007)
-  C. Bouchard et al. [HPQCD Collaboration], Phys. Rev. Lett. 111, no. 16, 162002 (2013) Erratum: [Phys. Rev. Lett. 112, no. 14, 149902 (2014)]
-  G. Hiller and F. Kruger, Phys. Rev. D 69, 074020 (2004)
-  R. Aaij et al. [LHCb Collaboration], JHEP 1406, 133 (2014)
-  R. Aaij et al. [LHCb Collaboration], JHEP 1302, 105 (2013)
-  C. Bobeth, G. Hiller and D. van Dyk, JHEP 1107, 067 (2011)
-  C. Bobeth, G. Hiller, D. van Dyk and C. Wacker, JHEP 1201, 107 (2012)
-  C. Bobeth, G. Hiller and D. van Dyk, Phys. Rev. D 87, no. 3, 034016 (2013)
-  R. Aaij et al. [LHCb Collaboration], JHEP 1509, 179 (2015)
-  R. Aaij et al. [LHCb Collaboration], JHEP 1602, 104 (2016)
-  A. Abdesselam et al. [Belle Collaboration], arXiv:1604.04042 [hep-ex].
-  J. P. Lees et al. [BaBar Collaboration], Phys. Rev. Lett. 109, 101802 (2012)
-  R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 115, no. 11, 111803 (2015) Addendum: [Phys. Rev. Lett. 115, no. 15, 159901 (2015)]
-  M. Huschle et al. [Belle Collaboration], Phys. Rev. D 92, no. 7, 072014 (2015)
-  V. Khachatryan et al. [CMS and LHCb Collaborations], Nature 522, 68 (2015)
-  B. Grinstein, M. J. Savage and M. B. Wise, Nucl. Phys. B 319 271 (1989)
-  A. J. Buras and M. Munz, Phys. Rev. D 52, 186 (1995)
-  C. Q. Geng, C. C. Lih and W. M. Zhang, Phys. Rev. D 62, 074017 (2000)
-  Y. Dincer and L. M. Sehgal, Phys. Lett. B 521, 7 (2001)
-  F. Kruger and D. Melikhov, Phys. Rev. D 67, 034002 (2003)
-  F. Kruger and L. M. Sehgal, Phys. Rev. D 55, 2799 (1997)
-  D. Melikhov, N. Nikitin and S. Simula, Phys. Lett. B 430, 332 (1998)
-  D. Melikhov, N. Nikitin and S. Simula, Phys. Rev. D 57, 6814 (1998)
-  D. Melikhov, Phys. Lett. B 516, 61 (2001)
-  D. Melikhov, O. Nachtmann, V. Nikonov and T. Paulus, Eur. Phys. J. C 34, 345 (2004)
-  J. J. Sakurai, Annals Phys. 11, 1 (1960)
-  M. Gell-Mann and F. Zachariasen, Phys. Rev. 124, 953 (1961)
-  G. J. Gounaris and J. J. Sakurai, Phys. Rev. Lett. 21, 244 (1968)
-  D. Melikhov, Phys. Rev. D 53, 2460 (1996)
-  D. Melikhov, Phys. Rev. D 56, 7089 (1997)
-  M. Neubert and B. Stech, Adv. Ser. Direct. High Energy Phys. 15, 294 (1998)
-  A. Kozachuk, D. Melikhov and N. Nikitin, Phys. Rev. D 93, no. 1, 014015 (2016)
-  [LHCb Collaboration], LHCb calorimeters: Technical design report (CERN, Geneva, 2000)
-  R. Voss and A. Breskin, The CERN Large Hadron Collider, accelerator and experiments (CERN, Geneva, 2009)
-  T. Aushev et al., arXiv:1002.5012 [hep-ex]
-  D. Melikhov and N. Nikitin, Phys. Rev. D 70, 114028 (2004)
-  C. Patrignani et al. (Particle Data Group), Chin. Phys. C, 40, 100001 (2016)
-  W. Y. Wang, Z. H. Xiong and S. H. Zhou, Chin. Phys. Lett. 30, 111202 (2013)
-  I. Balakireva, D. Melikhov, N. Nikitin and D. Tlisov, Phys. Rev. D 81, 054024 (2010)