Quasi-two-body decays in the perturbative QCD approach

In this study, we investigate quasi-two-body decays in the perturbative QCD approach. Two-meson distribution amplitudes are introduced to describe the final state interactions of the pair, which involve time-like form factors and Gegenbauer polynomials. We calculate the CP averaged branching ratios of the decays. Our results are in agreement with newly updated data measured by Belle II. This suggests that it is more appropriate to analyze these quasi-two-body B decays in the three-body framework than the two-body framework. We also predict direct CP asymmetries for the considered decay modes and find that is small and less than % in magnitude, whereas is larger and can reach a few percent. Our predictions can be tested in future B meson experiments.


I. INTRODUCTION qqqq qqqqq
Many experimental studies on three-body B meson decays [1][2][3][4][5][6][7][8][9] have been performed in recent years. This type of decay has gained increasing attention for the following reasons. (1) Many new resonance states are observed in the invariant mass distributions of three-body decays, which are difficult to understand in terms of a common meson or baryon and are known as exotic states. Researchers are puzzled by their inner structures and have proposed numerous assumptions, such as compact tetraquarks, pentaquarks, loosely bound hadronic molecules, and glueball and hybrid states. (2) These decays involve significantly more complicated QCD dynamics compared with two-body decay cases, which imposes a serious challenge to present theoretical frameworks. The hard b-quark decay kernels in three-body B decays contain two virtual gluons at the leading order, which are difficult to directly evaluate owing to the enormous number of Feynman diagrams. (3) Large direct CP asymmetries in the localized regions of phase spaces for three-body B decays have been observed in experiments. For three-body B decays, CP violation depends on the invariant mass, and regional CP asymmetries can be SU (3) measured through the Dalitz plot. However, such observables do not exist in two-body B decays. To study these three-body B decays, many approaches based on symmetry principles and factorization theorems have been proposed. Symmetry principles include U-spin [10][11][12][13], flavor symmetry [14][15][16][17], and the factorization-assisted topological-diagram amplitude (FAT) approach [18]. Factorization theorems include the QCD-improved factorization approach [19][20][21][22][23][24] and perturbative QCD (PQCD) approach [25][26][27][28][29][30][31][32][33][34], and it has been proposed that the factorization theorem of three-body B decays is approximately valid when two energetic final particles move collinearly and the bachelor particle recoils back. At the same time, the interactions between them can be neglected. Based on this quasi-two-body-decay approximation, two-meson distribution amplitudes (DAs) are introduced into the PQCD approach, where the strong dynamics between the two final hadrons in the resonant regions are included.
On the experimental side, the decays have recently been investigated by Belle II [35], and their branching ratios were obtained through different decay modes. (1) The two-body decays have been studied using different theories, including the PQCD approach [36][37][38]. Here, we study three-body radiative B decays in the quasi-two-body approximation using the PQCD approach. Verifying the consistency between the three-body and two-body analyses of the decays can support PQCD factorization for exclusive B meson decays. Furthermore, by studying these more subtle channels compared with two-body decays, we can deepen our understanding of isospin symmetry and the narrow width approximation involved in these decays. After introducing the new non-perturbative inputs, the two-meson DAs, the factorization formula for the three-body B decay can be written as [39,40] where denotes the meson DA, is the two-meson DA for that fly collinearly, and represents the convolution in parton momenta. Then, the hard kernel H for the b quark decay, similar to the two-body decay case, starts with the diagrams of single hard gluon exchange.

Kπ
This paper is organized as follows. In Sec. II, the kinematic variables for three-body radiative B decays are defined. The considered two-meson ( ) P-wave DAs are parametrized, whose normalization form factors are assumed to follow the relativistic Breit-Wigner (RBW) model. Then, Feynman diagrams and the total decay amplitudes for these channels are given. In Sec. III, the numerical results are presented and discussed, where we compare our predictions with other theoretical and experimental results. The summary is presented in the final section.

II. FRAMEWORK
We begin with the parametrization of the kinematic variables involved in the decays. In the rest frame of the B meson, we define the B meson momentum , K meson momentum , π meson momentum , meson momentum , and γ momentum in light-cone coordinates as is the B meson mass, the variable , ω is the invariant mass of the pair, and ζ is the momentum fraction for the K meson. The momenta of the light quarks in the B and mesons are denoted as and k, respectively, where and z are the momentum fractions.

Kπ
Because the massless photon is only transversely polarized, only the transversely polarized P-wave DAs of the pair contribute to the final results, which are defined as [30] Φ T with the functions [41] ϕ T Kπ (z, where , and the Gegenbauer moments associated with the transverse polarization , are determined in Ref. [34] and are listed in the next section. The strong interactions between the resonance and final-state meson pair can be factorized into a time-like form factor, which is guaranteed by the Watson theorem [42]. For narrow resonances, the RBW [43] function is a convenient model to effectively separate them from other resonant and nonresonant contributions with the same spin and has been widely used in experimental data analyses. Here, the time-like form factor is parametrized with the RBW line shape and can be expressed in the following form [44,45]: where and are the pole mass and width, respectively. The mass dependent width is defined as where is the magnitude of the momentum measured in the resonance rest frame, whereas is the value of corresponding to . is the orbital angular momentum in the system, and corresponds to the P-wave resonances. Owing to limited studies on the form factor , we use the two decay constants and of the intermediate particle to determine this factor through the ratio [27]. Note that this approximation is workable in the region with .
For the wave function of the heavy meson [46], we take Here, only the contribution of the Lorentz structure is considered because the contribution of the second Lorentz structure is numerically small and has previously been neglected [47]. For the distribution amplitude in Eq. (9), we adopt the following model: where the shape parameter GeV is well fixed using rich experimental data on the meson from many studies, and the coefficient is determined by the normalization . An updated discussion on the B meson wave function can be found in Ref. [48].

B. Analytic formulae
For quasi-two-body decays, the effective Hamiltonian relevant to the transition is given by [49] where the Fermi coupling constant [50], and and are the products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The scale μ separates the effective Hamiltonian into two distinct parts: the Wilson coefficients and local fourquark operators . The local four-quark operators are written as with the color indices i and j. Here, refer to the Lorentz structures . Note that the terms associated with the strange quark mass in the and operators are dropped.
The typical Feynman diagrams at the leading order for quasi-two-body decays are shown in Figs. 1-5, where we take the decay as an example. The contributions from the operators and the annihilation type diagrams are presented, and the analytic formulas for the decay amplitudes of each Feynman diagram can be found in our previous paper [32]. The wave functions and corresponding parameters must be replaced in the calculations.
Quasi-two-body decays in perturbative QCD approach Chin. Phys. C 46, 123105 (2022) By combining the amplitudes from the different Feynman diagrams, the total decay amplitude of the charged B meson decay is given as where correspond to contributions from the righthanded and left-handed photons, respectively, and the combinations of the Wilson coefficients are defined as Similarly, the total decay amplitudes for the decays are listed below.
Then, the differential decay rate can be described as where the squared amplitudes for the meson decays are summed in the helicity basis, is the mean lifetime of the meson, and the kinematic variables and denote the magnitudes of the K and γ momenta in the center-of-mass frame of the pair,

III. NUMERICAL RESULTS
The adopted input parameters in our numerical calculations are summarized below (the masses, decay constants, Gegenbauer moments, and QCD scale are in units As for the CKM matrix elements, we employ Wolfenstein parametrization with the inputs [50] λ = 0.22453 ± 0.00044, A = 0.836 ± 0.015, Kπ Using the differential branching ratio in Eq. (17) and the squared amplitudes in Eqs. (13) and (15) and integrating over the full invariant mass region for the resonant components, we obtain the branching ratios of the quasi-two-body decays as where the first source of errors originates from the shape parameter of the B meson DA, GeV, the second error is from the Gegenbauer coefficients in the kaon-pion DAs, GeV, GeV, and the final two errors are induced by nextto-leading-order effects in the PQCD approach, that is, changing the hard scale t from 0.80t to 1.2t and the QCD scale GeV. From our results, we can see that the dominant theoretical error arises from the uncertainty on , which is close to %. The error induced by the Gegenbauer coefficients in the pair DAs is smaller at approximately %. These four decays are mediated by the transition. In our calculations, the Feynman diagrams from the operator give the dominant contribution, which is proportional to .
Assuming isospin conservation for the strong decay , we can obtain the following relations: Under the narrow width approximation, the branching ratios of these quasi-two-body decays can be expressed as Using the experimental data given in Eq. (1), which were measured by Belle II in the last year, combined with isospin conservation (Eq. (22)) and the narrow width approximation (Eq. (23)), we can estimate the branching ratios of the quasi-two-body decays as follows: These estimates and our predictions are consistent with each other; therefore, it is reasonable to extend the PQCD approach to quasi-two-body B meson decays. Using the numerical results given with Eq. (21), we calculate the relative ratios between the branching ratios of the charged and neutral B decays, If we assume the branching ratio of the decay to be % and the narrow width approximation to be true, we can relate them with the ratio . Using data from PDG [50], the ratio can be obtained. Using the updated data measured by Belle II [35], the ratio . The values of theseratios once again support the usability and rationality of PQCD factorization for quasi-two-body B decays. Furthermore, we can relate these ratios with isospin asymmetry, which is defined as This can be rewritten using the formula The ratios suggest that isospin asymmetry between the decays and is small at approximately % in magnitude, which is supported by the present data. The branching ratios of the two-body decays have previously been calculated with the PQCD approach [37], where the results were . Two years later, they were updated to [53]. Comparing these two calculations, differences still exist. A comparison of all the theoretical results with the updated data measured by Belle II supports the idea that it is more appropriate to study these quasitwo-body B decays in the three-body framework than the two-body framework. Under the times of high precision measurement, these results should be carefully tested in LHCb experiments.

Kπ
B 0 s → Using the same two-meson DAs for the pair, we also calculate the branching ratios for the )γ decays and obtain the results as These two decays are induced by the transition, whose penguin contributions are proportional to and expected to be suppressed by one order of magnitude relative to those induced by the transition. From Eq. (28), we can obtain through the internal decay mode and through . Our predictions are consistent with the results given in Ref. [53] Br The consistency of these results indicates that the PQCD approach is applicable to two-body and three-body B decays at the same time. 76 We also predict the dependences of the and decay spectra shown in Fig. 6, which exhibit a maximum at the invariant mass at approximately 0.895 GeV. The curves for the other four B decay modes are similar because the same time-like form factors are used for the DAs. It is easy to see that the main contribution to the branching ratio originates from the region around the pole mass of the resonance, as expected. For example, the central values of the branching ratio are and when we integrate over ω by limiting the ranges of and , repectively. These amount to % and % of the total branching ratio , respectively.
The direct CP asymmetry of is defined as For our considered decays, direct CP asymmetry is mainly induced by interference between the contributions from the and tree operators, which is proportional to . We predict the direct CP asymmetries as where we find that the direct CP violations of the decays induced by the transition are less than %. For decays, our predictions are consistent with those of the two-body decays predicted using the PQCD approach, and [37], and recalculated to be and in Ref. [53]. On the experimental side, the direct CP asymmetries of the decays are given by the PDG [50] A CP ( where is larger than % and opposite in sign to the PQCD predictions. Nevertheless, large errors still exist. We hope that this divergence can be clarified in future LHCb and SuperKEKB experiments by measuring the three-body B decays and . If direct CP violation of more than a few percent is confirmed in the future, it can be considered that some new physics may contribute to these channels. Direct CP violation for decays induced by the transition are significantly larger than those of decays because the product of the CKM matrix ele- ment for the electro-magnetic penguin operator is proportional to , and that for the tree operator is either proportional to or . In other words, the tree contribution is not suppressed and can be comparative with the penguin contribution. Because we know that direct CP violation arises from the interference between the tree and penguin contributions, we can expect relatively large CP asymmetries for these two decays. In previous PQCD calculations [53], the authors obtained the following result: which is indeed significantly larger than those of decays. Though this result is larger than our predictions, it is clear that direct CP violation for the decay has a positive sign, which contrasts with those of decays.

IV. SUMMARY
In this study, we analyze the three-body radiative B decays with the pair originating from the intermediate state using the PQCD approach. Under the quasi-two-body-decay approximation, the pair DAs are introduced, which include final-state interactions between the pair in the resonant region. Both the resonant and nonresonant contributions are described by the time-like form factor , which is parametrized using the RBW formula for the P-wave resonance . Under the condition of the narrow width approximation and isospin conservation, the branching ratios of the decays are consistent with those of two-body decays calculated in previous PQCD calculations, which verifies that the PQCD approach can be extended to three-body B decays. More importantly, our predictions are closer to data recently measured by Belle II. This indicates that it is more appropriate to study quasi-two-body B decays in the three-body framework than the two-body framework. For the decays induced by the transition, their direct CP violations are small and less than %. If direct CP violation of more than a few percent is confirmed in the future, it can be considered that some new physics may contribute to these channels. For the decays induced by the transition, stronger interference exists between the tree and penguin contributions; therefore, relatively large CP asymmetries can be observed, which can be tested in the LHCb and Belle II experiments.