Abstract

We present the calculation of the form factor of the and decays in the frame work of QCD sum rules. We also find branching ratio = which is in agreement with the experimental data. Furthermore, we estimate the  MeV, where experimental bound for full width of is  MeV.

1. Introduction

Heavy quarkonium states (like and ) and their decay modes offer a laboratory to study the strong interaction in the nonperturbative regime. Charmonium in particular has served as a calibration tool for the corresponding techniques and models [1, 2]. Heavy quarkonium states can have many bound states and decay channels used to study and determine different parameters of standard model (SM) and QCD from the theoretical perspective. In particular, the calculation of bottomonium masses [3], total widths, coupling constants [47], and branching ratio can serve as benchmarks for the low energy predictions of QCD. In addition, the theoretical calculations on the branching ratio of radiative decays of heavy quarkonium states are relatively clean with respect to the hadronic or semileptonic decays, and their comparison with experimental data could provide important insights into their nature and hyperfine interaction. In this regard, the decay width and branching ratio of radiative decays of several and states analysed by using an effective Lagrangian approach valid for heavy quarkonia [8, 9].

The QCD sum rules have been used for the radiative transitions in charmonium and bottomonium cases [1015]. But exclusive and decays have not yet been studied in the framework of the QCD sum rules. The decay studied by CLEO detector operating at CESR [16, 17]. The next experimental studies of exclusive and decays are going to start with next operation of LHC.

We present the theoretical study on the form factor of exclusive and decays in the frame work of QCD sum rules. Note that in order to calculate the branching ratio we have to acquire information about the masses and decay constants of the participating particles. It is worth mentioning that the masses can be obtained either by means of the experimental results, that is, the particle data group or by the theoretical methods. The decay constants, on the other hand, can be calculated theoretically via different nonperturbative methods. In this respect, masses and decay constants and spectrum of heavy quarkonium states are calculated in the various approaches (see, e.g., [1, 1825]). Here, firstly, we calculate the form factor of and decays in the framework of three-point QCD sum rules; secondly, we calculate the branching ratio of exclusive and decays.

In Section 2, we introduce the QCD sum rules technique for the form factors of exclusive and decays. Last section is devoted to the numerical analysis and discussion.

2. QCD Sum Rules for the Form Factors

The three-point correlation function associated with the and vertex is given by where is the time ordering operator and is momentum of photon. Each meson and photon field can be described in terms of the quark field operators as follows:

We calculate the correlation function equation (1) in two different methods. In phenomenological or physical approach, it can be evaluated in terms of hadronic parameters such as masses, decay constants, and form factors. In theoretical or QCD side, on the other hand, it is calculated in terms of QCD parameters, which are quark and gluon degrees of freedom, by the help of the operator product expansion (OPE) in deep Euclidean region. Equating the structure calculated in two different approaches of the same correlation function, we get a relation between hadronic parameters and QCD degrees of freedom. Finally, we apply double Borel transformation with respect to the momentum of initial and final mesons ( and ). This final operation suppresses the contribution of the higher states and continuum.

2.1. Phenomenological Side

We insert the complete sets of appropriate vector meson ( ) and scalar meson ( ) states (regarding the conservation of the quantum numbers of corresponding interpolating currents) inside correlation functions equation (1). Here, vector state is either or and scalar state is state. After integrating over the and , we get the following result for the correlation function equation (1): where contains the contribution of the higher and continuum states with the same quantum numbers.

The matrix elements of the above equation are related to the hadronic parameters as follows: where is the form factor of transition and is the polarization vector associated with the vector meson. Using (4) in (3) and considering the summation over polarization vectors via the result of the physical side is as follows: We are going to compare the coefficient of structure for further calculation from different approaches of the correlation functions.

2.2. Theoretical (QCD) Side

Theoretical side consists of perturbative (bare loop; see Figure 1) and nonperturbative parts (the contributions of two gluon condensate diagrams, Figure 2). We calculate it in the deep Euclidean space ( and ). We consider this side as

2.2.1. Bare Loop

The perturbative part is a double dispersion integral as follows: where is called spectral density. We aim to evaluate the spectral density with the help of the bare loop diagram in Figure 1. One of the generic methods to calculate this bare loop integral is the Cutkosky method, where the quark propagators of Feynman integrals are replaced by the Dirac delta functions: Then, using the Cutkosky method we get spectral density as where and is the color number. Note that, since three functions of integrand must vanish simultaneously, the physical regions in the plane must satisfy the following inequality:

2.2.2. Two Gluon Condensates

We consider the two gluon condensate diagrams. Note that we do not include the heavy quarks condensate diagrams, since the heavy quark contributions are exactly reducible to the gluon condensate [26]. Now, as a nonperturbative part, we must add contributions coming from the gluon condensates presented in Figures 2(a), 2(b), 2(c), 2(d), 2(e), and 2(f).

These diagrams are calculated in the Fock-Schwinger fixed-point gauge [2729], where the vacuum gluon field is as follows: where is the gluon momentum and is the gluon field. In addition, the quark-gluon-quark vertex is used as

We come across the following integrals in calculating the gluon condensate contributions [30, 31]: where is the momentum of the spectator quark .

These integrals are calculated by shifting to the Euclidean space-time and using the Schwinger representation for the Euclidean propagator: This kind of expression is very easy for the Borel transformation since where is Borel parameter.

We perform integration over the loop momentum and over the two parameters which we use in the exponential representation of propagators [28]. As a final operation we apply double Borel transformations to and . We get the Borel transformed form of the integrals in (14) as where and and are the Borel parameters. The function is as follows: where where the circumflex of in the equations is used for the results after the double Borel transformation. As a result of the lengthy calculations we obtain the following expressions for the two gluon condensate:

Now, we can compare coefficient of (6) and (7). Our result related to the sum rules for the corresponding form factor is as follows: Note that finally we have to set for the real photon.

3. Numerical Analysis

In this section we calculate the value of form factors and the branching ratios. We use  GeV and  GeV [32], which correspond to the pole masses  GeV and  GeV [6]. Also,  MeV [32],  MeV [32],  MeV [32],  MeV [32],  MeV [33],  MeV [33],  MeV [3],  MeV [3], and the full width for :  MeV [32] are used.

To do further numerical analyses we have to know the value or range of the auxiliary parameters of QCD sum rules. Those are the continuum thresholds ( and ) and the Borel mass parameters ( and ). The physical results are required to be either weakly dependent on or independent of the aforementioned parameters. Therefore, we must consider the working regions of these auxiliary parameters where the dependence of the form factors is weak. We also consider the working regions for the Borel mass parameters and in a way that both the contributions of the higher states and continuum are sufficiently suppressed and the contributions coming from higher dimensions operators can be ignored. With the aforementioned conditions, we find the stable region for the form factor in the following intervals: and for decays (see also Figures 3 and 4). We also get and for decays.

The continuum thresholds, and , are fixed by the mass of the corresponding ground-state hadron. Note that they must not be greater than the energy of the first excited states with the same quantum numbers. In our numerical calculations the following regions for the continuum thresholds in and channels are used: and for and channels, respectively. Here, is the mass of either or meson and is the mass of either or meson. We studied the dependence of the form factor ( ) on continuum thresholds (see Figures 3 and 4). The figures show the weak dependence of the form factor for the chosen intervals.

Note that we follow the standard procedure in the QCD sum rules, where the continuum thresholds are supposed to be independent of Borel mass parameters and . However, this assumption is not free of uncertainties (see, e.g., [34]).

Considering the large negative enables us to evaluate the correlation function by means of OPE that gives better convergence property; therefore, we use the extrapolation of the form factor from negative into the physical region (positive region).

The best fit curve for is employed for the negative region and fit is extrapolated for the positive region (see Figure 5). The values , , and for and , , and for decays are obtained via (23).

Using in (23), we obtain the and the for and decays, respectively. It is worth mentioning that firstly, the contributions of the two gluon condensate in the value of the is about . Secondly, roughly 80% of the errors in our numerical calculation arise from the variation continuum thresholds, Borel mass parameter in the given intervals, and uncertainties of the input parameters and the remaining 20% occurs as a result of the quark masses when one proceeds from the to the pole-scheme mass parameters, the input parameters.

The matrix element for the decay of and is as follows: where and are the momentum and polarization of final state vector meson, that is, either or mesons, is the polarization of real photon, and is the momentum of initial scalar meson.

Using this matrix element, we get where is mass of either or meson and is either mass of or meson. The decay width for decays is

The branching ratio of can be evaluated with (26) and using the experimental total width that is This result is in good agreement with the experimental measurement [32] and the result of the potential model given in [8, 9], which are respectively.

We get for decays. Using this value we calculate the decay width as follows: This decay width and the measured branching ratio [17] allow us to evaluate the total width of . We estimate that the full width , which is consistent with the experimental results that indicate the full width [17].

Conflict of Interests

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

Acknowledgments

The author is grateful to T. M. Aliev and K. Azizi for their useful discussions.