Abstract

Exclusive decay is studied in the framework of the three-point QCD sum rules approach. The two gluon condensate contributions to the correlation function are calculated and the form factors of this transition are found. The decay width and total branching ratio for this decay are also calculated.

1. Introduction

The standard model (SM) Higgs boson which is one of the most important components of the SM has been discovered by the ATLAS [1] and CMS [2] collaborations. Nowadays, we aim to find out the new physics beyond the SM. Heavy mesons with the different flavors like and mesons can provide a good testing benchmark not only for the predictions of the SM but also for searching the new physics beyond SM. The LHCb experiment has aimed to test the SM predictions and discover the possible new physics signals. In this regard, a lot of the experimental data are released by the LHCb experiment [3].

The dominant decay mode of is [4]. Rare proceeds FCNC transitions. This decay is roughly of the same order as that of the [5]. In the SM framework, the rare decay is dominated by the Z-penguin and box diagrams involving top quark exchanges. The theoretical uncertainties related to the renormalization scale dependence of running quark mass can be essentially neglected after the inclusion of next-to-leading order corrections [6]. This decay is theoretically very clean process in comparison with the semileptonic decays like the decay and is also sensitive to the new physics beyond the SM [7]. Moreover, this decay is complementary to the decay. Note that the direct calculation of physical observables such as form factors suffers from sizable uncertainties. These can be greatly reduced through a combined analysis of the rare and [8] decays.

These decays have not yet been measured by the LHCb. There are no theoretical studies relevant to the form factors and branching ratios of decay. The form factors of these decays can be evaluated with the different approaches. Some of them are the light front, the constituent quark models [9], and the QCD sum rules. In this study the three-point QCD sum rules approach is used in the calculation of form factors. It is worth mentioning that the QCD sum rules have widely been utilized in calculation of the form factors (some of them can be found in [1017]).

The paper has 3 sections. In Section 2, the effective Hamiltonian and the three-point QCD sum rules approach are presented for completeness. In Section 3, the numerical values of form factors are given and the sensitivity of the branching ratio is studied and conclusion is presented.

2. Sum Rules for the Transition Form Factors

The FCNC decay is described within the framework of the SM at the quark level by the effective Hamiltonian [18, 19]: where is the Fermi constant, is the Weinberg angle, is the fine structure coupling constant, and

The is where . The explicit form of is given in [1820]. Note that gives about 3% contribution to the term [21].

The Wilson coefficients (in our case and ) can be calculated in any gauge and they are gauge independent and the results should be gauge invariant. The Wilson coefficients are calculated in gauge. It is worth mentioning that local operators in the considered problem have anomalous dimensions. We have checked that taking into account anomalous dimensions can change numerical results at most 10%.

The matrix element of the exclusive decays is found by inserting initial meson state and final meson state in (1): where is the polarization vector of meson, is the momentum of the , and is the momentum of meson. The matrix element of (4) is written in terms of the form factors as follows: Here, Lorentz invariant and parity conservation are considered. Also, , where are the dimensionless transition form factors. and is the transfer momentum or the momentum of the boson.

The matrix element in terms of the form factors is as where .

We try to calculate the aforementioned form factors by means of the QCD sum rules. The QCD sum rules begin with the following correlation functions: where the interpolating currents are and the and the meson states, respectively. consists of the vector () and axial vector () transition currents. After inserting the two complete sets of the and meson, the correlation functions in (7) are written as follows: where “” shows the contributions come from higher states and continuum of the currents with the same quantum numbers.

The and matrix elements are defined as follows: where and are the leptonic decay constants of and mesons, respectively. Using these equations and calculating the summation over the polarization of the vector meson , (8) is as follows: This correlation function is calculated in terms of the quarks and gluons parameters by means of the operator product expansion (OPE) as Each with and contains the perturbative and nonperturbative parts as in the following: The bare-loop diagram given in Figure 1(a) is the contribution of the perturbative part. The nonperturbative part consists of the two gluon condensates diagrams {see Figures 2(a)2(f)}. Hence, contributions of the light quark condensates {diagrams shown in Figures 1(b), 1(c), and 1(d)} vanish by applying the double Borel transformations [16].

The following double dispersion integrals are the contributions of the bare-loop diagrams in the correlation function: One of the basic methods to solve the Feynman Integrals in order to calculate the spectral densities is Cutkosky rules where the quark propagators are replaced by Dirac Delta Functions: , which indicates that all quarks are on-shell.

Three delta functions appear as a result of the applying Cutkosky rules. These delta functions have to vanish at the same time. Therefore, we get the following inequality from the arguments of the delta functions: where .

Following the standard calculations, the spectral densities are evaluated as where Now, it is aimed to calculate the nonperturbative part of (12) which consists of the gluon condensates diagrams shown in Figure 2. The gluon condensate contributions are calculated in Fock-Schwinger gauge because in this gauge the gluon field is expressed in terms of gluon field strength tensor directly. The following type of the integrals has to be calculated in order to get the results of the gluon condensate diagrams [15, 22]: where is the momentum of the spectator quark . The generic solutions for these integrals can be seen in [22, 23]. A part of our results for the contributions of the gluon condensate diagrams following similar methods shown in [22, 23] is given in the Appendix.

The Borel transformations are applied for both phenomenological and QCD side {see (11)} in order to suppress the contributions of higher states and continuum. The QCD sum rules for the form factors (, , , and ) are obtained by equalizing the Borel transformed forms of the physical side. The result is in the following formula: Note that the contributions of the gluon condensates () are already considered in the numerical analysis. However, each of these explicit expressions is extremely long; it is found unnecessary to show all of them in this study. Therefore, one of these expressions () is shown as a sample in Appendix. The and are the continuum thresholds in and channels, respectively. Also and .

3. Numerical Analysis

Having known the matrix element, that is, (6), the decay rate for decay is evaluated as follows:

The expression for the decay rate shows that we need to know the input parameters shown in Table 1, taken from [24].

Moreover, the values of the leptonic decay constants  GeV [25] and the gluon condensate  GeV4 [26] are necessary for the evaluation of the form factors. In addition, the form factors contain four auxiliary parameters: the Borel mass squares and and the continuum threshold and . The form factors are assumed to be independent of or weakly dependent on these auxiliary parameters in the suitable chosen regions named as “working regions.”

The contributions proportional to the highest power of are supposed to be less than about 30% of the contributions proportional to the highest power of . The lower bound of the and can be determined by the above condition. In addition, the contributions of continuum must be less than that of the first resonance. This helps us to fix the upper bound of the and . Therefore, we find the suitable region for the Borel mass parameters in the following intervals; and .

The numerical value of the and is supposed to be less than the mass squared of the first excited state meson with the same quantum numbers. In other words, the and are between mass squared of the ground sate meson and excited state meson with the same quantum numbers. The following regions for the and are chosen: and .

The form factors depend on the . The detail of the dependence is complicated. We fit them to the following function:

The , , and are given in Table 2.

The origin of the errors in Table 2 is the variation of , , and in the chosen intervals and the uncertainties of the input parameters.

In order to evaluate the branching ratio of the decay, the mean life time of the meson is needed. For the time being there is no experimental data on the mean life time of this meson. We follow the theoretical methods like Bethe-Salpeter model [27] and potential model [28] and estimate that the mean life time of the meson is in the order of the mean life time of the meson. We assume that the total life time  s [24]. Using the mean life time and the dependence of the form factors given by (20) in the kinematical allowed region we study the branching ratios for decay. Our results for three different values of the  GeV2 are presented in Table 3. In addition, Figure 3 depicts the dependence of the branching ratio on for full kinematical allowed region.

Finally, we calculate the integrated branching ratio for decay as follows:

To sum up, we investigated the branching ratio and decay rate of the decay. The form factors of this decay were found in the framework of the QCD sum rules. In addition, the contributions of the two gluon condensates diagrams to the correlations function were obtained.

Appendix

In this section, we present the explicit expression for the coefficients corresponding to the gluon condensates contributions of structure entering to the expression for the form factors in (18): where

Conflict of Interests

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