#### Abstract

The charmless decays (where and denote the light pseudoscalar and vector mesons, resp.) can occur only via the weak annihilation diagrams within the Standard Model. In this paper, we study these kinds of decays in the framework of QCD factorization, by adopting two different schemes: scheme I is similar to the method usually adopted in the QCD factorization approach, while scheme II is based on the infrared behavior of gluon propagator and running coupling. For comparison, in our calculation, we adopt three kinds of wave functions for meson. The branching ratios based on the two schemes are given. It is found that (a) the predicted branching ratios in scheme I are, however, quite small and almost impossible to be measured at the LHCb experiment and (b) in scheme II, by assigning a dynamical gluon mass to the gluon propagator, we can avoid enhancements of the contribution from soft endpoint region. The strength of annihilation contributions predicted in scheme II is enhanced compared to that obtained in scheme I.

#### 1. Introduction

The meson is the lowest-lying bound state of two heavy quarks with different flavors ( and ). Due to its flavor quantum numbers and being below the threshold, the meson is stable against strong and electromagnetic interactions and can decay only via weak interaction. Furthermore, the meson has a sufficiently large mass; each of the two heavy quarks can decay individually, resulting in rich decay channels [1]. Therefore, the meson is an ideal system to study weak decays of heavy mesons [2].

The experimental studies of meson properties started in 1998 when the Collider Detector at Fermilab (CDF) reported the first observation of meson through the semileptonic decay modes [3]. Thanks to the fruitful performance of the CDF, D0, and LHCb collaborations, both the mass [4–6] and the lifetime [7–9] of the meson have been measured quite accurately. At the Large Hadron Collider (LHC) with a luminosity of about , one could expect around events per year [10]. In addition, several hadronic decay channels, such as [11] and [12], have also been observed for the first time. In the following years, the properties of meson and the dynamics involved in decays will be further exploited through the precision measurements at the LHC with its high collision energy and high luminosity, opening therefore a golden era of physics [13].

The theoretical investigations have also been carried out on the properties of meson, such as its lifetime, its decay constant, and some of its form factors, based on different theoretical frameworks [2]. Due to its heavy-heavy nature and the participation of strong interaction, the hadronic decays are extremely complicated but, at the same time, provide great opportunities to study the perturbative and nonperturbative QCD and final-state interactions in heavy meson decays. Being weakly decaying and doubly heavy flavor meson, it also offers a novel window for studying the heavy-quark dynamics that is inaccessible through the and quarkonia [2]. These features have motivated an extensive study of decays in various theoretical approaches in the literature [14].

In this paper, we will focus on the two-body charmless hadronic decays, which can proceed only via the weak annihilation diagrams in the Standard Model (SM): the initial and quarks annihilate into and quarks, which form two light mesons by hadronizing with () pair emitted from a gluon. Detailed studies of these decays will be certainly helpful for further improving our understanding of the weak annihilation contributions, the size of which is currently an important issue in physics.

The recent measurements of the decays, especially of the pure annihilation processes and [15–17], indicate that the annihilation topologies can be significant, contrary to the common belief of their power suppression in the heavy-quark limit [18]. Although it was later noticed theoretically that the annihilation amplitudes may not be negligibly small in realistic meson decays [19], it is still very hard to make a reliable calculation of these diagrams, and quantitative predictions for them vary greatly between different approaches. In the QCD factorization (QCDF) approach [20], they can only be estimated in a model-dependent way due to the endpoint singularities [21]. In the soft-collinear effective theory (SCET) [22], they are argued to be factorizable and almost real with tiny strong phase [23], which is rather different from almost imaginary with large strong phase as predicted in the perturbative QCD approach (pQCD) [19]. In addition, the annihilation contributions in many decays usually involve both tree and penguin operators, and they interfere with many other different topologies, making it difficult to obtain an accurate value of annihilation by fitting the experimental data [24].

The charmless decays into two light mesons, coming only from a single tree operator, provide therefore an ideal testing ground for annihilation in heavy meson decays and deserve detailed studies using different theoretical approaches [25–28]. In this paper, we will revisit these decays in the QCDF framework, using two different schemes proposed to deal with the endpoint singularity and to avoid enhancements in the soft endpoint region: the divergence in scheme I is usually parameterized with at least two phenomenological parameters through the treatment [21], whereas in scheme II, one could use an infrared finite gluon propagator [29], to regulate the divergent integrals [30–34]. The different scenarios corresponding to different choices of and in scheme I have been thoroughly discussed in [21]. In scheme II, it is found that the hard spectator-scattering contributions are real and the annihilation corrections are complex with a large imaginary part [30, 31, 34]. These two different treatments used in decays could be further tested through the charmless decays.

The remaining parts of the paper are organized as follows. In Section 2, after recapitulating the theoretical framework for two-body charmless hadronic decays, we present the calculation of the annihilation diagrams in the QCDF framework with the two different schemes. The numerical results and discussions are given in Section 3. Finally, Section 4 contains the main conclusions and a short summary. The explicit expressions for the decay amplitudes and the relevant input parameters are collected in Appendices A and B, respectively.

#### 2. Theoretical Framework and Calculation

##### 2.1. The Effective Weak Hamiltonian and Hadronic Matrix Element

Using the operator product expansion and renormalization group (RG) equation, we can write the effective weak Hamiltonian for charmless ( denote the light pseudoscalar and vector mesons) decays as [35]where is the Fermi coupling constant and and () are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [36]. The four-quark operators arise from W-boson exchange and are defined, respectively, aswhere , are the color indices. The corresponding Wilson coefficients can be calculated using the RG improved perturbative theory [35].

To obtain the decay amplitude, the remaining work is to evaluate the hadronic matrix elements of the local operators , which is however quite difficult due to the participation of nonperturbative QCD effects. The Feynman diagrams for decays with the QCDF approach are shown in Figure 1, where (a), (b) and (c), (d) are nonfactorizable and factorizable topologies, respectively. Since the tree operators have the Dirac structure, the two factorizable diagrams (c) and (d) cancel each other exactly in the QCDF approach due to the conservation of the vector current and partial conservation of axial-vector current, plus the approximation that final states have the same asymptotic expressions for twist-2 and twist-3 distribution amplitudes (see (12)) in our calculation [21]. Moreover, due to the mismatch of the color indices, there are no contributions from diagrams (a) and (b) with the insertion of the color-singlet operator . Thus, there is only a single tree operator involved in the decay amplitudes, and the nonzero contribution comes only from diagrams (a) and (b).

**(a)**

**(b)**

**(c)**

**(d)**

In the QCDF framework and with the same hypotheses made for hadronic decays, the decay amplitude for charmless decays can be written as [21]where , are decay constants of the and mesons, respectively. The coefficient is defined as [21]where and , the superscript “” on refers to the gluon emission from the initial-state quarks, and the subscript “” on refers to the Dirac structure of the inserted four-quark operator . The basic building block can be expressed as the convolution of the hard kernels given by diagrams (a) and (b) in Figure 1 and the light-cone distribution amplitudes (LCDAs) of the initial- and final-state mesons, which will be detailed in the next two subsections.

##### 2.2. in Scheme I

In scheme I, the annihilation contributions to hadronic decays are evaluated by regularizing the divergent integrals on the basis of heavy-quark power counting [21]. Despite the fact that such a treatment is not entirely self-consistent in the context of a hard-scattering approach, it provides nevertheless a model to estimate the importance of annihilation, which, motivated by the first observation of the pure annihilation decay [15–17], has been revisited quite recently in [31, 37].

Following a similar treatment, we now estimate the annihilation topologies in charmless decays. In accordance with the convention adopted in [21], we find that the basic building block is given bywhen both mesons are pseudoscalar or when is a pseudoscalar and a vector meson. In the case when is a vector meson and a pseudoscalar, one has to change the sign of the second term in . When we take and , this result is in agreement with the expressions obtained in [21, 38]. In (5), and denote the relative size of the and quark masses withTheir appearance allows one to distinguish the origin of each term in the brackets: the ones involving must come from diagram (a), whereas those involving must come from diagram (b) in Figure 1. As always, and denote the leading-twist and twist-3 two-particle LCDAs of the final-state meson , respectively. The factor , once multiplied by , is used to normalize the twist-3 distribution amplitude; explicitly, we havewhere denote the running masses of the two valence quarks of a pseudoscalar, and is the scale-dependent transverse decay constant of a vector meson. Despite being formally suppressed by one power of in the heavy-quark limit, these terms are not always small numerically, especially in the case of pseudoscalar mesons [21].

In the calculation, we use three different types of distribution function for meson. The first one is the peak form (W-I) [39]The second one is the solution of the equation with the harmonic oscillator potential (W-II) [40]where , the reduced mass , and the quantum of energy GeV [41]. The third one is the quarkonium form (W-III) [42]where .

In (9) and (10), is normalization constant and the normalization condition isThe shape of the three distribution functions for meson are displayed in Figure 2.

To pursue the structure of the singularities of the building block , we take, for simplicity, the asymptotic expressions for the distribution amplitudes [21, 43]

The weak annihilation of exhibits endpoint singularities even at twist-2 order in the light-cone expansion for the final-state mesons. For decays, the endpoint singularities only at twist-3 level, the situation is the same as the hard spectator interactions of . For the twist-2 terms, the singularities are in the integration interval. It is found that the convolution integrals in (5) can be performed without problem as long as . There are, however, integrable singularities at or when , which can be dealt with using the prescription of Cauchy principal value integral. Taking the integralas an example, we show in Figure 3 its real and imaginary parts dependence on the parameter , and one can see clearly that the integral is finite as long as is different from and . In (5), the twist-3 terms are more complex and can not be expressed as polynomial of , so we make the integral interval of .

Rather than giving the explicit expressions for the convolution integrals, we present, with the default inputs GeV and GeV, the numerical results for the building block in the three different cases

*W-I*

*W-II*

*W-III*where the results are obtained with and GeV. Judging from the above expressions, the branching ratios obtained with W-I and W-III should be very close, and the W-II’s results will be smaller. It is noted that is identical to in our approximation and the annihilation contribution has a large imaginary part.

##### 2.3. in Scheme II

Instead of being parameterized with an* ad hoc* model-dependent cut-off, the endpoint divergences can also be regulated with an infrared (IR) finite gluon propagator that is characterized by a dynamical gluon mass, providing therefore a natural IR regulator [29]. This has been successfully applied to various hadronic decays in [30–34]. In this subsection, we will evaluate the building block in this scheme.

Instead of the perturbative expression that is IR divergent, the IR finite gluon propagator is obtained by solving an intricate set of coupled Dyson-Schwinger equations (DSE) for pure gauge QCD, under a systematic approximation and truncation [29]. It is also noted that any IR finite gluon propagator leads to a freezing of the IR coupling constant [44], meaning that the use of an IR finite gluon propagator must be accompanied by an IR finite coupling constant. The above information about the IR behavior of QCD has also been confirmed by the most recent lattice simulations [45, 46] (recent reviews, together with a list of references, on DSE solutions and lattice results about the infrared finite gluon propagator and running coupling constant could be found, e.g., in [47–50]). Here we adopt the gluon propagator derived by Cornwall many years ago [29]:where denotes the gluon momentum squared. The corresponding running coupling constant reads [29]where is the first coefficient of the QCD beta function, and the number of active quark flavors at a given scale. The dynamical gluon mass is given by [29]where MeV is the QCD scale, and the effective gluon mass with a typical value GeV [29]. It is interesting to note that similar values are found by fitting the experimental data on decays: GeV from decays [30] and GeV from decays [31]. In our calculation, we take GeV. As shown in Figure 3, both the gluon propagator (23) and the coupling constant (24) are IR finite and different from zero at the origin of momentum squared .

With the above prescription and the same convention used in scheme I, our final results for the building block can be expressed as ()when and . If , on the other hand, the sign of the second term in has to be changed. When taking the limits and , our results agree with the ones given in [30, 31, 34].

In (26), the time-like gluon momentum squared is given by and also depends on the longitudinal momentum fractions and , making the convolution integrals rather complicated. As shown in Figure 5, although the running coupling constant is rather large in the small region (see Figure 4(b)), the fact that only a small fraction comes from the region in the plane indicates that the annihilation contributions are still dominated by the region associated with a large imaginary part [30]. In scheme II, we also use three kinds of meson wave function in the calculation. With the default inputs GeV, GeV, GeV, and GeV, our numerical results for the building block read

**(a)**

**(b)**

*W-I*

*W-II*

*W-III*It is found that, compared to the values obtained in scheme I ((16)–(22)), the annihilation contributions predicted in scheme II are enhanced. These will apparently affect the predictions for charmless hadronic decays, which will be detailed in the next section.

#### 3. Numerical Results and Discussions

In the meson rest frame, the branching ratio for general charmless decays can be written as where is the meson lifetime and is the center-of-mass momentum of either of the two outgoing mesons, withThe decay amplitude can be obtained from the hadronic matrix element defined in (3); for convenience, we collect in Appendix A the explicit expressions of the decay amplitudes for the considered decay modes. The CP-violating asymmetries for all the considered decays are absent, since there involves only a single tree operator in the decay amplitudes, which can be clearly seen from (3).

With the theoretical expressions given above and the input parameters collected in Appendix B, we proceed to evaluate the CP-averaged branching ratios for these charmless decays. In our calculation, the default value of the renormalization scale is set at , which is approximately the averaged virtuality of the time-like gluon propagated in the annihilation diagrams. The numerical results based on the two schemes are collected in Table 2. In Tables 3, 4, and 5, we also present detailed error estimates induced by the theoretical uncertainties of input parameters. The first error shown corresponds to the variation of the CKM parameters and (named as “CKM”); the second error refers to the variation of the quark masses, decay constants, and the mixing angle (named as “hadronic”). The third error arises from the variation of the renormalization scale (named as “scale”). The last error reflects the uncertainty due to the dynamical gluon mass (named as “”).

Based on the results collected in Tables 2–5, we have the following observations and remarks:(i)The two-body charmless hadronic decays can be classified into two categories: the strangeness-conserving () and the strangeness-changing () processes. From the numerical results listed in Table 2, one can see that the branching rations of channels are generally much larger than those of ones. This is due to the large hierarchical structure between the two CKM matrix elements and , .(ii)In scheme I, the branching ratios obtained with W-I and W-III are very close, which vary in the ranges of to , being larger than the corresponding ones obtained with W-II. This is consistent with the wave functions; the sharp of W-III is very close to function as shown in Figure 2. In scheme II, the annihilation contributions are enhanced when we adopt the IR finite gluon propagator, and the branching ratios are not sensitive to the choice of wave function for meson.(iii)Among these charmless decays, only several decays modes, such as , , , , and , have relatively large branching ratios, being around in scheme II. All of these channels belong to the transitions that are CKM favored. It is found that branching ratio for decay is relatively large among decays into final states and , which is promisingly detected by experiments at the running Large Hadron Collider and forthcoming SuperKEKB.(iv)For , , and decays, on the other hand, since , the contributions from and components of the neutral mesons cancel each other exactly or almost, resulting in (approximate) zero branching ratios of these three channels. For , , and decays, due to the flavor decomposition , the interference between the two flavor components and of the neutral mesons is constructive, resulting in larger branching ratios. Taking into account the fact that , one can easily understand the pattern of their branching ratios. In particular, the decay modes and have similar branching ratios, because only the term is involved in the decay amplitudes.(v)For decays, the obtained branching ratios show a rather different pattern, , from that of . It is also observed that is much larger than , while is suppressed rather than enhanced compared to that of the mode. To understand the enhancement and suppression patterns, we should note that both the and terms contribute to these transitions, but with an opposite sign between them for the and final states, which is due to the fact that , , while . This results in a destructive interference for , but a constructive interference for the modes. Similar patterns have already been observed in the and decays [51].(vi)As discussed in [26], several relations among the charmless decay channels hold in the limit of exact SU() flavor symmetry. For decays, for example, one of such relations reads with the Cabibbo-suppressing factor . Similar relations could also be found for decays, with the replacements and/or . We find that the first equality holds exactly in both scheme I and scheme II, because the exact isospin symmetry is assumed in our calculation. The second equality is, however, violated by the differences between decay constants and light-quark masses, which account for the SU() breaking effect.(vii)As a comparison, the pQCD predictions [25] are shown in the last column of Table 2. On the whole, the central values for the branching fractions obtained in this paper are smaller than those obtained in the pQCD approach. The different choice of the renormalization scale may be the main reason leading to these discrepancies. There are two points to notice about the results: () As mentioned earlier, for , decays, their branching ratios are zero or almost zero in our result. But their branching ratios are not very small in the pQCD method, as shown in Table 2. () For , decays, their branching ratios are almost the same in our results. And in the pQCD approach, . The possible causes of the two differences are the same: In the pQCD approach, they use the LCDAs in the form of Gegenbauer polynomials which contain the SU() breaking effect.(viii)As the relevant CKM parameters have been measured quite precisely, the theoretical errors introduced by the CKM parameters are small. The uncertainty due to the variation of the dynamical gluon mass in scheme II is also found to be negligible. The main uncertainties are due to the variation of the renormalization scale, as well as the mixing parameters for and final states. The large scale dependence of the branching ratios is understandable, because only the leading order term in is taken into account in our calculation. Furthermore, the different choices of the renormalization scale also account for the main differences among our results and the ones presented in [25, 26].

Finally, we would like to point out that it is hard to estimate the systematical uncertainties coming from the hypothesis underlying our calculations, such as the one-gluon approximation for the annihilation mechanism, the use of asymptotic distribution amplitudes, and the neglect of -suppressed power corrections.

#### 4. Summary

Being the lowest-lying bound state of two heavy quarks with different flavors, the meson is an ideal system to study weak decays of heavy mesons. In this paper, we have carried out a detailed study of two-body charmless hadronic decays, which can proceed only via the weak annihilation diagrams within the SM and are, therefore, very suitable for further improving our understanding of the annihilation mechanism, the size of which is currently an important issue in physics. Explicitly, we have adopted two different schemes to deal with these decays: scheme I is similar to the usual method adopted in the QCDF approach, while scheme II is based on the infrared behavior of gluon propagator and running coupling. For comparison, we adopt three different kinds of distribution function for meson in our calculation. It is found that the strength of annihilation contributions predicted in scheme II is enhanced compared to that obtained in scheme I. The branching ratios are not sensitive to the choice of wave function for meson in scheme II. However, the predicted branching ratios are inconsistent with the corresponding ones obtained in the pQCD approach [25].

The large discrepancies among these theoretical predictions make it very necessary to make more detailed studies of these kinds of decays, especially from the experimental side. It is interesting to note that the LHCb experiment has the potential to observe the decays with a branching ratio of , which will certainly provide substantial information on these charmless decays and deepen our understanding of the annihilation mechanisms.

#### Appendix

#### A. Decay Amplitudes in the QCDF Approach

Starting with (3) and adopting the standard phase convention for the flavor wave functions of light and heavy mesons [21, 51, 56], one can easily write down the decay amplitude for a given decay mode. Firstly, there are eight charmless decays with the corresponding amplitude given, respectively, as (the exact isospin symmetry is assumed):The decay amplitudes for the charmless modes can be written, respectively, as:

#### B. Input Parameters

To get the Wilson coefficients at the lower scale , we adopt the following input parameters [54]:We also vary the renormalization scale in the region to assess the scale uncertainty.

For the CKM matrix elements, we use the Wolfenstein parameterization [57] and keep terms up to [35]:with the inputs and [58].

For the system, we adopt the Feldmann-Kroll-Stech (FKS) mixing scheme defined in the quark-flavor basis [56], where the physical states and are related to the flavor states and byThe decay constants and , as well as the other hadronic parameters related to and can then be expressed in terms of two decay constants and the mixing angle [51]. The values of these three parameters have been determined from a fit to experimental data, yielding [56]

Finally, a summary of the other input parameters entering our numerical analysis is given in Table 1. It is noted that the latest experimental determinations of and [54] compare positively within errors with the lattice results [59]. Our values of the vector-meson decay constants are taken from [52], which are an update of the ones extracted in [53]. The scale dependence of the transverse decay constants is taken into account via the leading-logarithmic running . The light quark masses given in the table are the running masses defined in the scheme; to get the corresponding pole and running masses at different scales, we use the NLO running formulae collected, for example, in [60]. The and quark masses are, however, defined as the pole masses.

#### Competing Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by CCNU-QLPL Innovation Fund (QLPL201411).