Research Article  Open Access
Decays with Perturbative QCD Approach
Abstract
Besides the traditional strong and electromagnetic decay modes, meson can also decay through the weak interactions within the standard model of elementary particle. With anticipation of copious data samples at the running LHC and coming SuperKEKB experiments, the twobody nonleptonic bottomchanging decays () are investigated with perturbative QCD approach firstly. The absolute branching ratios for and decays are estimated to reach up to about and , respectively, which might possibly be measured by the future experiments.
1. Introduction
The upsilon meson is the spintriplet wave state of bottomonium (bound state consisting of bottom quark and antibottom quark ) with wellestablished quantum number of [1]. The characteristic narrow decay widths of mesons for , 2, and 3 provide insight into the study of strong interactions (see Table 1, and note that, for simplicity, will denote , , and mesons in the following content if not specified definitely). The mass of meson is below meson pair threshold. meson decays into bottomed hadrons through strong and electromagnetic interactions are forbidden by the law of conservation of flavor number. The bottomchanging decays can occur only via the weak interactions within the standard model, although with tiny incidence probability. Both constituent quarks of upsilons can decay individually, which provide an alternative system for investigating the weak decay of heavyflavored hadrons. In this paper, we will study the nonleptonic ( and ) weak decays with perturbative QCD (pQCD) approach [2–4].
Experimentally, over data samples have been accumulated at Belle and BaBar experiments [5]. More and more upsilon data samples will be collected at the running hadron collider LHC and the forthcoming collider SuperKEKB (the SuperKEKB has started commissioning test run (http://www.kek.jp/en/NewsRoom/Release/)). There seems to exist a realistic possibility to explore weak decay at future experiments. Signals of decays should be easily distinguished with “charge tag” technique, due to the facts that the backtoback final states with different electric charges have definite momentum and energy in the rest frame of meson. meson has not been observed experimentally by now. meson production via the strong interaction is suppressed due to the simultaneous presence of two heavy quarks with different flavors and higher order in QCD coupling constant . , decays provide a novel pattern to study meson production. The identification of a single explicitly flavored meson could be used as an effective selection criterion to detect upsilon weak decays. Moreover, the radiative decay of meson provides a useful extra signal and a powerful constraint (the investigation on the radiative decay of meson can be found in, e.g., [6], with QCD sum rules). Of course, any discernible evidences of an anomalous production rate of single bottomed meson from upsilon decays might be a hint of new physics.
Theoretically, many attractive QCDinspired methods have been developed recently to describe the exclusive nonleptonic decay of heavyflavored mesons, such as the pQCD approach [2–4], the QCD factorization approach [7–9], and soft and collinear effective theory [10–13], and have been applied widely to vindicate measurements on meson decays. The upsilon weak decay permits one to further constrain parameters obtained from meson decay, and cross comparisons provide an opportunity to test various phenomenological models. The upsilon weak decay possesses a unique structure due to the CabibboKobayashiMaskawa (CKM) matrix properties which predicts that the channels with one meson are dominant. decay belongs to the favorable transition, which should, in principle, have relatively large branching ratio among upsilon weak decays. However, there is still no theoretical study devoted to decay for the moment. In this paper, we will present a phenomenological investigation on weak decay with the pQCD approach to supply a ready reference for the future experiments.
This paper is organized as follows. Section 2 focuses on theoretical framework and decay amplitudes for weak decays. Section 3 is devoted to numerical results and discussion. The last section is a summary.
2. Theoretical Framework
2.1. The Effective Hamiltonian
Theoretically, , weak decays are described by an effective bottomchanging Hamiltonian based on operator product expansion [19]: where [1] is the Fermi coupling constant; the CKM factors and correspond to and decays, respectively; with the Wolfenstein parameterization, the CKM factors are expanded as a power series in a small Wolfenstein parameter [1]: The local tree operators are defined as where and are color indices and the sum over repeated indices is understood.
The scale factorizes physics contributions into short and longdistance dynamics. The Wilson coefficients summarize the physics contributions at scale higher than and are calculable with the renormalization group improved perturbation theory. The hadronic matrix elements (HME), where the local operators are inserted between initial and final hadron states, embrace the physics contributions below scale of . To obtain decay amplitudes, the remaining work is to calculate HME properly by separating from perturbative and nonperturbative contributions.
2.2. Hadronic Matrix Elements
Based on LepageBrodsky approach for exclusive processes [20], HME is commonly expressed as a convolution integral of hard scattering subamplitudes containing perturbative contributions with universal wave functions reflecting nonperturbative contributions. In order to effectively regulate endpoint singularities and provide a naturally dynamical cutoff on nonperturbative contributions, transverse momentum of valence quarks is retained and the Sudakov factor is introduced within the pQCD framework [2–4]. Phenomenologically, pQCD’s decay amplitude could be divided into three parts: the Wilson coefficients incorporating the hard contributions above typical scale of , processdependent rescattering subamplitudes accounting for the heavy quark decay, and wave functions of all participating hadrons, which is expressed as where is the momentum of valence quarks and is the Sudakov factor.
2.3. Kinematic Variables
The light cone kinematic variables in rest frame are defined as follows: where and are the longitudinal momentum fraction and transverse momentum of valence quarks, respectively; and are the longitudinal and transverse polarization vectors, respectively, and satisfy relations and ; the subscript on variables , , , and corresponds to participating hadrons; namely, for meson, for the recoiled meson, and for the emitted pseudoscalar meson; and are positive and negative null vectors, respectively; , , and are the Lorentzinvariant variables; is the common momentum of final states. The notation of momentum is displayed in Figure 2(a).
2.4. Wave Functions
With the notation in [16, 21], wave functions are defined as where , , and are decay constants of , , and mesons, respectively.
Considering mass relations of and , it might assume that the motion of heavy valence quarks in and mesons is nearly nonrelativistic. The wave functions of and mesons could be approximately described with nonrelativistic quantum chromodynamics (NRQCD) [22–24] and timeindependent Schrödinger equation. For an isotropic harmonic oscillator potential, the eigenfunctions of stationary state with quantum numbers are written as [17] where parameter determines the average transverse momentum; that is, . Employing the substitution ansatz [25], where and are the longitudinal momentum fraction and mass of valence quark, respectively; then integrating out and combining with their asymptotic forms, the distribution amplitudes (DAs) for and mesons can be written as [17] where ; . According to NRQCD power counting rules [22], with and QCD coupling constant . The exponential function represents distribution. Parameters of , , , , , , , and are normalization coefficients satisfying the conditions
The shape lines of normalized DAs for and mesons are shown in Figure 1. It is clearly seen that DAs for and mesons fall quickly down to zero at endpoint , due to suppression from exponential functions; DAs for meson are symmetric under the interchange of momentum fractions , and DAs for meson are basically consistent with the feature that valence quarks share momentum fractions according to their masses.
(a)
(b)
(c)
(a)
(b)
(c)
(d)
Our study shows that only the leading twist (twist2) DAs of the emitted light pseudoscalar meson are involved in decay amplitudes (see Appendix). The twist2 DAs have the expansion [16],and are normalized as where are Gegenbauer polynomials: and each term corresponds to a nonperturbative Gegenbauer moment ; note that due to the normalization condition equation (12); parity invariance of the pion DAs requires Gegenbauer moment for .
2.5. Decay Amplitudes
The Feynman diagrams for weak decay are shown in Figure 2. There are two types. One is factorizable emission topology where gluon attaches to quarks in the same meson, and the other is nonfactorizable emission topology where gluon connects to quarks between different mesons.
With the pQCD master formula equation (4), the amplitude for decay can be expressed as [26] which is conventionally written as the helicity amplitudes [26]: where and the color number ; the subscript on corresponds to three different helicity amplitudes; that is, ; the superscript on denotes indices of Figure 2. The explicit expressions of building blocks are collected in Appendix.
3. Numerical Results and Discussion
In the centerofmass of meson, branching ratio for decay is defined as
The input parameters are listed in Tables 1 and 2. If not specified explicitly, we will take their central values as the default inputs. Our numerical results are collected in Table 3, where the first uncertainty comes from scale and the expression of is given in (A.4) and (A.5); the second uncertainty is from mass of and ; the third uncertainty is from hadronic parameters including decay constants and Gegenbauer moments; the fourth uncertainty is from CKM parameters. The followings are some comments:(1)Branching ratio for decay is about with pQCD approach, which is well within the measurement potential of LHC and SuperKEKB. For example, experimental studies have showed that production cross sections for meson in pp and pPb collisions are a few at the LHCb [27, 28] and ALICE [29, 30] detectors. Consequently, there will be more than data samples per data collected by the LHCb and ALICE, corresponding to a few hundreds of events. Branching ratio for decay, , is generally less than that for decay by one order of magnitude due to the CKM suppression, .(2)As it is well known, due to the large mass of , the momentum transition in decay may be not large enough. One might naturally wonder whether the pQCD approach is applicable and whether the perturbative calculation is reliable. Therefore, it is necessary to check what percentage of the contributions comes from the perturbative region. The contributions to branching ratio for decay from different region are shown in Figure 3. It can be clearly seen that more than 93% (97%) contributions come from (0.3) region, implying that decay is computable with the pQCD approach. As the discussion in [2–4], there are many factors for this, for example, the choice of the typical scale, retaining the quark transverse moment and introducing the Sudakov factor to suppress the nonperturbative contributions, which deserve much attention and further investigation.(3)Because of the relations among masses resulting in the fact that phase space increases with the radial quantum number in addition to the relations among decay widths , in principle, there should be relations among branching ratios for the same pseudoscalar meson . But the numerical results in Table 3 are beyond such expectation. Why? The reason is that the factor of in (16) has almost the same value for , so branching ratio is proportional to factor with the maximal value for . Besides, contributions from regions decrease with (see Figure 3), which enhance the decay amplitudes.(4)Besides the uncertainties listed in Table 3, other factors, such as the models of wave functions, contributions of higher order corrections to HME, and relativistic effects, deserve the dedicated study. Our results just provide an order of magnitude estimation.
 
The decay constant cannot be extracted from the experimental data because of no measurement on weak decay at the present time. Theoretically, the value of has been estimated, for example, in [18], with the QCD sum rules. From Table of [18], one can see that the value of is modeldependent. In our calculation, we will take the latest value given by the lattice QCD approach [15] just to offer an order of magnitude estimation on branching ratio for decays. 

(a)
(b)
(c)
4. Summary
decay via the weak interaction, as a complementary to strong and electromagnetic decay mechanism, is allowable within the standard model. Based on the potential prospects of physics at highluminosity collider experiment, decay into and final states is investigated with the pQCD approach firstly. It is found that the dominant contributions come from perturbative regions , which might imply that the pQCD calculation is practicable and workable; there is a promiseful possibility of searching for () decay with branching ratio about () at the future experiments.
Appendix
Building Blocks for Decays
The building blocks , where the superscript corresponds to indices of Figure 2 and the subscript relates to different helicity amplitudes, are expressed as follows: